弹性力学有限元法3PPT课件

1、2022-3-141第第3章章 弹性力学有限元法弹性力学有限元法2022-3-1423.1有限元法求解问题的基本步骤有限元法求解问题的基本步骤1.问题及求解域定义问题及求解域定义2.连续体离散化连续体离散化 即即有限元网格划分有限元网格划分，将连续体划分为有限个具，将连续体划分为有限个具有一定形状的单元组合体，相邻单元之间通过节点有一定形状的单元组合体，相邻单元之间通过节点相连接。相连接。3.单元分析单元分析（1）选择位移模式）选择位移模式位移法位移法：选择节点位移作为基本未知量。（应用较多）：选择节点位移作为基本未知量。（应用较多）力法力法：选择节点力为基本未知量。：选择节点力为基本未知量。

2、混合法混合法：取一部分力和一部分节点位移作为基本未知量。：取一部分力和一部分节点位移作为基本未知量。2022-3-143（2）分析单元的力学性质）分析单元的力学性质 列出单元节点和节点位移之间的关系式。应用几何方程和列出单元节点和节点位移之间的关系式。应用几何方程和物理方程来建立力和位移的方程式，导出单元刚度矩阵。物理方程来建立力和位移的方程式，导出单元刚度矩阵。节点载荷和节点位移之间的关系式为：节点载荷和节点位移之间的关系式为： 为单元刚度矩阵。为单元刚度矩阵。（3）计算等效节点力计算等效节点力：用等效的节点力来代替所有在单元：用等效的节点力来代替所有在单元上的力。上的力。eKeeeFK20

3、22-3-1444. 组成物体的整体方程组组成物体的整体方程组 由单元刚度矩阵构成整体刚度矩阵。对总体建立方程：由单元刚度矩阵构成整体刚度矩阵。对总体建立方程：5. 求解有限元方程和结果解释求解有限元方程和结果解释 根据边界条件和初始条件求解上式，得到节点位移。根据边界条件和初始条件求解上式，得到节点位移。FKu2022-3-1453.2 连续体离散化连续体离散化 结构的离散化也称为结构的离散化也称为有限元网格划分有限元网格划分，即将求解域近似，即将求解域近似为具有不同有限大小和形状且只在节点上彼此相连的有限个为具有不同有限大小和形状且只在节点上彼此相连的有限个单元组成的离散域。单元组成的离散

4、域。 常用的单元类型：常用的单元类型：1.杆单元杆单元 一维单元，位移仅是轴向座标的函数。一维单元，位移仅是轴向座标的函数。2.平面单元平面单元 二维单元，单元内任意点的应力、应变、位移仅与两个二维单元，单元内任意点的应力、应变、位移仅与两个座标方向的变量有关。座标方向的变量有关。2022-3-1462022-3-147 三角形单元：采用三角形单元：采用线性位移模式，线性位移模式，在整个单元内各点的在整个单元内各点的应变值为常数，所以也称为常应变单元或常应力单元。应变值为常数，所以也称为常应变单元或常应力单元。 矩形单元：采用矩形单元：采用双线性位移模式双线性位移模式，单元内的应力是线性，单元

5、内的应力是线性变化的。变化的。3. 薄板弯曲单元和薄板单元薄板弯曲单元和薄板单元ukxuxK2ukxmx2()kxmxx4. 多面体单元多面体单元2022-3-1482022-3-1495. 等参数单元等参数单元：单元内：单元内任一点的位移与节点位移之间任一点的位移与节点位移之间的关系的关系恰好和恰好和该点的坐标与节点坐标之间该点的坐标与节点坐标之间的关系相同。的关系相同。 任意四边形的边一般不平行于坐标轴，沿单元边的位任意四边形的边一般不平行于坐标轴，沿单元边的位移将按抛物线变化，而不是线性变化。移将按抛物线变化，而不是线性变化。 以直角坐标系以直角坐标系XOY下的任意直边四边形单元单元的形

6、心下的任意直边四边形单元单元的形心为坐标原点，用等分它四个边的两条直线为坐标轴，建立一为坐标原点，用等分它四个边的两条直线为坐标轴，建立一个非正交的局部座标系个非正交的局部座标系 ，使单元边界上的，使单元边界上的 、 是是 ，这样在局部坐标系中构成一个矩形单元。矩形单元的节点和这样在局部坐标系中构成一个矩形单元。矩形单元的节点和内部任一点都与原总体坐标系中的单元的节点和内部点形成内部任一点都与原总体坐标系中的单元的节点和内部点形成一一对应关系。一一对应关系。总体坐标系适用于整个结构，局部坐标系只总体坐标系适用于整个结构，局部坐标系只适用于具体某个单元。适用于具体某个单元。1o 12022-3-

