编了一个三维涡流场程序，大家帮忙看看哪里不对，谢了
fini/clear
/title,BKSC15000/38.5
keyw,magedg,1 !指定磁场分析
/prep7
/pnum,volu,1 !显示体积编号
/pnum,area,1 !显示面编号
et,1,solid117
et,2,solid117,1 !空气的单元特性
et,3,mesh200,7
mp,murx,1,1 !空气相对磁导率
mp,murx,2,1 !绕组相对磁导率
mp,murx,3,7240 !铁心相对磁导率
mp,rsvx,3,0.73e-6 !铁心电阻率
mp,murx,4,800 !夹件相对磁导率 Q235钢
mp,rsvx,4,1.43e-8 !夹件电阻率
!i=224.9 !额定电流A
i1=318 !选择了三相电流额定电流为224.9 224.9X根号2Xsin90
i2=-159 !224.9X根号2Xsin210
i3=-159 !224.9X根号2Xsin330 这样可以不？
acoila=184 !每匝导线截面积 平方毫米
acoilb=198.72
acoilc=160.8
!A项
jdena1=i1/acoila*1000000 !电流密度 A/m2 线圈在绕制过程中换线，有三种线规，截面积不同
jdena2=i1/acoilb*1000000
jdena3=i1/acoilc*1000000
!B项
jdenb1=i2/acoila*1000000
jdenb2=i2/acoilb*1000000
jdenb3=i2/acoilc*1000000
!C项
jdenc1=i3/acoila*1000000
jdenc2=i3/acoilb*1000000
jdenc3=i3/acoilc*1000000
!建立模型
local,12,1,0,0,0
wpcsys,,12
!铁心饼 V1-V6 1/4模型 电抗器有铁心饼之间有气隙
cyl4,0,0,0.04,0,0.11,180,0.05 !A内1 V1
cyl4,1.56,0,0.04,0,0.11,180,0.05 !B内1 V2
cyl4,3.12,0,0.04,0,0.11,180,0.05 !C内1 V3
cyl4,0,0,0.14,0,0.28,180,0.05 !A外1 V4
cyl4,1.56,0,0.14,0,0.28,180,0.05 !B外1 V5
cyl4,3.12,0,0.14,0,0.28,180,0.05 !C外1 V6
vsel,s,,,1 !复制A内饼18个 V1/V7-V23
vgen,9,1,,,,,0.05+0.01
vgen,2,14,,,,,0.05+0.022
vgen,9,15,,,,,0.05+0.01
vsel,s,,,2 !复制B内饼18个 V2/V24-V40
vgen,9,2,,,,,0.05+0.01
vgen,2,31,,,,,0.05+0.022
vgen,9,32,,,,,0.05+0.01
vsel,s,,,3 !复制C内饼18个V3/V41-V57
vgen,9,3,,,,,0.05+0.01
vgen,2,48,,,,,0.05+0.022
vgen,9,49,,,,,0.05+0.01
vsel,s,,,4 !复制A外饼18个 V4/V58-V74
vgen,9,4,,,,,0.05+0.01
vgen,2,65,,,,,0.05+0.022
vgen,9,66,,,,,0.05+0.01
vsel,s,,,5 !复制B外饼18个 V5/V75-V91
vgen,9,5,,,,,0.05+0.01
vgen,2,82,,,,,0.05+0.022
vgen,9,83,,,,,0.05+0.01
vsel,s,,,6 !复制C外饼18个 V6/V92-V108
vgen,9,6,,,,,0.05+0.01
vgen,2,99,,,,,0.05+0.022
vgen,9,100,,,,,0.05+0.01
!上铁轭 铁轭之间也有气隙 夹件是Q235钢
block,-0.28,3.4,0.03,0.11,1.106,1.426 !V109
block,-0.28,3.4,0.138,0.2825,1.106,1.426 !V110铁心铁轭总数量V1-V110
!线圈
cyl4,0,0,0.425,0,0.4445,180,0.755 !A内1 最内层线圈 体号 V111 每一相线圈分6层
cyl4,0,0,0.4695,0,0.491,180,0.723 !A2 V112
cyl4,0,0,0.516,0,0.5355,180,0.677 !A3 V113
cyl4,0,0,0.5605,0,0.58,180,0.646 !A4 V114
cyl4,0,0,0.605,0,0.6245,180,0.599 !A5 V115
cyl4,0,0,0.6495,0,0.679,180,0.565 !A6 V116
cyl4,1.56,0,0.425,0,0.4445,180,0.755 !B柱线圈 V117-V122
cyl4,1.56,0,0.4695,0,0.491,180,0.723
cyl4,1.56,0,0.516,0,0.5355,180,0.677
cyl4,1.56,0,0.5605,0,0.58,180,0.646
cyl4,1.56,0,0.605,0,0.6245,180,0.599
cyl4,1.56,0,0.6495,0,0.679,180,0.565
cyl4,3.12,0,0.425,0,0.4445,180,0.755 !C柱线圈 V123-V128
cyl4,3.12,0,0.4695,0,0.491,180,0.723
cyl4,3.12,0,0.516,0,0.5355,180,0.677
cyl4,3.12,0,0.5605,0,0.58,180,0.646
cyl4,3.12,0,0.605,0,0.6245,180,0.599
cyl4,3.12,0,0.6495,0,0.679,180,0.565
!夹件
block,-0.88,4,0.3125,0.324,1.116,1.416 !上夹件V129
!边界 空气
block,-2,5.12,0,1.679,0,2.426 !V131 干式的电抗器，可以加一个远场的边界么？怎么加？
vsel,s,,,all
vovlap,all
!赋属性
/pnum,mat,1
vsel,s,,,1,110,1 !铁心与铁轭
vatt,3,,1
allsel,all
vsel,s,,,111,128,1 !线圈
vatt,2,,1
allsel,all
vsel,s,,,129 !夹件
vatt,4,,2
allsel,all
vsel,s,,,131 !空气区域
vatt,1,,1
allsel,all
!划分网格
allsel,all
esize,0.1
vsel,s,,,1,129,1
mshkey,0
mshape,1,3d
vmesh,all
smrt,5
vsel,s,,,131
aesize,all,0.1
mshape,1,3D
mshkey,0
vmesh,all
!给线圈施加电流密度
vsel,s,,,111 !A相六层线圈
vsel,a,,,113,115,1
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdena1,0,0
vsel,s,,,112
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdena2,0,0
vsel,s,,,116
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdena3,0,0
vsel,s,,,117 !B相六层线圈
vsel,a,,,119,121,1
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenb1,0,0
vsel,s,,,118
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenb2,0,0
vsel,s,,,122
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenb3,0,0
vsel,s,,,123 !C相六层线圈
vsel,a,,,125,127,1
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenc1,0,0
vsel,s,,,124
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenc2,0,0
vsel,s,,,128
eslv,s
emodif,all,esys,12
bfe,all,js,1,,jdenc3,0,0
finish
/solu
antype,harmic,new
harfrq,50
allsel,all
nsel,s,loc,x,-2 !这样设定边界条件对不？
nsel,a,loc,x,5.12
nsel,a,loc,y,0
nsel,a,loc,y,1.679
nsel,a,loc,z,0
nsel,a,loc,z,2.426
d,all,az,0
!allsel,all !整个底面
!asel,s,,,781 !整个底面减去铁心中间剖面
!asel,a,,,5,6,1 !铁心剖平面
!asel,a,,,11,12,1
!asel,a,,,17,18,1
!asel,a,,,23,24,1
!asel,a,,,29,30,1
!asel,a,,,35,36,1
!asel,a,,,41,137,6
!asel,a,,,42,138,6
!asel,a,,,347,443,6
!asel,a,,,348,444,6
!asel,a,,,449,545,6
!asel,a,,,143,239,6
!asel,a,,,144,240,6
!asel,a,,,450,546,6
!asel,a,,,551,647,6
!asel,a,,,245,341,6
!asel,a,,,246,342,6
!asel,a,,,552,648,6
!asel,a,,,665,695,6
!asel,a,,,666,696,6
!asel,a,,,701,731,6
!asel,a,,,702,732,6
!asel,a,,,737,767,6
!asel,a,,,738,768,6
!da,all,az,0
allsel
solve
finish
/post1
esel,s,mat,,4 !选择材料特性为4的夹件
powerh !涡流损耗
plvect,b,,,,vect,elem,on !磁密等值线
/graphics,power
avres,2
plnsol,b,sum
/device,vectop,1 !绘出磁通密度矢量
plvect,b
/device,vectop,1 !绘出磁场强度矢量
plvect,h 补充一点：我想算夹件的涡流损耗，进而计算夹件的温升，但是计算出的涡流损耗只有130多W，肯定有地方不对，不知道是边界条件没有定义对还是别的什么原因，希望帮忙看看，谢谢了 我看了一下，发现电流密度施加的不对 不知道怎么改 每一相的线圈都用线面的命令建模 其中12可以改成13、14
local,12,1,0,0,0
wpcsys,,12 我现在正好也在算损耗和温升，看了楼主的帖子收获很大 ，表示感谢一下 呵呵呵 如果计算时间允许的话可以挨个加电流，然后看方向是不是对的。 剖分的粗细对结果影响很大，但是剖的太细，电脑无法计算了 总提示内存不足 电流没有加上去，要对每一相绕组设置各自的单元坐标系，才能保证电流沿圆周流动。
页:
[1]