清华大学弹性力学冯西桥FXQChapter微分提法.pptx

《清华大学弹性力学冯西桥FXQChapter微分提法.pptx》由会员分享，可在线阅读，更多相关《清华大学弹性力学冯西桥FXQChapter微分提法.pptx（86页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、弹性理论的微分提法 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第1页/共86页微分提法Chapter 6.1平衡方程(Navier)：几何方程(Cauchy)：应变协调方程：(Saint-Venant)本构方程：(1)应变应力公式：(Hooke)(2)应力应变公式：(Lam)第3页/共86页微分提法Chapter 6.1当选位移作基本量时只需考虑几何方程，协调方程将自动满足；当选应变作基本量时，只需满足协调方程，就能保证由几何方程积分出单值连续的位移场来。两个本构方程也是等价的，于是有两组基本方程组：第4页/共86页 第一组

2、 基本未知量：ij(6),ij(6),ui (3)平衡方程：(3)几何方程：(6)应力应变关系：(6)微分提法Chapter 6.1第5页/共86页 第二组 基本未知量：ij(6),ij(6)平衡方程：(3)协调方程：(6-3)应变应力关系：(6)微分提法Chapter 6.1第6页/共86页微分提法Chapter 6.1 边界条件：第7页/共86页微分提法Chapter 6.1 弹性理论中常见的三种边界条件：1.处处给定外部作用力 的力边界条件S。边界条件为：域内应力场的边界值应满足柯西公式 p不能消除刚体位移p要满足整体平衡条件。第8页/共86页微分提法Chapter 6.1分量形式为：当

3、 时称为自由表面，是力边界的特殊情况。集中力化为作用在微小面积上的均布表面力。集中力矩则化为非均布表面力。第9页/共86页微分提法Chapter 6.12.处处给定位移 约束的位移边界Su。域内位移场的边界值应等于给定边界值：有时也可指定边界位移的导数值(例如，转角为零)或应变值。在静力学问题中，所给的位移应足以防止物体的刚体运动。第10页/共86页微分提法Chapter 6.13.在部分边界S 上给定外力，部分边界Su上给定位移的混合边界S。这时要求 对于弹性动力学问题，还应给定初始条件。第11页/共86页4.混合型边界条件 在边界同一位置，给定部分位移分量和部分面力分量。51chapter

4、 3.6微分提法5.弹簧类边界条件第12页/共86页6.对称和反对称条件51chapter 3.6微分提法p对称载荷：在对称面上，所有对称场变量的一阶导数等于零，所有反对称场变量的值等于零。p反对称载荷：在对称面上，所有反对称场变量的一阶导数等于零，所有对称场变量的值等于零。第13页/共86页56chapter 3.6微分提法6.对称和反对称条件第14页/共86页微分提法Chapter 6.1 弹性力学问题微分提法的基本思想：从研究弹性体内的微元入手，导出描述微元静力平衡、变形几何及弹性关系的一组基本方程，加上相应的边界条件把弹性力学问题归结为求解偏微分方程组的边值问题。第15页/共86页微分

5、提法Chapter 6.1从求解的未知量方面考虑，可分为如下四类：位移为基本未知量 应变为基本未知量 应力（应力函数）为基本未知量 混合未知量第18页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第19页/共86页位移解法Chapter 6.2 位移解法是以位移分量ui作基本未知量的解法。即以位移分量的三个未知函数作为基本未知函数。这三个位移分量所对应的应力在物体内部应满足平衡微分方程。现经过下述步骤将平衡微分方程中的应力改用位移表示，从而得出用位移表示的平衡微分方程式。第20页/共86页位移

6、解法Chapter 6.2用位移表示的平衡方程（Lam-Navier方程）几何几何方程方程本构本构关系关系平衡平衡方程方程第21页/共86页位移解法Chapter 6.2 具体推导如下：先将几何关系代入广义虎克定律，可得 式中第22页/共86页位移解法Chapter 6.2第23页/共86页位移解法Chapter 6.2代入第一个以位移表示的平衡微分方程 第24页/共86页位移解法Chapter 6.2同样可得其余两个方程，即式中第25页/共86页位移解法Chapter 6.2 上式实质上是位移形式的平衡方程式，这就是位移法的基本方程式。综上指标形式第26页/共86页位移解法Chapter 6

