求解出现自由度的问题

崔佳乐

finish  
/clear  
/filename,EMAT_3D_A0,1                           !定义文件名
/title,EMAT_3D_A0                          !定义工作名称
/PREP7  
!定义单元
ET,1,SOLID97,1                               ！钢板

ET,2,SOLID117
                              ！线圈

ET,3,SOLID117                               ！空气

ET,4,SOLID97                              ！永磁体
ET,5,SOLID97
ET,6,SOLID97                           

MP,MURX,2,,1    !线圈
MP,RSVX,2,,0.175e-7   

TB,BH,1,,25  
TBPT,,100,0.46512
TBPT,,200,0.72993
TBPT,,300,0.90090
TBPT,,400,1.0204
TBPT,,500,1.1086
TBPT,,600,1.1765
TBPT,,700,1.2302
TBPT,,800,1.2739
TBPT,,900,1.3100
TBPT,,1000,1.3405
TBPT,,1400,1.4257
TBPT,,1800,1.4778
TBPT,,2200,1.5131
TBPT,,2600,1.5385
TBPT,,3000,1.5576
TBPT,,3400,1.5726
TBPT,,3800,1.5847
TBPT,,4200,1.5945
TBPT,,4600,1.6028
TBPT,,5000,1.6098
TBPT,,7000,1.6332
TBPT,,9000,1.6465
TBPT,,11000,1.6551
TBPT,,13000,1.6611
TBPT,,15000,1.6656
MP,RSVX,1,,1/5.0e6                !（电阻率）
MP,MURX,3,,1                      !空气

HC=895000

TB,BH,4,,30
TBPT,,130,0.10
TBPT,,170,0.20
TBPT,,197,0.30
TBPT,,218,0.40
TBPT,,250,0.50
TBPT,,290,0.60
TBPT,,338,0.70
TBPT,,400,0.80
TBPT,,472,0.90
TBPT,,570,1.00
TBPT,,682,1.10
TBPT,,810,1.20
TBPT,,975,1.30
TBPT,,1600,1.40
TBPT,,2520,1.50
TBPT,,3520,1.60
TBPT,,12000,1.90
TBPT,,17000,2.00
TBPT,,23000,2.10
TBPT,,32000,2.20
TBPT,,42500,2.30
TBPT,,44500,2.32
TBPT,,48200,2.35
TBPT,,52200,2.37
TBPLOT,BH,4,,,
TB,BH,5,,30
TBPT,,130,0.10
TBPT,,170,0.20
TBPT,,197,0.30
TBPT,,218,0.40
TBPT,,250,0.50
TBPT,,290,0.60
TBPT,,338,0.70
TBPT,,400,0.80
TBPT,,472,0.90
TBPT,,570,1.00
TBPT,,682,1.10
TBPT,,810,1.20
TBPT,,975,1.30
TBPT,,1600,1.40
TBPT,,2520,1.50
TBPT,,3520,1.60
TBPT,,12000,1.90
TBPT,,17000,2.00
TBPT,,23000,2.10
TBPT,,32000,2.20
TBPT,,42500,2.30
TBPT,,44500,2.32
TBPT,,48200,2.35
TBPT,,52200,2.37
TBPLOT,BH,5,,,
MP,MURX,6,,5.30504


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
ALLSEL
VGLUE,ALL

NUMCMP,VOLU
VSEL,S,,,1,23,1
VATT,2,,2,0  
!磁铁
BLOCK,-0.0450,-0.0050,,0.0150,-0.0070,0.0430   
VSEL,S,,,24
VATT,4,,4,0  
BLOCK,0.0385,0.0785,,0.0150,-0.0070,0.0430
VSEL,S,,,25
VATT,5,,5,0
BLOCK,-0.0450,0.0785,0.0150,0.0350,-0.0070,0.0430
VSEL,S,,,26
VATT,6,,6,0

BLOCK,-0.0600,0.0935,-0.0010,-0.0310,-0.0240,0.0600！钢板模型
VSEL,S,,,27
VATT,1,,1,0
ALLSEL
！空气模型
BLOCK,-0.25,0.25,-0.25,0.25,-0.25,0.25
VSEL,S,,,28
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,27


VSEL,S,,,1,26,1           !!!线圈，永磁体划分
MSHAPE,1,3D
MSHK,0
VMESH,ALL
                          ！！！空气划分
SMRT,4
VSEL,S,,,28
MSHAPE,1,3D         
MSHK,0
VMESH,ALL！！！

                          !加电压耦合
ASEL,S,,,2
NSLA,S,1
CP,1,VOLT,ALL



!求解
FINISH  
/SOL                      !进入求解器
ALLSEL
*do,i,1,34,1               !控制循环命令1-34，间隔1，循环34次
time,i*(1e-6)           !求解时间

a=MOD(i,2)             !奇偶函数，奇数为1，偶数为0     
curr=100*a             !电流
allsel                 !选择所有
!*  
ANTYPE,4               !选择瞬态分析
!*  
TRNOPT,FULL            !用full求解方法
!*  
LUMPM,0
!*
ASEL,S,EXT             !选择所有外表面
DA,all,AZ,0,0          !指定电磁为0边界条件


DA,68,VOLT,0         !线圈另一面指定电压为0

FK,6,AMPS,curr      !给线圈耦合面上的一点施加总电流

allsel
NSUBST,2,0,0           !每步指定2子步
OUTRES,ERASE
OUTRES,ALL,1           !输出所有
AUTOTS,0               !关闭自动荷载步
solve                  !求解
*enddo                 !结束循环，与*do结合使用
请问一下为什么求解时会出现下面的问题？