TCWAN 0100-2023 焊接数值模拟固有应变法

ICS25.160.01CCSJ33团 体 标 准T/CWAN0100—2023焊接数值模拟固有应变法Numericalweldingsimulationinherentstrainmethod2023-12-27发布 2024-02-01实施中国焊接协会 发布T/CWAN0100—2023T/CWAN0100—2023PAGE\*ROMANPAGE\*ROMANI目 次前 言 II范围 1规范性引用文件 1术语和定义 1基本流程 2一般要求 2焊接仿真分析模型建立规则 2仿真结果评估 6仿真分析报告 7附录A（规范性）基于固有应变法焊接数值模拟流程图 8附录B（资料性）固有应变的提取及数据修正 9附录C（资料性）文档模板 16附录D（资料性）对接接头MIG焊接的变形模拟仿真报告 17前 言本文件按照GB/T1.1—2020《标准化工作导则 第1部分：标准化文件的结构和起草规则》的规定起草。本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中国焊接协会提出并归口。本文件主要起草人：彭根琛、陆智俊、龙伟民、杨帆、张立平、魏世忠、张浩、冯吉才、张贵芝、黄瑞生、樊喜刚、纪昂、武汉琦、孟政宇、钟素娟、武鹏博、宋玉军、徐殿鑫、彭赫力、王耀、于华、方乃文、郭鹏。T/CWAN0100—2023T/CWAN0100—2023PAGEPAGE10焊接数值模拟固有应变法范围的焊接变形预测。规范性引用文件（包括所有的修改单）适用于本文件。GB/T31054机械产品计算机辅助工程有限元数值计算术语GB/T33582机械产品结构有限元力学分析通用规则ISO/TS18166焊接数值模拟 执行和文档(Numericalweldingsimulation—Executionanddocumentation)ISO/TR25901（所有部分）（Weldingandalliedprocesses0063—2023焊接数值模拟热弹塑性有限元方法术语和定义GB/T31054、GB/T33582、ISO/TR25901、ISO/TS18166和T/CWAN0063—2023界定的以及下列术语和定义适用于本文件。固有应变inherentstrain经历焊接热循环后，焊缝及其附近产生的永久应变。固有变形inherentdeformation固有应变在焊缝横截面上的积分结果。固有应变法inherentstrainmethod一种焊接变形数值模拟方法，其原理是将焊缝及其附近的固有应变作为初始值进行一次弹性分析，从而得到焊接结构的整体变形。基本流程固有应变法焊接数值模拟的基本流程主要包括有限元模型建立、模型求解、结果评估、仿真方案制定实施及报告编写，具体参见附录A。一般要求有限元模型分析前应符合下述要求：GB/T33582BSI单位制，质量单位选择吨，长度单位选择毫米，应力单位选择兆帕，时间单位选择秒，密度单位选择吨/立方毫米，力的单位选择牛，模型材料属性单位应与几何模型单位一致；——焊接CAD模型应结构完整，装配关系清晰明了；焊接仿真分析模型建立规则几何模型建立及简化几何模型简化应符合准确性和经济性两个原则。焊接仿真分析的几何模型简化原则：——应对模型进行干涉检查，避免模型中存在干涉；——应删除各个零件上的圆角、切角、退刀槽、凸台、无焊缝填充的倒角等细节特征；——应删除各零件上孔径与结构件尺寸比值小于1%的孔及安装孔；——应简化不规则的截面形状。有限元网格划分采用固有应变法进行焊接数值模拟的网格模型，可使用壳单元网格划分或实体单元网格划分。对于板厚与长/1:10的结构件，宜采用实体单元网格划分；反之，宜采用壳单元网格划分。基于壳单元的网格划分应遵守以下原则：——应对模型进行干涉检查，避免模型中存在干涉；——将模型导入前处理软件中，对几何模型抽取中面；——应调整需要对齐的平面，并去除中面上的圆倒角、切角等细节特征；——应去除中面间隙，使各中面相交；——应修补模型的破面，合并距离很近的边及节点；——应去除模型中存在的自由边；——应将对接接头、角接接头各组件中面的交线作为焊缝，如图1和图2所示；——应在搭接接头两组件中面之间建立一个面P，并将该面作为搭接接头P1组件中面的一部分，P1组件与P2组件的交线为搭接接头的焊缝，如图3所示；——应在端接接头两组件中面之间建立一个面，并将该面平均分为A，B两个面，分别赋给P1，P2两个组件的中面，两组件中面的交线为端接接头的焊缝，如图4所示；——建立焊接仿真分析模型时，宜使用四边形低阶单元；——应在焊缝位置划分足够多的网格来准确预测变形，小型结构焊缝长度方向宜不少于8个单元，大型结构宜不少于6个单元，超大型结构宜不少于4个单元。标引序号说明：

