新手上路

- 积分
- 37
- 极客币
- 55
- 主题
- 9
- 帖子
- 26
- 注册时间
- 2018-11-5
- 在线时间
- 1 小时
- 性别
- 保密
|
沙发

楼主 |
发表于 2018-11-16 15:45:54
|
只看该作者
由于磁场部分的参数与温度无关,我不想用循环求解,算温度场的时候直接读入对应部分的产热率。那么求解部分的命令流该如何修改?
ALLSEL,ALL
FINISH
/SOLU
ANTYPE,TRANS
TOFFST,273
TUNIF,300
d,nmax,temp,300
CNVTOL,HEAT,1
KBC,1
TRNOPT,FULL
AUTOS,ON
DELTIM,5,5,10,ON
OUTRES,ERASE
OUTRES,NSOL,LAST
PHYSICS,WRITE,THERMAL
SAVE
FINISH
ALLSEL,ALL
/SOL
EQSLV,PCG,1E-8
!!!! 采用循环进行求解,这是间接耦合方法
*DO,I,1,FTIME/TINC ! solution *do loop
*SET,TIME,TIME+TINC ! increment time
PHYSICS,READ,EMAG
/solu !read emag physics file
*IF,I,EQ,1,THEN
TUNIF,300 ! initial temperature
*ELSE
LDREAD,TEMP,LAST,,,,,RST ! read thermal analysis temperatures
*ENDIF
SOLVE ! solve harmonic analysis
FINISH
PHYSICS,READ,THERMAL ! read thermal physics file
/ASSIGN,ESAV,THERM,ESAV ! redirect files for use in thermal restart
/ASSIGN,EMAT,THERM,EMAT
/SOLU
*IF,I,GT,50,THEN
ANTYPE,TRANS,REST ! thermal restart
*ENDIF
TIME,TIME ! time at end of thermal run
ESEL,S,MAT,,2 ! select billet region
LDREAD,HGEN,,,,2,,RST ! apply coupled joule heating load from emag
ESEL,ALL
SOLVE
FINISH
/ASSIGN,ESAV ! reassign files to default
/ASSIGN,EMAT
*ENDDO ! end of solution looping
/post1
set,,,1,,10
PLNSOL,TEMP,,0
!gettemp
/post26
nsel,all
nsol,2,2,temp
plvar,2 |
|