为什么磁场求解不正确？

han_jj

请各位帮忙看一下，这是一个求线圈磁场的问题，使用的是边单元117，源选择的是36，进行谐态分析，磁场强度结果不正确，不知道是哪出了问题...是不是36不能用于谐态分析？（以下是命令流）

/prep7
hx=5
hy=5
hz=5
hyx=10
hyy=10
hyz=10   
r0=1   
r1=0.01
dz=0.02
dy=0.01
i=0.1
n=50
f=20000
/pnum,volu,1
/pnum,area,1
/pnum,line,1
/pnum,kp,1
n,51,hyx,hyy,hyz
n,52,hyx,-hyy,hyz
n,53,-hyx,hyy,hyz
n,54,-hyx,-hyy,hyz
n,61,hyx,hyy,-hyz
n,62,hyx,-hyy,-hyz
n,63,-hyx,hyy,-hyz
n,64,-hyx,-hyy,-hyz  
knode,,all
block,-hx,hx,-hy,hy,-hz,hz
v,13,14,15,16,4,2,1,3
v,9,10,12,11,7,8,6,5
v,9,10,13,16,7,8,4,3
v,11,12,14,15,5,6,2,1
v,9,11,15,16,7,5,1,3
v,10,12,14,13,8,6,2,4            ！建立模型
numcmp,all
/replot

et,1,solid117,6
et,2,solid117,0
et,3,sourc36
emunit,mks
mp,murx,1,1.0
mp,rsvx,1,1
mp,rsvy,1,1
mp,rsvz,1,4
mp,murx,2,1.0                !设置单元属性

vsel,s,volume,,1
vatt,1,,1
allsel
vsel,s,volume,,2,7
vatt,2,,2

smrtsize,1
vmesh,all

r,1,1,n*i,dy,dz
*get,maxn,node,,num,max
csys,1
n1=maxn+1
n2=maxn+2
n3=maxn+3
n,n1,r1,0,0
n,n2,r1,90,0
n,n3,0,0,0
type,3
real,1
mp,murx,4,1.0
mp,rsvx,4,3e-8
e,n1,n2,n3
/eshape,1                ！建立36单元

csys,0
nsel,s,loc,x,hyx
nsel,a,loc,x,-hyx
nsel,a,loc,y,hyy
nsel,a,loc,y,-hyy
da,all,az,0
allsel
nsel,s,loc,x,hyx
nsel,a,loc,x,-hyx
nsel,a,loc,y,hyy
nsel,a,loc,y,-hyy
nsel,a,loc,z,hyz
nsel,a,loc,z,-hyz
da,all,volt,0
allsel                    !边界条件

/solu
antype,harm
harfrp,f
sbctran
solve      !求解

han_jj

莫有人知道吗？

cqlyj2008

回复 1# han_jj
3D谐波分析，是基于单元边法的，不能用source36单元

han_jj

谢谢版主提醒，是我把117的自由度弄错了，source36只能用于自由度为AZ的分析中...