7、14106. 轴对称单元：轴对称单元： 几何形状是回转体，所受约束和外力对称于回转轴几何形状是回转体，所受约束和外力对称于回转轴的的机械机构称为轴对称问题。对此类问题一般采取柱坐标系机械机构称为轴对称问题。对此类问题一般采取柱坐标系来描述应力和变形。来描述应力和变形。 对于此类问题采用轴对称单元。对于此类问题采用轴对称单元。划分网格的基本原则：划分网格的基本原则：（1）网格数量：网格数量增加，计算精度会有所提高。）网格数量：网格数量增加，计算精度会有所提高。（2）网格疏密：在结构不同处采用不同的网格形式。）网格疏密：在结构不同处采用不同的网格形式。（3）单元阶次：网格数量较少时，计算精度差别较

8、大，采）单元阶次：网格数量较少时，计算精度差别较大，采用高阶单元。网格数量较多时，采用两种单元的精度相差不用高阶单元。网格数量较多时，采用两种单元的精度相差不大，采用低阶单元计算量降低。大，采用低阶单元计算量降低。2022-3-1411（4）网格质量：网格几何形状的合理性，网格质量）网格质量：网格几何形状的合理性，网格质量的好坏会影响计算精度，对于太差的网格形状程序的好坏会影响计算精度，对于太差的网格形状程序将会自动停止计算。将会自动停止计算。（5）网格分界面和分界点：结构中一些特殊位置的）网格分界面和分界点：结构中一些特殊位置的界面或特殊位置的点应分为网格边界或节点。界面或特殊位置的点应分为

9、网格边界或节点。（6）位移协调性：一个单元的节点必须也是相邻单）位移协调性：一个单元的节点必须也是相邻单元的节点，只有这样单元上的力和力矩才能够通过元的节点，只有这样单元上的力和力矩才能够通过节点传递到相邻单元。节点传递到相邻单元。（7）网格布局：对于对称结构应该划分对称单元。）网格布局：对于对称结构应该划分对称单元。（8）节点和单元编号：一般情况下程序自动编号。）节点和单元编号：一般情况下程序自动编号。2022-3-14123.3 单元分析单元分析1.单元的插值函数单元的插值函数 如果弹性体内的位移分量已知，则应变分量和如果弹性体内的位移分量已知，则应变分量和应力分量也可以确定了。应力分量也

10、可以确定了。 几何方程几何方程 虎克定律虎克定律FKuuE2022-3-1413 对于整体划分单元后，在每个单元的局部范围里可以采对于整体划分单元后，在每个单元的局部范围里可以采用比较简单的函数来近似地表达单元的真实位移，把各单元用比较简单的函数来近似地表达单元的真实位移，把各单元的位移函数连接起来，就可以近似表示整个区域的真实的位的位移函数连接起来，就可以近似表示整个区域的真实的位移函数。移函数。2022-3-1414 在离散体中任取一个单元，三个节点按逆时针方向顺序编在离散体中任取一个单元，三个节点按逆时针方向顺序编号为号为i，j，m。节点坐标分别表示为（。节点坐标分别表示为（xi，yi）

11、，（），（xj，yj），），（xm，ym）。）。2022-3-1415 对于弹性力学平面问题，一个三角形单元上的每对于弹性力学平面问题，一个三角形单元上的每个节点应有个节点应有2个位移分量，则三角形单元共有个位移分量，则三角形单元共有6个自个自由度：由度： 。三角形单元的节点位移矢量是：三角形单元的节点位移矢量是：单元节点力矢量是：单元节点力矢量是：,iijjmmu v u v uv( ,)eTiijjmmu v u v uv (,)eeeeeee TiijjmmFXYXYXY2022-3-1416 单元分析的基本任务单元分析的基本任务是建立单元节点力与节点位移之间的是建立单元节点力与节点位移

