弹塑性力学3弹性力学解题方法

1、第3章 弹性力学解题方法 按位移求解弹性力学问题 按应力求解弹性力学问题 平面问题和应力函数 半逆解法广义胡克定律1678年，年，R. Hooke发表了固体受力后应力和应变关发表了固体受力后应力和应变关系的定律系的定律胡克定律。胡克定律。“有多大伸长，就有多大力有多大伸长，就有多大力”xExEEzy)(1zyxE)(1zxyyE)(1yxzzEGxyxyGyzyzGzxzx各向同性体的胡克定律还可以用应变表示应力。各向同性体的胡克定律还可以用应变表示应力。zzyyxxGGG222zxzxyzyzxyxyGGG)21)(1 (ELam 弹性常数弹性常数3-1 按位移求解弹性力学问题弹性力学的一般

2、问题中，共包含弹性力学的一般问题中，共包含15个未知函数，个未知函数，将用将用15方程来求解。方程来求解。对于各向同性的弹性体：对于各向同性的弹性体： 3个平衡微分方程；个平衡微分方程； 6个几何方程（微分方程）；个几何方程（微分方程）； 6个物理方程（广义胡克定律）。个物理方程（广义胡克定律）。 边界条件（与上述方程组成封闭的定解问题）边界条件（与上述方程组成封闭的定解问题）弹性力学问题解法的分类：弹性力学问题解法的分类： 取取位移位移作为基本未知量。作为基本未知量。 位移法位移法 取取应力应力作为基本未知量。作为基本未知量。 应力法应力法按位移求解弹性力学问题时按位移求解弹性力学问题时，取

3、取u,v,w作为基本未知量。作为基本未知量。几何方程几何方程物理方程物理方程消去应变消去应变zwGyvGxuGzyx222)()()(xwzuGzvywGyuxvGzxyzxy平衡方程平衡方程消去应力消去应力0)(0)(0)(222zyxfwGzGfvGyGfuGxG21)21)(1 (GGEGLam位移方程位移方程力的边界条件力的边界条件zzzyzxyyzyyxxzxyxxFnmlFnmlFnmlzwGyvGxuGzyx222)()()(xwzuGzvywGyuxvGzxyzxy消去消去应力应力zyxFzwGnywzvmGxwzulGFzvywnGyvGmxvyulGFzuxwnGyuxvm

4、GxuGl)2()()()()2()()()()2(按位移求解弹性力学问题按位移求解弹性力学问题优点：优点：未知函数的个数比较少，即仅有三个未知量未知函数的个数比较少，即仅有三个未知量u,v,w。缺点：缺点：必须求解三个联立的二阶偏微分方程。必须求解三个联立的二阶偏微分方程。按位移求解问题是普遍适用的方法，特别是在数值按位移求解问题是普遍适用的方法，特别是在数值解中得到了广泛的应用，例如在有限元法，差分法解中得到了广泛的应用，例如在有限元法，差分法等数值计算方法中，得到了很好的应用。等数值计算方法中，得到了很好的应用。? 例1设有半空间体，单位体积的质量为设有半空间体，单位体积的质量为，在水平

5、边，在水平边界上受均布压力界上受均布压力q的作用，试用位移法求各位移分的作用，试用位移法求各位移分量和应力分量，假设在量和应力分量，假设在z=h处处z方向的位移方向的位移w=0。qh解：解：由于载荷和弹性体由于载荷和弹性体对对z轴对称，并且为轴对称，并且为半空间体半空间体可以假设可以假设)(0, 0zwwvu体积应变体积应变dzdwzwyvxu222222222)(dzwdwzyxw代入代入Lam 位移方程位移方程0)2(22gdzwdGGgdzwd222BAzgzGw2)1 (421力的边界条件力的边界条件1, 0nmlqFFFzyx, 0qdzdwGdzdwz0)2(qGdzdwz)1 (

6、221)(0qGAgzGz)1 (221)1 (2210qGA)1 (221位移边界条件位移边界条件0)(hzw2)1 (421)1 (221ghGqhGB)(2)()1 (4210, 022zhqzhgGwvu位移分量位移分量0),()(21zxyzxyzyxgzqgzq应力分量应力分量3-2 按应力求解弹性力学问题按应力求解弹性力学问题时按应力求解弹性力学问题时，取取6个应力分量个应力分量作为作为基本未知量。基本未知量。 变形协调变形协调 方程方程物理方程物理方程消去应变消去应变)(1)(1)(12)(12)(12)(1222222222222222222222222222222222zx

7、yzyxyxyxzyxzxzxzyxzyzyzxzxzxzyyzyzyxyxxyxyyzxzzzxyzxyyyzxyzxxxzxzyzzyxyyx平衡方程平衡方程改变形式改变形式xfzfzxyfzfzyxfyfyxzfyfxfxfxzfyfxfxfxzfyfxfxfxzxzxzyyzyxxyzyxxxzyxxxzyxxx222222222222222111111)(1211)(1211)(1211zyx相容方程相容方程2222222zyx0)1 (0)1 (0)1 (0)1 (0)1 (0)1 (222222222222222zxzyyxxxxzxyzxyxxx体积力为零或为常量体积力为零或为

8、常量02推导参照教材推导参照教材应力第一不变量应力第一不变量是调和函数是调和函数左式两边分别作左式两边分别作Laplace运算运算000222222zyx000222222zxyzxy应力分量是双调和函数应力分量是双调和函数按应力求解弹性力学问题按应力求解弹性力学问题优点：优点：边界条件比较简单，并且得到的应力表达式边界条件比较简单，并且得到的应力表达式在大多数具体问题中比位移表达式简单在大多数具体问题中比位移表达式简单。缺点：缺点：未知函数较多，所要求解的二阶偏微分方程未知函数较多，所要求解的二阶偏微分方程比较复杂。比较复杂。按应力求解比按位移求解一般来说容易些。按应力求解比按位移求解一般来

9、说容易些。但就解决弹性体问题的普遍性而言，按位移求解更但就解决弹性体问题的普遍性而言，按位移求解更具有普遍性。具有普遍性。对于实际问题，适当的选择求解方法。对于实际问题，适当的选择求解方法。3-3 平面问题和应力函数v 平面问题平面问题的分类：平面问题的分类：平面应力问题平面应力问题平面应变问题平面应变问题平面应力问题平面应力问题构件几何形状特征：构件几何形状特征：薄板薄板外力外力:平行于中面，沿厚度均匀平行于中面，沿厚度均匀分布，表面不受外力作用。分布，表面不受外力作用。xyzyzo表面面力边界条件：表面面力边界条件：0, 00222hzyzhzxzhzz由于薄板厚度很小，应力由于薄板厚度很

10、小，应力分量均匀分布分量均匀分布),(),(),(0, 0, 0321yxFyxFyxFxyyxxzyzz中面中面平面应变问题平面应变问题构件几何形状特征：构件几何形状特征：具有很长纵向轴的柱体具有很长纵向轴的柱体纵向轴纵向轴横截面的大小和形状沿轴横截面的大小和形状沿轴线不变；外力与轴线垂直线不变；外力与轴线垂直并且沿轴线不变；主体两并且沿轴线不变；主体两端受端受固定约束固定约束。z方向上位移方向上位移0w位移发生在位移发生在oxy平面内平面内),(),(021yxvyxuw根据几何方程根据几何方程0, 0yzxzz根据物理方程根据物理方程0)(1yxzzE()zxy ),(),(),(0,

11、0),(321yxFyxFyxFxyyxxzyzyxzv 应力函数在平面问题中，引进应力函数的概念，往往使求解在平面问题中，引进应力函数的概念，往往使求解问题变得简单。问题变得简单。0yxyxx0yxyyx无体力存在时无体力存在时假定假定yxxyxyyx22222,平衡方程将自然满足平衡方程将自然满足应力函数Airyyx),(只需求解以应力函数表示的协调方程只需求解以应力函数表示的协调方程平面应力问题：平面应力问题：yxGxyGxyGxyxyx2222222221)()1 (21)()1 (21变形协调方程变形协调方程024422444yxxx022边界条件边界条件yxFxmyxlFyxmyl

12、222222平面问题归结为平面问题归结为求解满足双调和求解满足双调和方程和给定边界方程和给定边界条件的函数条件的函数),(yx? 例2图示很长的矩形柱体，材料的比重为图示很长的矩形柱体，材料的比重为，将其放入，将其放入形状相同的刚性槽内若不考虑摩擦力，设应力函数形状相同的刚性槽内若不考虑摩擦力，设应力函数的形式为的形式为2232DxCyByyAx试求各应力分量、应变分量以及位移分量。试求各应力分量、应变分量以及位移分量。haaxyo解：解：xAxxyxDAyxCByyxyyx2222622222根据根据Airy应力函数可得应力函数可得应力边界条件应力边界条件处，hy 00yxy刚性槽的条件刚性

13、槽的条件0aaxdx无关变量和xx0 x2A2hDyx10)1(1yxxE)(hyy)(1hyx61B21hC)(1)21)(1 (1)1(1hyEExyy0dxuxKhyyEdyvy)2(1)21)(1 (120时，0vy0K针对求解的问题，根据材料力学已知解或弹性体的针对求解的问题，根据材料力学已知解或弹性体的边界形状和受力情况，假设部分应力为某种形式的边界形状和受力情况，假设部分应力为某种形式的函数，从而推断出应力函数；函数，从而推断出应力函数；然后用方程和边界条件确定尚未求出的应力分量，然后用方程和边界条件确定尚未求出的应力分量，或完全确定原来假设的尚未全部定下来的应力。或完全确定原来

14、假设的尚未全部定下来的应力。如能满足弹性力学的全部条件，则这个解就是正确如能满足弹性力学的全部条件，则这个解就是正确的解答，如不能满足全部条件，则需另外假定，重的解答，如不能满足全部条件，则需另外假定，重新求解。新求解。由于根据已有解或经验作了一定假设，使得问题的由于根据已有解或经验作了一定假设，使得问题的求解过程得以大大简化。求解过程得以大大简化。3-4 半逆解法? 例3图示立柱（厚度为单位厚度），在其侧面上，作用图示立柱（厚度为单位厚度），在其侧面上，作用有均布剪力有均布剪力，试用半逆解法求其应力分布规律。，试用半逆解法求其应力分布规律。hxyo解：解：假定纵向纤维互不挤压假定纵向纤维互不

15、挤压022yx)()(),(21xfyxfyx代入代入0220)()(424414dxxfddxxfdy0)(, 0)(424414dxxfddxxfd上式对于上式对于y取任何值均应成立取任何值均应成立KJxExDxxfICxBxAxxf232231)()(对应力分量无影响对应力分量无影响2323)(ExDxCxBxAxyCBxAxEDxBAxyxyyx2326)26(02边界条件：边界条件：在在x=0处，处，0, 0 xyx0C在在x=h处，处，xyx, 0BhAh232（主要边界条件，需精确满足）（主要边界条件，需精确满足）在在y=0处，处，00hxydx0023hBxAxhxyo22hB

16、hA在在y=0处，处，00hxydx0023hBxAx00hydx0232 EhDh（次要边界条件，使用圣维南原理建立）（次要边界条件，使用圣维南原理建立）00hxxdx0323 EhDh0 EDxhxhyhxhxyyx)23()31 (20应力分量：应力分量：hxyo3dy)0( YX0 xydyx60yxy12h2hll图示梁对应的边界条件：图示梁对应的边界条件：:2hy0, 0 xyy:lx0,6xyxdydh3mindh3maxMM3dy 对应于矩形截面梁的纯弯曲问题应力分布。对应于矩形截面梁的纯弯曲问题应力分布。(1)(1)例例:矩形梁的纯弯曲矩形梁的纯弯曲xy12h2hlldh3m

17、indh3maxMM常数常数 d 与弯矩与弯矩 M 的关系：的关系：220hhxdy(1)0622hhdydy(2)Mdyyhhx222226hhMdydyMhd32)2(3hMd 或yIMxyhMx312yhMx)12/(3此结果与材力中结果相同此结果与材力中结果相同.位移分量的求解位移分量的求解xyl1hMM12/3hMyyIMx0 xy0y（1）应变分量）应变分量)(1xyyE)(1yxxEGxyxy(a)IMyEyIMyEx10 xy（2）位移分量）位移分量0 xvyuxyIMyExux1IMyEyvy(c)0 xvyuxyIMyExux1IMyEyvy(c)(222xfyEIMv)(

18、1yfxyEIMu(d)将式将式 (d) 代入代入 (c) 中第三式，中第三式，得：得：)()(12yfxfxEIM整理得：整理得：0)()(21xfyfxEIM)(2xfxEIM)(1yf(e)01)(uyyf022)(vxxEIMxf将上式代入式（将上式代入式（d），得），得0uyxyEIMu02222vxxEIMyEIMv(f)(f)0uyxyEIMu02222vxxEIMyEIMvxyl1hMM（1）两端简支）两端简支02222vxxEIMyEIMv0uyxyEIMu(f)000yxu000yxv0202vlEIMl00u00vEIMl2ylxEIMu)2( 22)(2yEIMxxlE

19、IMv(3-3)梁的挠曲线方程：梁的挠曲线方程：xxlEIMvy)(20 与材力中结果相同与材力中结果相同00ylxv（2）悬臂梁）悬臂梁02222vxxEIMyEIMv0uyxyEIMu(f)边界条件边界条件0lxv0lxu22hyhh/2h/2由式（由式（f）可知，此边界条件无法满足。）可知，此边界条件无法满足。边界条件改写为：边界条件改写为：0, 000ylxylxvu（中点不动）（中点不动）00ylxxv（轴线在端部不转动）（轴线在端部不转动）00u0202vllEIM0lEIM得：得：00uEIMlv220EIMlyxlEIMu)( 222)(2yEIMxlEIMv(3-4)h/2h

20、/2挠曲线方程：挠曲线方程：20)(2|xlEIMvy与材料力学中结果相同与材料力学中结果相同讨论：讨论： 若取固定端边界条件为：若取固定端边界条件为：h/2h/20, 000ylxylxvu（中点不动）（中点不动）00ylxyu(中点处竖向线段转角为零中点处竖向线段转角为零)00u得到：得到：0202vlEIMl0EIMl02222vxxEIMyEIMv0uyxyEIMu求得：求得：00uEIMlv220EIMl此结果与前面情形相同。此结果与前面情形相同。（1）024422444yyxx(2-27)（2）xyyx,然后将然后将 代入式（代入式（2-26）求出应力分量：）求出应力分量：),(yx先由方程（先由方程（2-27）求出应