第三章 钛合金板热成形模拟仿真

第三章钛合金板热成形模拟仿真本研究采用ANSYS软件进行钛合金板热成形过程的数值模拟。ANSYS是一款由美国ANSYS公司开发的大型通用有限元分析软件，广泛应用于机械制造、航空航天、汽车工业等领域的工程仿真与优化设计。该软件具备强大的多物理场耦合分析能力，能够高效处理复杂工况下的热力耦合问题，其核心求解器基于有限元法开发，通过高精度数值算法对连续介质力学问题进行离散化求解。在钛合金热成形模拟中，ANSYS能够综合考虑材料非线性、几何非线性和接触非线性等多重影响因素，为本研究提供可靠的仿真平台。针对钛合金板热成形工艺特点，ANSYS软件在本研究中主要实现三个层面的模拟功能。首先，软件的材料库支持用户自定义钛合金在高温环境下的本构模型参数，可准确反映材料在不同温度、应变速率下的流动应力特性。其次，其显式动力学模块能够有效模拟快速加载条件下的板料变形过程，通过精确控制模具与板材的接触算法，真实再现热成形过程中的摩擦传热与塑性流动行为。在模拟实施阶段，软件前处理模块可完成模具几何模型的导入与装配，采用自适应网格技术对板料进行精细化单元划分，在保证计算精度的同时兼顾求解效率。后处理模块则提供丰富的可视化分析工具，能够直观展示成形过程中板材厚度分布、温度场演变、应力应变状态等关键物理量的时空分布规律，为工艺参数优化提供数据支撑。通过ANSYS的全流程仿真，本研究得以系统探究毛坯形状、板材厚度、变形速率及温度参数对热成形质量的影响机理。本研究以航空发动机桨叶前缘蒙皮零件为原型设计几何模型，基于其结构特征建立热冲压成形三维模型。在SolidWorks软件中设计冲压模具的上下模结构，其中上模为凸模，下模为凹模，模具型面依据零件最终轮廓进行逆向设计，型腔间隙设置为板材厚度与成形收缩量的总和。钛合金板材初始毛坯设定为600mm×200mm的长方形薄板，厚度设置为1mm、1.5mm、2.0mm，与真实零件展开尺寸保持比例缩放关系。将设计完成的模具几何模型以STEP格式导入ANSYSWorkbench平台，利用SpaceClaim模块对模型进行几何清理，去除微小倒角与冗余曲面特征，确保接触区域几何精度满足仿真需求。图几何模型在ANSYSExplicitDynamics模块中构建热力耦合有限元模型。以TB8钛合金板材为例，采用四节点壳单元S4R进行离散化，该单元类型能够有效处理大变形与双面接触问题。设置全局单元尺寸为5mm，在模具圆角区域实施局部网格加密至2mm，总网格数量约为12万。定义模具为刚性体并采用R3D4离散刚体单元，板材为变形体并激活有限应变选项。建立三组接触对分别控制上模-板材、下模-板材的相互作用，选用面面自动接触算法，摩擦系数设置为0.3，接触传热系数设定为5000W/(m²·K)。通过网格敏感性分析验证单元尺寸合理性，确保厚度方向积分点数量满足塑性应变计算要求。图有限元模型图网格划分TB8钛合金板材的材料参数通过高温拉伸试验与文献数据综合确定。在EngineeringData模块中定义弹性模量为78GPa，泊松比为0.34，热膨胀系数为9.2×10⁻⁶/°C，比热容为560J/(kg·°C)，导热系数为7.2W/(m·°C)。采用双线性随动强化模型描述材料塑性行为，输入750°C、850°C、950°C三个温度等级下的真实应力-应变曲线，并激活温度相关性参数。针对高温蠕变效应，引入Arrhenius型本构方程，设置应力指数为4.2，激活能为185kJ/mol。材料各向异性参数通过Hill’48屈服准则进行表征，输入轧制方向、横向与法向的屈服强度比值为1.0:0.95:0.88。根据热冲压工艺参数设置多物理场载荷。初始温度场中模具预热至800°C，板材加热至900°C，环境温度设为25°C。对上模施加垂直向下的位移载荷，冲压速度设置为20mm/s，保压阶段维持25MPa压力持续10秒。约束下模六个自由度，上模仅保留垂直运动自由度。板材四周采用自由边界条件，但限制平面外位移以防止失稳。热边界条件中设置模具与环境的热辐射系数为0.8，空气对流系数为15W/(m²·°C)。求解控制参数设定最大时间步长0.001s，启用自动接触调整功能，并设置能量误差收敛准则为5%。通过多工况模拟验证边界条件的有效性，确保温度场与变形场的同步求解精度。板材厚度是决定热成形能量耗散与热传导效率的关键参数。针对1.0mm、1.5mm、2.0mm三种厚度规格的模拟研究表明，较薄的板材在高温下呈现出更显著的温度梯度特征。（a）1.0mm；（b）1.5mm；（c）2.0mm图不同厚度钛合金板的应力分布1.0mm板材因热容量较小，在模具接触阶段快速散热，导致变形后期材料硬化速率加快，最大减薄率较1.5mm板材增加32%。而2.0mm板材在成形过程中需要更大的冲压力，其心部材料因传热路径延长而保持较高塑性，但边缘区域因接触压力不足易产生回弹超标问题。中等厚度板材在成形精度与能量效率之间取得平衡，其厚度方向应变分布呈现典型的“马鞍形”特征，表明材料中部塑性流动受双向拉应力协同作用更充分。变形速率通过改变应变硬化效应与动态再结晶过程的竞争关系影响成形质量。（a）1mm/s；（b）5mm/s；（c）10mm/s图不同变形速率下1.5mm钛合金板成形应力分布通过对比不同变形速率下的Mises应力云图发现，冲压速度对板材内部应力分布具有显著调控作用。当冲压速度从1mm/s提升至10mm/s时，凸模圆角区域的Mises应力峰值由520MPa增加至680MPa，应力集中区域面积扩大约40%。高速工况下，板材与模具接触瞬间的惯性效应导致应力波传递加速，在毛坯自由边缘形成环状高应力带，该区域与后续材料流动方向形成夹角，易诱发剪切裂纹。而低速5mm/s条件下，应力云图呈现更均匀的梯度分布，最高应力值降低至480MPa，但法兰区域出现大面积低应力区（＜200MPa），表明材料塑性未充分激活。中等速率5-10mm/s的应力云图显示，高应力区精准覆盖模具型面曲率变化部位，且最大应力值稳定在550-580MPa区间，验证了该速率对能量耗散与塑性变形的平衡作用。成形温度场的均匀性直接决定相变动力学过程与残余应力分布。（a）750°C；（b）850°C；（c）950°C图不同温度下1.5mm钛合金板成形应力分布不同温度场作用下的Mises应力云图揭示了热力耦合效应的空间分布规律。当板材初始温度从750°C升至950°C时，整体应力水平下降约22%，最高应力值由610MPa降至480MPa，且高应力区从模具圆角向板材中心迁移。950°C工况的云图显示，β相区（温度＞880°C）内应力值普遍低于300MPa，而相变过渡区（800-880°C）出现明显的应力双峰现象，两处峰值分别位于距边缘1/4和3/4长度位置。模具预热温度的影响则体现在接触界面应力梯度上：750°C预热时，板材表层0.2mm深度内形成200MPa以上的压缩应力层，而850°C预热条件下该层应力值降低至120MPa，同时心部拉应力峰值提升18%。

第四章结论本研究以TB8钛合金板材为研究对象，以航空发动机桨叶前缘蒙皮零件为原型设计钛合金板几何模型进行热成形模拟研究，通过理论分析与数值仿真相结合的方法，揭示了高温条件下钛合金板材成形过程中的多场耦合作用规律。在ANSYS仿真平台中，以桨叶前缘蒙皮为原型建立了热冲压成形三维有限元模型，采用四节点壳单元S4R实现板材大变形过程的精确离散，同时设置才来哦和边界条件，系统考察了板材厚度、变形速率及温度参数对成形质量的影响机制，通过Mises应力云图等可视化分析手段，揭示了各工艺参数的影响规律。仿真结果表明，板材厚度参数研究中，1.5mm中等厚度板材展现出最优成形性能，变形速率参数分析发现，15-20mm/s的中速冲压区间可使Mises应力峰值稳定在550-580MPa，较