由 fossiler 发表的一个涡流计算的例子,因为不太容易找,故单独列到这里来方便大家查找。
原文见: http://www.simwe.com/forum/thread-797200-2-1.html
- 以team21为例,分别采用自然计算与旋转节点,耦合自由度计算
- 拉板磁导率选为10000
- 损耗计算结果:
- 自然计算(默认情况):4.52386
- 旋转节点耦合自由度:0.139526E-01
- 下是命令流,大家看看是不是有问题
- a1=0.2
- a2=0.27
- r1=0.01
- r2=0.045
- h1=(0.458-0.024)/2
- h2=0.024
- g1=0.012
- b_h=0.82
- b_a=0.01
- b_b=0.36
- in=3000
- /prep7
- et,1,97
- et,2,97,1
- mp,murx,1,1
- mp,murx,2,10000
- mp,rsvx,2,1/(1.3889*1000000)
- rectng,-a1/2,a1/2,-a1/2,a1/2
- rectng,-a2/2,a2/2,-a2/2,a2/2
- aovlap,all
- adele,all
- LFILLT,3,2,r1, ,
- LFILLT,3,4,r1, ,
- LFILLT,4,1,r1, ,
- LFILLT,1,2,r1, ,
- LFILLT,6,7,r2, ,
- LFILLT,7,8,r2, ,
- LFILLT,8,5,r2, ,
- LFILLT,5,6,r2, ,
- x1=a1/2-r1
- y1=a1/2-r1
- wpoff,x1,y1
- CSWPLA,11,1,1,1,
- wpoff,-2*x1
- CSWPLA,12,1,1,1,
- wpoff,,-2*y1
- CSWPLA,13,1,1,1,
- wpoff,2*x1
- CSWPLA,14,1,1,1,
- wpoff,-x1,y1
- csys,0
- LSTR, 14, 22
- LSTR, 13, 21
- LSTR, 15, 23
- LSTR, 16, 24
- LSTR, 10, 17
- LSTR, 9, 18
- LSTR, 11, 19
- LSTR, 12, 20
- FLST,2,4,4
- FITEM,2,10
- FITEM,2,23
- FITEM,2,14
- FITEM,2,24
- AL,P51X
- FLST,2,4,4
- FITEM,2,24
- FITEM,2,8
- FITEM,2,4
- FITEM,2,18
- AL,P51X
- FLST,2,4,4
- FITEM,2,18
- FITEM,2,15
- FITEM,2,17
- FITEM,2,11
- AL,P51X
- FLST,2,4,4
- FITEM,2,1
- FITEM,2,17
- FITEM,2,5
- FITEM,2,19
- AL,P51X
- FLST,2,4,4
- FITEM,2,12
- FITEM,2,20
- FITEM,2,16
- FITEM,2,19
- AL,P51X
- FLST,2,4,4
- FITEM,2,2
- FITEM,2,6
- FITEM,2,20
- FITEM,2,21
- AL,P51X
- FLST,2,4,4
- FITEM,2,9
- FITEM,2,21
- FITEM,2,22
- FITEM,2,13
- AL,P51X
- FLST,2,4,4
- FITEM,2,7
- FITEM,2,3
- FITEM,2,22
- FITEM,2,23
- AL,P51X
- VEXT,all, , ,0,0,h1,,,,
- VGEN,2,all, , , , ,h1+h2, ,0
- wpoff,,,h1+h2/2
- wpoff,-(a2/2+g1+b_a/2),,,
- BLOCK,-b_a/2,b_a/2,-b_b/2,b_b/2,-b_h/2,b_h/2,
- wpoff,(a2/2+g1+b_a/2),,,
- BLOCK,-0.7*a2,0.65*a2,-0.65*b_b ,0.65*b_b ,-0.6*b_h,0.6*b_h ,
- vovlap,all
- numcmp,all
- asel,s,ext
- vsla,s
- cm,kq,volu
- allsel,all
- cmsel,u,kq
- lsel,s,length,,b_h
- asll,s
- vsla,r
- cm,jb,volu
- allsel,all
- cmsel,u,kq
- cmsel,u,jb
- cm,xq,volu
- vsel,r,loc,y,a1/2,a2/2
- vsel,r,loc,x,-a1/2,a1/2
- cm,xq_x1,volu
- cmsel,s,xq
- vsel,r,loc,y,-a1/2,-a2/2
- vsel,r,loc,x,-a1/2,a1/2
- cm,xq_x2,volu
- cmsel,s,xq
- vsel,r,loc,x,-a1/2,-a2/2
- vsel,r,loc,y,-a1/2,a1/2
- cm,xq_y1,volu
- cmsel,s,xq
- vsel,r,loc,x,a1/2,a2/2
- vsel,r,loc,y,-a1/2,a1/2
- cm,xq_y2,volu
- cmsel,s,xq
- cmsel,u,xq_x1
- cmsel,u,xq_x2
- cmsel,u,xq_y1
- cmsel,u,xq_y2
- cm,xq2,volu
- vsel,r,loc,x,0,10
- vsel,r,loc,y,0,10
- cm,xq_r1,volu
- cmsel,s,xq2
- vsel,r,loc,x,0,-10
- vsel,r,loc,y,0,10
- cm,xq_r2,volu
- cmsel,s,xq2
- vsel,r,loc,x,0,-10
- vsel,r,loc,y,0,-10
- cm,xq_r3,volu
- cmsel,s,xq2
- vsel,r,loc,x,0,10
- vsel,r,loc,y,0,-10
- cm,xq_r4,volu
- cmsel,s,kq
- vatt,1,,1,0
- cmsel,s,xq
- cmsel,u,xq2
- vatt,1,,1,0
- cmsel,s,xq_r1
- vatt,1,,1,11
- cmsel,s,xq_r2
- vatt,1,,1,12
- cmsel,s,xq_r3
- vatt,1,,1,13
- cmsel,s,xq_r4
- vatt,1,,1,14
- cmsel,s,jb
- vatt,2,,2,0
- allsel,all
- lsel,s,length,,b_a
- LESIZE,all, , ,8, , , , ,1
- lsel,s,length,,b_b
- LESIZE,all, , ,20, , , , ,1
- lsel,s,length,,b_h
- LESIZE,all, , ,50, , , , ,1
- cmsel,s,jb
- vsweep,all
- allsel,all
- lsel,s,length,,(a2-a1)/2
- LESIZE,all, , ,3, , , , ,1
- lsel,s,length,,h1
- LESIZE,all, , ,6, , , , ,1
- cmsel,s,xq
- vsweep,all
- MSHAPE,1,3D
- MSHKEY,0
- allsel,all
- cmsel,s,kq
- aslv,s,1
- lsla,s,1
- LESIZE,all,0.03, , , , , , ,1
- vmesh,all
- allsel,all
- !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!(这一部分可选可不选,用来比较)
- cmsel,s,jb
- asel,s,ext
- lsla,s,ext
- nsll,s,1
- cm,mubdry2n,node
- nsla,s,all
- norl,all,,1 !这一条是我加上的,因为相对于下面nora的那一指令
- cmse,u,mubdry2n
- cm,mubdry1n,node
- nora,all,1
- alls
- *get,nmax,node,,num,max
- cmsel,s,jb
- eslv,s,1
- nsle
- cm,scrap_e,elem
- cm,scrap_n,node
- ngen,2,nmax,all
- egen,2,nmax,all
- cmse,s,scrap_n
- cmse,s,scrap_e
- VCLEAR, all
-
- vdele,all,,,1
- alls
- numc,node
- alls ! UNSELECT "FLAT SURFACE" NODES, LEAVING CORNER
- cmse,u,mubdry2n ! AND EDGE NODES. COUPLE AX/Y/Z ?????这一条为什么不是cmse,s,mubdry2n
- cpintf,ax,1e-8
- cpintf,ay,1e-8
- cpintf,az,1e-8
- alls ! UNSELECT EDGE/CORNER NODES, LEAVING ONLY FLAT
- cmse,u,mubdry1n ! SURFACE NODES. COUPLE AY/Z ONLY?????这一条为什么不是cmse,s,mubdry1n
- cpintf,ay,1e-8
- cpintf,az,1e-8
- !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
- allsel,all
- /sol
- csys,4
- ANTYPE,3
- HARFR,50
- allsel,all
- ASEL,S,EXT
- DA,all,ASYM
- s=(a2-a1)*h1/2
- cmsel,s,xq_x1
- vsel,r,loc,z,0,10
- BFV,all,JS,-in*1.4141414/s ,,,,
- cmsel,s,xq_y1
- vsel,r,loc,z,0,10
- BFV,all,JS,,-in*1.4141414/s,,,
- cmsel,s,xq_x2
- vsel,r,loc,z,0,10
- BFV,all,JS,in*1.4141414/s,,,,
- cmsel,s,xq_y2
- vsel,r,loc,z,0,10
- BFV,all,JS,,in*1.4141414/s,,,
- cmsel,s,xq2
- vsel,r,loc,z,0,10
- BFV,all,JS,,in*1.4141414/s,,,
- cmsel,s,xq_x1
- vsel,r,loc,z,0,-10
- BFV,all,JS,in*1.4141414/s ,,,,
- cmsel,s,xq_y1
- vsel,r,loc,z,0,-10
- BFV,all,JS,,in*1.4141414/s,,,
- cmsel,s,xq_x2
- vsel,r,loc,z,0,-10
- BFV,all,JS,-in*1.4141414/s,,,,
- cmsel,s,xq_y2
- vsel,r,loc,z,0,-10
- BFV,all,JS,,-in*1.4141414/s,,,
- cmsel,s,xq2
- vsel,r,loc,z,0,-10
- BFV,all,JS,,-in*1.4141414/s,,,
- allsel,all
- solve
- /post1
- set,last
复制代码 |