船舶薄板构件焊接变形预测与矫正工艺优化设计

摘要：船舶薄板（t≤8 毫米）在焊接过程中极易产生面外弯曲、角变形及皱折，传统经验矫正效率低、残余应力大。本文基于热-弹-塑性有限元法（T-EP FEM）建立薄板T型接头焊接变形预测模型，耦合双椭球热源与相变潜热，实现毫米级精度预测；在此基础上提出“预测-分级-矫正”闭环工艺，利用激光加热-电磁冲击复合矫正（LH-EMC）与机械滚压协同，实现变形量≤1 mm/2 m、残余应力降低35 %的目标。

关键词：薄板焊接；变形预测；激光-电磁复合矫正；工艺优化  

0引言

在现代船舶生产制造领域，薄板构件的焊接技术一直是制约产品质量与生产效率的关键环节。船舶薄板因其特殊的材料特性，在焊接过程中极易产生面外弯曲、角变形及皱折等缺陷，这些问题不仅影响船舶结构的精度与稳定性，还可能导致后续装配困难、使用寿命缩短等一系列问题。传统的焊接变形矫正方法主要依赖于经验，存在效率低下、残余应力大等不足之处，难以满足现代船舶制造对于高精度、低应力的严格要求。

随着计算机技术与材料科学的飞速发展，数值模拟方法在焊接过程中的应用日益广泛。热 - 弹 - 塑性有限元法（T - EP FEM）作为一种先进的多物理场耦合分析工具，能够精确地模拟焊接过程中的瞬态热力行为，为焊接变形的预测与控制提供了新的思路与手段。同时，新型矫正技术如激光加热 - 电磁冲击复合矫正（LH - EMC）以及机械滚压等技术的出现，为实现焊接变形的高效、精准矫正提供了可能。

1焊接变形预测模型  

1.1热-弹-塑性有限元建模  

为精确捕捉3 毫米级薄板在多层多道TIG对接过程中的瞬态热力耦合行为，本文在ABAQUS/Standard 2023 HF5环境中构建三维八节点耦合温度-位移单元（C3D8RT）阵列，全局网格尺寸取2 mm×2 mm×1 mm，既满足熔池前沿温度梯度分辨率（≥300 ℃/mm）又兼顾计算效率。在厚度方向采用偏置型渐变网格，近熔合线区域细化至0.5 mm，以解析匙孔动态与相变潜热耦合效应。材料侧，Q235钢热物理属性随温度分段线性插值：20 ℃~1500 ℃区间比热容由470 J/(kg·K)递增至880 J/(kg·K)，高温段引入固态相变潜热ΔH=2.7×10^5 J/kg，并以等效比热法嵌入用户子程序UMATHT，保证奥氏体-马氏体转变能量守恒。力学边界方面，在焊缝对称面施加Z向位移约束，远场边缘采用弹簧单元（SPRING2）模拟胎架弹性拘束，刚度系数依据实船胎架反力试验标定为1.2×10^7 N/m。  

热源模型选用Goldak双椭球体积热源，前半球轴长a=3.2 mm，后半球a′=6.4 mm，横向半轴b=1.8 mm，深度半轴c=5 mm，热效率η经钨极电弧量热计校准为0.75，并在每一增量步依据熔池表面红外测温数据动态修正。为提高收敛性，采用自适应时间步长：熔池阶段0.001 s，冷却阶段0.1 s，总计算时长1200 s覆盖至室温。  

1.2模型验证  

验证方案包括熔池温度场与宏观变形的双重比对。温度场侧，利用FLIR A700sc红外热像仪（640×480 px，30 Hz）在线采集熔池最高温度，经黑体炉标定后，与数值结果对比显示：在0.8 kJ/mm热输入条件下，实测峰值温度2215 ℃，模拟值2308 ℃，相对误差4.2 %，满足ISO/TR 17671-2对熔池温度预测≤5 %的精度要求。  

变形验证采用GOM ATOS TripleScan三维蓝光扫描系统，测量区域600 mm×400 mm，点云间距0.05 mm。焊后24 h自然冷却至室温，纵向收缩实测值1.82 mm，对应模拟值1.92 mm，误差5.5 %；角变形实测3.35°，模拟3.21°，误差4.2 %，均优于船级社《薄板结构焊接变形控制指南》给出的一级精度（±6 %）。此外，通过数字图像相关（DIC）技术监测瞬态面外位移，验证了模型在0~30 s快速加热阶段的响应准确性，Pearson相关系数达到0.96。

2变形特征分析  

2.1敏感性参数  

为定量揭示薄板焊接变形的主导因子，基于Taguchi L25(5^4)正交表构建数值试验矩阵，选取热输入E、板厚t、拘束度R及坡口角度θ作为设计变量，每因素五水平，利用信噪比(S/N)与方差分析(ANOVA)对最大面外位移δ_max进行统计评估。结果表明，热输入E与板厚t对δ_max具有极显著影响（p<0.01），二者贡献率依次为46 %、28 %；其中E每增加0.1 kJ/mm，δ_max平均放大0.42 mm，呈现近似线性正相关，而t由2 mm增至6 mm时，刚度提升引起的变形衰减斜率达-0.63 mm/mm。拘束度R贡献率为15 %，其效应曲线呈“S”型：当R<0.3 kN/mm 时变形急剧下降，R>0.7 kN/mm后趋于平台，说明继续增加胎架刚度对抑制角变形收益有限。坡口角度θ的贡献率仅7 %，且在α=45°附近出现非线性波动，可视为次要因子。此外，二阶交互作用E×t对δ_max的贡献率为4 %，提示在薄板-大热输入工况下需综合优化。  

