弹性力学的求解方法和一般性原理

1、For personal use only in study and research; not for commercial use第五章 弹性力学的求解方法和一般性原理蚁知识点肅芅莅学习思路:芄 通过应力状态、 应变状态和本构关系的讨论， 已经建立了一系列的弹性 力学基本方程和边界条件。 本节的主要任务是将基本方程和边界条件作综合总结， 并且对求解方法作初步介绍。肀 弹性力学问题具有 15 个基本未知量， 基本方程也是 15 个，因此问题求 解归结为在给定的边界条件下求解偏微分方程。蚀 由于基本方程与 15 个未知量的内在联系，例如已知位移分量，通过几 何方程可以得到应变分量， 然后通过物

2、理方程可以得到应力分量； 反之，如果已 知应力分量， 也可通过物理方程得到应变分量， 再由几何方程的积分求出位移分 量，不过这时的应变分量必须满足一组补充方程， 即变形协调方程。 基于上述的 理由，为简化求解的难度，可以选取部分未知量作为基本未知量求解。膇 根据基本未知量， 弹性力学问题可以分为应力解法、 位移解法和混合解 法。肃 上述三种求解方法对应于偏微分方程的三种边值问题。膀学习要点：1.薄2.袁3.芀4.腿首先将弹性力学基本方程综合如下:芆1.平衡微分方程袄用张量形式描述%讥二°莀2.几何方程岛.=(1，- +% J薈用张量形式描述2'?变形协调方程蚃当然，具体求解弹

3、性力学问题时，并不需要同时求解十五个基本未知量，可以 而且必须做出必要的简化。 根据几何方程和本构方程可见， 位移、应力和应变分 量之间不是相互独立的。蒀 假如已知位移分量， 通过几何方程可以得到应变分量， 然后通过物理方 程可以得到应力分量。 反之，如果已知应力分量， 也可通过物理方程得到应变分 量，再由几何方程的积分求出位移分量， 不过这时的应变分量必须满足一组补充 方程，即变形协调方程。罿 基于上述的理由，为简化求解的难度，选取部分未知量作为基本未知 量。蒆 若以位移函数作为基本未知量求解，称为位移解法；蒂 若以应力函数作为基本未知量，称为应力解法；蕿 若以部分位移分量和部分应力分量作为

4、基本未知量，称为混合解 法。蒀 在给定的边界条件下，求解偏微分方程组的问题， 数学上称为偏微分 方程的边值问题。袄 按照不同的边界条件，弹性力学有三类边值问题。蒅 第一类边值问题： 已知弹性体内的体力Fbx，Fby，Fbz 和其表面的面力Fs< , Fsy , Fsz,求平衡状态的弹性体内各点的应力分量和位移分量，这时的边界 条件为面力边界条件。第二类边值问题：已知弹性体内的体力分量Fbx, Fby, Fbz以及表面的位移分量码匕W，求平衡状态的弹性体内各点的应力分量和位移分量，这时的边界条件为位移边界条件。薇第三类边值冋题：已知弹性体内的体力分量 Fbx， Fby， Fbz，以及物体表

5、面的部分位移分量和部分面力分量，求平衡状态的弹性体内各点的应力分量和位 移分量。这时的边界条件在面力已知的部分， 用面力边界条件，位移已知的部分 用位移边界条件，称为混合边值问题。蚆以上三类边值问题，代表了一些简化的实际工程问题。若不考虑物体的刚体位移，则三类边值问题的解是唯一的。芄4.边界条件：虿如果物体表面的面力Fsx, Fsy，Fsz为已知，则边界条件应为:羈称为面力边界条件，用张量符号表示为如果物体表面的位移见叱和已知，则边界条件应为羃称为位移边界条件。除了面力边界条件和位移边界条件，还有混合边界条件蝿 综上所述，弹性力学的基本未知量为三个位移分量， 六个应力分量和六 个应变分量，共计

6、十五个未知量。基本方程为三个平衡微分方程，六个几何方程 和六个物理方程，也是十五个基本方程。荿 这里没有考虑变形协调方程，原因是位移已经作为基本未知量。 对于任 意的单值连续的位移函数，如果设其有三阶的连续导数，则变形协调方程仅仅是 几何方程微分的结果，自然地满足，所以位移作为基本未知量时，不需要考虑变 形协调方程。螆要使基本方程有确定的解，还要有对应的面力或位移边界条件X,/5帝程就是在给定的边界条件下，就十五个未知量求解十五个基本方程。衿位移解法是以位移函数作为基本未知函数求解的，所以需要通过几何方程将位 移函数表达为应变分量，再通过物理方程将其表达为应力分量， 代入平衡微分方 程即可得到

7、位移解法的基本方程。螀首先，根据物理方程和几何方程，可以得到由位移分量表达的应力分量,即薇其中螅将上述位移表示的应力分量代入平衡微分方程，整理后可得V = + 罿这里是拉普拉斯运算符号，即於卽'曲祎上述方程是以位移表示的平衡微分方程，称为 拉梅（Lam）方程，它可以表示为张量形式羅或表达为矢量形式v =（F+ 7+）薃上式中 创 比 为拉普拉斯算符矢量聿对于边界条件，如果物体表面的位移已知，则直接由位移形式给定，即使用位 移边界条件。芇如果给定的边界条件是物体表面的面力，则面力边界条件式需用位移分量表示，将应力分量代入物理方程，整理可得位移分量表示的面力边界条件：蚇或表达为张量形式莂显

8、然，如果给定的边界条件是面力边界条件， 那么位移解法的边界条件表达式十分复杂，因此求解的难度将是比较大的。莃总之，如果以位移函数作为基本未知函数求解弹性力学问题，归结为在给定的边界条件下求解位移表示的平衡微分方程，即拉梅方程。蚈位移分量求解后，则可通过几何方程和物理方程求出相应的应变分量和应力分量。膅以应力作为基本未知函数求解弹性力学问题时，应力分量必须满足平衡微分方 程和面力边界条件。莅 但是仅此还不够，仅仅满足上述条件的应力分量并不是真正的应力。 因 为这组应力分量求出的应变分量代入几何方程， 将可能得到一组矛盾方程，不可 能求出单值连续的位移分量。要使这组方程不矛盾，则要求应力分量不仅满

9、足平 衡微分方程和面力边界条件，而且应力分量对应的应变分量必须满足变形协调方 程。蒂 这个问题也可以从物理上解释，应力分量满足平衡微分方程和面力边界 条件，只能保证物体的平衡，但是不能保证物体的连续。只有这组应力分量求出 的应变分量满足变形协调方程时，才能保证变形后的物体是连续的。聿 当位移分量作为基本未知函数求解时， 变形协调方程是自然满足的。如 果位移表示基本未知量，只有应力作为基本未知函数求解时，变形协调方程作为 一组补充方程是必须的。袇因此，对于应力解法，应力分量必须满足平衡微分方程和变形协调方程。膄由于变形协调方程是由应变分量表达的， 在应力解法中，需要将其转换为由应力分量表达。薂将

10、物理方程改写为蒀其中莅将上式代入变形协调方程的第一，四两式，可得轮换x，y，z可得其余四个方程。由此可得应力表达的变形协调方戈-IL”k 8a0蚂下面我们对应力分量表示的东少1 +询氐张叙卽氐作 简化蚇首先对平衡微分方程的第二和第三两式分别对 z, y求偏导数，然后相加可以得到肇将上式与变形协调方程的第二式相加后并整理，可得蚂上式为简化后的方程，轮换x，y，z以后，可得另外两个类似的公式螂综上所述，我们一共得到以下六个关系式:肇上述方程即为应力分量表达的变形协调方程，通常称为贝尔特拉米-米切尔方程。蒄如果弹性体体力为常量，则应力分量表达的变形协调方程可以简化为蚅上述方程为应力分量表达的变形协调

11、方程，通常简称为应力协调方程。但是应该注意：应力是不需要协调的，其实质仍为应变分量所满足的变形协调关 系。总而言之， 在以应力函数作为基本未知量求解时， 归结为在给定的 边界条件下，求解平衡微分方程和应力表达的变形协调方程所组成的偏微分方程 组。膆混合解法以六个应力分量和三个位移分量作为基本未知量求解弹性力学问题。 通过物理方程中消去应变分量， 其基本方程为平衡微分方程和由应力分量表达的 几何方程，即蒃 这里有三个平衡微分方程和六个几何方程， 共计九个方程对应九个未知 函数，加上给定的边界条件，则可得到唯一的解。袂 弹性力学的基本求解方法的应用要根据问题性质， 主要是根据边界条件 选择使用。衿

12、 对于面力边界条件问题，使用应力解法；蚄 位移边界条件应用位移解法；节 混合解法主要应用于混合边界条件， 即弹性体的部分边界位移已知， 部 分边界面力已知的问题。羂本节将从位移表达的平衡微分方程和应力表达的变形协调方程入手，推导体力 为常量时的应力分量， 应变分量， 位移分量，以及体积应力和体积应变所遵循的 规律，为进一步分析和理解弹性力学问题作必要的准备。羆 将位移分量表示的平衡微分方程的三个公式分别对 x， y，z 求偏导数， 然后相加可得莆由于肁所以（无十十“可玄日=0肂即莇由体积应力和体积应变的关系，可得v2e = o袄由上述公式可知，如果体力为常量，体积应力和体积应变均满足拉普拉斯（

13、Laplace）方程，即体积应力函数和体积应变函数均为调和函数。肄为了证明弹性力学解的唯一性定理，首先证明一个重要的定理，即应变能定理賺应变能定理是指：弹性体在外力作用下处于平衡状态时， 物体内存储的弹性势能，即应变能，等于外力由原始位置到平衡位置所做的功。 假如外力是由 零连续变化到其最终数值的，则在加载的过程中，物体始终是处于平衡状态 的。螈以下证明弹性体的应变能定理。设弹性体处于体力Fbx, Fby, Fbz和面力Fsc,Fsy,Fs.的作用下，弹性体内产生位移u，v，w。则外力在位移过程中作功为将面力边界条件代入上式的第二个积分，并利用高斯积分公式，可得袃因此芁由此可以证明，外力所做的

14、功等于弹性体存储的弹性势能。腿对于一般的工程构件，即弹性体，由于偏微分方程边值问题在数学上求解的困 难，因此直接根据给定的边界条件求解弹性力学的基本方程是十分困难的。羄为了避开偏微分方程边值问题直接求解的困难，在弹性力学问题的求解中，经常采用的方法是逆解法和半逆解法。薂逆解法就是根据研究问题的性质和研究对象特点， 确定基本未知量，写出相应的基本方程并且假设一组满足全部基本方程的应力函数或位移函数。然后在确定的坐标系下，考察具有确定的几何尺寸和形状的物体， 根据边界条件确定 表面作用面力或者已知位移。由此确定假设函数可以求解的弹性力学问 题。莁半逆解法就是对于给定的弹性力学问题，根据弹性体的几何

15、形状，受力特征和变形的特点或者已知的一些简单结论， 如材料力学得到的初等结论，假设 部分应力分量或者部分位移分量的函数形式为已知，由基本方程确定其他的未知 量，然后根据边界条件确定未知函数中的待定系数。薀逆解法和半逆解法的应用将在以后的章节中作详细介绍螆逆解法和半逆解法的求解过程带有"试算"的性质，显然弹性力学解的唯一性定理是逆解法和半逆解法的理论依据。蚅利用变形能定理，容易证明弹性力学 解的唯一性定理蒁 唯一性定理认为：假如弹性体内受已知体力的作用，表面受已知面力作用，或者表面位移为已知； 或者部分表面面力已知，部分表面位移已知。则弹性体处于平衡状态时，弹性体内任一点的应

16、力分量和应变 分量都是唯一的。如果弹性体表面有部分或全部位移已知，则位移分量也是唯一的。螇 下面我们证明解的唯一性定理蒈 假设在同一体力 Fbx, Fby, Fbz 的作用下， 并在同一边界条件下有两组不同的弹性力学解，有两组 不同的位移分量，应变分量和应力分量，即莄第一组为膈蒁第一组为袅为了证明这两组解相同，假设这两组解的差为一组新的解答，有膂由于第一组应力和第二组应力均为弹性力学的解，其应力应满足平衡微分方程。因此，两组平衡微分方程相减可得薁 因此，第三组应力满足体力为零的平衡微分方程。薈解的迭加原理 ：弹性力学解的迭加原理是指在线弹性条件下，对于满足小变形 条件的弹性体，在两组不同的外力

17、作用下所得到的弹性力学解相加等于这两组外 力共同作用于弹性体的解答。设弹性体在第一组外力，体力Fbx , Fby, Fbz和面力Fsx,Fsy, Fsz作用下，产生的应力分量为匚j；蚁同一弹性体在第二组外力，体力F'bx， F'by，F'bz和面力F'sx， F'sy，F'sz作用下，产生的应力为Jj。罿显然上述两组外力所产生的两组应力均应满足平衡微分方程和应力表示的变形协调方程，以及对应的边界条件。由于基本方程和边界条件都是线性的， 因此迭加后的应力显然满足两组外力共同作用时的平衡微分方程和边界条件。即迭加后的应力满足全部基本方程和边界条件，是

18、两组外力共同作用时的弹性力学 解答。因此可以证明弹性力学解满足迭加原理。肅应该注意的是由于全部基本方程和边界条件是由变形前的坐标描述的，因此只有在小变形的条件下才可以使用迭加原理。即变形对外力作用点位置的改 变可以忽略不计。羄弹性力学问题的求解是在给定的边界条件下求解基本方程。由于偏微分方程边 值问题的性质，弹性力学的解必然要求物体表面的外力或者位移是按一定的规律 分布的。对于工程实际问题，构件表面面力或者位移是很难满足这个要求的。 这 使得弹性力学解的应用将受到极大的限制。螁为了扩大弹性力学解的适用范围，放宽这种限制，圣维南提出了局部影响原理。肆 圣维南局部影响原理 其主要内容为：物体表面某

19、一小面积上作用的外力 力系，如果被一个静力等效的力系所替带， 那么物体内部只能导致局部应力的改 变。而在距离力的作用点较远处，其影响可以忽略不计。螇根据圣维南局部影响原理，假如我们用一静力等效力系取代弹性体上作用的原外力，则其影响仅在力的作用区域附近。离此区域较远处，几乎不受影响'螃 局部影响原理是在实验中被多方证明的。 例如用一个钳子夹住铁杆， 钳 子对铁杆的作用相当于一组平衡力系。 实验证明， 无论作用力多大， 在距离力的 作用区域比较远处，几乎没有应力产生。袀 这一结论也可以通过有限元程序计算分析得到。对于矩形板， 在等效力系一个集中力F；两个集中力F/2和分布载荷F/I的作用，

20、计算数据表明：静力 等效的力系只能导致弹性体局部应力的改变。 而在距离力的作用点较远处， 其影 响可以忽略不计。蒇 圣维南局部影响原理也可以表述为： 如果物体的任一部分作用一个平衡 力系，则此平衡力系在物体内产生的应力分布， 仅局限于该力系作用的局部区域， 在离该区域比较远处，这种影响便急剧减少。芅学习思路:蒂 通过应力状态、 应变状态和本构关系的讨论， 已经建立了一系列的弹性 力学基本方程和边界条件。 本节的主要任务是将基本方程和边界条件作综合总结， 并且对求解方法作初步介绍。羀弹性力学问题具有 15 个基本未知量， 基本方程也是 15个，因此问题求解归结为在给定的边界条件下求解偏微分方程。

21、袈由于基本方程与 15 个未知量的内在联系，例如已知位移分量，通过几何方程可以得到应变分量， 然后通过物理方程可以得到应力分量； 反之，如果已 知应力分量， 也可通过物理方程得到应变分量， 再由几何方程的积分求出位移分 量，不过这时的应变分量必须满足一组补充方程， 即变形协调方程。 基于上述的 理由，为简化求解的难度，可以选取部分未知量作为基本未知量求解。羇根据基本未知量， 弹性力学问题可以分为应力解法、 位移解法和混合解法。上述三种求解方法对应于偏微分方程的三种边值问题。肀学习要点：艿 1.肆 2.羅 3.螁 4.莁学习思路:袈 以位移函数作为基本未知量求解弹性力学问题的方法称为位移法。螄

22、位移解法的基本方程是位移表示的平衡微分方程。 位移分量求解后， 则 可以通过几何方程和物理方程求出相应的应变分量和应力分量。袁 如果问题的边界条件为位移边界条件， 边界条件描述比较简单。 如果问 题为面力边界条件，由于边界条件是通过位移函数的导数描述的， 因此应用困难。蒈 总之若以位移为基本未知函数求解时，归结为在给定的边界条件下求解位移 表示的平衡微分方程，即拉梅方程。芆学习要点：薃 1.羁 2.衿 3. 。羇学习思路 :薆 如果选用应力分量或者应力函数作为基本未知量求解弹性力学问题称 为应力解法。肁 应力解法的基本方程不仅有平衡微分方程， 而且有变形协调方程。 因为 仅仅满足平衡微分方程的

23、应力分量并不一定是真实应力， 这组应力分量求出的应 变分量代入几何方程， 将可能得到一组矛盾方程， 这就不可能求出单值连续的位 移分量。芀 由于变形协调方程是应变表示的，在应力解法中，需要转化为基本未知量应 力分量表示。蒅 利用平衡微分方程的求导形式简化变形协调方程， 可以得到应力分量表示的 变形协调方程。莄 总之，在以应力函数作为基本未知量求解时，归结为在给定的边界条件下， 求解平衡微分方程和应力表达的变形协调方程所组成的偏微分方程。膁学习要点：蚀 1.膇 2.膃 3.芁 4.肁学习思路:薅 如果选取应力分量和位移分量作为基本未知量求解弹性力学问题， 称为 混合解法。膆 基本方程为平衡微分方

24、程和应力分量表达的几何方程。芁 混合解法三个平衡微分方程和六个几何方程， 共计九个方程对应九个未知函 数，加上给定的边界条件，则可得到唯一的解。芈学习要点：1.羅学习思路 :莀 本节讨论体力为常量条件下， 弹性力学的基本未知量的特性。 通过这些 物理量性质的研究，将有助于今后进一步探讨弹性力学问题。虿 从位移表达的平衡微分方程和应力表达的变形协调方程入手， 可以得到 体积应力函数和体积应变函数均为调和函数； 而位移分量， 应变分量和应力分量 均为双调和函数。聿学习要点：蚄 1.蒀 2.肀 3. 。蒇学习思路 :蒃 本节通过应变能原理推导弹性力学的解的唯一性定理薀 弹性力学解的唯一性定理： 假如

25、弹性体内受已知体力的作用， 物体表面 面力已知，或者表面位移已知；或者部分表面面力已知，部分表面位移已知。则 弹性体处于平衡状态时， 弹性体内任一点的应力分量和应变分量都是唯一的。 对 于表面有部分或全部位移已知的，则位移分量也是唯一的。解的唯一性定理证明可以根据不同的方法证明艿 解的唯一性定理的意义是为弹性力学问题的求解提供了重要的理论依据。 由 于偏微分方程边值问题求解的困难， 因此在弹性力学问题分析中， 经常需要使用 逆解法或半逆解法。而解的唯一性定理为这些方法奠定了基础。学习要点：蚀 1.2.蚆 3.芅学习思路 :螀 由于偏微分方程边值问题求解困难， 因此直接由给定的边界条件求解弹 性

26、力学的基本方程几乎是不可能的。 所以对于弹性力学问题的求解， 经常采用的 方法是逆解法和半逆解法。逆解法 就是根据问题的性质，确定基本未知量，写出相应的基本方程并且假 设一组满足全部基本方程的应力函数或位移函数。 然后在确定的坐标系下， 考察 具有确定的几何尺寸和形状的物体， 其表面将受什么样的面力作用或者将有什么 样的位移，然后确定假设的函数。半逆解法 就是对于给定的弹性力学问题，根据弹性体的几何形状，受力特征 和变形的特点或已知的一些简单结论， 如材料力学得到的初等结论， 假设部分应 力分量或者部分位移分量的函数形式为已知， 由基本方程确定其他的未知量， 然 后根据边界条件确定未知函数中的待定系数。逆解法和半逆解法将在以后的章节中作介绍。其求解过程带有“试算”的性质，显然