7、.2 边界条件 若给定的是位移边界条件，则直接用位移表示，即若给定的是表面力的边界条件，则可将其表面力以位移表示(以x方向为例)，第27页/共86页位移解法Chapter 6.2代入第28页/共86页位移解法Chapter 6.2 用位移表示的外力边界条件：第29页/共86页位移解法Chapter 6.2u 微分方程的解：齐次方程通解特解（易得）u 齐次的Lam-Navier方程(即fi0的无体力情况)：将齐次方程对xi求导，并对指标i迭加后得而第30页/共86页位移解法Chapter 6.2是非零常数，故第一应变不变量应满足调和方程 其中 称为调和算子或拉普拉斯算子。第31页/共86页 根据

8、 ，其中K为常数。故第一应力不变量 (或平均正应力)也满足调和方程：上式作调和运算得：位移解法Chapter 6.2即即第32页/共86页位移解法Chapter 6.2由连续性条件及式(1)式得又其中 称为重调和算子。上式说明位移分量ui应满足重调和方程。第33页/共86页位移解法Chapter 6.2连续性条件 代入于是说明应力及应变分量也都满足重调和方程。第34页/共86页位移解法Chapter 6.2 弹性力学问题解的性质u 调和函数的性质（见下册p43）(1)调和函数的各阶导数均为调和函数(2)若 为调和函数，则 也是调和函数u 调和函数和双调和函数的关系(1)若 为调和函数，则 也是

9、双调和函数(2)若 为调和函数，则 是双调和函数(3)若 为调和函数，则 是双调和函数第35页/共86页位移解法Chapter 6.2 综上所述，在无体力情况下，第一应变不变量、第一应力不变量和平均正应力0都是调和函数。位移分量ui，应变分量ij和应力分量ij都是重调和函数。于是弹性力学的无体力问题在数学上归结为调和方程和重调和方程的边值问题。对于常体力情况ficonst，不难验证这个结论同样适用。对于变体力情况，可先找一个特解(不必满足边界条件)，然后与上述齐次解迭加，使全解满足全部边界条件。第36页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法

10、 叠加原理 解的唯一性定理 圣维南原理Chapter 6第37页/共86页应力解法Chapter 6.3 平衡方程：(3)协调方程：(6-3)应变应力关系：(6)应力解法是以应力分量作基本未知量的解法。第38页/共86页应力解法Chapter 6.3Beltrami-Michell方程：本构关系本构关系代入协调方程代入协调方程利用平衡方程利用平衡方程u 思路：消掉其余基本量，仅用应力表示：这就是应力解法的定解方程，称为应力协调方程或贝尔脱拉密-密乞尔方程，简称B-M方程，共含六个二阶椭圆方程。E.Beltrami(1835-1900)第39页/共86页应力解法Chapter 6.3分 量 形

11、式第40页/共86页应力解法Chapter 6.3 前面曾指出，六个应变协调方程并不完全独立，不能由它们独立解出六个应变分量。以此类推，六个应力协调方程也不可能完全独立，所以用应力解法解题时通常要求在域内同时满足六个B-M方程和三个平衡方程，且在边界上满足三个力边界条件。对于全部边界给定外力的边值问题，应力解法可以避开几何关系直接解出工程中关心的应力分量。但应力解法处理位移边界条件相当困难。应力解法涉及六个二阶B-M方程，三个一阶平衡方程和三个力边界条件，对于几何形状或载荷分布较复杂的问题求解比较困难。第41页/共86页应力解法Chapter 6.3 对于动力学问题，应把惯性力纳入体力项，B-

12、M方程中的fi,j应改为 于是应变协调方程变为：第42页/共86页应力解法Chapter 6.3 最后提一下以应变分量ij作基本未知量的应变解法。由于应力与应变间的虎克定律是代数方程，应变解法的求解难度不会比应力解法有实质性的改善，而边界条件用应力表示则方便得多，所以很少采用应变解法。第43页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第44页/共86页应力函数解法Chapter 6.4 在位移解法中，引进三个单值连续的位移函数，使协调方程自动满足，问题被归结为求解三个用位移表示的平衡方程。

13、应变分量可由位移偏导数的组合来确定。与此类似，在应力解法中也可以引进某些自动满足平衡方程的函数，称为应力函数，把问题归结为求解用应力函数表示的协调方程。应力分量可由应力函数偏导数的组合来确定。第45页/共86页应力函数解法Chapter 6.4u 应力解法比较复杂u 与应力解法相比，此方法独立的应力分量只有3个事先满足几何方程，只需满足平衡方程：事先满足平衡方程，只需满足协调方程：几何方程几何方程本构方程本构方程满足平衡方程（满足平衡方程（3 3个）个）？方程？方程满足应力协调方程（满足应力协调方程（3 3个）个）满足平衡方程的满足平衡方程的第46页/共86页Bianchi恒等式：无体力平衡方

