弹性力学总结

弹性力学关于应力变分法问题

一、起源及发展

1687年，Newton在《自然哲学的数学原理》中提出第一个变分问题一一定轴转动阻力最小的旋转曲面形状问题；1696年，Bernoulli提出了著名的最速降线问题；到18世纪，经过Euler, Lagrange等人的努力，逐渐形成变分法。古典变分法的基本容是确定泛函的极值和极值点，它为许多数学、物理、科技、工程问题提供了强有力地数学工具。现代理论证明，微分方程（组）中的变分法是把微分方程（组）化归为其对应泛函的临界点（即化为变分问题），以证明其解的存在性及解的个数。讨论对应泛函临界点的存在性及其个数的基本方法是Morse理论与极小极论（Minimax Theory） <>变分法有着深刻的物理背景，某种意义上，自然界一切物质运动均可以用某种形式的数理方程表示，一般数理方程又与一定的泛函相对应，所以一切物质运动规律都遵从“变分原理”。

由于弹性力学变分解法，实质上就是数学中的变分法应用于解弹性力学问题，虽然在讨论的近似解法中使用变分计算均甚简单（类似微分），但"变分” 的概念却极为重要，它关系到我们队一系列力学变分原理中“虚”的概念的建立与理解。以下，就应力变分法进行讨论。

二、定义及应用

（1）、应力变分方程

设有任一弹性体，在外力的作用下处于平衡。命円为实际存在的应变分量, 它们满足平衡微分方程和应力边界条件，也满足相容方程，其相应的位移还满足位移边界条件。现在，假想体力和应变边界条件上给定的面力不变而应力分量发生了微小的改变茨％,即所谓虚应力或应力的变分，使应力分量成为5j+叫假定他们只满足平衡微分方程和应力边界条件。

既然两组应力分量都满足同样体力和面力作用下的平衡微分方程和应力边界条件，应力分量的变化必然满足无体力时的平衡微分方程。即

(a)

在位移给定的边界上，应力分量的变分必然伴随着面力分量的变分

根据应力边界条件的要求，应力分量的变分在边界上必须满足

㈣+〃％+,% =3几 mda s +ndT xz +/Jr n . = >

zzJcr. + Idt^ + mdr^ = 3f z a 则应变余能的变分应为 6V C = J J J 5 v c dxdydz, = Jf J (色-^cr v + …+2- + • • Jdxdydz,。

将上式代入，得

SV C = JUG® + …+ 7V Q 兀 + ・・ Jcbafyd 乙。

再将几何方程代入，得

SV c — Ilf[ — &7丫 + ・• • + (―— + —)^r v . + - • -]dxdydzo ex cy cz

9r

根据分部积分和奧一高公式，对上式右边进行处理:

.1 ? d(J x dxdydz = U lud

ex

ex

最后可得 5 匕=f J[“(/5v + 恥 j + 加r J + …]dS -

\ \ 込 + 二

+ — 5 g) + …]〃xdydz 。

ex oy cz 再将(a)、(b)代入，即得 6-av

(b)

^V=U {u8f x + v3f y + w3f z)dS°

这就是所谓应力变分方程，有的文献把它叫做卡斯蒂利亚诺变分方程。最小余能原理：

北("兀+诂兀+加兀)dS=O。

上式也可以改写为：

世 _ J i(uf x + v7y + wf z)dS ] = Oo

(2)、应力变分法

由推到出的应力变分方程，使其满足平衡方程和应力边界条件，但其中包含若干待定系数，然后根据应力变分方程解决这些系数，应力分量一般可设为：b%Jo+工九(bj (C)

m tn

其中观是互不依赖的m个系数，(％)“是满足平衡微分方程和应力边界条件的设定函数，(qj,”是满足“没有体力和面力作用时的平衡微分方程和应力边界条件”的设定函数。这样，不论系数A m如何取值，(爲丄总能满足平衡微分方程和应力边界条件。

注意：应力的变分只是由系数Am的变分来实现。

如果在弹性体的每一部分边界上，不是面力被给定，便是位移等于零，则应力变分方程得也=0 ,

即：仝丄=o (d)

込

应变余能匕是％的二次函数，因而方程(d)将是Am的一次方程。这样的方程共有m个，恰好可以用来求解系数,Am从而由表达式(c)求得应力分量。

如果在某一部分边界上，位移是给定的，但并不等于零，则在这一部分边界上须直接应用变分方程（11-18）,即

我.=J f x +v5/s 4-wdf.）dSo

在这里，u, V、\v是已知的，积分只包括该部分边界，面力的变分与应力的变

分两者之间的关系即：

5/、= u^a- + n 3r y , +"春 >

J/. = 11&6 + l&J + mdTo

带入方程的右边积分后，将得岀如下的结果：

+ 诂兀 + W 万 JdS = £ B 加

m

其中Bm 是常数，另一方面，我们有:

因而得:

这将仍然是九的一次方程而且总共有m 个,仍然可以用来求解系数 从而由表达式（C ）求得应力。

（3）、应力函数方法

由于应力分量的数量有点多，确定起来较为困难，通常用应力函数方法。 在平面应力问题中，如果体力分量为常数，则存在应力函数。将应力函数设 为：

①二①0 +工九①加’

m

其中每为互不依赖的m 个系数。这样就只需使①（）给出的应力分量满足实 际的应力边界条件，并使①川给岀的应力分量满足无面力时的应力边界条件。

在平面应力问题中，有=T v ： =T ：x =0 ,而且b*、b v 、和不随坐标Z 而

变。在Z 方向取一个单位厚度，则用应力分量表示的应变余能表达式为 K =袒 J +

cr v 2- + 2（1 + //）r n .2Wyo

2匕

（加

=12…）