有限元平面问题

[]
()1
,2
1
1112
1
=∴---++==
i i i j
k i j k i i k k j j i k
k
j j i i y x N y x y x y x y x y x y x y x y x y x A
即：()()()1,,,===k k k j j j i i i y x N y x N y x N （由i,j,k 轮换性知） 同理可证：()()0,,==k k i j j i y x N y x N （作业：证明：()()k j i j i y x N j j i ,,0,=≠=）
因此
()()()()()()()()()()⎪⎪
⎩⎪⎪
⎨⎧===⎩⎨⎧≠===⇒======1
,,0,,0,,1,00,,1,,0,0
,,0,,1,k k k j j k i
i k j
i
i k k j j j j i i j k k i j j i i i i y x N y x N y x N j i j i j N y x N y x N y x N y x N y x N y x N δ （2-12）
即形函数在自己节点上为1，在其余节点上为0。

2. 在单元上任意一点，三个形函数之和为1，即
()()()1,,,=++y x N y x N y x N k j i 。

证明：
()()()()()()()[]
y c c c x b b b a a a A
y c x b a y c x b a y c x b a A
y x N y x N y x N k j i k j i k j i k k k j j j i i i k j i ++++++++=
++++++++=++21
21
,,,
⎪
⎩⎪
⎨⎧-=-=-=j i k i k j k j i y y b y y b y y b ⎪
⎩⎪
⎨⎧-=-=-=i
j k k i j j k i x x c x x c x x c ()()()()()()()()()()()()1
,,,0
02=++∴=-+-+-=++=-+-+-=++=-+-+-=++y x N y x N y x N x x x x x x c c c y y y y y y b b b A y x y x y x y x y x y x a a a k j i i j k i j k k j i j i i k k j k j i i j j i k i j k j k k j k j i （2-13）
由此可见，三个形函数中只有2个是独立的，即第三个可由其余两个表示。

3. ij 边上的形函数()k j i i N i ,,=与节点k 的坐标无关（i ，j, k 轮换），即在ij 边上有：
()()()⎪⎪⎪⎪
⎩
⎪
⎪⎪⎪⎨⎧
=--=---=0,,1,y x N x x x x y x N x x x x y x N k i
j i j
i
j i i
（i ，j, k 轮换） （2-14）
证明：设 节点i 坐标：()i i y x ,,节点j 坐标：()j j y x ,。

求：ij
边的直
线方程。

i j i i j i x x x x y y y y --=
--
()i j i
j i
i y y x x x x y y ---+=∴ ()
()()()
()()(){}()()()()()k
k j j k i i k j k i j i k j k i j i j i j j i j i j i j j j i k k i j j j j j j k j i i k
k k i j k i j c c b c b x x A x x c c b x x b A x x c c b x x b y c x b a A x b x b y c x b a A x x c b y c x b a A y c x b a A y x N N ij N y x x x c b
b y y
c x x +--=
⎥⎦⎤⎢⎣⎡
--
-=⎥⎦⎤⎢⎣⎡---+++=-+++=⎭
⎬⎫⎩⎨⎧
⎥⎦⎤⎢⎣⎡--++=++--=∴-=-=-222222,N ,:,,N y y ,j j i １
１１１
１１＝计算边上的值在得到代入直线方程为：
()()0,==++i i j i j i j j
y x N y c x b a
k
k
j j i i y x y x y x A 1112
1
=
()()()()A
y x y x y x y x y x y x x x y y x x y y c b c b j k i j k i i k k j j i k i j i i j i k j k k j 2=---++=-----=- ()k
i
j c x x y x N -=
∴, i j k x x c -=∴ ()i
j i
j x x x x y x N ij --=
∴,边上有：在
在j i 边上：
()()()(){}[]()0,2222,==++=-+++=⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡--++=++=
i i k i k k k k i k i k k k i k k i k k k k k k k y x N y c x b a A x x b y c x b a A x x c b y c x b a A y c x b a A y x N １１
１１
由性质2 ：
()()()i
j i
k j i x x x x y x N y x N y x N ---
=--=∴1,,1,
即在i ,j 边上有：
()()()⎪⎪⎪⎪
⎩
⎪
⎪⎪⎪⎨⎧
=--=---=0,,1,y x N x x x x y x N x x x x y x N k i
j i j
i
j i i
（2-15） 证毕。