14、程：引进Beltrami应力函数张量：应力函数解法Chapter 6.4 可见两者形式上相似。只有3个独立分量定义应力分量：第47页/共86页应力函数解法Chapter 6.4 通过定义上述应力分量，可验证平衡方程将自动满足。上式便是用应力函数mn表示的应力公式。将上式代入应力协调方程，并考虑无体力情况得：这是应力函数解法的定解方程，称为应力函数协调方程。第48页/共86页应力函数解法Chapter 6.4 对于三维弹性力学问题，可选mn六个分量中的三个作为应力函数，共有17种选择方案，其中最常用的是：Maxwell应力函数：第49页/共86页应力函数解法Chapter 6.4对协调方程作如下

15、替换：于是可以写出用麦克斯韦应力函数表示的应力公式：第50页/共86页应力函数解法Chapter 6.4分量形式第51页/共86页应力函数解法Chapter 6.4 Morera应力函数：对协调方程作如下替换：第52页/共86页应力函数解法Chapter 6.4用莫雷拉(Morera)应力函数表示的应力公式：第53页/共86页应力函数解法Chapter 6.4分量形式第54页/共86页 则 就是平面问题中的艾瑞(Airy,G.B.)应力函数。若令莫雷拉应力函数为 则(x,y)就是柱形杆扭转问题中的普朗特应力函数应力函数解法Chapter 6.4 对于二维弹性力学问题，可进一步简化。若令麦克斯韦

16、应力函数为第55页/共86页应力函数解法Chapter 6.4 综上所述，应力函数解法既保留了应力解法的优点(能直接求解应力分量)，又吸收了位移解法的思想(能自动满足平衡方程，基本未知量降为3个)，所以是弹性理论中最常用的解法之一。第56页/共86页位移应力应变几何方程本构方程应力公式应力函数协调方程平衡方程红线：位移解法蓝线：应力函数解法Chapter 6.4 位移解法与应力函数解法的求解思路 圈表示物理量，框表示关系式，双框表示最后导出的定解方程，实箭头表示推导过程，虚箭头表示自动满足。第57页/共86页应力函数解法Chapter 6.4几何方几何方程程本构方本构方程程平衡方程平衡方程满满

17、足足应应力力公公式式积分积分（位移单值条件）（位移单值条件）协调方协调方程程满满足足解解出出解解出出自自动动满满足足自自动动满满足足第58页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第59页/共86页叠加原理Chapter 6.5 描述 作用在物体上的两组外力(包括表面力和体积力)的总和在物体内部所产生的效果(应力、应变及位移等)，等于此两组外力分别作用效果的总和。或者说物体受两组载荷共同作用时的应力或位移场就等于每组载荷单独作用时的应力或位移场之和，且与加载顺序无关。第60页/共86页叠

18、加原理Chapter 6.5设第一组载荷为体力 和面力 ，第二组为体力和面力 ，它们引起的应力和位移场分别为 和及 和 ，且仅考虑线弹性小变形情况，则联合载载引起的应力和位移场为具体阐述为：第61页/共86页叠加原理Chapter 6.5下面以应力解法为基础证明叠加原理的正确性：将叠加后的载荷代入应力解法各个方程后有：代入平衡方程 应力协调方程 力边界条件第62页/共86页叠加原理Chapter 6.5由于上述各式是线性微分方程或是代数方程，根据线性方程的性质将上述各式改写成：第63页/共86页由前提假设，和 分别是载荷 和单独作用时的解，根据它们满足的平衡方程，B-M方程和力边界条件可直接判

19、断上述各式的等号是成立的。因而证明了叠加后的应力场能满足应力解法的全部方程和边界条件，由此叠加原理得证。同理可以类似证明位移叠加原理的正确性。叠加原理Chapter 6.5第64页/共86页叠加原理Chapter 6.5小结 叠加原理用于位移边界条件时要求总位移满足给定的边界条件，而 和 单独不一定要满足位移边界条件。对于大变形情况，几何方程将出现非线性项，平衡方程也受到变形的影响，因为叠加原理不再适用。例如：同时受轴向和横向力的梁的纵横弯曲问题，薄壁构件的弹性稳定性问题，板壳结构的大挠度问题。第65页/共86页叠加原理Chapter 6.5 对于非线性弹性材料或弹塑性材料，本构方程是非线性的

