弹性力学用差分法及变分法解平面问题课件

用差分法和变分法解平面问题用差分法和变分法解平面问题章章弹性力学弹性力学弹性力学弹性力学一、差分公式的推导一、差分公式的推导二、弹性体的形变势能和外力势能二、弹性体的形变势能和外力势能三、位移变分方程三、位移变分方程四、位移变分法四、位移变分法五、位移变分法例题五、位移变分法例题第五章用差分法和变分法解平面问题第五章用差分法和变分法解平面问题内容提要内容提要弹性力学简明教程(第三版)徐芝纶院士(1911-1999)弹性力学弹性力学 弹性力学的基本解法是，根据静力平衡条件，形变与位移弹性力学的基本解法是，根据静力平衡条件，形变与位移之间的几何条件和形变与应力之间的物理条件，建立微分方程之间的几何条件和形变与应力之间的物理条件，建立微分方程和边界条件。和边界条件。因此，弹性力学问题属于微分方程的边值问题。通过求解，得出函数表示的精确解答。对于工程实际问题，由于荷载和边界较复杂，难以求出函数式的解答。为此，人们探讨弹性力学的各种近似解法，主要有差分法、变分法和有限单元法。差分公式的推导差分公式的推导一一弹性力学弹性力学3用差分法和变分法解平面问题用差分法和变分法解平面问题fxo差分法差分法是微分方程的一种数值解法，是微分方程的一种数值解法，它不是去求解函数它不是去求解函数f(x)，而是求函数在一些结点上的值，而是求函数在一些结点上的值 。将导数用有限差商来代替将导数用有限差商来代替将微分用有限差分来代替将微分用有限差分来代替将微分方程用差分方程（代数方程）代替，求解微分方程问将微分方程用差分方程（代数方程）代替，求解微分方程问题化为求解差分方程问题。题化为求解差分方程问题。差分公式的推导差分公式的推导一一弹性力学弹性力学4用差分法和变分法解平面问题用差分法和变分法解平面问题在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格。设()为弹性体内的某一个连续函数，该函数在平行于x轴的一根网线上，如在3-上，它只随x坐标的改变而变化。在邻近结点处，函数f可展为泰勒级数如下：差分公式的推导差分公式的推导一一(a)弹性力学弹性力学5用差分法和变分法解平面问题用差分法和变分法解平面问题 只考虑离开结点充分近的那些结点，即（0）充分小。于是可不计（0）的三次及更高次幂的各项，则上式简写为：在结点，0；在结点1，0。代入(b)得：联立（c）、（d），解得差分公式：差分公式的推导差分公式的推导一一(b)(c)(d)弹性力学弹性力学6用差分法和变分法解平面问题用差分法和变分法解平面问题同理，在网线4-0-2上可得到差分公式：以上（）（）是基本差分公式，从而可导出其它的差分公式如下：差分公式的推导差分公式的推导一一弹性力学弹性力学7用差分法和变分法解平面问题用差分法和变分法解平面问题差分公式()及()是以相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称为端点导数公式。应当指出：中点导数公式与端点导数公式相比，精度较高。因为前者反映了结点两边的函数变化，而后者却只反映了结点一边的函数变化。因此，我们总是尽可能应用前者，而只有在无法应用前者时才不得不应用后者。差分公式的推导差分公式的推导一一弹性力学弹性力学8用差分法和变分法解平面问题用差分法和变分法解平面问题一、差分公式的推导一、差分公式的推导二、弹性体的形变势能和外力势能二、弹性体的形变势能和外力势能三、位移变分方程三、位移变分方程四、位移变分法四、位移变分法五、位移变分法例题五、位移变分法例题第五章用差分法和变分法解平面问题第五章用差分法和变分法解平面问题内容提要内容提要弹性力学简明教程(第三版)徐芝纶院士(1911-1999)弹性力学弹性力学弹性力学变分法，因其泛函就是弹性体的能量（如形变势能、弹性力学变分法，因其泛函就是弹性体的能量（如形变势能、外力势能），又称为能量法。外力势能），又称为能量法。泛函是以函数为自变量的一类函数。泛函是以函数为自变量的一类函数。变分法，是研究泛函及其极值的求解方法。变分法，是研究泛函及其极值的求解方法。弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二弹性力学变分法，是区别于微分方程边值问题的另一种独弹性力学变分法，是区别于微分方程边值问题的另一种独立解法，分为：立解法，分为：位移变分法：取位移函数为自变量，并以势能极小值条件导出位移变分法：取位移函数为自变量，并以势能极小值条件导出变分方程。变分方程。应力变分法：取应力函数为自变量，并以余能极小值条件导出应力变分法：取应力函数为自变量，并以余能极小值条件导出变分方程。变分方程。弹性力学弹性力学10用差分法和变分法解平面问题用差分法和变分法解平面问题（2）因应力和应变均从）因应力和应变均从0增长到增长到 ，故单位体积上，故单位体积上，应力所做的功是应力所做的功是 非线性非线性 关系关系 线线 性性 关系关系（1）作用于微小单元上的应力，是邻近部分物体对它的作用力，）作用于微小单元上的应力，是邻近部分物体对它的作用力，可看成是作用于微小单元上的可看成是作用于微小单元上的“外力外力”。、应力的功和形变势能（内力势能）、应力的功和形变势能（内力势能）弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二弹性力学弹性力学11用差分法和变分法解平面问题用差分法和变分法解平面问题线性的应力与应变关系非线性的应力与应变关系弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二、应力的功和形变势能（内力势能）、应力的功和形变势能（内力势能）弹性力学弹性力学12用差分法和变分法解平面问题用差分法和变分法解平面问题（3）对于平面应力问题 或平面应变问题 单位体积上应力所做的功都是 (c)弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二、应力的功和形变势能（内力势能）、应力的功和形变势能（内力势能）弹性力学弹性力学13用差分法和变分法解平面问题用差分法和变分法解平面问题弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二、应力的功和形变势能（内力势能）、应力的功和形变势能（内力势能）（4）假设没有转化为非机械能和动能，则应力所做的功全部）假设没有转化为非机械能和动能，则应力所做的功全部转化为弹性体的内力势能，又称为形变势能，或应变能，转化为弹性体的内力势能，又称为形变势能，或应变能，存贮于物体内部。存贮于物体内部。单位体积的形变势能（形变势能密度）。单位体积的形变势能（形变势能密度）。（5）整个弹性体的形变势能）整个弹性体的形变势能(d)弹性力学弹性力学14用差分法和变分法解平面问题用差分法和变分法解平面问题（6）将物理方程代入，平面应力问题的形变势能密度）将物理方程代入，平面应力问题的形变势能密度 ，可，可用形变表示为用形变表示为对于平面应变问题，对于平面应变问题，将将再将几何方程代入，再将几何方程代入，可用位移表示为可用位移表示为整个弹性体的形变势能为整个弹性体的形变势能为弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二、应力的功和形变势能（内力势能）、应力的功和形变势能（内力势能）(5-16)弹性力学弹性力学15用差分法和变分法解平面问题用差分法和变分法解平面问题（1 1）U U 是应变或位移的二次泛函，故不能应用叠加原理。是应变或位移的二次泛函，故不能应用叠加原理。（2 2）应变或位移发生时，）应变或位移发生时，U U 总是正的，即总是正的，即 （3 3）U U 的大小与受力次序无关。的大小与受力次序无关。（4 4）对应变的导数，等于对应的应力：对应变的导数，等于对应的应力：(5-15)弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二2 2、形变势能的性质、形变势能的性质弹性体每单位体积中的形变势能对于任一形变分量的弹性体每单位体积中的形变势能对于任一形变分量的改变率，等于相应的应力分量。改变率，等于相应的应力分量。弹性力学弹性力学16用差分法和变分法解平面问题用差分法和变分法解平面问题外力势能外力势能外力做了功，必然消耗了相同外力做了功，必然消耗了相同 值的势能。当取值的势能。当取 时的外力功和能时的外力功和能为零，则：为零，则：(b)外力功：外力功：弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二3 3、弹性体上的外力功和外力势能、弹性体上的外力功和外力势能 弹性力学弹性力学17用差分法和变分法解平面问题用差分法和变分法解平面问题弹性体的总势能，是外力势能和内力（形变）势能之和，(h)弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二4 4、弹性体的总势能、弹性体的总势能弹性力学弹性力学18用差分法和变分法解平面问题用差分法和变分法解平面问题弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二例题例题1 试证明，在同样的应变分量试证明，在同样的应变分量 下，平面应变情况下下，平面应变情况下单位厚度的形变势能大于平面应力情况下的形变势能。单位厚度的形变势能大于平面应力情况下的形变势能。对于平面应变情况，只需将上式中对于平面应变情况，只需将上式中 ，变换为变换为解：平面应力情况下，单位厚度的形变势能：解：平面应力情况下，单位厚度的形变势能：(a)弹性力学弹性力学19用差分法和变分法解平面问题用差分法和变分法解平面问题代入，得显然，方括号内将式(a)中的，都作为式(b)的变换，整理后得平面应变情况下的形变势能公式，(c)弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二弹性力学弹性力学20用差分法和变分法解平面问题用差分法和变分法解平面问题从式(c)可见，在平面应变情况下，形变势能中的第1，2，3项均大于平面应力情况下的值，而第4项不变。因此，平面应变的形变势能大于平面应力的形变势能U。弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二弹性力学弹性力学21用差分法和变分法解平面问题用差分法和变分法解平面问题CDEFAB弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二例题例题2图示一板块，在铅直方向均布拉力作用下发生拉伸变形，并图示一板块，在铅直方向均布拉力作用下发生拉伸变形，并使之两端固定下来，若在其中切开一小口时，试说明板的形变势使之两端固定下来，若在其中切开一小口时，试说明板的形变势能将发生什么变化？能将发生什么变化？解：解：当线切开时，线上的应力趋于当线切开时，线上的应力趋于0，而形变势能是正定，而形变势能是正定，当应力，当应力 时，相应的形变势能也失去。因此，时，相应的形变势能也失去。因此，板的总的形变势能减少。板的总的形变势能减少。弹性力学弹性力学22用差分法和变分法解平面问题用差分法和变分法解平面问题弹性体的形变势能和外力势能弹性体的形变势能和外力势能二二 当线切开后，边界和仍是固定的，我们可以比较两种状态：当线切开后，边界和仍是固定的，我们可以比较两种状态：切开后，切开后，线仍然处于闭合状态，不发生张开，这是不稳定线仍然处于闭合状态，不发生张开，这是不稳定的平衡状态。的平衡状态。线张开，出现裂纹，这是稳定的平衡状态。由于系统的稳线张开，出现裂纹，这是稳定的平衡状态。由于系统的稳定平衡状态与邻近的状态相比，总势能处于极小值，而定平衡状态与邻近的状态相比，总势能处于极小值，而(a)，(b)两种状态的外力势能不变，因此，两种状态的外力势能不变，因此，(b)的形变势能小于的形变势能小于(a)，即形变势能将减少。，即形变势能将减少。弹性力学弹性力学23用差分法和变分法解平面问题用差分法和变分法解平面问题一、差分公式的推导一、差分公式的推导二、弹性体的形变势能和外力势能二、弹性体的形变势能和外力势能三、位移变分方程三、位移变分方程四、位移变分法四、位移变分法五、位移变分法例题五、位移变分法例题第五章用差分法和变分法解平面问题第五章用差分法和变分法解平面问题内容提要内容提要弹性力学简明教程(第三版)徐芝纶院士(1911-1999)弹性力学弹性力学1.1.实际平衡状态的位移实际平衡状态的位移 ，必须满足，必须满足用位移表示的平衡微分方程（在A中）;用位移表示的应力边界条件（在上);位移边界条件（在上）。(a)其中，属于静力平衡条件，属于约束条件。对于实际位移，可将看成是必要条件，而，是充分条件。位移变分方程位移变分方程 三三弹性力学弹性力学25用差分法和变分法解平面问题用差分法和变分法解平面问题（在 上）。2.2.虚位移状态虚位移状态 虚位移（数学上称为位移变分）虚位移（数学上称为位移变分），表示在约束条件允许下，平衡状态附近的微表示在约束条件允许下，平衡状态附近的微小位移增量，虚位移应满足小位移增量，虚位移应满足 上的约束边上的约束边界条件，即界条件，即 (b)位移变分方程位移变分方程 三三弹性力学弹性力学26用差分法和变分法解平面问题用差分法和变分法解平面问题虚位移不是实际外力作用下发生的，而是假想由其他干扰产生的。因此，虚位移状态就构成实际平衡状态附近的一种邻近状态。(c)位移变分方程位移变分方程 三三弹性力学弹性力学27用差分法和变分法解平面问题用差分法和变分法解平面问题(d)变分与微分的比较变分与微分的比较位移变分方程位移变分方程 三三微分是在同一状态下，研究由于位置微分是在同一状态下，研究由于位置(坐标坐标)改变而引起函数改变而引起函数的改变。其中的自变量为坐标变量，而因变量为函数，如的改变。其中的自变量为坐标变量，而因变量为函数，如位移，有位移，有变分是在同一点位置上，由于状态改变而引起泛函的改变。变分是在同一点位置上，由于状态改变而引起泛函的改变。其中的自变量为状态函数，如位移；而因变量为泛函，如其中的自变量为状态函数，如位移；而因变量为泛函，如 ，有，有(e)弹性力学弹性力学28用差分法和变分法解平面问题用差分法和变分法解平面问题由于微分和变分都是微量，所以 a.它们的运算方式相同，如式(d)，(e)；b.变分和微分可以交换次序，如 (f)位移变分方程位移变分方程 三三弹性力学弹性力学29用差分法和变分法解平面问题用差分法和变分法解平面问题当发生虚位移（位移变分）时，由于虚位移引起虚应变，外力势能的变分：外力势能的变分：外力的虚功（外力功的变分）：外力的虚功（外力功的变分）：3.3.在虚位移上弹性体的功和能在虚位移上弹性体的功和能 位移变分方程位移变分方程 三三弹性力学弹性力学30用差分法和变分法解平面问题用差分法和变分法解平面问题形变势能的变分，即实际应力在虚应变上的虚功,由于实际应力在虚应变之前已存在,所以作为常力计算，故无系数。(j)位移变分方程位移变分方程 三三弹性力学弹性力学31用差分法和变分法解平面问题用差分法和变分法解平面问题（1）在封闭系统中，假设没有非机械能的改变，也没有动能的改变，则按照能量守恒定律，在虚位移过程中形变势能的增加 应等于外力势能的减少（即等于外力所做的虚功 ）。所以4.4.弹性力学中位移变分方程的导出弹性力学中位移变分方程的导出 位移变分方程位移变分方程 三三弹性力学弹性力学32用差分法和变分法解平面问题用差分法和变分法解平面问题（2）位移变分方程将式(g)的代入上式，得它表示，在实际平衡状态发生位移的变分时，所引起的形变势能的变分，等于外力功的变分。位移变分方程位移变分方程 三三(5-22)弹性力学弹性力学33用差分法和变分法解平面问题用差分法和变分法解平面问题（3）虚功方程 将式(j)的 代入上 式，得位移变分方程位移变分方程 三三虚功方程表示：如果在虚位移发生之前，弹性体处于虚功方程表示：如果在虚位移发生之前，弹性体处于平衡状态，则在虚位移过程中，外力在虚位移上所做的虚平衡状态，则在虚位移过程中，外力在虚位移上所做的虚功等于应力在虚应变上所做的虚功。功等于应力在虚应变上所做的虚功。(5-24)弹性力学弹性力学34用差分法和变分法解平面问题用差分法和变分法解平面问题其中 形变势能的变分，如式(j)所示，外力功的变分,如式(g)所示。（4）最小势能原理式(k)可写成其中U弹性体的形变势能，如5-4式(d),W弹性体的外力功，如5-4式(a)。可以证明，式(n)可以写成为 位移变分方程位移变分方程 三三弹性力学弹性力学35用差分法和变分法解平面问题用差分法和变分法解平面问题由于弹性体的总势能为故式(o)可以表示为再将总势能对其变量（位移或应变）作二次变分运算，可得(p)(q)位移变分方程位移变分方程 三三极小势能原理：在给定的外力作用下，在满足位移边界条极小势能原理：在给定的外力作用下，在满足位移边界条件的所有各组位移状态中，实际存在的一组位移对应于总势能件的所有各组位移状态中，实际存在的一组位移对应于总势能为极小值。为极小值。弹性力学弹性力学36用差分法和变分法解平面问题用差分法和变分法解平面问题一、差分公式的推导一、差分公式的推导二、弹性体的形变势能和外力势能二、弹性体的形变势能和外力势能三、位移变分方程三、位移变分方程四、位移变分法四、位移变分法五、位移变分法例题五、位移变分法例题第五章用差分法和变分法解平面问题第五章用差分法和变分法解平面问题内容提要内容提要弹性力学简明教程(第三版)徐芝纶院士(1911-1999)弹性力学弹性力学 位移变分法是取位移为基本未知函数的。位移函数应预先满足 上的位移边界条件，然后再满足位移变分方程。位移变分法位移变分法 四四弹性力学弹性力学38用差分法和变分法解平面问题用差分法和变分法解平面问题(a)（1）因位移函数是未知的，在变分法中采用设定位移试函数的方法，令位移变分法位移变分法 四四弹性力学弹性力学39用差分法和变分法解平面问题用差分法和变分法解平面问题其中 和 均为设定的x，y的函数，并在边界 上，令 （在 上）（在 上）(c)(b)位移变分法位移变分法 四四弹性力学弹性力学40用差分法和变分法解平面问题用差分法和变分法解平面问题 所以 已满足了 上的位移边界条件。而 ，用来反映位移状态的变化，故位移的变分为(d)位移变分法位移变分法 四四弹性力学弹性力学41用差分法和变分法解平面问题用差分法和变分法解平面问题位移的变分通过，的变分来反映，故形变势能的变分为（2）位移(a)还必须满足位移变分方程(d)将式(d)，(f)代入(e)得位移变分法位移变分法 四四弹性力学弹性力学42用差分法和变分法解平面问题用差分法和变分法解平面问题因虚位移（位移变分）中的，是完全任意的，独立的，为了满足上式，必须:位移变分法位移变分法 四四弹性力学弹性力学43用差分法和变分法解平面问题用差分法和变分法解平面问题式(g)是瑞利-里茨变分方程。它是关于，的线性代数方程组，由上式可解出，,从而得到位移的解答。位移变分法位移变分法 四四弹性力学弹性力学44用差分法和变分法解平面问题用差分法和变分法解平面问题一、差分公式的推导一、差分公式的推导二、弹性体的形变势能和外力势能二、弹性体的形变势能和外力势能三、位移变分方程三、位移变分方程四、位移变分法四、位移变分法五、位移变分法例题五、位移变分法例题第五章用差分法和变分法解平面问题第五章用差分法和变分法解平面问题内容提要内容提要弹性力学简明教程(第三版)徐芝纶院士(1911-1999)弹性力学弹性力学例例1 1 图示矩形图示矩形板板abab，在上，在上边及右边受有边及右边受有均布压力均布压力 及及 ，而左边和下，而左边和下边受有法向连边受有法向连杆的约束。杆的约束。位移变分法例题位移变分法例题 五五弹性力学弹性力学46用差分法和变分法解平面问题用差分法和变分法解平面问题位移边界(a)应用瑞利应用瑞利-里茨法里茨法，设定位移，设定位移只取 和位移变分法例题位移变分法例题 五五弹性力学弹性力学47用差分法和变分法解平面问题用差分法和变分法解平面问题 在平面应力状态下可得即位移变分法例题位移变分法例题 五五弹性力学弹性力学48用差分法和变分法解平面问题用差分法和变分法解平面问题由可得即解得位移变分法例题位移变分法例题 五五弹性力学弹性力学49用差分法和变分法解平面问题用差分法和变分法解平面问题(e)例例2 2本题全部为位移边界条件：位移变分法例题位移变分法例题 五五弹性力学弹性力学50用差分法和变分法解平面问题用差分法和变分法解平面问题本题以y轴为对称轴，所以u应为x的奇函数,v应为x的偶函数。(f)设定位移势函数为设定位移势函数为位移变分法例题位移变分法例题 五五弹性力学弹性力学51用差分法和变分法解平面问题用差分法和变分法解平面问题位移(g)已满足对称性条件(f)和全部边界条件(e)。因全部为位移边界条件且均已满足，所以从55式(u)可见，也可应用伽辽金变分法。位移变分法例题位移变分法例题 五五弹性力学弹性力学52用差分法和变分法解平面问题用差分法和变分法解平面问题将位移(g)代入上式，求出得出的位移解答与书中用瑞利-里茨法给出的结果相同。因，故伽辽金变分方程伽辽金变分方程为(h)位移变分法例题位移变分法例题 五五弹性力学弹性力学53用差分法和变分法解平面问题用差分法和变分法解平面问题