ANSYS第二讲专用课件.pptVIP

第二讲 有限元基本原理 有限单元法的基本思路 结构的离散化 用有限元法对结构进行应力分析时，首先要将结构进行离散化。即将一个连续体看成由有限个单元组成的体系。弹性力学平面问题中最常见的单元是三角形单元。 所有作用在单元上的载荷都按静力等效的原则移置到结点上，并在受几何约束的结点处设置相应的铰支座。这样就得到了用以代替原来弹性体的有限单元计算模型。 位移模式 取一个典型的三角形单元进行力学分析。在有限单元位移法中，假设结点上的位移是基本未知量。为了能用单元的结点位移表示单元中的应变和应力分量，必须假定一个位移模式，也就是说根据单元的结点位移去构造单元上的位移插值函数。 位移插值函数 采用线性插值，即假定单元上的位移分量是坐标的线性函数： 它们可以由结点位移确定如下： 联立求解上述方程，可得： 其中： 而： 是三角形ijm的面积。 于是可以得到： 其中： 同理得： 可以将位移模式改写为矩阵模式： 有了单元的位移模式，就可以借助平面问题的几何和物理方程，导出用单于的结点位移表示单元中的应变和应力分量的公式。 由： 得到： 或简写为： 将应变代入物理方程： 可得： 即为用单元中的结点位移表示单元中应力的关系式。 式中[D]为弹性矩阵，对于平面应力问题，矩阵为： 单元的总势能 我们已经知道由各个单元的位移模式就形成了整个结构的位移模式。按弹性力学最小势能原理，结构中最接近于真实解的位移应该是使结构总势能取得最小值的那组位移函数。 由于在位移函数公式中，结点位移为自变量，这样就使一个泛函的极值问题变为一个多元函数的极值问题。为此我们来讨论单元的总势能关于结点位移的表达式。 每一个单元的总势能由该单元的应变能以及此单元上所有外力的势能组成。 单元的应变能 平面应力状态下，设物体厚度为h，则单元中的应变能为： 将{ε}和[Bi]代入上式，应用矩阵相乘的转置的逆序法则,注意到弹性矩阵[D]的对称性,有： 因为矩阵[B]及[D]的元素都是常量，所以可记： 从而单元的应变能可写为： 利用{ε}=[B]{δ}e，有： 注意到[B]=[Bi Bj Bm]，记子矩阵 单元上体积力的势能 物体中常见的体力为旋转离心体力和重力。在平面问题中，体积力在z轴方向的分力为零，设单元体积中的体积力为： 单元上体积力具有的势能为： 单元上表面力的势能 设物体边界上一单元某边上受到表面力的作用，单位长度上所受到的表面力为： 则单元上表面力的势能为： 单元节点上集中力的势能 如果弹性物体受到集中力{R}e 的作用，通常划分单元网格时都在集中力的作用点设置结点。设某单元3个结点上所受到的集中力为： 于是该单元上集中力的势能是： 单元中的总势能 综合前面的几种情况，可以得到单元中的总势能为： 单元中的总势能 分别引进单元体积力，表面力，集中力向量如下： 单元中的总势能 则单元中的总势能可以表示为： 物体中的总势能 把各单元的总势能叠加起来，就可得到整个弹性体的总势能。为了便于叠加和归并，需将单元刚度矩阵表达式作适当的改写。 假设结构离散化后共有n个结点，将编号为 l的结点位移记为： 则结构的结点位移向量： 是一个2n维的列向量。 物体中的总势能（续） 可将单元刚度矩阵式用补零的办法由6X6的矩阵扩大到2nX2n的矩阵 物体中的总势能（续） 如果在物体上划分的单元总数是e0,再引进 结构的总刚度阵： 物体总势能就可写为： 物体中的总势能（续） 代入约束条件后的弹性体总势能可以写为： 空间问题的有限单元法 用有限单元法求解弹性力学空间问题，首先也要将连续的空间物体用一系列的单元离散化。 空间问题中，最简单的是四面体单元。离散的空间结构是这些单元只在节点处以空间铰相互连接的集合体。 空间问题的有限单元法（续） 位移模式 空间问题中，每一个结点有3个位移分量，单元结点位移向量由12个分量组成，分别表示为： 位移模式（续） 假定单元内的位移分量为坐标的线性函数: 位移模式（续） 将上式中的第一式应用于4个结点，则有： 位移模式（续） 由上式可解出a1,a2,a3和a4再代回位移分量的表达式，可得： 式中： 为形函数，其中： 位移模式（续） 位移模式（续） 用同样的方法，可以得到： 合并υ,ν,ω的表达式，可以将单元内任一点的位移写为： 单元中的应变和应力 在空间问题中，每点有6个应变分量，由几何关系： 单元中的应变和应力(续) 单元中的应变和应力(续) 可以看出，应变矩阵[B]中的元素都是