考虑SD效应及初始缺陷影响的钢压杆弹塑性稳定分析
李斌1,2, 李传习2, 刘雪松2, 张玉平2
1.湖南理工学院 土木建筑工程学院,湖南 岳阳 414006
2.长沙理工大学 土木与建筑学院,长沙 410076

第一作者:李斌(1981—),男,浙江绍兴人,副教授,博士生.研究方向为钢结构设计与理论.email: Leebin2009@126.com.

摘要

考虑钢压杆一般存在初变形和初偏心等初始缺陷及钢材拉压强度的差异,基于Ježek法推导在轴向偏心荷载作用下钢压杆弹塑性失稳破坏时的极限荷载公式.通过探讨SD效应、长细比及初始缺陷等因素对钢压杆弹塑性失稳破坏极限荷载的影响规律,拟合得到了两端简支的矩形截面钢压杆弹塑性极限荷载近似计算公式;所得解析理论与近似计算公式可以蜕化求解带初弯曲的轴压构件、带初偏心的压杆、理想的中心受压杆及压弯构件考虑SD效应时的弹塑性极限荷载.分析结果表明:考虑SD效应时钢压杆的弹塑性稳定极限荷载提高10.5%~12.2%.因此,在实际工程中考虑SD效应可以使钢压杆的强度得以充分发挥.

关键词: 钢压杆; 极限荷载; 弹塑性稳定; SD效应; 初始缺陷; 公式拟合
中图分类号:TU391 文献标志码:A 文章编号:1673-0291(2017)03-0090-06
Elastic-plastic stability analysis of steel compression members considering SD effect and initial defects
LI Bin1,2, LI Chuanxi2, LIU Xuesong2, ZHANG Yuping2
1.College of Civil Engineering & Architectural, Hunan Institute of Science and Technology, Yueyang Hunan 414006,China;
2.School of Civil Engineering and Architectural, Changsha University of Science & Technology, Changsha 410076,China;
Abstract

Considering the initial imperfections of steel compression members,such as initial deformation and initial eccentricity, and the difference of tension and compression strength of steel, the formula for calculating ultimate load of steel compression members under elastic-plastic buckling failure caused by axial eccentric load is derived by the Ježek method. The approximate calculation formula of elastic-plastic ultimate load of a rectangular section steel compression member with two ends simply supported is obtained through studying the influences of strength difference(SD)effect, slenderness and initial defects on ultimate load of elastic-plastic buckling failure of steel compression members. The analytical solution and approximate calculation formula in this paper can be used to solve the elastic-plastic limit load of the steel compression members with initial bending, steel compression members with initial eccentric, ideal center compression members or beam-columns considering SD effect. The results show that the elastic-plastic ultimate load of steel compression members increases by 10.5%~12.2% when SD effect is considered. Therefore, the strength of steel compression members can be fully played when SD effect is considered in practical engineering.

Keyword: steel compression member; ultimate load; elastic-plastic stability; SD effect; initial defect; formular fitting

钢结构受压杆件和受压柱被广泛应用于土木工程建设的各个领域, 如桁架桥的弦杆、网架中的杆件、工业厂房及海上钻井平台的立柱等, 这些杆件由于长细比较大, 容易引起弹性或弹塑性失稳.其中, 以小挠度弯曲近似理论为基础的Euler公式, 很好地解决了各种边界条件下等截面理想压杆的弹性稳定问题[1], 并得到了广泛的工程应用[2].基于梁段挠曲线近似微分方程解析解的传递矩阵法[3], 基于变截面压杆与等截面压杆两者弯曲应变能、顶点位移等效的等效刚度法[4], 基于压杆临界荷载、失稳模态与轴向振动主频率、主模态一致, 进行Ritz展开的模态摄动法[5]及基于几何非线性有限元理论的数值分析法[6]等方法也为变截面压杆弹性稳定临界荷载的求解提供了多种有效途径.而切线模量理论、割线模量理论、Jež ek法[7, 8, 9]和Shanley模型[10]等近似方法和力学分析模型的提出及应用使得压杆弹塑性失稳时极限荷载的计算分析得以开展和简化.其中, 基于失稳截面静力平衡的Jež ek法具有力学概念清晰和简便易行的优点, 将其用于求解受压、压弯等受力构件的弹塑性稳定性可以得到精确度较高的近似理论解.文献[7, 8, 9]利用Jež ek法计算矩形、工字型及空心方管等截面压弯构件碳纤维片加固前后弹塑性失稳破坏时的极限荷载, 并考虑了初弯曲[9]的影响.尹一平等[11]根据Jež ek的思想结合中国铁路桥梁钢结构设计规范计算了无初始残余应力箱形截面杆件压溃荷载的近似计算公式.这些研究成果证明了该方法的简便性、有效性和精确性.