a）几何模型 b）中面模型1——几何模型中的组件P1；2——几何模型中的焊缝；3——几何模型中的组件P2；4——P1的中面；5——P2的中面；6——中面模型中的组件P1；7——中面模型中的焊缝；8——中面模型中的组件P2。

图1对接接头的简化 标引序号说明：1——2——几何模型中的焊缝；

a）几何模型 b）中面模型3——几何模型中的组件P2；4——P1的中面；5——P2的中面；6——中面模型中的组件P1；7——中面模型中的焊缝；8——中面模型中的组件P2。

图2角接接头的简化标引序号说明：

a）几何模型 b）中面模型1——几何模型中的组件P1；2——几何模型中的焊缝；3——几何模型中的组件P2；4——P1的中面；5——P2的中面；6——中面模型中的组件P1；7——中面模型中的焊缝；8——中面模型中的组件P2;9——中面模型中的面P。

图3搭接接头的简化标引序号说明：

a）几何模型 b）中面模型——P1；——几何模型中的焊缝；——P2；——P1的中面；——P2的中面；——P1；——A；——中面模型中的焊缝；——中面模型中的面B；10——中面模型中的组件P2。

图4端接接头的简化T/CWAN0063—2023的规定。焊接接头固有变形计算（B，主要包括以下内容：a）附录B.1给出了采用实验方法提取固有变形的步骤；b）附录B.2给出了采用热弹塑性法提取固有变形的步骤；c）附录B.3给出了采用经验公式计算固有变形的步骤。计算条件设置力学边界条件的施加模型结构边界条件的设定应符合以下原则：——焊接仿真分析模型的结构边界条件设定应符合实际工况；——当在自由状态下进行焊接变形计算时，应限制仿真模型的6个自由度，以防止工件发生刚体位移和转动，宜使用如图5所示的约束方式；——当在装夹条件下进行焊接变形计算时，应在夹具作用位置施加位移约束，限制仿真模型的6个自由度，以防止工件发生不符合夹具限制的位移和转动，宜使用如图6所示的约束方式。图5限制刚体位移的约束方式图6限制装夹条件下刚体位移的约束方式固有变形的施加基于壳单元模型的固有变形的施加应符合以下原则：——单元法向应保持一致；——固有变形的施加位置应符合6.2.2中规定的焊缝位置。基于实体单元模型的固有应变的施加应符合以下原则：——宜施加到焊缝热影响区；——宜通过修改材料热膨胀系数的方法施加固有应变，应保证固有应变的方向与坐标系的方向一有限元模型的检查网格模型的检查网格划分完成后，应对网格进行检查，检查内容及要求如下：——模型中不应有重复、折叠、畸变和连接问题的网格；1的规定；——应优先保证焊缝位置的网格质量；——应检查相连组件间是否存在未耦合的节点，并进行节点耦合操作。表1 单元质量检查控制参数类别四边形单元三角形单元六面体单元楔形单元四面体单元长宽比≤5.0≤5.0≤5.0≤5.0≤5.0翘曲度≤16°—≤18°≤18°—偏斜度≤60°≤60°≤60°≤60°≤60°内角40°～135°20°～120°40°～135°20°～120°20°～120°雅克比≥0.3≥0.3≥0.3≥0.3≥0.3模型设定的检查模型设定完成后，应对模型进行检查，检查内容及要求如下：——力学边界检查，依据模型在实际焊接过程中的约束状态，检查模型的约束。仿真结果评估表象评估结果表象评估，具体方法如下：——检查模型的奇异比；——检查模型的收敛性；——分析焊接变形的合理性。变形趋势评估通过三维扫描的形式确定物理样件的变形趋势，具体方法如下：——采用三维扫描仪对焊接后的物理样件进行整体扫描；——将有限元分析的网格模型导出为stl文件；——采用最小二乘法进行两个文件拟合，比对模型的重合度。关键位置数据评估——对于对接接头，关键位置可设置为焊缝上表面中心线上的等间距点；——对于角接头，关键位置可设置为角焊缝中心线、角焊缝焊趾线上的等间距点；——对于搭接接头，关键位置可设置为焊缝上下表面中心线上的等间距点；——对于端接接头，关键设置可设置为左右两条焊缝，内外表面中心线上的等间距点。仿真方案制定与实施基于7.1~7.3仿真分析报告根据具体的分析对象、分析目的等编写分析报告，模板见附录C，仿真实施过程可参见附录D示出的对接接头MIG焊接的变形模拟仿真实例，报告至少包括但不限于以下内容：——仿真对象，详细描述组件或者结构、焊接工艺和参数；——仿真目标，介绍分析问题背景，说明采取的分析手段、阐述预期结果细节；——物理模型，描述需要模拟的物理过程，给定的边界条件及简化假设；——分析过程，描述几何模型的简化、网格划分、材料属性、边界条件、求解参数及求解方式；/——模型校核，应给出保证模拟结果正确的措施，主要有模型核查、与物理样件对比验证。附 录 A（规范性）基于固有应变法焊接数值模拟流程图基于固有应变法焊接数值模拟基本流程如图A.1所示。A.1基于固有应变法焊接数值模拟流程图附 录 B（资料性）实验方法

