运行之后总是出现这样的错误:
*** 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 |