Q460 高强钢焊接T 形截面纵向残余应力分布

熊晓莉，卢娅囡，卢梦丹，管欣旺，李璐

(河南工业大学,郑州,450001)

0 序言

高强钢因具有强度高，延性、韧性好等特点，在船舶、桥梁和钢结构等工程领域中逐渐得到推广应用[1-2].钢构件在制造过程中的钢板切割和焊接将不可避免地在构件中产生残余应力和变形，残余应力作为一种重要的初始力学缺陷会使得高强钢压杆受力后某些部位提前进入塑性状态，降低刚度，并影响其稳定承载能力.

近些年，国内外学者在高强钢焊接截面残余应力分布试验研究方面取得了丰硕的成果，但多数针对焊接工字形截面和焊接箱形截面的残余应力分布规律进行研究，并提出了相应的残余应力分布模型[3-6]，专门针对焊接T 形截面的研究工作开展的并不多.熊晓莉等人[7]采用分割法测量了经等离子切割的Q460 高强钢板焊接成T 形截面的残余应力；Cao 等人[8]采用分割法测量经火焰切割的800 MPa 高强钢板焊接T 形截面的残余应力，探讨了宽厚比的影响，提出两种不同的残余应力分布模型.

基于有限元的数值模拟研究作为一种重要的预测截面残余应力的手段，被越来越多的学者用来预估钢板切割和截面焊接产生的残余应力.一方面，钢板切割过程的数值模拟研究主要聚焦于切割参数对截面残余应力的影响.Jokiaho 等人[9]基于ABAQUS 研究了切割速度和预热对低合金耐磨厚钢板火焰切割残余应力分布的影响.另一方面，截面焊接过程的数值模拟研究主要聚焦于材料属性、约束条件、焊接工艺参数等因素对焊接截面残余应力的影响[10-14]；杜宝帅等人[15]基于ANSYS 研究了Q460 高强钢多层多道焊平板对接焊的残余应力分布；魏雷等人[16]基于ANSYS 研究了Q460 高强钢T 形接头单边开坡口非对称焊的残余应力分布.上述研究也是仅考虑了焊接过程产生的残余应力，未考虑钢板切割初始应力对焊接截面残余应力分布的影响.Wei 等人[17]基于ANSYS 研究了低碳钢平板火焰切割和对接焊连续加工过程后的残余应力分布，研究发现钢板切割初始应力对焊接截面残余应力的影响不容忽视，除此之外还未见有专门针对高强钢切割和截面焊接连续加工过程后截面残余应力分布的相关研究.

采用试验研究与数值模拟相结合的方法，探讨经过钢板切割和截面焊接连续加工过程后的Q460 高强钢焊接T 形截面纵向残余应力分布.采用等离子切割方法制作4 个Q460 高强钢焊接T 形截面残余应力试件，并采用分割法测量上述4 个试件截面的残余应力大小及分布，借助通用有限元程序ANSYS，运用热-结构耦合分析方法模拟钢板切割及截面焊接的过程，最后将试验结果与数值模拟结果进行对比，验证有限元分析过程的正确性.基于正确的有限元分析方法，获得更多不同截面尺寸的Q460 高强钢焊接T 形截面纵向残余应力分布，基于丰富的试验和数值模拟数据分析截面几何参数对纵向残余应力分布的影响，提出Q460 高强钢焊接T 形截面纵向残余应力分布模型，为后续钢压杆整体稳定承载力数值模拟研究提供参考.

1 钢材性能和残余应力测量

1.1 钢材性能

选用8 和12 mm 厚度的Q460 低合金高强度结构钢板，依照国家标准GB/T 2975《钢及钢产品力学性能试验取样位置及试样制备》要求制备试件，并进行室温单向拉伸试验，测量钢材力学性能参数，具体结果如表1 所示.钢材弹性模量、抗拉强度和屈强比都能满足国家标准GB/T 1591《低合金高强度结构钢》和GB50011《建筑抗震设计规范》要求，能够将表1 中数据直接应用于后续研究.

