第
29
卷第
3
期
2008
年
6
月
河南科技大学学报:自然科学版
Journal
of
Henan
University
of
Science
and
Technology:
Natural
Science
Vo
l.
29
No.3
Jun.
2008
文章编号:
1672 - 6871
(2008)
03
- 0091 -
04
无网格伽辽金法在静电场中的应用
霍允,夏茂辉,李俊贤
(燕山大学理学院,河北秦皇岛
066004)
摘要
t
无网格伽辽金法(
EFGM)
是近几年发展起来的与有限元相似的数值算法,并在电磁场分析中得到初步
的应用。本文采用移动的最小二乘法构造形函数,从能量泛函的弱变分形式中得到控制方程,并用
Lagrange
乘子处理本质边界条件,从而得到数值解。基于
MATLAB
平台实现了一维静电场中的
EFGM
算法,并与解析
解进行比较,结果表明了
EFGM
算法的正确性和有效性。
关键词:元网格;移动最小二乘法;
Lagrange
乘子;静电场数值分析
;MATLAB
中图分类号
;TM153.1
文献标识码
:A
。
前言
无网格法最早出现在
1977
年,此方法不需要划分单元,只需要节点信息,且具有求解精度高、收敛
快、后处理方便等优点,
1990
年以来在力学领域发展迅速,目前代表性的无网格方法
[1
-
2J
有;光滑粒子
法
(SPH)
,扩散单元法
(DEM)
,无网格伽辽金法
[3
-4J
(EFGM)
,单元分解法
(PUM)
,再生核粒子法
(RKPM)
等。其中,无网格伽辽金法是美国西北大学的
Belytschko
教授提出的,目前在国内外电磁场领
域有关此方法的研究仍处在起步阶段
[5-7]o
本文以力学领域的研究成果为基础,探讨了静电场中的无
网格伽辽金法。首先,介绍元网格伽辽金法的理论基础,然后,给出二维静电场离散方程的实现过程,最
后,为简单起见,以一维静电场为例说明该方法在电场中是可行且有效的。
1
移动最小二乘近似
无网格伽辽金法采用移动最小二乘近似
(MLS)
构造形函数。
考虑域。中一场变量
u(X)
的未知标量函数,定义在
X
处
u(X)
的
MLS
近似表达式为
uh(X)=;
二
Pj(X)
α
:/X)
= p
T
(x)
α
(X)
式中
p(X)
为对于三维空间坐标
X
T
=
[耳,
yJ
的基函数,该基函数
p(X)
为完备多项式。
式(1)中的
α
(X)
为一系数向量,可表示为
α
T
(X)
= 1 a
1
(X)
α2
(X)
…
am(X)
f
(2)
式
(2)
中的系数向量
α
(X)
是
X
的函数,该系数
α
可通过下列加权离散
L
2
范数取极小值而得到
、‘』
J
咽'且
'''飞
J
=;二
ω
(X
-
X
,)
[pT(X
,)
α
(X)
- U
.J
2
(3)
式中
n
为包含在权函数
ω
(X
-
X
,)
笋
0
的支持域中的节点数
,
U
i
为
U
在
X=X
i
处的节点参数。式
(3)
为一加权残量的泛函,它利用未知场变量节点参数值和其近似而构成
,
J
关于
α
(X)
求驻点值可得到
aJ
=
0
θα
(4)
该运算将产生如下线性关系
A(X)
α
(X)
=
B(X)
U
(5)
基金项目:凉山大学博士基金项目
(B272)
作者简介:霍
允
(1983
-
)
.女,河北石家庄人,硕士;夏茂辉(1
963
-
)
.男,黑龙江鹤岗人,教授,博士,硕士生导师,主要研究方向
为有限元、元网格算法.
收稿日期
:2007
-
10
-
28