12、之间的关系式：关系式：式中式中 是是6*6的矩阵，称为的矩阵，称为单元刚度矩阵单元刚度矩阵。 将单元的位移分量将单元的位移分量u，v取为坐标取为坐标x，y的多项式，且位移的多项式，且位移场函数场函数u，v在三个节点处的数值应该等于这三个节点处的六在三个节点处的数值应该等于这三个节点处的六个位移分量。个位移分量。即有：即有：eeeFKeK2022-3-1417123456( , )( , )u x yaa xa yv x yaa xa y在在i，j，m三点应该有：三点应该有：123456123456123456iiiiiijjjjjjmmmmmmuaa xa yvaa xa yuaa xa yv

13、aa xa yuaa xa yvaa xa y(1)2022-3-1418由上式可以确定由上式可以确定 的值。将其带入的值。将其带入（1）式就可以得到用）式就可以得到用单元节点位移表示的单单元节点位移表示的单元位移模式元位移模式。126,a aa( , )0( , )0( , )0( , )0( , )0( , )0( , )( , )jimjimNx yN x yNx yu x yNx yN x yNx yv x yeuN2022-3-1419 N称为形函数矩阵或插值函数矩阵。称为形函数矩阵或插值函数矩阵。插值函数具有如下性质：插值函数具有如下性质：（1）在节点上插值函数的值有：）在节点上插

14、值函数的值有：（2）在单元内任一点各插值函数的和等于一。）在单元内任一点各插值函数的和等于一。1( , )0jiiijN x y当j=i（i,j,m)当1ijmNNN2022-3-1420 线性单元线性单元 平面单元平面单元 立体单元立体单元 总体、局部和自然坐标总体、局部和自然坐标 数值积分：高斯勒让德多项式数值积分：高斯勒让德多项式 ANSYS实例实例2022-3-14213.1 一维单元一维单元1. 形函数形函数形函数在有限元分析中形函数在有限元分析中,扮演非常主要的扮演非常主要的角色角色. 除作为元素除作为元素(单元单元)的内插函数的内插函数,将元将元素内的位移或温度分布素内的位移或温

15、度分布,以节点位移或节以节点位移或节点温度表示外点温度表示外,在余量法中的迦辽金法中在余量法中的迦辽金法中,亦可作为加权函数来用亦可作为加权函数来用. 此外此外,亦可将分布亦可将分布载荷转换为集中力与弯矩载荷转换为集中力与弯矩,分别施加在各分别施加在各节点上节点上.形函数根据其多项式的幂次形函数根据其多项式的幂次,分为一次、分为一次、二次、三次与高次等。二次、三次与高次等。2022-3-1422位移沿着单元的分布可以用一个线性函数近似。如图所示。位移沿着单元的分布可以用一个线性函数近似。如图所示。u( ) eixxjxjuiul2022-3-1423 一维一次元素的形函数中，函数值沿单一坐标轴

16、一维一次元素的形函数中，函数值沿单一坐标轴以线性变化。假设位移函数沿以线性变化。假设位移函数沿x轴线性变化，位轴线性变化，位移函数移函数uu(x)可写成可写成: u=a1+a2x向量形式向量形式: 假设在假设在i、j节点的位移值分别为节点的位移值分别为ui和和uj , 有有: u=ui 在在X=Xi处处 u=uj 在在X=Xj处处 211aaxu2022-3-1424将节点的值带入线性方程将产生两个方程和两个未知量：将节点的值带入线性方程将产生两个方程和两个未知量：求解未知量求解未知量a1和和a2得到：得到：12iiuaa X12jjuaa X1ijjijiu Xu XaXX2jijiuuaX

17、X2022-3-1425由节点的值表示的单元的位移分布为：由节点的值表示的单元的位移分布为：改写一下形式得到：改写一下形式得到：定义形函数：定义形函数： 其中其中l为单元长度。为单元长度。( )ijjijiejijiu Xu XuuuXXXXX( )jeiijjijiXXXXuuuXXXXjjijiXXXXsXXliijjiXXXXsXXl2022-3-1426由形函数表示的单元的位移分布为：由形函数表示的单元的位移分布为：写成矩阵形式：写成矩阵形式：可以使用形函数和相应的节点值来表示给定单元上的任意的未知量的可以使用形函数和相应的节点值来表示给定单元上的任意的未知量的变化。变化。( )eii

