TCWAN 0063-2023 焊接数值模拟热弹塑性有限元方法

ICS25.160.01CCSJ33团 体 标 准T/CWAN0063—2023焊接数值模拟热弹塑性有限元方法Numericalweldingsimulation-Thermalelastic-plasticfiniteelementmethod2023-12-27发布 2024-02-01实施中国焊接协会 发布T/CWAN0063—2023T/CWAN0063—2023PAGE\*ROMANPAGE\*ROMANI目 次前 言 II范围 1规范性引用文件 1术语和定义 1基本流程 2一般要求 2焊接仿真分析模型建立规则 2求解器设置 7仿真结果评估 7仿真分析报告 8附录A（规范性）基于热弹塑性法焊接数值模拟流程图 9附录B（资料性）文档模板 10附录C（资料性）对接接头MAG焊接的变形模拟仿真报告 11前 言GB/T1.1—20201本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中国焊接协会提出并归口。本文件起草单位：江苏徐工工程机械研究院有限公司、郑州机械研究所有限公司、哈尔滨职业技术学院、河南科技大学、中国机械总院集团哈尔滨焊接研究所有限公司、中信重工机械股份有限公司、青岛海尔空调器有限总公司、哈尔滨电机厂有限责任公司、上海航天精密机械研究所、上海工程技术大学、北京星航机电装备有限公司、南昌航空大学。本文件主要起草人：张立平、蹤雪梅、龙伟民、张贵芝、郝双双、魏世忠、孟政宇、滕彬、武汉T/CWAN0063—2023T/CWAN0063—2023PAGEPAGE10焊接数值模拟热弹塑性有限元方法范围本文件适用于基于热弹塑性有限元方法的焊接结构件熔化焊焊接变形及焊接残余应力的焊接数值模拟。规范性引用文件（包括所有的修改单）适用于本文件。GB/T3651金属高温导热系数测量方法GB/T43384339金属材料热膨胀特征参数的测定GB/T22315金属材料弹性模量和泊松比试验方法GB/T25846工业用γ射线密度计GB/T31054机械产品计算机辅助工程有限元数值计算术语GB/T33582机械产品结构有限元力学分析通用规则NB/SH/T0632比热容的测定差示扫描量热法0009焊接术语熔化焊ISO/TS18166焊接数值模拟 执行和文档(Numericalweldingsimulation-Executionanddocumentation)术语和定义

