有限单元法：杆系结构有限元[详细]

1、有限单元法,第1章 杆系结构有限元,P-2/62,第1章 杆系结构有限元,1.1 有限单元法及其发展概况 1.2 有限元方法及解题步骤 1.3 单元分析 1.4 整体分析 1.5 按单元定位向量形成总刚度方程 1.6 约束处理及求解,P-3/62,第1章 概论,1.1 有限单元法及其发展概况 1.2 有限元方法及解题步骤 1.3 单元分析 1.4 整体分析 1.5 按单元定位向量形成总刚度方程 1.6 约束处理及求解,P-4/62,1.1.1 力学问题的解法,力学问题的求解方法可分为：解析法、数值法、半解析法 解析法（连续、精确）：问题的求解可以得到具体的数学表达式，可以计算任意位置的问题的精

2、确解答。如悬臂梁求挠度等。 数值法（离散、近似）：得不到具体表达式，只能得到某些离散点处的近似值。 半解析法：是将解析法和数值法结合起来形成的一种方法。,P-5/62,连续问题一般都是通过数学方程式的运算和求解获得精确解。但在实际工程中，可用的数学方法通常使得能用连续化处理的问题极为有限，或者使其过于简化。 相反，由于计算机的出现，使得求解离散问题比较容易。,计算力学：是根据力学中的理论，利用现代电子计算机和各种数值方法，解决力学中的实际问题的一门新兴学科。它横贯力学的各个分支，不断扩大各个领域中力学的研究和应用范围，同时也逐渐发展自己的理论和方法。 包括：有限元单元法、有限差分法、边界元法等

3、,P-6/62,常用的数值解法有下列几类： （1）有限差分法将微分方程化为差分形式，求近似解。 （2）加权残值法将微分方程化为加权积分形式，求近似解（加权残值法的五种常用做法是配点法，子域法，加辽金法，最小二乘法，矩法）。 （3）有限单元法将微分方程问题化为能量驻值问题，采用插值函数，求近似解。 （4）边界元法只在边界上进行离散。 常用的半解析法主要有：有限条法和有限元线法、薄层法等方法。,其中应用最广泛的方法就是有限单元法。几乎可以说，一切弹性力学问题都可以利用有限元法求解。,P-7/62,有限单元法(Finite Element Method)(FEM)是计算力学的重要分支，是一种将连续体

4、离散化，以求解各种力学问题的数值方法。,有限元方法的基本思想：将结构物看成由有限个人为划分的单元组成的整体，以单元结点上的解代表整体的解。它是一种化整为零、集零为整、化未知为已知的方法。不同的学科，所求解的参数不同。在计算力学中，主要有以下三种： 位移型：以结点位移为未知量。 力型：以结点力为未知量。 混合型：某些地方以结点位移为未知量，另外一些以结点力为未知量。 我们主要就“位移型”有限元进行讲解。,P-8/62,结构力学位移法,单跨超静定梁,利用已知的单跨超静定梁的内力计算公式的成果（固定求固端力，放松求单元的杆端力）,结点的平衡条件附加约束上的反力=0,建立求解未知位移的平衡方程组,P-

5、9/62,结构力学位移法,位移法：建立求解未知位移的平衡方程组,矩阵位移法，并利用计算机求解,1956年，Turner，Clough，Martin和Toop ，刚架的矩阵位移法推广用于弹性力学平面问题作近似分析,有限单元法的刚度法,P-10/62,有限单元法的刚度法,单元，结点,单元，结点,分析各单元的特性，得到物理条件（节点力与节点位移关系）,利用在结点相交的平衡条件、变形协调条件,建立求解各结点位移(ui、vi)的联立方程,P-11/62,1.1.2 有限单元法发展概况,产生构想：20世纪40年代，麦克亨利（McHenry）、雷尼柯夫（Horenikoff）、纽马克（Newmark）等首次

