正弦激励的瞬态场为什么不收敛？（附命令流）

kt266a

我做了一个简单模型：中心是截面积不变的铜柱，外边包了一层铁。最外边当然是空气包。铜柱通过幅值1A的交流电。

大家不要问为什么不用时谐场，因为我在做同步发电机，这是做一个简单模型试一试三维瞬态。

边界条件：
空气包侧面AX=0，AY=0，AZ=0
上下圆面AX=0，AY=0

按照ansys帮助文件VM186的做法施加电流激励。

我反复调整参数，初始的第一步长T从1E-7    一直改到  1E-9，CURRENT还是不收敛。CUR SEG倒是很快就收敛。

经检查剖分密度满足趋肤深度

另外，求解时出现警告：
there are 1 small equation  solver  pivot  terms. this maybe a  transitional effect caused  by the newton-raphson  logic.


命令流
/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,,,3/1000,3/1000,20/1000
con4,,,8/1000,8/1000,20/1000
cylind,,16/1000,,20/1000,
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,200                !磁导率
mp,rsvx,3,1.12e-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
!LSEL,S,LOC,Z,0,0
!LSEL,S,LOC,x,8,8
!LESIZE, all, 1,
CSYS,1
MSHKEY,0
MSHAPE,1,3D
allsel,all
VMESH,ALL


/SOLU   
ANTYPE,TRANS                  ! NONLINEAR TRANSIENT DYNAMIC ANALYSIS
T=1E-9                                  ! 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,,,10
NSLA,S,1
NSEL,U,LOC,x,3/1000,3/1000
!CP,2,VOLT,ALL
D,all,volt,0
Asel,S,,,9
NSLA,S,1
NSEL,U,LOC,x,3/1000,3/1000
!CP,3,VOLT,ALL
D,all,volt,0
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,16/1000,16/1000
NSLA,S,1
D,all,ax,0
D,all,ay,0
D,all,az,0
ASEL,S,LOC,z,0,0
ASEL,A,LOC,z,20/1000,20/1000
NSLA,S,1
D,all,ax,0
D,all,ay,0
!NEQIT,1
!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(360*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,3
FINISH

kt266a

另外，我搞不懂：
[img]为什么我输入的都是线性材料，还出现了非线性分析的迭代过程？[/[/img]color]

huahudie

按照你的这个步骤能出来集肤效应的效果？？？我看你的步骤跟我的是一样的，只不过我用的是117单元，加的是电流密度载荷，其他的都一样了，但是我的没有集肤效应的效果。
   你的这个不收敛是不是每一个载荷步迭代次数太少了，我看你的好像是一个载荷步迭代了26次，你改大点，改成50试试。我也不太懂，刚做这个，我就是这样猜想的。

余光耀

有高导磁性材料的模型最好不要用97单元，因为高导磁材料模型的边界条件相对于97单元来说比较差

kt266a

改收敛判据