请教：涡流场+solid97+circu124 一直调不对

biomzll

运行之后总是出现这样的错误：

*** WARNING ***                         CP =      66.753   TIME= 15:13:13
Coefficient ratio exceeds 1.0e8 - Check results.                        
*** WARNING ***                         CP =      71.199   TIME= 15:13:23
Pivoting has been activated for the Sparse Matrix Solver.  The memory   
required may greatly exceed predicted amount.  In that event rerun      
with increased memory using -m option (refer to the Basic Analysis      
Procedures Guide, Ch.  19).                                             
*** ERROR ***                           CP =    9932.880   TIME= 17:24:18
Solution not converged at time 1.E-04 (load step 1 substep 1).         
  Run terminated.   


请各位高手不吝赐教...........   谢谢了。
下面我把命令流贴上了。
!将容积导体放置在激励线圈与检测线圈之间，通有交流电的线圈产生交变磁场，即激励磁场B，穿过导体时,由于电磁感应作用使导
体内产生涡流，该涡流同时会干生出磁场B1，并能改变原激励磁场的强弱。在检测线圈上可以检测到B1.
finish  $ /clear  

/title,transient analysis
/filname,lzq,1
!!!!!!!模型建立第一部分!!!!!!
/prep7
LOCAL,11,1,0,0,0.128,0,0,0,1,1
LOCAL,12,1,0,0,-0.113,0,0,0,1,1

WPCSYS,-1,0
CSYS,0
wprota,,90
cylind,,0.01,-0.04,0.04,,,         !容积导体
cylind,0.01,0.098,-0.04,0.04       !盐溶液

WPCSYS,-1,11
CSYS,11
LOCAL,11,1,0,0,0.128,0,0,0,1,1                           
!3激励线圈壳部分 正中间
cylind,0,0.03884,-0.025,0.025,,,
cylind,0.03884,0.0394,-0.025,0.025,,,
!cylind,0,0.03884,0.103,0.153,,,
!cylind,0.03884,0.0394,0.103,0.153,,,   !4激励线圈线部分 距0.5厘米

WPCSYS,-1,12
CSYS,12
LOCAL,12,1,0,0,-0.113,0,0,0,1,1
!检测线圈壳部分
cylind,0,0.014,-0.01,0.01,,,
cylind,0.014,0.01456,-0.01,0.01,,, !检测线圈线部分

WPCSYS,-1,0
CSYS,0
block,-0.3,0.3,-0.3,0.3,-0.2,0.32    !外部空气

vsel,all
vovlap,all                    !布尔运算
numcmp,volu                   !体积压缩号
!!!!!!!!!!!!!!!!!模型建立完毕!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!定义单元!!!!!!!!!!!!!!!

!vsel,s,volu,,1    !容积导体
!vsel,s,volu,,2    !激励线圈柱
!vsel,s,volu,,3    !检测线圈柱
!vsel,s,volu,,4    !盐溶液
!vsel,s,volu,,5    !激励线圈线
!vsel,s,volu,,6    !检测线圈线
!vsel,s,volu,,7    !空气包

ET,1,SOLID97,1,0       !容积导体
ET,2,SOLID97       !激励柱
ET,3,SOLID97       !检测柱
ET,4,SOLID97,1,0       !盐溶液
ET,5,SOLID97,3
!KEYOPT,5,1,3      !激励线圈线
ET,6,SOLID97,3
!KEYOPT,6,1,3      !检测线圈线（这里用的也是第一关键选项设置为3）
ET,7,SOLID97      !空气包

!!!!!!!!!!!!定义材料属性!!!!!!!!!!!!!!!!!

mp,murx,1,1           !容积导体
mp,rsvx,1,1.75e-7   
mp,murx,2,1
mp,murx,3,1
mp,murx,4,1           !盐溶液
mp,rsvx,4,2.7

mp,murx,5,1
mp,rsvx,5,1.75e-7     !激励线圈线
mp,murx,6,1
mp,rsvx,6,1.75e-7     !检测线圈线

mp,murx,7,1           !空气包

R,1,0.000028,150,0.0000065,0,1,0
R,2,0.0000112,60,0.0000009,0,1,0

vsel,s,,,1
vatt,1,,1,0
vsel,s,,,2
vatt,2,,2,11
vsel,s,,,3
vatt,3,,3,12
vsel,s,,,4
vatt,4,,4,0
vsel,s,,,5
vatt,5,1,5,11
vsel,s,,,6
vatt,6,2,6,12
vsel,s,,,7
vatt,7,,7,0

!!!!!!!!!!!!!!!!!材料属性定义完毕!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!网格划分!!!!!!!!!!!!!!!!!!!!!!!!
!vsel,s,volu,,1    !容积导体
!vsel,s,volu,,2    !激励柱
!vsel,s,volu,,3    !检测柱
!vsel,s,volu,,4    !盐溶液
!vsel,s,volu,,5    !激励线
!vsel,s,volu,,6    !检测线
!vsel,s,volu,,7    !空气包

vsel,s,volu,,2
mshape,1,3d
mshkey,0
smrtsize,4
vmesh,all

