弹性力学关于应力变分法问题一、起源及发展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 cz9r根据分部积分和奧一高公式，对上式右边进行处理:.1 ? d(J x dxdydz = U lud<y x dS - g (Sbjdxdydz,exex最后可得 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 .2Wyo2匕（加=12…）对于平面应变问题，［（1 一”）（b, + b；） - 2〃bs + 2 叮｝dxdy.如果所考虑的弹性体是单连体，体力为常量,应力分量£、5.应当与<1 y -V“无关，可以取“二0,于是平面应力情况下的表达式和平面应力情况下的表达式都简化为K=^-f J（b「+ 町 + 2 叮）dxdy.即得用应力函数表示应变余能的表达式K需川（21-加+（竽-小+2（竽）2畑。
2E 0- ax cxoy在应力边界问题中，因为面力不能有变分,旳=0。
应为应力分量以及应变余能的变分是通过系数Am的变分来实现的，所以上式归结为将将应力函数表达式代入，即得川（云丁 _/>）=（尹）+（=_/』）丁（尹）+ 內。
九dy dx % dx 6?①d6’①2 ----------- （——）］dxdy = ^dxdy dA m dxdy= 12 …）可以用来决定系数A叫从而确定应力函数°,再由应力函数°求得应力分量。
由于是近似解，应力分量不能精确满足相容条件，由应力分量求得的应变分量也不能精确满足变形协调条件，不能根据几何方程求得位移分量。
应力函数法的要点是要找到满足全部边界条件的应力函数，二这种函数一般任然难以找到，尤其在边界不规整的情况下。
所以应力方法的应用在这一点上受到极大的限制。
(4)、典型例题：例1：设有宽度为2a,高度为b的矩形薄板，左右两边和下边被固定约束，2上边的位移被给定为"=o v = -/Z(l-4).不计体力。
试求薄版的位移分量和a~应力分量。
解：取坐标系底部为x轴，对称轴为y轴，则该问题是一个轴对称问题一一及约束情况，几何形状以及所受的外来因素都对称于某个坐标轴。
本题中，对称轴显然是y轴。
这样，位移u,v关于y轴对称。
首先考察位移u：薄板左右两边：= 0 (说明u中含有(x2-a2)项或(a2 -x2)项)薄板下边：(")冋=0 (说明u中含有(y-O)项)薄板上边：(“)心=0 (说明u中含有(y-b)项或(b-y)项)所以u所以表达成:u =Ai(a2-x2)y(b-y)(这里m=l,即取一个系数州)由此可得u,v的表达式为：w = A l(l-—)—(1--)cf a a a►V =_7(I_4)2+B1(I-4)2(I-Z)cr b CT b bW = °可以满足位移边界条件：(叽o =0(V).v=O = 09(5=0cT由于U是X的奇函数，V是X的偶函数, 对称条件满足。
此外，由(i)得："]=( 疋y y2 万一护万一訐vl=(l_M-p) cr b少du at/ r- f 由瓦引心血丽訂几吨du dA{= -c h ab^ = -q2abEab2(1-v2)(2Aj + 2vB[) = -q l abEab2(1 — /)(2Bj + 2vA,) = -q2abyE E直杆弯曲的应变能为，在不考虑剪切效应时下面用最小势能原理来确定参数,(2a 2 +6«3) dx-F(a 2l3 +a y C) 由最小势能原理三、总结与思考所谓弹性力学的变分解法就是基于力学能量原理求解弹性力学的变分方法， 这种方法从其本质而言，是要把原来在给定的边界条件下求解的微分方程组的问 题变为泛函求极值的问题，而在求问题的近似解时，泛函的极值问题又可变成函 数的极值问題，因而最终把问題归结为求解线性代数方程组。
变分法在理论物理中非常重要：在拉格朗日力学中，以及在最小作用原理在 量子力学的应用中。
变分法提供了有限元方法的数学基础，它是求解边界值问题 的强力工具。
它们也在材料学中研究材料平衡量使用。
而在纯数学中的例子有， 黎曼在调和函数中使用狄力克雷原理。
应力变分法在力学领域同样拥有很髙的地位，这正说明了力学在学术界的重 要地位，通过应力变分法地学习，许多难题将更容易得到解答，所以，在以后的 学习生活中，我们将不会停止对力学的探究和学习，相信力学对我们的影响将是 巨大的。
参考文献：【1】弹性力学第四版徐芝纶髙等教育【2】弹性力学复习解题指导致 王俊民 同济大学[3] 弹性力学理论概要与典型题解王光钦西南交通大学[4] 弹性力学容精要与典型题解章军水利水电6禺 4(2f/2 + 6a 3)clx- FL: = 0 6& ■ 12(2f/2+66/3)Jx-FZ?=0E,。