但是文献[12, 13]等研究成果表明, 塑性材料的压缩和拉伸屈服极限存在明显差异, 对结构的应力和应变等力学性能有一定影响.压缩与拉伸屈服极限之比γ 一般为1~4/3, 这种强度差异简称SD(strength difference)效应.文献[14]考虑SD效应利用对称性对两端固支超静定梁的弹塑性全过程进行了分析, 获得了极限荷载.文献[15]分析了拉压强度不同的厚壁球壳的加载应力、残余应力及安定极限.文献[16]对内外压力作用下的厚壁筒进行了应力分析, 证明了材料拉压性能的差异对厚壁筒的轴向应力和环向应力有一定的影响, 拉压强度不同的材料不能简化成无体积变化的材料.上述研究主要针对SD效应对结构强度及刚度的影响进行分析, 对弹塑性稳定性能影响的研究及应用目前鲜有报道.

本文作者基于Jež ek法, 建立了考虑SD效应影响时失稳截面受压区单侧塑性屈服、受压受拉双侧塑性屈服2种失稳模式, 利用静力平衡条件、变形协调条件及极值条件推导弹塑性失稳极限荷载.在对SD效应、长细比、初始缺陷、屈服强度等参数影响规律分析的基础上, 实现了两端简支矩形截面钢压杆弹塑性极限荷载近似计算公式的拟合.

1 考虑SD效应的弹塑性极限荷载
1.1 基本假定

具有初始缺陷的理想弹塑性矩形截面钢压杆, 两端铰支, 长为 l.杆件初始变形按照一阶屈曲模态假设, 即初始几何变形曲线 y0=v0sin(πx/l), v0为初始变形.偏心荷载F, 偏心距e作用下的变形曲线为 y1=v1sin(πx/l), 如图1所示.

图1 偏心钢压杆变形曲线Fig.1 Deformation curves of an eccentric steel compression member

1.2 受压侧出现塑性区

轴向偏心荷载作用下, 图1所示钢压杆弹塑性失稳时截面的应变沿截面高度分布符合平截面假定, 应力、应变分布情况如图2所示.

图2 受压区单侧屈服时的应变和应力Fig.2 Strain and stress of the section yielded in compression

图中 bh分别为截面的高和宽; 阴影部分表示塑性受压区, εyfy分别为构件受拉屈服时的应变和屈服强度, fy=Eεy; εtσt分别为截面受拉区下缘的应变和应力, σt=Eεt; E为材料的弹性模量; γ 为压缩与拉伸屈服强度之比; φ 为失稳截面的曲率; 为截面的弹性核高度.

由静力平衡条件可得失稳截面x=l/2处轴力F和弯矩M分别为

F=bhγfy-bd2(γfy+σt)(1)

M=bd(γfy+σt)(h4-d6)=F(e+v0+v1)(2)

由截面的曲率

-y″(x)=v1k(x)E=γfy+σtEd(3)

解得:

σt=k(x)v1d-γfy(4)

式中: k(x)=Eπ2l2sin(πxl)=k0sin(πxl).

将式(4)分别代入式(2)和式(1), 且 x=l/2, 得:

F(e+v0+v1)=-bk0v1d36+hbk0v1d24(5)

d=2(bhfy-F)bk(x)v1(6)

式(5)两边对 v1求导, 并注意极值条件 dFdv1=0, 可得:

F=bk0d22(-d3+h2)+bk0v1d2(h-d)dddv1(7)

联立式(5)、式(6)和式(7)可求得钢压杆中点截面处单侧出现塑性区时弹塑性失稳的极限荷载 Fu.

1.3 受压、受拉双侧出现塑性区

当截面下表面应力 σt逐渐增大到与 fy相等时, 截面受拉侧也将达到屈服极限, 即受压与受拉侧均出现塑性区.此时, x=l/2截面处的应变和应力示意图如图3所示.其中, c表示截面受拉侧塑性区高度.

