平面问题有限元解法公式推导讲解PPT课件.ppt

《平面问题有限元解法公式推导讲解PPT课件.ppt》由会员分享，可在线阅读，更多相关《平面问题有限元解法公式推导讲解PPT课件.ppt（78页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于平面问题有限元解法公式推导讲解第一张，PPT共七十八页，创作于2022年6月有限单元法基本思想n有限单元法的思想是将物体（连续的求解域）离散成有限个且按一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，从而将一个连续的无限连续的无限自由度自由度问题简化为离散的有限自由度问题有限自由度问题求解的一种数值分析法数值分析法。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。n有限单元法的分析步骤如下：q物体离散化q单元特性分析q单元组集，整体分析q求解未知节点的位移q由节点的位移求解各单元的位移和应力第二张，PPT共七十八页，创作于2022年6月06.04.2023

2、有限元单元模型中几个重要概念n单元q网格划分中每一个小的块体n节点q确定单元形状、单元之间相互联结的点n节点力q单元上节点处的结构内力n载荷q作用在单元节点上的外力（集中力、分布力）n约束p限制某些节点的某些自由度n弹性模量（杨式模量）En泊松比（横向变形系数）n密度单元单元单元单元载荷节点节点力约束第三张，PPT共七十八页，创作于2022年6月06.04.20231.研究内容研究内容内容：内容：弹性体在外力或温度作用下的应力、变形、位弹性体在外力或温度作用下的应力、变形、位移等分布规律。移等分布规律。任务：任务：解决弹性体的强度、刚度、稳定性问题。解决弹性体的强度、刚度、稳定性问题。弹性力学

3、的内容及基本假定弹性力学的内容及基本假定2.研究对象研究对象一般弹性实体结构：一般弹性实体结构：三维弹性固体、板状结构、杆件等三维弹性固体、板状结构、杆件等第四张，PPT共七十八页，创作于2022年6月06.04.2023弹性力学的内容及基本假定弹性力学的内容及基本假定3.研究方法研究方法由平衡方程、几何方程、物理方程三方面分析由平衡方程、几何方程、物理方程三方面分析4.数学理论基础数学理论基础 偏微分方程（高阶，二、三个变量）偏微分方程（高阶，二、三个变量）数值解法数值解法：能量法（变分法）、差分法、：能量法（变分法）、差分法、有限单元法等。有限单元法等。第五张，PPT共七十八页，创作于20

4、22年6月06.04.2023弹性力学的内容及基本假定弹性力学的内容及基本假定5.基本假定基本假定（1）.连续性假定连续性假定整个物体的体积都被组成物体的介质充满，不留下任何空隙。整个物体的体积都被组成物体的介质充满，不留下任何空隙。作用：使得使得、u 等量表示成坐标的连续函数。等量表示成坐标的连续函数。第六张，PPT共七十八页，创作于2022年6月06.04.2023弹性力学的内容及基本假定弹性力学的内容及基本假定（2）.完全弹性假定完全弹性假定 假定物体完全服从虎克（假定物体完全服从虎克（HookeHooke）定律，应力与应变间成线性比）定律，应力与应变间成线性比例关系。例关系。脆性材料脆

5、性材料 一直到破坏前，都可近似为线弹性的；一直到破坏前，都可近似为线弹性的；塑性材料塑性材料 比例阶段，可视为线弹性的。比例阶段，可视为线弹性的。（3）.均匀性假定均匀性假定 假定整个物体是由同一种材料组成的，各部分材料性质相同。假定整个物体是由同一种材料组成的，各部分材料性质相同。作用：作用：弹性常数（弹性常数（E、）等）等不随位置坐标而变化；不随位置坐标而变化；取微元体分析的结果可应用于整个物体。取微元体分析的结果可应用于整个物体。第七张，PPT共七十八页，创作于2022年6月06.04.2023弹性力学的内容及基本假定弹性力学的内容及基本假定（4）.各向同性假定各向同性假定（5）.小变形

6、假定小变形假定 假定物体内一点的力学性质在所有各个方向都相同。假定物体内一点的力学性质在所有各个方向都相同。作用：作用：弹性常数（弹性常数（E、）不随坐标方向而变化；不随坐标方向而变化；假定位移和形变是微小的，即物体受力后物体内各点位移假定位移和形变是微小的，即物体受力后物体内各点位移远远小于物体的原来的尺寸。远远小于物体的原来的尺寸。作用：作用：建立方程时，可略去高阶微量；建立方程时，可略去高阶微量；可用变形前的尺寸代替变形后的尺寸。可用变形前的尺寸代替变形后的尺寸。使求解的方程使求解的方程线性化线性化。第八张，PPT共七十八页，创作于2022年6月06.04.2023基本概念：基本概念：外