固有应变的提取及数据修正利用实验方法提取对接接头固有变形的主要步骤如下：——根据图B.1及表B.1所示的细节参数截取两块板材，并根据板材厚度填充表B.1第一行数据；——拼装两块平板并点焊固定，使用冲头和锤子在平板上下表面进行标记，位置如图B.1中的黑点所示；——宜采用三坐标仪测量标记点位置坐标，计算表B.1第二行相关数据并填充；——测量后在自由状态下焊接平板；——焊后测量接头相关尺寸，填充表B.1第三行数据；——将相关数据带入公式（B.1）～（B.3）计算接头纵向固有收缩、横向固有收缩和横向固有弯曲，其中ΔL或Δb为负值时表示收缩；——一般情况下纵向固有弯曲（θx*）取值0rad。注1：OL为该对接接头总长度。注2：B为接头总宽度。注3：L、b为参考点位置，其中L参考点仅在下表面标记，b参考点在上、下表面均标记。4：ΔL为焊接后上下表面的平均收缩量。5：Δbtb的收缩量。6：Δbb为焊接后接头下表面b7：e1L参考点的上翘量。8：e2L参考点的上翘量。注9：d1、d2分别为左、右L参考点与焊缝中心的距离。图B.1 对接焊接头的实验测量表B.1 对接接头固有变形提取中的相关尺寸 单位为mm尺寸 OL L B bt bb b=(bt+bb)/2

d1 d2 e1 e2推荐 100t 60t 50t 30t 30t 20t 20t 0 0定位焊后焊后注1：t为板厚。2：bt为b在焊后接头上表面的对应值。3：bbb在焊后接头下表面的对应值。式中：δx*——纵向固有收缩，单位为毫米δy*——横向固有收缩，单位为毫米θy*——横向固有弯曲，单位为弧度（rad）；θy1——接头左侧横向固有弯曲，单位为弧度（rad）；θy2——接头右侧横向固有弯曲，单位为弧度（rad）；Δb——上表面收缩（Δbt）与下表面收缩（Δbb）的平均值。

……（B.1）……………………（B.3）对于由一块平板和腹板组成的T形接头（如图B.2所示），可用类似于对接接头的方式来测量接头的固有变形。当难以准确量取平板上表面b参考点间距时，可只在平板背面做好标记测量距离变化。平（和公式（B.7）～（B.9）计算，通常认为平板和腹板的纵向固有弯曲（θx*）为0rad。注1：H表示腹板高度。注2：ΔH注3：c表示焊接后腹板顶端横向位移。式中：

图B.2 T形接头的实验测

…………（B.4）（B.5）…（B.6）δxp*——平板纵向固有收缩，单位为毫米（mm）；δyp*——平板横向固有收缩，单位为毫米（mm）；……θyp*——平板横向固有弯曲，单位为弧度（rad）；t——平板板厚，单位为毫米（mm）。……式中：δxw*——腹板纵向固有收缩，单位为毫米δyw*——腹板横向固有收缩，单位为毫米腹板横向固有弯曲，单位为弧度（rad）。热弹塑性有限元法测量法

