求助：ANSYS与ANSOFT励磁分析结果不同（后附APDL命令流）

fivev5

同样励磁电流密度,ANSYS与ANSOFT不同结果，在ANSYS里面使用的单位全部为MKS，请问各位大大，这个是原因？B差两个数量级
   
/prep7
/unit,mks
et,1,plane13
mp,murx,1,1
mp,murx,2,1000
mp,murx,3,1
mp,murx,4,2000
!主轴
axislenth=260e-3
axisr=2.5e-3
axisr1=20e-3
axisl1=160e-3
axisr2=15e-3
axisl2=50e-3
axisl3=50e-3
wp1=axisl2
wp2=axisl2+axisl1
wp3=axisl2+axisl1+axisl3
rectng,axisr,axisr2,0,axisl2
rectng,axisr,axisr1,wp1,wp2
rectng,axisr,axisr2,wp2,wp3
asel,all
aadd,all
numcmp,area
numcmp,line
!线圈架
xqr1=20.5e-3
xqr2=52.5e-3
xqr3=22.5e-3
xql1=wp1+3.5e-3
xql2=xql1+2e-3
xql3=xql2+150e-3
xql4=xql3+2e-3
k,100,xqr1,xql1
k,101,xqr2,xql1
k,102,xqr2,xql2
k,103,xqr3,xql2
k,104,xqr3,xql3
k,105,xqr2,xql3
k,106,xqr2,xql4
k,107,xqr1,xql4
*do,i,100,106
l,i,i+1
*enddo
l,100,107
*do,i,11,18
*if,i,eq,11,then
lsel,s,,,i
*else
lsel,a,,,i
*endif
*enddo
al,all
!线圈
xqnr=xqr3
xqwr=xqr3+20e-3
rectng,xqnr,xqwr,xql2,xql3
!外壳
wkr1=105e-3/2
wkr2=75e-3
wkr3=wkr1+2e-3
wkl1=xql1-13.5e-3
wkl2=wkl1+5e-3
wkl3=wkl2+170e-3
wkl4=wkl3+5e-3
k,200,wkr1,wkl1
k,201,wkr2,wkl1
k,202,wkr2,wkl2
k,203,wkr3,wkl2
k,204,wkr3,wkl3
k,205,wkr2,wkl3
k,206,wkr2,wkl4
k,207,wkr1,wkl4
*do,i,200,206
l,i,i+1
*enddo
l,200,207
*do,i,23,30
*if,i,eq,23,then
lsel,s,,,i
*else
lsel,a,,,i
*endif
*enddo
al,all
!下端盖
dgr1=axisr2
dgr2=22e-3
dgr3=55e-3/2
dgr4=dgr3+5e-3
dgr5=wkr2
dgl1=wkl1-2e-3
dgl2=wkl1
dgl3=dgl2+10e-3
dgl4=wkl1-8e-3
k,300,dgr1,dgl1
k,301,dgr2,dgl1
k,302,dgr2,dgl2
k,303,dgr3,dgl2
k,304,dgr3,dgl3
k,305,dgr4,dgl3
k,306,dgr4,dgl2
k,307,dgr5,dgl2
k,308,dgr5,dgl4
k,309,dgr1,dgl4
*do,i,300,308
l,i,i+1
*enddo
l,300,309
*do,i,31,40
*if,i,eq,31,then
lsel,s,,,i
*else
lsel,a,,,i
*endif
*enddo
al,all
!上端盖
dgr1=axisr2
dgr2=22e-3
dgr3=55e-3/2
dgr4=dgr3+5e-3
dgr5=wkr2
dgl11=wkl4+2e-3
dgl22=wkl4
dgl33=dgl22-10e-3
dgl44=wkl4+8e-3
k,400,dgr1,dgl11
k,401,dgr2,dgl11
k,402,dgr2,dgl22
k,403,dgr3,dgl22
k,404,dgr3,dgl33
k,405,dgr4,dgl33
k,406,dgr4,dgl22
k,407,dgr5,dgl22
k,408,dgr5,dgl44
k,409,dgr1,dgl44
*do,i,400,408
l,i,i+1
*enddo
l,400,409
*do,i,41,50
*if,i,eq,41,then
lsel,s,,,i
*else
lsel,a,,,i
*endif
*enddo
al,all
!轴承
rectng,dgr1,dgr3,dgl2,dgl3
rectng,dgr1,dgr3,dgl22,dgl33
!工具
tnr=2.5e-3
twr=10e-3
tl=50e-3
rectng,tnr,twr,0,-tl
asel,all
arsym,x,all
!工件
wpl=100e-3
wph=20e-3
wjx=2e-3
wpp1=-tl-wjx
rectng,wpl/2,-wpl/2,wpp1,wpp1-wph
!空气块
rectng,dgr5,-dgr5,wpp1-wph-20e-3,wp3+20e-3
!赋予材料特性
aovlap,all
numcmp,area
asel,all
aatt,1
asel,s,,,1
asel,a,,,7
asel,a,,,13,17
asel,a,,,28,29
aatt,2
asel,all
aesize,all,0.01
amesh,all
/solu
asel,s,,,2
esla,s
BFE,all,JS,1,,,3000/(xqwr-xqnr)*(xql3-xql2)
asel,s,,,8
esla,s
BFE,all,JS,1,,,-3000/(xqwr-xqnr)*(xql3-xql2)
asel,all
esel,all
/STATUS,SOLU
solve

                               
登录/注册后可看大图

沙隆巴斯

整差100倍？是不是什么地方参数漏了。米搞成了厘米

fivev5

回复 2# 沙隆巴斯

模型应该没有出现单位上的问题
但是，这里就只设置相对磁导率，
加载的时候就只加了电流密度
其他就没有设置什么参数了，
也不知道还要什么参数了
电磁分析的菜鸟

沙隆巴斯

电流密度也是有单位的

fivev5

回复 4# 沙隆巴斯
谢谢楼上，这里电流密度的单位为A/m2，是MKS系吧

xiaolu

回复 1# fivev5

请问，这个是apdl语言吗，能否在anosft中这样运行参数呢，谢谢！

沙隆巴斯

楼主写的是apdl，ansoft里也有相应的语言，好像是vb之类的

fivev5

现在发现无论三维还是二维
ansoft和ANSYS的结果差别都比较大

zengxiaodong

虽然说Ansys比Ansoft要精确得多，但是差别这么大就不对了，肯定哪里出错了！

保夫鲁沙

很不错，谢谢楼主分享

缘尽怀念

楼主用的是古典的吗