wi (P, Pi )
(P, Pi )Rd 0

R(Pi ) 0
（ i 1,2, , n ）
配点法又叫点匹配法。
（2）子域法
将求解区域划分成 n 个子域，
每次选取权函数在一个子域上为 1， 其他子域上为零。
wi 01（（PP不在在子子域域i内i内））
即权函数为
wi

R ci
（ i 1,2, , n ）
（4）伽辽金法
选取权函数序列与基函数序列相同。
wi ui
ui Rd 0

m
ui L u jd ui fd

j 1

（ i 1,2, , n ）
在上述几种加权余量法中， 伽辽金法应用最广泛。 有限元法基于伽辽金法
➢计算系数阵
Ki, j Ni L(N j ) d bi Ni f d
这些积分可以分单元进行。例如对右
图所示的局部编码，K01、K00以及b0 的计算公式为：
K00
N L(N )d 1 2 3 4 5 6 0
0
K01 16 N0L(N1)d
7.3 二维泊松方程的有限元法
Ki, j Ni L(N j ) d bi Ni f d
➢ 以二维静电场泊松方程的求解 为例。
Ku b
2u 2u L(u) f
x2 y2 u g

目标：依据加权余量法，利用分域基，建立离散的代数 方程组，即确定系数{Kij} 和{bi}。
Rd 0
i
（ i 1,2, , n ）
（3）最小二乘法
按使方程余量平方积分最小选取权函数。
令 I (c1, c, , cm ) R2d

使 I 最小的条件为
I 0 ci
（ i 1,2, , n ）
得
R Rd 0
ci
u(x, y) a bx cy
代入三个顶点的坐标和函数值， 可以解出a、b、c。得到
u(x,
y)

u1
1 ( x,
y)

u2
2
(x,
y)

u3
3 ( x,
y)
111
其中，
1 2
x1
x2
x3
y1 y2 y3
11 1
1

1 2
x
x2
x3
x1
x
x3
y1 y y3
单元节点的编号按 逆时针方向排列！
1 11
3

1 2
x1
x2
x
y1 y2 y
u(x,
y)

u1
1 ( x,
y)

u2
2 (x,
y)

u3
3 ( x,
y)
记住我们的任务 —寻找基函数
u (x, y) 1N1 2 N2 3N3
对比
u(x,
要保持对 称性；有 更简便的 做法
f 为已知函数。
为求 u ，设有一组完备、线性无关的函数 u1, u2 , , uk , ， 取其前 m 项的线性组合作为 u 的近似解 u 。 若当 m 时，有 u u ， 则称 u1, u2 , , uk , 为基函数序列， uk 为基函数。 c1, c2 , , cm 为待定系数。
Ni是连续的，从而通过Ni构造的逼近函数也是连续的。
➢计算系数阵 Ki, j Ni L(N j ) d bi Ni f d
在积分 Kij NiL(N j )d 中，对于确定的 i，j的有效取值为i
本身以及与节点i相联的周围节点，积分的有效区域为以i、 j 为公共节点的所有三角形单元 ，在这些单元中Ni、Nj才有 交叠。
7 二维泊松方程的有限元法
有限元法 可以从变分原理导出， 也可以从加权余量法导出。 前者需要补充泛函、变分法、欧拉方程、 泛函极值等数学知识，推导过程比较复杂。 后者相对比较直观，而且应用范围更广， 推导过程简单。
7.1 加权余量法
1、加权余量概念
假定边值问题方程
Lu f
式中， u 为未知函数， L 是算符（算子），表示对 u 的一种运算，

M
M O M M M


Kmi L Kmj L Kmm um bm
L
O M M
[K
(e)
]


K (e) ii
K (e) ji

K (e) mi
K (e) ij
K (e) jj
K (e) mj
K (e) im
➢ 场域离散
二维问题常使用三角形单元离散，便于处理复杂的场域形 状，容易实现。
单元：互不重叠，覆盖全部场域；每个单元内介质是 单一、均匀的。
节点：网格的交点，待求变量的设置点。 该步骤需要记录的信息： 节点编号、节点坐标 节点属性(激励源、是否边界等) 单元编号 单元节点编号 单元介质
➢ 三角形单元内的基函数 设三角形三个顶点处待求函数值 分别为u1, u2, u3。如果单元足够小， 可以采用线性近似，将单元内任 意p点的u(x,y)表示为
K (e) jm

K (e) mm

b(e) i

f (e)
e
N (e) i
d

3
f
(e)
➢ 第一类边界条件（强加边界条件）
第一类边界节点是指边界上函数值 i fi 已知。因此处理
方法是，合成整体系数阵之后，将该节点所在行的主元素置 1，其它元素均置零，同时将右端项中对应元素设为已知函 数值。
位置是第i行、j列；因此
K
(e) ij
必须合成到整体矩阵的第i行、
j列元素上。
整体矩阵合成：
O
L M M

Kii L Kij L Kim


ui


bi

MO M
M M M

K ji
K jj
K mj

u
j



b
j
b N fd 0
1 2 3 4 5 6 0
以下把单元e的贡献记为
K (e) ij

e
N (e) i
L(
N
(e) j
)d
b(e) i

e
N (e) i
f
(e)d
这样，就有
K00

K (1) 00

K (2) 00

K (3) 00

K (4) 00

K (5) 00
u 的近似解表示为
m
u c ju j
j 1
将近似解代入方程，得余量
m
R Lu f c j Lu j f
j 1
如果余量为零，说明已经满足方程，
即 u 是方程的精确解。
一般情况下余量不为零。 只能放松约束， 强制余量的加权积分为零。 即
wi Rd 0 （ i 1,2, , n ）

式中 wi 为权函数， w1, w2 , , wk , 为权函数序列，
权函数之间要求线性无关。 权函数的不同选择导致不同的近似方法。
2、几种加权余量法
（1）配点法
在求解区域中选取 n 个点 P1, P2 , , Pn ， 让方程的余量在这 n 个点上为零。
即选权函数为
f
(e)d
i
i
由于单元很小，做单元分析时 通常可以取 f (e) 为常数值（可以 认为等于三个顶点上的平均 值）。因此
b(e) i

f (e)
e
Ni(e)d

3
f (e)
单元矩阵：
➢ 上述以节点为序的分析过程对于有限元原理的说明是易于
理解的。而在实际编程中，更有效率的是以单元为序，逐

K (6) 00
K01

K (1) 01

K (6) 01
b0

b(1) 0

b(2) 0

b(3) 0

b(4) 0

b(5) 0

b(6) 0
每个
K
(e ij
)
或
b(e) i
的计算都在具体的单元内单独考虑（称
为单元分析）。
b N f d 右端项元素：
b(e) i

e
N (e) i
个计算单元系数阵[K(e)]，然后合成整体系数阵[K]。单元
系数阵[K(e)]定义为
[K
(e)
]


K (e) ii
K (e) ji

K (e) mi
K (e) ij
K (e) jj
K (e) mj
K (e) im
K (e) jm

K (e) mm

设 i, j, m 是节点的整体编号，元素Kij在整体矩阵中的实际
y)

u1
1 ( x,
y)

u2
2 (x,
y)

u3
3 ( x,
y)
可得
Ni

i (x,
y)
( i 1, 2, 3 )
基函数Ni常被称为插值函数或者形状函数，具有以下性质：
（1）是插值的；
1 (i j) （2）Ni (xj , y j ) 0 (i j) （3）在相邻单元的公共边界上，