弹性力学经典变分原理资料下载.pdf

1、图 3.1 物质坐标系和空间坐标系矢量PQ表示了质点P的位移，记为u。为简单和方便起见，一般取两个参照系相重合，这时位移矢量u的分量iu可以用下式来表示,（1,2,3）iiiuxXi（3.1.1）其中变形后质点的坐标）3,2,1（ixi与变形前的坐标）3,2,1（iXi存在着确定的关系。我们可以把变形后质点的坐标看成是变形前质点物质坐标的函数，即123（,）,（1,2,3）iixx XXXi（3.1.2）也可以用其逆变换（数学上要求Jacobi 行列式不为零）来表述，也就是从变形后空间坐标描述的质点，来追涉变形前这一质点的坐标123（,）,（1,2,3）iiXX x x xi（3.1.3）如果

2、把位移u看作是变形前坐标、即物质坐标的函数123（,）,（1,2,3）iiuu XXXi（3.1.4）称之为Lagrange 描述描述。如果把位移u看作是变形后坐标、即空间坐标的函数123（,）,（1,2,3）iiuu x x xi（3.1.5）称之为Euler 描述描述。我们取变形前P点）,（321XXX及相邻112233（d,d,d）P XXXXXX，它们之间的长度平方为3201dddiiisXX（3.1.6）它们变形后相应的Q点）,（321xxx及相邻112233（d,d,d）Q xx xx xx，其长度平方为321dd diiisx x（3.1.7）根据变形前后的坐标关系有3311dd

3、,ddiiijjjjjjxXxXXxiXx从而有33220,11dd（）ddijiji jijxxssXXXX（3.1.8）或者33220,11dd（）d dijiji jijXXssx xxx（3.1.9）如果定义3121ijijijxxEXX（3.1.10）及3121ijijijXXxx（3.1.11）则有220dd2ddijijssEXX（3.1.12）220dd2d dijijssx x（3.1.13）上述表达式中，有重复下标的,i j，已省略了相应的求和记号3311,ij，称为 Einstein 约定。我们称ijEE为 Lagrange-Green 应变张量（用 Lagrange 坐

4、标系来描述），把ij 称作为Euler-Almansi 应变张量（用 Euler 坐标系来描述）。如果我们在 Lagrange 坐标系中,沿着某一个特定的坐标方向取一个微分单元1123d（d0,dd0）XXXR,其变形前长度为01ddsX而变形后的长度为110d12dsEs因此，该微段变形前后的相对伸长量为01110dd121dssEEs（3.1.14）可见11E与线元的相对伸长有关。当111E时，111EE。如 果 在 Lagrange 坐 标 系 中 沿 坐 标 轴 方 向 取 两 个 相 互 垂 直 的 微 元，分 别 为11d（d,0,0）XR和22d（0,d,0）XR，它们的长度分别

5、为011022dd,ddsXsX那么在变形后它们长度1ds和2ds分别为1111d12dsEX（3.1.15）2222d12dsEX（3.1.16）变形后两个微段对应向量的内积为1212121212cos d ddd2ddkkxxs sXXEXXXX112212kkxxEXX（3.1.17）其中为变形后两个微段之间的夹角。所以1212121211222dd2cosd d1212EXXEs sEE（3.1.18）如果记变形前后两个微元之间夹角的变化（减少）为，也就是说2（3.1.19）那么1211222sincos1212EEE（3.1.20）当11221,1EE时，可以表示为12122,/2E

6、E（3.1.21）所以说，12E是与剪切变形有关的量。如果用空间坐标系来描述变形，也就是说，位移矢量u的分量iu用变形后的坐标来描述1,23（,）,（1,2,3）iiiiiuxXxX x x xi1,23（,）,（1,2,3）iiiiiXxuxu x x xi那么312131213121ijijijijijijjijiijXXxxuuxxuuuuxxxx（3.1.22）在小变形情况下，如果忽略高阶小量后，那么有12jiijjiuuuu（3.1.23）我们称之为 Cauchy 微小应变。在工程上描述的应变为xux，yvy，zwzywzvyz，xvyuxy，zxuwzx把他们写成矩阵的形式为000

7、000000Txyzyzzxxyxzyuvyzxwzyx （3.1.24）也就是（）TEu（3.1.25）其中TuvwuTxyzyzzxxy 000（）（,）000000 xzyxyzyzxzyxE=E 式中 代表梯度算子xyzijk i,j,k代表zyx,方向的单位向量。3.1.2 应力分析应力分析图 3.2 物体受力如图所示,通常作用于物体的外力可以分为两种：一种是分布在物体表面的作用力，例如一个物体对另一物体作用的压力，象水压力等，我们称之为面力（surface traction）；另一种是分布在物体体积内部的力，象重力、磁力或运动物体的惯性力等，我们称之为体力（bodyforce）。图

8、 3.3 内力和应力当一个物体处于平衡状态时，假如我们设想从中分离出一部分B，其表面用S表示。S上任意一点，其邻域S面上作用的合力为F，应力0limSS Fp正应力剪应力截面上应力,（）p与截面法向有关。当取定坐标系统xoy后,可以用每个坐标面上的沿坐标轴的三个应力分量来表示应力状态。根据剪应力互等定律,其中独立的分量有 6 个,我们记为应力张量（满足坐标变换规律）xxyxzijyxyyzzxzyz,zyyzzxxzyxxy,（3.1.26）应力的符号规则:外法线方向与坐标轴方向一致的截面上,沿坐标轴正方向的应力为正,沿坐标轴负方向的应力为负；反之,外法线方向与坐标轴方向相反的截面上,沿坐标轴

