有限元计算原理与方法..

1.有限元计算原理与方法

有限元是将一个连续体结构离散成有限个单元体，这些单元体在节点处相互铰结，把荷载简化到节点上，计算在外荷载作用下各节点的位移，进而计算各单元的应力和应变。用离散体的解答近似代替原连续体解答，当单元划分得足够密时，它与真实解是接近的。

1.1.有限元分析的基本理论

有限元单元法的基本过程如下：

1.1.1.连续体的离散化

首先从几何上将分析的工程结构对象离散化为一系列有限个单元组成，相邻单元之间利用单元的节点相互连接而成为一

个整体。单元可采用各种类型，对于三

维有限元分析，可采用四面

体单元、五西体单元和六面体

单元等。在Plaxis 3D Foundation

程序中，土体和桩体主要采用包含6个高

斯点的15节点二次楔形体单元，该单元

由水平面为6 节点的三角形单元和竖直面

为四边形8节点组成的，其局部坐标下

的节点和应力点分布见图3.1，界面单元

采用包含9个高斯点的8个成对节点四边

形单元

图3.1 15 节点楔形体单元节点和应力点分布在可能出现应力集中或应力梯度较大

的地方，应适当将单元划分得密集些; 若连续体只在有限个点上被约束，则应把约束点也取为节点：若有面约束，则应把面约束简化到节点上去，以便对单元组合体施加位移边界条件，进行约束处理；若连续介质体受有集中力和分布荷载，除把集中力作用点取为节点外，应把分布荷载等效地移置到有关节点上去。

最后，还应建立一个适合所有单元的总体坐标系。

由此看来，有限单元法中的结构已不是原有的物体或结构物，而是同样材料

的由众多单元以一定方式连接成的离散物体。因此，用有限元法计算获得的结果

只是近似的，单元划分越细且又合理，计算结果精度就越高。与位移不同，应力和应变是在Gauss积分点（或应力点）而不是在节点上计算的，而桩的内力则可通过对桩截面进行积分褥到。

1.1.

2.单元位移插值函数的选取

在有限元法中，将连续体划分成许多单元，取每个单元的若干节点的位移作为未知量，即{ }e =[U i，V i，W i，...]T，单元体内任一点的位移为{f}二[u,v,w]T。

引入位移函数N（x,y,3表示场变量在单元内的分布形态和变化规律，以便用场变量在节点上的值来描述单元内任一点的场变量。因此在单元内建立的位移模

式为：

