急！环形线圈中加载电流后 产生的磁通路径不对
本帖最后由 环佩叮当 于 2012-4-15 19:32 编辑finish
/clear
!可选参数!
/prep7
pi=2*asin(1)
csys,1
motor_len=0.413867!电机轴向长度
zz_nj=2.177000 !转子内径
zz_wj=2.257000!转子外径
jh_xs=0.800000!极弧系数
zz_qx=0.003000 !气隙厚度
yct_k=0.120000!永磁体宽
yct_h=0.040000!永磁体高
yct_hpb=0.080285!永磁体周向长度
bt=0.120000 !定子齿宽
js=72 !极对数
dz_c_h=0.085000 !定子槽高
dz_c_b=0.162000!定子槽宽
dz_wj=2.475000!定子外径
dz_nj=2.300000!定子内径
par_a=1
current=303.868563 !电流
current_dens=303*110/dz_c_b/(dz_c_h)!电流密度
ET,1,solid117
mp,murx,1,1
sbh_num=24
sb1=0.000000
sh1=0.000000
sb2=0.200000
sh2=33.440000
sb3=0.300000
sh3=42.200000
sb4=0.400000
sh4=49.200000
sb5=0.500000
sh5=55.330000
sb6=0.600000
sh6=60.510000
sb7=0.700000
sh7=68.470000
sb8=0.800000
sh8=78.820000
sb9=0.900000
sh9=94.750000
sb10=1.000000
sh10=114.730000
sb11=1.100000
sh11=141.720000
sb12=1.150000
sh12=163.220000
sb13=1.200000
sh13=192.680000
sb14=1.250000
sh14=234.870000
sb15=1.300000
sh15=298.570000
sb16=1.350000
sh16=390.130000
sb17=1.400000
sh17=517.520000
sb18=1.450000
sh18=796.180000
sb19=1.500000
sh19=1233.800000
sb20=1.550000
sh20=1990.450000
sb21=1.600000
sh21=3025.480000
sb22=1.650000
sh22=4299.360000
sb23=1.700000
sh23=6050.960000
sb24=1.750000
sh24=7961.780000
TB,BH,2,,sbh_num
*do,i,1,sbh_num,1
tbpt,,sh%i%,sb%i%
*enddo
mp,murx,3,1.090782
mp,mgxx,3,875367.000000
mp,murx,4,1.090782
mp,mgxx,4,-875367.000000
ibh_num=39
ib1=0.100000
ih1=130.000000
ib2=0.200000
ih2=170.000000
ib3=0.300000
ih3=197.000000
ib4=0.400000
ih4=218.000000
ib5=0.500000
ih5=250.000000
ib6=0.600000
ih6=290.000000
ib7=0.700000
ih7=338.000000
ib8=0.800000
ih8=400.000000
ib9=0.850000
ih9=435.000000
ib10=0.900000
ih10=472.000000
ib11=0.950000
ih11=520.000000
ib12=1.000000
ih12=570.000000
ib13=1.050000
ih13=620.000000
ib14=1.100000
ih14=682.000000
ib15=1.150000
ih15=750.000000
ib16=1.200000
ih16=810.000000
ib17=1.250000
ih17=885.000000
ib18=1.300000
ih18=975.000000
ib19=1.350000
ih19=1200.000000
ib20=1.400000
ih20=1600.000000
ib21=1.450000
ih21=2000.000000
ib22=1.500000
ih22=2520.000000
ib23=1.550000
ih23=3020.000000
ib24=1.600000
ih24=3520.000000
ib25=1.650000
ih25=4040.000000
ib26=1.700000
ih26=4760.000000
ib27=1.750000
ih27=6600.000000
ib28=1.800000
ih28=8300.000000
ib29=1.850000
ih29=10000.000000
ib30=1.900000
ih30=12000.000000
ib31=1.950000
ih31=14000.000000
ib32=2.000000
ih32=17000.000000
ib33=2.050000
ih33=20000.000000
ib34=2.100000
ih34=23000.000000
ib35=2.150000
ih35=27300.000000
ib36=2.200000
ih36=32000.000000
ib37=2.250000
ih37=37000.000000
ib38=2.300000
ih38=42500.000000
ib39=2.350000
ih39=48200.000000
TB,BH,5,,ibh_num
*do,i,1,ibh_num,1
tbpt,,ih%i%,ib%i%
*enddo
!!!!!!!!!!绕组布置计算
/prep7
csys,1
zz1= (motor_len-2* bt- dz_c_b)/2
k,1,0,0,zz1
k,2,zz_wj,360*(jh_xs/2-2)/2/js, zz1
k,3,zz_wj,360*(-jh_xs/2-2)/2/js, zz1
aa=zz_wj*sin(360*jh_xs*pi/4/js/180)
k,4,( zz_wj +yct_h),asin(aa/( zz_wj +yct_h))*180/pi-360*2/2/js, zz1
k,5,( zz_wj +yct_h),-asin(aa/( zz_wj +yct_h))*180/pi-360*2/2/js, zz1
larc,2,3,1,( zz_wj) !line:1
lstr,2,4 !line:2
larc,4,5,1,( zz_wj +yct_h) !line:3
lstr,3,5
zz2= motor_len-zz1-yct_k
k,10,0,0, zz2
k,6, zz_wj,360*(jh_xs/2-2)/2/js, zz2
k,7, zz_wj,360*(-jh_xs/2-2)/2/js, zz2
k,8,( zz_wj +yct_h),asin(aa/( zz_wj +yct_h))*180/pi-360*2/2/js, zz2
k,9,( zz_wj +yct_h),-asin(aa/( zz_wj +yct_h))*180/pi-360*2/2/js, zz2
larc,6,7,10, zz_wj
lstr,6,8
larc,8,9,10,( zz_wj +yct_h)
lstr,7,9
k,15,0,0,0
k,11,(zz_wj+yct_h+zz_qx),asin(aa/(zz_wj+yct_h+ zz_qx))*180/pi-360*2/2/js,0
k,12,( zz_wj+yct_h+zz_qx),-asin(aa/( zz_wj+yct_h+zz_qx))*180/pi-360*2/2/js,0
k,13, dz_wj,asin(aa/ dz_wj)*180/pi-360*2/2/js,0
k,14, dz_wj,-asin(aa/ dz_wj)*180/pi-360*2/2/js,0
larc,11,12,15,( zz_wj+yct_h+zz_qx) !line:1
lstr,11,13 !line:2
larc,13,14,15, dz_wj !line:3
lstr,12,14
al,1,2,3,4
voffst,1, yct_k
al,5,6,7,8
voffst,7, yct_k
al,9,10,11,12 !体编号3
voffst,13, motor_len
!al,13,14,15,16 !体编号4
!vdrag,19,,,,,,17
cylind,(dz_nj),(dz_nj+dz_c_h),zz1+bt,zz1+bt+dz_c_b!体编号4
cylind,dz_nj,(dz_nj+dz_c_h),0, zz1!体编号5
cylind,dz_nj,(dz_nj+dz_c_h),motor_len-zz1, motor_len !体编号6
vsbv,3,4
vsbv,7,5
vsbv,3,6
numcmp,volu
cylind, zz_nj, zz_wj,0, motor_len,-360*(jh_xs/2+2)/2/js,360*(jh_xs/2+4)/2/js !转子轭部
allsel,all
vsel,s,loc,x, zz_wj,( zz_wj+ yct_h)
vgen,7,all,,,0,360/2/js,0
allsel,all
vglue,all
numcmp,volu
allsel,all
vsel,s,loc,x,( zz_wj+ yct_h + zz_qx),dz_wj
vgen,4,all,,,0,360/js,0
allsel,all
cylind,(dz_nj+dz_c_h1+dz_c_h2),(dz_nj+dz_c_h),zz1+bt,zz1+bt+dz_c_b,-360*(jh_xs/2+2)/2/js,360*(jh_xs/2+4)/2/js!线圈
cylind,( dz_nj - zz_qx), dz_nj,0,motor_len,-360*(jh_xs/2+2)/2/js,360*(jh_xs/2+4)/2/js!气隙
allsel,all
vglue,all
numcmp,volu
cylind,zz_nj-0.1,dz_wj+0.1,-0.1,motor_len+0.1,-360*(jh_xs/2+2)/2/js-0.5,360*(jh_xs/2+4)/2/js+3 !空气包围圈
numcmp,volu
allsel,all
VOVLAP,all
!把气隙单独出来
asel,s,loc,z, motor_len+0.1
vsla,s,0
cm,air,volu
yc=7
*do,i,1,yc,1
vsel,s,loc,x, zz_wj,( zz_wj+ yct_h)
vsel,r,loc,z,zz1,zz1+ yct_k
ccc=(i-3)*360/2/js
vsel,r,loc,y,-360*jh_xs/4/js+ccc,360*jh_xs/4/js+ccc
cm,a%i%,volu
*enddo
cmsel,s,a1
*do,i,3,7,2
cmsel,a,a%i%
*enddo
cm,pmaN,volu
cmsel,s,a2
cmsel,a,a%4%
cmsel,a,a%6%
!cmsel,a,a%8%
cm,pmaS,volu
*do,j,1,yc,1
vsel,s,loc,x, zz_wj,( zz_wj+ yct_h)
vsel,r,loc,z,motor_len-zz1-yct_k, motor_len-zz1
ddd=(j-3)*360/2/js
vsel,r,loc,y,-360*jh_xs/4/js+ddd,360*jh_xs/4/js+ddd
cm,b%j%,volu
*enddo
cmsel,s,b1
*do,i,3,7,2
cmsel,a,b%i%
*enddo
cm,pmbS,volu
cmsel,s,b2
cmsel,a,b%4%
cmsel,a,b%6%
!cmsel,a,b%8%
cm,pmbN,volu
cmsel,s,pmaN !选择永磁体
cmsel,a,pmbN
cm,pmN,volu
cmsel,s,pmaS
cmsel,a,pmbS
cm,pmS,volu
allsel,all
cmsel,u,air
asel,s,loc,z,0 !选择线圈
asel,a,loc,z, motor_len
vsla,s
cm,temp,volu
allsel,all
cmsel,u,air
cmsel,u,temp
cmsel,u,pmS
cmsel,u,pmN
cm,conduct,volu
allsel,all
vsel,s,loc,x, dz_nj, dz_wj
cmsel,u,conduct
cmsel,u,air
cm,guigang2,volu !定子铁心
allsel,all !选择转子上的硅钢
lsel,s,loc,x, zz_nj
lsel,r,loc,z,0
asll,s,0
vsla,s,0
cm,temp2,volu
allsel,all
cmsel,s,temp2
cmsel,u,air
cm,guigang1,volu
allsel,all !选择气隙
cmsel,u,air
cmsel,u,guigang2
cmsel,u,conduct
cmsel,u,guigang1
cmsel,u,pmS
cmsel,u,pmN
cm,qixi,volu
local,11,1,0,0,0
allsel,all
cmsel,s,qixi
vatt,1,,1,0
cmsel,s,guigang2
vatt,2,,1,0
cmsel,s,guigang1
!vatt,5,,1,0
vatt,2,,1,0
cmsel,s,conduct
vatt,1,,1,0
cmsel,s,pmS
!vatt,3,,1,11
vatt,1,,1,0
cmsel,s,pmN
!vatt,4,,1,11
vatt,1,,1,0
allsel,all
cmsel,s,air
vatt,1,,1,0
allsel,all
cmsel,s,qixi
aslv,s
lsla,s
LESIZE,all,20e-3, , , , , , ,1
cmsel,s,pmS,volu
cmsel,a,pmN,volu
aslv,s
lsla,s
LESIZE,all,60e-3, , , , , , ,1
cmsel,s,conduct,volu
aslv,s
lsla,s
LESIZE,all,60e-3, , , , , , ,1
cmsel,s,guigang2,volu
cmsel,a,guigang1,volu
aslv,s
lsla,s
LESIZE,all,60e-3, , , , , , ,1
cmsel,s,air,volu
aslv,s
lsla,s
LESIZE,all,100e-3, , , , , , ,1
allsel,all
MSHAPE,1,3D
vmesh,all
csys,1
!boundary
allsel,all
ASEL,S,EXT
da,all,all,0,0
allsel,all
!save,tfpm,db
cmsel,s,conduct
bfv,all,JS,0,-current_dens,0
!BFA,all,JS, , ,-current_dens*cos((a_fif)*pi/180)
allsel,all
/solu
!solution
solve
/post1
csys,1
rsys,1
allsel,all
!!!!!!!!!!定子轭平均磁密，最大，最小磁密计算
cmsel,s,guigang2,volu
vsel,r,loc,y,-asin(aa/(zz_wj+yct_h+ zz_qx))*180/pi,asin(aa/(zz_wj+yct_h+ zz_qx))*180/pi
eslv,s
wpoffs,0,0,motor_len/2
sucr,phi2,cplane
sumap,,b,z
supl,phi2,
sueval,flux2,BZ,INTG
sta_cm_avg=abs(flux2/yct_hpb/(dz_wj-dz_nj-dz_c_h))
wpoffs,0,0,-motor_len/2
*cfopen,Noload_Field_Result,txt,,append
*cfwrite,sta_e_cm,flux2
*cfclose
!!!!!!!!!!定子轭平均磁密，最大，最小磁密计算
cmsel,s,guigang2
cmsel,a,pmS
cmsel,a,pmN
vsel,r,loc,y,-asin(aa/(zz_wj+yct_h+ zz_qx))*180/pi,asin(aa/(zz_wj+yct_h+ zz_qx))*180/pi
eslv,s
/ANG,1,90,XS,1
/ANG,1,90,zS,1
PLVECT,B, , , ,VECT,ELEM,ON,0
大家好！我想在环形线圈中加载电流，永磁体的属性被我幅值为跟空气一样，只观察线圈电流产生的磁场，发现线圈电流产生的磁场没有经过转子，直接从线圈槽口穿过去了，很迷惑，求助！ 直接贴图，上程序一时看不明白。
页:
[1]