第6章弹性薄板小挠度弯曲问题的基础变分原理16K.docx

第6章弹性薄板小挠度弯曲问题的基础变分原理16K.docx

《第6章弹性薄板小挠度弯曲问题的基础变分原理16K.docx》由会员分享，可在线阅读，更多相关《第6章弹性薄板小挠度弯曲问题的基础变分原理16K.docx（27页珍藏版）》请在冰豆网上搜索。

第6章弹性薄板小挠度弯曲问题的基础变分原理16K

第6章弹性薄板小挠度弯曲问题的基础变分原理

平分板厚度的平面称为板的中面，一般地，当板的厚度不大于板中面最小尺寸的时的板称为薄板，薄板的中面是一个平面。

薄板在垂直于中面的载荷作用下发生弯曲时，中面变形所形成的曲面称为弹性曲面或挠度面，中面内各点在未变形中面垂直方向的位移称为板的挠度。

薄板弯曲的精确理论应是满足弹性力学的全部基本方程，但这在数学上将会遇到很大的困难。

1850年，G.R.基尔霍夫（KirchhoffGustavRobert，基尔霍夫古斯塔夫·罗伯特，德国物理学家，1824-1887年）除采用弹性力学的基本假设外，还提出了一些补充的假设，从而建立起了薄板小挠度弯曲的近似理论。

这些假设是：

第一，变形前垂直于板中面的直线，在板变形后仍为直线，并垂直于变形后的中面，而且不经受伸缩；第二，与中面平行的各面上的正应力与应力,和相比属于小量；第三，在横向载荷作用下板发生弯曲时，板的中面并不伸长，这也就是说，薄板中面内各点都没有平行于中面的位移分量。

用变分法可以导出薄板弯曲问题的平衡微分方程和边界条件。

当板的形状和边界条件较复杂时，直接求解偏微分方程时比较困难的，以变分法为基础的各种近似解是求解这类问题的一个重要途径。

本章讨论了用于薄板小挠度弯曲问题的一些基础变分原理，这包括虚功原理、最小位能原理、最小余能原理、两类自变量广义变分原理并推广到三类自变量广义变分原理。

§6.1基本方程与边界条件回顾

取坐标平面与中面重合，轴垂直于中面，,和轴构成一个右手直角笛卡儿坐标系。

变形后的板内各点沿,和轴方向的位移分别用,和表示。

由Kirchhoff假设，可以得到

，，（6-1）

并利用弹性力学中位移与应变之间的关系式，可以得到薄板中任意点的应变分量为

，，（6-2）

其余3个应变分量,和根据假设都等于零，即

，，（6-3）

由薄板的平衡关系，可以确定板的横向分布载荷与剪力,以及弯矩,和扭矩（,,统称为内力矩）与,之间的关系式。

这里要注意，,,是单位中面宽度内的内力矩，它们的因次是千克力，,是单位中面宽度内的内力，它们的因次是千克力/米。

弯矩、扭矩和剪力的正方向如图6-1所示。

图6-1弯矩、扭矩和剪力的正方向

平衡方程为

（6-4）

在薄板弯曲理论中，剪力,不产生应变，因而也不作功，因此可以从（6-4）式中消去,，得到

（6-5）

以后凡提到薄板弯曲平衡方程，都是指（6-5）式而言。

而内力,不再作为独立的量看待。

上面两组方程仅仅是力的平衡方程，它们未涉及到板的材料性质。

与内力矩相对应的广义应变是挠度面的曲率，在小挠度弯曲理论中，它们与挠度的关系为

，，（6-6）

内力矩与曲率的关系可以通过应变能密度表示出来，若将表示为的函数，则有

，，（6-7）

这种关系式对于线性或非线性材料都成立。

对于线性的弹性体，是的正定的二次齐次函数。

在各向同性的情况下，的算式为

（6-8）

将（6-8）式代入（6-7）式，然后再将（6-6）式代入，得到内力矩与挠度的关系式为

（6-9）

以上各式中称为板的弯曲刚度，其中为板的厚度，为材料的泊松系数。

如果我们定义为广义应变，为广义应力，即

