图1.2 轴心压杆的平衡状态2
用上述理想轴心压杆的情况来描述钢结构的失稳现象是不够的,钢结构的失稳现象就其性质而言,可以分为三类稳定问题。
1.2.1 分支点失稳
理想的(即无缺陷的、笔直的)轴心受压杆件和理想的中面内受压的平板的失稳(屈曲)都属于分支点失稳。也称平衡分岔失稳,或称第一类失稳。
图1.3a为一理想轴心受压构件,当轴向压力P< Pcr时,压杆沿轴向只被压缩Δc,杆始终处于直线平衡状态,称为原始平衡状态。此时如果在其横向施加微小干扰,杆件会呈微弯曲状态而偏离原平衡位置,但是撤去此干扰后,压杆立即恢复到原直线平衡状态。可见,原始平衡状态具有唯一的平衡形式。
当P= Pcr时,压杆会突然弯曲,该现象称为丧失稳定,或称为屈曲。如图1.3b所示,构件由原来挺直的平衡状态转变到微弯曲的平衡状态。从图1.3c表示的荷载(P)—位移(δ)曲线中可以看出,当荷载到达A点后,杆件可能有两个平衡路径,即直线AC和水平线AB(AB’), A点称为两个平衡路径的分支点,或分岔点。由于在同一个荷载点出现了平衡分支现象,所以将此种失稳现象称为分支点失稳。
(a)原始平衡 (b)临界平衡 (c)P—δ曲线
图1.3 理想轴心受压构件
分支点失稳又可以分为稳定分支点失稳和不稳定分支点失稳两种。 1. 稳定分支点失稳
图1.3c所示荷载—位移曲线是根据小挠度理论分析得到的,如按大挠度理论分析,轴心受压构件屈曲后,荷载随横向位移加大而略有增加,但横向位移的增长速度远大于轴向力的提高速度,如图1.4b所示。轴心压杆屈曲后,荷载—位移曲线是AB或AB’,这种平衡状态是稳定的,属于稳定分支点失稳。由于压杆因弯曲变形而产生弯矩,在压力和弯矩的共同作用下,杆件最大弯矩作用截面边缘纤维先屈服,随着塑性发展,压杆很快就达到承载能力极限状态,即极限荷载Pu与屈曲荷载Pcr相差很小,因此,轴心受压构件屈曲后强度并不能被利用。
对图1.5a所示四边有支撑的薄板,当中面均匀压力P达到屈曲荷载Pcr后,板发生凸曲,同时在板中面产生横向薄膜拉应力,牵制了板的变形,使板屈曲后仍能承受较大的荷载增量,屈曲后板仍处于稳定平衡状态,该板的失稳属于稳定分支点失稳。薄板屈曲后荷载—位移曲线如图1.5b中的AB或AB’所示,由于薄板的极限荷载Pu远超过屈曲荷载Pcr,所以可以利用板屈曲后的强度。
3
(a)轴心受压构件 (b)P—δ曲线
图1.4 大挠度弹性理论分析的荷载—位移关系
(a)中面均匀受压的四边支承薄板 (b)P—w曲线
图1.5 中面均匀受压的四边支承薄板的荷载—位移关系
2. 不稳定分支点失稳
如果结构或构件发生分支点失稳后,只能在远比临界荷载低的条件下维持平衡状态,则称此类失稳为不稳定分支点失稳。图1.6a所示承受均匀压力的圆柱壳的失稳就是不稳定分支点失稳,荷载—位移曲线如图1.6b中的OAB或OAB’所示。
(a)均匀受压圆柱壳 (b)荷载—位移曲线
图1.6 不稳定分支点失稳
4
1.2.2 极值点失稳
图1.7a所示偏心受压构件,作用力P的偏心距为e,其失稳过程的压力(P)—挠度(Δ)曲线见图1.7b。随着压力P的增加,偏心压杆的挠度Δ也随之增长,形成曲线的上升段OA,压弯构件处于稳定平衡状态;但是到达曲线的最高点A时,构件的抵抗力开始小于外力作用,即A点为压弯构件承载力的极限点,表示压弯构件开始丧失整体稳定,Pu为偏心压杆的最大承载力,也称为偏心压杆的极限荷载或压溃荷载;A点之后出现了曲线的下降段AB,为了维持构件的平衡状态必须不断降低端部压力P,构件处于不稳定平衡状态。从压弯构件的失稳过程可知,其荷载—位移曲线只有极值点,没有出现由直线平衡状态向弯曲平衡状态过渡的分岔点,构件弯曲变形的性质始终不变,称这种失稳为极值点失稳,也称为第二类失稳。
(a)偏心受压构件 (b)荷载(P)—挠度(δ)曲线
图1.7 极值点失稳
1.2.3 跃越失稳
对两端铰接的坦拱结构(图1.8a),在均布荷载q作用下产生挠度w,其荷载—挠度曲线(图1.8b)也有稳定的上升段OA,但是到达曲线的最高点A时会突然跳跃到一个非临近的具有很大变形的C点,即由向上拱起的位形突然跳到下垂的位形,与A点对应的荷载qcr为坦拱的临界荷载;下降段AB不稳定,BC段虽然稳定上升,但是因为结构已经破坏而不能被利用。这种结构由一个平衡位形突然跳到另一个非临近的平衡位形的失稳现象称为跃越失稳。跃越失稳既无平衡分支点,又无极值点,但与不稳定分支失稳又有相似之处,都在丧失稳定平衡后经历一段不稳定平衡,然后达到另一个稳定平衡状态。钢结构油罐、扁球壳顶盖等的失稳也属此种类型。
qw (a)均布荷载作用下的坦拱 (b)荷载—挠度曲线
图1.8 跃越失稳
5
1.3 临界力的计算方法
结构由稳定平衡到不稳定平衡的界限状态称为临界状态。结构处于临界状态时的荷载值称为临界荷载值,稳定计算的主要目的在于确定临界荷载值。求临界荷载值的方法很多,可分为精确计算方法和近似计算方法两大类,其中静力法、能量法分别是两类方法中常用的计算方法。
1.3.1 静力法
静力法即静力平衡法,也称中性平衡法,此法是求解临界荷载的最基本方法。对第一类弹性稳定问题,在分支点存在两个临近的平衡状态:原始直线平衡状态和产生了微小弯曲变形的平衡状态。静力法就是根据已发生了微小弯曲变形后结构的受力条件建立平衡微分方程,而后解出临界荷载。下面以图1.9a所示两端铰接轴心受压直杆说明静力法的原理和计算步骤。
当荷载达到临界荷载(P= Pcr)时,压杆会突然弯曲,由原来的直线平衡状态转变到图1.9a中实线表示的微弯的曲线平衡状态。此时杆件除弯曲外,还受压缩及剪切作用,由于压缩和剪切的影响很小,一般忽略不计,则任一截面(图1.9b)内力矩与外力矩的平衡关系为
MPy (1.1)
由挠曲线的近似微分方程
EIyM (1.2)
可得
EIyPy0 (1.3)
式中:E为材料弹性模量,I为杆件截面惯性矩。令2P,式(1.3)为一常系数微分方程
(EI)y2y0 (1.4)
其通解为
yC1sinxC2cosx (1.5)
当两端铰接时,边界条件为
(1.6)
xl时,y0将边界条件代入式(1.5),得如下齐次方程组
C10C210 (1.7) CsinlCcosl021当C1C20时,满足式(1.7),但由式(1.5)知,此时y0,表示杆件处于直线平衡状态,与图1.9b不符。对应杆件曲线平衡状态,要求y0,即C1、C2有非零解,为此要求方程组(1.7)的系数行列式必须等于零,即
D()0sinl10 (1.8)
coslx0时,y0上式D()0为稳定特征方程,解之得
sinl0 (1.9)
则有
ln (n=0,1,2, ┄) (1.10)
即
Pn22EIl2 (1.11)
当n=1时,得到P的最小值Pcr,即分支屈曲荷载,又称欧拉(Euler)临界荷载
Pcr2EI/l2 (1.12)
6
(a)轴心受压 (b)任一截面平衡关系
图1.9 两端铰接轴心受压构件
由上述可见,静力法求临界荷载首先假定杆件已处于新的平衡状态,并据此列出平衡微分方程,然后解此方程并结合边界条件得到一组与未知常数数目相等的齐次方程;对于新的平衡形式要求齐次方程组的系数行列式必须等于零,即D0,从而解出临界力Pcr。稳定特征方程D0通常简称为稳定方程。
1.3.2 能量法
静力法通过建立轴心受压构件微弯状态时的平衡方程求出临界荷载的精确解,但是对于有些轴心受压构件,如变截面的或者压力沿轴线变化的构件,静力法得到的是变系数微分方程,求解十分困难,有时甚至无法求解,这时就需要采用其它方法,如近似计算方法中的能量法求解。能量法已广泛应用于轴心受压构件、压弯构件、受弯构件和板壳结构的稳定计算。
用能量法求解临界荷载的途径主要有能量守恒原理和势能驻值原理。
1. 能量守恒原理求解临界荷载
用能量守恒原理解决结构弹性稳定问题的方法是铁摩辛柯(Timoshenko)首先提出的,故又称为铁摩辛柯能量法[10]。保守体系处在平衡状态时,贮存在结构体系中的应变能等于外力所做的功,此即能量守恒原理。
当作用着外力的弹性结构偏离原始平衡位置而产生新的微小位移时,如果应变能的增量U大于外力功的增量W,即此结构具有恢复到原始平衡位置的能力,则结构处于稳定平衡状态;如果UW,则结构处于不稳定平衡状态而导致失稳;临界状态的能量关系为
UW (1.13)
式(1.13)是铁摩辛柯能量法计算临界力的基本方程。
仍以图1.9a所示两端铰接轴心受压直杆说明能量守恒原理求解临界力的具体过程。当轴向力P=Pcr时,压杆发生横向挠曲,杆件中产生弯曲应变能增量
1lM2U0dx (1.14)
2EI以MEIy代入后,有
1lU0EIy2dx (1.15)
2由直线平衡状态过渡到曲线平衡状态过程中,外荷载P所做的功为
WP (1.16)
7
式中,是力P作用点下降的距离。
在压杆上任取一微段dx,变形后与轴x的夹角为θ,微段dx弯曲前后在轴x上投影的长度差为
dx1cos2dxsin2因杆件处于微弯状态,θ角很小,故有tg22 (1.17)
dyy,则可推得 dx21122dxsin2dxtg2dxtgydx (1.18)
22222故有
01cosdx则
l1l2ydx (1.19) 02Pl2WP0ydx (1.20)
2则式(1.13)可表达为
ll220EIydxPcr0ydx (1.21)
求出临界力为
Pcrl2EIydx0l2ydx0 (1.22)
式中y(x)是满足位移边界条件的任一可能曲线位移方程。具体应用将在第二章详细说明。
2. 势能驻值原理求解临界荷载
势能驻值原理指:受外力作用的结构,当位移有微小变化而总势能不变,即总势能有驻值时,结构处于平衡状态。其表达式
0 (1.23)
式中为结构总势能的一阶变分,有
UVUW (1.24)
其中U是虚位移引起的结构内应变能的变化,它总是正值;W表示外力因虚位移而作的功,且外力势能的变化V等于外力虚功的负值,即VW。这样,势能驻值原理还可以表述为:弹性变形体对每一个和约束相容的虚位移,其总势能的一阶变分为零,则该体系处于平衡状态。
势能驻值原理与平衡方程是等价的,用该原理可以解决复杂结构的弹性稳定问题。如很多结构很难直接建立平衡方程,则可以先写出结构总势能,然后利用0,即可得到平衡方程。还可以先假定构件挠曲线形状,给出挠曲线方程,将其代入总势能,通过0解出临界荷载。若给出的挠曲线方程满足几何边界条件,称求解临界荷载的方法为里兹(Ritz)法 [11];若给出的挠曲线方程不仅满足几何边界条件,而且满足自然边界条件,则称其为迦辽金(Galerkin)方法[12]。
8