《有限元基本理论》PPT课件.ppt

1、ANSYS软件及其工程应用,四川大学水利水电学院 费文平,主要内容,第1章 有限元基本理论 第2章 ANSYS功能简介 第3章 ANSYS基本过程 第4章 ANSYS入门与准备 第5章 模型输入及修复 第6章 坐标系 第7章 选择、组件与部件 第8章 实体建模技术 第9章 布尔操作 第10章 单元属性,第11章 网格划分 第12章 加载求解技术 第13章 后处理技术 第14章 结构非线性分析 第15章 模态分析 第16章 耦合和约束方程 第17章 APDL基础 第18章 子模型 第19章 热分析 第20章 热-应力耦合分析,第一章 有限元基本理论,平衡方程,几何方程,物理方程,边界条件,物理系

2、统,有限元离散,单元的位移场,(假定单元内位移函数),单元节点关系,求解区域的位移场、应力场,简单化,1.1 有限元分析 (FEA),有限元分析 是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。它利用简单而又相互作用的元素，即单元，用有限数量的未知量去逼近无限未知量的真实系统。,定义,1.2 有限单元法的基本思想,将连续的结构离散成有限个单元，并在每一单元中设定有限个节点，将连续体看作只在节点处相连接的一组单元的集合体。 选定场函数的节点值作为基本未知量，并在每一单元中假设一近似插值函数，以表示单元中场函数的分布规律。 利用力学中的某种变分原理去建立用以求节点未知量的有限单元法方

3、程，将一个连续域中有限自由度问题化为离散域中有限自由度问题。,1.3 物理系统举例,几何体 载荷 物理系统,1.3.1 平衡方程,1.3.2 几何方程,1.3.3 物理方程(本构方程),拉梅系数,体积应变,剪切模量,1.3.4 边界条件,应力边界条件,位移边界条件,1.4 有限元模型,真实系统,有限元模型,有限元模型 是真实系统理想化的数学抽象。,定义,1.5 自由度(DOFs),自由度(DOFs) 用于描述一个物理场的响应特性。,结构 DOFs,ROTZ,UY,ROTY,UX,ROTX,UZ,定义,1.6 节点和单元,节点:空间中的坐标位置，具有一定自由度和存在相互物理作用。,单元: 一组节

4、点自由度间相互作用的数值、矩阵描述（称为刚度或系数矩阵)。单元有线、面或实体以及二维或三维的单元等种类。,有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。,载荷,定义,1.6 节点和单元 (续),信息是通过单元之间的公共节点传递的。,.,.,.,A,B,.,.,.,.,.,.,.,.,A,B,.,.,.,2 nodes,1.6 节点和单元 (续),节点自由度是随连接该节点 单元类型 变化的。,J,I,I,J,J,K,L,I,L,K,I,P,O,M,N,K,J,I,L,三维杆单元 (铰接) UX, UY, UZ,三维梁单元 UX, UY, UZ, ROTX, ROTY,

5、 ROTZ,二维或轴对称实体单元 UX, UY,三维四边形壳单元 UX, UY, UZ, ROTX, ROTY, ROTZ,三维实体热单元 TEMP,J,P,O,M,N,K,J,I,L,三维实体结构单元 UX, UY, UZ,1.7 单元形函数,FEA仅仅求解节点处的DOF值。 单元形函数是一种数学函数，规定了从节点DOF值到单元内所有点处DOF值的计算方法。 因此，单元形函数提供出一种描述单元内部结果的“形状”。 单元形函数描述的是给定单元的一种假定的特性。 单元形函数与真实工作特性吻合好坏程度直接影响求解精度。,1.7 单元形函数(续),1.7 单元形函数(续),DOF值可以精确或不太精确