图3 双侧屈服时的应变和应力Fig.3 Strain and stress of the section yielded in compression and tension

同理, 由静力平衡条件可得:

F=γfybh-bd2(γ+1)fy+(γ+1)fybc(8)

M=F(e+v0+v1)=bd2(γ+1)fy(h2-d3-c)+bc(γ+1)fy(h2-c2)(9)

由式(3)可知, 此时失稳截面的曲率满足条件

v1k0E=γfy+fyEd(10)

d=(γ+1)fyk0v1(11)

将式(11)代入式(8)中得:

c=a-d2(12)

其中 a=γfybh-F(γ+1)fyb.

将式(11)和式(12)代入式(9)可得:

F(e+v0+v1)=(γ+1)bfy(-124d2+ha2-a22)(13)

式(13)两边对 v1求导, 注意极值条件 dFdv1=0, 可得:

F=-(γ+1)dbfy12dddv1(14)

由式(11)得:

dddv1=-(γ+1)fyk(x)v12(15)

将式(11)和式(15)带入式(14)中可得:

F=(γ+1)3bfy312v13k02(16)

联立式(11)、式(12)、式(13)和式(16)可以求得钢压杆中点截面处双侧出现塑性区时弹塑性失稳的极限荷载值 Fu.

1.4 钢压杆的弹塑性极限荷载

计算钢压杆弹塑性极限荷载时, 首先假设钢压杆单侧出现塑性区, 联立式(5)、式(6)和式(7)求得钢压杆中点截面处单侧出现塑性区时弹塑性失稳的极限荷载 Fu和挠度v1.

Fuv1回代式(4)和式(6), 若受拉区下边缘应力满足

σt=k(x)v1d-γfyfy(17)

k(x)v1d/fy1+γ(18)

则钢压杆仅受压侧出现塑性区.反之, 钢压杆受拉侧已经进入塑性状态, 应联立式(11)~式(13)及式(16)求解钢压杆弹塑性失稳的极限荷载 Fu.

2 参数分析与公式拟合

某理想弹塑性钢压杆长l=8m, 横截面宽b=0.1m, 高h=0.2m, 偏心距为e, 初始变形最大挠度为v0, 弹性模量E=206GPa.为表述方便, 令p1=v0/l× 103, p2=(e+v0)/l× 103.

2.1 极限荷载影响参数分析

由解析理论可知, 初偏心距e与最大初始变形v0 的大小对钢压杆弹塑性稳定极限荷载的影响作用一致.在下面的分析中, 取偏心距e=0.038 m, 分别计算不同最大初始变形形v0及压缩与拉伸屈服极限之比γ 的钢压杆弹塑性稳定极限荷载、中点处挠度、弹性核高度、受拉侧出现塑性区分界点及失稳截面塑性区分布情况, 令 p1=v0/l×103.计算结果分别列于表1~表4.

表1 极限荷载Fu的计算结果 Tab.1 Calculation results of ultimate load Fu kN
表2 中点挠度v1的计算结果 Tab.2 Calculation results of deflection v1 at the midpoint mm
表3 弹性核高度d的计算结果 Tab.3 Calculation results of hight d of the elastic core mm

表1~表3可知, 随着初始变形 v0增大, 初始缺陷引起的弯矩对截面内力的贡献增大, 使得失稳截面更容易进入双侧塑性状态.导致极限荷载值和失稳截面的弹性核高度减小, 中点挠度增大, 这与材料的工作性能相符合.

表4中, γ从1增加到1.3的过程中, 失稳截面由单侧塑性状态进入双侧塑性状态的分界点分别为5、4.125、3.375、2.625、1.875、1.25、0.625.由表1~表4分析可得, 考虑SD效应的影响时, 受压区的屈服强度提高, 相对 γ=1.0而言, 受拉区进入塑性状态的分界点提前, 说明材料更容易进入双侧塑性状态.此时, 构件极限荷载提高10.5%~12.2%, 中点挠度增大10.8%~28.6%, 相应的弹性核高度的增幅稳定在3.8%左右, 可见材料潜能得以进一步发挥.

表4 失稳截面塑性区分布情况 Tab.4 Distribution of plastic zone in instability section

