Abaqus帮助文档翻译（一）快速介绍有限元

1.介绍
没什么好讲的，大部分情况是在吹自己有多厉害，只翻译了关于有限元的知识。
1.1有限元方法的快速回顾
任何有限元分析的第一步都是通过一个有限元的集合来描述一个实体几何的结构。每个有限元表示了物质结构的离散部分。有限单元由共享节点连接，节点和有限单元的集合称为网格（mesh）。每单位长度、面积或网格内的单元数称为网格密度。在应力分析中，节点位移（the displacements of the nodes）是Abaqus计算的基本变量。一旦节点位移已知，每个有限元中的应力和应变就可以很容易地确定。

1.1.1使用隐式的方法获得结点位移
我们从一个简单的桁架开始，如图1所示，一个桁架的一端受约束，一端受加载。

分析的目的是求出桁架自由端位移、桁架内部应力和桁架约束端反力。在这种情况下，如图1所示的杆将用两个桁架单元建模。在Abaqus桁架中，单元只能承受轴向载荷。图2显示了离散化的模型以及节点和单元标签。

模型中每个节点的自由体图如图3所示。一般来说，每个节点将承受一个施加于模型的外荷载P和一个由附加在该节点上的单元的应力引起的内部荷载I。对于一个处于静态平衡的模型，作用在每个节点上的合力必须为零;也就是说，每个节点的内部和外部负载必须相互平衡。对于节点a，则平衡方程为:

假设杆的长度变化很小，单元1中的应变为：

其中，分别为节点a，b的位移，L为单元的原始长度。
假设材料是线弹性的，杆中的应力由应变乘以杨氏模量E给出：

作用在端节点上的轴向力等于杆内的应力乘以杆的截面积a，由此得到内力、材料性能和位移之间的关系:

.
节点a保持平衡，可以得到：

.
节点b处的平衡必须考虑在该节点连接的两个单元所作用的内力。单元1的内力，现在作用方向相反，所以是负的。得到节点b的平衡方程：

.
节点c的平衡方程：

.
对于隐式方法，需要同时求解平衡方程才能得到所有节点的位移。将平衡方程写成矩阵的形式。如果两个单元的性质和维数相同，则平衡方程可简化为：

.
通常，每个单元的单元刚度项EA/L都不同，因此我们把这个模型中两个单元的单元刚度写成K1与K2。我们需要求得是外力P与内力I平衡的平衡方程的解。当考虑收敛性和非线性（convergence and nonlinearity）时，我们把它写成：

.
因此，对于完整的二元三节点结构，我们修改符号，将平衡方程改写为：

.
在隐式方法中，这个方程组可以求解得到三个未知变量的值$u^b$