20、；对于载荷随变形而变化的非保守力系情况或边界用非线性弹簧支撑的约束情况，边界条件是非线性的；叠加原理对这些情况都将失效。第66页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第67页/共86页解的唯一性定理Chapter 6.6u 无初应力的线弹性问题的解是唯一的u 给“试凑”解法提供了理论基础第68页/共86页解的唯一性定理Chapter 6.6下面给出定理的证明（反证法）：先假设存在两种不同的解，它们的位移场和应力场 及 都满足基本微分方程和给定边界条件。然后证明，对线性弹性问题两解之差必

21、为零，因而只能有唯一解。第69页/共86页解的唯一性定理Chapter 6.6 由前提假设，及 是同一物体在同一载荷 及相同边界条件下的两个解，它们都满足平衡方程 力边界条件 位移边界条件(在V内)(在S上)(在Su上)第70页/共86页解的唯一性定理Chapter 6.6 把以上左右两式对应相减，由叠加原理可知，两解之差 必然满足无体力平衡方程 和齐次边界条件第71页/共86页解的唯一性定理Chapter 6.6将(1)式两边乘ui，对体积V积分，并利用高斯公式得 其中第一项面积分的积分域为S=S+Su,根据(2)和(3)式，被积函数在边界上总有一个因子ui或ijj为零，所以第一项等于零。再

22、利用ij的对称性和线弹性应变能公式，上式可化为第72页/共86页解的唯一性定理Chapter 6.6 对于线弹性问题应变能处处正定，故上式要求W=0，即两解之差是ij0和ij0的无变形状态，因而 于是证明了应力场和应变场的解是唯一的。第73页/共86页解的唯一性定理Chapter 6.6小结一般说，位移场 和 之间还可能差一个刚体位移，但是绝大多数弹性力学问题都给定足以限制刚体运动的位移约束条件，因而位移场的解也是唯一的。以上证明的前提是叠加原理、应变能正定性和应力张量对称性。线弹性理论能自动满足这些条件，因为线弹性问题的解是唯一的。无论用什么方法求得的解，只要能满足全部基本方程和边界条件，就

23、一定是问题的真解。第74页/共86页解的唯一性定理Chapter 6.6 不满足唯一性定理的情况 (1)材料非线性 (2)几何非线性 (3)载荷与变形耦合问题（边界条件非线性）(4)有初应力的情况第75页/共86页微分提法、解法及一般原理 弹性力学问题的微分提法 位移解法 应力解法 应力函数解法 叠加原理 解的唯一性定理 圣维南原理Chapter 6第76页/共86页圣维南原理Chapter 6.7u 局部影响原理 由作用在物体局部表面上的自平衡力系（合力与合力矩为零）所引起的应力和应变，在远离作用区（距离远大于该局部作用区的线性尺寸）的地方将衰减到可以忽略不计的程度。第77页/共86页圣维南

24、原理Chapter 6.7u 静力等效原理 若把作用在物体局部表面上的外力，用另一组与它静力等效（合力与合力矩与它相等）的力系来代替，则这种等效处理对物体内部应力应变状态的影响将随远离该局部作用区的距离增加而迅速衰减。u 没有严格的数学证明，仅为事实所验证第78页/共86页圣维南原理Chapter 6.7 圣维南原理的意义u 可以估算一组自平衡力系的影响范围 r O(l)yxozrPyxozrPPyxozrPPPP第79页/共86页圣维南原理Chapter 6.7l第80页/共86页圣维南原理Chapter 6.7u 当无法严格给出边界上作用力的着点分布规则时，可以用一组静力等效力系代替u 在数学上和椭圆型方程的性质有关。弹性动力学问题（双曲型方程）不满足圣维南原理第81页/共86页圣维南原理Chapter 6.7u 利用圣维南原理可将位移边界转化为等效的力边界,如图所示:第82页/共86页圣维南原理Chapter 6.7 圣维南原理的适用条件u 只有当力的作用区为结构最小的特征尺度时，圣维南原理才成立。双力矩薄壁杆件：h b(R)l第83页/共86页圣维南原理Chapter 6.7薄壁杆件在工程结构中有广泛的应用机翼轨道舱第84页/共86页圣维南原理Chapter 6.7谢谢！第85页/共86页感谢您的观看！第86页/共86页