用solid97和circu124求解导体上的电感，结果不对，求赐教！
本帖最后由 zhengfengfen 于 2011-10-10 17:59 编辑用solid97和circu124求解导体上的电感，电路和有限元模型都很简单：电路仅为一个电压源；而有限元模型为一段导线，边界定义为无限远边界。但是，系统先是在模型中给出：
Large deflection effects are not included for element type SOLID97.
Large deflection effects are not included for element type INFIN111.
Large deflection effects are not included for element type CIRCU124.
的警告。然后在求解所得的结果也不正确：电流很小，趋势也不对。请诸位帮忙，看看问题到底出在什么地方，先谢谢了！
所附为代码：
finish
/clear
/filname,inductor
/title, test conductor inductance
/prep7
emunit, mks
Pi=acos(-1)
my_volt =1000
MP, RSVX, 2, 4e-6 !the resistivity of the conductor
MP, MURX, 2, 1!the permibility of the conductor
MP, RSVX, 3, 1e6
MP, MURX, 3, 1 !the permibility of air
MP, MURX, 4, 1
r, 10, Pi*6.4e-3*6.4e-3 !define the cross sectional area of coil
et, 1, solid97 ! air
et, 2, solid97, 4 ! current forced coil
et, 3, infin111, 1, 0
et, 6, circu124, 4, 4
et, 7, circu124, 7
n, 1
n, 2
n, 3
wpoffs, 0, 0, -15e-3
blc5, 0, 0, 30e-3, 30e-3, 30e-3
wpoffs, 0, 0, 10e-3
cyl4, 0, 0, , , 6.4e-3, , 10e-3
wpoffs, 0, 0, 5e-3
vsel,all
vovlap, all
!end of solid modeling
allsel, all
vsel, s, , , 2
vatt, 2, 10, 2
allsel, all
vsel, s, , , 3
vatt, 3, , 1
smrtsize, 3
mshape, 1, 3D
mshkey, 0
allsel, all
vmesh, all
!Construct the inifinite region
allsel, all
asel, s, loc, z, 15e-3
asel, a, loc, z, -15e-3
asel, a, loc, x, 15e-3
asel, a, loc, x, -15e-3
asel, a, loc, y, 15e-3
asel, a, loc, y, -15e-3
cm, aouter, area
extopt, esize, 1
MAT, 4
type, 3
vext, all, , , , , , 1.1, 1.1, 1.1
aslv
cmsel,u, aouter
nsla, s, 1
sf,all,inf
!get the upper surface node on the conductor
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, 5e-3
*get, ppower, node, 0, num, min
!get the lower surface node on the conductor
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, -5e-3
*get, npower, node, 0, num, min
!end of finding surface node
TYPE,7 ! mass conductor
r, 1, 1
real, 1
E, 1, 2, ppower, npower
TYPE,6 ! VOLTAGE SOURCE
R,4, 0, 0, 900e-6, my_volt! PIECEWISE LINEAR LOAD
REAL,4
E, 1, 2, 3 ! VOLTAGE SOURCE ELEMENT
!assign the boundary conditions
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, 5e-3
sf, all, MCI, -1
nsel, a, , , 1
cp, 1, volt, all
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, -5e-3
sf, all, MCI, +1
cp, 3, volt, all
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, -5e-3
cm, my_c1, node
allsel, all
vsel, s, mat, , 2
nslv, s, 1
nsel, r, loc, z, 5e-3
cm, my_c2, node
allsel, all
cmsel, s, my_c1
cmsel, a, my_c2
cp, 4, CURR, all
!end of assign
save
finish
/solu
ANTYPE,TRANS
NROPT, FULL
TRNOPT,FULL
NLGEOM, ON
neqit, 50
CNVTOL, CSG, , 0.5e-8
allsel, all
OUTRES,ALL,ALL
DELTIM, 300e-6
TIME, 900e-6
SOLVE
FINISH
/post26
esol, 2, 15019, , smisc, 2
plvar, 2,
页:
[1]