finish /clear
/filename,EMAT_3D_A0,1 !定义文件名 /title,EMAT_3D_A0 !定义工作名称 /PREP7 !定义单元
ET,1,SOLID117 !钢板
ET,2,SOLID117
!线圈
ET,3,SOLID117 !空气
ET,4,SOLID117 !永磁体 ET,5,SOLID117
ET,6,SOLID117 KEYOPT,1,1,1 KEYOPT,2,1,1 EMUNIT,MKS
MP,MURX,2,1 !线圈
MP,RSVX,2,0.175e-7
MP,MURX,1,1
MP,RSVX,1,0.283e-7 !电阻率 MP,MURX,3,1 !空气
HC=895000
MP,MURX,4,11 MP,MURX,5,11
MP,MURX,6,1000
MP,MGXX,4,0 MP,MGYY,4,HC MP,MGZZ,4,0
MP,MGXX,5,0 MP,MGYY,5,-HC MP,MGZZ,5,0
BLOCK,,0.0005,,0.0005,,0.0355 !!!线圈模型 VGEN,12,ALL,,,0.0030,,,,0
VSEL,NONE
BLOCK,,0.0035,,0.0005,-0.0005,0 VGEN,6,ALL,,,0.0060,,,,0 VSEL,NONE
BLOCK,0.0030,0.0065,,0.0005,0.0355,0.0360 VGEN,5,ALL,,,0.0060,,,,0 VSEL,NONE
BLOCK,0.0035,0.0060,,0.0005,-0.0005,0 VGEN,5,ALL,,,0.0060,,,,0 VSEL,NONE
BLOCK,,0.0005,,0.0005,0.0355,0.0360 VGEN,2,ALL,,,0.0330,,,,0
VSEL,NONE
BLOCK,0.0005,0.0030,,0.0005,,0.0355 VGEN,11,ALL,,,0.0030,,,,0 VSEL,NONE
BLOCK,0.0005,0.0030,,0.0005,0.0355,0.0360 VGEN,6,ALL,,,0.0060,,,,0
ALLSEL VGLUE,ALL NUMCMP,VOLU
!磁铁
BLOCK,-0.0450,-0.0050,,0.0150,-0.0070,0.0430 BLOCK,0.0385,0.0785,,0.0150,-0.0070,0.0430
BLOCK,-0.0450,0.0785,0.0150,0.0350,-0.0070,0.0430
!钢板模型
BLOCK,-0.0600,0.0935,-0.0010,-0.0310,-0.0240,0.0600
!空气模型
BLOCK,-0.0050,0.0385,,0.0150,-0.0070,0.0430 !线圈上空气 BLOCK,-0.0600,0.0935,-0.0010,0,-0.0240,0.0600!钢板上空气 BLOCK,-0.25,0.25,-0.25,0.25,-0.25,0.25
VOVLAP,ALL !布尔搭接操作 NUMCMP,VOLU
VSEL,S,,,54
VSEL,INVE VSEL,U,,,53 VSEL,U,,,51
!为实体分配单元属性 !线圈
VSEL,S,,,1,13,1 VSEL,A,,,25
VSEL,A,,,37,47,1 VATT,2,,2,0
!永磁体 VSEL,S,,,48 VATT,4,,4,0 VSEL,S,,,49 VATT,5,,5,0 VSEL,S,,,51 VATT,6,,6,0 !钢板 VSEL,S,,,50 VATT,1,,1,0 !空气
VSEL,S,MAT,,1 VSEL,A,MAT,,2 VSEL,A,MAT,,4 VSEL,A,MAT,,5 VSEL,A,MAT,,6 VSEL,INVE VATT,3,,3,0 !网格划分 !钢板
LSEL,S,,,211 LSEL,A,,,215
LESIZE,ALL,,,40,,,,,0 LSEL,S,,,210
LESIZE,ALL,,,10,,,,,0 ALLSEL
MSHAPE,0,3D MSHK,1 VMESH,50
VSEL,S,,,1,49,1 VSEL,A,,,51,53,1 MSHAPE,1,3D MSHK,0
VMESH,ALL
VSEL,S,,,54 SMRT,4 VMESH,ALL
!加电压耦合 ASEL,S,,,170 NSLA,S,1
CP,1,VOLT,ALL
!求解
FINISH
/SOL !进入求解器 ALLSEL
ANTYPE,TRANS OUTRES,ALL,, TRNOPT,FULL ALLSEL
TM_START=1E-6 TM_END=1.25E-2
TM_INCR=0.069995E-2 PHI=0
AUTOTS,ON
*DO,TM,TM_START,TM_END,TM_INCR TIME,TM PHI=PHI+1 I=4*COS(PHI) AS=6e-7 JS=I/AS
allsel !选择所有 !*
ANTYPE,4 !选择瞬态分析 LUMPM,0 !*
ASEL,S,EXT !选择所有外表面
DA,all,AZ,0,0 !指定电磁为0边界条件
DA,176,VOLT,0 !线圈另一面指定电压为0
FK,230,AMPS,JS !给线圈耦合面上的一点施加总电流
allsel
NSUBST,2,0,0 !每步指定2子步 OUTRES,ERASE
OUTRES,ALL,1 !输出所有
AUTOTS,0 !关闭自动荷载步 solve !求解
SAVE,EMAT_3D_A0_resu.db
*enddo FINISH
结束循环,与*do结合使用 !