有限元分析思想

,单击此处编辑母版标题样式,*,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,有限元分析,Finite Element Analysis,课程目标,了解什么是有限单元法、有限单元法的基本思想。,学习有限单元法的原理，主要结合弹性力学问题来介绍有限单元法的基本方法，包括单元分析、整体分析、载荷与约束处理、等参单元等概念。,初步学会使用商用有限元软件分析简单工程问题。,了解有限元软件的基本结构和有限单元法当前的进展情况。,课程评估,作业,25%,随堂测验,25%,专题分析与报告,50%,课程进度安排,第,16,周，课堂教学,第,710,周，学习,ANSYS,与上机练习,第,1215,周，课堂教学,第,16,周，课程总结,/,交流,参考书目,王勖成,邵敏编著,.,有限单元法基本原理和数值方法,.,北京,:,清华大学出版社,1997,朱伯芳著,.,有限单元法原理与应用（第,2,版）,.,北京,:,中国水利水电出版社,1998,Saeed Moaveni.Finite Element Method Theory and Application with ANSYS.New Jersey:Prentice-Hall,Inc.,1999,O.C.Zienkiewicz,R.L.Taylor.The finite element method(5th ed).Oxford;Boston:Butterworth-Heinemann,2000,郭和德编,.,有限单元法概论，清华大学，,1998,1,有限单元法简介,1.1,有限单元法的形成,1.2,有限单元法的基本思路,1.3,有限单元法的计算步骤,1.4,有限单元法的进展与应用,1.1,有限单元法的形成,两类典型的工程问题,第一类问题，可以归结为有限个已知单元体的组合。,例如，材料力学中的连续梁、建筑结构框架和桁架结构。,平面桁架结构，由,6,个承受轴向力的“杆单元”组成。,1889,年建成的,Effiel,塔，由,18036,个部件组成,大型编钟“中华和钟”的振动分析及优化设计，,由曾攀教授完成。,第二类问题，通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。例如弹性力学问题，热传导问题，电磁场问题等。,热传导问题的控制方程与换热边界条件如下：,300MW,汽轮机低压转子淬火,两类问题的对比,第一类问题的研究对象称为离散系统。离散系统是可解的，但是求解复杂的离散系统，要依靠计算机技术。,第二类问题的研究对象称为连续系统。可以建立描述连续系统的基本方程和边界条件，通常只能得到少数问题的解析解。对于许多实际的工程问题，需要用近似算法求解。,有限单元法形成的背景,在寻找近似解法的过程中，工程师和数学家从两个不同的路线得到了相同的结果，即有限单元法（,Finite Element Method),。,有限单元法的形成可以回顾到二十世纪,50,年代，它的形成直接得益于土木结构分析中的矩阵位移法和在飞机结构分析中所获得的成果。,1954-1955,年，在航空工程杂志上发表了一组能量原理和结构分析论文。,1956,年，在纽约举行的航空学会年会上介绍了一种新的计算方法，将矩阵位移法推广到求解平面应力问题。,1960,年，,R.W.Clough,在他的名为“,The finite element in plane stress analysis”,的论文中首次提出了有限元（,Finite Element,）这一术语,从固体力学的角度来看，桁架结构与分割成有限个分区后的连续体在结构上存在相似性。,有限单元法的数学基础（,1,）,数学家们则发展了微分方程的近似解法，包括有限差分方法，变分原理和加权余量法。,在,1963,年前后，经过,J.F.Besseling,R.J.Melosh,R.E.Jones,R.H.Gallaher,T.H.H.Pian,（卞学磺）等许多人的工作，认识到有限单元法就是变分原理中,Ritz,近似法的一种变形，发展了用各种不同变分原理导出的有限元计算公式。,有限单元法的数学基础（,2,）,1965,年和（张佑启）发现只要能写成变分形式的所有场问题，都可以用与固体力学有限单元法的相同步骤求解。,1969,年和指出可以用加权余量法特别是,Galerkin,法，导出标准的有限元过程来求解非结构问题。,我国学者的贡献,陈伯屏（结构矩阵方法）,钱令希（余能原理）,钱伟长（广义变分原理）,胡海昌（广义变分原理）,冯康（有限单元法理论）,20,世纪,60,年代初期，冯康等人在大型水坝应力计算的基础上，独立于西方创造了有限元方法并最早奠定其理论基础。,-,数学辞海,第四卷,1.2,有限单元法的基本思路,将连续系统分割成有限个分区或单元,用标准方法对每个单元提出一个近似解,将所有单元按标准方法组合成一个与原有系统近似的系统。,自重作用下等截面直杆的解,受自重作用的等截面直杆如图所示，杆的长度为,L,，截面积为,A,，弹性模量为,E,，单位长度的重量为,q,，杆的内力为,N,。,试求：杆的位移分布，杆的应变和应力。,自重作用下等截面直杆的材料力学解答,自重作用下等截面直杆的有限单元法解答,1,）离散化,如图所示，将直杆划分成,n,个有限段，有限段之间通过一个铰接点连接。,称两段之间的连接点为结点，称每个有限段为单元。,第,i,个单元的长度为,L,i,，包含第,i,，,i+1,个结点。,2,）用单元节点位移表示单元内部位移,第,i,个单元中的位移用所包含的结点位移来表示,。,第,i,结点的位移,第,i,结点的坐标,第,i,个单元的应变,应力,内力,3,）把外载荷集中到节点上,把第,i,单元和第,i+1,单元重量的一半，集中到第,i+1,结点上,4,）建立结点的力平衡方程,对于第,i+1,结点，由力的平衡方程可得,：,令,根据约束条件，,对于第,n+1,个结点，第,n,个单元的内力与,第,n+1,个结点上的外载荷平衡，,建立所有结点的力平衡方程，再加上约束条件可以得到由,n+1,个方程构成的方程组，可解出,n+1,个结点的位移。,1.3,有限单元法的计算步骤,有限单元法的计算步骤归纳为以下三个基本步骤：,网格划分（离散化）,单元分析,整体分析,网格划分,对弹性体进行必要的简化，再将弹性体划分为有限个单元组成的离散体。,单元之间通过单元节点相连接。,由单元、结点、结点连线构成的集合称为网格。,网格划分,通常把三维实体划分成四面体（,Tetrahedron,）或六面体（,Hexahedron,）单元的网格,四面体,4,结点单元,六面体,8,结点单元,网格划分,三维实体的四面体单元划分,三维实体的六面体单元划分,网格划分,通常把平面问题划分成三角形（,triangular),或四边形,(quadrilateral),单元的网格,三角形,3,节点单元,四边形,4,节点单元,网格划分,平面问题的三角形单元划分,平面问题的四边形单元划分,单元分析,弹性力学问题的单元分析，就是建立各个单元的节点位移和节点力之间的关系式。,由于将单元的节点位移作为基本变量，单元分析首先要为单元内部的位移确定一个近似表达式，然后计算单元的应变、应力，再建立单元中节点力与节点位移的关系式。,单元分析,单元有三个结点,I,、,J,、,M,，每个结点有两个位移,u,、,v,和两个结点力,U,、,V,。,单元分析,单元的所有结点位移、结点力，可以表示为结点位移、结点力向量（,vector,）,单元的结点位移和结点力之间的关系用张量（,tensor,）来表示,整体分析,对由各个单元组成的整体进行分析，建立节点外载荷与结点位移的关系，以解出结点位移，这个过程为整体分析。,在边界结点,i,上受到集中力作用。结点,i,是三个单元的结合点，因此要把这三个单元在同一结点上的结点力汇集在一起建立平衡方程。,整体分析,i,结点的结点力：,i,结点的平衡方程：,1.4,有限单元法的进展与应用,1.4.1,主要研究领域与有限元软件,从二十世纪,60,年代中期以来，进行了大量的理论研究，不但拓展了有限单元法的应用领域，还开发了许多通用或专用的有限元分析软件。,主要的研究领域有：,计算方法：大型线性方程组的解法，非线性问题的解法，动力问题计算方法。,高精度单元,多物理场耦合与复杂材料模型,目前使用较多的通用有限元软件如下表所列,软件名称,简介,MSC/Nastran,著名结构分析程序，最初由,NASA,研制，由,MSC,维护开发,MSC/Dytran,动力学分析程序,MSC/Marc,非线性分析软件,ANSYS,通用结构分析软件,ADINA,非线性分析软件,ABAQUS,非线性分析软件,一些专用的有限元软件如下表所列,软件名称,简介,Deform,金属体积成形分析,Autoform,金属板料成形分析,DYNAFORM,金属板料成形分析,SysWeld,焊接与热处理分析,MSC,中国,ANSYS,中国,http,：,ABAQUS,ADINA,DEFORM,AUTOFORM,SYSWELD,飞箭公司的有限元程序自动生成系统,FEPG,应用实例,有限单元法已经成功地应用在以下一些领域：,固体力学,包括强度、稳定性、震动和瞬态问题的分析；,传热学,电磁场,流体力学,转向机构支架的强度分析（刘道勇，东风汽车工程研究院动，用,MSC/Nastran,完成）,金属体积成形过程的分析（用,Deform,软件完成）,型材挤压成形模拟,T,形锻件的成形模拟,焊接残余应力分析（用,Sysweld,完成）,焊接过程的温度分布与轴向残余应力,结构与焊缝布置,复杂形状工件的组织转变预测（用,NSHT3D,完成）,淬火,3.06 min,时的温度、组织分布,