一个计算涡流损耗的例子

lntu_lxm

/UNITS,SI
/prep7
et,1,117
emunit,mks
afexin=3000
pcu1=0.00000003
pcu2=0.00000003
abi=20
tiee=2000
shiyou=1
you=1
emunit,mks
mp,murx,1,shiyou
mp,murx,2,afexin
mp,murx,3,1
mp,rsvx,3,pcu1
mp,murx,4,1
mp,rsvx,4,pcu2
mp,murx,5,1
mp,rsvx,5,pcu1
mp,murx,6,1
mp,rsvx,6,pcu1
mp,murx,7,1
mp,rsvx,7,pcu1
mp,murx,8,1
mp,rsvx,8,pcu1
mp,murx,9,1
mp,rsvx,9,pcu1
mp,murx,10,1
mp,rsvx,10,pcu1
mp,murx,11,abi
mp,murx,12,tiee
mp,murx,13,you
a=1.510
a1=0.62
a2=0.61
b=0.088
c=3.000
d=3.720
e=4.280
d1=0.787
d2=0.912
g1=1.022
g2=1.192
f=0.170
g=0.085
h=0.100
m=0.06
n=0.440
o=0.02
p=0.080
s=3.490
r=1.660
v1=2.050
v2=2.080
v=2.188
k=0.010
l=0.01
gy1=0.0762
gy2=0.5131
gy3=0.2153
gy4=0.2204
dy1=0.0891
dy2=0.521
dy3=0.2849
dy4=0.145
t=-1.00
yl=0.128125
wpoffs,f+a1/2+t,c/2,b+a1+v2/2+g
cylind,0,a1/2,0,s/2,0,360
wpoffs,r,0,0
cylind,0,a/2,0,s/2,0,360
wpoffs,r,0,0
cylind,0,a1/2,0,s/2,0,360
wpoffs,-r-r,0,s/2-a1/2
wprota,0,0,90
cylind,0,a2/2,0,r+r,0,360
vadd,1,2,3,4
wprota,0,0,-90
wpoffs,r,0,-s/2+a1/2
cylind,d1,d2,0,dy4+dy3+dy2+dy1,0,360
cylind,g1,g2,0,gy4+gy2+gy1+gy3,0,360
cylind,a/2,d1,0,v2/2,0,360
cylind,d2,g1,0,v/2,0,360
wpoffs,-r-f-a1/2,0,0
block,0,e,-c/2,c/2,0,d-b-a1-v2/2-g
block,-0.010,e+0.010,-c/2-0.010,c/2+0.010,0,d-b-a1-v2/2-g+0.010
vovlap,all
numcmp,volu
local,11,1,f+a1/2+t+r,c/2,b+a1+v2/2+v2/4+g
vsel,s,volu,,1
vatt,3,1,1
esys,11
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,2
vatt,4,1,1
esys,11
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,3
vatt,13,1,1
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,4
vatt,13,1,1
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,6
vatt,2,1,1
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,5
vatt,11,1,1
esize,0.1
MSHAPE,1,3D
MSHKEY,0
vmesh,all
vsel,s,volu,,7
vatt,13,1,1
esize,0.5
MSHAPE,1,3D
MSHKEY,0
vmesh,all
/solu
antype,harmonic
harfrq,50
asel,s,area,,48                       !加同量平行边界条件
asel,a,area,,49
asel,a,area,,50
asel,a,area,,51
asel,a,area,,52
da,all,az,0
allsel,all
vsel,s,,,1
eslv,s
bfe,all,js,,0,-1.57e+6,0
vsel,s,,,2
eslv,s
bfe,all,js,,0,0.864e+6,0
  
allsel,all
solve
finish
/post1
ESEL,S,MAT,,5
POWERH

ponderadams

计算算得是哪部分涡流损耗？
单单powerh命令就可以了么？
117单元不选keyopt（1）=1？