6、提出用框架方法求解力学问题，用简单弹性杆排列代替连续体的各个小部分，能够得到连续问题的相当好的解答。 1943年，柯兰特（Courant）第一次假设饶曲函数在一个划分的三角形单元集合体的每个单元上为简单线性函数。由于计算机尚未出现，这篇论文没有引起应有的注意。 (1945年，世界上出现了第一台电子数字计算机“ENIAC”，用于计算弹道。1956年，晶体管电子计算机。1959年，集成电路计算机。）,一、起源,P-12/62,3）1955年，德国斯图加特大学的J.H. Argyris教授发表了一组能量原理与矩阵分析的论文，奠定了有限元方法的理论基础。 4）第一个尝试：1956年，特纳（Turner

7、）、克拉夫（Clough）等将刚架分析中的位移法扩展到弹性力学平面问题，并用于飞机的结构分析和设计，系统研究了离散杆、梁、三角形的单元刚度表达式，并求得了平面应力问题的正确解答。 第一次命名：1960年，克拉夫（Clough），建立在虚位移原理或最小势能原理的基础上。在处理剖面弹性问题时，第一次提出并使用“有限元方法”的名称,P-13/62,O.C.Zienkiewize在他的“The Finite Element Method”一书中（参考书6），一开头便称： “人类大脑是有限的，以致不能一次就弄清周围许多（自然存在的和创造出的）复杂事物的特性。因此，我们先把整个系统分成特性容易了解的单个元

8、件或单元，然后由这些元件重建原来的系统以研究其特性，这是工程师、科学家甚至经济学家都采用的一种自然的方法。许多经典的数学近似方法以及工程师们用的直接近似法都属于这一范畴。” 因此，从这一意义上说确定有限单元法的准确起源时间是困难的。尽管如此, 有限单元法从应用意义上讲，它的发展始于20世纪60年代。,P-14/62,二、 三个黄金时代,第一个黄金时代： 1960年起。最重要的工作来自结构工程师，第一次解决了诸如汽车、飞机、水坝等复杂结构的力学分析。尽管当时有限元法的数学基础尚未完全建立（尽管与今天相比,当时的成就十分有限），但该方法获得了巨大成功。论文统计：1961年10篇、1965年67篇、

9、1968年303篇。,P-15/62,第二个黄金时代： 始于1965年前后，属于数值分析研究者。1963年开始出现数值分析家的论文，他们终于认识了有限元法的基本原理，事实上是逼近论、偏微分方程及变分形式和泛函分析的巧妙结合。 并得出结论：直接刚度法（即有限元法）的基础是变分原理，它是基于变分原理的一种新型里兹法（采用分区插值方案的新型里兹法）。这样就使数学界与工程界得到沟通，获得共识。从而使有限元法被公认为既有严密理论基础、又有普遍应用价值的一种数值方法。 在有限元的这一发展时期，还有必要说明两件事：一个是有限单元法与变分原理之间的关系，另一个是我国学者对有限单元法的贡献。,P-16/62,从

10、有限元法的创立过程看两者的关系：,1956年，Turner，Clough，Martin和Toop ，刚架的矩阵位移法推广用于弹性力学平面问题作近似分析,1963年的交汇融合(直接刚度法，即有限元法的基础是变分原理 )，这正是有限元创立的真实和曲折过程。通过这种曲折，正好看出两者关系的密切难分。,1943年，Courant ，变分原理和分片插值方法来求扭转问题的近似解，“变分解法”,P-17/62,各种变分原理以及其对应的有限元,不同性质的单元对应于不同的变分原理。 （1）协调位移元（采用的位移插值函数在单元间精确协调）最小势能原理。 （2）非协调位移元（采用的位移插值函数在单元间不精确协调）分

11、区势能原理。 （3）广义协调位移元（采用的位移插值函数在单元间广义协调）分区势能原理的退化形式。 （4）应力杂交元（采用应力试函数，满足平衡微分方程）最小余能原理。 （5）混合元（采用混合试函数，含位移、应力和应变）广义变分原理。 （6）分区混合元（部分单元采用位移试函数，其余单元采用应力试函数）分区混合能量原理。,P-18/62,我国学者对有限单元法的贡献,我国数学家冯康等人从1960年前后开始，也创造了系统化的有限元算法（1965年文“基于变分原理的差分格式”）编写了程序，解决了当时国防和经济建设中的一些重大课题，并奠定了数学理论基础。因此可以说，有限元法是在欧、美和中国被独立发展的。 有

12、限元及变分原理的研究领域是我国学者的研究强项。胡海昌于1954年提出的弹性力学广义变分原理为有限元法的发展提供了理论基础。冯康提出的基于变分原理的差分格式实质上就是今天的有限元法。龙驭球提出的分区和分项能量原理(1980)，分区混合有限元 (1982)，样条有限元(1984)，广义协调元(1987)和四边形面积坐标理论(1997)等，使有限元方法的分析能力和应用领域得到很大提升。,P-19/62,在学术界影响广泛的国内著作 钱伟长，变分法与有限元. 北京: 科学出版社, 1980 胡海昌，弹性力学的变分原理及其应用. 北京: 科学出版社, 1981 朱伯芳，有限单元法原理与应用. 北京: 中国

13、水利水电出版社, 第1版1979, 第2版1998）是兼备科学性和实用性的巨著。 龙驭球，有限元法概论. 北京: 高等教育出版社, 第1版, 1978, 第2版, 1991, 龙驭球，新型有限元引论. 北京: 清华大学出版社, 1992等,P-20/62,第三个黄金时代： 有限单元法的广泛应用。自1970年代开始,有限元法被迅速应用到各领域。,钢结构接头的应力分析,汽车碰撞分析,P-21/62,三、有限单元法的现状与发展趋势,1 变分原理与数值方法 针对有限元发展需要而提出的变分原理的新形式，例如分区势能、余能、混合能量变分原理，分区变分原理的退化形式及其应用，含参数变分原理及其应用，压电复合

14、材料结构变分原理，基于应变梯度理论的细观力学变分原理。 2 新型单元构造方法 对单元构造的现有模式作进一步开拓，如杂交元、混合元、拟协调元、应变元、样条元、无限元、 高精度单元（带转角元）、大单元、超（巨型）单元 ,子结构、土壤结构流体的相互作用等。 发展新的构造模式，如基于广义协调理论的广义协调元、基于分区混合变分原理的分区混合元，理性有限元，基于四边形面积坐标的四边形元，基于解析试函数的有限元等。,P-22/62,3疑难现象及破解对策 有限元学科发展中，还遗留一些疑难现象和问题，有的长期未得到破解。这些尚待破解的疑难问题自然就成为关注的焦点。例如 多种闭锁现象（剪切闭锁、薄膜闭锁、不可压缩

15、闭锁）； 网格畸变敏感现象； 有的非协调元不收敛现象； 虚假零能模式现象； 解的晃动现象； 精度损失现象（位移元的应力，层合板的层间应力）； 应力奇点现象； 数值计算病态现象。,P-23/62,4复杂深层问题 材料非线性和几何非线性有限元分析； 壳体结构屈曲稳定性分析； 塑性成型有限元分析； 撞击破坏过程数值模拟； 基于应变梯度理论的有限元法。 解的精度的研究：数值方法的误差估计、收敛性、可靠性、自适应性和优化； 自适应有限元方法 。,P-24/62,有限元的数学理论原理表明，若单元构造适当，当网格无限加密时其解应收敛到精确解，但对于一个特定网格（含单元）在缺乏精确解的情况下很难估计其精度，此

16、问题近十来年所进行的研究已取得重要进展。其方法思路是:在求得给定网格下的有限元位移解和应力解后，采用应力恢复法得到结点应力的改进值，然后利用形函数求改进估计，或称为最佳猜测应力。当按此求得的误差大于允许值时，可根据误差预测出结构各区域中单元的合理尺寸（网格尺寸）。据此，重新划分网格（由计算机自动完成），可反复多次近似于迭代过程。,P-25/62,5 耦合交叉问题 耦合问题有流固耦合、气液固耦合、土结构流体耦合、力电耦合等。 学科交叉问题有生物力学、微电子科学、材料科学中的数值模拟与优化设计等。 6与其他方法的联合沟通 例如：有限元边界元；有限元有限差分；有限元法无网格法；数值法解析法。 7.

17、向新领域的扩展 已由弹性力学平面问题扩展到空间问题和板、壳问题等整个固体力学，能对原子能反应堆、堤坝、飞机、船体、齿轮叶片等复杂结构进行应力分析。 由静力平衡问题扩展到稳定与动力问题，由弹性力学扩展到弹塑性、粘弹性、断裂力学等，由结构分析扩展到结构优化。,P-26/62,8 软件开发与CAD/CAE技术 有限元软件与有限元理论几乎是同时诞生的。只有得到工程应用，理论才会具有生命力，这就是软件的作用往往大于论文的理由。计算力学要不断吸收计算机科学技术的成果，与CAD/CAE技术共同发展。 （1） 通用程序的进一步完善 增加(强)前、后处理功能，重新包装(用户界面,多媒体)，扩充单元库、动力及非线

18、性等。 （2） 大型有限元程序微机化、网格化。 （3） 由专业研究转向普及应用（科学成果转化成生产力） （4） 向高层次发展,结构分析通用程序简介,自1980年代开始，世界各国，特别是发达国家，都花费巨大的人力和物力开发大型、通用的结构分析程序。,P-28/62,ANSYS是由美国ANSYS公司开发的大型通用有限元分析软件，ANSYS公司自1970年成立以来，不断吸取世界最先进的计算方法和计算机技术，引导世界有限元分析软件的发展，以其先进性、可靠性、开放性等特点，被全球工业界广泛认可，拥有全球最大的用户群。1995年，在分析设计类软件中第一个通过ISO9001国际质量体系认证。 ANSYS采用

19、三维实体描述法建立几何模型，几十种图素库可以模拟任意复杂的几何形状，强大的布尔运算实现模型的精雕细刻；提供多种网格划分方法，可以实现网格密度及形态的精确控制。具体划分方法有拉伸网格划分，智能自由网格划分，映射网格划分，自适应网格划分等。,ANSYS,P-29/62,SAP,SAP是Structural Analysis Program 的英文缩写，SAP程序由美国加里福尼亚大学伯克利分校的威尔逊（Wilson）、巴瑟(Bath)和彼特森（Pertson)等人在1970年代初开始研制完成的，经过二十多年的不断发展、完善，成为一个在国际上普遍受欢迎的通用结构分析软件。（80年代由北大几位老师引入我

