弹性力学有限元法.ppt

《弹性力学有限元法.ppt》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、第三章弹性力学有限元法第三章弹性力学有限元法3.1有限元法求解问题的基本步骤3.2连续体离散化3.3单元分析3.4整体分析3.5边界约束条件处理3.6求解、计算结果的整理和有限元后处理3.1有限元法求解问题的基本步骤工程结构的几何简化载荷简化边界条件的简化建立力学模型连续体离散化用合适的单元将连续体划分为有限个具有规则形状的的单元集合，单元的选取应视所分析问题的性质、规模和精度要求而定。3.1有限元法求解问题的基本步骤单元分析位移模式(位移函数、插值函数）选取单元分析单元的位移模式一般采用多项式，因为多项式计算简便，

2、并且随着项数的增加，可以逼近任何一段光滑的函数曲线。建立单元刚度矩阵建立单元节点力列阵弹性力学几何方程和物理方程静力等效原则3.1有限元法求解问题的基本步骤整体分析和有限元方程求解建立整体刚度矩阵建立整体节点力列阵代入边界条件选择适当的代数方程求解高斯消元法三角分解法波前法雅克比迭代法3.1有限元法求解问题的基本步骤结果后处理和分析应力误差的减小结果输出方式结果分析3.2连续体离散化杆状单元平面杆单元只能承受轴向的拉压载荷，平面杆单元每个节点有两个自由度。空间杆单元只能承受轴向的拉压载荷，平面杆单元每个节点有三个自由

3、度。轴力单元3.2连续体离散化梁单元平面梁单元平面梁单元每个节点有三个自由度，两个线位移一个角位移，可承受平面内的体力，集中力、分布力和垂直平面的弯矩的作用。空间梁单元空间梁单元每个节点有六个自由度，三个线位移三个角位移，可承受各个方向的体力，集中力、分布力和弯矩的作用。杆状单元3.2连续体离散化平面单元平面单元属于二维单元，只能承受单元平面内的分布力和集中力，不能承受面外载荷。三节点三角形单元每个节点有两个自由度，因此只能采用线性模式，应变值为常量，也称为常应变单元或常应力单元。3.2连续体离散化四节点四边形单元四

4、边形单元有四个节点，每个节点也有两个自由度，采用双线性位移模式，计算精度较高。平面单元3.2连续体离散化薄板弯曲单元和薄板单元薄板弯曲单元四边形弯曲单元四边形单元有四个节点，每个节点有三个自由度，主要承受横向载荷和绕水平轴的弯矩。3.2连续体离散化薄板弯曲单元三角形弯曲单元三角形单元有三个节点，每个节点有三个自由度，主要承受横向载荷和绕水平轴的弯矩。薄板弯曲单元和薄板单元3.2连续体离散化薄板弯曲单元和薄板单元薄板单元四边形薄板单元四边形薄板单元有四个节点，每个节点有五个自由度，可承受各个方向的载荷和绕水平轴的弯矩。

5、3.2连续体离散化薄板弯曲单元和薄板单元薄板单元三角形薄板单元三角形单元有三个节点，每个节点有五个自由度，可承受各个方向的载荷和绕水平轴的弯矩。3.2连续体离散化多面体单元四面体单元四面体单元有四个节点，每个节点有三个自由度。3.2连续体离散化六面体单元六面体单元有八个节点，每个节点有三个自由度。多面体单元3.2连续体离散化等参单元求解实际问题时人们总希望用最少的单元实现比较高的计算精度，而且所选用的单元对复杂结构也有比较好的适应性。计算期望：手段：单元形状的变化和单元内位移函数的变化用相同数目的结点参数和相同的插值

6、函数进行变换。四边形四节点等参元3.2连续体离散化四边形四节点等参元插值函数注意:总体坐标系适用于整体结构，局部坐标系只适用于具体某个单元。常用的对于平面问题还有八节点等参元，空间问题有八节点空间等参元，二十节点等参元等。等参单元3.2连续体离散化轴对称单元对于回转结构，如果约束条件和载荷都对称于回转轴，其应力、应变和位移也都对称于回转轴线，这类应力应变问题称为轴对称问题，通常用柱坐标来描述应力、应变和位移，单元为实心圆环体，仅截面不同三角形环单元3.2连续体离散化四边形环单元回转圆锥薄壳单元轴对称单元3.3单元分析

7、单元的插值函数（各种多项式）3.3单元分析单元的插值函数（各种多项式）四节点矩形单元的插值多项式3.3单元分析单元的插值函数（各种多项式）令：矩阵形式：3.3单元分析单元的插值函数（各种多项式）形函数特点：本点为1，他点为0在单元内任一点各形函数之和等于1单元任意一条边上的形函数，仅与该边两端节点的坐标有关，而与其他节点无关3.3单元分析单元的插值函数（各种多项式）单元形函数必须满足的条件：位移模式在单元内连续，在单元的公共边界处协调.位移模式必须包括单元的刚体位移.位移模式还必须包括单元的常应变状态.3.3单元分析

8、单元分析弹性力学平面问题的几何方程：B称为应变矩阵，其分块子矩阵为：3.3单元分析单元分析代入无量刚插值函数，应变矩阵为：令：由平面问题物理方程单元内任意一点的应力可表示为S为应力矩阵3.3单元分析单元分析D为弹性矩阵其表达式为：E是杨氏模量，是泊松比应力矩阵S的分块子矩阵为3.3单元分析单元分析代入无量刚插值函数，应力矩阵为：对于平面应变问题