有限元理论与方法.docVIP

第一章绪论

有限元发展过程:

有限元法在西方起源于收音机和导弹的结构设计，发表这方面文章最早而且最有影响的是西德J.H.Argyrb教授，于1954—1955年间分阶段在《AircraftEngineering》上发表上许多有关这方面的论文，并在此基础上写成了《能量原理与结构分析》，此书内容提供了有限元法的理论基础。美国的M.T.Turner、R.W.cloagh、H.C.martin和L.J.Topp等人于1956年发表了了篇题为《复杂结构的刚度和挠度分析》一文，此文提出了计算复杂结构刚度影响系数的方法，并说明了如何利用计算机进行分析。美国于1960年在一篇介绍平面应力分析的论文中，首先提出了有限元的名字。1965年英国及其合作者解决了将有限元法应用于所有场的问题，使有限元法的应用更加广泛。

线性有限元

有限元几何非线性

非线性有限元

材料非线性

有限元法的基本思路:

有限元法的基本思路和基本原理以结构力学中的位移法为基础，把复杂的结构或连续体看成为有限个单元的组合，各单元彼此在节点处连续而组成整体，把连续体分成有限个单元和节点，称之为离散化，

先对单元进行特性分析，然后根据各单元在节点处的平衡协调条件建立方程，综合后作整体分析。

这样一分一合，先离散再综合的过程，就把复杂结构或连续体的计算问题转化为简单单元的分析与综合问题。

有限元分析中可采取三种方法:

位移法——取节点位移作为基本未知数

力法——取节点力作为基本未知数

混合法——

有限元法分析过程:

1、结构离散化(单元划分)

2、选择位移模式

为了能用节点位移表示单元体的位移、应变和应力，在分析连续体时，必须对单元中位移的分布做出一定的假定，也就是假定位移是坐标的某种简单函数，这种函数称为位移模式或位移函数(形函数)。

e(1)uN,,,,,,,,

3、分析单元的力学特性

(1)利用几何方程:由位移表达式导出用点位移表示单元应变的关系式

e,,,,为单元内任一点的应变列阵(2),,,,,,,,,

(2)利用物理方程，由应变的表达式导出用节点位移表示单元应力的关系式

e,,,,,,DD(3),,,,,,,,,,,,

是单元内任一点的应力列阵是材料的弹性矩阵,D,,,,

(3)利用虚功原理建立作用于单元上的节点力和节点位移之间的关系式，即单元的刚度方程(平衡方程)

eeKR,,,,,,,,

4、计算等效节点力

弹性体经过离散化后，假定力是通过节点从一个单元传递到另一个单元，但是作为实际的连续体，力是从单元的公共边界传递到另一个单元的，因而，这种作用在单元边界上的表面力、体积力、集中力等都需要等效移置到节点上去，所用方法虚功等效。

5、组装总刚度阵，建立结构的平衡方程

有两方面内容:?组装总刚?组装总的载荷列阵

得到:KR,,,,,,,,

6、求解结点的位移和计算单元应力

第二章有限元法的理论基础—加权余量

法和变分原理

本章要点

?微分方程的等效积分形式及其“弱”形式的实质和构造方法，任意函数和场

函数应满足的条件。

?不同形式的加权余量法中权函数的形式和近似解的求解步骤，以及伽辽金

(Galerkin)方法的特点。

?线性自伴随微分方程变分原理的构造方法和泛函数的性质，以及自然边界条

件和强制边界条件的区别。

?经典里兹(Ritz)方法的求解步骤、收敛性及其局限性。

?两种形式的虚功原理(虚位移原理和虚应力原理)的实质和构造方法。?从虚功原理导出最小位能原理和最小余能原理的途径和各自的性质，以及场

函数事先应满足的条件。

1.1引言

在工程和科技领域内，对于许多力学问题和物理问题，人们可以给出它们的数学模型，即应遵循的基本方程(常微分方程和偏微分方程)和相应的定解条件。但能用解析的方法求出精确解的只是少数方程性质比较简单，且几何形状相当规则的情况。对于大多数问题，由于方程的非线性性质，或由于求解域的几何形状比较复杂，则只能采用数值方法求解。20世纪60年代以来，随着电子计算机的出现，特别是最近20年来软、硬件技术的飞

速发展和广泛应用，数值分析方法已成为求解科学技术问题功能强大的有力工具。

已经发展的偏微分方程数值分析方法可以分为两大类。一类是以有限差分法为代表，其特点是直接求解基本方程和相应定解条件的近似解。一个问题的有限差分法的求解步骤归纳为:首先将求解域划分为网格，然后在网格的节点上用差分方程来近似微分方程。当采用较密的网格，即较多的节点时，近似解的精度可以得到改进，借助于有限差分法，能够求解相当复杂的问题，特别是求解方程建立于固结在空间坐标(欧拉(Euler)坐标系)的流体力学问题，有限差分法有自身的优势。因此在流体力学领域内，至今仍占支