第五章弹性薄板小挠度弯曲问题的变分原理(16K)

89第五章弹性薄板小挠度弯曲问题的变分原理平分板厚度的平面称为板的中面，一般地，当板的厚度t不大于板中面最小尺寸的5/1时的板称为薄板，薄板的中面是一个平面。薄板在垂直于中面的载荷作用下发生弯曲时，中面变形所形成的曲面称为弹性曲面或挠度面，中面内各点在未变形中面垂直方向的位移称为板的挠度。薄板弯曲的精确理论应是满足弹性力学的全部基本方程，但这在数学上将会遇到很大的困难。1850年，G.R.Kirchhoff除采用弹性力学的基本假设外，还提出了一些补充的假设，从而建立起了薄板小挠度弯曲的近似理论。这些假设是：第一，变形前垂直于板中面的直线，在板变形后仍为直线，并垂直于变形后的中面，而且不经受伸缩；第二，与中面平行的各面上的正应力z与应力x,y和xy相比属于小量；第三，在横向载荷作用下板发生弯曲时，板的中面并不伸长，这也就是说，薄板中面内各点都没有平行于中面的位移分量。用变分法可以导出薄板弯曲问题的平衡微分方程和边界条件。当板的形状和边界条件较复杂时，直接求解偏微分方程时比较困难的，以变分法为基础的各种近似解是求解这类问题的一个重要途径。本章讨论了用于薄板小挠度弯曲问题的一些基础变分原理，这包括虚功原理、最小位能原理、最小余能原理、两类自变量广义变分原理并推广到三类自变量广义变分原理。§5.1基本方程与边界条件回顾取坐标平面oxy与中面重合，z轴垂直于中面，x,y和z轴构成一个右手直角笛卡儿坐标系。变形后的板内各点沿x,y和z轴方向的位移分别用u,v和w表示。由Kirchhoff假设，可以得到xwzzyxu),,(，ywzzyxv),,(，),(),,(yxwzyxw（5-1）并利用弹性力学中位移与应变之间的关系式，可以得到薄板中任意点的应变分量为22xwzx，22ywzy，yxwzxy22（5-2）其余3个应变分量z,xz和yz根据假设都等于零，即0z，0xz，0yz（5-3）由薄板的平衡关系，可以确定板的横向分布载荷),(yxq与剪力xQ,yQ以及弯矩xM,yM和扭矩xyM（xM,yM,xyM统称为内力矩）与xQ,yQ之间的关系式。这里要注意，xM,yM,xyM是单位中面宽度内的内力矩，它们的因次是千克力，xQ,yQ是单位中90面宽度内的内力，它们的因次是千克力/米。弯矩、扭矩和剪力的正方向如图5-1所示。平衡方程为),(yxqyQxQQyMxMQyMxMyxyyxyxxyx（5-4）在薄板弯曲理论中，剪力xQ,yQ不产生应变，因而也不作功，因此可以从（5-4）式中消去xQ,yQ，得到0),(22222yxqyMyxMxMyxyx（5-5）以后凡提到薄板弯曲平衡方程，都是指（5-5）式而言。而内力xQ,yQ不再作为独立的量看待。上面两组方程仅仅是力的平衡方程，它们未涉及到板的材料性质。与内力矩相对应的广义应变是挠度面的曲率xyyxkkk,,，在小挠度弯曲理论中，它们与挠度w的关系为22xwkx，22ywky，yxwkxy2（5-6）内力矩与曲率的关系可以通过应变能密度U~表示出来，若将U~表示为xyyxkkk,,的函数，则有xxkUM~，yykUM~，xyxykUM~21（5-7）这种关系式对于线性或非线性材料都成立。对于线性的弹性体，U~是xyyxkkk,,的正定的图5-1弯矩、扭矩和剪力的正方向91二次齐次函数。在各向同性的情况下，U~的算式为)])(1(2)[(21~22xyyxyxkkkkkDU（5-8）将（5-8）式代入（5-7）式，然后再将（5-6）式代入，得到内力矩与挠度的关系式为yxwDMxwywDMywxwDMxyyx222222222)1()()(（5-9）以上各式中)1(1223EtD称为板的弯曲刚度，其中t为板的厚度，为材料的泊松系数。如果我们定义}{为广义应变，M为广义应力，即xyyxxyyxMMMMyxwywxwkkk}{22}{22222（5-10）则有}]{[}{DM（5-11）式中的][D为弯曲刚度矩阵。（5-8）式可以写为}]{[}{21~DUT（5-12）余应变能密度*~U看作是内力矩xM,yM,xyM的函数，其值定义为UkMkMkMUxyxyyyxx~2~*（5-13）并且有xxMUk*~，yyMUk*~，xyxyMUk*~2（5-14）同样，对于线性的弹性体，*~U是xM,yM,xyM的正定的二次齐次函数。如果以广义应力}{M表示余应变能密度，则有}]{[}{21~*MCMUT（5-15）式中1][][DC。92（5-12）式与（5-15）式都是以后经常要用到的表达式。注意，对于线弹性薄板，应变能密度与余应变能密度在数值上是相等的，即*~~UU。将（5-9）式代入（5-5）式，得到以挠度表示的各向同性薄板的平衡方程为),()2(4422444yxqywyxwxwD（5-16）或),(22yxqwD（5-16/）在处理具体问题时，经常遇到坐标旋转而引起的变换。如果坐标由oxy转变为o，如图5-2所示，则两个坐标系中坐标的关系为cossin,sincoscossin,sincosyxyxyx（5-17）对于挠度w，有),(),(wyxw，从而cossinsincosywxwwywxww（5-18）及二阶偏导为sincos)sin(cossincoscoscossin2sinsincossin2cos22222222222222222222222222ywyxwxwwywyxwxwwywyxwxww（5-19）弯矩、扭矩的变换公式为sincos)sin(cossincoscoscossin2sinsincossin2cos222222yxyxyxyxyxyxMMMMMMMMMMMM（5-20）剪力的变换公式为图5-2坐标转换93cossinsincosyxyxQQQQQQ（5-21）在板的弯曲问题中，有三种典型的边界条件，简述如下。设为板在xy平面上的定义域，板的边界为C，令n为沿边界外向法线的方向，s为边界的切线，（n,s）的转向与（x,y）的转向是一致的，如图5-3所示。第一种边界为固支边界1C，在这种边界上，其挠度与法向斜率均为给定的，即有nn,（在1C上）（5-22）第二种边界为简支边界2C，在这种边界上，其挠度与法向弯矩为给定的，即有nnMMww,（在2C上）（5-23）第三种边界为自由边界3C，在自由边界上，作用在边界上的力为给定的。从内力和力矩看，在边界上共有三个，即nnsnQMM,,，但其中并不完全独立，因为从作功角度来看，nsM和nQ并不完全独立。事实上，若边界上的挠度有一变分wδ，则nnsQM,在wδ上所作之功wδ是swQswMwCnnsd]δδ[δ3（5-24）利用分部积分，上式又可以写成33|δdδ)(δCnsCnnswMswQsMw（5-25）由（5-25）式可见，切向扭矩nsM可以分解为沿着周边边界3C的分布载荷sMns及作用于3C两端的集中力||nsM，而3C两端是支座（不是固支边便是简支边）。从实际板的受力来分析，可以看到集中力||nsM为作用在角点上，一般是影响到支座上的力，而对板的变形无影响。因此，分布载荷sMns与剪力nQ构成沿自由边界3C上的分布力，这部分边界力的虚功为swQsMCnnsdδ)(3与wδ相对应的广义力为nnsQsM，自由边的边界条件应取为)(,sqQsMMMnnsnn（在3C上）(5-26))(sq为已知的作用在3C上的线分布载荷。图5-3板的边界94§5.2虚功原理和功的互等定理力学上，可能位移是指满足位移连续条件的位移。在薄板弯曲问题中只有一个广义位移),(yxw，因此，),(yxw可能作为可能位移的条件是：ywxww,,是x,y的连续可导函数，并且在边界上满足连续条件：上）（在上）（在21,CwwCn（5-27）同样，由可能位移w也可得到相应的可能曲率。可能内力是指与某种外力保持平衡关系的内力。在薄板弯曲问题中，内力有xM,yM,xyM，这三个内力组成一组可能内力的条件是：在板的内部满足平衡方程（5-5）式，在板的边界上满足条件上）（在上）（在32)(,CsqQsMMMCMMnnsnnnn（5-28）根据能量守恒定理，外力在可能位移上所作的功等于可能内力在可能应变上所作的功，通常把这一关系叫做虚功原理。在薄板弯曲问题中，若把支座反力也看作外力，则虚功原理的数学形式是CnCnnsyxyxsnwMswQsMyxqwyxywMyxwMxwMdd)(dddd)2(22222（5-29）上式中，w为可能挠度，xyyxMMM,,是可能内力，它们之间可以完全独立而彼此无任何联系。下面给出（5-29）式的数学证明。为了书写简单，引入下面符号：),cos(),,cos(ynmxnl现在将取为n的方向，取为s的方向，则可以利用（5-18）、（5-20）、（5-21）式等将swnwQMMnnsn,,,,等用ywxwQQMMMyxxyyx,,,,,,等表示出来，下面证明中将用到这些公式。从（5-29）式中，等号右边两个线积分可作如下化简（并引用（5-22）、（5-23）式的边界条件），并得到CCnnnssnwMswQsMdd)(95sywmMywlxwmMxwlMswmQlQsywlxwmMmlMMlmywmxwlMmlmMMlswmQlQsswMnwMswQyCxyxCyxxyxyCyxyxCyxCCnsnnd})({d)(d)}(])()([)()2{(d)(d)(d2222（5-30）再将（5-4）式的关系代入（5-29）式右边第一个积分项里的q中，展开后为yxywMyxwMxwMsywmMywlxwmMxwlMswmQlQyxywyMxMxwyMxMswmQlQyxywQxwQswmQlQyxwyQxQyxqwyxxyxyCxCyxyxyxyxCyxyxCyxyxdd)2(d])([d)(d]d)()[(d)(dd)(d)(dd)(dd22222（5-31）将（5-30）式和（5-31）式代入（5-29）式的右端，可以证明其左端等于右端。对于虚功原理方程（5-29）式，还可以表示为以下恒等式CnnnsyxyxyxsnwMwQsMyxywMyxwMxwMwyQxQd}){(dd}2){(22222（5-32）式中的yxQQ,代表（5-4）式前两个方程的缩写。这里所谓恒等式，是指公式（5-32）中96的wMMMxyyx,,,是四个可以任意选取的函数。该式要求wMMMxyy