音圈参数设置请教（ansys计算出来的阻抗结果为理论值的2倍）

soy546

在ansys中进行谐响应分析求音圈阻抗时，阻抗为15.7861，阻抗结果与7.8相差一倍 ，请

帮忙看看哪里出错了！

线圈外径：9.0mm
内径：8.75mm
匝数：34
层数：2
线径：0.055mm
电阻为：7.8欧
ρ=RS/L   

ρ=1.946e-8

S ＝3.14*0.055^2/4*1e-6=2.37e-9m^2
为单匝线圈截面积

R＝7.8Ω
L =3.14*8.875*34=0.95m
截面积为：1.625e-7
填充因子：s/(0.125*1.3*1e-6)=0.4959

/prep7
/title,stranded coil in free space
et,1,53,,,1 ! air
et,2,53,2,,1! voltage forced coil
et,3,110,,,1
emunit,mks
mp,murx,1,1
mp,murx,2,1
mp,rsvx,2,1.946e-8! resistivity of coil
a=4.375e-3
b=1.3e-3
s=0.25e-3
rectng,a,a+s,-b/2,b/2
pcirc,0,6*b,-90,90
pcirc,0,12*b,-90,90
aovlap,all
numcmp,area
/REPLOT
asel,s,area,,1
aatt,2,1,2
asel,s,area,,2
aatt,1,1,3
asel,all

FLST,2,1,5,ORDE,1  ！网格划分

FITEM,2,3

FLST,5,2,4,ORDE,2

FITEM,5,2

FITEM,5,4

CM,_Y,LINE

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y

LESIZE,_Y1, , ,50, , , , ,1

FLST,5,2,4,ORDE,2

FITEM,5,1

FITEM,5,3

CM,_Y,LINE

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y

LESIZE,_Y1, , ,10, , , , ,1

FLST,5,2,4,ORDE,2

FITEM,5,5

FITEM,5,8

CM,_Y,LINE

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y

LESIZE,_Y1, , ,30, , , , ,1

FLST,5,2,4,ORDE,2

FITEM,5,11

FITEM,5,-12

CM,_Y,LINE

LSEL, , , ,P51X

CM,_Y1,LINE

CMSEL,,_Y

LESIZE,_Y1, , ,1, , , , ,1

FLST,2,1,5,ORDE,1

FITEM,2,3

AESIZE,P51X,1e-3,

FLST,5,3,5,ORDE,2

FITEM,5,1

FITEM,5,-3

CM,_Y,AREA

ASEL, , , ,P51X

CM,_Y1,AREA

CHKMSH,'AREA'

CMSEL,S,_Y

MSHKEY,2

AMESH,_Y1

MSHKEY,0

CMDELE,_Y

CMDELE,_Y1

CMDELE,_Y2

r,1, 1.625e-7,34, , ,0.4959
esel,s,mat,,2 ! get coil elements
nsle,s
cp,1,curr,all ! couple curr dof in coil
allsel,all
csys,1
nsel,s,loc,x,12*b
sf,all,inf! set infinite surface flag
csys,0
nsel,s,loc,x,0
d,all,az,0 ! flux-parallel condition
allsel,all
finish
/solu
antype,harm
esel,s,mat,,2
bfe,all,vltg,,1.26! 12 volts load (real)
esel,all
HARF,100,
solve
finish
/post1
ESEL,S,MAT,,2
ETABLE,inductance,NMISC,9
ETABLE,impendance,NMISC,8
SSUM