表1 钢材力学性能Table 1 Test results of steel material property

1.2 残余应力测量试件设计及加工

为测量Q460 高强钢焊接T 形截面的残余应力，采用等离子切割方法获取所需截面尺寸的钢板，切割后钢板采用角焊缝连接，焊接时焊丝采用ER55-D2(直径1.2 mm)，采用自动实心焊丝CO2气体保护焊，其焊接电流为220 A，电弧电压为25 V，焊接速度为250 mm/min，电流极性为直流反接，导电嘴与工件距离为20 mm ± 7 mm，保护气体流量为30 L/min，焊脚尺寸为8 mm.使用上述方法共设计制作4 个焊接T 形截面试件，图1 为试件截面几何形状示意图，表2 为试件截面几何尺寸.

图1 试件截面几何形状示意图Fig.1 Schematic diagram of the geometry of the crosssection of the specimen

表2 试件截面几何尺寸Table 2 Sectional dimensions of specimens

1.3 残余应力测量方法

分割法是将截面切割分条使其应力释放，测量切割前后条带的长度，并计算截面残余应力，此测量方法可得到较为准确的截面残余应力及分布，具体测量方法及过程详见文献[7].

1.4 残余应力测量结果

截面纵向残余应力测量结果可详见文献[7]，根据测量结果可以看出，Q460 高强钢焊接T 形截面纵向残余应力有如下特点.

(1) 翼缘残余应力呈对称分布，翼缘焊缝周边区域和外伸端为残余拉应力，翼缘外伸段中部为残余压应力.

(2) 腹板焊缝周边区域和外伸端为残余拉应力，腹板中部为残余压应力.

(3) 截面焊缝处的最大残余拉应力小于钢材实测屈服强度，也小于钢材名义屈服强度.

2 焊接残余应力数值模拟试验

2.1 有限元模型建立

钢板切割的有限元模型尺寸与试件尺寸一致(图2)，割缝宽度设为1 mm，钢板切割几何形状示意图和切割方向如图2a 所示；根据表2 中T 形截面几何尺寸建立有限元模型，T 形截面焊接几何形状示意图和焊接方向如图2b 所示；有限元模型网格划分时，为了能够在得到准确计算结果的同时尽量减少计算时间成本，在割缝(焊缝) 和热影响区内采用较细的网格划分，距离割缝(焊缝) 和热影响区越远，网格划分越稀疏.图3 为钢板切割和T 形截面焊接有限元模型的网格划分示意图.

图2 钢板切割几何形状和T 形截面焊接示意图Fig.2 Schematic diagram of steel plate cutting and Tsection welding.(a) steel plate cutting; (b) Tsection welding

图3 有限元模型网格划分示意图Fig.3 Schematic diagram of meshing of finite element model for steel plate cutting and T-section welding.(a) steel plate cutting; (b) T-section welding

2.2 材料属性定义

国家标准GB50017—2017《钢结构设计标准》[18]并未对不同温度下的钢材热物理性能参数进行描述，因此参考欧洲标准Eurocode 3[19]确定相关参数.密度不随温度变化而变化，取恒定值ρ=7 850 kg/m3；泊松比也与温度变化无关，可根据第1.1 节材料性能试验测量结果取值，板厚为8 mm 的钢材泊松比取恒定值0.275，板厚为12 mm的钢材泊松比取恒定值0.285.当施加热源时，温度发生变化，钢材的热物理性能和力学性能也会随之改变，变化曲线如图4 所示，图4 中钢材常温下的屈服强度和弹性模量采用第1.1 节材料性能试验的实测值.后续3.1 节将补充板厚为10 mm 的有限元模型进行参数分析，考虑到材料性能试验并未测量板厚为10 mm 的钢材力学性能参数，因此采用线性内插法得到板厚为10 mm 的相关性能参数.

图4 钢材热物理性能和力学性能随温度变化Fig.4 Variation curve of thermal properties and mechanical properties of steel with temperature