（6-10）

则有

（6-11）

式中的为弯曲刚度矩阵。

（6-8）式可以写为

（6-12）

余应变能密度看作是内力矩,,的函数，其值定义为

（6-13）

并且有

，，（6-14）

同样，对于线性的弹性体，是,,的正定的二次齐次函数。

如果以广义应力表示余应变能密度，则有

（6-15）

式中。

（6-12）式与（6-15）式都是以后经常要用到的表达式。

注意，对于线弹性薄板，应变能密度与余应变能密度在数值上是相等的，即。

将（6-9）式代入（6-5）式，得到以挠度表示的各向同性薄板的平衡方程为

（6-16）

或

（6-16/）

在处理具体问题时，经常遇到坐标旋转而引起的变换。

如果坐标由转变为，如图6-2所示，则两个坐标系中坐标的关系为

（6-17）

对于挠度，有，从而

图6-2坐标转换

（6-18）

及二阶偏导为

（6-19）

弯矩、扭矩的变换公式为

（6-20）

剪力的变换公式为

（6-21）

在板的弯曲问题中，有三种典型的边界条件，简述如下。

设为板在平面上的定义域，板的边界为，令为沿边界外向法线的方向，为边界的切线，（,）的转向与（,）的转向是一致的，如图6-3所示。

第一种边界为固支边界，在这种边界上，其挠度与法向斜率均为给定的，即有

（在上）（6-22）

第二种边界为简支边界，在这种边界上，其挠度与法向弯矩为给定的，即有

（在上）（6-23）

第三种边界为自由边界，在自由边界上，作用在边界上的力为给定的。

从内力和力矩看，在边界上共有三个，即，但其中并不完全独立，因为从作功角度来看，和并不完全独立。

事实上，若边界上的挠度有一变分，则在上所作之功是

（6-24）

利用分部积分，上式又可以写成

（6-25）

由（6-25）式可见，切向扭矩可以分解为沿着周边边界的分布载荷及作用于两端的集中力，而两端是支座（不是固支边便是简支边）。

从实际板的受力来分析，可以看到集中力为作用在角点上，一般是影响到支座上的力，而对板的变形无影响。

因此，分布载荷与剪力构成沿自由边界上的分布力，这部分边界力的虚功为

与相对应的广义力为，自由边的边界条件应取为

（在上）（6-26）

为已知的作用在上的线分布载荷。

§6.2虚功原理和功的互等定理

力学上，可能位移是指满足位移连续条件的位移。

在薄板弯曲问题中只有一个广义位移，因此，可能作为可能位移的条件是：

是,的连续可导函数，

并且在边界上满足连续条件：

（6-27）

同样，由可能位移按式（6-10）也可得到相应的可能曲率。

可能内力是指与某种外力保持平衡关系的内力。

在薄板弯曲问题中，内力有,,，这三个内力组成一组可能内力的条件是：

在板的内部满足平衡方程（6-5）式，在板的边界上满足条件

（6-28）

根据能量守恒定理，外力在可能位移上所作的功等于可能内力在可能应变上所作的功，通常把这一关系叫做虚功原理。

在薄板弯曲问题中，若把支座反力也看作外力，则虚功原理的数学形式是

（6-29）

上式中，为可能挠度，是可能内力，它们之间可以完全独立而彼此无任何联系。

下面给出（6-29）式的数学证明。

为了书写简单，引入下面符号：

现在将取为的方向，取为的方向，则可以利用（6-18）、（6-20）、（6-21）式等将等用等表示出来，下面证明中将用到这些公式。

从（6-29）式中，等号右边两个线积分可作如下化简（并引用（6-22）、（6-23）式的边界条件），并得到

（6-30）

再将（6-4）式的关系代入（6-29）式右边第一个积分项里的中，展开后为

（6-31）

将（6-30）式和（6-31）式代入（6-29）式的右端，可以证明其左端等于右端。

对于虚功原理方程（6-29）式，还可以表示为以下恒等式

（6-32）

式中的代表（6-4）式前两个方程的缩写。

这里所谓恒等式，是指公式（6-32）中的是四个可以任意选取的函数。