7、力、应力、形变、位移。外力、应力、形变、位移。1.外力：外力：体力、面力体力、面力(1)体力体力 分布在物体分布在物体体积体积内的力内的力 体力分布集度体力分布集度（矢量）（矢量）xyzO单位：单位：N/m3kN/m3说明：说明：f 是坐标的连续分布函数是坐标的连续分布函数；弹性力学中的几个基本概念弹性力学中的几个基本概念p第九张，PPT共七十八页，创作于2022年6月06.04.2023(2)面力面力 分布在物体表面的力分布在物体表面的力 面力分布集度（矢量）面力分布集度（矢量）xyzO单位：单位：1N/m2=1Pa(帕)1MN/m2=106Pa=1MPa(兆帕)说明：说明：弹性力学中的几个

8、基本概念弹性力学中的几个基本概念是坐标的连续分布函数是坐标的连续分布函数；p第十张，PPT共七十八页，创作于2022年6月06.04.20232.应力应力(1)一点应力的概念一点应力的概念AF内力内力(1)物体内部分子或原子间的相互作用力物体内部分子或原子间的相互作用力;(2)由于外力作用引起的相互作用力由于外力作用引起的相互作用力.(不考虑不考虑)P截面上截面上P点的应力点的应力应力矢量应力矢量.的极限方向的极限方向应力分量应力分量n(法线法线)应力的法向分量应力的法向分量 正应力正应力应力的切向分量应力的切向分量 切应力切应力单位单位:MPa(兆帕)应力关于坐标连续分布应力关于坐标连续分布

9、弹性力学中的几个基本概念弹性力学中的几个基本概念第十一张，PPT共七十八页，创作于2022年6月06.04.2023(2)一点的应力状态一点的应力状态通过一点通过一点P 的各个面上应力状况的集合的各个面上应力状况的集合 称为一点的应力状态称为一点的应力状态x面的应力：面的应力：y面的应力：面的应力：z面的应力：面的应力：弹性力学中的几个基本概念弹性力学中的几个基本概念第十二张，PPT共七十八页，创作于2022年6月06.04.2023用矩阵表示：用矩阵表示：其中，只有其中，只有6个量独立。个量独立。切应力互等定理切应力互等定理应力应力正负号正负号的规定：的规定：正应力正应力 拉为正，压为负。拉

10、为正，压为负。切应力切应力 坐标坐标正面正面上，与坐标正向一致时为正；上，与坐标正向一致时为正；坐标坐标负面负面上，与坐标正向相反时为正。上，与坐标正向相反时为正。xyzO弹性力学中的几个基本概念弹性力学中的几个基本概念第十三张，PPT共七十八页，创作于2022年6月06.04.20233.形变形变形变形变 物体的形状改变物体的形状改变xyzO（1）线段长度的改变）线段长度的改变（2）两线段间夹角的改变。）两线段间夹角的改变。PBCA用正应变用正应变度量度量切应变切应变度量度量（切应变（切应变两垂直线段夹角两垂直线段夹角（直角）（直角）的改变量）的改变量）三个方向的正应变：三个方向的正应变：三

11、个平面内的切应变：三个平面内的切应变：(1)一点形变的度量一点形变的度量应变的正负：应变的正负：正应变正应变：伸长时为正，缩短时为负；伸长时为正，缩短时为负；切应变切应变：以直角变小时为正，变大时为负；以直角变小时为正，变大时为负；弹性力学中的几个基本概念弹性力学中的几个基本概念第十四张，PPT共七十八页，创作于2022年6月06.04.2023(2)一点应变状态一点应变状态其中其中应变无量纲；应变无量纲；4.位移位移注：注：一点的位移一点的位移 矢量矢量S应变分量均为位置坐标的函数应变分量均为位置坐标的函数xyzOSwuvP位移分量：位移分量：u x方向的位移方向的位移 分量；分量；v y方

12、向的位移方向的位移 分量；分量；w z方向的位移方向的位移 分量。分量。量纲：量纲：m 或 mm弹性力学中的几个基本概念弹性力学中的几个基本概念第十五张，PPT共七十八页，创作于2022年6月06.04.2023 工程力学问题建立力学模型的过程中，一般从三方面进行简化：结构简化 如空间问题向平面问题的简化，向轴对称问题的简化，实体结构向板、壳结构的简化。受力简化 如：根据圣维南原理，复杂力系简化为等效力系等。材料简化根据各向同性、连续、均匀等假设进行简化。第十六张，PPT共七十八页，创作于2022年6月06.04.2023平面问题的基本理论平面问题的基本理论n任何一个实际的弹性力学问题都是空间