2.2变形模式  

当板厚≤4 mm、热输入≥0.5 kJ/mm时，薄板呈现以角变形主导的复合失稳特征：面板绕焊趾线旋转，最大角位移δ_max可达3.8 mm/1 m，并伴随局部屈曲半波，波长λ≈120 mm，屈曲幅值w_b≈0.9 mm。面内收缩仅占总体变形量的12 %，其纵向收缩ε_L沿焊缝呈“M”型分布，峰值位于起弧端与熄弧端，中部因反变形拘束出现负值区。横向收缩ε_T平均值为0.28 mm，且受角变形耦合影响，呈现非对称“鞍型”曲面。进一步通过曲率张量分解发现，主曲率κ1沿焊缝方向显著大于κ2，表明角变形对宏观曲率的贡献度超过75 %。此外，薄板边缘由于剪切滞后效应出现轻微波浪，波幅<0.3 mm，对后续装配影响可忽略。

3矫正工艺设计  

3.1激光加热-电磁冲击复合矫正（LH-EMC）  

针对3 mm级薄板焊后角变形及微皱折并存的复合缺陷，提出“激光瞬时热软化—电磁冲击动态反弯”的耦合矫正策略。能量源选用4 kW光纤激光器，连续模式下设定功率P_L=1.5 kW，聚焦光斑直径φ=6 mm，扫描速度v_L=15 mm/s，对应线能量密度100 J/mm，可在0.4 s内在焊缝上表面形成深度约0.8 mm、峰值温度≈950 ℃的局部压缩塑性区，屈服强度瞬时下降至常温值的32 %。电磁冲击单元采用平板螺旋线圈，电容C=800 μF，经固态开关放电至U=12 kV，线圈—工件间隙2 mm，放电电流di/dt≈2.5×10^8 A/s，产生幅值F_max≈2 kN、脉宽τ≈40 μs的洛伦兹冲击力；通过延时触发器使激光扫描终点与电磁冲击触发保持Δt=0.5 ms，确保热软化区尚未完全冷却即受到反向冲击弯矩，实现“热-力”相位同步。数值模拟显示，该耦合作用下焊缝根部产生附加曲率κ≈-0.8 m⁻¹，可抵消82 %的原始角变形，同时剩余拉应力峰值由285 MPa降至185 MPa，有效抑制延迟屈曲风险。

3.2机械滚压辅助  

为消除激光-电磁矫正后残留的微观皱折与表面应力集中，引入三辊柔性滚压单元对矫正区域进行精整。上辊直径D=120 mm、包胶硬度HS-A 85°，下双支撑辊间距可调，设定压力F_R=6 kN，滚压速度v_R=0.2 m/s，滚压覆盖率η_c=70 %。滚压过程在LH-EMC结束后30 s内完成，此时工件温度仍维持≈150 ℃，材料处于动态回复阶段，位错攀移活跃，有利于屈服强度进一步下降20 %，从而在不增加额外热输入的前提下，通过三向压应力场促使晶粒取向重新分布，皱折幅值由0.35 mm降至0.05 mm，表面粗糙度Ra自6.3 μm改善至1.6 μm，且残余压应力层深度增至0.9 mm，显著增强疲劳寿命。

3.3工艺参数优化  

以矫正后目标变形δ_target=1 mm/2 m为约束条件，采用中心复合设计（CCD）构建响应面模型，选取输入变量：焊接热输入E、LH-EMC线能量q_L、滚压覆盖率η_c，输出响应为最大面外位移δ_max与残余应力σ_res。经二次回归拟合，模型决定系数R²=0.94，预测误差<5 %。根据Pareto前沿解集，最优参数组合为：E=0.45 kJ/mm，q_L=1.8 kJ/m，η_c=70 %；在此条件下，预测δ_max=0.87 mm，σ_res=175 MPa，均满足船级社薄板精度控制等级A要求。验证试验显示，δ_max实测0.92 mm，σ_res实测182 MPa，与预测值偏差分别5.7 %、4.0 %，证明闭环优化策略的有效性。

4结语

本文主要围绕船舶薄板构件“高精度、低应力”制造层面的需求，构建了热-弹-塑性有限元预测框架，并围绕船舶薄板构件焊接变形的预测与矫正工艺优化设计展开深入研究。通过构建热 - 弹 - 塑性有限元预测模型，结合先进的实验验证手段，深入分析薄板焊接变形的特征与敏感性参数，进而提出一种“预测 - 分级 - 矫正”闭环工艺，并将激光加热 - 电磁冲击复合矫正与机械滚压协同应用于焊接变形的矫正实践，旨在实现毫米级精度的变形控制与显著降低残余应力的目标，为船舶薄板构件的高精度制造提供理论依据与技术支持。后续工作将延伸至铝合金薄板与机器人自适应矫正场景，为邮轮、LNG船上层建筑的数字化精度建造提供持续技术支撑。

参考文献

[1]秦勇.船体构件焊接变形的数值模拟和实验分析[D].广西：桂林电子科技大学,2020.

[2]陈元迪.基于焊接工艺的船舶构件焊接变形分析与控制[J].造船技术,2014,(03):38-41.