20 3 Vol． 20 No． 3

2011 9 Computer Aided Engineering Sept． 2011

CAD / CAE / CAM 1006 － 0871 2011 03-0092-05

Independent CAD / CAE / CAM Innovation Forum

FEPG
1 2

1． 100190
2． 100190

FEPG．

FEPG 5
． FEPG 、
、
、 ． FEPG
、 ．
FEPG
TB115 B

Finite element analysis software platform FEPG

LIANG Guoping1 TANG Juzhen2
1． Academy of Mathematics and Systems Science Chinese Academy of Sciences Beijing 100190 China
2． Beijing FEGEN Software Co． Ltd． Beijing 100190 China

Abstract As to the fact that the most finite element softwares are generally applied to certain fields and
types of problems which cannot meet the requirement of workers in the other fields especially in the
fields of research and teaching finite element analysis software platform FEPG is developed which can
automatically generate finite element programs by the techniques such as finite element language． Finite
element language is introduced which are used to describe finite element problems． By the language
expressions of partial differential equations and finite element algorithms which are used to solve
differential equations can be described by simply inputting the desired expressions and formulas in finite
element language users can get all the finite element computation source programs by automatic
generation function of FEPG five examples are indicated to introduce the way of preparing partial
differential equations and algorithms in finite element language． The main techniques used in FEPG are
introduced including generator componentization formula library Compared with the library of

2011-05-19 2011-05-27
“ ” 2009AA044501
1938— E-mail linjilingcn@ yahoo． com． cn

http / / www． chinacae． cn

 3 FEPG 93

traditional finite element software formula library is easy to be read modified and maintained by users
and can decrease software complexity and improve software extensibility． Compared with the domestic and
foreign major general-purpose finite element softwares the development with FEPG can be achieved with
less code and investment and short development cycle but the interfaces of FEPG are not friendly and
users are required with a higher starting point．
Key words FEPG finite element finite element language software development software platform
generator componentization formula library

0 FEPG
Partial Differential
、 Equations PDE
、

． ． FEPG

、 ．
．
FEPG

．
． 5
FEPG 2 ．
1
、
①
FEPG － · k  u = q in Ω
FORTRAN

