编了一个三维涡流场程序，大家帮忙看看哪里不对，谢了

binbin831120

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

binbin831120

补充一点：我想算夹件的涡流损耗，进而计算夹件的温升，但是计算出的涡流损耗只有130多W，肯定有地方不对，不知道是边界条件没有定义对还是别的什么原因，希望帮忙看看，谢谢了

binbin831120

我看了一下，发现电流密度施加的不对 不知道怎么改

binbin831120

每一相的线圈都用线面的命令建模   其中12可以改成13、14
local,12,1,0,0,0
wpcsys,,12

changjiang111

我现在正好也在算损耗和温升，看了楼主的帖子收获很大 ，表示感谢一下 呵呵呵

沙隆巴斯

如果计算时间允许的话可以挨个加电流，然后看方向是不是对的。

binbin831120

剖分的粗细对结果影响很大，但是剖的太细，电脑无法计算了 总提示内存不足

aaa19871988

电流没有加上去，要对每一相绕组设置各自的单元坐标系，才能保证电流沿圆周流动。