第一章 绪论
1.1 引言
有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题。事实上，有限单元法已经成为在已知边界条件和初始条件下求解偏微分方程组的一般数值方法。
对于有限元的线性分析，我们已经比较熟悉，它已作为设计工具被广泛采用。但绝大多数实际问题属于非线性问题，根据产生非线性的原因，非线性问题主要有三种类型：
1． 材料非线性(简称材料非线性或物理非线性)
其特点是应力σ与应变ε之间为非线性关系，通常与加载历史有关，加载和卸载不是同一途径（图1.1），因而其物理方程
shl是谁D εσ=
中的弹性矩阵D 是应变ε的函数。但材料非线性问题仍属于小变形问题，位移和应变是微量，其几何方程是线性的。土、岩石、混凝土等具有典型的材料非线性性质，所以，土坝、岩土地基的稳定性和加固，地下洞室和边坡的稳定性都应当按材料非线性问题处理。
2．几何非线性
几何非线性属于大变形问题，位移和应变或者它们中一个是有限量。可能会有三种情况：大位移(包括线位移和角位移)、小应变，小位移、大应变和大位移、大应变。此时反映应变和位移关系的几何方程是非线性方程，例如，正应变
z ε可表示为
ΛΛ+∂∂+∂∂+∂∂+∂∂=
](()[(21222x
w x v
x u x u x ε 剪应变xy γ表示为
ΛΛ+∂∂∂∂+∂∂∂∂+∂∂∂∂+∂∂+∂∂=
y
w x w y v x v y u x u y u x v xy γ 如果应力和应变之间的关系也是非线性的，就变成了更复杂的双重非线性问题。不过，在几何非线性问题中一般都认为应力在弹性范围内，σ与ε之间呈线性关系。工程中
的实体结构和板壳结构都存在几何非线性问题，例如弹性薄壳的大挠度分析，压杆或板壳在弹性屈曲后的稳定性问题。
在采用有限元方法分析非线性问题时，以上两类都表现为结构的整体劲度矩阵K 不再是常量矩阵，而是结点位移δ的函数，还有一类问题是结点荷载R 与δ有关，这就是边界非线性问题，又称接触非线性。
3．接触非线性
由于接触体的变形和接触边界的摩擦作用，使得部分边界条件随加载过程而变化，且不可恢复。这种由边界条件的可变性和不可逆性产生的非线性问题，称为接触非线性。工程中有许多接触非线性的问题，如混凝土坝纵缝和横缝缝面的接触，面板堆石坝中钢筋混凝土面板与垫层之间的接触，岩体节理面或裂隙面的工作状态等。
目前，研究工程非线性问题比较有效的工具是非线性有限元方法。要使这一方法实用化，有两个问题必须解决。第一，由于非线性问题的数值计算工作量大大增加，需要有相当大计
算能力的工具。近10年来，高速计算机的发展已基本上满足了这一需要，同时计算费用也在继续减小。第二，非线性求解方法的精度和收敛性必须被验证。由于发展和改进了多种类型的单元，以更好
地模拟结构的工作，找到更有效的非线性求解方法，并且积累了许多经验可应用于实际工程问题，现在已经能够比较有把握地完成非线性有限元分析。非线性有限元方法正在成为一种强有力的计算工具被研究人员和工程人员所采用。
peak1.2 线弹性有限元方法的回顾
一般来说，非线性有限元法可归结为一系列线弹性问题。因而，线弹性有限元是非线性有限元法的基础。两者不仅在分析方法和求解步骤上有相似之处，而且后者要不断调用前者的结果。为此，本节将扼要地回顾一下线弹性有限元方法（以位移有限元为例）。
线弹性有限元的分析步骤大致可分： 1． 结构的离散化
(a) 轴对称圆筒 (b)取1/4结构离散
图1.2 圆筒有限元离散
在有限元分析中，首先要将连续体剖分成有限个结点上相联的单元（图1.2），单元之间的力依靠结点传递，铰接结点只能传递力，刚接结点可以传递力和力矩。根据结构的形状和计算要求，可以选择各种单元类型，图1.2采用了八结点六面体的空间等参单元。
在位移有限元中，以结点的位移δ为基本未知量。设结构离散后，总的自由度数为n ，则整体的结点位移向量为
T n ][21δδδΛΛ=δ (1-1)
单元的结点位移向量为
T d d d e w v u w v u w v u ][222111ΛΛ=δ (1-2)
d 是单元的结点数。
连续体离散后，还需要对位移边界进行约束处理，将连续的位移边界条件离散为结点的位移边界条件。例如，位移为零的边界可处理为该边界上的结点具有铰支座（图1.2）。
2． 单元分析
根据所选择的单元类型，确定单元的位移模式，将外荷载转化为等效结点荷载列阵，导出单元的应变、应力矩阵及劲度矩阵。
(1) 单元位移模式
旅游胜地英文采用广义坐标法或Seredipity 方法[1]导出位移插值函数N i 的表达式，就可以建立单元内任一点的位移f = [u v w ]T 与单元结点位移之间的关系：
e
δe N δf = (1-3)
N i 反映了单元内位移的分布形状，所以又称形函数。对于d 个结点的三维单元，
N = (1-4) ⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡d d d N N N N N N N N N 0
0000000000002
1
21
21ΛΛi N 是整体坐标（x y z ）的函数，对于图1-3
所示的等参单元，由于需要进行从实际单元（图1-3(a)）到母单元（图1-3(b)）的坐标变换，将是局部坐标（i N ξ η ζ）的函数，
坐标变换可表示为
e Nx x = (1-5)
式中
x = [x y z ]T (1-6)
为单元内任意点的整体坐标。
e x =[x 1 y 1 z 1 x 2 y 2 z 2……x d y d z d ]T (1-7)
为单元各结点的整体坐标。
（2）单元的结点荷载e
R
作用在单元上的集中力[
]
T
z
y x
P P P P =，体力[]T
Z Y
X
=p 及面力
[
T
Z Y
X
=p ]
必须转换成等效的单元结点荷载列阵e R ，
[]T
d d
d
e Z Y X Z Y X Z Y X Λ
22
2
1
1
1=R (1-8)
根据虚功原理可导出
∫∫∫∫∫++=ds dv T T T e p N p N P N R (1-9)
（3）应变和应力矩阵
英文金额大写对于小变形问题，应变和位移之间满足以下的几何方程
Lf ε= (1-10)
式中
swiftcode是什么意思(1-11)
[
T zx yz
xy
z
y x
γγγεεε=ε]⎥⎥⎥⎥⎥
⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎢⎣
⎡∂∂∂∂∂∂∂∂∂∂∂∂
∂∂∂∂∂∂
=x y
z
z x y z y x 0
0000
000L (1-12) 为微分算子矩阵。将式（1-3）代入（1-10）可得
e B δε= (1-13)
其中
LN B = (1-14)
为单元的应变矩阵。
对于线弹性体，应力和应变遵循以下的物理方程
D εσ= (1-15)
D 为弹性矩阵，可表示为
⎥⎥
⎥⎥
⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+++=G G G G Sym G G 000000000000222λλλλλλD (1-16)
λ和G 为拉密系数。将式（1-13）代入（1-15）可得
e S δσ= (1-17)oursky
其中
DB S = (1-18)
称单元的应力矩阵。
（4）单元劲度矩阵
在有限元方法中，各单元之间的相互作用力表现为结点力，设单元的结点力为
[]T
d d
d e W V U W V U W V U Λ
222111=F (1-19)
则虚功方程可表示为
∫∗∗=V
T e T e dv σεF δ)()( (1-20)
式中，是单元的结点虚位移，e
∗δ
∗
ε是虚应变。将式（1-13）和（1-15）代入式（1-20）
，可得出出单元结点力与结点位移的关系
e e k δF = (1-21)
式中
∫=V
T dv DB B k (1-22)
为单元劲度矩阵。
3． 整体分析
弹性力学问题在按位移求解时，与基本微分方程对应的泛函是弹性体的势能p π，它是位移的函数。与微分方程等价的变分方程就是最小势能原理，即
0=p δπ (1-23)
在有限元方法中，弹性体的总势能为各单元势能的集合，单元势能由单元的形变势能和外力势能组成
e
p πe U e
V e
e e p V U +=π (1-24)
式中
∫
=V
T
e dv U σε2
1 (1-25) ∫∫−−−=σ
foreign是什么意思S T V
T T e ds dv V p f p f P f (1-26)
将式（1-13）和（1-15）代入（1-25），式（l-3）代入（l-26），并注意式（1-22）和（1-9），可得
e T e e T e e p R δk δδ)()(2
1
−=
π (1-27) 引入单元选择矩阵，使
e上海市回国留学人员服务中心
C δC δe e =
则
R δK δδT T p −=2
1
π (1-28)
根据式（1-23），有
0=∂∂δ
p π
由（1-28）式导出位移有限元的支配方程
R K δ= (1-29)
其中
∑=e
e T e kC C K )( (1-30)
为整体劲度矩阵。
∑=e
e T e R C R )( (1-31)
是整体结点荷载列阵。整体分析的目的就是将k 集合成K ，e
R 集合成R 。
4． 解方程，求位移和应力
求解支配方程（1-29），得到结点位移，然后根据式（1-18）和（1-15）求出单元的应变和应力。
δ
1.3 学习内容及课程安排
值此之际本门课程主要为土木、水利类专业的本科生和研究生开设，根据土木、水利工程的特点，限于学时，课程内容主要研究材料非线性问题的有限元方法。与其它一些以理论为主的有限元教材不同，为了使本教材更实用，重点放在非线性有限元方法的计算机实践上。当然，作为一门完整的课程，对材料的本构关系及非线性问题的求解也要作必要的介绍。
工程中出现的材料非线性来自多方面的原因，例如，可能来自材料的弹塑性或某种形式的超弹性；也可能由于时间因素引起，如材料的粘性或动态瞬时特性。无论是一维、二维还是三维结构都可能存在以上类型的非线性，限于篇幅，本教材只能选其中几种进行讨论，并以弹塑性为主。|
第一章回顾了弹性有限元方法，并简单介绍了几类非线性问题。
第二章对非线性方程组的求解方法进行研究，列出一些基本的求解方法，并讨论了非线性计算的迭代收敛准则。
第三章主要研究材料的本构关系，包括非线性弹性、弹塑性、粘弹性、粘塑性等，其中以弹塑性本构关系为主。
第四章对材料非线性有限元方法进行了介绍，重点研究了迭代求解的过程。本章还列出几种在非线性有限元分析中常遇到的特殊单元。
职场英语培训班
第五章给出一个二维的非线性有限元程序，对程序结构、主要变量和主要的子程序作了说明。