18、jjusus u( )ieijjuussu2022-3-14272. 形函数的性质形函数的性质1. 在相应的节点上值为在相应的节点上值为1,而在相邻节点上值为而在相邻节点上值为0. 和和 和和2.形函数的和为形函数的和为1。1iijjiiXXXXXXXXsll0jjjjjiXXXXXXXXsll0iiiiijXXXXXXXXsll1jjjiijXXXXXXXXsll1jiijjijiXXXXssXXXX2022-3-14283. 形函数对于形函数对于x导数的和为导数的和为011()()0jijijijijiXXXXdddXXXdXXXXXXX 2022-3-14293. 实例实例（a）悬臂梁在

19、）悬臂梁在X4cm处的温度由单元（处的温度由单元（2）来表示：）来表示：(2)(2)(2)32223323XXXXTSTSTTTll5442413436.333TC1234(1)(2)(3)2cm3cm5cmx123450413420TTCTT2022-3-1430（b）悬臂梁在）悬臂梁在X=8cm处的温度由单元（处的温度由单元（3）来表示：）来表示：对于这个例子，注意对于这个例子，注意 和和 的区别。的区别。(3)(3)(3)34334434XXXXTSTSTTTll10885342025.655TC(2)3S(3)3S2022-3-1431 假设承受的是轴向负荷假设承受的是轴向负荷, 应用

20、线性单元应用线性单元, 柱体的柱体的垂直位移由下式确定垂直位移由下式确定.AB2022-3-1432（a）应用总体坐标）应用总体坐标Y，点，点A的位移由单元（的位移由单元（1）表示：）表示：（b）点）点B的位移由单元（的位移由单元（4）表示：）表示：(1)(1)(1)21112212YYYYuS uS uuull15 1010000.032830.021881515uin (4)(4)(4)54445545YYYYuSuSuuull605252450.075040.084420.079411515uin2022-3-14333.2 一维三节点单元一维三节点单元 一维三节点单元一维三节点单元用二

21、次函数代替线性用二次函数代替线性函数要求使用三个节函数要求使用三个节点来定义一个单元点来定义一个单元,这这是因为至少要有三个是因为至少要有三个点才能确定一个二次点才能确定一个二次函数函数.第三个点可以取第三个点可以取在单元的中点在单元的中点.2321xaxaau 32121aaaxxu2022-3-1434 如单元的温度分布可如单元的温度分布可以表示以表示为为:并且节点的值为并且节点的值为:T=Ti 在在X=Xi处处T=Tk 在在X=Xk处处T=Tj 在在X=Xj处处 2321XcXccTe 2022-3-1435 产生的三个方程和三个未知量产生的三个方程和三个未知量:求解求解c1，c2和和c

22、3，整理后得到由节点的值和形函数表示的单，整理后得到由节点的值和形函数表示的单元温度分布：元温度分布：将以上表达式写成矩阵形式为：将以上表达式写成矩阵形式为：212321232123iiikkkjjjTcc Xc XTcc Xc XTcc Xc X( )eiijjkkTsTs Ts T2022-3-1436这里形函数为这里形函数为：一般来说，对于给定单元，由节点的值表示的参数变化一般来说，对于给定单元，由节点的值表示的参数变化 可以写为：可以写为：( )ieijkjkTTsssTT2222()()2()()2()()ijkjikkijsXXXXlsXXXXlsXXXXl( )ieijkjkss

23、s2022-3-1437 对于给定单元节点的对于给定单元节点的位移变化位移变化: 在形函数的性质中在形函数的性质中,二二次形函数关于次形函数关于X的导数的导数之和不为零之和不为零. kjikjieuuuSSSu2022-3-14383.3 一维四节点单元一维四节点单元 在有限元公式中，二次插值提供了较为精确的结果。然在有限元公式中，二次插值提供了较为精确的结果。然而，如果需要更高的精度，可以使用更高阶的插值函数，例而，如果需要更高的精度，可以使用更高阶的插值函数，例如三次多项式。这样可以使用三次函数表示给定变量的空间如三次多项式。这样可以使用三次函数表示给定变量的空间变化。用三次函数代替二次函