同理知：（轮换）
在jk 边上有： 在ki 边上有：
()()()⎪⎪⎪⎪
⎩
⎪
⎪⎪⎪
⎨⎧--=
---==i j i k i j i j
i x x x x y x N x x x x y x N y x N ,1,0, ()()()⎪⎪⎪
⎩
⎪
⎪⎪⎨⎧---==--=i j i k j i
j i
i
x x x x y x N y x N x x x x y x N 1,0
,, 几何表示：
五、三角形单元位移函数的收敛性
（要点提示：单元位移函数的三条收敛准则及意义）
下面我们来验证所设的位移函数⎩⎨
⎧++=++=y
x v y
x u 654321αααααα 满足收敛准则（三条）。

1、 单元的位移函数解反映单元的刚体位移（包含有）
由几何方程：⎪⎪⎪⎩
⎪⎪⎪⎨⎧∂∂+∂∂=∂∂=∂∂=
x v y u y v
x u xy
y x γεε 寻找物体发生刚体位移的条件。

若物体发生刚体位移，则有：
()()x f v y f u x v y u y v
x u
xy
y
x 21,000==⇒⎪⎪⎪⎩
⎪⎪⎪⎨⎧=∂∂+∂∂==∂∂==∂∂=γεε ()的任意函数为y u 由0=xy γ得：
()()021=∂∂+∂∂x x f y y f ()()dx
x df dy y df 21=-∴
等式两侧分别为x 和y 的函数，要使其相等只有：
()()const dx
x df dy y df ===-
∴ω21
积分：()()⎩⎨
⎧+=-=x
v x f y
u y f ωω0201 式中00,v u 为积分常数
故位移：()()00201===⇒⎩
⎨⎧-==-==xy y x x v x f v y
u y f u γεεωω
即：（不难证明）以上两项是发生刚体位移的充要条件。

因为这是xy y x γεε==的情形。

故：
刚体转动的转角
方向的刚体位移方向的刚体位移---ωy x u 00v
事实上，将位移函数改变形式为：
⎪⎪⎩
⎪⎪⎨
⎧
+-+++=-++++=y
x x v y y x 635534535321222
2u αααααααααααα
显然可看出：
⎪⎪⎩⎪
⎪⎨⎧
----转的刚体转动体现绕无关）方向的刚体位移（与体现无关）方向的刚体位移（与体现z y x y y x x 2,,3
541α
ααα （其它系数意义后述）
2、
单元位移函数解反映单元的常应变
由：⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂+∂∂=∂∂=
∂∂=
x v
y u y v
x u
xy y x γεε 可以得到：
⎪⎪⎪⎩
⎪⎪⎪⎨⎧+=∂∂+∂∂==∂∂=
=∂∂=5362ααγαεαεx v
y u y v
x u
xy y x ()()53532
1
21ααααβαγ+++=+=xy
显然：⎪⎩⎪
⎨⎧-+--常量
体现了单元的剪应变为方向的常应变体现了单元沿方向的常应变体现了单元沿53
62ααααy x
由此看出，但单元的各应变均为常量。

故三角形单元在位移函数：]
⎩
⎨
⎧++=++=y x v y x u 654321αααααα 下个典的各个应变量均为常量。

故称ε为常应变单元。

3、 单元的位移函数在单元内部连续，在边界与相邻单元协调。

显然，设⎩
⎨⎧++=++=y x v y x u 654321αααααα 是单元内部的连续函数。

下面考察下边界上协调（一致）
的问题。

