有限元计算原理与方法

《有限元计算原理与方法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、1.有限元计算原理与方法有限元是将一个连续体结构离散成有限个单元体，这些单元体在节点处相互铰结，把荷载简化到节点上，计算在外荷载作用下各节点的位移，进而计算各单元的应力和应变。用离散体的解答近似代替原连续体解答，当单元划分得足够密时，它与真实解是接近的。1.1.有限元分析的基本理论有限元单元法的基本过程如下：1.1.1.连续体的离散化首先从几何上将分析的工程结构对象离散化为一系列有限个单元组成，相邻单元之间利用单元的节点相互连接而成为一个整体。单元可采用各种类型，对于三维有限元分析，可采用四面体单元、五西体单元和六面体单元等。在Plaxis3DFoundation程序中，土体

2、和桩体主要采用包含6个高斯点的15节点二次楔形体单元，该单元由水平面为6节点的三角形单元和竖直面为四边形8节点组成的，其局部坐标下的节点和应力点分布见图3.1，图3.115节点楔形体单元节点和应力点分布界面单元采用包含9个高斯点的8个成对节点四边形单元。在可能出现应力集中或应力梯度较大的地方，应适当将单元划分得密集些；若连续体只在有限个点上被约束，则应把约束点也取为节点：若有面约束，则应把面约束简化到节点上去，以便对单元组合体施加位移边界条件，进行约束处理；若连续介质体受有集中力和分布荷载，除把集中力作用点取为节点外，应把分布荷载等效地移置到有关节点上去。最后，还应建立一个适

3、合所有单元的总体坐标系。由此看来，有限单元法中的结构已不是原有的物体或结构物，而是同样材料的由众多单元以一定方式连接成的离散物体。因此，用有限元法计算获得的结果只是近似的，单元划分越细且又合理，计算结果精度就越高。与位移不同，应力和应变是在Gauss积分点(或应力点)而不是在节点上计算的，而桩的内力则可通过对桩截面进行积分褥到。1.1.1.单元位移插值函数的选取在有限元法中，将连续体划分成许多单元，取每个单元的若干节点的位移作为未知量，即，单元体内任一点的位移为。引入位移函数N(x,y,z)表示场变量在单元内的分布形态和变化规律，以便用场变量在节点上的值来描述单元内任一点的场

4、变量。因此在单元内建立的位移模式为：（3-1）其中：，I为单位矩阵。按等参元的特性，局部坐标到整体坐标的坐标转换也采用与位移插值类似的表达式。经过坐标变化后子单元与母单元(局部坐标下的规则单元)之间建立一种映射关系。不管内部单元或边界附近的单元均可选择相同的位移函数，则为它们建立单元特性矩阵的方法是相同的。因此，对于15节点楔形体单元体内各点位移在整体坐标系下一般取：上式中的为整体坐标系下节点处的位移值，为在局部坐标系下节点相应的形函数。1.1.1.单元特性分析利用几何方程、本构方程、虚功原理或位能变分方程求解单元节点力与节点位移关系的表达式，即单元刚度矩阵。根据几何方程可建

5、立单元内的应变矩阵：（3-3）其中，对于小变形线性弹性问题，根据物理方程建立单元内的应力矩阵：（3-5）其中,为几何矩阵，为弹性矩阵，为应力矩阵，。根据虚功原理求出单元中的节点力：（3-6）其中为单元的劲度矩阵，对于整体结构上的任一点i，建立平衡方程：为i节点上的外荷。上式表示与围绕i点的各单元在i点上的节点力之和相平衡。1.1.1.总体特性分析对每一个位移未知的节点，都可写出3-7式的方程，利用结构力的平衡条件和边界条件把各个单元按原来的结构重新联接起来，形成分析对象的整体有限元平衡方程组：（3-8）其中，为整体劲度矩阵，；为整个结构的节点位移矩阵，为整个结构的节点荷载矩阵

6、，是已知的。由式(3-8)求出节点位移,由式(3-3)、式(3-5)求出各单元的应变和应力。1.2.非线性有限元分析非线性现象是在实际的结构分析中经常遇到的问题。与线性分析相比，非线性分析中荷载与位移之间的关系已不是直线关系，而是曲线关系。土体的非线性分析一般来说采用非线性的分析方法，选用适当的土体本构系，进行有限元计算。非线性问题一般有材料非线性和几何非线性两种。几何非线性即存在大变形，其变化的几何形状可能引起结构的非线性响应，即应变与位移的关系不里线性，应变不仅包括位移对坐标的一阶导数，还要包括高阶导数。在进行小应变或者小变形分析时，假定位移和变形总是足够小(这种假定取决

7、于特定分析要求中的精度等级)可以忽略结构变形对系统刚度的影响，即基于最初几何形状的结构刚度的一次迭代足以计算出分析结果。随着变形位移增长，一个有限单元的已移动的坐标可以多种方式改变结构的刚度，进行多次迭代来获得一个有效的解，这就是几何非线性。除了结构大变形引起剐度变化以外。许多与材料有关的参数同样可以改变结构刚度。材料的非线性即是材料的应力—应交关系是非线性的。主要有弹性非线性模型和弹塑性模型两大类。弹性非线性理论是以弹性理论为基础，在微小的荷载增量范围内，把土看作弹性材料，从一个荷载增量变化到另一个荷