13、问题，但是如任何一个实际的弹性力学问题都是空间问题，但是如果所考察的弹性体具有某种果所考察的弹性体具有某种特殊的形状，特殊的形状，并且承受的并且承受的是某些是某些特殊的外力和约束特殊的外力和约束，就可以把空间问题简化为，就可以把空间问题简化为近似的近似的平面问题平面问题。n两种典型的平面问题两种典型的平面问题q平面应力问题平面应力问题q平面应变问题平面应变问题第十七张，PPT共七十八页，创作于2022年6月06.04.2023平面应力问题平面应力问题(1)几何特征几何特征xyyztba 一个方向的尺寸比另两个方一个方向的尺寸比另两个方向的尺寸小得多。向的尺寸小得多。平板平板如：板式吊钩，旋转圆

14、盘，工字形梁的腹板等如：板式吊钩，旋转圆盘，工字形梁的腹板等(2)受力特征受力特征外力外力（体力、面力）和（体力、面力）和约束约束，仅，仅平行于板面作用平行于板面作用，沿，沿 z 方方向不变化。向不变化。第十八张，PPT共七十八页，创作于2022年6月06.04.2023xyyztba(3)应力特征应力特征如图选取坐标系，以板的中面为如图选取坐标系，以板的中面为xy 平面，垂直于中面的任一直线为平面，垂直于中面的任一直线为 z 轴。轴。由于板面上不受力，有由于板面上不受力，有因板很薄，且外力沿因板很薄，且外力沿 z 轴方向不变。轴方向不变。可认为可认为整个薄板的各整个薄板的各点点都有：都有：由

15、切应力互等定理，有由切应力互等定理，有结论：结论：平面应力问题只有三个应力分量：平面应力问题只有三个应力分量：xy应变分量、位移分量也仅为应变分量、位移分量也仅为 x、y 的函数，与的函数，与 z 无关。无关。第十九张，PPT共七十八页，创作于2022年6月06.04.2023平面应变问题平面应变问题(1)几何特征几何特征水坝水坝滚柱滚柱厚壁圆筒厚壁圆筒 一个方向的尺寸比另两个一个方向的尺寸比另两个方向的尺寸方向的尺寸大得多，且沿长度大得多，且沿长度方向几何形状和尺寸不变化。方向几何形状和尺寸不变化。近似认为无限长近似认为无限长(2)外力特征外力特征 外力外力（体力、面力）（体力、面力）平行于

16、横截面平行于横截面作用，且作用，且沿长度沿长度 z 方向不变化方向不变化。约束约束 沿长度沿长度 z 方向不变化方向不变化。(3)变形特征变形特征 如图建立坐标系：以任一横截面为如图建立坐标系：以任一横截面为 xy 面，任一纵线为面，任一纵线为 z 轴。轴。设设 z方向为无限长，则方向为无限长，则沿沿 z 方向都不变化，方向都不变化，仅为仅为 x，y 的函数。的函数。任一横截面均可视为对称面任一横截面均可视为对称面第二十张，PPT共七十八页，创作于2022年6月06.04.2023水坝水坝任一横截面均可视为对称面，则有任一横截面均可视为对称面，则有所有各点的位移矢量都平行于所有各点的位移矢量都

17、平行于 x y 平面。平面。平面位移问题平面位移问题 平面应变问题平面应变问题注：注：平面应变问题中平面应变问题中但是，但是，第二十一张，PPT共七十八页，创作于2022年6月06.04.2023 如图所示三种情形，是否都属平面问题？是平面应如图所示三种情形，是否都属平面问题？是平面应力问题还是平面应变问题？力问题还是平面应变问题？平面应力问题平面应力问题平面应变问题平面应变问题非平面问题非平面问题第二十二张，PPT共七十八页，创作于2022年6月06.04.2023三大基本方程三大基本方程n根据根据静力学静力学、几何学几何学和和物理学物理学三方面条件，建立三套方程。三方面条件，建立三套方程。

18、q平面问题中，根据微分体的平衡条件，建立平面问题中，根据微分体的平衡条件，建立平衡微分方程平衡微分方程：(1-1)q根据微分线段上形变与位移之间的几何关系，建立根据微分线段上形变与位移之间的几何关系，建立几何方程几何方程：(1-2)q根据应力与形变之间的物理关系，建立根据应力与形变之间的物理关系，建立物理方程物理方程：(1-3)(1-3)第二十三张，PPT共七十八页，创作于2022年6月06.04.2023平衡微分方程平衡微分方程n从弹性体中取出一个微分体，根据平衡条件导从弹性体中取出一个微分体，根据平衡条件导出应力分量与体力分量之间的关系式，也就是出应力分量与体力分量之间的关系式，也就是平面

19、问题的平面问题的平衡微分方程平衡微分方程。n从弹性体中取出一个微小的正平行六面体，从弹性体中取出一个微小的正平行六面体，它在它在x和和y方向的尺寸分别为方向的尺寸分别为dx和和dy，在，在z方向的尺方向的尺寸为一个单位长度。寸为一个单位长度。n以以x为投影轴，列出投影的平衡方程：为投影轴，列出投影的平衡方程：n约简以后，两边除以约简以后，两边除以dxdy，得：，得：n同理，以同理，以y为投影轴，列出投影的平衡方程，化简得为投影轴，列出投影的平衡方程，化简得：第二十四张，PPT共七十八页，创作于2022年6月06.04.2023几何方程几何方程n经过弹性体内的任意一点经过弹性体内的任意一点P，沿