2.3 热源模型选择

现有相关研究中常用的热源模型有高斯热源模型、双椭球热源模型、基于生死单元的焊接热源加载模型及表面高斯热源和生热率组合热源模型等.基于生死单元的焊接热源加载模型无需构造复杂的热流密度函数，能够在保证计算精度的前提下节省计算时间，因此将采用基于生死单元的焊接热源加载模型对钢板切割过程和T 形截面焊接过程进行数值模拟.

钢板切割是将热源施加在割缝处，使割缝处钢材熔融形成熔渣，高速热流使熔渣脱落形成切口.因此在钢板割缝处施加生热率模拟钢板切割过程，生热率Q的计算公式为

式中：η为钢板切割时输入热量的效率，取0.35；I为切割电流，取220 A；U为切割电压，取25 V；VH为单位荷载步割缝的体积.切割速度为250 mm/min，两道割缝之间设置冷却时间为400 s.

截面焊接时热源施加在焊缝处，使焊条熔化与钢材熔融，并冷却后形成焊缝.可在焊缝处直接施加生热率模拟T 形截面焊接过程，此时式(1) 中η为焊接时输入热量的效率，取0.85；I为焊接电流，取220 A；U为焊接电压，取25 V；VH为单位荷载步焊缝的体积.焊接速度为250 mm/min，两道焊缝之间设置冷却时间为400 s，两道焊缝的焊接参数保持不变，焊接结束后自然冷却至室温.有限元模型未考虑切割顺序和焊接顺序对残余应力的影响，每个有限元模型均采用同一切割和焊接顺序.

2.4 热分析

钢板切割过程与T 形截面焊接过程均采用热-结构耦合分析方法.在热分析阶段，有限元模型选用能够适用于瞬态热分析的三维热实体单元SOLID70.为模拟割缝和焊缝与空气对流换热，钢板切割过程和T 形截面焊接过程设置对流换热系数为10 W/(m2·K)，环境温度为20 °C.

采用 “生死单元” 法模拟钢板割缝脱落和T 形截面焊缝形成过程.切割时，对割缝单元依次施加热源，为模拟切缝处熔渣脱落，热源施加完毕后杀死割缝单元；焊接时，先将焊缝单元杀死，模拟焊接准备过程，然后依次激活被杀死单元并施加热源，在焊接过程中，假定焊缝与钢材具有相同的性能参数.

2.5 结构分析

在结构分析阶段，采用ETCHG 命令将三维热实体单元SOLID70 转化为结构单元SOLID185.塑性选项使用双线性等向强化模型BISO，定义随温度变化的屈服强度和切线模量，为防止模型发生位移，设置边界条件，钢板切割的边界条件和T 形截面焊接的边界条件分别如图2 所示.热分析对应时刻的节点温度在结构分析中作为载荷施加到相应节点进行求解，利用ANSYS 后处理命令读取相应位置的纵向残余应力.

2.6 有限元分析结果

为模拟钢板切割和T 形截面焊接的连续加工过程，可将钢板切割和T 形截面焊接分别得到的纵向残余应力值叠加，即可获得考虑钢板切割初始应力的焊接T 形截面纵向残余应力值.考虑到钢板若在同一位置先切割后焊接，切焊耦合应力场由焊接产生的残余应力决定[20]，因此在焊接T 形截面腹板的焊缝处，纵向残余应力由焊接产生的残余应力决定，其它位置的纵向残余应力则为钢板切割和截面焊接叠加后获得.对残余应力测量试验截面ST-1～ST-4 的焊接过程进行数值模拟，分别获取图2b 中横截面A-A 位置处钢板切割、T 形截面焊接和二者叠加的截面纵向残余应力分布，结果汇总于图5，从图5 可以得出如下结论.

图5 钢板切割、T 形截面焊接和切焊叠加的截面纵向残余应力分布Fig.5 Longitudinal residual stress distribution of steel plate cutting,T-section welding and superimposed cutting and welding.(a) ST-1;(b) ST-2;(c)ST-3;(d) ST-4

