急需帮忙，ansys高斯分布温度场求解时总是自动退出
本帖最后由 desire123 于 2011-7-8 10:56 编辑各位高手们大家好，我现在在做一个高斯分布温度场的仿真，在求解时总是会自动退出，检查了很多遍也没有找出问题，希望大神们能帮小弟一把，不胜感激
finish
/clear,START
/FILNAME,3D THERMAL,1
/prep7
/units,si
R0=0.008
H=0.002
R=0.002
PI=3.1415926
P=10
ABSORD=0.96
et,1,mesh200
keyopt,1,1,6
ET,2,SOLID70
KEYOPT,2,2,1
ET,3,SURF152
R,1
keyopt,3,4,1
keyopt,3,5,0
keyopt,3,8,3
MAT,1
MPTEMP,1,20,50,100,160,200,260,300
MP,KXX,1,2.3
MPDATA,C,1,,0.83,8.17,23.52,36.40,42.68,50.29,54.76
MP,DENS,1,3.14
K,1,0,0,0
K,2,0.008,0,0
K,3,0.008,,0.002
K,4,0,0,0.002
A,1,2,3,4
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
APLOT
LSEL,S,LINE,,1,3,2
LESIZE,ALL,,,8
ALLSEL
LSEL,S,LINE,,2,4,2
LESIZE,ALL,,,10
ALLSEL
AATT,1,1,1
AMESN,ALL
NSEL,S,LOC,Z,0.002
NSEL,R,LOC,X,0
nrefine,all,,,3,1,off,off
TYPE,2
ESIZE,,9
VROTAT, ALL,,,,,,1,4,180, NSEG
NUMCMP,ALL
NUMMRG,ALL,1e-6, , ,LOW
eplot
TYPE,3
NSEL,S,LOC,Z,0.002
ESURF
EPLOT
FINISH
!高斯激光分布的表函数
Qm=P*ABSORD/(PI*R**2)
*DIM,LASTER,TABLE,5,19,1,X,Y,,1
*DO,I,1,5,1
LASTER(I,0)=(I-1)*0.00025
*ENDDO
*DO,J,1,19,1
LASTER(0,J)=J-1
*ENDDO
*DO,I,1,5,1
*DO,J,1,19,1
XCENTER=0
YCENTER=0
DISTANCE=(I-1)*0.00025
*IF,DISTANCE,LE,R,THEN
LASTER(I,J)=Qm*exp(-DISTANCE**2/(R**2))
*ELSE
LASTER(I,J)=0
*ENDIF
*ENDDO
*ENDDO
FINISH
/solu
antype,trans
trnopt,full,,,,,
tunif,20
tref,20
kbc,1
timint,on !打开瞬态效应timint,key,lab
autots,on
TIME,2
nsubst,10,10,1
outres,all,all
esel,s,type,,3
SFE,ALL,5,CONV,,20,20
ALLSEL
EPLOT
NSEL,S,LOC,Z,0
D,ALL,TEMP,20
ALLSEL
CSYS,1
NSEL,S,LOC,X,0.008
SF,ALL,CONV,20,20
ALLSEL
CSYS,0
NSEL,S,LOC,Z,0.002
SF,ALL,HFLUX,%LASTER%
SOLVE
FINISH
solve之前加个allsel就不会自动退出了 2# woomii 谢谢自动退出已经解决，但是求解结果不对，我是打算在圆形中心出加载高斯激光，但是求解结果是如附件中所示结果，能否帮小弟看看是哪里有问题了，万分感激啊 我曾经也遇到过这样的问题，后来发现跟mesh size，还有加载的face no有关系
页:
[1]