20、，沿x轴和轴和y轴轴的正方向取两个微小长度的线段的正方向取两个微小长度的线段PAdx和和PBdy。假定弹性体受力后，。假定弹性体受力后，P,A,B三点分别移动到三点分别移动到P,A,B.n线段线段PA的正应变是：的正应变是：注注:由于位移微小，由于位移微小，y方向的位移方向的位移v引起的引起的PA的伸缩，是高一阶微量，略去不计。的伸缩，是高一阶微量，略去不计。n线段线段PB的正应变是：的正应变是：n线段线段PA与与 PB之间的直角的改变，即切应变之间的直角的改变，即切应变n线段线段PA的转角的转角是是：n线段线段PB的转角的转角是：是：第二十五张，PPT共七十八页，创作于2022年6月06.0

21、4.2023物理方程物理方程n在理想的弹性体中，形变分量和应力分量之间的关系，在材料力在理想的弹性体中，形变分量和应力分量之间的关系，在材料力学根据胡克定律导出如下：学根据胡克定律导出如下：n在平面应力问题中，在平面应力问题中，式变为：式变为：n在平面应变问题中，在平面应变问题中，只要将上式中的只要将上式中的E换为换为 ，换为换为 就得到平面应变问题的物理方程。就得到平面应变问题的物理方程。第二十六张，PPT共七十八页，创作于2022年6月06.04.2023n假定已知任一点假定已知任一点P处坐标面上的应力分量处坐标面上的应力分量x，y，x y=y x。求经过该点的，平行于。求经过该点的，平行

22、于z轴而倾斜于轴而倾斜于x轴和轴和 y轴的轴的任何倾斜面上应力。任何倾斜面上应力。n在在P点附近取一个平面点附近取一个平面AB，它平行于上述斜面，并经过，它平行于上述斜面，并经过P点划出一个微小的三棱柱点划出一个微小的三棱柱PAB。当。当AB无限小而趋于无限小而趋于P点时，平面点时，平面AB上的应力就成为斜面上的应力。上的应力就成为斜面上的应力。平面问题中一点的应力状态平面问题中一点的应力状态n设斜面设斜面AB 的长度为的长度为ds，则，则PB面及面及PA面的长度分别为面的长度分别为 lds及及mds，而，而PAB的面积为的面积为 ldsmds/2，棱柱的厚度设为，棱柱的厚度设为1。n由由x轴

23、平衡条件，得：轴平衡条件，得：n其中，其中，fx为体力分量。将上式除以为体力分量。将上式除以ds，并令，并令ds趋于趋于0（斜面（斜面AB趋于趋于P点），点），即得：即得：n由由y轴平衡条件，得：轴平衡条件，得：n用用n表示斜面表示斜面AB的的外法线方向，其方向余弦为：外法线方向，其方向余弦为：第二十七张，PPT共七十八页，创作于2022年6月06.04.2023边界条件边界条件q若在若在su部分部分边界上给定了约束位移分量边界上给定了约束位移分量 和和 ，则对于此边界上的，则对于此边界上的每一点，位移函数每一点，位移函数u和和v应满足条件：应满足条件：q其中其中（u）s 和和（v）s 是位移