20、国，SAP）。 SAP2000,P-30/62,ADINA,ADINA是由国际上著名的美国麻省理工学院K.J. Bathe教授领导的，ADINA R 4.3 薄板弯曲问题程序设计,P-45/62,第5章 稳定问题有限元 第6章 动力有限元法 6.1 结构振动微分方程 6.2 单元质量矩阵和单元阻尼矩阵 6.3 结构的自由振动和特征值问题 6.4 地震反应分析的计算机方法 6.5 线性地震反应时程分析法 自学： 第7章 非线性分析有限元法 第8章 钢筋砼结构非线性分析有限元法 第9章 高层建筑结构分析的有限元法,P-46/62,主要参考书,1有限元法及其应用，江见鲸，何放龙，何益斌等，机械社20

21、06版（土木工程研究生系列教材） 2有限元法概论 龙驭球, 高教社1978第一版，1991第二版 3结构矩阵分析原理 赵超燮, 人教社1982版 4有限单元法原理与应用 朱伯芳, 水电社1979版, 71.211/ZBF 5有限元法基本原理与数值方法王勖成,邵敏,清华大学社1988版 6The Finite Element Method Zeinkiewicy O C，Taylor R L,O xford: Butterworth-Heinemann, 2000,P-47/62,7The Finite Element Method (third edition) Zeinkiewicy O C