该式要求具有一定连续可导性质，例如要求的一阶导数应该是连续而且是可导的。

利用上面说明的虚功方程（6-29）式，我们很容易导出功的互等定理。

在（6-29）式中，应再次指明内力与挠度是彼此独立的，它们之间是无任何联系的。

现在有两组载荷对同一块板作用，形成两组解，分别为

第一组载荷作用下，产生的内力与挠度为

第二组载荷作用下，产生的内力与挠度为

分别形成的虚功方程为

第一组外载及内力取第二组的位移（）为虚位移，有

（6-33）

第二组外载及内力取第一组的位移（）为虚位移，有

（6-34）

（6-33）式等号右边可以引用（6-11）式，得到下式

（6-35）

考虑到，，则（6-35）式可写成

（6-36）

（6-36）式即是薄板的功的互等定理，还可以写成

（6-37）

由于采用了线性的应力应变关系，所以无论是外力功的互等定理（6-37）式，还是内力功的互等定理（6-36）式，都是能量守恒原理和线性性质的后果。

§6.3最小位能原理

考虑板在横向分布载荷作用下处于平衡，并假定在板的边界上三种支持都存在的情况。

整个板的总位能包括两部分，一部分为板的应变能，它的算式为

（6-38）

为板的应变能密度，其算式如（6-12）式。

另一部分为外力包括分布载荷及边界力的位能，可写为

（6-39）

于是，整个板的总位能为

（6-40）

在最小位能原理中，挠度为唯一经受变分的自变函数，称这种变分为“一个自变函数的变分问题”。

令是精确解，与相应的弯矩、剪力为等，它们满足方程（6-4）式、（6-11）式和边界条件（6-22）、（6-23）及（6-26）式。

令为一个可能挠度，则最小位能原理指出：

与精确解相应的总位能小于任何其它可能挠度相应的总位能。

现在令

（6-41）

满足下面的边界条件

（在边界上）

（在边界上）（6-42）

与相应的总位能为

（6-43）

式中是把（6-41）式代入（6-38）式，以代替所得到的结果，即

（6-44）

其中

（6-45）

而

（6-46）

根据（6-29）式的虚功方程，可以证明

这样便有

（6-47）

从（6-44）式，不论为任何不全为零的组合，恒有。

因此有

（6-48）

这便是最小位能原理。

若将最小位能原理写成变分的形式，则有

（6-49）

利用分部积分，参考（6-30）的推导，由（6-40）式可以得到

（6-50）

§6.4最小余能原理

考虑与上节相同的薄板弯曲问题。

令为精确解。

再命为一组可能内力，它们满足下列方程

（6-51）

和在边界、上的边界条件

在上：

（6-52a）

在上：

，（6-52b）

系统的余能包括两部分。

一部分为余应变能，它的算式为

（6-53）

式中为余应变能密度。

对于线性的应力与应变关系，它可以表示为（6-15）式。

另一部分为已知的边界位移的余功，它的算式为

（6-54）

整个板的总余能为

（6-55）

总余能为自变函数的泛函。

现在取

，，，

，

且满足下列方程和边界条件

（6-56）

在上：

（6-57a）

在上：

（6-57b）

以上两式表示内力增量在边界上对应为零的外载荷。

并有

（6-58）

于是有

（6-59）

式中为内力增量相应的余应变能，而中间一项代表下列算式：

（6-60）

根据（6-29）式的虚功方程，可以证明

这样（6-59）式可以化为

（6-61）

如果不全为零，那么由（6-15）式可知

于是可得到

（6-62）

这便是最小余能原理。

将最小余能原理表达成变分形式，为

（6-63）

最小余能原理是一种条件变分原理，因为可能内力必须满足平衡条件（6-51）式。

Southwell（索斯韦尔）指出，利用应力函数方法可以把以上条件变分问题化为无条件变分问题。

齐次方程的解可以用两个应力函数与表示之，如

（6-64）

再命是平衡方程的一组特解。

于是可将内力表达为

（6-65）

将以上算式代入（6-63）式，可将余能表示成自变量和的泛函。