6、地等于在节点处的真实解，但单元内的平均值与实际情况吻合得很好。 这些平均意义上的典型解是从单元DOFs推导出来的（如：结构应力、热梯度）。,1.7 单元形函数(续),如果单元形函数不能精确描述单元内部的DOFs，就不能很好地得到导出数据，因为这些导出数据是通过单元形函数推导出来的。 当选择了某种单元类型时，也就十分确定地选择并接受该种单元类型所假定的单元形函数。 在选定单元类型并随之确定了形函数的情况下，必须确保分析时有足够数量的单元和节点来精确描述所要求解的问题。,1.8 直杆受自重作用的拉伸问题,1.8 直杆受自重作用的拉伸问题(续),就整个直杆来说，位移函数U(x)是未知的，但对每一单元

7、可以近似地假设一位移函数，它在结点上等于结点位移。此处，假设单元中的位移按线性分布 ，即：,1.8 直杆受自重作用的拉伸问题(续),有了位移插值函数，就可以按材料力学公式求出应变和应力用节点位移表示的公式：,1.8 直杆受自重作用的拉伸问题(续),外载荷与结点的平衡方程,为第i个结点上承受的外载荷,1.8 直杆受自重作用的拉伸问题(续),假定将直杆分割成3个单元，每个单元长为a=L/3，则对结点2，3，4列出的平衡方程为：,1.8 直杆受自重作用的拉伸问题(续),1.8 直杆受自重作用的拉伸问题(续),联立求解线性代数方程组得：,1.9 有限单元法解题的一般步骤,结构的离散化 选择位移模式 建

8、立平衡方程 求解节点位移 计算单元中的应力和应变,1.9.1 结构的离散化,将分析的结构物分割成有限个单元体，使相邻的单元体仅在节点处相连接，而以如此单元的结合体去代替原来的结构。,1.9.2 选择位移模式(形函数),首先对单元假设一个位移差值函数，或称之为位移模式，得到用节点位移表示单元体内任一点的唯一的关系式 有了位移模式，就可利用几何关系和应力-应变关系表出用单元节点位移表示单元中应变和应力的表达式,1.9.3 三角形单元的形函数,基本假定：假定单元内的位移可以用一个比较简单的函数来表示，如线性插值函数。这在单元划分比较密的情况下是合理可行的。,1.9.3 三角形单元的形函数(续),将三

9、角形单元的3个顶点的2个方向位移代入位移函数可求出6个待定系数。即可用节点的位移表示内部任意一点的位移：,1.9.4 建立平衡方程,可利用最小势能原理建立结构的节点载荷和节点位移之间的关系式，即结构的平衡方程,1.9.5 求解结点位移,将边界条件代入线性代数方程组 后，经解算可求得所有未知的结点位移。,1.9.6 计算单元中的应变和应力,依据求得的结点位移，由 可求得单元中任一点的应变和应力。,平面问题的有限单元法,结构的离散化,用有限元法对结构进行应力分析时，首先要将结构进行离散化。即将一个连续体看成由有限个单元组成的体系。弹性力学平面问题中最常见的单元是三角形单元。 所有作用在单元上的载荷

10、都按静力等效的原则移置到结点上，并在受几何约束的结点处设置相应的铰支座。这样就得到了用以代替原来弹性体的有限单元计算模型。,位移模式,取一个典型的三角形单元进行力学分析。在有限单元位移法中，假设结点上的位移是基本未知量。为了能用单元的结点位移表示单元中的应变和应力分量，必须假定一个位移模式，也就是说根据单元的结点位移去构造单元上的位移插值函数。,位移模式（续）,位移插值函数,采用线性插值，即假定单元上的位移分量是坐标的线性函数： 它们可以由结点位移确定如下：,位移模式(续),联立求解上述方程，可得：,位移模式(续),其中： 而： 是三角形ijm的面积。,位移模式(续),于是可以得到： 其中：