(1) 仅钢板切割时，钢板沿厚度的纵向残余应力分布相同；除切割边缘处为残余拉应力外，钢板其余部分均为残余压应力.究其原因，钢板切割时基于生死单元的焊接热源模型是均匀施加在割缝处所有单元的，因此沿钢板厚度方向的纵向残余应力不发生改变；近热影响区为纵向残余拉应力，反之则为纵向残余压应力.

(2) 仅T 形截面焊接时，焊缝周边区域和腹板外伸端为残余拉应力，翼缘外伸端和腹板中部为残余压应力.焊接时焊缝处的高温需要通过钢板由近及远耗散，近焊缝区域一般产生残余拉应力，远离焊缝的区域一般为残余压应力，最终达到整个截面上残余应力自平衡状态.由于焊缝位于翼缘的上表面且与腹板相交位置处，热量沿板厚方向的传导因受到周边高温单元的影响，会出现近焊缝处翼缘和腹板的纵向残余应力沿板厚方向分布不均匀的现象，且板厚越大，这种不均匀现象越明显.此外还需要注意的是焊缝周边的残余拉应力有时会超过钢材的名义屈服强度，甚至会超过钢材的实测屈服强度，这主要是因为热源集中在焊缝处，导致部分单元提前进入了材料的强化阶段.然而考虑到超过屈服强度的部分一般在焊接T 形截面的近焊缝处，该位置的纵向残余应力其对钢压杆的整体稳定承载力影响不大，因此在后续残余应力分布模型的选取时可偏安全地将焊缝处的最大残余拉应力按钢材的名义屈服强度考虑.

(3) 将钢板切割和T 形截面焊接的纵向残余应力进行叠加后，焊缝周边区域、翼缘和腹板外伸端为残余拉应力，翼缘外伸段中部和腹板中部为残余压应力.与仅考虑T 形截面焊接时的纵向残余应力相比，叠加钢板切割加工过程产生的残余应力后，T 形截面焊缝周边区域的残余拉应力值略有减少，腹板外伸端的残余拉应力值明显变大，翼缘外伸段中部和腹板中部的残余压应力明显变大，翼缘外伸端的残余应力则由压应力变为了拉应力.上述变化表明，钢板切割产生的残余应力对Q460 高强钢焊接T 形截面残余应力分布的影响不能忽略.

2.7 数值模拟结果与试验结果对比

为了验证有限元分析结果的正确性，选取图2b 中横截面A-A 位置处的数值模拟结果与1.4 节中的试验结果进行对比，具体情况如图6 所示，从图6 可以得出如下结论.

图6 截面纵向残余应力分布的试验结果和数值模拟结果对比Fig.6 Comparison of test results and numerical simulation results of longitudinal residual stress distribution in the section.(a) ST-1;(b) ST-2;(c)ST-3;(d) ST-4

(1) 焊接T 形截面纵向残余应力分布情况的数值模拟结果与试验结果一致，即焊缝周边区域、翼缘和腹板外伸端为残余拉应力，翼缘外伸段中部和腹板中部为残余压应力.

(2) 翼缘焊缝附近的残余拉应力的数值模拟结果与试验结果有差异，但这种差异主要是由有限元方法与试验方法的精度不同所导致的.有限元方法可以精细地描述焊接热影响区的残余应力分布情况，具体表现为有限元分析所得翼缘上、下表面的残余拉应力明显不同.主要原因是翼缘上表面在焊缝位置处与腹板相连，焊接时热量输入集中，且降温速度较慢导致该处产生较大的残余拉应力，翼缘下表面相对降温速度较快导致该处产生相对较小的残余拉应力.用分割法进行残余应力测量时，在翼缘和腹板相交处，切条的难度较大，甚至出现条带上、下表面的残余应力实测值相同的情况，这也是由分割法在焊缝处使用受限决定的，但是翼缘焊缝附近上、下表面残余拉应力平均值的数值模拟结果与试验结果却相当接近.