vsel,s,volu,,5
mshape,1,3d
mshkey,0
smrtsize,3
vmesh,all

vsel,s,volu,,4      !smrtsize划分
mshape,1,3d
mshkey,0
smrtsize,2
vmesh,all
vsel,s,volu,,1      !smrtsize划分
mshape,1,3d
mshkey,0
smrtsize,5
vmesh,all
vsel,s,volu,,3
mshape,1,3d
mshkey,0
smrtsize,4
vmesh,all

vsel,s,volu,,6
mshape,1,3d
mshkey,0
smrtsize,3
vmesh,all
vsel,s,volu,,7
mshape,1,3d
mshkey,0
vmesh,all




!!!!!!!!!!!!!!!!!!网格划分完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!场路耦合!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!耦合自由度!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

esel,s,mat,,5
nsle,s
cp,1,curr,all
cp,2,emf,all
*get,n1,node,0,num,min     !取最小编号激励
cm,k1,elem  
esel,s,mat,,6
nsle,s
cp,3,curr,all
cp,4,emf,all
*get,n2,node,0,num,min     !取最小编号检测

cm,k2,elem  
allsel,all


!!!!!!!!!!!!!!耦合完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


!!!!!!!!!!!!!场路连接!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!激励源
ALLSEL,ALL
CSYS,0
ET,9,CIRCU124,4,1     !交流电压源
ET,10,CIRCU124,5        !绞线圈  
ET,11,circu124,0       !电阻抗
R,9,,10,1000,,0
R,10,1                  !对称因子等于1
R,11,5                   !电阻5

!激励电压源
N,,0,0.1,0.7              !NMAX-5
N,,0,0,0.7               !NMAX-4
N,,0,-0.1,0.7                !NMAX-3
N,,0,0.1,0.5              !NMAX-2

!检测电阻
N,,0,0.25,-0.5              !NMAX-1
!N,,0,0.15,-0.5              
N,,0,0.05,-0.5               !NMAX

!激励场路耦合
*Get,Nmax,Node,0,Num,Max
Type,9
Real,9
E,Nmax-5,Nmax-3,Nmax-4  

Type,10
Real,10
E,Nmax-3,Nmax-2,n1      

Type,11
Real,11
E,Nmax-5,Nmax-2

D,Nmax-3,volt,0          !零线端接地

!检测端
allsel,all
csys,0
ET,12,circu124,0     !电阻抗
ET,13,CIRCU124,5     !绞线圈

R,12,10               !5欧姆的电阻
R,13,1

Type,12
Real,12
E,Nmax-1,Nmax  

Type,13
Real,13
E,Nmax-1,Nmax,n2
D,Nmax,volt,0


!!!!!!!!!!!!!!!设置边值条件!!!!!!!!!!!!!!!!!!
csys,0
nsel,s,loc,x,-0.3
nsel,a,loc,x,0.3
D,ALL,AY,0
D,ALL,AZ,0
D,ALL,AX,0
nsel,a,loc,y,-0.3
nsel,a,loc,y,0.3
D,ALL,AX,0
D,ALL,AZ,0
D,ALL,AY,0
nsel,a,loc,z,-0.2
nsel,a,loc,z,0.32
D,ALL,AX,0
D,ALL,AY,0
D,ALL,AZ,0
FINISH

!!!!!!!!!!!!!!!!!!!!求瞬态分析!!!!!!!!!!!!!!
/solu
antype,trans
outres,all,all        !控制数据库和结果输出文件
trnopt,full           !选择求解方法试试全波方法
nropt,full,,on         !牛顿拉普逊法
eqslu,sparse       !选择稀疏矩阵求解器
time,0.001      !一个周期
deltime,0.001/10 !设置载荷步长
outres,all,all
allsel,all
solve
save
finish

沙隆巴斯

*** WARNING ***                         CP =      66.753   TIME= 15:13:13
Coefficient ratio exceeds 1.0e8 - Check results.                        
*** WARNING ***                         CP =      71.199   TIME= 15:13:23
Pivoting has been activated for the Sparse Matrix Solver.  The memory   
required may greatly exceed predicted amount.  In that event rerun      
with increased memory using -m option (refer to the Basic Analysis      
Procedures Guide, Ch.  19).                                             
*** ERROR ***                           CP =    9932.880   TIME= 17:24:18
Solution not converged at time 1.E-04 (load step 1 substep 1).         
  Run terminated

针对警告和错误，
第一个，有不合适的剖分单元，这个到还好说
第二个，选择其他的求解器
第三个，这个很致命，求解不收敛，这个是由前面一系列警告带来的最终结果。
建议，重新检查网格剖分，查找不合适参数。

biomzll

谢谢沙隆巴斯的关注和回复。
你的意思是还是剖分部分有问题?嗯，我再改改看吧。
剖分的时候采用体分割剖分的话，体分割之后还需要体粘合吗？
请不吝赐教...

hchangduan

你不粘合就不是一个整体了

沙隆巴斯

剖分完了再粘还管用吗？最好还是整体剖

shagf

这个算例不错，俺也学习一下

ycc1111

学习了，持续学习中