1、2001年10月1日,La-1,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有限元分析 (FEA) 方法,Definition,2001年10月1日,La-2,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有限元分析 (FEA),有限元分析是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。,定义,历史典故 结构分析的有限元方法是由一批学术界和工业界的研究者在二十世纪五十年代到二十世纪六十年代创立的。 有限元分析理论已有100多年的历史，是悬索桥
2、和蒸汽锅炉进行手算评核的基础。,2001年10月1日,La-3,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有限单元法（续）,将连续的结构离散成有限个单元，并在每一单元中设定有限个节点，将连续体看作只在节点处相连接的一组单元的集合体。 选定场函数的节点值作为基本未知量，并在每一单元中假设一近似差值函数已表示单元中场函数的分布规律。 利用力学中的某种变分原理去建立用以求节点未知量的有限单元法方程，将一个连续域中有限自由度问题化为离散域中有限自由度问题。,2001年10月1日,La-4,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),物理系统举例,几
3、何体 载荷 物理系统,2001年10月1日,La-5,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有限元模型,真实系统,有限元模型,有限元模型 是真实系统理想化的数学抽象。,定义,2001年10月1日,La-6,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),自由度（DOFs）,自由度(DOFs) 用于描述一个物理场的响应特性。,结构 DOFs,ROTZ,UY,ROTY,UX,ROTX,UZ,2001年10月1日,La-7,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),节点和单元,节点: 空间中的坐标位置，具有一定自由度和
4、 存在相互物理作用。,单元: 一组节点自由度间相互作用的数值、矩阵描述（称为刚度或系数矩阵)。单元有线、面或实体以及二维或三维的单元等种类。,有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。,载荷,2001年10月1日,La-8,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),节点和单元 (续),每个单元的特性是通过一些线性方程式来描述的。作为一个整体，单元形成了整体结构的数学模型。尽管梯子的有限元模型低于100个方程（即“自由度”），然而在今天一个小的 ANSYS分析就可能有5000个未知量，矩阵可能有25，000，000个刚度系数。,历史典故
5、 早期 ANSYS是随计算机硬件而发展壮大的。ANSYS最早是在1970年发布的，运行在价格为1，000，000的CDC、由Univac和IBM生产的计算机上，它们的处理能力远远落后于今天的PC机。一台奔腾PC机在几分钟内可求解50005000的矩阵系统，而过去则需要几天时间。,2001年10月1日,La-9,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),节点和单元 (续),信息是通过单元之间的公共节点传递的。,.,.,.,A,B,.,.,.,.,.,.,.,.,A,B,.,.,.,2 nodes,2001年10月1日,La-10,ANSYS培训教程 版本 5.5 XJ
6、TU MSSV(001128),节点和单元 (续),节点自由度是随连接该节点 单元类型 变化的。,J,I,I,J,J,K,L,I,L,K,I,P,O,M,N,K,J,I,L,三维杆单元 (铰接),UX, UY, UZ,三维梁单元,二维或轴对称实体单元,UX, UY,三维四边形壳单元,UX, UY, UZ,三维实体热单元,TEMP,J,P,O,M,N,K,J,I,L,三维实体结构单元,ROTX, ROTY, ROTZ,ROTX, ROTY, ROTZ,UX, UY, UZ,UX, UY, UZ,2001年10月1日,La-11,ANSYS培训教程 版本 5.5 XJTU MSSV(001128)
7、,单元形函数,FEA仅仅求解节点处的DOF值。单元形函数是一种数学函数，规定了从节点DOF值到单元内所有点处DOF值的计算方法。因此，单元形函数提供出一种描述单元内部结果的“形状”。单元形函数描述的是给定单元的一种假定的特性。 单元形函数与真实工作特性吻合好坏程度直接影响求解精度。,2001年10月1日,La-12,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),单元形函数(续),2001年10月1日,La-13,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),单元形函数(续),遵循:DOF值可以精确或不太精确地等于在节点处的真实解，但单元内的平均值
8、与实际情况吻合得很好。这些平均意义上的典型解是从单元DOFs推导出来的（如，结构应力，热梯度）。如果单元形函数不能精确描述单元内部的DOFs，就不能很好地得到导出数据，因为这些导出数据是通过单元形函数推导出来的。,2001年10月1日,La-14,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),单元形函数(续),遵循原则:当选择了某种单元类型时，也就十分确定地选择并接受该种单元类型所假定的单元形函数。 在选定单元类型并随之确定了形函数的情况下，必须确保分析时有足够数量的单元和节点来精确描述所要求解的问题。,2001年10月1日,La-15,ANSYS培训教程 版本 5.5
9、 XJTU MSSV(001128),直杆受自重作用的简单拉伸问题,2001年10月1日,La-16,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),就整个直杆来说，位移函数U(x)是未知的，但对每一单元可以近似地假设一位移函数，它在结点上等于结点位移。此处，假设单元中的位移按线性分布 ，即：,2001年10月1日,La-17,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有了位移插值函数，就可以按材料力学公式求出应变和应力用节点位移表示的公式：,2001年10月1日,La-18,ANSYS培训教程 版本 5.5 XJTU MSSV(001128)
10、,Definition,外载荷与结点的平衡方程,为第i个结点上承受的外载荷,2001年10月1日,La-19,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),假定将直杆分割成3个单元，每个单元长为a=L/3，则对结点2，3，4列出的平衡方程为：,2001年10月1日,La-20,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),2001年10月1日,La-21,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),解线性代数方程组,联立求解方程组得：,Definition,2001年10月1日,La-22,ANSYS培训教程 版本 5.
11、5 XJTU MSSV(001128),有限单元法解题的一般步骤,Lesson Objectives,2001年10月1日,La-23,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),结构的离散化 选择位移模式 建立平衡方程 求解节点位移 计算单元中的应力和应变,2001年10月1日,La-24,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),结构的离散化,将分析的结构物分割成有限个单元体，使相邻的单元体仅在节点处相连接，而以如此单元的结合体去代替原来的结构。,Definition,2001年10月1日,La-25,ANSYS培训教程 版本 5.5
12、XJTU MSSV(001128),选择位移模式,首先对单元假设一个位移差值函数，或称之为位移模式，得到用节点位移表示单元体内任一点的唯一的关系式,Definition,2001年10月1日,La-26,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),有了位移模式，就可利用几何关系和应力-应变关系表出用单元节点位移表示单元中应变和应力的表达式,选择位移模式（续）,2001年10月1日,La-27,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),建立平衡方程,可利用最小势能原理建立结构的节点载荷和节点位移之间的关系式，即结构的平衡方程,Definition,2001年10月1日,La-28,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),求解结点位移,将线性代数方程组 代入边界条件后，经解算可求得所有未知的结点位移。,Definition,2001年10月1日,La-29,ANSYS培训教程 版本 5.5 XJTU MSSV(001128),计算单元中的应变和应力,依据求得的结点位移，由可求得单元中任一点的应变和应力。,Definition,