(3) 腹板外伸端残余拉应力的数值模拟结果与试验结果吻合较好.考虑到该处的残余拉应力对钢压杆的整体稳定承载力有利，偏于安全地取平均值进行对比.以ST-3 为例，腹板外伸端平均残余拉应力的有限元分析结果为249.99 MPa，试验值为245.96 MPa，前者仅比后者高出1.6%，且未超过钢材名义屈服强度的5%，然而翼缘外伸端残余拉应力的数值模拟结果与试验结果差异较大，这可能与现场试验过程中的约束施加及人工测量的误差有关，但考虑到该处的残余拉应力对钢压杆的整体稳定承载力有利，后续研究过程中可偏安全地在该处取较小的残余拉应力值.

(4) 翼缘外伸段中部和腹板中部的残余压应力的数值模拟结果与试验结果吻合较好.以模型ST-3 为例，翼缘左、右两个外伸段中部的最大残余压应力的有限元分析结果分别为−146.92 和−121.73 MPa，忽略两道焊缝先后焊接次序不同引起翼缘左右两外伸段残余应力的差异，两者的平均值为−134.33 MPa；试验值为−145.20 和−125.60 MPa，平均值为−135.4 MPa；对比有限元分析结果的平均值和试验结果的平均值可知，前者绝对值仅比后者的低0.79%.同样腹板中部的最大残余压应力的有限元分析结果为−119.09 MPa，试验值为−113.8 MPa，前者绝对值也仅比后者的高4.64%，翼缘和腹板的残余压应力均值的有限元分析结果与试验结果相差均未超过钢材名义屈服强度的5%.

综上所述，翼缘和腹板纵向残余应力分布的数值模拟结果与试验结果总体吻合较好.虽然在焊缝周边区域，出现了按照两种不同方法获得的残余拉应力不一致的情况，但考虑到这部分区域位于截面中间位置，相应的残余拉应力对焊接T 形截面钢压杆整体稳定承载力的影响不大，而对高强钢焊接T 形截面钢压杆整体稳定承载力影响较大的翼缘和腹板外伸端残余拉应力、中部残余压应力的数值模拟结果与试验测量值均吻合较好，因此可以认为提出的有限元分析方法正确合理，可用于后续补充模型的计算，为参数分析提供更多的数据支撑.

3 参数分析

3.1 补充模型设计

为分析截面几何参数(板件宽厚比和板厚) 对Q460 高强钢焊接T 形截面纵向残余应力的影响，考虑到前述试验及有限元分析获得的残余应力数据有限，特补充设计6 种不同截面尺寸的有限元模型，具体截面尺寸如表3 所示.

表3 补充模型的截面几何尺寸Table 3 Section geometry of supplementary models

3.2 纵向残余应力分布汇总

按照前述的焊接残余应力数值模拟试验过程，对钢板切割和焊接连续加工过程进行数值模拟，获得各截面的纵向残余应力分布，如图7 所示.

图7 补充模型的纵向残余应力分布Fig.7 Longitudinal residual stress distribution for supplementary model.(a) STA-1;(b) STA-2;(c) STA-3;(d) STA-4;(e) STA-5;(f) STA-6

为了便于后续问题研究，将数值模拟得到的Q460 高强钢焊接T 形截面纵向残余应力分布进行简化，简化的纵向残余应力分布模型如图8 所示，图中的各截面纵向残余应力代表值汇总于表4.图8 和表4 中，σfrt1为翼缘焊缝区域的最大残余拉应力；σfrt2和σfrt3为翼缘外伸端的最大残余拉应力；σfrc1和σfrc2为翼缘外伸段中部的最大残余压应力；σwrt1为腹板焊缝区域的最大残余拉应力；σwrt2为腹板外伸端的最大残余拉应力；σwrc为腹板中部的最大残余压应力；a～j分别表示纵向残余应力各线性分布区段的长度.

图8 简化的纵向残余应力分布模型Fig.8 Simplified longitudinal residual stress distribution model

