摘 要：为了研究不锈钢工字形轴心受压构件的整体稳定承载力，采用ANSYS软件建模对不锈钢轴压构件进行非线性有限元模拟，将模拟结果与试验结果进行对比，验证了所建立的有限元模型的准确性.采用经试验验证的有限元模型对不同几何初始缺陷、截面残余应力、材料力学性能、截面宽厚比以及长细比的不锈钢工字形构件整体稳定承载力进行了参数分析，通过对比可确定材料力学性能、构件长细比为主要影响承载力的因素.在参数分析的基础上通过稳定承载力的数据拟合提出了整体稳定系数的三段式计算方法，并将该计算方法与试验数据进行对比，表明此计算方法可以较为准确地计算不锈钢工字形轴心受压构件整体稳定承载力.
关键词：不锈钢；有限元法；整体稳定性；参数分析；计算方法
中图分类号：TU391 文献标识码：A
文章编号：1674-2974（2016）03-0055-11
不锈钢材料具有易维护和全生命周期成本低等优势，在建筑结构中得到广泛的应用[1].目前有关不锈钢结构的研究比较广泛，但研究尚处于初期阶段.对于不锈钢材料的应力应变关系，Mirambell和Real[2]提出了两阶段的Ramberg-Osgood模型，此外，Quach等[3]通过引入σ2.0，提出了一个三阶段修正模型.Gardner等[4]和Quach等[5]分别对不同情况构件提出了残余应力模型.王元清[6]和杨璐[7]等进行了相应的试验和理论研究，提出了纯弯作用下焊接工字形不锈钢梁的整体稳定计算表达式.Yuan等[8]对不锈钢焊接截面轴心受压构件的相关稳定性能进行了试验研究，并提出了相应的设计方法.舒赣平等[9-10]对冷弯不锈钢轴心受压构件的稳定承载力以及压弯构件平面稳定承载力进行了研究.Gardner等[11-12]根据若干不锈钢试件的试验研究结果进行分析，总结并提出了连续强度法[13].在不锈钢受压构件研究方面，Gardner[14-15]等对不锈钢的短柱试件进行了试验，对比了现行欧洲不锈钢设计规范的计算结果，提出了新的设计方法和设计建议.此外还有一些学者对不锈钢构件的变形性能和不锈钢管混凝土进行了研究.
不锈钢的应力应变曲线特征以及构件受力性能特征与普通钢有很大的不同，不锈钢材料具有明显的非线性特征以及明显的应变硬化特性，因此在进行不锈钢整体稳定性能分析时不宜使用钢结构规范.目前，中国不锈钢设计规范正在编制中，对不锈钢轴心受压构件整体稳定性能的研究对规范的编纂有一定的价值.
1 试验概况
1.1 试件设计
本文对奥氏体型和双相体型焊接工字型不锈钢共计22根构件进行了轴压试验，根据试验设定其中12个构件绕弱轴出现失稳，10个构件绕强轴出现失稳.所有构件在设计前均进行了试算，确保施加荷载能够使构件发生整体失稳，同时限制构件的截面尺寸以防构件出现局部屈曲.构件的材料属性通过材性实验获得.
1.2 试验装置
试验过程中使用的加载装置如图1所示，采用500T液压式长柱压力试验机进行加载，加载过程中，试件两端各布置一个单刀铰，使得加载装置与柱子端部实现单向铰接，单刀铰的转动轴线与试件弯曲失稳平面垂直.通过单刀铰约束方向来控制构件绕强轴或弱轴的失稳，单刀铰转动中心至柱端面距离为170 mm，试件的铰接长度Lt=L+340.
1.3 测量内容
位移计架设示意图如图2所示，在柱中设置2个位移计LVDT-5和LVDT-6用于测量试件失稳平面内柱中截面的水平位移，同时在柱中失稳平面外设置一个位移计LVDT-7用于测量试件失稳时失稳平面外柱中截面的水平位移.通过LVDT-3和LVDT-4两个位移计可以测量试件在受压时的竖向变形，即试验仪器加载点的竖向位移.此外在柱两端截面各布置了4个应变片，用于根据试件在弹性受力阶段端部截面的应变分布推算荷载初偏心值.在试验开始前，采用光学测量设备通过测量沿柱长方向四分点位置处截面中心偏离柱两端截面中心连线的距离来对每个试件的整体几何初始弯曲进行了测量.
2 有限元方法及试验验证
对工字形构件的整体稳定性能的研究需要考虑不同的影响因素，并分别进行参数化分析，由于试验本身的局限性，需引入有限元软件进行分析.在建立有限元模型及计算分析的过程中考虑了材料的非线性、构件几何初始缺陷以及截面残余应力的因素，并用试验结果对有限元模型进行了验证.
2.1 有限元模型的建立
2.1.1 单元的选择及边界约束
本文主要对轴压构件的整体稳定性能进行研究，因此采用BEAM188单元.为了与试验的柱端约束情况取得一致，首先在有限元模型的两端采用固定约束，在此基础上释放特定方向的扭转.对于工字形截面构件，应对构件绕强轴与弱轴失稳分别进行考虑[16].
2.1.2 不锈钢材料的本构模型
不锈钢的本构模型中比较准确的有两阶段的R-O模型和三阶段模型，其中两阶段模型相对简洁，使用较多.为能更好地反映试验的实际情况，本文采用试验测得的3条应力应变曲线的平均值模型.有限元分析采用多线性等向强化本构模型进行模拟.应力应变试验曲线与平均值模型曲线的对比如图3所示.
2.1.3 构件的初始缺陷
目前，在钢结构稳定计算中，各国规范都考虑了构件的初始几何缺陷.在进行有限元模拟时可偏于安全地取一阶整体屈曲模态作为几何初始缺陷的变形状态，用总的初始缺陷作为构件的一阶模态的最大位移，并对模型节点的坐标进行更新以实现对初始缺陷的模拟.初始几何缺陷模型如图4所示.
2.1.4 构件残余应力的分布
焊接构件中普遍存在残余应力，且残余应力的存在会对构件的极限承载力产生影响.本文采用袁焕鑫[17]测得的残余应力分布图[图5（a）]和他提出的残余应力分布模型[图5（b）].根据残余应力自平衡特点将残余应力分布模型进行简化，如图5（c）所示，并将其施加在有限元模型中，如图5（d）所示[17]. 2.2 有限元模型与试验对比
2.2.1 荷载位移曲线对比
采用上述方法对22根不锈钢焊接工字形截面轴压构件进行有限元模拟，有限元模拟得到各构件的荷载与端部的竖向位移曲线与试验曲线的对比如图6所示，模拟得到各构件的荷载与构件中点处水平位移曲线与试验曲线的对比如图7所示，其中取构件内侧与外侧的残余应力测量值的平均值作为实测值.
由图6和图7可知，试验值与模拟值的荷载位移曲线匹配得较好，有限元模型能够较准确地模拟奥氏体型与双相体型不锈钢构件整体失稳的真实受力情况.由于试验需要克服单刀铰的摩擦力，无法达到理想状态；另一方面，在进行有限元模拟时，构件初始缺陷是按构件最大初始缺陷的一阶模态来取用的，因此，试验和数值模拟曲线之间存在一定的差别.此外，模拟还受到了构件的计算长度以及残余应力模型简化的影响.
2.2.2 极限荷载的对比
目前，欧洲不锈钢规范以及美国不锈钢规范是不锈钢设计的主要依据.本文对比了模拟和试验得到的极限承载力，同时也分别按欧洲不锈钢规范和美国不锈钢规范对稳定承载力进行了计算，通过对比验证了有限元模拟的准确性，并对两种规范的合理性进行了评估.
由表1可知，双相体型的模拟值与试验值的差别比奥氏体型模拟值与试验值的差别小，这主要是由于双相体型构件在极限荷载状态下未达到屈服状态，构件均属于弹性失稳；构件绕强轴失稳时模拟值与试验值的差别比构件绕弱轴失稳时的差别小.总体来看，有限元模拟值与试验值的误差控制在10%以内，平均误差在5%以内，本文的有限元方法能够很好地模拟不锈钢轴心受压构件整体稳定承载力.
对于不锈钢构件，根据欧洲不锈钢规范得到的设计值普遍低于试验值和模拟值，随长细比的增加两者的差异逐渐变小，欧洲不锈钢规范相对较保守且符合实际情况.根据美国不锈钢规范得到的设计值整体上高于试验值与模拟值，美国规范相对较激进不适用于焊接不锈钢构件.当构件长细比较小时，根据美国规范得到的设计值与构件极限承载力差别较小，而对于长细比较大的构件按照欧洲规范得到的设计值与构件的极限荷载差别较小.
综合来看，试验与模拟之间的差别在可接受的范围之内，故此有限元方法比较准确.
3 参数分析
本节使用以上经试验确定的有限元方法进行参数化分析，分别讨论构件几何初始缺陷、截面残余应力、材料力学性能、截面宽厚比以及长细比对构件极限承载力的影响.
3.1 构件几何初始缺陷
为研究构件几何初始缺陷对构件极限承载力的影响，本文在不改变其他参数情况下，采用改变初始缺陷大小的方法进行验证.本文选取初始缺陷系数分别为0.000 5L，0.001L和0.002L时构件绕强轴失稳及绕弱轴失稳2种情况进行有限元分析，并将初始缺陷值为0.000 5L与0.002L时的极限荷载值与初始缺陷值为0.001L时的极限荷载值进行对比.对比结果如图8所示.
图8中横坐标表示构件长细比，纵坐标表示初始缺陷分别为0.002L，0.000 5L构件的极限荷载F0.002，F0.000 5与初始缺陷为0.001L构件的极限荷载F0.001的比值.由图8可知，初始缺陷为0.000 5L时的极限荷载比初始缺陷为0.001L的极限荷载高5%左右，且随正则化长细比的变化有轻微的波动；初始缺陷为0.002L时的极限荷载比初始缺陷为0.001L的极限荷载低7%左右，且随正则化长细比的变化有轻微的波动.初始缺陷大小对构件的极限承载力影响较大，但对不同长细比的构件，初始缺陷对于绕弱轴失稳和绕强轴失稳2种情况的极限承载力影响不大.此外，通过奥氏体型与双相体型对比可知，不同初始缺陷对两者的影响基本相同，但初始缺陷对双相体型构件影响较小且对绕弱轴失稳构件影响更大.
3.2 构件截面残余应力的影响
为研究构件截面残余应力的影响，本文在不改变其他参数情况下，分别对考虑残余应力与不考虑残余应力2种情况进行了有限元分析，并对考虑残余应力时构件的极限荷载与不考虑残余应力时的极限荷载进行了对比.对比结果如图9所示.
图9中横坐标表示构件长细比，纵坐标表示不考虑残余应力时构件的极限荷载Fw与考虑残余应力的极限荷载Fy的比值.从构件的残余应力影响来看，随长细比的变化，残余应力对构件的极限承载力的影响也随之发生改变；残余应力对构件绕弱轴失稳时的极限承载力影响较大，表明绕弱轴失稳时构件对残余应力更为敏感.此外，残余应力对奥氏体型不锈钢构件的极限承载力的影响较大，而对双相体型不锈钢构件的极限承载力的影响较小，主要是由于残余应力的存在使得奥氏体型不锈钢构件更容易达到屈服应力fy.
3.3 材料力学性能的影响
在探讨材料力学性能的影响时，在保证其他影响因素不变的前提下，分别取用2种奥氏体型与2种双相体型不锈钢材料进行有限元分析，并对有限元模拟得到的极限荷载进行对比，对比结果如图10所示.
图10中横坐标表示构件长细比，纵坐标表示2种牌号的材料极限荷载F1和F2的比值.由图10可知，绕强轴失稳和弱轴失稳时材料力学特性对整体稳定承载力的影响均表现为随构件长细比的增加先变大后变小.材料的力学性能对双相体型不锈钢构件的极限承载力影响较大，同一种材料使用正则化长细比考虑材料特性后2种牌号的材料稳定系数与正则化长细比的关系基本相同.
3.4 截面宽厚比的影响
考虑宽厚比的影响时，在不改变其他参数的情况下，通过改变截面宽厚比，使用有限元分析获得极限荷载并对比按照中国钢结构规范中φ=F/fyA计算得到稳定系数的变化.如图11所示.
图11中横坐标表示构件宽厚比，纵坐标表示不考虑残余应力时与考虑残余应力时的整体稳定系数.由图11可知，腹板的宽厚比和翼缘的宽厚比对绕弱轴失稳和绕强轴失稳的极限承载力的影响可以忽略，且宽厚比对奥氏体型与双相体型不锈钢构件极限承载力的影响相同. 3.5 构件长细比的影响
当讨论长细比对构件的极限荷载的影响时，在不改变其他因素的前提下，通过改变构件的长度来改变构件的长细比，通过有限元分析确定当构件长细比改变时构件极限荷载按照中国钢结构规范中φ=F/fyA计算得到整体稳定系数变化，如图12所示.
图12中横坐标表示构件长细比，纵坐标表示构件的整体稳定系数.由图12可知，奥氏体型与双相体型不锈钢的极限荷载随长细比的变化均呈现出较明显的改变，构件长细比越大极限荷载值越小.
4 计算方法的提出
4.1 计算方法
本文采用正则化长细比构建整体稳定系数，将柱子曲线分为3段分别进行计算.第1段，由于正则化长细比较小，构件失稳时由于截面屈服应力已超过fy，边缘屈服准则已不再成立，因此对于此类问题可采用Gardner提出的连续强度法.由于此种方法计算比较复杂且在实际中此类长细比较小的构件较少，因此为了使用方便以及曲线的完整可采用一段保守的函数曲线来代替.第2段，考虑到材料的非线性特性，构件失稳形式为弹塑性失稳，综合美国钢结构规范以及冷弯构件的直接强度法，可采用美国冷弯不锈钢规范中给出的公式形式进行计算.第3段，此时构件正则化长细比较大，构件失稳形式为弹性失稳，可采用基于构件边缘屈服准则的perry公式的形式进行计算.
4.2 数据拟合
通过对奥氏体型与双相体型不锈钢进行参数化有限元分析，将得到的有限元参数分析结果使用Matlab提出的稳定系数公式进行拟合，最终确定公式中系数的取值，并考虑到系数的简便以及曲线的连续性对系数进行简单处理，如表2所示.
将通过计算得到的柱子曲线与欧洲不锈钢规范和美国不锈钢规范得到的曲线进行对比可以看出：计算得到的柱子曲线高于欧洲规范曲线，且长细比越大两者柱子曲线的差别越小；当正则化长细比较小时，计算得到的柱子曲线高于规范曲线，当正则化长细比较大时，计算得到的柱子曲线低于规范曲线.当正则化长细比较小时，计算得到的柱子曲线与美国规范曲线间的差距较小，当正则化长细比较大时，计算得到的柱子曲线与欧洲不锈钢规范的柱子曲线间的差距较小.
此外，对于奥氏体型不锈钢构件，模拟得到的整体稳定系数与欧洲规范曲线差值较小，对于双相体型不锈钢构件，模拟得到的整体稳定系数与欧洲规范曲线差值较大，表明不同材料对构件整体稳定承载力影响较大，因此两者应分别进行讨论.构件绕强轴失稳时模拟得到的整体稳定系数与欧洲规范曲线差值小于构件绕弱轴失稳时模拟得到的整体稳定系数与欧洲规范曲线差值，因此应分别讨论工字形构件绕强轴失稳和绕弱轴失稳2种情况.
本文提出的三段式计算方法适用于奥氏体型和双相体型2类不锈钢构件.由图14可知，模拟值明显高于欧洲不锈钢规范曲线，表明欧洲不锈钢规范较保守；模拟值与美国不锈钢柱子曲线相差较大，表明美国不锈钢规范不适用于焊接不锈钢构件.三段式与试验和有限元数据点的分布形态吻合较好，对工字形截面构件的承载力计算值与试验和有限元的平均比值均近似为1，表明建议公式能够对奥氏体型、双相体型不锈钢的工字形截面构件的整体稳定承载力进行合理的计算.
5 结 论
1）通过考虑了不锈钢材料力学性能、焊接残余应力、整体几何初始缺陷等因素的有限元模型对试验进行验证，表明了有限元模型的可靠性和适用性.
2）通过参数分析可确定不锈钢构件的材料特性、几何初始缺陷与长细比对构件整体稳定承载力影响较大，截面残余应力对构件整体稳定承载力影响较小，截面宽厚比对构件整体稳定承载力影响可忽略.
3）欧洲不锈钢规范对于构件整体稳定性能的预测较保守，美国不锈钢规范对于焊接构件可能不适用.
4）计算了742个焊接工字形不锈钢构件的数值算例，根据算例结果提出了三段式计算方法，建议公式可很好地预测不锈钢构件的整体稳定承载力，可以为工程设计应用和相关设计规范编制提供参考.
参考文献
[1] 王元清， 袁焕鑫， 石永久，等. 不锈钢结构及其应用和研究现状[J]. 钢结构， 2010（2）： 1-13.
WANG Yuan-qing， YUAN Huan-xin， SHI Yong-jiu， et al. A review of current applications and research of stainless steel structure [J]. Steel Structure， 2010（2）： 1-13.（In Chinese）
[2] MIRAMBELL E， REAL E. On the calculation of deflections in structural stainless steel beams： an experimental and numerical investigation [J]. Journal of Constructional Steel Research， 2000， 54（1）： 109-133.
[3] QUACH W M， TENG J G， CHUNG K F.Three-stage full-range stress-strain model for stainless steels[J]. Journal of Structural Engineering，2008，134（9）： 1518-1527.
[4] GARDNER L， CRUISE R B. Modeling of residual stresses instructural stainless steel sections[J]. Journal of Structural Engineering， 2009， 135（1） ： 42-53. [5] QUACH W M， TENG J G， CHUNG K F. Residual stresses in press-braked stainless steel sections I： coiling and uncoiling of sheets[J]. Journal of Constructional Steel Research， 2009， 65（8）：1803-1815.
[6] 王元清，高博，戴国欣，等.焊接不锈钢工字形截面受弯构件整体稳定与设计方法[J].土木建筑与环境工程， 2013， 35（1）：57-62.
WANG Yuan-qing， GAO Bo， DAI Guo-xin，et al.Analysis on overall buckling behavior of stainless steel welded I-beams and design methods[J]. Journal of Civil Architectural &Environmental Engineering， 2013， 35（1）：57-62.（In Chinese）
[7] YANG Lu， WANG Yuan-qing， GAO Bo，et al.Two calculation methods for buckling reduction factors of stainless steel welded I-section beams[J]. Thin-Walled Structures， 2014， 83： 128-136.
[8] YUAN H X， WANG Y Q， GARDNER L，et al.Local-overall interactive buckling of welded stainless steel box section compression members [J]. Engineering Structures， 2014， 67（5）：62-76.
[9] 舒赣平， 郑宝锋， 沈晓明. 不锈钢压弯构件平面内稳定承载力计算方法研究[J].工业建筑， 2012， 42（5）： 41-44.
SHU Gan-ping， ZHENG Bao-feng，SHEN Xiao-ming. In-plane stability design method of stainless steel beam-columns[J]. Industrial Construction， 2012，42（5）：41-44.（In Chinese）
[10]舒赣平， 郑宝锋， 沈晓明. 不锈钢轴心受压构件稳定承载能力计算方法研究[J].工业建筑， 2012， 42（5）：21-28.
SHU Gan-ping， ZHENG Bao-feng，SHEN Xiao-ming. New design method for stainless steel cold-formed tubular columns under axial load [J]. Industrial Construction， 2012，42（5）：21-28. （In Chinese）
[11]GARDNER L， NETHERCOT D A. Experiments on stainless steel hollow sections-part 1： material and cross-sectional behavior [J]. Journal of Constructional Steel Research， 2004， 60（9）：1291-1318.
[12]GARDNER L， NETHERCOT D A. Experiments on stainless steel hollow sections-part 2： member behavior of columns and beams [J]. Journal of Constructional Steel Research， 204， 60（9）： 1319-1332.
[13]GARDNER L. The continuous strength method [J]. Structures & Buildings， 2008， 161（S3）： 127-133.
[14]GARDNER L， NETHERCOT D A.Experiments on stainless steel hollow sections―part 2： member behavior of columns and beams[J]. Journal of Constructional Steel Research， 2004， 60（9）：1319-1332.
[15]GARDNER L，NETHERCOT D A.Experiments on stainless steel hollow sections―part 1： material and cross-sectional behaviour[J]. Journal of Constructional Steel Research， 2004， 60（9）：1291-1318.
[16]班慧勇， 施刚， 石永久. Q420 高强度等边角钢轴压构件整体稳定性能设计方法研究[J]. 工程力学，2014， 31（3）：63-70.
BAN Hui-yong， SHI Gang， SHI Yong-jiu. Investigation on design method of overall buckling behaviour for Q420 high strength steel equal-leg angle members under axial compression [J]. Engineering Mechanics，2014， 31（3）：63-70. （In Chinese）
[17]袁焕鑫. 焊接不锈钢轴心受压构件局部稳定和相关稳定性能研究[D].北京：清华大学土木工程系，2014.
YUAN Huan-xin.Local and local-overall buckling behaviour of welded stainless steel member under axial compression[D]. Beijing：Department of Civil Engineering，Tsinghua University， 2014. （In Chinese）
[18]EN 1993-1-4 Eurocode 3： design of steel structures，part 1-4： general rules-supplementary rules for stainless steels [S]. London：BSI，2006.
[19]ASCE 8-02 Specification for the design of cold-formed stainless steel structural members [S]. Virginia ： American Society of Civil Engineers， 2001.