ansys计算电感

ansys计算电感
模型文件chock.lgw如下,如何求解线圈电感
/title static analysis of the chock
keyw,magedg,1
/prep7
define the elements and materials
et,1,117
mp,murx,1,the air
mp,murx,2,the gap
mp,murx,3,the coil
mp,rsvx,3,3e-8
mp,rsvx,4,60e-4 the iron
mp,murx,4,12000
ix=0.0457
iy=0.1524
iz=0.127
gz=0.0254
fz=0.0762
fy1=0.0254
fy2=0.0308
ch=0.0229
cw=0.018
gap1=0.005
gap2=0.002
i=200
n=36
jdens=i*n/(ch*cw)
*dim,cur,array,2
*dim,coils,char,2
cur(1)=200
coils(1)='coil1'
cur(2)=200
coils(2)='coil2'
/pnum,volu,1
block,0,ix/2,0,fy2/2,-iz/2,-iz/2+gz
block,0,ix/2,fy2/2,iy/2-fy1,-iz/2,-iz/2+gz
block,0,ix/2,iy/2-fy1,iy/2,-iz/2,-iz/2+gz
block,0,ix/2,iy/2-fy1,iy/2,-iz/2+gz,0
vsymm,x,1,4
vsymm,y,1,8
vsymm,z,1,16
block,0,ix/2,0,fy2/2,-fz/2,-fz/4-ch/2
block,0,ix/2,0,fy2/2,-fz/4-ch/2,-fz/4+ch/2
block,0,ix/2,0,fy2/2,-fz/4+ch/2,fz/4-ch/2
block,0,ix/2,0,fy2/2,fz/4-ch/2,fz/4+ch/2
block,0,ix/2,0,fy2/2,fz/4+ch/2,fz/2-gap1
block,0,ix/2,0,fy2/2,fz/2-gap1,fz/2
vsymm,x,33,38
vsymm,y,33,44
vglue,all
numcmp,all
block,ix/2,ix/2+gap2,0,fy2/2,-fz/4-ch/2,-fz/4+ch/2
block,0,ix/2,fy2/2,fy2/2+gap2,-fz/4-ch/2,-fz/4+ch/2
cyl4,ix/2,fy2/2,0,0,gap2,90,-fz/4-ch/2
cyl4,ix/2,fy2/2,0,0,gap2,90,-fz/4+ch/2
vsbv,59,60
numcmp,all
vsymm,x,57,59
vsymm,y,57,62
vsymm,z,57,68
vsel,s,volu,57,80
vglue,all
numcmp,all
block,ix/2+gap2,ix/2+gap2+cw,0,fy2/2,-fz/4-ch/2,-fz/4+ch/2
block,0,ix/2,fy2/2+gap2,fy2/2+gap2+cw,-fz/4-ch/2,-fz/4+ch/2
cyl4,ix/2,fy2/2,gap2,0,gap2+cw,90,-fz/4-ch/2
cyl4,ix/2,fy2/2,gap2,0,gap2+cw,90,-fz/4+ch/2
vsbv,83,84
numcmp,all
vsymm,x,81,83
vsymm,y,81,86
vsymm,z,81,92
vsel,s,volu,81,104
vglue,all
numcmp,all
plot the air
block,-ix*1.5,ix*1.5,-iy*1.5,iy*1.5,-iz*1.5,iz*1.5
vsel,all
vovlap,all
numcmp,all
mesh
vsel,s,volu,14,36,36-14
vsel,a,volu,40,52,52-40
vatt,2,1
vsel,s,volu,1,13
vsel,a,volu,15,35
vsel,a,volu,37,39
vsel,a,volu,41,51
vsel,a,volu,53,56
vatt,4,1
vsel,s,volu,57,80
vatt,1,1
vsel,s,volu,81,84
vsel,a,volu,89,92
vsel,a,volu,97,100
vatt,3,1,1
vsel,s,volu,85,88
vsel,a,volu,93,96
vsel,a,volu,101,104
vatt,3,2,1
vsel,s,volu,105
vatt,1,1
vsel,s,volu,1,104
vmesh,all
vsel,s,volu,105
mshkey,0
mshape,1,3d
smrtsize,4
vmesh,all
esel,s,real,1
cm,coils(1),elem
esel,s,real,2
cm,coils(2),elem
allsel,all
finish
/solu
allsel,all
nsel,s,loc,x,ix*1.5
nsel,a,loc,y,-iy*1.5
d,all,az,0
vsel,s,volu,84,100,100-84
vsel,a,volu,88,104,104-88
eslv,s
bfe,all,js,jdens
vsel,s,volu,81,97,97-81
vsel,a,volu,85,101,101-85
eslv,s
bfe,all,js,0,jdens
vsel,s,volu,82,98,98-82
vsel,a,volu,86,102,102-86
eslv,s
bfe,all,js,-jdens
vsel,s,volu,83,99,99-83
vsel,a,volu,87,103,103-87
eslv,s
bfe,all,js,0,-jdens
csys,1
vsel,s,volu,92,96,96-92
eslv,s
bfe,all,js,0,jdens
vsel,s,volu,91,95,4
eslv,s
bfe,all,js,0,jdens
vsel,s,volu,89,93,4
eslv,s
bfe,all,js,0,jdens
vsel,s,volu,90,94,4
eslv,s
bfe,all,js,0,jdens
csys,0
allsel,all
solve
finish
该源码过于冗繁,现在,不管是谁,我就给他分数
命令流修改部分我没有收到,能否将修改的命令流和求解电感的命令流发到我的邮箱里个,我的邮箱是gfgf_99@163.com
请注意,悬赏分数已经提高到150