表4 截面纵向残余应力代表值汇总Table 4 Summary of representative values of longitudinal residual stress in section

3.3 板件宽厚比对纵向残余应力的影响

图9 绘制了板厚为8 和12 mm 截面的宽厚比与纵向残余应力之间的关系.从图中可以看出，翼缘与腹板焊缝周边区域的残余拉应力和翼缘与腹板外伸端残余拉应力均与宽厚比无关；翼缘外伸段中部和腹板中部的残余压应力与宽厚比成反比.究其原因，翼缘焊缝处的残余拉应力与截面焊接时的热输入能量有关，与板宽厚比无关；翼缘和腹板外伸端为切割边，该处的残余拉应力仅与钢板切割时的热输入能量有关，与宽厚比无关；但因较小宽厚比的板件外伸端更易受到焊接后板件中部产生的残余压应力的影响，因而出现了翼缘外伸端残余拉应力明显与宽厚比成正比、腹板外伸端残余拉应力略微与宽厚比成正比的现象，因此从根本上讲，翼缘和腹板外伸端残余拉应力与宽厚比关系不大.翼缘外伸段中部和腹板中部的残余压应力主要是由钢板切割和T 形截面焊接后钢材的不均匀冷却造成的，热输入能量大小不变，宽厚比越大，残余压应力的分布范围越大，相应的残余压应力越小，当宽厚比较大时(H/tw＞20)，残余压应力的分布范围对应力值的影响不明显.

图9 板件宽厚比与纵向残余应力的关系Fig.9 Relationship between width-thickness ratio of plate and longitudinal residual stress.(a) plate thickness is 8 mm;(b) plate thickness is 12 mm

3.4 板厚对纵向残余应力的影响

为方便对比，选取表3 中宽厚比相同而板厚不同的两组截面，第一组截面为STA-1，STA-3 和STA-5，第二组截面STA-2，STA-4 和STA-6.图10为第一组和第二组截面的板厚与纵向残余应力之间的关系，从图中可以看出，焊接工艺参数不变的情况下，翼缘和腹板焊缝周边区域的最大残余拉应力与板厚成正比；翼缘和腹板外伸端的残余拉应力与板厚无关；翼缘外伸段中部和腹板中部的残余压应力与板厚成反比.

图10 板厚与纵向残余应力之间的关系Fig.10 Relationship between plate thickness and longitudinal residual stress.(a) the first set of sections;(b) the second set of sections

结合图10、图7 和表4 中的数据可知，板厚越大，虽然翼缘和腹板焊缝周边区域的最大残余拉应力越大，但翼缘和腹板沿板厚方向残余拉应力的分布差异也越大，且部分板厚较大的截面近焊缝位置附近的钢材已进入了强化阶段；考虑到焊缝处的残余拉应力仅与焊接热输入能量有关，在实际取值时可以认为焊缝处的残余拉应力上限为钢材名义屈服强度，与板厚无关；此外随着板厚的增大，翼缘外伸段中部和腹板中部的不均匀冷却速度变慢，最终形成的残余压应力变小.

3.5 截面纵向残余应力的自平衡状态

翼缘、腹板和全截面的不平衡应力σerr计算式为

式中：n为截面划分单元的数量；σri为各单元对应的残余应力，拉应力取正值，压应力取负值；Ai为各单元的面积；A为截面的总面积.

按照式(2) 计算得到各有限元模型的板件不平衡应力，并计算出不平衡应力系数，将两者汇总于图11 中.不平衡应力系数为板件不平衡应力与实测屈服强度的比值.从图中可以看出，所有模型的各板件不平衡应力系数均未超过5%，可以认为截面纵向残余应力满足自平衡条件.

图11 截面纵向残余应力不平衡应力Fig.11 Section longitudinal residual stress unbalanced stress

4 纵向残余应力分布模型

4.1 分布模型形状