24、的边界值，是位移的边界值，和和 在边界上是坐标的已在边界上是坐标的已知函数。知函数。n边界条件表示在边界上位移与约束，或应力与面力之间的关系式。它可以分为边界条件表示在边界上位移与约束，或应力与面力之间的关系式。它可以分为位移位移边界条件边界条件、应力边界条件应力边界条件和和混合边界条件混合边界条件。q位移边界条件位移边界条件：q应力边界条件应力边界条件：q若在若在su部分部分边界上给定了面力边界上给定了面力 和和 ，则由平衡条件得出平面应，则由平衡条件得出平面应力问题的应力力问题的应力(或面力）边界条件为：或面力）边界条件为：其中，其中，l,m是边界面外法线的方向余弦。是边界面外法线的方向余

25、弦。第二十八张，PPT共七十八页，创作于2022年6月06.04.2023圣维南原理圣维南原理n在求解弹性力学问题时，应力分量、形变分量和位移分量必须在求解弹性力学问题时，应力分量、形变分量和位移分量必须满足区域内的满足区域内的三套基本方程三套基本方程，还必须满足边界上的，还必须满足边界上的边界条件边界条件。但。但是，要使边界条件得到完全满足，往往遇到很大的困难。是，要使边界条件得到完全满足，往往遇到很大的困难。n圣维南原理圣维南原理可为简化局部边界上的应力边界条件提供很大方便。可为简化局部边界上的应力边界条件提供很大方便。n圣维南原理圣维南原理表明，如果把物体的一小部分边界上的面力，变换为分

26、表明，如果把物体的一小部分边界上的面力，变换为分布不同但布不同但静力等效静力等效的面力（的面力（主矢相同，对同一点的主矩也相同主矢相同，对同一点的主矩也相同），），那么，近处的应力分布将有显著的改变，但是远处所受的影响那么，近处的应力分布将有显著的改变，但是远处所受的影响可以不计。可以不计。第二十九张，PPT共七十八页，创作于2022年6月06.04.2023圣维南原理的应用圣维南原理的应用n例，设有柱形构件，在两端截面的形心例，设有柱形构件，在两端截面的形心受到大小相等而方向相反的拉力受到大小相等而方向相反的拉力F（a）。如果把一端或两端的拉力变换为静力等效如果把一端或两端的拉力变换为静力等

27、效的力，则只有虚线划出的部分的应力分布的力，则只有虚线划出的部分的应力分布有显著的改变，而其余部分所受影响是可有显著的改变，而其余部分所受影响是可以不计的。以不计的。n由于（由于（d）图中，面力连续分布，边界条件简单，应力容易求得。其它三种情况，应力难以）图中，面力连续分布，边界条件简单，应力容易求得。其它三种情况，应力难以求得。把求得。把d情况下的应力解答应用到其它三个情况，虽不能满足两端的应力边界条件，情况下的应力解答应用到其它三个情况，虽不能满足两端的应力边界条件，但仍然可以表明离杆端较远处的应力状态，没有显著的误差。但仍然可以表明离杆端较远处的应力状态，没有显著的误差。n图图e，构件右

28、端有位移边界条件，构件右端有位移边界条件，d情况的解答，不能满足位移情况的解答，不能满足位移边界条件，但边界条件，但e图右端的面力，一定是合成为经过截面形心的力图右端的面力，一定是合成为经过截面形心的力F。所以把图。所以把图d情况的解情况的解答应用于图答应用于图e时，仍然只是在靠近两端处有显著的误差，而在离两端较远之处，误差可时，仍然只是在靠近两端处有显著的误差，而在离两端较远之处，误差可以不计。以不计。第三十张，PPT共七十八页，创作于2022年6月06.04.2023圣维南原理的应用圣维南原理的应用n例，厚度例，厚度=1的梁中，左右两端的梁中，左右两端x=l，的边界面是正、负，的边界面是正

29、、负x面，其上作用有一般分面，其上作用有一般分布的面力布的面力 。按照严格的应力边界条件，应力分量在边界上满足：。按照严格的应力边界条件，应力分量在边界上满足：n上式要求在边界上上式要求在边界上y值不同的各点，应力分量与对应的面力分量必须处处相等，这种严值不同的各点，应力分量与对应的面力分量必须处处相等，这种严格的条件是较难满足的。格的条件是较难满足的。n当当lh时，时，x=l 是梁的边界的一小部分，可以应用是梁的边界的一小部分，可以应用圣维南原理，利用静力等效条件圣维南原理，利用静力等效条件来代替，即，使来代替，即，使应力的主矢量应力的主矢量和和主矩主矩分别等于分别等于对应的对应的面力的主矢

30、量面力的主矢量和和主矩主矩。第三十一张，PPT共七十八页，创作于2022年6月06.04.2023圣维南原理的应用圣维南原理的应用n应力的主矢量和主矩的绝对值分别等于面力的主矢量和主矩的绝对值；应力的主矢量和主矩的绝对值分别等于面力的主矢量和主矩的绝对值；n面力的主矢量和主矩的方向就是应力的主矢量和主矩的方向。面力的主矢量和主矩的方向就是应力的主矢量和主矩的方向。第三十二张，PPT共七十八页，创作于2022年6月06.04.2023按位移法求解平面问题按位移法求解平面问题n以上几节已经建立了弹性力学平面问题的以上几节已经建立了弹性力学平面问题的 基本方程和边界条件，基本方程和边界条件，即：即：

31、平衡微分方程、几何方程平衡微分方程、几何方程和和物理方程物理方程，以及，以及位移的边界条件位移的边界条件和和应力的边界条件应力的边界条件。n求解弹性力学平面问题即求解求解弹性力学平面问题即求解3个应力分量个应力分量、3个形变分量个形变分量及及2个个位移分量位移分量的未知函数。通常采用类似于代数方程中消元法进行求的未知函数。通常采用类似于代数方程中消元法进行求解。解。n按位移求解的方法，称为按位移求解的方法，称为位移法位移法。它以位移分量为基本未知函数。它以位移分量为基本未知函数。n按应力求解的方法，称为按应力求解的方法，称为应力法应力法。它以应力分量为基本未知函数。它以应力分量为基本未知函数。

32、第三十三张，PPT共七十八页，创作于2022年6月06.04.2023按位移法求解平面问题按位移法求解平面问题n平面问题中，取位移分量平面问题中，取位移分量u和和v为基本未知函数。为基本未知函数。n从方程中消去形变分量和应力分量：从方程中消去形变分量和应力分量：n将几何方程代入上式将几何方程代入上式n利用平衡微分方程和边界条件，导出用位移表示的平衡微分方程：利用平衡微分方程和边界条件，导出用位移表示的平衡微分方程：第三十四张，PPT共七十八页，创作于2022年6月06.04.2023按位移法求解平面问题按位移法求解平面问题n利用应力边界条件利用应力边界条件n得到用位移表示的应力边界条件得到用位

33、移表示的应力边界条件其中：其中：n位移边界条件如（位移边界条件如（1-4）不变）不变n按位移法求解平面应力问题时，要使位移分量在区域内满足平衡微分方程，在按位移法求解平面应力问题时，要使位移分量在区域内满足平衡微分方程，在边界上满足位移边界条件或应力边界条件。边界上满足位移边界条件或应力边界条件。第三十五张，PPT共七十八页，创作于2022年6月06.04.2023按位移法求解平面问题（例题）按位移法求解平面问题（例题）n设有如图所示的杆件，在设有如图所示的杆件，在y方向的上端为固定，而方向的上端为固定，而下端为自由，受自重体力下端为自由，受自重体力fx0，fyg的作用。试的作用。试用位移法求

34、解此问题。用位移法求解此问题。n解：将这个问题简化为一维问题处理。解：将这个问题简化为一维问题处理。n设设u=0,v=v(y)，泊松比，泊松比0。代入位移表示的平衡微分方程，。代入位移表示的平衡微分方程，得：得：n第一式自然满足，第二式成为：第一式自然满足，第二式成为：n解出：解出：第三十六张，PPT共七十八页，创作于2022年6月06.04.2023按位移法求解平面问题（例题）按位移法求解平面问题（例题）n设有左图所示的杆件，在设有左图所示的杆件，在y方方向的上端为固定，而下端为向的上端为固定，而下端为自由，受自重体力自由，受自重体力fx0，fyg的作用。试用位移法求的作用。试用位移法求解此

35、问题。解此问题。n解出：解出：n上下边的边界条件分别要求：上下边的边界条件分别要求：n将（将（a）式代入（）式代入（b）式得：）式得：B0，n再代入（再代入（c）式，即得：）式，即得：n得到解答：得到解答：第三十七张，PPT共七十八页，创作于2022年6月06.04.2023有限元单元法分析步骤（一）n结构离散化结构离散化q 将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点连接。结构的将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点连接。结构的离散化是有限元法分析的第一步，关系到计算精度和效率，包括以下三个方面：离散化是有限元法分析的第一步，关系到计算精度和效率，包

36、括以下三个方面：n单元类型的选择单元类型的选择。选定单元类型，确定单元形状、单元节点数、。选定单元类型，确定单元形状、单元节点数、节点自由度数等。节点自由度数等。n单元划分单元划分。网格划分越细，节点越多，计算结果越精确，但计算。网格划分越细，节点越多，计算结果越精确，但计算量越大。网格加密到一定程度后计算精度提高就不明显，对应应量越大。网格加密到一定程度后计算精度提高就不明显，对应应力变化平缓区域不必要细分网格。力变化平缓区域不必要细分网格。n节点编码节点编码。注意：注意：有限元分析的结构有限元分析的结构已不是已不是原有的物体或结构物原有的物体或结构物，而是由同样材料、众多，而是由同样材料、

37、众多单元以一定方式连接成的离散物体。用有限元分析计算所获得的结果是单元以一定方式连接成的离散物体。用有限元分析计算所获得的结果是近似的（满足工程要求即可）。近似的（满足工程要求即可）。平面问题有限单元法基本概念平面问题有限单元法基本概念第三十八张，PPT共七十八页，创作于2022年6月06.04.2023有限元单元法分析步骤（二）有限元单元法分析步骤（二）n单元特性分析单元特性分析q 选择未知量模式选择未知量模式n选择节点位移作为基本未知量时，称为位移法；选择节点位移作为基本未知量时，称为位移法；n选节点力作为基本未知量时，称为力法；选节点力作为基本未知量时，称为力法；n取一部分节点位移和一部

38、分节点力作为未知量，称为混合法。取一部分节点位移和一部分节点力作为未知量，称为混合法。q分析单元力学性质分析单元力学性质n根据单元材料性质、形状、尺寸、节点数目、位置等，找出单元根据单元材料性质、形状、尺寸、节点数目、位置等，找出单元节点力和节点位移关系式，应用几何方程和物理方程建立力和位节点力和节点位移关系式，应用几何方程和物理方程建立力和位移的方程式，从而导出单元刚度矩阵。移的方程式，从而导出单元刚度矩阵。q计算等效节点力计算等效节点力n作用在单元边界上的表面力、体积力或集中力都需要等效地移到作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，即用等效力来替代所有作用在单元上

39、的力。节点上去，即用等效力来替代所有作用在单元上的力。第三十九张，PPT共七十八页，创作于2022年6月06.04.2023有限元单元法分析步骤（三）n整体分析整体分析q集成整体节点载荷矢量集成整体节点载荷矢量 F。结构离散化后，单元之间通过节点传递。结构离散化后，单元之间通过节点传递力，作用在单元边界上的表面力、体积力或集中力都需要等效地移到力，作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体节点编码顺节点上去，形成等效节点载荷。将所有节点载荷按照整体节点编码顺序组集成整体节点载荷矢量。序组集成整体节点载荷矢量。q组成整体刚度矩阵组成

40、整体刚度矩阵K，得到总体平衡方程：，得到总体平衡方程：q引进边界约束条件，解总体平衡方程求出节点位移。引进边界约束条件，解总体平衡方程求出节点位移。通过上述分析可以看出有限单元法的基本思想是通过上述分析可以看出有限单元法的基本思想是“一分一合一分一合”，分是为了进行单，分是为了进行单元分析，合是为了对整体的结构进行综合分析。元分析，合是为了对整体的结构进行综合分析。第四十张，PPT共七十八页，创作于2022年6月06.04.2023有限单元法中基本量的矩阵表示有限单元法中基本量的矩阵表示n有限单元法有限单元法(FEM)中，为了简洁清晰地表示各个基本量以及它们之中，为了简洁清晰地表示各个基本量以

41、及它们之间的关系，也为了便于编制程序利用计算机进行计算，广泛采用矩间的关系，也为了便于编制程序利用计算机进行计算，广泛采用矩阵表示和矩阵运算。阵表示和矩阵运算。q平面问题中平面问题中，物体受体力，可用，物体受体力，可用体力列阵体力列阵表示：表示：(1)q物体受面力，可用物体受面力，可用面力列阵面力列阵表示：表示：(2)q3个应力分量的个应力分量的应力列阵应力列阵表示：表示：(3)q3个形变分量的个形变分量的应变列阵应变列阵表示：表示：(4)q2个位移分量的个位移分量的位移列阵位移列阵表示：表示：(5)第四十一张，PPT共七十八页，创作于2022年6月06.04.2023弹性力学中基本方程的矩阵

42、表示弹性力学中基本方程的矩阵表示q几何方程的几何方程的矩阵表示为：矩阵表示为：(6)q物理方程物理方程矩阵表示为：矩阵表示为：(7)q 利用应力列阵和应变列阵（利用应力列阵和应变列阵（3）、（）、（4）得：）得：(8)其中矩阵其中矩阵 (9)只与弹性常数只与弹性常数E及及有关，称为有关，称为平面问题的弹性矩阵平面问题的弹性矩阵。第四十二张，PPT共七十八页，创作于2022年6月06.04.2023虚位移原理虚位移原理n用用u*和和v*表示虚位移，用表示虚位移，用 表示与该虚位移相应的虚应变。表示与该虚位移相应的虚应变。n根据根据虚功方程虚功方程：处于平衡状态的变形体，：处于平衡状态的变形体，外

43、力外力在在虚位移虚位移上所做的上所做的虚功虚功等于等于应力应力在在虚应变虚应变上所做的上所做的虚功虚功。n对于厚度为对于厚度为t的薄板，虚功方程可用矩阵表示为：的薄板，虚功方程可用矩阵表示为：n其中，其中，分别为分别为体力列阵体力列阵，面力列阵面力列阵和和应力列阵应力列阵。为为虚位移列阵虚位移列阵n有限单元法中，作用于弹性体的各种外力常以作用于某些点的等效集中力来有限单元法中，作用于弹性体的各种外力常以作用于某些点的等效集中力来代替。在厚度为代替。在厚度为t的薄板上，设作用于的薄板上，设作用于i点的集中力沿点的集中力沿x及及y方向的分量为方向的分量为Fix,Fiy,作用作用于于j点的力为点的力

44、为Fjx,Fjy等。这些集中力以及它们相应的虚位移用列阵表示为：等。这些集中力以及它们相应的虚位移用列阵表示为：为为虚应变列阵虚应变列阵第四十三张，PPT共七十八页，创作于2022年6月06.04.2023虚位移原理虚位移原理(续续)n代入虚功方程，得：代入虚功方程，得：n 上式为上式为集中力作用下的虚功方程集中力作用下的虚功方程。n集中力列阵集中力列阵 （13）n虚位移列阵虚位移列阵 （14）n外力在虚位移上所做的功为：外力在虚位移上所做的功为：第四十四张，PPT共七十八页，创作于2022年6月06.04.2023（1）取三角形单元的节点位移为基本未知量：）取三角形单元的节点位移为基本未知量

45、：(a)其中，其中，称为单元的称为单元的节点位移列阵节点位移列阵；（2）应用插值公式，由）应用插值公式，由单元节点位移单元节点位移求出单元的求出单元的位移函数位移函数：(b)其中，其中，N 称为称为形函数矩阵形函数矩阵；（3）应用几何方程，由）应用几何方程，由单元的节点位移单元的节点位移求出求出单元的应变单元的应变：(c)其中，其中，B是表示是表示 与与 之间关系的矩阵；之间关系的矩阵；三角形单元离散化结构分析步骤三角形单元离散化结构分析步骤第四十五张，PPT共七十八页，创作于2022年6月06.04.2023 (f)其中，其中，Fe 是单元的节点力，是单元的节点力，k称为称为单元劲度列阵单元

46、劲度列阵；对三角形板单元，节点力为：对三角形板单元，节点力为：(e)（5）应用）应用虚功方程虚功方程，导出单元节点力与，导出单元节点力与节点位移之间的关系。对右图中的节点位移之间的关系。对右图中的i节点：节点：节点对单元的作用力为节点力，节点对单元的作用力为节点力，作用于单元上。作用于单元上。三角形单元离散化结构分析步骤三角形单元离散化结构分析步骤(续）续）（4）应用）应用物理方程物理方程，由单元的节点位移，由单元的节点位移求出单元的应力：求出单元的应力：(d)其中，其中，S称为称为应力转换矩阵应力转换矩阵；Fe是作用于单元的外力，此外，单元内部还作用有应力。根据虚功方程，从而得是作用于单元的

47、外力，此外，单元内部还作用有应力。根据虚功方程，从而得到节点力的公式：到节点力的公式：第四十六张，PPT共七十八页，创作于2022年6月06.04.2023（7）列出各节点的平衡方程，组成整个结构的平衡方程组。由于节点列出各节点的平衡方程，组成整个结构的平衡方程组。由于节点i受有环绕节点的受有环绕节点的单元移置而来的节点载荷单元移置而来的节点载荷 和节点力和节点力 因而因而i节点的平衡节点的平衡方程为：方程为：（i=1,2,n）(h)三角形单元离散化结构分析步骤三角形单元离散化结构分析步骤(续）续）（6）应用虚功方程，将单元中的外力载荷向节点移置，化）应用虚功方程，将单元中的外力载荷向节点移置

48、，化为为节点载荷节点载荷（即求出单元的节点载荷）：（即求出单元的节点载荷）：(g)将（将（f）代入（）代入（h），整理得：），整理得：（j)其中，其中，K称为称为整体刚度矩阵整体刚度矩阵，FL是是整体节点载荷列阵整体节点载荷列阵，是是整体节点位移列阵整体节点位移列阵。在上述求解步骤中，（在上述求解步骤中，（2）至（）至（6）是针对每个单元进行的，称为）是针对每个单元进行的，称为单元分析单元分析；（7）是针对整个结构进行的称为）是针对整个结构进行的称为整体分析整体分析。第四十七张，PPT共七十八页，创作于2022年6月06.04.2023n对三角形对三角形i,j,m三个节点，位移函数应当等于该节

49、点的位移值，即：三个节点，位移函数应当等于该节点的位移值，即：三角形单元的位移模式三角形单元的位移模式n对每个单元，只要求得单元中的位移函数，就可以应用对每个单元，只要求得单元中的位移函数，就可以应用几何方程几何方程求得应变，再应求得应变，再应用用物理方程物理方程求得应力。有限单元法中常取节点位移为基本未知量，由求得应力。有限单元法中常取节点位移为基本未知量，由单元的节点单元的节点位移位移求出求出单元中的位移函数单元中的位移函数是首先必须解决的问题。是首先必须解决的问题。n可以假定一个位移模式，来表示单元中的位移函数（即在单元中做出位移插值可以假定一个位移模式，来表示单元中的位移函数（即在单元

50、中做出位移插值函数）。三角形单元中，可以假定位移分量只是坐标的线性函数，即假定：函数）。三角形单元中，可以假定位移分量只是坐标的线性函数，即假定：n6个方程解出个方程解出1-6，代入，代入u，v式整理得：式整理得：其中：第四十八张，PPT共七十八页，创作于2022年6月06.04.2023三角形单元的位移模式三角形单元的位移模式nNi也可以写成为：也可以写成为：其中系数其中系数ai,bi,ci是：是：其中其中A就等于三角形就等于三角形ijm的面积：的面积：n 按照解析几何学，在图示的坐标系中，为了得出的面积按照解析几何学，在图示的坐标系中，为了得出的面积A不致成为负值，节点不致成为负值，节点i