9、正方向的应力为负,沿坐标轴负方向的应力为正。3.1.3 截面上应力截面上应力在某一个外方向（,）Txyzn n nn的截面上，根据力的平衡关系，截面上应力p沿三个坐标轴上的应力分量为31,1,2,3iijjjpni也就是说xxxxyyxzzpnnnyyxxyyyzzpnnnzzxxzyyzzpnnn图 3.4 应力张量与截面上应力写成矩阵形式为000000000 xyxxzyzyyzxyzzzyxzxxypnnnpnnnpnnn（3.1.27）也就是说（）pE n（3.1.28）式中TxyzppppTxyzyzzxxy（）E n就是将（）E 中的梯度矢量替换成截面的法向单位矢量n，即000（）

10、000000 xzyyzxzyxnnnnnnnnnE n（3.1.29）3.1.4 平衡方程平衡方程应力分量在物体内部的平衡方程为3,10,1,2,3ij jijfi（3.1.30）写成分量的形式为0 xzxxyxfzyx0yyzyxyfzyx0zzyzzxfzyx其中zyxfff,分别是体积力在zyx,轴上的分量。如果把平衡方程表示成矩阵的形式为0000000000 xyxzyyzzzxxyxzyffyzxfzyx也就是（）0Ef （3.1.31）式中Txyzffff3.1.5 应变能、余应变能及应力与应变关系应变能、余应变能及应力与应变关系物体发生弹性变形时，外力所做的功等于物体中所储存的

11、应变能。而这种应变能与物体的变形过程无关，只同物体的最终变形状态有关，也就是说只与最终的应变有关。我们在物体中隔离出一个微元d d dx y z。该微元上的应变分量为xyzxyzzyx,，作用微元表面上的应力分量为xyzxyzzyx,。记物体的应变能密度为U（也就是单位体积的应变能），那么储存在该微元上的应变能为d d dU x y z。根据前面的说明，应变能密度U应该是应变分量xyzxyzzyx,的函数。如果此时微元的应变有一个微小变化xyzxyzzyx,，相应的应变能密度也有了一个微小的变化U，根据能量守衡，有xyxyzxzxyzyzzzyyxxU从中我们可以得到xxU，yyU，zzUyz

12、yzU，zxzxU，xyxyU写成矩阵的形式为TU（3.1.32）TU 用积分形式表示为dddddddTxxyyzzyzyzzxzxxyxyU 通过下式定义的是余应变能密度VUVxyxyzxzxyzyzzzyyxx（3.1.33）也就是说xyxyzxzxyzyzzzyyxxVU（3.1.34）写成矩阵的形式有TVU （3.1.35）用积分形式表示为dddddddTxxyyzzyzyzzxzxxyxyV利用应力与应变之间的关系，可以把上式右边表示成应力的形式，也就是说把V表示成应力分量xyzxyzzyx,的函数。对上式取变分，有TTTTVU （3.1.36）因此有TV （3.1.37）图 3.5

13、 应变能和余应变能密度上面这些关系对于线弹性变形和非线性弹性变形都是适用的。对于非线性的弹性变形，U和V不仅在数学形式上不一样，而且在数值也不相等。对于线弹性变形，应力和应变之间关系是线性的，应变能密度U和应变余能密度V数值上相等12TUV （3.1.38）如果应力和应变之间的关系表示为 A a那么U和V可以表示成12TU A（3.1.39）12TV a（3.1.40）由于能量的正定性,A和a都必须是对称正定的六阶矩阵,而且它们之间互为逆矩阵，也就是说AaI这里I是六阶单位矩阵。3.1.6 边界条件边界条件在弹性力学的定解问题中，除了必要的微分方程外，还需要给定合适的条件。这种边界条件是多种多

14、样的，我们这里只讨论两种典型的情况，即给定位移的位移边界和给定面力的应力边界。记B为物体的总边界，我们可以把总边界分成两部分12B=B+B。其中1B上有位移边界条件u=u（3.1.41）2B上有应力边界条件（）E n=p（3.1.42）其中TuvwuTxyzpppp分别是边界1B上给定的位移向量和2B上给定的单位面积上外力向量。图 3.6 边界条件3.1.7 几何可能位移和静力可能应力几何可能位移和静力可能应力几何可能位移几何可能位移:在位移边界1B上满足位移边界条件u=u，且在整个区域内满足连续条件（可以得到相应的应变）的位移称为几何可能位移几何可能位移，一般用ku来表示。静力可能应力静力可

15、能应力:在应力边界2B上满足应力边界条件（）E n=p，且在整个区域内满足应力平衡条件（）0E+f=的一组应力称为静力可能应力静力可能应力，一般用s 来表示。3.1.8 弹性力学精确解弹性力学精确解弹性力学的精确解u,应满足下列微分方程和边界条件（1）几何关系（）T Eu ，内（2）平衡方程（）0Ef ，内（3）本构关系TU 或者TV ，内（4）边界条件（）E n=p，2B上u=u，1B上3.2 一个重要的恒等式一个重要的恒等式对于三维空间上任意一个连通区域，始终成立下面的恒等关系恒等关系（）d（）d（）dTTTTBB EuE nuEu（3.2.1）其中B是区域的边界，（,）Txyzn n nn是边界B上的外法线方向，3uR和6R 是两组任意独立的函数，式中000（）000000 xzyyzxzyxE 000（）000000 xzyyzxzyxnnnnnnnnnE n写成分量的形式为333,1,1,1dddiji jijjiij jii ji ji jBun u Bu 证明证明:123（）xyzEEEE 123（）xyznnnE nEEE其中123,E EE是三个常数矩阵。12310000000000100001