24、数，要求使用四个节点来定义变化。用三次函数代替二次函数，要求使用四个节点来定义一个单元，这是因为至少要有四个节点才能确定一个三阶多一个单元，这是因为至少要有四个节点才能确定一个三阶多项式。单元被分成等长的三段。四个节点的取法如图所示。项式。单元被分成等长的三段。四个节点的取法如图所示。应用三次近似考虑上例，典型单元的温度分布可以表示为：应用三次近似考虑上例，典型单元的温度分布可以表示为：2022-3-1439iXjXkXmXmTkTjTiTTX3l3l3l2022-3-1440并且节点的值为：并且节点的值为： 在在 处处 在在 处处 在在 处处 在在 处处由节点的值和形函数表示的单元温度分布：

25、由节点的值和形函数表示的单元温度分布：写成矩阵的形式为：写成矩阵的形式为：( )231234eTcc Xc Xc XikmTTTTTTTTikmjXXXXXXXX( )eiijjkkmmTSTS TS TS T( )jeijkmkmTTTSSSSTT2022-3-1441 形函数为形函数为:)()(227)()(227)()(29)()(293333kjimmjikmkijmkjiXXXXXXlSXXXXXXlSXXXXXXlSXXXXXXlS 2022-3-1442 当插值函数的阶数增加时当插值函数的阶数增加时,可以用拉格朗日可以用拉格朗日插值代替以上的方法来得到形函数插值代替以上的方法来得

26、到形函数: 对于三次插值函数对于三次插值函数,每个节点相关的形函数每个节点相关的形函数可以用三个函数的乘积表示可以用三个函数的乘积表示.函数的乘积在函数的乘积在给定的节点上为给定的节点上为1,而在其他节点上为而在其他节点上为0. 如考虑节点如考虑节点i:2022-3-1443 若将若将X=Xj X=Xk X=Xm代入方程代入方程, 函数函数Si的值为零的值为零. 当在给定节点上计算形函数时当在给定节点上计算形函数时,即即X=Xi时时,函数的函数的值为值为1:)32)(3)()()(111lllaXXXXXXamikiji )()(1mkjiXXXXXXaS 3129la 2022-3-1444

27、 由节点值表示的三次插值函数任意参数由节点值表示的三次插值函数任意参数的变化可的变化可以表示为以表示为: 三次形函数性质三次形函数性质: 形函数在相应节点上值为形函数在相应节点上值为1,而在另一个相邻节点上值为而在另一个相邻节点上值为0; 如果对形函数求和如果对形函数求和,结果为结果为1; 对三次形函数的求导将得到二次的结果对三次形函数的求导将得到二次的结果. mkjimkjieSSSS2022-3-14453.4整体、局部和自然坐标整体、局部和自然坐标 在有限元建模中，可以使用几个参考系。整体坐标用来表示每个在有限元建模中，可以使用几个参考系。整体坐标用来表示每个节点的位置和每个单元的方向，

28、并用来施加边界条件和负荷。另一方节点的位置和每个单元的方向，并用来施加边界条件和负荷。另一方面，需要使用局部和自然坐标系，以简化计算。对于一维单元，整体面，需要使用局部和自然坐标系，以简化计算。对于一维单元，整体坐标坐标X和局部坐标和局部坐标x的关系为：的关系为：iXXx2022-3-1446 在一维一次形函数的表达式中带入由局部坐标在一维一次形函数的表达式中带入由局部坐标x表示的表示的X有：有：()1jjiiXXXXxxllls ()iiijXXXxXxllls0 xl2022-3-14473.4.1 一维线性自然坐标一维线性自然坐标 自然坐标是局部坐标的无量纲形式。使用自然坐标容易在上限自

29、然坐标是局部坐标的无量纲形式。使用自然坐标容易在上限1和下限和下限1之间积分。令：之间积分。令：这里这里x是局部坐标，局部坐标与自然坐标的关系如图所示。是局部坐标，局部坐标与自然坐标的关系如图所示。21xl2022-3-1448 通过将由通过将由 表示的表示的x带入形函数的表达式能够带入形函数的表达式能够得到自然线性形函数。自然线性形函数具有线性得到自然线性形函数。自然线性形函数具有线性形函数相同的性质。形函数相同的性质。1(1)2is1(1)2js2022-3-14493.4.2 一维自然二次和三次形函数一维自然二次和三次形函数将将 带入形函数的表达式得到二次自然形带入形函数的表达式得到二次自然形函数为：函数为：iXXx(1)2lx1(1)21(1)2(1)(1)ijksss 2022-3-1450 一维三次自然形函数为：一维三次自然形函数为：1(1) 31 31161(1) 31 31169(1)1 31169(1) 13116ijkmsss