平面问题有限元解法公式推导讲解.pptx

平面问题有限元解法公式推导讲解平面问题有限元解法公式推导讲解2023/3/162023/3/16有限单元法基本思想有限单元法的思想是将物体（连续的求解域）离散成有限个且按有限单元法的思想是将物体（连续的求解域）离散成有限个且按一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，从而将一个从而将一个连续的无限自由度连续的无限自由度连续的无限自由度连续的无限自由度问题简化为离散的问题简化为离散的有限自由度问题有限自由度问题有限自由度问题有限自由度问题求解的一种求解的一种数值分析法数值分析法数值分析法数值分析法。物体被离散后，通过对其中各个单元进。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。行单元分析，最终得到对整个物体的分析。有限单元法的分析步骤如下：有限单元法的分析步骤如下：l l物体离散化物体离散化l l单元特性分析单元特性分析l l单元组集，整体分析单元组集，整体分析l l求解未知节点的位移求解未知节点的位移l l由节点的位移求解各单元的位移和应力由节点的位移求解各单元的位移和应力第1页/共78页2023/3/162023/3/16有限元单元模型中几个重要概念有限元单元模型中几个重要概念单元单元l l网格划分中每一个小的块体网格划分中每一个小的块体节点节点l l确定单元形状、单元之间相互联确定单元形状、单元之间相互联结的点结的点节点力节点力l l单元上节点处的结构内力单元上节点处的结构内力载荷载荷l l作用在单元节点上的外力作用在单元节点上的外力（集中力、分布力）（集中力、分布力）约束约束pp限制某些节点的某些自由度限制某些节点的某些自由度弹性模量弹性模量（杨式模量）（杨式模量）E E泊松比泊松比（横向变形系数）（横向变形系数）密度密度单元单元单元单元载荷节点节点力约束第2页/共78页2023/3/162023/3/161.研究内容研究内容内容：内容：弹性体在外力或温度作用下的应力、弹性体在外力或温度作用下的应力、变形、位移等分布规律。变形、位移等分布规律。任务：任务：解决弹性体的强度、刚度、稳定性问题。解决弹性体的强度、刚度、稳定性问题。弹性力学的内容及基本假定弹性力学的内容及基本假定2.研究对象研究对象一般弹性实体结构：一般弹性实体结构：三维弹性固体、板状结构、杆件等三维弹性固体、板状结构、杆件等第3页/共78页2023/3/162023/3/16弹性力学的内容及基本假定弹性力学的内容及基本假定3.研究方法研究方法由平衡方程、几何方程、物理方程三方面分析由平衡方程、几何方程、物理方程三方面分析4.数学理论基础数学理论基础 偏微分方程（高阶，二、三个变量）偏微分方程（高阶，二、三个变量）数值解法数值解法：能量法（变分法）、差分：能量法（变分法）、差分法、有限单元法等。法、有限单元法等。第4页/共78页2023/3/162023/3/16弹性力学的内容及基本假弹性力学的内容及基本假定定5.基本假定基本假定（1）.连续性假定连续性假定整个物体的体积都被组成物体的介质充满，不留下任何空隙。整个物体的体积都被组成物体的介质充满，不留下任何空隙。作用：使得使得、u 等量表示成坐标的连续函数。等量表示成坐标的连续函数。第5页/共78页2023/3/162023/3/16弹性力学的内容及基本假定弹性力学的内容及基本假定（2）.完全弹性假定完全弹性假定 假定物体完全服从虎克（假定物体完全服从虎克（HookeHooke）定律，应力与应变间）定律，应力与应变间成线性比例关系。成线性比例关系。脆性材料脆性材料 一直到破坏前，都可近似为线弹性的；一直到破坏前，都可近似为线弹性的；塑性材料塑性材料 比例阶段，可视为线弹性的。比例阶段，可视为线弹性的。（3）.均匀性假定均匀性假定 假定整个物体是由同一种材料组成的，各部分材料性质相同。假定整个物体是由同一种材料组成的，各部分材料性质相同。作用：作用：弹性常数（弹性常数（E、）等）等不随位置坐标而变化；不随位置坐标而变化；取微元体分析的结果可应用于整个物体。取微元体分析的结果可应用于整个物体。第6页/共78页2023/3/162023/3/16弹性力学的内容及基本假定弹性力学的内容及基本假定（4）.各向同性假定各向同性假定（5）.小变形假定小变形假定 假定物体内一点的力学性质在所有各个方向都相同。假定物体内一点的力学性质在所有各个方向都相同。作用：作用：弹性常数（弹性常数（E、）不随坐标方向而变化；不随坐标方向而变化；假定位移和形变是微小的，即物体受力后物体内各点位移假定位移和形变是微小的，即物体受力后物体内各点位移远远小于物体的原来的尺寸。远远小于物体的原来的尺寸。作用：作用：建立方程时，可略去高阶微量；建立方程时，可略去高阶微量；可用变形前的尺寸代替变形后的尺寸。可用变形前的尺寸代替变形后的尺寸。使求解的方程使求解的方程线性化线性化。第7页/共78页2023/3/162023/3/16基本概念：基本概念：外力、应力、形变、位移。外力、应力、形变、位移。1.外力：外力：体力、面力体力、面力(1)体体力力 分布在物体分布在物体体积体积内的力内的力 体力分布集度体力分布集度（矢量）（矢量）xyzO单位：单位：N/m3kN/m3说明：说明：f 是坐标的连续分布函数是坐标的连续分布函数；弹性力学中的几个基本概念弹性力学中的几个基本概念p第8页/共78页2023/3/162023/3/16(2)面面力力 分布在物体表面的力分布在物体表面的力 面力分布集度（矢量）面力分布集度（矢量）xyzO单位：单位：1N/m2=1Pa(帕)1MN/m2=106Pa=1MPa(兆帕)说明：说明：弹性力学中的几个基本概念弹性力学中的几个基本概念是坐标的连续分布函数是坐标的连续分布函数；p第9页/共78页2023/3/162023/3/162.应力应力(1)一点应力的概念一点应力的概念AF内力内力(1)物体内部分子或原子间的相互物体内部分子或原子间的相互作用力作用力;(2)由于外力作用引起的相互作用力由于外力作用引起的相互作用力.(不考虑不考虑)P截面上截面上P点的应力点的应力应力矢量应力矢量.的极限方向的极限方向应力分应力分量量n(法线法线)应力的法向分应力的法向分量量 正应力正应力应力的切向分应力的切向分量量 切应力切应力单位单位:MPa(兆帕)应力关于坐标连续分布应力关于坐标连续分布弹性力学中的几个基本概念弹性力学中的几个基本概念第10页/共78页2023/3/162023/3/16(2)一点的应力状态一点的应力状态通过一点通过一点P 的各个面上应力状况的集合的各个面上应力状况的集合 称为一点的应力状态称为一点的应力状态x面的应力：面的应力：y面的应力：面的应力：z面的应力：面的应力：弹性力学中的几个基本概念弹性力学中的几个基本概念第11页/共78页2023/3/162023/3/16用矩阵表示：用矩阵表示：其中，只有其中，只有6个量独立。个量独立。切应力互等定切应力互等定理理应力应力正负号正负号的规定：的规定：正应力正应力 拉为正，压为负。拉为正，压为负。切应力切应力 坐标坐标正面正面上，与坐标正向一致时为正；上，与坐标正向一致时为正；坐标坐标负面负面上，与坐标正向相反时为正。上，与坐标正向相反时为正。xyzO弹性力学中的几个基本概念弹性力学中的几个基本概念第12页/共78页2023/3/162023/3/163.形变形变形变形变 物体的形状改物体的形状改变变xyzO（1）线段长度的改变）线段长度的改变（2）两线段间夹角的改变。）两线段间夹角的改变。PBCA用正应变用正应变度量度量切应变切应变度量度量（切应变（切应变两垂直线段夹角两垂直线段夹角（直角）（直角）的改变量）的改变量）三个方向的正应变：三个方向的正应变：三个平面内的切应变：三个平面内的切应变：(1)一点形变的度量一点形变的度量应变的正负：应变的正负：正应变正应变：伸长时为正，缩短时为负；伸长时为正，缩短时为负；切应变切应变：以直角变小时为正，变大时为负；以直角变小时为正，变大时为负；弹性力学中的几个基本概念弹性力学中的几个基本概念第13页/共78页2023/3/162023/3/16(2)一点应变状态一点应变状态其中其中应变无量纲；应变无量纲；4.位移位移注：注：一点的位移一点的位移 矢量矢量S应变分量均为位置坐标的函数应变分量均为位置坐标的函数xyzOSwuvP位移分量：位移分量：u x方向的位移方向的位移 分量；分量；v y方向的位移方向的位移 分量；分量；w z方向的位移方向的位移 分量。分量。量纲：量纲：m 或 mm弹性力学中的几个基本概念弹性力学中的几个基本概念第14页/共78页2023/3/162023/3/16 工程力学问题建立力学模型的过程中，一般从三方面进行简化：结构简化 如空间问题向平面问题的简化，向轴对称问题的简化，实体结构向板、壳结构的简化。受力简化 如：根据圣维南原理，复杂力系简化为等效力系等。材料简化根据各向同性、连续、均匀等假设进行简化。第15页/共78页2023/3/162023/3/16平面问题的基本理论平面问题的基本理论任何一个实际的弹性力学问题任何一个实际的弹性力学问题都是空间问题，但是如果所考都是空间问题，但是如果所考察的弹性体具有某种察的弹性体具有某种特殊的形特殊的形状，状，并且承受的是某些并且承受的是某些特殊的特殊的外力和约束外力和约束，就可以把空间问，就可以把空间问题简化为近似的题简化为近似的平面问题平面问题。两种典型的平面问题两种典型的平面问题l l平面应力问题平面应力问题平面应力问题平面应力问题l l平面应变问题平面应变问题平面应变问题平面应变问题第16页/共78页2023/3/162023/3/16平面应力问题平面应力问题(1)几何特征几何特征xyyztba 一个方向的尺寸比另两个一个方向的尺寸比另两个方向的尺寸小得多。方向的尺寸小得多。平板平板如：板式吊钩，旋转圆盘，工字形梁的腹板等如：板式吊钩，旋转圆盘，工字形梁的腹板等(2)受力特征受力特征外力外力（体力、面力）和（体力、面力）和约束约束，仅，仅平行于板面作用平行于板面作用，沿，沿 z 方向不变化。方向不变化。第17页/共78页2023/3/162023/3/16xyyztba(3)应力特征应力特征如图选取坐标系，以板的中面如图选取坐标系，以板的中面为为xy 平面，垂直于中面的任一直线平面，垂直于中面的任一直线为为 z 轴。轴。由于板面上不受力，有由于板面上不受力，有因板很薄，且外力因板很薄，且外力沿沿 z 轴方向不变。轴方向不变。可认为可认为整个薄板的整个薄板的各点各点都有：都有：由切应力互等定理，有由切应力互等定理，有结论：结论：平面应力问题只有三个应力分量：平面应力问题只有三个应力分量：xy应变分量、位移分量也仅为应变分量、位移分量也仅为 x、y 的函数，与的函数，与 z 无关。无关。第18页/共78页2023/3/162023/3/16平面应变问题平面应变问题(1)几何特征几何特征水水坝坝滚滚柱柱厚壁圆厚壁圆筒筒 一个方向的尺寸比另一个方向的尺寸比另两个方向的尺寸两个方向的尺寸大得多，大得多，且沿长度方向几何形状和且沿长度方向几何形状和尺寸不变化。尺寸不变化。近似认为无限长近似认为无限长(2)外力特征外力特征 外力外力（体力、面力）（体力、面力）平行于横截面平行于横截面作作用，且用，且沿长度沿长度 z 方向不变化方向不变化。约束约束 沿长度沿长度 z 方向不变化方向不变化。(3)变形特征变形特征 如图建立坐标系：以任一横截面为如图建立坐标系：以任一横截面为 xy 面，任一纵线为面，任一纵线为 z 轴。轴。设设 z方向为无限长，则方向为无限长，则沿沿 z 方向都不变化，方向都不变化，仅为仅为 x，y 的函数。的函数。任一横截面均可视为对称面任一横截面均可视为对称面第19页/共78页2023/3/162023/3/16水水坝坝任一横截面均可视为对称面，则有任一横截面均可视为对称面，则有所有各点的位移矢量都平行于所有各点的位移矢量都平行于 x y 平面。平面。平面位移问题平面位移问题 平面应变问题平面应变问题注：注：平面应变问题中平面应变问题中但是，但是，第20页/共78页2023/3/162023/3/16 如图所示三种情形，是否都属平面问题？是平如图所示三种情形，是否都属平面问题？是平面应力问题还是平面应变问题？面应力问题还是平面应变问题？平面应力问题平面应力问题平面应变问题平面应变问题非平面问题非平面问题第21页/共78页2023/3/162023/3/16三大基本方程三大基本方程三大基本方程三大基本方程根据根据根据根据静力学静力学静力学静力学、几何学几何学几何学几何学和和和和物理学物理学物理学物理学三方面条件，建立三套方程。三方面条件，建立三套方程。三方面条件，建立三套方程。三方面条件，建立三套方程。l l平面问题中，根据微分体的平衡条件，建立平面问题中，根据微分体的平衡条件，建立平面问题中，根据微分体的平衡条件，建立平面问题中，根据微分体的平衡条件，建立平衡微分方程平衡微分方程平衡微分方程平衡微分方程：(1-1)(1-1)l l根据微分线段上形变与位移之间的几何关系，建立根据微分线段上形变与位移之间的几何关系，建立根据微分线段上形变与位移之间的几何关系，建立根据微分线段上形变与位移之间的几何关系，建立几何方程几何方程几何方程几何方程：(1-2)(1-2)l l根据应力与形变之间的物理关系，建立根据应力与形变之间的物理关系，建立根据应力与形变之间的物理关系，建立根据应力与形变之间的物理关系，建立物理方程物理方程物理方程物理方程：(1-3)(1-3)(1-3)(1-3)第22页/共78页2023/3/162023/3/16平衡微分平衡微分方程方程从弹性体中取出一个微分体，根据平衡从弹性体中取出一个微分体，根据平衡从弹性体中取出一个微分体，根据平衡从弹性体中取出一个微分体，根据平衡条件导出应力分量与体力分量之间的关条件导出应力分量与体力分量之间的关条件导出应力分量与体力分量之间的关条件导出应力分量与体力分量之间的关系式，也就是平面问题的系式，也就是平面问题的系式，也就是平面问题的系式，也就是平面问题的平衡微分方程平衡微分方程平衡微分方程平衡微分方程。从弹性体中取出一个微小的正平行六面从弹性体中取出一个微小的正平行六面从弹性体中取出一个微小的正平行六面从弹性体中取出一个微小的正平行六面体，它在体，它在体，它在体，它在x x和和和和y y方向的尺寸分别为方向的尺寸分别为方向的尺寸分别为方向的尺寸分别为d dx x和和和和d dy y，在，在，在，在z z方向的尺寸为一个单位长度。方向的尺寸为一个单位长度。方向的尺寸为一个单位长度。方向的尺寸为一个单位长度。n以以x为投影轴，列出投影的平衡方程：为投影轴，列出投影的平衡方程：n约简以后，两边除以约简以后，两边除以dxdy，得：，得：n同理，以同理，以y为投影轴，列出投影的平衡方程，化简得为投影轴，列出投影的平衡方程，化简得：第23页/共78页2023/3/162023/3/16几何方程几何方程几何方程几何方程经过弹性体内的任意一点经过弹性体内的任意一点经过弹性体内的任意一点经过弹性体内的任意一点P P，沿，沿，沿，沿x x轴和轴和轴和轴和y y轴的正方向取两个微小长轴的正方向取两个微小长轴的正方向取两个微小长轴的正方向取两个微小长度的线段度的线段度的线段度的线段PAPAdxdx和和和和PBPBdydy。假。假。假。假定弹性体受力后，定弹性体受力后，定弹性体受力后，定弹性体受力后，P P,A A,B B三点分三点分三点分三点分别移动到别移动到别移动到别移动到PP,AA,B.B.n线段线段PA的正应变是：的正应变是：注注:由于位移微小，由于位移微小，y方向的位移方向的位移v引起的引起的PA的伸缩，是高一阶微量，略去不计。的伸缩，是高一阶微量，略去不计。n线段线段PB的正应变是：的正应变是：n线段线段PA与与 PB之间的直角的改变，即切应变之间的直角的改变，即切应变n线段线段PA的转角的转角是是：n线段线段PB的转角的转角是：是：第24页/共78页2023/3/162023/3/16物理方程物理方程在理想的弹性体中，形变分量和应力分量之间的关系，在材在理想的弹性体中，形变分量和应力分量之间的关系，在材在理想的弹性体中，形变分量和应力分量之间的关系，在材在理想的弹性体中，形变分量和应力分量之间的关系，在材料力学根据胡克定律导出如下：料力学根据胡克定律导出如下：料力学根据胡克定律导出如下：料力学根据胡克定律导出如下：n在平面应力问题中，在平面应力问题中，式变为：式变为：n在平面应变问题中，在平面应变问题中，只要将上式中的只要将上式中的E换为换为 ，换为换为 就得到平面应变问题的物理方程。就得到平面应变问题的物理方程。第25页/共78页2023/3/162023/3/16n假定已知任一点假定已知任一点P处坐标面上的应力分量处坐标面上的应力分量x，y，x y=y x。求经过该点的，平行于。求经过该点的，平行于z轴而倾斜于轴而倾斜于x轴和轴和 y轴的任何倾斜面上应力。轴的任何倾斜面上应力。n在在P点附近取一个平面点附近取一个平面AB，它平行于上述斜面，并经过，它平行于上述斜面，并经过P点划出一个微小的三棱柱点划出一个微小的三棱柱PAB。当。当AB无限小而趋于无限小而趋于P点时，平面点时，平面AB上的应力就成为斜面上的应力。上的应力就成为斜面上的应力。平面问题中一点的应力状态平面问题中一点的应力状态n设斜面设斜面AB 的长度为的长度为ds，则，则PB面及面及PA面的长度分别为面的长度分别为 lds及及mds，而，而PAB的面积为的面积为 ldsmds/2，棱柱的厚度设为，棱柱的厚度设为1。n由由x轴平衡条件，得：轴平衡条件，得：n其中，其中，fx为体力分量。将上式除以为体力分量。将上式除以ds，并令，并令ds趋于趋于0（斜面（斜面AB趋于趋于P点），点），即得：即得：n由由y轴平衡条件，得：轴平衡条件，得：n用用n表示斜面表示斜面AB的的外法线方向，其方向余弦为：外法线方向，其方向余弦为：第26页/共78页2023/3/162023/3/16边界条件边界条件q若在若在su部分部分边界上给定了约束位移分量边界上给定了约束位移分量 和和 ，则对于此边界上，则对于此边界上的每一点，位移函数的每一点，位移函数u和和v应满足条件：应满足条件：q其中其中（u）s 和和（v）s 是位移的边界值，是位移的边界值，和和 在边界上是坐标的在边界上是坐标的已知函数。已知函数。n边界条件表示在边界上位移与约束，或应力与面力之间的关系式。边界条件表示在边界上位移与约束，或应力与面力之间的关系式。它可以分为它可以分为位移边界条件位移边界条件、应力边界条件应力边界条件和和混合边界条件混合边界条件。q位移边界条件位移边界条件：q应力边界条件应力边界条件：q若在若在su部分部分边界上给定了面力边界上给定了面力 和和 ，则由平衡条件得出平面，则由平衡条件得出平面应力问题的应力应力问题的应力(或面力）边界条件为：或面力）边界条件为：其中，其中，l,m是边界面外法线的方向余是边界面外法线的方向余弦。弦。第27页/共78页2023/3/162023/3/16圣维南原理圣维南原理圣维南原理圣维南原理在求解弹性力学问题时，应力分量、形变分量和位移分量在求解弹性力学问题时，应力分量、形变分量和位移分量在求解弹性力学问题时，应力分量、形变分量和位移分量在求解弹性力学问题时，应力分量、形变分量和位移分量必须满足区域内的必须满足区域内的必须满足区域内的必须满足区域内的三套基本方程三套基本方程三套基本方程三套基本方程，还必须满足边界上的，还必须满足边界上的，还必须满足边界上的，还必须满足边界上的边边边边界条件界条件界条件界条件。但是，要使边界条件得到完全满足，往往遇到很。但是，要使边界条件得到完全满足，往往遇到很。但是，要使边界条件得到完全满足，往往遇到很。但是，要使边界条件得到完全满足，往往遇到很大的困难。大的困难。大的困难。大的困难。圣维南原理圣维南原理圣维南原理圣维南原理可为简化局部边界上的应力边界条件提供很大可为简化局部边界上的应力边界条件提供很大可为简化局部边界上的应力边界条件提供很大可为简化局部边界上的应力边界条件提供很大方便。方便。方便。方便。圣维南原理圣维南原理圣维南原理圣维南原理表明，如果把物体的一小部分边界上的面力，表明，如果把物体的一小部分边界上的面力，表明，如果把物体的一小部分边界上的面力，表明，如果把物体的一小部分边界上的面力，变换为分布不同但变换为分布不同但变换为分布不同但变换为分布不同但静力等效静力等效静力等效静力等效的面力（的面力（的面力（的面力（主矢相同，对同一点主矢相同，对同一点主矢相同，对同一点主矢相同，对同一点的主矩也相同的主矩也相同的主矩也相同的主矩也相同），那么，近处的应力分布将有显著的改变，），那么，近处的应力分布将有显著的改变，），那么，近处的应力分布将有显著的改变，），那么，近处的应力分布将有显著的改变，但是远处所受的影响可以不计。但是远处所受的影响可以不计。但是远处所受的影响可以不计。但是远处所受的影响可以不计。第28页/共78页2023/3/162023/3/16圣维南原理的应用圣维南原理的应用例，设有柱形构件，在两端截面的例，设有柱形构件，在两端截面的例，设有柱形构件，在两端截面的例，设有柱形构件，在两端截面的形心受到大小相等而方向相反的拉形心受到大小相等而方向相反的拉形心受到大小相等而方向相反的拉形心受到大小相等而方向相反的拉力力力力F F（a a）。如果把一端或两端的。如果把一端或两端的。如果把一端或两端的。如果把一端或两端的拉力变换为静力等效的力，则只有拉力变换为静力等效的力，则只有拉力变换为静力等效的力，则只有拉力变换为静力等效的力，则只有虚线划出的部分的应力分布有显著虚线划出的部分的应力分布有显著虚线划出的部分的应力分布有显著虚线划出的部分的应力分布有显著的改变，而其余部分所受影响是可的改变，而其余部分所受影响是可的改变，而其余部分所受影响是可的改变，而其余部分所受影响是可以不计的。以不计的。以不计的。以不计的。n由于（由于（d）图中，面力连续分布，边界条件简单，应力容易求得。其它）图中，面力连续分布，边界条件简单，应力容易求得。其它三种情况，应力难以求得。把三种情况，应力难以求得。把d情况下的应力解答应用到其它三个情况，情况下的应力解答应用到其它三个情况，虽不能满足两端的应力边界条件，但仍然可以表明离杆端较远处的应力虽不能满足两端的应力边界条件，但仍然可以表明离杆端较远处的应力状态，没有显著的误差。状态，没有显著的误差。n图图e，构件右端有位移边界条件，构件右端有位移边界条件，d情况的解答，情况的解答，不能满足位移边界条件，但不能满足位移边界条件，但e图右端的面力，一定是合成为经过截面形图右端的面力，一定是合成为经过截面形心的力心的力F。所以把图。所以把图d情况的解答应用于图情况的解答应用于图e时，仍然只是在靠近两端处时，仍然只是在靠近两端处有显著的误差，而在离两端较远之处，误差可以不计。有显著的误差，而在离两端较远之处，误差可以不计。第29页/共78页2023/3/162023/3/16圣维南原理的应用圣维南原理的应用例，厚度例，厚度例，厚度例，厚度=1=1的梁中，左右两端的梁中，左右两端的梁中，左右两端的梁中，左右两端x=x=l l，的边界面是正、负，的边界面是正、负，的边界面是正、负，的边界面是正、负x x面，其上作用面，其上作用面，其上作用面，其上作用有一般分布的面力有一般分布的面力有一般分布的面力有一般分布的面力 。按照严格的应力边界条件，应力分。按照严格的应力边界条件，应力分。按照严格的应力边界条件，应力分。按照严格的应力边界条件，应力分量在边界上满足：量在边界上满足：量在边界上满足：量在边界上满足：n上式要求在边界上上式要求在边界上y值不同的各点，应力分量与对应的面力分量必须处值不同的各点，应力分量与对应的面力分量必须处处相等，这种严格的条件是较难满足的。处相等，这种严格的条件是较难满足的。n当当lh时，时，x=l 是梁的边界的一小部分，可以应用是梁的边界的一小部分，可以应用圣维南原理，利用圣维南原理，利用静力等效条件来代替，即，使静力等效条件来代替，即，使应力的主矢量应力的主矢量和和主矩主矩分别等于分别等于对应的对应的面面力的主矢量力的主矢量和和主矩主矩。第30页/共78页2023/3/162023/3/16圣维南原理的应用圣维南原理的应用n应力的主矢量和主矩的绝对值分别等于面力的主矢量和主矩的绝对值；应力的主矢量和主矩的绝对值分别等于面力的主矢量和主矩的绝对值；n面力的主矢量和主矩的方向就是应力的主矢量和主矩的方向。面力的主矢量和主矩的方向就是应力的主矢量和主矩的方向。第31页/共78页2023/3/162023/3/16按位移法求解按位移法求解平面问题平面问题以上几节已经建立了弹性力学平面问题的以上几节已经建立了弹性力学平面问题的以上几节已经建立了弹性力学平面问题的以上几节已经建立了弹性力学平面问题的 基本方程和边基本方程和边基本方程和边基本方程和边界条件，即：界条件，即：界条件，即：界条件，即：平衡微分方程、几何方程平衡微分方程、几何方程平衡微分方程、几何方程平衡微分方程、几何方程和和和和物理方程物理方程物理方程物理方程，以，以，以，以及及及及位移的边界条件位移的边界条件位移的边界条件位移的边界条件和和和和应力的边界条件应力的边界条件应力的边界条件应力的边界条件。求解弹性力学平面问题即求解求解弹性力学平面问题即求解求解弹性力学平面问题即求解求解弹性力学平面问题即求解3 3个应力分量个应力分量个应力分量个应力分量、3 3个形变分个形变分个形变分个形变分量量量量及及及及2 2个位移分量个位移分量个位移分量个位移分量的未知函数。通常采用类似于代数方程的未知函数。通常采用类似于代数方程的未知函数。通常采用类似于代数方程的未知函数。通常采用类似于代数方程中消元法进行求解。中消元法进行求解。中消元法进行求解。中消元法进行求解。按位移求解的方法，称为按位移求解的方法，称为按位移求解的方法，称为按位移求解的方法，称为位移法位移法位移法位移法。它以位移分量为基本。它以位移分量为基本。它以位移分量为基本。它以位移分量为基本未知函数。未知函数。未知函数。未知函数。按应力求解的方法，称为按应力求解的方法，称为按应力求解的方法，称为按应力求解的方法，称为应力法应力法应力法应力法。它以应力分量为基本。它以应力分量为基本。它以应力分量为基本。它以应力分量为基本未知函数。未知函数。未知函数。未知函数。第32页/共78页2023/3/162023/3/16按位移法求解平面问题按位移法求解平面问题平面问题中，取位移分量平面问题中，取位移分量平面问题中，取位移分量平面问题中，取位移分量u u和和和和v v为基本未知函数。为基本未知函数。为基本未知函数。为基本未知函数。从方程中消去形变分量和应力分量：从方程中消去形变分量和应力分量：从方程中消去形变分量和应力分量：从方程中消去形变分量和应力分量：n将几何方程代入上式将几何方程代入上式n利用平衡微分方程和边界条件，导出用位移表示的平衡微分方程：利用平衡微分方程和边界条件，导出用位移表示的平衡微分方程：第33页/共78页2023/3/162023/3/16按位移法求解平面问题按位移法求解平面问题利用应力边界条件利用应力边界条件利用应力边界条件利用应力边界条件n得到用位移表示的应力边界条件得到用位移表示的应力边界条件其中：其中：n位移边界条件如（位移边界条件如（1-4）不变）不变n按位移法求解平面应力问题时，要使位移分量在区域内满足平衡微分按位移法求解平面应力问题时，要使位移分量在区域内满足平衡微分方程，在边界上满足位移边界条件或应力边界条件。方程，在边界上满足位移边界条件或应力边界条件。第34页/共78页2023/3/162023/3/16按位移法求解平面问题按位移法求解平面问题（例题）（例题）设有如图所示的杆件，在设有如图所示的杆件，在设有如图所示的杆件，在设有如图所示的杆件，在y y方向的上端为固定，方向的上端为固定，方向的上端为固定，方向的上端为固定，而下端为自由，受自重体力而下端为自由，受自重体力而下端为自由，受自重体力而下端为自由，受自重体力f fx x0 0，f fy ygg的作用。的作用。的作用。的作用。试用位移法求解此问题。试用位移法求解此问题。试用位移法求解此问题。试用位移法求解此问题。n解：将这个问题简化为一维问题处理。解：将这个问题简化为一维问题处理。n设设u=0,v=v(y)，泊松比，泊松比0。代入位移表示的平衡微。代入位移表示的平衡微分方程，得：分方程，得：n第一式自然满足，第二式成为：第一式自然满足，第二式成为：n解出：解出：第35页/共78页2023/3/162023/3/16按位移法求解平面问题（例题）按位移法求解平面问题（例题）设有左图所示的杆件，在设有左图所示的杆件，在设有左图所示的杆件，在设有左图所示的杆件，在y y方向的上端为固定，而下方向的上端为固定，而下方向的上端为固定，而下方向的上端为固定，而下端为自由，受自重体力端为自由，受自重体力端为自由，受自重体力端为自由，受自重体力f fx x0 0，f fy ygg的作用。试用位的作用。试用位的作用。试用位的作用。试用位移法求解此问题。移法求解此问题。移法求解此问题。移法求解此问题。n解出：解出：n上下边的边界条件分别要求：上下边的边界条件分别要求：n将（将（a）式代入（）式代入（b）式得：）式得：B0，n再代入（再代入（c）式，即得：）式，即得：n得到解答：得到解答：第36页/共78页2023/3/162023/3/16有限元单元法分析步骤（一）有限元单元法分析步骤（一）结构离散化结构离散化结构离散化结构离散化l l 将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点连接。结构的离散化是有限元法分析的第一步，关系到计算精度和效率，连接。结构的离散化是有限元法分析的第一步，关系到计算精度和效率，连接。结构的离散化是有限元法分析的第一步，关系到计算精度和效率，连接。结构的离散化是有限元法分析的第一步，关系到计算精度和效率，包括以下三个方面：包括以下三个方面：包括以下三个方面：包括以下三个方面：单元类型的选择单元类型的选择单元类型的选择单元类型的选择。选定单元类型，确。选定单元类型，确。选定单元类型，确。选定单元类型，确定单元形状、单元节点数、节点自由定单元形状、单元节点数、节点自由定单元形状、单元节点数、节点自由定单元形状、单元节点数、节点自由度数等。度数等。度数等。度数等。单元划分单元划分单元划分单元划分。网格划分越细，节点越多，。网格划分越细，节点越多，。网格划分越细，节点越多，。网格划分越细，节点越多，计算结果越精确，但计算量越大。网计算结果越精确，但计算量越大。网计算结果越精确，但计算量越大。网计算结果越精确，但计算量越大。网格加密到一定程度后计算精度提高就格加密到一定程度后计算精度提高就格加密到一定程度后计算精度提高就格加密到一定程度后计算精度提高就不明显，对应应力变化平缓区域不必不明显，对应应力变化平缓区域不必不明显，对应应力变化平缓区域不必不明显，对应应力变化平缓区域不必要细分网格。要细分网格。要细分网格。要细分网格。节点编码节点编码节点编码节点编码。注意：注意：注意：注意：有限元分析的结构有限元分析的结构有限元分析的结构有限元分析的结构已不是已不是已不是已不是原有的物体或结构物原有的物体或结构物原有的物体或结构物原有的物体或结构物，而是由同样，而是由同样，而是由同样，而是由同样材料、众多单元以一定方式连接成的离散物体。用有限元分析计材料、众多单元以一定方式连接成的离散物体。用有限元分析计材料、众多单元以一定方式连接成的离散物体。用有限元分析计材料、众多单元以一定方式连接成的离散物体。用有限元分析计算所获得的结果是近似的（满足工程要求即可）。算所获得的结果是近似的（满足工程要求即可）。算所获得的结果是近似的（满足工程要求即可）。算所获得的结果是近似的（满足工程要求即可）。平面问题有限单元法基本概念平面问题有限单元法基本概念第37页/共78页2023/3/162023/3/16有限元单元法分析步骤（二）有限元单元法分析步骤（二）有限元单元法分析步骤（二）有限元单元法分析步骤（二）单元特性分析单元特性分析单元特性分析单元特性分析l l 选择未知量模式选择未知量模式选择未知量模式选择未知量模式 选择节点位移作为基本未知量时，称选择节点位移作为基本未知量时，称选择节点位移作为基本未知量时，称选择节点位移作为基本未知量时，称为位移法；为位移法；为位移法；为位移法；选节点力作为基本未知量时，称为力选节点力作为基本未知量时，称为力选节点力作为基本未知量时，称为力选节点力作为基本未知量时，称为力法；法；法；法；取一部分节点位移和一部分节点力作取一部分节点位移和一部分节点力作取一部分节点位移和一部分节点力作取一部分节点位移和一部分节点力作为未知量，称为混合法。为未知量，称为混合法。为未知量，称为混合法。为未知量，称为混合法。l l分析单元力学性质分析单元力学性质分析单元力学性质分析单元力学性质 根据单元材料性质、形状、尺寸、节根据单元材料性质、形状、尺寸、节根据单元材料性质、形状、尺寸、节根据单元材料性质、形状、尺寸、节点数目、位置等，找出单元节点力和点数目、位置等，找出单元节点力和点数目、位置等，找出单元节点力和点数目、位置等，找出单元节点力和节点位移关系式，应用几何方程和物节点位移关系式，应用几何方程和物节点位移关系式，应用几何方程和物节点位移关系式，应用几何方程和物理方程建立力和位移的方程式，从而理方程建立力和位移的方程式，从而理方程建立力和位移的方程式，从而理方程建立力和位移的方程式，从而导出单元刚度矩阵。导出单元刚度矩阵。导出单元刚度矩阵。导出单元刚度矩阵。l l计算等效节点力计算等效节点力计算等效节点力计算等效节点力 作用在单元边界上的表面力、体积力作用在单元边界上的表面力、体积力作用在单元边界上的表面力、体积力作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，或集中力都需要等效地移到节点上去，或集中力都需要等效地移到节点上去，或集中力都需要等效地移到节点上去，即用等效力来替代所有作用在单元上即用等效力来替代所有作用在单元上即用等效力来替代所有作用在单元上即用等效力来替代所有作用在单元上的力。的力。的力。的力。第38页/共78页2023/3/162023/3/16有限元单元法分析步骤（三）有限元单元法分析步骤（三）整体分析整体分析整体分析整体分析l l集成整体节点载荷矢量集成整体节点载荷矢量集成整体节点载荷矢量集成整体节点载荷矢量 F F。结构离散化后，单元之间通过节点传。结构离散化后，单元之间通过节点传。结构离散化后，单元之间通过节点传。结构离散化后，单元之间通过节点传递力，作用在单元边界上的表面力、体积力或集中力都需要等效递力，作用在单元边界上的表面力、体积力或集中力都需要等效递力，作用在单元边界上的表面力、体积力或集中力都需要等效递力，作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体节点编码顺序组集成整体节点载荷矢量。节点编码顺序组集成整体节点载荷矢量。节