……（B.7）…（B.8）………………（B.9）利用测量法从热弹塑性有限元计算结果中提取固有变形的主要步骤如下：——利用热弹塑性有限元法计算简单焊接接头变形；——参考B.1中提及的方法标记有限元模型中的参考点，并测量相关数据；——根据公式（B.1）～（B.3）计算对接接头固有变形，根据公式（B.4）～（B.9）计算T形接头固有变形，依然认为纵向固有弯曲为0rad。积分法利用积分法从热弹塑性有限元计算结果中提取固有变形的主要步骤如下：——利用热弹塑性有限元法计算简单焊接接头变形；——从计算结果中提取焊缝截面上沿厚度方向（Z向）x方向上的εyεy和X-Yy.3——根据塑性应变分布图，利用公式（B.10）～（B.13）计算接头固有变形：……………（B.11）………（B.12）………（B.13）式中：t——焊接板材厚度，单位为毫米（mm）；z——以板中面为参考的厚度方向坐标。图B.3 塑性应变分量εpx，εpy,εpxy在焊缝横截面上的分布图经验公式方法基于壳单元的经验公式法对接接头固有变形计算经验公式如下：………………（B.14）…………（B.15）…………（B.16）（B.18）0…………（B.19）式中：Q——有效热输入，单位为焦每毫米（J/mm）；——焊接电压，单位为伏特（V）；——焊接电流，单位为安培（A）；v——焊接速度，单位为毫米每秒（mm/s）；n——焊接道次；h——等效板厚，单位为毫米（mm）；t1——对接接头其中一板板厚，单位为毫米（mm）；t2——对接接头另一板板厚，单位为毫米（mm）；δx*——纵向固有收缩，单位为毫米（mm）；——弹性模量，单位为兆帕斯卡δy*——横向固有收缩，单位为毫米（mm）；θy*——横向固有弯曲，单位为弧度（rad）；θx*——纵向固有弯曲，单位为弧度（rad）；t——平板板厚，单位为毫米（mm）。角接接头固有变形计算经验公式如下：

………（B.20）…………（B.21）…………（B.22）…………………（B.23）…………（B.24）…………（B.25）………………（B.26）……（B.27）…………（B.28）式中：Qp——角接接头底板有效热输入，单位为焦每毫米（J/mm）；Qw——角接接头腹板有效热输入，单位为焦每毫米t1——腹板板厚，单位为毫米（mm）；t2——底板板厚，单位为毫米（mm）；——室温时引起工件纵向焊接变形的纵向收缩力，单位为牛顿底板纵向固有收缩，单位为毫米底板横向固有收缩，单位为毫米θyp*——底板横向固有弯曲，单位为弧度（rad）；θxp*——底板纵向固有弯曲，单位为弧度（rad）；δxw*——腹板纵向固有收缩，单位为毫米δyw*——腹板横向固有收缩，单位为毫米腹板横向固有弯曲，单位为弧度（rad）；腹板纵向固有弯曲，单位为弧度（rad）。基于实体单元的经验公式法基于实体单元的固有应变计算经验公式如下：

…………（B.29）…………（B.30）…………（B.31）…………（B.32）…………（B.33）…………（B.34）…………（B.35）…………（B.36）………（B.37）………（B.38）式中：Q——焊接线能量，单位为焦每毫米（J/mm）；——焊接电压，单位为伏特（V）；——焊接电流，单位为安培（A）；v——焊接速度，单位为毫米每秒（mm/s）；Wx——焊缝纵向固有应变的总和，单位为平方毫米（mm2）;Wy——焊缝横向固有应变的总和，单位为F——焊缝横截面积，单位为平方毫米（mm2）；εx——纵向固有应变；εy——横向固有应变；K——纵向固有应变系数，单位为立方毫米每焦（mm3/J）；ξ——横向固有应变系数，单位为立方毫米每焦（mm3/J）；h——板材厚度，单位为毫米（mm）。对接接头固有变形修正.4（.3～（B.48)计算四个固有变形参数（纵向固有收缩，横向固有收缩，纵向固有弯曲，横向固有弯曲），记为（a0，b0，c0，d0）。选择的测量点应遵守以下原则：——对于横向收缩和横向弯曲，选择中间横截面上的点；——对于纵向收缩，选择位于边缘和焊缝之间的点；——对于纵向弯曲，选择靠近焊缝的点；注：Pi（xi，yi，zi）表示板（壳模型）平面上的任意点。图B.4 对接接头实验测量或FEM取点示意图纵向固有收缩a0可以通过焊缝长度方向的平均收缩评价,较大的纵向弯曲使两个测量点之间的直线距离比实际收缩短，故应使用曲线距离：式中：LC——两个测量点间的曲线距离，单位为毫米（mm）；B——对接接头宽度,单位为毫米（mm）；dz——两个测量点间的曲线对应弦长，单位为毫米（mm）；x0,i——点Pi的初始横坐标，单位为毫米（mm）；y0,i——点Pi的初始纵坐标，单位为毫米（mm）；xi——点Pi的变形后的横坐标，单位为毫米（mm）；yi——点Pi的变形后的纵坐标，单位为毫米（mm）；Uzi——Piz方向上的位移，单位为毫米（mm）