图4~图7分别表示压缩与拉伸屈服极限之比γ 、长细比λ 、初始缺陷 e+v0及拉伸屈服极限 fy对构件极限荷载值的影响.构件的极限荷载值随γ 增大呈线性增大; 随初始缺陷及长细比λ 的增大都呈幂指数级大幅减小; 随拉伸屈服极限fy的增大呈幂指数增大.

图4 γ 值对极限荷载的影响Fig.4 Influence of γ on ultimate load

图5 初始缺陷对极限荷载的影响Fig.5 Influence of initial defects on ultimate load

图6 长细比对极限荷载的影响Fig.6 Influence of slenderness ratios on ultimate load

图7 fy对极限荷载的影响Fig.7 Influence of fy on ultimate load

2.2 极限荷载公式拟合

基于上述各种参数的影响规律分析, 对长细比、压缩与拉伸屈服极限之比、初始缺陷及拉伸屈服强度等进行多元回归, 拟合得到两端简支矩形截面偏心钢压杆极限荷载 Fu的近似计算公式

Fu=A×10-3(γ+B)λC(e+v0l×10-3+1)DEbhfy235(19)

参数 ABCD可由表5确定.

2.3 计算理论及近似公式的参数分析

1)设ev0e+v0等于0, 本文计算理论及近似公式可以分别求解带初弯曲的轴压构件、带初偏心的压杆及理想的中心受压杆考虑SD效应时的弹塑性极限荷载.

2)设极限荷载作用下, 由初偏心距和初始变形引起的弯矩 F(e+v0)=M0, 即令式(5)、式(13)或式(19)中 e+v0=M0/F, 则本文计算理论及近似公式同样适用于考虑SD效应时压弯构件弹塑性失稳极限荷载的计算.

以文献[7]某承受初始弯矩M0的矩形压弯钢柱为例, 计算弹塑性失稳时的极限荷载, 计算结果列于表6.由表6可见, 当不计SD效应影响时, 本文解析理论计算结果与文献[7]完全一致; 近似公式精度较高, 计算得到的极限荷载与解析解相比最大仅差2.82%.因此, 可以验证本文结果的正确性.

表5 参数ABCD取值 Tab.5 Values of parameter A, B, C and D
表6 极限荷载Fu的计算结果 Tab.6 Calculation results of ultimate load Fu kN
3 结论

1)考虑SD效应对提高钢压杆的极限承载力有着明显的作用.随着压缩与拉伸屈服极限之比γ 值的增大, 钢压杆双侧相对提前出现塑性区, 极限荷载值呈线性增长.因此, 在工程实际中考虑SD效应有助于充分发挥材料的强度.

2)初偏心距与最大初始变形的大小对钢压杆工作性能的影响作用一致.两种初始缺陷和长细比对钢压杆的极限承载力的影响很大.初始缺陷的增大使得构件的极限承载力和弹性核高度明显的下降, 构件中点的挠度值明显增大.另外, 构件的极限荷载值随着钢压杆长细比λ 值的增大急剧下降, 随着拉伸屈服极限 fy的增大小幅度增大.

3)本文解析理论及近似公式具有广泛的适用性, 可以求解带初弯曲的轴压构件、带初偏心的压杆、理想的中心受压杆及压弯构件考虑SD效应时的弹塑性极限荷载.

The authors have declared that no competing interests exist.