由形函数的第3条性质，我们证明： 对于相邻的两个单元j i e e ,,21为公共边界。

ij 边上的N:
()()()⎪⎪⎪⎪
⎩
⎪
⎪⎪⎪⎨⎧
=--=
---=0,,1,y x N x x x x y x N x x x x y x N k i
j i j i
j i i
分别写出两个单元在公共边上的位移表达
式。

对于单元1e ，其位移函数为：⎪⎩⎪
⎨⎧++=++=k
k j j i i e k k j j i i e v N v N v N v u N u N u N u 1
1
（*）
对于单元2e ，其位移函数为：⎪⎩⎪
⎨⎧++=++=m
m j j i i e m m j j i i e v N v N v N v
u N u N u N u 2
2 （**） Ij 为单元1e ，2e 的公共边界。

由形函数的性质3我们知道：
⎪⎪⎪⎪
⎩
⎪
⎪⎪⎪⎨⎧
==--=
---=0
1m
k i
j i j i
j i i N N x x x x N x
x x x N 仅与节点 i 有关。

因此，对于1e ：⎪⎪⎩
⎪⎪⎨
⎧--+⎪⎪⎭⎫ ⎝⎛---=--+⎪⎪⎭
⎫ ⎝⎛---=j i j i i i j i e j i j i i i j i e v x x x x v x x x x v u x x x x u x x x x u 1111
对于2e ：⎪⎪⎩
⎪⎪⎨⎧--+⎪⎪
⎭⎫ ⎝⎛---=--+⎪⎪⎭
⎫ ⎝⎛---=j
i j i i i j i e j i j i i i j i e v x x x x v x x x x v u x x x x u x x x x u 1122 与节点k,m 无关，仅与i, j 节点坐
标及i i v u ,有关。

j i x x ,-- 已知常数
i i v u ,-- 节点位移唯一
边界上x 唯一确定u,v
由1e
u 和2e
u 比较及1e
v 和2e
v 比较知：在公共边界上各点，ij 上位移u,v 是唯一的。

由上知：三角形单元的位移函数⎩
⎨⎧++=++=y x v y
x u 654321αααααα 满足收敛性条件。

Note: 用三角形单元计算则位移是连续的。

而应力、应变是阶梯的。

位移法（假
设位移）的结果位移要好（比应力准确）。

2-5 三角形单元的刚度矩阵（单刚）
提示：我们已经建立了三角形单元的位移函数；导出了三角形单元的形函数；并用形函数来表示其位移函数；最后，我们证明了三角形单元位移函数的收敛性。

下面我们要推导三角形单元的单元刚度矩阵。

在推导单刚前我们还有些准备工作要做。