………………（B.39）……（B.40）……（B.41）……（B.42）……（B.43）……………（B.44）……………（B.45）b0，可用三个垂直于焊缝方向的截面上的平均收缩率近似表示：…………（B.46）c0，可用焊缝处或焊缝附近两条纵线处的平均弯曲曲率表示：…………（B.47）d0，可用焊缝处或焊缝附近两条纵线处的弯曲角表示,弯曲角单位应为弧度（rad）：………（B.48）当仿真结果与实际变形差距较大时，应对接头固有变形参数进行修正，具体步骤如下：——提取焊接试板上各标记点的实际焊接变形(Ux0(Pi)，Uy0(Pi)，Uz0(Pi))B.4.2～B.4.5计算初始固有变形参数(a0，b0，c0，d0)；——将初始固有变形参数(a0，b0，c0，d0)JWRIAN软件中进行焊接变形预测，提取预测结果中各标记点变形(Ux1(Pi)，Uy1(Pi)，Uz1(Pi))B.4.2～B.4.5计算固有变形参数c1，d1)ratio_a1=a1/a0，ratio_b1=b1/b0，ratio_c1=c1/c0，ratio_d1=d1/d0评估准确性，1.0，说明初始固有变形参数(a0，b0，c0，d0)较为准确，可将该参数作为实ERa1=a1-a0，ERb1=b1-b0，ERc1=c1-c0，ERd1=d1-d0；a1=a0/ratio_a1，b1=b0/ratio_b1，c1=c0/ratio_c1，d1=d0/ratio_d1更新固有变形参数JWRIAN中输入更新后的B.4.2～B.4.5计算固有变形参数ERc2=c2-c0，ERd2=d2-d0，如果误差足够小，则将(a1，b1，c1，d1)作为实际接头的固有变形；(a1，b1，c1，d1)和(a2，b2，c2，d2)B.5JWRIAN并计算其误差，如果误差足够小，则将其作为实际接头固有变形。图B.5 固有变形参数线性估计T/CWAN0100—2023附 录 CT/CWAN0100—2023附 录 C（资料性）PAGEPAGE16表C.1给出了仿真报告文档模板。

文档模板表C.1 文档模板编制：XXXX公司名称：XXXX部门：XXXX依据：焊接数值模拟-固有应变法版本：V1日期：JJJ-MM-TT第X页共X页封面简要说明仿真对象仿真目标物理和数学模型求解方法和应用软件结果和结论摘要确保模拟结果质量采取的措施摘要措施备注/说明已检查[]是[]否模型检查已验证[]是[]否仿真结果与实验结果对比，例如:激光三维扫描变形测量其他：备注（可选）T/CWAN0100—2023T/CWAN0100—2023PAGEPAGE17附 录 D（资料性）对接接头MIG焊接的变形模拟仿真报告表D.1给出了对接接头MIG焊接的变形模拟仿真报告。表D.1 对接接头MIG焊接的变形模拟仿真报告公司名称：XXXX部门：XXXX依据：焊接数值模拟固有应变法编制：XXXX版本：V1日期：JJJ-MM-TT1页共4页封面简要说明仿真对象：201不锈钢对接接头的MIG脉冲电弧焊，板厚4mm，间隙为2.0mm的I型坡口。仿真目标：模拟201不锈钢MIG焊接变形和残余应力。物理和数学模型：/求解方法和应用软件：固有应变法焊接变形仿真，Jwrian结果和结论摘要：计算出的变形与光学变形测量变形趋势一致，数值相近。确保模拟结果质量采取的措施摘要措施备注/说明已检查[]是[]否模型检查已验证[]是[]否仿真结果与实验结果对比，激光三维扫描变形测量。其他：备注（可选）公司名称：XXXX部门：XXXX依据：焊接数值模拟固有应变法编制：XXXX版本：V1日期：JJJ-M