{
． u x y = g on Γ1
u
－k = a u －b on Γ2
． FEPG n

、 ∫ k ( ux δxu + uy δyu ) dΩ = ∫ qδudΩ + ∫ k unδudΩ

Ω Ω Ω
．
PDE
1 disp u
coor x y
shap q 4
gaus q
mate ek eq 1． 0d0 0． 0
． stif

① 2003 “FEPG FEPG ”．

http / / www． chinacae． cn

94 2011

dist = + u / x u / x * ek + u / y u / y * ek func ex ey exy

load = + u * eq shap % 1 % 2
end gaus % 3
GCN mate pe pv fu fv
defi S| c6 fact = pe / 1． + pv / 1． － 2． * pv
a ell
vect u u v
Startsin a
vect ed ex ey
Solvsin a
vect fs fu fv
2 PDE
matr de 2 2

{
1． － pv pv
σ xx σ xy pv 1． － pv
+ + fx = 0
x y func
σ xy σ yy ex = + u / x
+ + fy = 0
x y ey = + v / y
exy = + u / y + v / x
u v u v stif
ε xx = ε = ε = +
x yy y xy y x dist = + ed_i ed_j * de_i_j* fact + exy exy *
0． 5 － pv * fact
 σ xx  load = + u_i * fs_i
σ  = E
· end
 yy  1 +ν 1 － 2ν
 σ xy  3
N-S
1 － ν ν 0   ε xx 
 ν

 0
1 －ν 0   ε yy 
   {－ μΔ u + u· u +  p = f
·u = 0
0 0． 5 － ν   ε xy 

PDE
u = u0 v = v0 defi
disp u v p
T x = f1 T y = f2 coor x y
coef un vn
T x = f1 u v w T y = f2 u v w func div
shap q 9

∫ [ ( σx )
σ xy shap q 4 p 1 2 3 4
xx
+ + f x δu + gaus q
V y

( ) ]
σ xy σ yy mate emu rou fx fy \1． 0d － 3 1． 0d0 0 0
+ + f y δv dV = 0 gvs q4
x y
func
div = + u / x + v / y
∫
V
σ xx δε xx + σ yy δε yy + σ xy δε xy dV = stif
dist = + u / x u / x * emu + u / y u / y * emu +
∫V
f x δu + f y δv dV + ∫
Γ
T x δu + T y δv dΓ
v / x v / x * emu + v / y v / y * emu + div p －
p div + p / x p / x * 1． 0d － 10 +
． PDE p / y p / y * 1． 0d － 10
① load = + u * fx* rou + v * fy* rou
disp u v End
coor x y 4

① 2003 “FEPG FEPG ”．

http / / www． chinacae． cn

 3 FEPG 95

{
－ · K  φ = Q A = 0 φ = 0
φ K K = 1． 0 Q n· － σjωA － σ  φ = 0
Q = 0．
． PDE
defi

{
－ · K  φ + mφ = Q disp Ax Ay Az u
coor x y z
φ≥y K = 0 m = ε
func curlx curly curlz div
φ ＜ y K = 1． 0 m = 0
shap % 1 % 2
PDE
gaus % 3
defi
mass % 1 1． 0d0
disp u
mate sigma epsilon omega
coor x y z
fmu fjx fjy fjz
coef un
vect A Ax Ay Az
shap w 4
vect x x y z
gaus % 3
vect fj fjx fjy fjz
mate pk eq 1． 0 0． 0
vect curl curlx curly curlz
stif
func
S| c6 pk1 = pk curlx = + Az / y － Ay / z
S| c6 em = 0． 0d0 curly = + Ax / z － Az / x
curlz = + Ay / x － Ax / y
S| c6 if un． lt． y + 1． 0d － 11 then
div = + Ax / x + Ay / y + Az / z
S| c6 pk1 = 0． 0d0 stif
S| c6 em = 1． 0d － 8 S| c6 eu1 = omega* sigma

S| c6 endif S| c6 eu2 = 1 / fmu

dist = + u / x u / x * pk1 + u / y u / y * pk1 + dist = + curl_i curl_i * | eu2 0． 0 | + div div *
u / z u / z * pk1 + u u * em | eu2 0． 0 | + A_i A_i * | 0． 0 eu1 | +
load = － u * eq u / x_i A_i * | sigma 0． 0 | + A_i u / x_i *
end | 0． 0 eu1 | + u / x_i u / x_i * | sigma 0． 0 |
1． load = + A_i * fj_i* | 1． 0 0． 0 |
end

2
FEPG 2 、
3
．

1
Fig． 1 Analysis on 3D unconfined seepage flow
field of earth dam

5 A－φ PDE 2

( 1μ  × A ) －  ( 1μ ·A ) +
Fig． 2 Generator technology
 ×

 σjωA + σ  φ = 0 、 、
 FORTRAN
 · ( － σjωA － σ  φ ) = 0

( 1μ  × A ) －  ( 1μ ·A ) = J
90% ．
 × s
FEPG
、

http / / www． chinacae． cn

96 2011

． ． FEPG 、

． ． FEPG 1．
FEPG 1 FEPG
． Tab． 1 Comparison of FEPG with general finite

、 ． element software
FEPG

FEPG

、 、PDE 、
、
、 、 、 “ ” “ ”

N-S 、Dacy
A-φ ．
FEPG
4
、
FEPG
FORTRAN 90% ．
．
PDE
、
、 、 ．

FEPG
CPU CPU
．

3 ．

、
ADINA ANSYS Abaqus MSC Marc
．
COMSOL ．
FEPG
．
FEPG
FORTRAN 、 、 ．
．
FEPG
．
、

1 ZIENKIEWICZ Q C TAYLOR R L ZHU J Z． The finite element method its basis and fundamentals M ． 6th ed． Oxford Butterworth-
Heinemann 2008 733．
2 ． M ． 2009 1-80．
3 ． J ． 1985 2 2 121-125．
LIANG Guoping FU Zizhi． The component programming method and its language J ． Comput Struct Mech ＆ Applications 1985 2 2
121-125．

http / / www． chinacae． cn