有限元 1-概论有元 1-概论.doc

《有限单元法》讲义 《有限单元法》讲义 第1章 概论 计算力学是根据力学中的理论，利用现代电子计算机和各种数值方法，解决力学中的实际问题的一门新兴学科。它横贯力学的各个分支，不断扩大各个领域中力学的研究和应用范围，同时也逐渐发展自己的理论和方法。 有限单元法是计算力学的重要分支，是一种将连续体离散化，以求解各种力学问题的数值方法。在大学里学习过四大力学（理、材、结、弹）。众所周知，结构力学的研究对象是杆件结构(如框架、桁架)；弹性力学的主要研究对象是板、壳、实体结构(如平板、水坝)。 在结构力学的位移法中，是将计算对象转换成若干已知的单跨超静定梁(如图1-1)。 图 1-1 位移法计算简图 在有限单元法的刚度法中,是将计算对象划分成有限个容易求解的单元(“Elements”，又称元素“Components”或元件)。如图1-2所示深梁,可划分成若干个矩形单元、三角形单元,单元与单元之间用结点联结,常以结点的位移分量作为基本未知量。根据结构受力特性的不同，确定相应的位移分量，如图示深梁，则取每个结点的u,v位移为基本未知量。而图示平板则每个结点取3个位移分量（） 图 1-2 有限元法的单元划分 在结构力学位移法中,是利用已知的单跨超静定梁的内力计算公式的成果(固定求固端力,放松求单元的杆端力),由结点的平衡条件─附加约束上的反力＝0,建立求解未知位移的平衡方程组）。 在有限元刚度法中,先分析各元素的自身特性，然后利用在结点相交的平衡条件、变形协调条件、物理条件,建立求解各结点位移(ui、vi)的联立方程。 在工程领域,许多情况下都是利用有限个已知的元件得到一个合理的模型,以此求解原来的复杂的问题,这种求解问题的方法就叫作“离散”。 利用“离散(discrete)”方法求得的解答是 近似解 反之,若问题只有利用无穷小(无限个元件)才能定义,便称它为连续问题， 利用“连续(continuous)”方法求得的解答是 精确解。 连续问题一般都是通过数学方程式的运算和求解获得精确解。但在实际工程中，可用的数学方法通常使得能用连续化处理的问题极为有限，或者使其过于简化。 相反，由于计算机的出现，使得求解离散问题比较容易，利用离散法求解复杂问题最有效的方法就是有限单元法。几乎可以说，一切弹性力学问题都可以利用有限元法求解。 1.1 有限单元法发展概况 一、起源 有限单元法(The Finite Element Method)一词是R.W.Clough 于1960年首次提出(R.W.Clough, The Finite Element in plane stress an analysis,proc.2nd ASCE conf. on Electronic Computation, pittsburgh,pa., sept. 1960). 但有的学者认为，有限元法的基本思想在1941年就已产生(《有限元法概论》,邹仲康等译，湖南科技社，1988年版)。 1941年, Horenikoff提出所谓 网格法──将平面弹性体看成一批链杆和梁。 1943年 报告了 Courant著作中 使用的一组三角形单元，提出了单元的概念。1943年出版的Courant 1941年作的报告“Variational methods for the solution of problems of equilibrium and vibration”。他用变分原理和分片插值方法来求扭转问题的近似解。他在论文题目上把后来称为有限元的这种解法归结为“变分解法”。由于计算机尚未出现，这篇论文没有引起应有的注意。 其次，工程技术界第一篇有限元论文是1956年Turner，Clough，Martin和Toop的论文“Stiffness and deflection analysis of complex structures”。他们把刚架的矩阵位移法推广用于弹性力学平面问题作近似分析，在题目上把这种解法归结为复杂结构的刚度法（或直接刚度法），随后Clough定名为有限元法。 而O.C.Zienkiewize在他的“The Finite Element Method”一书中（参考书5），一开头便称： The limitation of the human mind are such that it can not grasp the behaviour of its complex surroundings and creations in one operation. Thus the process of subd