文档介绍：有限元分析理论基础
有限元分析概念
有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所虚位移原理还是虚应力原理，他们所依赖的几何方程和平衡方程都是基于小变形理论的，他们不能直接应用于基于大变形理论的力学问题。
3、最小总势能法
应变能：作用在物体上的外载荷会引起物体变形，变形期间外力所做的功以弹性能的形式储存在物体中，即为应变能。
由n个单元和m个节点组成的物体的总势能为总应变能和外力所做功的差：

最小势能原理：对于一个稳定的系统，相对于平衡位置发生的位移总会使系统的总势能最小，即：
，i=1,2,3,……,n
有限元法的收敛性
有限元法是一种数值分析方法，因此应考虑收敛性问题。
有限元法的收敛性是指：当网格逐渐加密时，有限元解答的序列收敛到精确解；或者当单元尺寸固定时，每个单元的自由度数越多，有限元的解答就越趋近于精确解。
有限元的收敛条件包括如下四个方面：
1）单元内，位移函数必须连续。多项式是单值连续函数，因此选择多项式作为位移函数，在单元内的连续性能够保证。
2）在单元内，位移函数必须包括常应变项。每个单元的应变状态总可以分解为不依赖于单元内各点位置的常应变和由各点位置决定的变量应变。当单元的尺寸足够小时，单元中各点的应变趋于相等，单元的变形比较均匀，因而常应变就成为应变的主要部分。为反映单元的应变状态，单元位移函数必须包括常应变项。
3）在单元内，位移函数必须包括刚体位移项。一般情况下，单元内任一点的位移包括形变位移和刚体位移两部分。形变位移与物体形状及体积的改变相联系，因而产生应变；刚体位移只改变物体位置，不改变物体的形状和体积，即刚体位移是不产生变形的位移。空间一个物体包括三个平动位移和三个转动位移，共有六个刚体位移分量。
由于一个单元牵连在另一些单元上，其他单元发生变形时必将带动单元做刚体位移，由此可见，为模拟一个单元的真实位移，假定的单元位移函数必须包括刚体位移项。
4）位移函数在相邻单元的公共边界上必须协调。对一般单元而言，协调性是指相邻单元在公共节点处有相同的位移，而且沿单元边界也有相同的位移，也就是说，要保证不发生单元的相互脱离开裂和相互侵入重叠。要做到这一点，就要求函数在公共边界上能由公共节点的函数值唯一确定。对一般单元，协调性保证了相邻单元边界位移的连续性。
但是，在板壳的相邻单元之间，还要求位移的一阶导数连续，只有这样，才能保证结构的应变能是有界量。
总的说来，协调性是指在相邻单元的公共边界上满足连续性条件。
前三条又叫完备性条件，满足完备条件的单元叫完备单元；第四条是协调性要求，满足协调性的单元叫协调单元；否则称为非协调单元。完备性要求是收敛的必要条件，四条全部满足，构成收敛的充分必要条件。
在实际应用中，要使选择的位移函数全部满足完备性和协调性要求是比较困难的，在某些情况下可以放松对协调性的要求。
需要指出的是，有时非协调单元比与它对应的协调单元还要好，其原因在于近似解的性质。假定位移函数就相当于给单元施加了约束条件，使单元变形服从所加约束，这样的替代结构比真实结构更刚一些。但是，这种近似结构由于允许单元分离、重叠，使单元的刚度变软了，或者形成了（例如板单元在单元之间的绕度连续，而转角不连续时，刚节点变为铰接点）对于非协调单元，上述两种影响有误差相消的可能，因此利用非协调单元有时也会得到很好的结果。在工程实践中，非协调元必须通过“小片试验后”才能使用。
应力的单元平均或节点平均处理方法
最简单的处理应力结果的方法是取相邻单元或围绕节点各单元应力的平均值。

这种方法最常用于3节点三角形单元中。这种最简单而又相当实用的单元得到的应力解在单元内是常数。可以将其看作是单元内应力的平均值，或是单元形心处的应力。由于应力近似解总是在精确解上下振荡，可以取相邻单元应力的平均值作为此两个单元合成的较大四边形单元形心处的应力。
如2单元的情况下，取平均应力可以采用算术平均，
即平均应力=（单元1的应力+单元2的应力）/2。
也可以采用精确一些的面积加权平均，
即平均应力=[单元1应力× 单元1的面积+单元2应力× 单元2面积]/（单元1面积+单元2面积）
当相邻两单元面积相差不大时，两者的结果基本相同。在单元划分时应避免相邻两单元的面积相差太多，从而使求解的误差相近。
一般而言，3节点三角形单元的最佳应力点是单元的中心点，此点的应力具有1阶的精度。

首先