参考文献
[1] TIMOSHENKO S P, GERE J M. Theory of elastic stability[M]. Tokyo: McGrawHill-Kogakusha Ltd, 1961: 46-100. [本文引用:1]
[2] 钢结构设计规范: GB 50017—2015[S][S]. 北京: 中国计划出版社, 2015.
Code for design of steel structures: GB 50017—2015[S]. Beijing: China Planning Press, 2015. (in Chinese) [本文引用:1]
[3] 刘庆潭. 含锥形变截面压杆稳定计算的传递矩阵法[J]. 计算结构力学及应用, 1996, 13(3): 364-368.
LIU Qingtan. The method of transfer matrix of calculation on stability of compressed bar with variable section[J]. Computational Structural Mechanics and Applications, 1996, 13(3): 364-368. (in Chinese) [本文引用:1]
[4] 吴亚平. 变截面压杆稳定计算的等效刚度法[J]. 力学与实践, 1994, 16(1): 58-60.
WU Yaping. The equivalent stiffness method for stability of variable cross-section compression bars[J]. Mechanics in Engineering, 1994, 16(1): 58-60. (in Chinese) [本文引用:1]
[5] 楼梦麟, 李建元. 变截面压杆稳定问题半解析解[J]. 同济大学学报(自然科学版), 2004, 32(7): 857-860.
LOU Menglin, LI Jianyuan. Semi-analytical approach for stability of variable cross-section compression bars[J]. Journal of Tongji University (Natural Science), 2004, 32(7): 857-860. (in Chinese) [本文引用:1]
[6] LI Guoqiang, LI Jinjun. Effects of shear deformation on the effective length of tapered columns with I-section for steel portal frames[J]. Structural Engineering and Mechanics, 2000, 10(5): 479-489. [本文引用:1]
[7] 杨刚, 王海滨, 张爱锋, . 碳纤维增强压弯构件弹塑性失稳解析解[J]. 工程力学, 2007, 24(5): 53-57.
YANG Gang, WANG Haibin, ZHANG Aifeng, et al. Analytical solution of elastic-plastic buckling of member wrapped by carbon fibre under the axial compressive load and bending moment[J]. Engineering Mechanics, 2007, 24(5): 53-57. (in Chinese) [本文引用:4]
[8] NIU Peng, YANG Gang, GONG Benqi, et al. Analytical solution of elastic-plastic buckling of a H-shaped steel column under axial compressive load and bending moment[J]. Key Engineering Material, 2010, 417/418: 857-860. [本文引用:2]
[9] 钮鹏, 金春福. 几何缺陷影响下的CFRP-方钢管极限承载力解析解[J]. 工程力学, 2015, 32(增1): 322-326.
NIU Peng, JIN Chunfu. Analytical solutions on ultimate bearing capacity of a square CFRP-steel tube member with initial imperfection[J]. Engineering Mechanics, 2015, 32(S1): 322-326. (in Chinese) [本文引用:3]
[10] 李国强, 陈琛. 火灾下轴向约束钢柱性能的Shanley理论模拟[J]. 同济大学学报(自然科学版), 2013, 41(4): 491-495.
LI Guoqiang, CHEN Chen. Theoretical simulation of axial restrained steel column in fire with Shanley model[J]. Journal of Tongji University(Natural Science), 2013, 41(4): 491-495. (in Chinese) [本文引用:1]
[11] 尹一平, 方秦汉. 箱形杆件压溃荷载近似计算[J]. 华中科技大学学报(城市科学版), 2008, 25(3): 128-131.
YIN Yiping, FANG Qinhan. Approximate calculation on ultimate carrying capacity of column with box section[J]. Journal of Huazhong University of Science and Technology (Urban Science Edition), 2008, 25(3): 128-131. (in Chinese) [本文引用:1]
[12] CASEY J, JAHEDMOTLAGH H. The strength-differential effects in plasticity[J]. International Journal of Solids and Structures, 1984, 20(4): 377-393. [本文引用:1]
[13] LIANG Yaping, WANG Huizhen, REN Xingmin. Elastic-plastic limit analysis ofcombination cylinders made of SD effect materials based on unified strength theory[J]. Journal of University of Science and Technology of China, 2008, 38(4): 364-368. [本文引用:1]
[14] 李斌, 韦成龙, 李传习. 集中荷载作用下两端固支梁考虑SD效应的极限荷载分析[J]. 重庆交通大学学报(自然科学版), 2014, 33(3): 6-10.
LI Bin, WEI Chenglong, LI Chuanxi. Limit load of the beam fixed at two ends under concentrated load considering SD effect[J]. Journal of Chongqing Jiaotong University (Natural Science), 2014, 33(3): 6-10. (in Chinese) [本文引用:1]
[15] 徐栓强, 俞茂宏. 拉压强度不同材料厚壁球壳的安定性分析[J]. 机械设计与制造, 2005(1): 36-37.
XU Shuanqiang, YU Maohong. Shakedown analysis of thick spherical shell of material with different strength in tension and compression[J]. Machinery Design & Manufacture, 2005(1): 36-37. (in Chinese) [本文引用:1]
[16] 阮澍铭. 拉压性能不同材料全量型本构关系及厚壁筒的应力分析[J]. 力学季刊, 2003, 24(3): 423-427.
RUAN Shuming. Constitutive equation of deformation theory for materials with different behaviors in tension and compression and stress analysis of thick-walled tube[J]. Chinese Quarterly of Mechanics, 2003, 24(3): 423-427. (in Chinese) [本文引用:1]