一、 三角形单元的应变矩阵[B]
将位移函数写出来：⎪⎩⎪⎨⎧++=++=k
k j j i i k
k j j i i v N v N v N v u N u N u N u
其中：
()y c x b a A
N j j j i ++=2１
),,(k j i i = 把位移函数u,v 代入几何方程：
()()()()()()[]
⎪⎪⎪
⎪⎪⎩⎪⎪⎪⎪
⎪⎨
⎧
+++++=+++++=∂∂+∂∂=++=∂∂=++=∂∂=k k k k j j j j i i i i k
k j j i i k k j j i i xy k k j j i i y k k j j i i x v b u c v b u c v b u c A
v b v b v b u c u c u c A x v y u v c v c v c A y v u b u b u b A x u 21212121γεε
写成矩阵的形式就是：（单元上任一点的应变）
{}⎪⎪⎪⎪⎭
⎪⎪⎪⎪⎬
⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧=k k j j i i k k
j j i i k j i k j i
z y x v u v u v u b c b c b c
c c c b b b A 000
000
21εεεε （2-16） 或{}[]{}e
B ∆=ε （2-17）
式中：[]⎥⎥
⎥⎦
⎤⎢⎢⎢⎣
⎡=k k
j
j
i
i k j i k j i
b c b c b c c c c b b b A
B 00000021 （2-18） 或分块：
[][][][][]k j i B B B B = （2-19）
式（2-16）表示单元节点位移{}e
A 与单元应变{}ε的关系。

矩阵[]
B 称为应变矩阵。

式（2-18）表示应变矩阵为常数矩阵()
j k i k i i x x c y y b -=-=,，再次证明三节点三角
形单元为常应变单元。

二、 三角形单元的应力矩阵[]S
由物理方程知：()()⎪⎪
⎪
⎩
⎪
⎪
⎪⎨⎧
-∙-=+-=+-=xy xy y x y y x x u E u E u E γμτεμεσμεεσ2111122
2 ()xy xy E γμτ+=12 用矩阵表示：⎪
⎭⎪⎬⎫
⎪⎩⎪⎨⎧⎥⎥⎥⎥⎦
⎤
⎢⎢⎢⎢⎣⎡--=⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧xy y x xy y x u E γεεμμμτσσ210
0010112 （2-20）
或缩写为：{
}[]{}εσD =（2-21） 其中：[]⎥⎥⎥
⎥⎦
⎤
⎢⎢⎢
⎢⎣⎡
--=210
0010112μμμu
E D （2-22） 称为弹性矩阵（仅与弹性常数有关）。

把{}[]{}∆=B ε代入物理方程{
}[]{}εσD =，得到：{}[]{}[][]{}∆==B D D εσ 令[][][]B D S = （2-23）
则有：{
}[]{}∆=S σ （2-24） 式（2-24）表示应力与节点位移的关系。

[]S 由式（2-23）给出，称为三角形单元的应力矩阵。

显然，弹性矩阵[]D 及应变矩阵[]B 都是常量矩阵。

故应力矩阵[]S 也是一个常量矩阵。

因此三节点三角形单元的应变和应力都是常量。

三、三角形单元的单刚
建立了应力与节点位移的关系式（2-24），我们就可以推导单刚了。

我们用虚功原理来推导。

一般来说，有限元的单刚最普通的方法是用变分原理来推导。

求泛函的变分（functional 泛函的函数）。

（在力学上就是最小泛解的变分原理）。

由于我们尚未解除变分原理且对于弹
力问题，用虚功原理推导就可以了。

有人证明了用虚功原理推导和用最小泛解的变分原理来推导单刚，对于弹性力学的问题结果是一致的。

⎩⎨
⎧接受；偏于使用概念直观、清楚；容易虚功原理来推导：力学
偏于纯理论研究严密；数学基础扎实；变分原理来推导：推理
下面我们用虚功原理来推导三节点三角形单元的单刚。

1. 单刚的推导（单元是平衡的：对其应用虚功原理）
如图所示，三角形单元的节点位移和节点力为：
节点位移：
{}{
}
T
k
k
j j i i k j i e
v u v u v u =⎪⎭
⎪
⎬⎫
⎪⎩⎪⎨⎧∆∆∆=∆
节点力：
{}{
}
T
yk
xk
yj xj yi xi k j i e F F F F F F F F F F =⎪⎭
⎪
⎬⎫
⎪⎩⎪⎨⎧=
给定一组虚位移：（每个单元都可能有虚位移）（虚位移是人为假设的任意位移，其唯一的条件就是约束所允许）
{}{
}
T
k
k
j j i i k j i e
v u v u v u δδδδδδδδδδ=⎪⎭
⎪
⎬⎫⎪⎩⎪⎨⎧∆∆∆=∆
产生虚应变：
{}e
xy y x e
⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧=δγδεδεδε
则单元的外虚功为：（节点力）
{}
{}
e
e k yk k xk j yj j xj i yi i xi F v F u F v F u F v F u F T
∆=+++++δδδδδδδ
单元的内力虚功为：
(){}{}tdxdy
tdxdy e
eT
A
xy xy y y x x
A
⎰⎰
⎰⎰=++σδεδγτδεσδεσ
由虚功原理知：
{}{}
{}{}tdxdy
F eT
A
e
eT
⎰⎰
=∆σδεδ
我们设法把等式右侧的应力和虚应变换成位移和虚位移表示：
{}[]{}
εσD = {}[]{}
∆=B ε
{}[][]{}
e
e
B D ∆=∴σ
{}
[]{}
e
e
B ∆⋅=δδε
代入虚功方程右侧：
{}{}[]{}()[][]{}{}[][][]{}tdxdy B D B tdxdy B D B tdxdy e
T
eT A e T
e A e
eT A ⎰⎰
⎰⎰⎰⎰∆∆=∆∆=δδσδε{}eT ∆δ及{}e ∆都与
x,y 无关。

在有限元中当我们研究一个单元时单元内的任一点位移可由
节点位移表示，是x,y 及节点位移的函数节点位移我们认为是已知量。

故有：
{}{}{}[][][](){}e T A eT e eT tdxdy B D B F ∆∆=∆⎰⎰δδ [][][]B D B T 可能是x,y 的函数
令：[][][][]⎰⎰=
tdxdy B D B k T
A e （2-25）
则：{}[]{}e
e
e
k F ∆= （2-26）
式（2-25）为三角形单元的单刚。

式（2-26）为三角形单元的单元刚度方程。

由于[][]D B ,均为常数矩阵。

故有：[][][][]B D B t A k T
e
⋅⋅= （2-27）
书上P72已将[]e
k 各项展开。

（把k m →）大家可看一下。

2. 单刚[]e k 的物理意义
把单刚分块，则单元刚度方程可写成：
e
k j i e kk
e kj
e ki
e jk e jj e ji e ik e ij e
ii e
k j i k
k
k k k k k k k F F F ⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧∆∆∆⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧⨯66 （2-28）
展开得：{}e
k e ik e j e ij e i e ii e
i k k k F ∆+∆+∆=
{}
k jk j jj i ji e
j
k k k F ∆+∆+∆=
{}k kk j kj i ki e
k k k k F ∆+∆+∆=
显然：e ii k 表示当0,1=∆=∆=∆k j i 时在i 节点产生的节点力e i F e
ij k 表示当0,1=∆=∆=∆k i j 时在i 节点产生的节点力e i F
3. 单元[]e k 的性质（与外力无关）
1）[]e
k 是66⨯的对称矩阵，即e ji e
ij k k = （互等定理）且主元非负，且e ii k >0
2）[]e
k 是奇异矩阵：（最简单的想法，最笨的做法是证明|[]e
k |=0）
[]e k 是66⨯方阵，我们这样做： 把1）、3）、5）三个加在一起。

（见P123）
看
第
一
列
相
加
的
结
果
：
()()k j i i k j i i i k i k i j i j i i c c c c b b b b c c b b c c b b c b ++-+++=-++-++-+
=∑2
12
12121221μ
μ
μμ
⎪
⎩⎪
⎨⎧-=-=-=j
i k i k j k j i y y b y y b y y b ⎪
⎩⎪
⎨⎧-=-=-=i
j k k i j j
k i x x c x x c x x c
∑=∴=++=++∴0
0,01k j i k j i c c c b b b
同理可证：
∑
=0i
()6,5,4,3,2=i
∴[]e
k 是奇异的。

3）[]e
k 的影响因素
A. 单元的几何参数：大小Ａ（平缓过渡问题），厚度ｔ，方位（节点坐标差）()i i y x ,
B. 单元的材料特性：[]0,⇒μE
()21A A =
至此，我们已经推出了单刚，并对[]u 进行了讨论。

有了单刚后我们就可以利用平衡条件建立总刚了。

3—6 结构刚度矩阵—总刚
提示：推导出单元的刚度矩阵，就意味着我们有了单元上节点力与节点位移的关系。

与一维的问题相同，我们下一步工作就是要找到结构的总刚度矩阵。

建立以结构节点位移为未知数的结构刚度方程。

一、节点的平衡方程（内力与外力的平衡）
我们仍然用一个简单的直观的例子来推导总刚度方程，然后不失一般性的推广到一般的结构。

结构离散如图所示，取出节点3来研究节点的平衡。

首先写出各单元的单元刚度方程。

单元（1）)
1(431)
1(4443
41
3433311413
11
)
1(431⎪⎭
⎪⎬⎫
⎪⎩⎪⎨⎧∆∆∆⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=⎪
⎭
⎪⎬⎫⎪⎩⎪⎨⎧k k k k k k k k k F F F
单元（2）)
2(321)
2(3332
31
232221131211)2(321⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧∆∆∆⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧k k k k k k k k k F F F
单元（3）)
3(352)
3(3335325355522325
22)
3(352⎪⎭⎪⎬⎫
⎪⎩⎪⎨⎧∆∆∆⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧k k k k k k k k k F F F
单元（4）)
4(453)4(4445
43
545553343533)4(453⎪⎭
⎪⎬⎫⎪⎩⎪⎨⎧∆∆∆⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=⎪
⎭
⎪⎬⎫
⎪⎩⎪⎨⎧k k k k k k k k k F F F
对节点3列出平衡方程
外力：{}⎭
⎬⎫⎩⎨⎧=333y x P P P 内力：{}⎪
⎭
⎪⎬⎫
⎪⎩⎪⎨⎧=⎭⎬⎫⎩⎨⎧=∑∑e y e x y x F F F F F 33333
由平衡条件：{}{}33F P = ——平衡方程
二、总刚的形成（先写方程，再定义总刚）
结构的总刚度矩阵可由结构全部节点的平衡方程写出。

我们仅以例中的结构，第3 节点
的平衡方程说明如何建立该结构的第三个方程（子块）
按节点形式展开节点3 的平衡方程：
{}{}{}()
()()
()
)
4(434535333)
3(333535232)
2(333232131)1(5)1(35)1(3)1(33)1(1)1(31333∆+∆+∆+∆+∆+∆+∆+∆+∆+∆+∆+∆===∑k k k k k k k k k k k k F F P e
e
注意到：i e
i ∆=∆ ()5,1-=i 则（按节点重排）：
{}()()()()()5)4(35
)3(354)4(34)1(343)4(33)3(33)2(33)1(332)3(32)2(321)2(31)1(313∆++∆++∆++++∆++∆+=k k k k k k k k k k k k P
把该结构的总刚中第3个子块写出：
⎪⎪⎪⎭
⎪⎪⎪⎬⎫
⎪⎪⎪⎩⎪⎪⎪⎨⎧∆∆∆∆∆⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢
⎢⎢⎢⎢⎣
⎡=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∑∑∑∑∑-543214
,3354
,1344
1333
,2322,13154321e e
e
e
e k k k k k P P P P P
在总刚度方程中，总刚矩阵的第3行（子块）的元素为第3 个节点全部相关单元的单刚中对应下标的元素（子块）之和。

（相关节点—若i,j 同属于e 则i ,j 为相关节点；相关单元—与节点i 相连的单元为i 的相关单元）
同理，由结构其余节点的平衡方程，可以得到总刚的全部内容：
⎪⎪⎪⎭⎪⎪⎪
⎬⎫
⎪⎪⎪⎩⎪⎪⎪⎨⎧∆∆∆∆∆⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥
⎦
⎤⎢⎢⎢
⎢⎢⎢⎢⎢⎢⎢⎣
⎡=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∑∑∑∑∑∑∑∑∑∑∑∑∑-543214,355)
4(54
4
,353
)
3(52
)4(454
,144
4,143
)1(414
,3354,134
4133
3
,232
2,131253,223
3,222
21
142,11312
2
,11154321000
0k k k
k k k
k k k k k k
k k k k
k k k k k P P P P P 杆系结构的总刚中
ij k 只解来自某一单刚。

而平面问题ij k 可能来自两个单
刚。

更一般地，对于一个结构，若将其离散为m 个单元，n 节点则有：
⎪⎪⎪⎪⎪⎭
⎪
⎪
⎪
⎪⎪
⎬⎫
⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧∆∆∆∆∆⎥⎥
⎥⎥⎥
⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣
⎡=⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧n j i nn nj
ni
n n jn jj ji j j in ij
ii i i n j i n j i n j i k k k k k k k k k k k k k k k k k k k k k k k k k P P P P P
212
121212222221111121121
其22⨯为ij k 子块 [
]
T
yi xi
i p P P = ； {}{}T
i i
i v u =∆
或者写为：
{}[]{}
∆=k P
总刚形成方法：
1.对角线子块ii k 为：节点i 相关单元（ｅ）的单刚中)(e ii k 求和。

即
ii k =∑)
(e ii
k
2.非对角线子块ij k :
()
()⎪⎩⎪⎨⎧==∑0
,2,ij e ij ij k j i k k j i 为非相关节点时：若—对平面问题—个单元求和至多为相关节点时：若因此，通俗一点就是：⎩⎨⎧副子块看线主子块看点
实际上，程序是按单刚中对应下标求和进行的。

三、总刚度矩阵的性质
1. 总刚[k]为对称矩阵。

[][]
T
ji
ij
k k =（子块）
可由单刚对称性及总刚形成的方法看
ｅｘ．)2(13)1(1313k k k += )
2(31)1(3131k k k += T
k k 3113= ()2,1=e
２．总刚是一个奇异矩阵
明确的物理解释
３．总刚是一个稀疏矩阵⇒必然的有许多非相关节点 0=ij k
有条件的：带状稀疏（节点编号满足螺旋法则）
{}[]{}∆=k P
总刚度矩阵是奇异的，要进行约束处理（代入已知的边界条件） 荷载是节点载荷，不具一般性
３－７约束处理
结构总刚度方程中，[]k 是奇异的，无法直接求解{}[]{}∆=k P ，需要进行约束处理 处理方法与一维问题相同
⎪⎩
⎪
⎨⎧划零置一法置大数法消行消列法 ３－８载荷处理（载荷移置）－等效节点载荷
一、问题的提出
在有限元法中，结构离散后，单元之间的联系及单元间力的传递都是通过节点实现的。

在我们推导单元刚度方程及结构总刚方程时，我们也都隐含了一个假设，那就是：作用于单元上及结构上的载荷为节点载荷。

那么到目前为止，对于一个弹性力学的平面问题，如果作用于结构上的所有荷载都是节点载荷（集中力），我们已经可以用有限元法求解问题了。

但是，在实际的工程问题中，作用一个结构上的载荷是多种
多样的，也是比较复杂的。

我们只有对任何种类的荷载都能用节点荷载来表示，有限元法才有生命力。

对于一维杆系结构，它只是非常简单的一种情况，我们用求解单跨超静定梁的杆端力的办法，就可以进行荷载处理了。

但二维以上的问题（弹力）就不能作类似的简单处理。

例子如下（图）
二、荷载移置的原则
对于弹性力学问题而言，在处理荷载的移置时，必须满足如下的条件，才能保证计算结果达到较高的精度。

１．圣维南原理
()
()⎩
⎨
⎧
尺寸
单一尺寸远远小于结构
小范围移置
由虚功等效来保证
静力等效
２．载荷移置的唯一性（由位移函数唯一来保证）
三、载荷移置的一般公式
在进行荷载移置之前，我们首先对结构上各种的载荷进行分类载荷的种类：
WORD格式编辑整理（１）非节点的集中力（单元上某一点）（２）体积力（单元内分布）
（３）面力（单元边界上分布）
以下逐步介绍载荷移置
专业知识分享。