GB/T31054、GB/T33582、T/CWAN0009和ISO/TS18166界定的以及下列术语和定义适用于本文热源模型heatsourcemodel熔化焊接热能特征及其与工件作用后在工件上分布的数学表示。热弹塑性有限元方法thermalelastic-plasticfiniteelementmethod双椭球热源模型doubleellipsoidalheatsourcemodel焊接热源在工件上呈前后两个半椭球分布的数学表示。生死单元birth-deathelement一种模拟焊缝填充的方法，其原理是在焊道未填充之前，将代表焊道的单元材料性能乘以一个较小的比例系数，当焊道填充之后，单元材料属性恢复到真实值。基本流程热弹塑性法焊接数值模拟的基本流程主要包括有限元模型建立、模型求解、结果评估、仿真方案制定实施及报告编写，具体参见附录A。一般要求有限元模型分析前应符合下述要求：GB/T33582B——焊接结构件CAD模型应结构完整，装配关系清晰明了；焊接仿真分析模型建立规则几何模型简化几何模型简化应符合准确性和经济性两个原则。准确性是指应保留分析对象主要特征，以保证受力特性不变。经济性是指应用较小分析成本取得足够精度的计算结果。焊接仿真分析的几何模型简化原则：——应对模型进行干涉检查，避免模型中存在干涉；——应忽略对焊接变形影响较小且不重要的零件，如吊耳、安装板、固定板等；——应删除各个零件上的圆角、退刀槽、凸台、无焊缝填充的倒角等细节特征；——应删除各零件上孔径与结构件尺寸比值小于1%的孔及安装孔；——应简化不规则的截面形状；——在钝边不大于2mm或者要求单面焊双面成型的情况下，带钝边倒角应修改成尖角，如图1所示；——将模型导入前处理软件中，应消除模型的自由边，并检查模型的干涉情况；——应修补模型的破面，合并距离很近的边及节点。1钝边的简化有限元网格划分建立焊接仿真分析模型时宜选择三维六面体单元。建立焊接仿真分析模型时宜选用低阶单元。2mm3m～421:31:2的网格过渡形式，降低远离焊缝位置的网格数量。3ⅠⅡa1取9mm)，网格可多次过渡。图2网格过渡形式图标引序号说明：①——坡口区域；Ⅰ——热影响区；Ⅱ——网格过渡区；Ⅲ——母材区；a1——热影响区宽度；a2——过渡区宽度；a3——母材宽度；b——板厚。图3模型分块情况图2层以上的实体单元网格。5。材料模型建立用于焊接仿真分析的材料热物理性能参数应包括热导率、比热容、密度、线膨胀系数、弹性模是关于温度的函数。Von-Mises屈服准则，材料进入塑性区后服从塑性流动准则和等向强化准则。材料低温区间(0℃～600℃)1规定的方法和设备进行测量。材料高温间(≥600℃)的热物理性能参数可通过数据外插拟合获得，也可通过具备材料性能计算功能的软件获得。材料室温屈服强度≥890MPa的情况下，应考虑相变/热焓转变对焊接变形及残余应力的影响。参数名 测试方法参考标准测试设备参数名 测试方法参考标准测试设备热导率 直接通电纵向热流法GB/T3651管式防热炉、热电偶密度 射线吸收量测定法GB/T25846放射性同位素密度计线膨胀系数 接触式示差法GB/T4339膨胀仪比热容/热焓转变 差示扫描量热法NB/SH/T0632差示扫描量热仪弹性模量 动态法GB/T22315共振测量装置屈服强度 试样拉伸法GB/T4338高温拉伸试验机W/(m•K)Kg/m31/KJ/(kg·C)N/m2N/m2热源参数定义电弧焊宜选用双椭球体/双椭圆面热源模型，激光焊宜选用高斯旋转体热源/高斯面热源模型，激光-电弧复合焊宜选用高斯旋转体热源/高斯面热源与双椭球/双椭圆热源的叠加。1e-3～1e-5。应采用试验与仿真对比的方法确定焊接热源的参数，包括以下两种方法。通过对比热循环曲线确定焊接热源参数：5mm处的热循环曲线；4a），测试该位置的热循环曲线；4b)，反复调整10%以内。热电偶放置位置 b)热循环曲线对比图4热循环曲线法通过对比焊缝形貌确定焊接热源参数：在焊接试板上截取试样，打磨、抛光、腐蚀，获得焊接接头的宏观金相；2倍熔深，进行焊接热分析，提取焊接仿真分析的热源形貌；5，判断其吻合性，反复调整焊接热源参数，直至吻合良好，完成热源校核。图5焊接接头宏观金相对比图计算条件设置模型热边界条件设定模型热边界条件的设定应符合以下原则：——根据焊接操作指导文件上的焊接方向、焊接顺序，及确定的焊接热源参数设定焊接热边界条件；5W/(m2•K)～25W/(m2•K)，也可根据牛顿定律测定对流换热系数。模型结构边界条件设定模型结构边界条件的设定应符合以下原则：——焊接仿真分析模型的结构边界条件设定应符合实际工况；366；6自由焊接时模型约束图——结构件在工装约束状态下焊接时，可采用刚体接触或者固定约束的方法模拟焊接工装的装夹，并在工装撤去后，释放接触/约束，并施加辅助约束；——可采用接触的方法考虑工艺拉筋对焊接变形的影响；模型边界条件的有效简化如计算机算力不足或计算规模较大，应对焊接边界条件进行简化，简化方法如下：——焊道集中方法，对于多层多道焊，可将焊道合并，在等效热效应前提下，一次输入焊接能量，设定焊接热源，保证焊接区域内单位体积的能量相对于简化前不变；——分段移动热源方法，将一条焊道分为多段，每段采用热循环曲线或者带状热源进行加载；——热机解耦计算方法，将多步的热分析计算结果用于一步的焊接变形及应力计算分析。有限元模型的检查网格模型的检查网格划分完成后，应对网格进行检查，检查内容及要求如下：——模型中不应有重复、连接问题、折叠及畸变网格；2的规定；——优先保证焊接及热影响区位置的网格质量。2单元质量检查控制参数参数长宽比雅克比内角翘曲度偏斜度数值5～10≥0.540°～135°≤18°≤60°模型设定的检查模型设定完成后，应对模型进行检查，检查内容及要求如下：——材料参数检查，根据焊接仿真分析模型设定的材料热物理性能参数，复核材料各个参数，检查材料各热物理性能参数曲线平滑性，无突变点，单位应统一；——力学边界检查，依据模型在实际焊接过程中的约束状态，检查模型的约束；求解器设置有限元分析应合理设置分析环境、时间参数、收敛准则、存储步数、需求分析结果，应符合以下要求：——焊接热力模拟可采用直接耦合，也可采用间接耦合。推荐计算效率高的间接耦合方式；0.1s～0.5s60%焊接熔池重合。焊后冷却阶段，温度梯度较小，5s～10s；——板材厚度不大于5mm的结构件进行焊接仿真分析时，应采用大应变模型；——应采用集中质量热容矩阵克服瞬态温度场分析中的阶跃现象；——应采用并行计算的方式进行焊接仿真分析，每个CPU计算的网格数应控制在1.5万～3万；——应输出温度、应力和位移结果三种场变量。仿真结果评估表象评估结果表象评估，具体方法如下：——检查模型的奇异比；——检查模型的收敛性；——分析焊接变形及焊接残余应力的合理性。变形趋势评估通过三维扫描的形式确定物理样件的变形趋势，具体方法如下：——采用三维扫描仪对焊接后的物理样件进行整体扫描；——将有限元分析的网格模型导出为stl格式的文件；——采用最小二乘法进行两个文件拟合，比对模型的重合度。关键位置数据评估进行物理样件试验，对比关键位置的有限元结果和试验结果，以试验结果为依据修正有限元模型，重新计算评估，保证变形结果偏差不大于20%。仿真方案制定与实施基于8.1~8.3仿真分析报告根据具体的分析对象、分析目的等编写分析报告，模板见附录B，仿真实施过程可参见附录C示出的对接接头MAG焊接的变形模拟仿真实例，报告至少包括但不限于以下内容：——仿真对象，详细描述组件或者结构、焊接工艺和参数；——仿真目标，介绍分析问题背景，说明采取的分析手段、阐述预期结果细节；——物理模型，描述需要模拟的物理过程，给定的边界条件及简化假设；——分析过程，描述几何模型的简化、网格划分、材料模型、边界条件、求解参数及求解方式；——结果分析，应给出关键位置的焊接变形值及典型的图标结果，如变形、应力/应变云图。图表应简明易懂且不应有无关的信息；——模型校核，应给出保证模拟结果正确的措施，主要有模型核查、与物理样件对比验证。附 录 A（规范性）基于热弹塑性法焊接数值模拟流程图基于热弹塑性法焊接数值模拟基本流程如图A.1所示。A.1基于热弹塑性法焊接数值模拟流程图B.1所示。