11、同理得：,位移模式(续),可以将位移模式改写为矩阵模式：,单元中的应变和应力,有了单元的位移模式，就可以借助平面问题的几何和物理方程，导出用单于的结点位移表示单元中的应变和应力分量的公式。 由：,单元中的应变和应力(续),得到： 或简写为：,单元中的应变和应力(续),将应变代入物理方程： 可得： 即为用单元中的结点位移表示单元中应力的关系式。,单元中的应变和应力(续),式中D为弹性矩阵，对于平面应力问题，矩阵为：,单元的总势能,我们已经知道由各个单元的位移模式就形成了整个结构的位移模式。按弹性力学最小势能原理，结构中最接近于真实解的位移应该是使结构总势能取得最小值的那组位移函数。 由于在位移函

12、数公式中，结点位移为自变量，这样就使一个泛函的极值问题变为一个多元函数的极值问题。为此我们来讨论单元的总势能关于结点位移的表达式。 每一个单元的总势能由该单元的应变能以及此单元上所有外力的势能组成。,单元的应变能,平面应力状态下，设物体厚度为h，则单元中的应变能为：,单元的应变能（续）,将和Bi代入上式，应用矩阵相乘的转置的逆序法则,注意到弹性矩阵D的对称性,有：,单元的应变能（续）,因为矩阵B及D的元素都是常量，所以可记：,单元的应变能（续）,从而单元的应变能可写为： 利用=Be，有：,单元的应变能（续）,注意到B=Bi Bj Bm，记子矩阵,单元上体积力的势能,物体中常见的体力为旋转离心体

13、力和重力。在平面问题中，体积力在z轴方向的分力为零，设单元体积中的体积力为： 单元上体积力具有的势能为：,单元上表面力的势能,设物体边界上一单元某边上受到表面力的作用，单位长度上所受到的表面力为： 则单元上表面力的势能为：,单元节点上集中力的势能,如果弹性物体受到集中力Re 的作用，通常划分单元网格时都在集中力的作用点设置结点。设某单元3个结点上所受到的集中力为： 于是该单元上集中力的势能是：,单元中的总势能,综合前面的几种情况，可以得到单元中的总势能为：,单元中的总势能,分别引进单元体积力，表面力，集中力向量如下：,单元中的总势能,则单元中的总势能可以表示为：,物体中的总势能,把各单元的总势

14、能叠加起来，就可得到整个弹性体的总势能。为了便于叠加和归并，需将单元刚度矩阵表达式(2-18)作适当的改写。 假设结构离散化后共有n个结点，将编号为 l的结点位移记为： 则结构的结点位移向量： 是一个2n维的列向量。,物体中的总势能（续）,可将单元刚度矩阵式用补零的办法由6X6的矩阵扩大到2nX2n的矩阵,物体中的总势能（续）,如果在物体上划分的单元总数是e0,再引进 结构的总刚度阵： 物体总势能就可写为：,物体中的总势能（续）,代入约束条件后的弹性体总势能可以写为：,空间问题的有限单元法,空间问题的有限单元法,用有限单元法求解弹性力学空间问题，首先也要将连续的空间物体用一系列的单元离散化。

15、空间问题中，最简单的是四面体单元。离散的空间结构是这些单元只在节点处以空间铰相互连接的集合体。,空间问题的有限单元法（续）,位移模式,空间问题中，每一个结点有3个位移分量，单元结点位移向量由12个分量组成，分别表示为：,位移模式（续）,假定单元内的位移分量为坐标的线性函数:,位移模式（续）,将上式中的第一式应用于4个结点，则有：,位移模式（续）,由上式可解出a1,a2,a3和a4再代回位移分量的表达式，可得： 式中： 为形函数，其中：,位移模式（续）,位移模式（续）,用同样的方法，可以得到： 合并,的表达式，可以将单元内任一点的位移写为：,单元中的应变和应力,在空间问题中，每点有6个应变分量，由几何关系：,将,的表达式代入上式，得到: 式中：,单元中的应变和应力(续),单元中的应变和应力(续),可以看出，应变矩阵B中的元素都是常量，从而单元中的应变都是常量，故线性位移模式的四面体单元是常应变单元。 由应力-应变关系，得到单元中的应力为