22、. 52.54/ZOC(3)( 中译本: 有限元法 科学社1985版, 52.4/JKW) 8变分法及有限元上册 钱伟长著，科学出版社1980版 9Concepts and Applications of Finite Element AnalisisCook R D， 10弹性 塑性 有限元 欧阳鬯、马文华, 湖南科技社1983版 11结构分析的有限条法 Y.K.CHEUNG, 王贻荪等译, 交通社1982版 12工程结构抗震动力学 李国豪, 上海科技社1980版, 86.213/LGH 13有限元法新论原理、程序、进展龙志飞等，中国水利水电社，2001,P-48/62,14新型有限元论龙驭

23、球，龙志飞，岑松，清华大学社2004版 15有限元法与板壳分析袁驷，崔京浩主编，清华大学社2005版 16钢筋混凝土有限元与板壳极限分析 沈聚敏等, 清华大学社1993版 17固体力学发展趋势 黄克智等, 北京理工大学社1995版 18工程结构计算机仿真分析 江见鲸,贺小岗 清华大学社1997版 19钢筋混凝土房屋结构计算机辅助设计尚守平,何放龙,刘光栋，建工社1999版 20广义协调元理论与四边形面积坐标法，龙志飞等，中国矿业大学社，2000版,P-49/62,第1章 杆系结构有限元,1.1 有限单元法及其发展概况 1.2 有限元方法及解题步骤 1.3 单元分析 1.4 整体分析 1.5 按

24、单元定位向量形成总刚度方程 1.6 约束处理及求解,P-50/62,一、有限元的力学分析方法,1.方法特点 从方法论看，有限元法是分析综合法的一种应用。先将结构分解为单元，再将单元合成结构，在一分一合中求得结构问题的解答。,结构,单元,化整为零难化易 分解 积零为整复原型 合成,P-51/62,从力学渊源看，有限元法是由刚架计算的矩阵位移法演变而来的。由刚架分析移植到弹性力学，矩阵位移法就变成了有限元法。矩阵位移法与有限元法共同的特点就是分解、合成法 。,刚架分析的 矩阵位移法,弹性力学的有限元法,移植,从数学角度看，有限元法是连续问题的一种离散化近似解法。把原来属于无限自由度的问题近似地按有

25、限自由度的问题来处理。把原来的微分方程问题转化为代数方程问题。,连续问题 无限自由度问题 微分方程问题,离散问题 有限自由度问题 代数方程问题,离散化,P-52/62,三角形面积,圆面积,P-53/62,3 有限元推导方法,从采用的力学工具分:直接法和变分法 (1)直接法：不拘泥于严格的数学形式，而侧重于易懂的物理概念。典型代表： 刚度法,P-54/62,(2)变分法：是把有限元法归结为求泛函的极值问题，使有限元建立在更加坚实的数学基础上（变分法及有限元钱伟长，科学社，1980版） 从选择基本未知量的角度分: 位移法 取结点位移作为基本未知量(displacement/stiffness method) 力 法 取结点力作作为基本未知量 (Force/Flexibility method) 混合法 同时取力和位移作为基