电路-磁场-结构耦合分析中结果不收敛

镜湖

最近在做一个电磁推进的电路-磁场-结构的耦合分析，可是求解时老不能收敛，求楼主和各位朋友给些指点啊！
问题：就是一个电容器给一个外线圈放电，线圈内部的圆筒感应出涡流，我想分析涡流对内部圆筒径向结构上可能造成的影响，Y方向我都固定了。电路如图1，有限元模型如图2.划分网格如图3.程序我列在下面了。弄了好几天了，各位帮帮忙了。

                  图1

                         图2
[localimg=180,147]3[/localimg]
                     图3

程序！
/COM,3504
/PREP7
/TITLE,VZ,WORK_2

ET,1,53,3,,1    !线圈
ET,2,13,,,1     !空气
ET,3,13,4,,1    !抛体
ET,4,110,,,1   !远场单元
EMUNIT,MKS
MP,MURX,1,1
MP,MURX,2,1
MP,MURX,3,1      !磁导率
MP,MURX,9,1
MP,RSVX,1,3E-8
MP,RSVX,3,2.8E-8 !电导率
MP,EX,3,7.6E10
MP,NUXY,3,0.35   !泊松比
R,1,0.025,50,,1,0.9
S=0.07
R=3*S/2
N,1              !电路单元建模
*REPEAT,6,1
RECTNG,0.008,0.02,0.15,0.4    !抛体
RECTNG,0.07,0.12,0,0.25   !线圈
PCIRC,0,15*S,0,90         !远场单元
PCIRC,0,12*S,0,90         !空气

AOVLAP,ALL
NUMCMP,ALL
/PNUM,AREA,1
APLOT
ALLSEL,ALL
!分配材料
ASEL,S,AREA,,1               !抛体
AATT,3,,3
ASEL,S,AREA,,2               !线圈
AATT,1,1,1
ASEL,S,AREA,,4               !空气
AATT,2,1,2
ASEL,S,AREA,,3               !远场
AATT,9,,4
ASEL,ALL
CSYS,1
!划分网格
LSEL,S,LOC,X,13.5*S
LESIZE,ALL,,,1  
ESIZE,,15
AMESH,3
ASEL,ALL
MSHAPE,0,2D
SMRTSIZE,2
MSHKEY,0
ESIZE,S/5
AMESH,2
AMESH,1
SMRTSIZE,2
ESIZE,S/4
MSHAPE,1,2D
MSHKEY,0
AMESH,4
LSEL,S,LINE,,1
NSLL,S
D,ALL,UY,0
LSEL,S,LINE,,3
NSLL,S
D,ALL,UY,0
ESEL,S,MAT,,1
NSLE,S
CP,1,CURR,ALL
CP,2,EMF,ALL
ALLSEL,ALL
CSYS,1
NSEL,S,LOC,X,15*S
SF,ALL,INF
CSYS,0
NSEL,S,LOC,X,0
D,ALL,AZ,0
!电路建立
ET,5,124,2      !电容
ET,6,124,0      !电阻
ET,7,124,1      !电感
ET,8,124,5      !耦合区
ET,9,125        !二极管

R,4,592E-6,500
R,5,0.2
R,6,2E-5,0
R,7,1
R,8

N1=NODE(S,0,0)
TYPE,5
REAL,4
E,2,1
TYPE,6
REAL,5
E,2,3
TYPE,7
REAL,6
E,3,4
TYPE,8
REAL,7
E,4,1,N1
TYPE,9
REAL,8
E,1,2
D,1,VOLT,0
NSEL,ALL
ESEL,ALL
FINISH
allsel,all
!求解
/SOLU
ANTYPE,TRANS
OUTRES,ALL,ALL
TIME,0.001
DELTIM,0.001/1000,0.001/1000,0.001/500
SOLVE
finish

镜湖

       网格划分

镜湖

各位朋友帮帮忙了，实在是弄得头晕了！

whatxshj

把线圈的自由度keyopt(1)=2试下

镜湖

4# whatxshj


谢谢你的提醒，我回到实验室改一下试试。