文献[8] 提出的800 MPa 高强钢焊接T 形截面残余应力分布模型，结合前述试验及有限元分析结果，Q460 高强钢焊接T 形截面纵向残余应力分布模型仍采用图8 所示的分布模型形状，翼缘纵向残余应力对称分布(σfrt2=σfrt3，σfrc1=σfrc2).

由3.5 节分析可知，各板件纵向残余应力满足自平衡条件，因此翼缘和腹板纵向残余应力分别满足如下要求，即

式中：σfr和σwr分别为翼缘和腹板的残余应力；Af和Aw分别为翼缘和腹板的面积.

由图8 可知，翼缘和腹板的纵向残余应力分布范围应分别满足下式要求

式中：a～e为翼缘纵向残余应力各线性分布区段的长度；f～j为腹板纵向残余应力各线性分布区段的长度.

4.2 残余拉应力取值

根据前述分析可知，翼缘和腹板焊缝周边区域的残余拉应力对钢压杆的稳定承载力有利，忽略其沿板厚方向分布的不均匀性以及部分钢材进入强化阶段的影响，可偏安全地取钢材名义屈服强度fy，翼缘残余应力分布范围e取hf+tw/2；腹板残余应力分布范围f取hf.

翼缘和腹板外伸端的残余拉应力与宽厚比无关，考虑到板件外伸端的残余拉应力对钢压杆的整体稳定承载力有利，偏于安全地将取值进行如下简化：翼缘外伸端残余拉应力取0.15fy，分布范围a取(bf−hf)/10，翼缘残余应力过渡段b的分布范围取(bf−hf)/10；腹板外伸端残余拉应力0.65fy，分布范围j取(H−hf)/10，腹板残余应力过渡段i的分布范围取(H−hf)/10.

4.3 残余压应力取值

前述分析表明，翼缘和腹板上的最大残余压应力与宽厚比和板厚均成反比，因此对前述试验和有限元分析结果进行非线性拟合，可得翼缘外伸段中部和腹板中部的残余压应力为

式中：σfrc为翼缘外伸段中部的残余压应力；σwrc为腹板中部的残余压应力.

纵向残余应力分布模型中的残余应力分布范围计算公式汇总于表5，残余应力分布模型及相应的计算公式汇总于图12.根据表5 和图12 计算可得拟合曲线并将其绘制于图7.将拟合曲线与图7 中的数值模拟结果进行对比，可以看出两者吻合较好.因此，可将图12 中的Q460 高强钢焊接T 形截面纵向残余应力分布模型应用于后续钢压杆整体稳定承载力研究.

图12 Q460 高强钢焊接T 形截面纵向残余应力分布模型Fig.12 Model of longitudinal residual stress distribution in welded T-section of Q460 high-strength steel

表5 残余应力分布范围计算公式汇总Table 5 Summary of calculation formulas for residual stress distribution range

5 结论

(1) 采用分割法测量了高强钢焊接T 形截面试件的纵向残余应力，采用等离子切割方法下料并焊接而成的T 形截面，其焊缝周边区域、翼缘和腹板外伸端为残余拉应力，翼缘外伸段中部和腹板中部为残余压应力.

(2) 借助通用有限元程序，运用热-结构耦合分析方法模拟钢板切割及截面焊接的过程，获得高强钢焊接T 形截面纵向残余应力的数值模拟结果，该结果与试验结果吻合较好，从而验证了有限元分析方法的正确性.

(3) 基于合理的有限元分析方法，改变截面几何参数，获得更多不同截面尺寸的高强钢焊接T 形截面的纵向残余应力分布.以丰富的试验及有限元分析数据为基础，分析截面几何参数对纵向残余应力的影响，得出结论，残余拉应力与板件宽厚比和板厚无关；残余压应力随板件宽厚比和板厚的增大而减小；翼缘、腹板和全截面均满足自平衡条件.

(4) 对试验数据和数值模拟结果进行分段线性拟合，建立Q460 高强钢焊接T 形截面纵向残余应力分布模型，为后续钢压杆整体稳定承载力研究过程中的初始缺陷取值提供重要参考.