变压器中绕组受力形变问题

js2010

最近仿照Vm172做了一个变压器模型，想观察一下绕组所受应力及形变情况，但是结果不理想，老是有错误提示且结果不收敛，如图所示
命令流如，望大侠指教：
/PREP7  
RECTNG,0,0.515,0,1.768              
RECTNG,0.26,0.515,0.184,1.584      
RECTNG,0.278,0.343,0.304,1.464  !低压绕组
RECTNG,0.385,0.465,0.304,1.464  !高压绕组
AOVLAP,all  
numcmp,area
ET,1,13,4,,1   !低压绕组材料
ET,2,13,4,,1   !高压绕组材料
ET,3,13,,,1    !铁心
ET,4,13,,,1    !空气

MP,EX,1,10.76E10                    ! SOLENOID MODULUS OF ELASTICITY
MP,NUXY,1,.35                       ! SOLENOID POISSON RATIO
MP,EX,2,10.76E10                    ! SOLENOID MODULUS OF ELASTICITY
MP,NUXY,2,.35                       ! SOLENOID POISSON RATIO
MP,MURX,1,1
MP,MURX,2,1
MP,MURX,3,1000
MP,MURX,4,1
asel,s,,,1            !赋予材料属性
aatt,1,,1  
asel,s,,,2  
aatt,2,,2  
asel,s,,,3  
aatt,3,,3   
asel,s,,,4  
aatt,4,,4   
allsel,all  
lsel,s,,,10,16,2                !高低压绕组划分
LESIZE,all, , ,30, , , , ,1
lsel,s,,,9,15,2
LESIZE,all, , ,1, , , , ,1  
asel,s,,,1,2
MSHAPE,0,2D
MSHKEY,1
amesh,all   
aSEL,INVE                       !空气和铁心划分  
SMRT,2  
MSHAPE,1,2D
MSHKEY,0
amesh,all
  
allsel,all  
lsel,s,,,1,3,2                  
DL,all, ,ASYM                    

ESEL,S,MAT,,1
NSLE,S
NSEL,R,LOC,x,0.278
CP,1,Ux,ALL                         ! COUPLE SOLENOID NODAL Ux DISP.
allsel
ESEL,S,MAT,,2
NSLE,S
NSEL,r,LOC,X,0.465                      ! SELECT NODES AT OUTER RADIUS
CP,2,ux,ALL
allsel
      
FINISH
/SOLU
asel,s,,,1  
bfa,all,js,,,-1438615.21684  
allsel,all  
asel,s,,,2  
bfa,all,js,,,1245831.1568   
allsel,all  
NSEL,ALL
ESEL,ALL
KBC,1                               ! STEP BOUNDARY CONDITIONS
OUTRES,,LAST
SOLVE
FINISH