关闭瞬态效应仿真永磁体充磁过程

崔佳乐

选择瞬态分析后关闭瞬态效应仿真永磁体，为什么会出现下面的这种情况呢？永磁体下铝板上的电流为0，但铝板所受洛伦兹力却不是0，命令流如下：
finish

/clear  
/filename,EMAT_3D_A0,1                           !定义文件名
/title,EMAT_3D_A0                          !定义工作名称
/CONFIG,FSPLIT,2000 !1=4mb 把.rst文件分割
/PREP7  
!定义单元
ET,1,SOLID117                               !钢板

ET,2,SOLID117
                              !线圈

ET,3,SOLID117                               !空气

ET,4,SOLID117                              !永磁体

KEYOPT,1,1,1
KEYOPT,2,1,1
EMUNIT,MKS

MP,MURX,1,1                       !铝板相对磁导率
MP,RSVX,1,0.283e-7                !铝板电阻率


MP,MURX,2,1                           !线圈相对磁导率
MP,RSVX,2,0.172e-7                     !铝板电阻率

MP,MURX,3,1                      !空气相对磁导率

HC=880000                        !永磁铁矫顽力
MP,MURX,4,5.30504
mp,mgxx,4,0
mp,mgyy,4,HC                    !定义Y方向为磁化方向
mp,mgzz,4,0

BLOCK,-0.01075,-0.01025,,0.00025,-0.011,0.011           !!!线圈模型
VGEN,7,ALL,,,0.0035,,,,0               !复制7个体
VSEL,NONE

BLOCK,-0.01075,-0.00675,,0.00025,-0.011,-0.01150
VGEN,3,ALL,,,0.0070,,,,0
VSEL,NONE

BLOCK,-0.00725,-0.00325,,0.00025,0.011,0.01150
VGEN,3,ALL,,,0.0070,,,,0
VSEL,NONE

BLOCK,-0.00675,-0.00375,,0.00025,-0.011,-0.01150
VGEN,3,ALL,,,0.0070,,,,0
VSEL,NONE

BLOCK,-0.00725,-0.01025,,0.00025,0.011,0.01150
VGEN,3,ALL,,,0.0070,,,,0
VSEL,NONE

BLOCK,-0.01025,-0.00725,,0.00025,-0.011,0.011
VGEN,6,ALL,,,0.0035,,,,0

BLOCK,-0.01075,-0.01025,,0.00025,0.011,0.01150

BLOCK,0.01075,0.01025,,0.00025,-0.011,-0.01150


ALLSEL
VGLUE,ALL
NUMCMP,VOLU
!磁铁
BLOCK,-0.012,0.012,0.00125,0.01625,-0.015,0.015
!铝板模型
BLOCK,-0.022,0.662,-0.00075,-0.00225,-0.025,0.025

!空气模型
BLOCK,-0.022,0.662,-0.00075,0,-0.025,0.025
BLOCK,-0.012,0.012,0,0.00125,-0.015,0.015
BLOCK,-0.0700,0.7000,-0.0500,0.0500,-0.0500,0.0500  !远场空气模型

VOVLAP,ALL                                  !布尔搭接操作
NUMCMP,VOLU

VSEL,S,,,32  !远场空气模型
VSEL,INVE                              
VSEL,U,,,28   !磁铁
VSEL,U,,,31   !导线上空气

!为实体分配单元属性
!线圈

VSEL,S,,,1,8,1
VSEL,A,,,15
VSEL,A,,,22,27,1
VATT,2,,2,0  

!永磁体
VSEL,S,,,28
VATT,4,,4,0  

!铝板
VSEL,S,,,29
VATT,1,,1,0
!空气
VSEL,S,MAT,,1                          
VSEL,A,MAT,,2   
VSEL,A,MAT,,4         
VSEL,INVE
VATT,3,,3,0
!网格划分
!铝板
LSEL,S,,,115
LESIZE,ALL,,,440,,,,,0
LSEL,S,,,119
LESIZE,ALL,,,30,,,,,0
LSEL,S,,,114
LESIZE,ALL,,,4,,,,,0
ALLSEL
MSHAPE,0,3D
MSHK,1
VMESH,29

VSEL,S,,,1,28,1
VSEL,A,,,31
MSHAPE,1,3D
MSHK,0
VMESH,ALL         

VSEL,S,,,30
SMRT,5
VMESH,ALL

VSEL,S,,,32
SMRT,5
VMESH,ALL
!加电压耦合
ASEL,S,,,167         !选择线圈输入端端面
NSLA,S,1        
CP,1,VOLT,ALL        !耦合输入端端面电压
ASEL,S,,,196
NSLA,S,1
D,ALL,VOLT,0         !线圈另一面（输出端）指定电压为0
!*
ALLSEL

allsel,all
finish

!求解  
/SOL
NCNV,2
allsel,all  
ANTYPE,4
TRNOPT,FULL   !指定载荷步中所需要的子步数
TIMINT,OFF
TIME,1E-5
ASEL,S,,,72,73,1              !选择所有外表面
ASEL,A,,,75,78,1                  
DA,ALL,ASYM
NSUBST,7,,
AUTOTS,ON
KBC,0      !
FK,205,AMPS,0
allsel,all  
/STATUS,SOLU
OUTRES,ALL,ALL
SOLVE   
SAVE                         !给线圈耦合面上(输入端)的一点施加总电         
SAVE,EMAT_3D_A0_resu,d