第六章薄板的失稳大型梁、柱等结构构件,通常都由板件组合而成,为了节省材料,板件通常宽而薄,薄板在面内压力作用下就可能失稳,并由此导致整个构件的承载力下降因此,对板件失稳和失稳后性态的研究也是结构稳定的重要问题。6.1板的稳定微分方程板按照其厚度可分为厚板、薄板和薄膜三种。将板的厚度t与板幅面的最小宽度b相比:如果t/b>1/5~1/8时,板被称为厚板;当1/80~1/100<t/b<1/5~1/8,这个范围内的板称为薄板;而当t/b<1/80~1/100时,板被称为薄膜。薄膜没有抗弯刚度,完全靠薄膜张力来支承横向荷载作用。薄板不仅具有抗弯刚度,还可能存在薄膜张力。薄板在横向荷载作用下,如果产生的挠度w≤t/5,则属于薄板的小挠度弯曲问题。薄板中的物理方程和内力表达式与弹性力学中平面应力问题的物理方程和内力表达式相同。薄板中于薄板上下表面等距离的面称为中面。当薄板在中面内承受平行于中面的荷载而失稳时,也可以根据静力平衡准则来确定板的临界荷载。下面根据小挠度理论给出薄板的弹性稳定微分方程。根据小挠度理论给出薄板的弹性稳定微分方程2222244224442)2(ywNyxwNxwNywyxwxwDyxyx式中:,为单位宽度板的抗弯刚度。上式是一个以挠度w为未知量的常系数线性四阶偏微分方程。板的边界条件表达式:1)简支边:w=0,2)固定边:w=0,3)自由边:022yw022xw00ywxw,022xw022yw0)2(2333yxwxw,,,)1(1223EtD6.2受压简支板的弹性失稳如图所示四边简支矩形板,板的中面上作用有荷载Nx=-Px、Nx=0、Nxy=0;因此,板的弹性稳定微分方程式为板的弹性稳定微分方程式为:根据板的边界条件:当x=0和x=a时,w=0、、;当y=0和y=b时,w=0、、。符合这些边界条件的板的挠曲面可用二重三角级数表示为式中:m和n分别是板失稳时,在x和y方向的半波数,为各项的待定常数。0)2(224422444xwPywyxwxwDx022xw022yw022xw022yw11sinsinmnmnbynaxmAwmnAnm;,3,2,1,,3,2,1对w微分两次和四次后代入偏微分方程,得由于和均不为零,Amn也不为零,否则板仍然为平面平衡状态,所以解得式中:。只有当n=1时,上式有最小值,所以有意义的临界荷载为:11222444224224440sinsin2mnxmnbynaxmamDPbnbanmamAwaxmsinbymsin0222244422422444amDPbnbanmamx2222mbanambbDPx)1(1223EtD222mbaambbDPcrx或;式中:K—为稳定系数,由,可得,K=4,则。要使上式成立,m=a/b必须是整数,如果m=a/b不是整数,则或计算临界荷载时,m的取值应使K值最小。当n=1时,即在y方向成一个半波的条件下,K随m和a/b成曲线关系,如前面图所示。22bDKPcrx2)(mbaambK0dmdKbam224bDPcrx222)(mbaambbDPcrx22bDKPcrx板的长宽比在时,m=1,板以一个半波的形式失稳;长宽比在与之间时,m=2,板以两个半波的形式失稳;长宽比在与之间时,m=3,余类推。而当长宽比时,K值已非常接近于最小值。板的临界应力:由上式可知:单向均匀受压板的临界应力与板的宽厚比的平方成反比,而与板的长度无关。对于单向均匀受压的矩形板,当加载边为简支,而非加载边为各种不同的支承条件时,稳定系数K的最小值如下表所示。2/ba266120.4/ba222)/()1(12tbEKtPcrxcrx稳定系数K的最小值序号12345非加载边的支承条件一边简支一边自由一边固定一边自由两边简支一边简支一边固定两边固定稳定系数K0.4251.2804.0005.4206.970只有当时,同时边界条件由简支变为固定,稳定系数K才会有较大提高。通过上述讨论可知,对于单向均匀受压的狭长板,用增加横向加劲肋来改变a/b,从而提高稳定系数的做法并无明显的效果;如果把加劲肋的间距取得小于2b又很不经济。而对于很宽的薄板,如果采用纵向加劲肋以减少板的宽度b倒是有效的。例如在板的纵向中心加一条加劲肋时。2/ba2216bDPcrx屈曲系数与板件长宽比的关系失稳系数与板件长宽比的关系6.3能量法计算板的弹性失稳对于四边简支的均匀受压板,计算公式中每一项都有,因此,可以从各项中将其分离出来,这样用平衡法求解很方便;而当板的支承条件不是简支时,三角函数则无法分离,这时就需要用能量法来求解。1)板的总势能已知Π=U+V,则弹性应变能外力势能bynaxmsinsindxdyyxwywxwywxwDUab00222222222212abxyyxdxdyywxwPywPxwPV00222212)瑞利—李兹法假设符合板的几何边界条件的挠曲面函数为将此式代入总势能公式中,经积分后,根据势能驻值原理建立一组的线性代数方程组。其非零解的条件是方程组的系数行列式为零,即可得板的稳定方程。3)伽辽金法已知板的平衡偏微分方程为,需假定符合板的几何与自然边界条件的挠曲面函数,现假定挠曲面函数11),(mnmnyxfAw0,,0,01211mnAAA0)(wL为可组成伽辽金方程组上面方程组经积分后可得到A1,A2,……An的线性方程组,为得到它们的非零解,其系数行列式︱D︱=0应为零,则可得稳定方程,由稳定方程可解的临界荷载。niiiyxAw1),(abnababdxdyyxwLdxdyyxwLdxdyyxwL000020010),()(0),()(0),()(6.4均匀受剪简支板的弹性失稳均匀受剪的四边简支板如图所示,在其对角线方向因受压而失稳,失稳时板的波长与另一对角线方向的拉力有关。对于长板,失稳时的半波长度约为板宽的1.25倍。当采用能量法求解剪切临界荷载时,板的挠曲面函数可用二重三角函数表示。但是对于均匀受剪的四边简支板,可以利用均匀受剪四边简支的正方形板来求解临界荷载的近似值。如采用伽辽金法求解时,板的中面力而。板的平衡偏微分方程为设满足几何和自然边界条件的挠曲面函数为可得伽辽金方程组为将挠曲面函数w(x,y)的微分代入平衡偏微分方程式中,经积分后得到yxxyyxxyppNN0,0yxNN022)(24422444yxwDpywyxwxwwLxyayaxAayaxAw2sin2sinsinsin210sinsin)(00dxdyayaxwLaa02sin2sin)(00dxdyayaxwLaa板的失稳条件是解得这个均匀受剪简支正方形板的剪切失稳临界荷载计算结果与精确解相比,误差为19%;如果采用更多项的挠曲面函数,则可提高解的精确度。09322124ADpAaxy0169322241AaADpxy0169329322424aDpDpaxyxy22241.1189aDaDpxy经过对矩形板更精确的理论分析,可得式中:ks为剪切失稳系数。均匀受剪矩形板的剪切稳定系数如下:对于四边简支的受剪板当a≥b时,当a≤b时,对于四边固定的受剪板当a≥b时,当a≤b时,22bDkpsxy2)/(0.434.5abks2)/(34.50.4abks2)/(6.598.8abks2)/(98.86.5abks6.5单向受压简支板的失稳后强度以上研究的板的失稳都是建立在小挠度理论基础上的,即认为板在失稳时的挠度远小于其厚度,忽略了板失稳时中面上的薄膜拉力。如果板的支承构件刚度较大,板的临界应力虽然不高,但失稳后并不破坏。板中的应力将重新分布,并产生薄膜拉力,使板的承载能力远远超过其临界荷载,这种现象被称为失稳后强度。此时,板的挠度与板的厚度相比已不是一个小量,所以需按大挠度理论来求解。6.5.1平衡偏微分方程薄板失稳后,板的中面产生了数值远大于其厚度的挠度,外荷载作用下中面力已不再是常量,在非荷载作用的方向也同时产生了中面力,因此需按大挠度理论研究薄板的失稳后强度。如图所示,从板中取出的微元体dxdyt上作用着中面力Nx、Ny和Nxy,由这些中面力在x方向的平衡条件,忽略其中的高阶微量后可得同理,由这些中面力在y方向的平衡条件可得由这些中面力在z方向的平衡条件可得由上式可见大挠度理论的平衡方程与小挠度理论平衡方程式的形式相同,但它是变系数的。已知中面力Nx、Ny和Nxy都包括了作用于中面的外荷载和因为板挠曲而产生的薄膜力,所以都是变量。这三个方程中有了四个因变量,属于变系数偏微分方程,并需要根据板的变形条件补充一个变形协调方程。0xNyNxyy0yNxNxyx2222244224442)2(ywNyxwNxwNywyxwxwDyxyx6.5.2变形协调方程薄板微元体中面的应变与挠度的变形协调方程为6.5.3薄板的大挠度方程组为了简化计算过程,应设法减少未知量,可引入满足中面力平衡方程的应力函数F(x,y)。则将上面表达式代入平衡方程和变形协调方程中,则可得到以挠度w和应力函数F(x,y)为变量的力平衡方程和变形协调方程2222220202202ywxwyxwyxxyxyyx22yFtNx22yFtNyyxFtNxy2-薄板的大挠度方程组为6.5.4单向受压简支板的失稳后强度单向均匀受压四边简支矩形板如下图所示,现研究其失稳后的强度。yxwyxFywxFxwyFDtywyxwxw222222222244224442)2(22222244224442ywxwyxwEyFyxFxF矩形板在平面外的边界条件为:当x=0和x=a时,w=0和;当y=0和y=b时,w=0和。板在发生失稳后会产生应力重新分布,故在研究板的失稳后强度时,还要考虑在板自身平面内的边界条件。为此假设:板失稳后其外形不变,及板的边缘仍保持直线;沿板的四周不产生剪应力,即Nxy=0;y=0和y=b的两条边在y方向可以自由移动。在上述假设条件下来求解板的失稳后强度,假设符合板边界条件的挠曲面函数为022xw022ywbyaxmfwsinsin将挠曲面函数代入变形协调方程式,即上式的全解由特解为Fp和通解为Fc两部分组成,其特解为将上式代入变形协调方程式,可得,。由此可得通解Fc由平衡偏微分方程式的齐次方程求得,即相当于w=0,处在板失稳前的平衡状态。byaxmbamEyFyxFxF2cos2cos21222244422444byCaxmBFp2cos2cos222232fbmEaB222232fabEmCbyabmaxmbmafEFp2cos2cos322222222此时板的中面力Nx=px,Ny=0和Nxy=0;可得积分得。则变形协调方程式的全解为应力函数F(x,y)与