波形钢腹板梁焊接残余应力分布及试验验证.pdf

文章编号：0258-2724(2024)02-0289-09DOI:10.3969/j.issn.0258-2724.20220054波形钢腹板梁焊接残余应力分布及试验验证冀伟1,2，刘勇2（1.湖南大学土木工程学院，湖南长沙410082；2.兰州交通大学土木工程学院，甘肃兰州730070）摘要：为研究波形钢腹板梁焊接残余应力的分布规律，通过有限元软件建立其三维热弹塑性模型，利用热-力耦合分析技术对其焊接温度场和应力场进行有限元数值仿真，采用双椭球体热源和修改单元材料属性的方法实现能量输入和焊缝填充，并将模拟计算结果与实测值进行对比分析.结果表明：有限元预测的波形钢腹板梁残余应力分布与实测结果具有相同趋势，在波形钢腹板梁焊缝的弯折角处，残余应力发生一定幅度的连续应力波动；底板和腹板的残余应力峰值均出现在焊缝中心区域，其值分别为材料屈服强度的 1.30 倍和 1.26 倍；底板纵向残余拉应力在焊缝中心线两侧 78mm 范围内急速下降后缓慢过渡为压应力，在底板较窄一侧压应力线性增大，最大值约为材料屈服强度的 0.61 倍；在底板较宽一侧压应力线性减小，并在边缘处转化为拉应力；焊接速度对残余应力分布扰动不大，而对残余应力峰值影响较显著；当焊接速度从 150mm/min 增加至 250mm/min 时，横向和纵向的残余应力最大值分别增加了 27.11%和 5.88%.关键词：波形钢腹板梁；焊接残余应力；热力学；数值模拟；焊接速度中图分类号：U441.5文献标志码：AWelding Residual Stress Distribution and Experimental Verification ofCorrugated Steel Web GirdersJI Wei1,2,LIU Yong2(1.CollegeofCivilEngineering,HunanUniversity,Changsha410082,China;2.CollegeofCivilEngineering,LanzhouJiao-tongUniversity,Lanzhou730070,China)Abstract:To study the distribution law of welding residual stress of corrugated steel web girders,a three-dimensionalthermo-elastoplasticmodelwasestablishedbyfiniteelementsoftware,andtheweldingtemperaturefieldandstressfieldwerenumericallysimulatedbythefiniteelementmethodusingthermalandmechanicalcouplinganalysistechnology.Doubleellipsoidalheatsourceandmodifiedelementmaterialpropertieswereusedforenergyinputandweldfilling,respectively.Thesimulationresultswerecomparedwiththemeasuredvalues.The results show that the residual stress distribution of corrugated steel web girders predicted by the finiteelementmethodhasthesametrendasthemeasuredvalues.Atthebendingangleoftheweldofthecorrugatedsteelwebgirder,theresidualstressfluctuatescontinuouslyinacertainrange.Thepeakvalueoftheresidualstressinthebottomplateandwebappearsinthecentralareaoftheweld,whichis1.30timesand1.26timestheyieldstrengthofthematerial.respectively.Thelongitudinalresidualtensilestressofthebottomplatedecreasesrapidlywithin78mmonbothsidesoftheweldcenterlineandthenslowlytransitionstocompressivestress.Inaddition,thecompressivestressonthenarrowsideofthebottomplateincreaseslinearly,andthemaximumvalueisabout0.61timestheyieldstrengthofthematerial.Thecompressivestressdecreaseslinearlyonthewidersideofthebottomplateandisconvertedtothetensilestressattheedge.Theanalysisindicatesthattheweldingspeedhaslittleeffectontheresidualstressdistributionbuthasasignificanteffectonthepeakvalueofresidualstress.When收稿日期：2022-01-18修回日期：2022-04-18网络首发日期：2022-05-20基金项目：国家自然科学基金（52168019，51868039）；甘肃省重点研发计划（23YFGA0043）；中央高校基本科研项目（541109030099）第一作者：冀伟（1982），男，教授，博士，研究方向为组合箱梁桥设计理论，E-mail：引文格式：冀伟，刘勇.波形钢腹板梁焊接残余应力分布及试验验证J.西南交通大学学报，2024，59(2)：289-297JIWei,LIUYong.WeldingresidualstressdistributionandexperimentalverificationofcorrugatedsteelwebgirdersJ.JournalofSouthwestJiaotongUniversity,2024，59(2)：289-297第59卷第2期西南交通大学学报Vol.59No.22024年4月JOURNALOFSOUTHWESTJIAOTONGUNIVERSITYApr.2024theweldingspeedincreasesfrom150mm/minto250mm/min,themaximumresidualstressintransverseandlongitudinaldirectionsincreasesby27.11%and5.88%,respectively.Key words:corrugatedsteelwebgirder;weldingresidualstress;thermodynamics;numericalsimulation;weldingspeed波形钢腹板梁在制造过程中，顶底板和腹板通过焊接连接，期间产生复杂的焊接残余应力，对结构的机械性能造成一定不利影响.因此，分析波形钢腹板梁焊接残余应力大小及分布对其力学性能的研究有重要意义1.目前，国内外学者已对直腹板工字钢梁的焊接残余应力展开了研究2-5，而关于波形钢腹板工字钢梁焊接残余应力的研究文献相对较少.Kubo 等6通过分割法实测了工字梁和波形钢腹板梁焊接残余应力大小及分布，并给出了两者的焊接残余应力分布模型，但试验方法过程较复杂，且经济性与可操作性不高.随着计算机软硬件技术的快速发展，有限元软件已成为当前国内外众多学者研究焊接残余应力的重要手段.Kollr 等7提出一种模拟波形钢腹板梁热切割和焊接的有限元模型，通过有限元结果与实测结果的对比验证了有限元模型的正确性.但其采用的生死单元加载热源方法并不能准确地模拟出实际焊接中不同位置的热弧现象8.以上研究的对象均为缩尺比例的波形钢腹板梁，腹板和底板的厚度均有所减小，其残余应力的大小及分布模式与真实值存在差异9.此外，实际制造中，波形钢腹板梁所用钢板较厚，目前的试验方法很难完全测量焊接残余应力，数值模拟则是一种有效的解决途径10.本文采用单向顺序耦合法对波形钢腹板梁各组成部分的温度和残余应力分布进行数值模拟，通过试验结果验证数值模拟结果的正确性，为考虑焊接残余应力的波形钢腹板梁的研究与设计提供参考依据.1 数值模拟分析方法本文采用的 SYSWELD 焊接有限元软件综合考虑了温度对相变和金属相变潜热对温度等多物理场的耦合，其自带的多种焊接模块和材料库可以满足不同焊接材料多种焊接方式的仿真要求.由于焊接模拟时，应力对温度场计算的影响很弱，面对耗时费力的计算过程，采用计算效率高且对结果影响小的顺序耦合法是一个合理的解决方案11.所谓顺序耦合法，即先根据给定的焊接和热边界条件进行非线性瞬态热分析，获得整个焊接过程的温度场，然后，将各载荷步的温度历程作为热荷载读入至模型的应力应变场计算中，完成热-应力耦合计算.焊接模拟具有计算量大、耗费时间长等特点.在进行模拟分析时需对模型进行必要的简化，弱化或忽略对热-力耦合计算结果影响不大的因素，故本文作出以下假设：1）材料服从 VonMises 屈服准则；2）不考虑结构的初始变形和焊接过程中产生的缺陷；3）焊接过程中，焊接参数趋于恒定；4）材料热物性随温度变化，忽略熔池金属的蒸发；5）视焊缝填充材料与母材具有相同的热物性质参数；6）根据结构对称性，取波形钢腹板标准梁的其中一个波段进行建模分析.2 波形钢腹板梁有限元模型 2.1 网格划分以甘肃博睿重装制造的波形钢腹板梁为工程背景，选取一个整波进行建模分析.波形钢腹板梁采用 Q345qDNH 耐候钢，腹板由钢板冷弯轧制而成，波型为 1200 型，高度为 1144mm，厚度为 16mm，底板由 500mm28mm 的矩形平钢板制作而成.波形钢腹板梁的几何构造如图 1 所示.1652703501505003301652701 200R2625001 1441 200281628R90(a)平面图(b)剖面图图1波形钢腹板梁的几何尺寸Fig.1Geometricdimensionsofcorrugatedsteelwebgirder图 2 为波形钢腹板梁的网格划分示意，选用8 节点六面体实体单元对模型进行离散化.在应力变化较大的焊缝区域内网格划分较密集，网格尺寸约为 3mm；远离焊缝处网格划分较稀疏，最大网格290西南交通大学学报第59卷尺寸约为 20mm，有限元模型共包含 199450 个节点、239015 个单元.5001 2001 156腹板底板zxyO(a)总体网格划分腹板内侧底板上侧底板下侧腹板外侧焊缝网格稀疏区网格过渡区网格密集区(b)焊接接头细化(c)底板细化Oxy图2有限元模型网格划分Fig.2Meshingoffiniteelementmodel 2.2 材料的热-力特征高温热源作用于焊缝位置处，材料的热学参数和力学参数均会随温度的变化而发生改变.本文Q345qDNH 耐候钢热物理属性采用 SYSWELD 软件公共材料库里与其牌号相对应的 S355J2G3 钢材进行计算，相关属性如图 3、4 所示12-14.03006009001 2001 50000.20.40.60.81.01.21.4热物理性能参数温度/导热系数/(101 J(mms)1)比热/(102 J(g)1)密度/(102 gmm3)图3S355J2G3 钢热物理属性Fig.3ThermophysicalpropertiesofS355J2G3steel本文采用修改单元材料属性的方法来间接模拟焊缝填充的过程，即假定未焊接部分的材料参数为一较小且合理的值，模拟其为“空相”，待热源作用于相应节点时，其单元材料属性被重新激活为设定值，完成材料的“填充”.相比“生死单元法”，其优点是不会因为刚度矩阵发生突变，而导致求解过程中出现收敛困难的情况.03006009001 2001 500050100150200250300350400力学性能参数温度/杨氏模量/GPa泊松比/102屈服强度/MPa热膨胀系数/(107 1)图4S355J2G3 钢力学属性Fig.4MechanicalpropertiesofS355J2G3steel 3 焊接温度场模拟与分析 3.1 热源模拟frab为更好反映焊接热源的能量分布，并解决平面热源无法考虑熔深的问题，本次模拟采用 Goldak提出的双椭球热源模型15，其几何形状如图 5 所示.图中：、分别为热源前、后两椭球的熔池长度，为热源的熔深，为热源模型的 1/2 融宽.arafzObxay图5双椭球热源模型Fig.5Doubleellipsoidalheatsourcemodel沿着焊接方向，双椭球热源分为前、后两部分，其能量分布满足：qf(x，y，z)=63ffQabfexp3(x2a2f+y2b2+z2a2),（1）qr(x，y，z)=63frQabrexp3(x2a2r+y2b2+z2a2),（2）Q=UIUIfffr式中：为有效输入功率，为焊接电压，为焊接电流，为热效率系数；、分别为热源前、第2期冀伟，等：波形钢腹板梁焊接残余应力分布及试验验证291ff+fr=2后两椭球内任一点处的热流密度分数，.在移动热源作用下，焊件整体温度传导满足式（3）.Tt=QVc+c(2Tx2+2Ty2+2Tz2),（3）TctQV式中：、分别为温度、密度、比热容、时间、热传导系数、内热源.波形钢腹板梁焊接加工属于典型的 T 型接头焊接问题，本文中构件的初始温度取 28，采用 CO2气体保护焊，焊接速度为 200mm/min，焊接电压为3545V，电流为 350450A，热效率为 0.85.3.2 边界条件在进行温度场计算时，主要能量损失来自于试件表面与周围介质的热交换.焊接数值模拟中，通过牛顿冷却公式和 Stefan-Boltzmann 定律考虑焊件表面的热交换和热辐射，其数学表达式为qc=hc(TsTe),（4）qr=(T4sT4e),（5）qcqrTsTehcStefan-Boltzmann式中：、分别为焊接表面的热交换和热辐射，为试件表面温度，为周围环境温度，和为对流和辐射的换热系数，为常数.为简化计算过程，将对流换热和辐射传热整合为等效热对流换热13，如式（6）所示.qt=qc+qr=heq(TsTe)，（6）heq式中：为综合传热系数.3.3 温度场分布当电弧作用于焊缝处，相应位置温度迅速达到材料的熔点并形成熔池.随着热源的移动，熔池前段汲取能量不断熔化，熔池后端释放能量相继冷却.模型各点温度随时间不断发生变化，腹板和底板沿焊缝方向均呈现明显的温度梯度.图 6（a）（c）分别为波形钢腹板梁在焊接过程中 3 个代表时刻的温度场分布图，即起弧、焊接中和收弧阶段；图 6（d）为第 10000.0s 时的温度场分布.(a)t=3.7 s261503004005006507509001 0001 1501 2501 4001 500261503004005006507509001 0001 1501 2501 4001 500261503004005006507509001 0001 1501 2501 4001 50027.9928.0628.1228.1828.2528.3228.3828.4428.5128.5828.6428.7028.77T/T/T/T/(b)t=195.0 s(c)t=394.0 s(d)t=10 000.0 s图6温度场变化云图Fig.6Changenephogramoftemperaturefield从图 6（a）可以看出：在 3.7s 时，焊缝处的温度达到材料的固液相线，焊缝处金属熔化形成熔池；随着时间的推移，焊件温度场进入准稳态，并持续至焊接结束；在熔池前段，温度场梯度线间距较小，在熔池后端，温度场梯度线间距较宽，整个焊接过程持续了 394.0s；冷却时，焊件的整体温度缓慢下降.从图 6（d）可以看出：当冷却 9606.0s 后，焊件整体温度趋于 28，可认为此时焊件已经完全冷却.图 7 为波形钢腹板梁 1/2 截面处，距离焊缝 0、10、20、30、40、50mm 各点的温度历程曲线.从图 7中可以看出：当电弧作用在这一截面时，各个节点的温度迅速上升，焊缝处温度超过材料的熔点，最高温度达到 1548.0；当电弧离开后，各节点温度开始下降，向常温方向发展，距离焊缝中心越近，其温度峰值越高，冷却速率也越快.总体看来，各节点的升温速率明显高于冷却速率，焊件的冷却速率随时间推移而愈发缓慢，有限元的温度场分析结果与实际温度变化过程相符.04008001 2001 60004008001 2001 600T/t/s 0 30 mm 10 mm 40 mm 20 mm 50 mm 图7距离焊缝不同距离的温度历程曲线Fig.7Temperaturechangecurvesatdifferentdistancesfromtheweld292西南交通大学学报第59卷 4 焊接残余应力场模拟与分析 4.1 边界条件在进行应力-应变场计算时，所加位移约束既要保证计算过程中焊件不发生刚体位移，又不能阻碍其自由形变.由于在实际焊接加工过程中，焊件整体没有外部的拘束，因此，本文在波形钢腹板梁的底板上建立了 3 个节点集合，并对其进行必要位移约束，如图 8 所示.图中：L1为焊缝曲线，L2和 L3分别为结构 1/2 横截面处的底板和腹板曲线.L1L2L3zxyO位移图8位移边界条件Fig.8Displacementboundarycondition 4.2 应力-应变本构关系为保证节点间的数据传输，采用与热学计算相同的有限元网格进行力学计算，将上节温度场计算结果按相同的载荷步依次施加至应力-应变场计算模型的每个节点上.在应力分析中，采用各向同性硬化准则，通过定义钢材在不同温度下的应力-应变曲线来考虑材料的应变硬化效应.材料处于弹性或塑性状态的应力（d）-应变（d）关系16为d=DdCdT,（7）d=de+dp+dth,（8）DCdTdedpdth式中：为弹塑性矩阵，为与温度有关的向量，为温度变化的增量，、分别为弹性应变、塑性应变和热应变的增量.4.3 残余应力场分析分别选取 L1、L2、L33 条路径上的纵、横向残余应力进行分析.定义 x 轴为残余应力纵向，y 轴为横向.4.3.1纵向分布图 9 为波形钢腹板梁焊缝处（L1路径）的残余应力分布.从图 9 可以看出：纵向和横向残余应力具有相同的分布规律；在波形钢腹板梁弯折角处，应力分布集中，发生一定幅度的连续应力变化，幅值高达100MPa；在距离焊缝始、末端 0.09m 范围内纵向和横向残余应力均发生突变，这是因为在起弧和灭弧阶段的能量输入较大，从而造成此阶段的温度梯度也随之变大，最终导致波形钢腹板梁的应力在焊缝始、末两端发生突变；纵向和横向残余应力峰值都出现在焊接收弧阶段，分别为 443.6、206.6MPa.00.20.40.60.81.01.21.44002000200400600800残余应力/MPa焊缝长度/m纵向应力横向应力图9焊缝处的残余应力分布Fig.9Residualstressdistributionattheweld4.3.2横向分布根据圣维南原理可知，焊缝始、末处小范围的残余应力突变对焊件整体残余应力的扰动不大，所以可忽略两端边界条件的影响，取沿焊缝方向中段横截面的残余应力进行研究9.图 10 给出了波型钢腹板梁 1/2 横截面处，沿底板（L2路径）厚度方向 3 条路线的纵向和横向残余应力的分布曲线.由图 10（a）可知：在距离焊缝中心5050mm内，底板横向残余应力变化大，分布集中，在焊缝处达到峰值，远离焊缝中心处缓慢降低并趋于 0；沿底板厚度 3 条不同路径的横向应力均呈现为拉压交替.焊缝中心区域处，横向残余应力在底板上、中侧为压应力，下侧为拉伸应力，其中，最大残余应力是位于底板上侧的压应力，最高达到 107MPa.由图 10（b）可以看出：底板纵向残余应力符合一般理论分布规律，即焊缝及其周围区域受拉，远离焊缝处受压；残余拉应力沿焊缝两边总计约 78mm的范围分布，在焊接中心处达到最大值，其数值约为材料屈服强度的 1.30 倍；底板较窄一侧的残余压应力自焊缝向底板边缘处线性增大，最大值约为材料屈服强度的 0.61 倍；另一侧残余压应力则呈线性减小，其值接近于 0；沿底板厚度方向，纵向残余应力的分布基本一致，上、中、下表面应力最大值分别为第2期冀伟，等：波形钢腹板梁焊接残余应力分布及试验验证293450.3、425.2、401.9MPa；在底板较宽侧，从下表面至上表面，其受力状态从受压过渡为受拉，这可能会对结构的整体受力性能产生影响.00.10.20.30.40.512090603003060横向残余应力/MPa底板宽度/m底板上侧底板中性层底板下侧焊缝(a)横向残余应力00.10.20.30.40.5300100100300500纵向残余应力/MPa底板宽度/m底板上侧底板中性层底板下侧焊缝78 mm0(b)纵向残余应力图10底板处的残余应力分布Fig.10Residualstressdistributionatthebottomplate4.3.3竖向分布图 11 为波形钢腹板梁 1/2 截面处，腹板（L3路径）内、外两侧沿梁高方向的残余应力分布曲线，其中外侧为焊缝所在一侧.由图 11（a）可知：腹板内、外侧的横向残余应力在焊缝附近处分布集中且复杂；距离焊缝中心约 3 倍板厚范围内，横向残余应力呈锯齿状的连续变化，远离焊缝区域，应力表现为一段平稳的长平台，其值趋于 0.由图 11（b）可知：腹板内、外侧的纵向残余应力峰值均出现在焊缝中心处，分别为 431MPa 和 436MPa；在距离焊缝中心50240mm 为腹板内、外侧的受压域，残余压应力最大值为 71MPa，约为材料屈服强度的 0.20 倍；沿腹板高度方向，纵向残余应力逐渐稳定并在腹板端部存在约 6MPa 的应力.4.4 残余内力自平衡分析取本文中 1200 型波形钢腹板梁的 1/2 横截面（即x=600mm）处为研究对象，对图 10（b）和图 11（b）的焊接残余应力分布曲线进行积分，得到腹板和底板焊接残余应力的内力平衡情况，如表 1 所示17.00.20.40.60.81.01.220020406080横向残余应力/MPa腹板高度/m腹板内侧腹板外侧焊缝(a)横向残余应力00.20.40.60.81.01.21000100200300400500腹板纵向应力/MPa腹板高度/m腹板内侧腹板外侧焊缝(b)纵向残余应力图11腹板处的残余应力分布Fig.11Residualstressdistributionattheweb表 1 焊件残余应力合力分析Tab.1CompositionofresidualstressonweldingsectionkN位置拉应力合力压应力合力小计 合计底板上（外）侧220.6201.918.7下（内）侧186.6176.510.1平均203.6189.214.4腹板上（外）侧90.1104.614.54.2下（内）侧67.380.212.9 2.8平均78.792.413.70.7由表 1 可知：底板和腹板的合力分别表现为残余拉力和残余压力；通过内力叠加，焊件在 1/2 横截面处的总合力仅为 0.7kN，故在此横截面处的残余拉、压应力可看作自相平衡，同理可知，焊件整体的残余应力也应处于自相平衡状态.4.5 残余应力试验验证本次焊接残余应力测试试验采用的设备有294西南交通大学学报第59卷CM500 电焊机、BE120-3AA 型应变片、笔记本电脑和东华 DH3816N 静态应变测试仪.为尽可能减小试验误差，测点位置选择在距离焊缝较远处区域.将测点位置处打磨光滑并采用耐高温胶对应变片进行粘贴，通过数据线和网线将应变片、应变测试仪和搭载应变测试系统的电脑相联18.对甘肃博睿重装制造的 NY-XHL1 型波形钢腹板梁 1/8 横截面附近处的底板和腹板相应测点焊接全过程的应变数据进行采集，具体测点布置如图 12 所示.5+15+42035+1.6+25应变片12015120(a)正视图(b)剖面图图12应变片布置（单位：cm）Fig.12Arrangementofstraingauges（unit:cm）波形钢腹板梁完全冷却后，理论上由温度引起的热应变为 0.由于测点位置距离焊缝较远，此区域的塑性应变较小，故假定材料处于弹性状态，按式（9）、（10）反算焊接残余应力6.x=E12(x+y)，（9）y=E12(y+x)，（10）Exyxy式中：、为材料在常温状态下的弹性模量和泊松比，（）和（）分别为纵（横）方向上的应力与应变.由于在焊接过程中受到高温的影响，个别应变片发生部分剥落，对应位置的测点数据误差明显，故舍去.图 13 为波形钢腹板梁 1/8 横截面处，底板上表面和腹板外侧的有限元计算结果与实测值对比.从图 13 可以看出：在靠近焊缝的测量结果与理论计算结果有一定的偏差.究其原因是：在前期热学计算时，由于各种原因导致温度场计算结果并不完美，后续应力-应变场计算结果必然存在偏差；焊缝周围温度过高，导致应变片工作性能下降，因此，实验结果精度有所下降，但上述误差是不可避免的.总体看来，波形钢腹板梁焊接残余应力模拟值与实测值数据基本吻合，大致趋势相同，证明了有限元预测焊接残余应力可行性.00.10.20.30.40.5300100100300500腹板高度/m应力/MPa理论计算值实测值(a)底板上表面00.20.40.60.81.01.21000100200300400500应力/MPa腹板高度/m理论计算值实测值(b)腹板上外侧图13残余应力实测值与理论计算值的比较Fig.13Comparisonbetweenmeasuredandtheoreticallycalculatedvaluesofresidualstress 5 焊接速度对残余应力的影响在实际焊接中，不同的焊接速度对应不同的线能量密度，焊件的温度梯度也会随之改变，最终表现为焊件在不同的焊接速度下，焊接温度场和应力场的差异19.由 4.3 节分析可知，波形钢腹板梁 1/2 横截面底板处的残余应力值最大，应力集中现象最突出.因此，本节提取此处的残余应力，对其在不同焊接速度下的分布与大小进行分析讨论.采用控制变量法，揭示焊接速度对波形钢腹板梁的底板残余应力影响规律.焊接移动速度变化对底板上侧的残余应力分布影响如图 14 所示.图中：200mm/min 为工厂现采用的焊接速度.从图 14 可以看出：横向和纵向的残余应力分布规律不会随着焊接速度的变化而改变；保持其他焊接参数不变，当焊接速度从 150mm/min 增加到250mm/min 时，横向残余压应力和纵向残余拉应力最大值分别增加了 27.11%和 5.88%；从应力峰值的第2期冀伟，等：波形钢腹板梁焊接残余应力分布及试验验证295增加量上来看，横向残余应力对焊接速度的敏感度要高于纵向残余应力.00.10.20.30.40.5150100500500.130.140.150.161201008060横向残余应力/MPa底板宽度/m 150 mm/min 200 mm/min 250 mm/min(a)横向残余应力00.10.20.30.40.53001001003005000.12 0.14 0.16 0.18420435450465纵向残余应力/MPa底板宽度/m 150 mm/min 200 mm/min 250 mm/min(b)纵向残余应力图14在不同焊接速度下的残余应力Fig.14Residualstressatdifferentweldingspeeds 6 结论1）焊接过程中，波形钢腹板梁整体形成准稳态温度场，波形钢腹板和底板沿焊缝前进方向均呈现相同的温度梯度.冷却过程中，冷却速率随着时间的增加而逐渐放缓.2）焊缝处纵向和横向残余应力具有相同变化趋势，且都在波形钢腹板弯折角处发生一定幅度的连续应力变化，幅值高达 100MPa.纵、横向残余应力在灭弧阶段达到最大值，分别为 443.6、206.6MPa.3）底板纵向残余拉应力在焊缝及热影响区内达到最大值，约为材料屈服强度的 1.30 倍；在距离焊缝稍远处呈受压状态，最大压应力出现在钢底板较窄一侧的边缘处，其值约为材料屈服强度的0.61 倍.与完全对称结构不同，波形钢腹板梁在钢底板较宽一侧的边缘处，纵向残余应力由压应力转变为拉应力，这可能会对结构的受力性能产生一定影响.4）焊接速度对波形钢腹板梁的底板残余应力分布扰动不大，但对应力的峰值大小具有较大的影响.当焊接速度从 150mm/min 提升至 250mm/min，横向和纵向残余应力峰值分别提高了 27.11%和5.88%，横向残余应力较纵向残余应力对于焊接速度具有更高的敏感度.参考文献：卫星，邹修兴，姜苏，等.正交异性钢桥面肋板焊接残余应力的数值模拟J.桥梁建设，2014，44（4）：27-33.WEIXing,ZOUXiuxing,JIANGSu,etal.Numericalsimulationofresidualstressinrib-to-topplateweldingoforthotropicsteelbridgedeckJ.BridgeConstruction，2014，44(4)：27-33.1WANGYB,LIGQ,CHENSW.Theassessmentofresidual stresses in welded high strength steel boxsectionsJ.Journal of Constructional Steel Research，2012，76：93-99.2丁阳，郭鹏.高强钢焊接工字梁整体稳定性能分析J.建筑结构，2015，45（21）：25-29.DINGYang,GUOPeng.Investigationonoverallbuck-lingbehaviorofhighstrengthsteelweldedI-beamsJ.BuildingStructure，2015，45(21)：25-29.3PASTERNAK H,LAUNET B,KRAUSCHE T.Welding of girders with thick plates-fabrication,measurementandsimulationJ.JournalofConstructionalSteelResearch，2015，115：407-416.4于海丰，周建伟，张岩，等.小截面焊接工字钢残余应力分布试验研究J.建筑结构学报，2016，37（增 1）：388-392，398.YU Haifeng,ZHOU Jianwei,ZHANG Yan,et al.Experimentalstudyonresidualstressesdistributiononsmall welded I-shaped sectionJ.Journal of BuildingStructures，2016，37(S1)：388-392，398.5KUBOM,WATANABEK.Residualstressmeasure-ment of corrugated steel web sectionsJ.Journal ofJapanSocietyofCivilEngineers，2017，73：248-258.6KOLLR D,KOVESDI B.Welding simulation ofcorrugatedwebgirderspart1:effectofmanufacturingonresidualstressesandimperfectionsJ.Thin-WalledStructures，2020，146：462-477.7王会利，王犇，谢常领，等.钢桥面板 U 肋焊接残余应力影响因素分析J.世界桥梁，2019，47（1）：53-58.WANG Huili,WANG Ben,XIE Changling,et al.AnalysisoffactorsinfluencingweldingresidualstressesinUribsofsteeldeckplateJ.WorldBridges，2019，47(1)：53-58.8曹宝雅，丁幼亮.板件厚度对钢桥面板顶板纵肋焊接9296西南交通大学学报第59卷残余应力的影响分析J.东南大学学报（自然科学版），2016，46（3）：565-571.CAO Baoya,DING Youliang.Influence analysis ofplatethicknessonweldingresidualofsteeldeckdeck-ribJ.JournalofSoutheastUniversity(NaturalScienceEdition)，2016，46(3)：565-571.HUMJ,LIKJ,CAIZP,etal.Anewweldmaterialmodel used in welding analysis of narrow gap thick-walled welded rotorJ.Journal of ManufacturingProcesses，2018，34：614-624.10张建勋，刘川.焊接应力变形有限元计算及其工程应用M.北京：科学出版社，2017.11DENGD,MURAKAWAH.Predictionofweldingdis-tortionandresidualstressinathinplatebutt-weldedjointJ.Computational Materials Science，2008，43(2)：353-365.12顾颖.U 肋加劲钢桥面板焊接残余应力与变形研究D.成都：西南交通大学，2016.13QIANGB,LIYD,YAOCR,etal.Through-thicknessweldingresidualstressanditseffectonstressintensityfactors for semi-elliptical surface cracks in a butt-weldedsteelplateJ.EngineeringFractureMechanics，2018，193：17-31.14周灿丰，陈智，焦向东，等.APIX65 管道深水铺设GMAW 横向焊接温度场J.焊接学报，2020，41（9）：60-68，100.ZHOU Canfeng,CHEN Zhi,JIAO Xiangdong,et al.Study on temperature field of GMAW horizontal15welding for deep water laying of API X65 pipeJ.TransactionsoftheChinaWeldingInstitution，2020，41(9)：60-68，100.秦荣.大跨度桥梁结构M.北京：科学出版社，2008.16赵秋，吴冲.U 肋加劲板焊接残余应力数值模拟分析J.工程力学，2012，29（8）：262-268.ZHAOQiu,WUChong.NumericalanalysisofweldingresidualstressofU-ribstiffenedplateJ.EngineeringMechanics，2012，29(8)：262-268.17强斌，李亚东，顾颖，等.钢桥对接焊缝残余应力及变形场数值分析与试验验证J.铁道学报，2017，39（9