附 录 B（资料性）文档模板表B.1文档模板公司名称：部门：依据：焊接数值模拟-热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT第X页共X页封面简要说明仿真对象仿真目标物理和数学模型求解方法和应用软件结果和结论摘要确保模拟结果质量采取的措施摘要措施备注/说明已校核[]是[]否热源校核已检查[]是[]否模型检查已验证[]是[]否仿真结果与实验结果对比，例如:激光三维扫描变形测量、盲孔法残余应力测试其他：备注（可选）T/CWAN0063—2023T/CWAN0063—2023PAGEPAGE11附 录 C（资料性）对接接头MAG焊接的变形模拟仿真实例对接接头MAG焊接变形模拟仿真报告如表C.1所示。表C.1对接接头MAG焊接的变形模拟仿真报告公司名称：XXXX部门：XXXX依据：焊接数值模拟热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT1页共7页封面简要说明仿真对象：Q345低合金钢对接接头的MAG脉冲电弧焊，板厚6mm，30°V型坡口。仿真目标：模拟Q345钢MAG焊接变形和残余应力。物理和数学模型：热-力间接耦合模型；热源模型—具有容积热源的热控制模型（Goldak）；材料塑性本构模型—具有与应变速率有关的各向同性硬化的弹塑性模型。求解方法和应用软件：热-机联合仿真，Sysweld结果和结论摘要：仿真模拟获得的温度场与试验确定的温度场一致，该温度场由热循环测试校核；模拟获得的横向、纵向残余应力与相应的X射线测量结果十分吻合；计算的变形与光学变形测量变形趋势一致，数值相近。确保模拟结果质量采取的措施摘要措施备注/说明已校核[]是[]否热源校核已检查[]是[]否模型检查已验证[]是[]否X射线残余应力测试。其他：备注（可选）公司名称：XXXX部门：XXXX依据：焊接数值模拟热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT第2页共7页模拟对象工件尺寸：300mm×120mm×6mm母材：Q345低合金钢 填充材料：G3Si1保护气体：M21 保护气体流量：25L/min工艺：半自动化MAG焊接，135 电源：逆变直流电焊机PHOENIX521（EWM）电弧类型：脉冲电弧 平均功率：6783.34W焊接速度：9.167mm/s 焊丝直径：1.2mm干伸长：13mm 焊枪角度：0°焊缝开口角度：30°(2×15°) 接头形式，拼焊位置：V型对接，0mm，150mm，300mm预热温度：20°C 夹紧：静态确定单位为毫米标引序号说明：1——热电偶；2——陶瓷垫板；3——装配平台；4——拼焊点。图C.1工件示意图模拟的目的考虑材料相变行为，构建焊接模拟模型，进行焊接变形和残余应力分析，并进行合理实验验证。公司名称：XXXX部门：XXXX依据：焊接数值模拟热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT第3页共7页物理模型做以下几点假设：——热-力”向“热”——假设部件没有预应力存在；——根据焊缝几何形状生成三维模型，根据焊接要求划分网格；——忽略组件尺寸公差；——3D网格模型中填充焊缝单元，计算过程采用生死单元技术；——采用三点约束施加力学边界条件。数学模型和求解方法热源模型(.1.2)(.3。6(C.1)q 1 ex 1 1 1 6 (C.2)q 2 e 2 2 2 12 (.3)式中：a1,a2,b,c——双椭球形体热源形状参数Φ ——有效功率；x ——与热源中心x方向上的距离y ——与热源中心y方向上的距离z ——与热源中心z方向上的距离；f1，f2——前半部分、后半部分椭球的能量分配系数，f1+f2=2；q1，q2——前半部分、后半部分椭球的热流密度。公司名称：XXXX部门：XXXX依据：焊接数值模拟热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT第4页共7页图C.2双椭球形体热源模型C.4.2材料塑性本构模型具有与应变速率有关的各向同性硬化的弹塑性模型。对于相变塑性，假设相变塑性应变εtp与现有应力σ相对于单轴应力状态之间是线性关系，因此以下适用：(C.4)式中：K——转换塑性常数，单位MPa-1。实施模型网格划分和约束设置考虑到焊缝形状和部件几何形状，使用六面体单元进行3D网格划分。焊接方向上的单元长为1mm（见图C.3）。标引序号说明：1-焊接方向。a焊缝金属单元大小：ex=1mm，ey=ez=0.5mm。图C.3试板网格和约束热源功率热源功率Q/WC1mmC2mmabmmx0mmy0mmz0mmMA110603.06.04.52.00.00.00.0-3公司名称：XXXX部门：XXXX依据：焊接数值模拟热弹塑性有限元方法编制：XXXX版本：V1日期：JJJ-MM-TT第5页共7页表C.1给出了工况和时间步的设置。表C.1计算时间步长概述工况 热分析 机械分加热t∈[0s;32.7s] 固定:dt=0.109s 固定:dt=0.0545s冷却t∈[32.918s;100s] 固定:dt=1.0s 固定:dt=1.0s冷却t∈[201s;2000s] 自适应:dtMax=50s 自适应:dtMax=50sC.5.2热源描述热输入通过双椭球热源模型实现。焊接轨迹由作为焊接开始的点P1（0,60,0）和作为焊接结束的P2（300,60,0）限定。热源沿轨迹以vw=9.167mm/s的速度均匀地移动。态高温范围受热源的影响。双椭球热源的形状参数列于表C.2中。参数X0，Y0，Z0描述热源到假定的实际弧的基点的几何P1和P21.2倍。表C.2双椭球热源参数