{f}=[N]{}e（3-1）其中：[N^[IN1,IN2,IN3??…1N15]，I 为单位矩阵。

按等参元的特性，局部坐标（，，）到整体坐标（x,y,z）的坐标转换也采用与位移插值类似的表达式。经过坐标变化后子单元与母单元（局部坐标下的规则

单元）之间建立一种映射关系。不管内部单元或边界附近的单元均可选择相同的位移函数，则为它们建立单元特性矩阵的方法是相同的。因此，对于15节点楔形体单元体内各点位移在整体坐标系（x, y, z）下一般取：

15

u 二人

N

i (

, ,

)u

i

i =1 15

v 八 N i (, i 丄 15

w 八 N i (

, , )W i

i J

上式中的（U i ，V i ，wJ 为整体坐标系下节点i 处的位移值，N i 「，「）为在局部 坐标系下节点相应的形函数。

1?1?3?单元特性分析

利用几何方程、本构方程、虚功原理或位能变分方程求解单元节点力与节 点位移关系的表达式，即单元刚度矩阵。

(3-3)

其中[B1=[B ，B ??…BJ ,

对于小变形线性弹性问题，根据物理方程建立单元内的应力矩阵:

「S'

(3-2)

根据几何方程可建立单元内的应变矩阵｛；｝=｛；

z z y Y Y ｝:

x y Z xy yZ zx -

1岔M /点x

0 0 1 0

刑/ Ey 0 0 0 cN i / cz

科/￡y

EN i /&x

0 cN i / cz cN i /&y

| 5Ni / & 0 cNi / ex J

(3-4)

｛；「｝珂D]｛ ；｝珂D][B]｛、｝e

(3-5)

其中,[B]为几何矩阵，[D]为弹性矩阵，[S]为应力矩阵, [S]=[D][B]

o

｛F ｝e

. 根据虚功原理求出单元中的节点力

｛F ｝_[k]

(3-6)

其中

[k]

为单元的劲度矩阵，

[k]

=

[B]T[D][B]dxdxdz

e

[B i ]

对于整体结构上的任一点i ,建立平衡方程:

' {F}二{R}

e

｛R ｝

为i 节点上的外荷。上式表示

｛R ｝

与围绕i 点的各单元在i 点上的节 点力之和相

平衡。

对每一个位移未知的节点，都可写出 3-7式的方程，利用结构力的平 衡条件和边界条件把各个单元按原来的结构重新联接起来，形成分析对象 的整体有限元平衡方程组:

(3-8)

其中，［k ］为整体劲度矩阵，K ij = k j .｛ ｝为整个结构的节点位移矩 阵，｛ R ｝为整个结构的节点荷载矩阵，是已知的。由式 (3-8)求出节点位

移｛「? ｝ 由式(3-3)、式(3-5)求出各单元的应变和应力。

1.2. 非线性有限元分析

非线性现象是在实际的结构分析中经常遇到的问题。与线性分析相 比，非线性分析中荷载与位移之间的关系已不是直线关系，

而是曲线关系

土体的非线性分析一般来说采用非线性的分析方法，选用适当的土体本构 系，进行有限元计算。

非线性问题一般有材料非线性和几何非线性两种。

几何非线性即存在大变形，其变化的几何形状可能引起结构的非线性

(3-7)

1.1.4.

总体特性分析

响应，即应变与位移的关系不里线性，应变不仅包括位移对坐标的一阶导数，还要包括高阶导数。在进行小应变或者小变形分析时，假定位移和变形总是足够小（这种假定取决于特定分析要求中的精度等级）可以忽略结构变形对系统刚度的影响，即基于最初几何形状的结构刚度的一次迭代足以计算出分析结果。随着变形位移增长，一个有限单元的已移动的坐标可以多种方式改变结构的刚度，进行多次迭代来获得一个有效的解，这就是几何非线性。

除了结构大变形引起剐度变化以外。许多与材料有关的参数同样可以

改变结构刚度。材料的非线性即是材料的应力一应交关系是非线性的。主

要有弹性非线性模型和弹塑性模型两大类。弹性非线性理论是以弹性理论为基础，在微小的荷载增量范围内，把土看作弹性材料，从一个荷载增量变化到另一个荷载增量，土体的弹性常数发生变化，以考虑非线性；弹塑性模型理论认为土体的变形包括弹性和塑性变形两部分，把弹性理论和塑性理论结合起来建立的本构模型。土体中的弹塑性本构关系都是用增量形

式表示的，因此，计算方法也宜用增量法。某级荷载增量PR］作用下,

各单元的应力状态不同。有些可能处于弹性区，则刚度矩阵要用弹性矩阵

［D］，有些可能产生塑性屈服，则须运用屈服准则、硬化规律和流动法则

建立的弹塑性刚度矩阵［D ep］来代替［D］。反映到式（3-5），其中的矩阵

［D ］不是常量其随应力或应变改变，由此推导的劲度矩阵［K］也随应力或变形而变。对于相适应流动法则g = f，贝y：

Cf歼T

[D]{—}{—}T[D] 式中A为塑性硬化模量，

是硬化参数函数。因此，不管是材料非线性还是几

何非线性，推出的劲度矩阵将随位移而变。

[D epP[D]- (3-8)

A {—}T[D]{—}

因此，不管是材料非线性还是几何非线性，推出的劲度矩阵将随位移而变。

[KC ）]{「} ={R} （3-10）这是位移的非线性方程组。直接解这样的方程组是困难的，因此简化为一系列