请问solid97加amps载荷后为什么不收敛（附命令流）？

kt266a

模型很简单：一根载有交流电的铜棒，一头粗一头细；外层包有一定厚度的铁皮，考虑涡流。最外面是空气包。
——————————————————————————————————————
请问solid97加amps载荷后，如果NEQIT,1    为什么只算出电密没算磁场？

但是，改成NEQIT,5以后，确实出现H，B的结果。但是，似乎不均匀，与模型、载荷的对称性不相符。
——————————————————————————————————————

上面的问题查清楚了。是因为我把设置平衡迭代次数的指令错记为子步数的指令！

下面还有很多问题，请各位接着看下去。

——————————————————————————————————————

我目前计算一台同步发电机遇到困难，所以，先在一些简化的小模型上试验操作。

因为发电机是转子加直流，定子加交流，不能用谐波分析。

以后，大家回帖时就不要问我为什么要用瞬态场算正弦激励问题。

建模、剖分的命令流：

/filname,test,1
finish
!/uis控制GUI操作行为msgpop控制那些从ansys错误子程序中来的信息在消息对话框中显示3只有错误现实出来
/uis,msgpop,3
!clear清除数据标识
/clear
!/prep7进入模型创建预处理器
/prep7
!*afun在参数表达式中，为角度函数指定单位。rad输入输出中使用弧度单位deg使用角度单位stat显示该命令当前设置
*afun,deg
btol,defa
con4,,,2,5,30
con4,,,3,6,30
cylind,,10,,30,
allsel,all
vovlap,all
ALLSEL,ALL
NUMMRG,ALL
NUMCMP,ALL
!定义材料参数
mp,murx,1,1                   !空气导磁率
mp,murx,2,1                   !铜的磁导率
mp,rsvx,2,0.0434e-6      !铜的电阻率
mp,murx,3,1000             !虚构导磁材料磁导率
mp,rsvx,3,0.0434e-6      !电阻率
!单元类型
!涡流区
et,11,solid97
keyopt,11,1,1
!非涡流区
et,21,solid97            !air
keyopt,21,1,0
et,31,solid97            !dc97
keyopt,31,1,5
!ASEL,S,LOC,Z,0,0

vsel,S,,,1
vatt,2,,31
vsel,S,,,2
vatt,3,,11
vsel,S,,,3
vatt,1,,21
smrtsize,1
allsel,all
MSHKEY,0
MSHAPE,1,3D
VMESH,ALL

kt266a

已经知道为什么求解不出磁场了。
我把平衡迭代次数和子步数的指令记错了。呵呵……

现在是这个
/SOLU   
ANTYPE,TRANS                  ! NONLINEAR TRANSIENT DYNAMIC ANALYSIS
T=1E-8                                  ! INITIALIZE TIME PARAMETER
C=0                                       ! INITIALIZE COUNTER PARAMETER
N=3600                                     ! NUMBER OF TIME INCREMENTS PER TURN


CSYS,1

Asel,S,,,1
NSLA,S,0
*get,NA1,node,0,num,min


Asel,S,,,1
NSLA,S,1
CP,1,VOLT,ALL

Asel,S,,,2
NSLA,S,1
D,all,volt,0

vsel,S,,,2
NSLv,S,0
*get,NA2,node,0,num,min

Nsel,S,,,NA2
D,all,volt,0

ASEL,S,LOC,x,20,20
NSLA,S,1
D,all,az,0


!NEQIT,50
NCNV,0
PI=3.14159265                       ! VALUE OF PI
CON=0.02/N                         ! SET TIME INCREMENT


*DO,I,1,3                               ! REPEAT MACRO EXECUTION

TIME,T
  IA=SIN(2*PI*50*T)            ! CALCULATE CURRENT

  Nsel,S,,,NA1
  F,NA1,AMPS,IA                   ! APPLY CURRENT TO A NODE IN CONDUCTOR
  NEQIT,26                              
  NSUBST,5                                            !  ITERATION PER TIME STEP
  KBC,0
  OUTRES,ALL,ALL
  LSWRITE,I

  T=T+CON                            ! INCREMENT TIME
  C=C+1                                  ! INCREMENT COUNTER

*ENDDO


allsel,all
LSSOLVE,1,3
FINISH

kt266a

既然开了这个帖子，希望高手进来讨论一下：
哪些参数明显影响计算精度？剖分网格尺寸，平衡迭代次数，子步数？
应该怎么设置呢？



NSUBST,5
还是出错

现改成AUTOTS,ON

求解中还是有不收敛。另外，H、B图像较正常，JT、JS电密好像不对劲，“铁筒”的JHEA发热密度根本不像样，简直体现不出模型和载荷的几何对称性。

kt266a

/SOLU   
ANTYPE,TRANS                  ! NONLINEAR TRANSIENT DYNAMIC ANALYSIS
T=1E-8                                  ! INITIALIZE TIME PARAMETER
C=0                                       ! INITIALIZE COUNTER PARAMETER
N=360                                     ! NUMBER OF TIME INCREMENTS PER TURN


CSYS,1

Asel,S,,,1
NSLA,S,0
*get,NA1,node,0,num,min


Asel,S,,,1
NSLA,S,1
CP,1,VOLT,ALL

Asel,S,,,2
NSLA,S,1
D,all,volt,0


allsel,all
vsel,S,,,2
NSLv,S,1
*get,NA2,node,0,num,min

Nsel,S,,,NA2
D,all,volt,0

ASEL,S,LOC,x,20,20
NSLA,S,1
D,all,az,0


!NEQIT,50
NCNV,0
PI=3.14159265                       ! VALUE OF PI
CON=0.02/N                         ! SET TIME INCREMENT


*DO,I,1,5                               ! REPEAT MACRO EXECUTION

TIME,T
  IA=SIN(2*PI*50*T)            ! CALCULATE CURRENT

  Nsel,S,,,NA1
  F,NA1,AMPS,IA                   ! APPLY CURRENT TO A NODE IN CONDUCTOR
  !NEQIT,26                              
  NSUBST,1                                            !  ITERATION PER TIME STEP
  !AUTOTS,ON
  KBC,0
  OUTRES,ALL,ALL
  LSWRITE,I

  T=T+CON                            ! INCREMENT TIME
  C=C+1                                  ! INCREMENT COUNTER

*ENDDO


allsel,all
LSSOLVE,1,5
FINISH

正在运行中，还是迭代到26次仍不收敛

我想不明白：我这个模型没有非线性材料，ansys为什么一个劲迭代?

kt266a

附图如下，可以看到不收敛:

perch88

怀疑是循环的问题
把循环去掉，直接复制粘贴，反正只有5次
对了以后再加循环