小弟在用ANSYS分析温度场时不同载荷步下温度差别很大,而且载荷步子步数不一样,这是什么原因。。。
这是我的命令流:
/SOLU
/uis,msgpop,3
/NERR,5,10000000
!DA,9,UX,0
!DA,9,UY,0
!DA,9,UZ,0
!DA,10,UX,0
!DA,10,UY,0
!DA,10,UZ,0
!DA,11,UX,0
!DA,11,UY,0
!DA,11,UZ,0
!DA,12,UX,0
!DA,12,UY,0
!DA,12,UZ,0
ESEL,all
*GET,NEMAX,ELEM,,NUM,MAX
*GET,NEMIN,ELEM,,NUM,MIN
T=0
DT0=1e-6
V=8e-3
L=100e-3
Q=90/8
h1=0.7e-3
LNUM=(5-1)/2
LNUM1=32
r2=r0+h1
TINC=2*PI*r2/V
w=2*PI/TINC
C=V*360/2/PI/r2
DT1=2*PI*r2/V/LNUM1
C1=360/LNUM1
dt2=0.4*DT1
TINC1=2*PI/(LNUM1*w)
ALLS
ANTYPE,TRANS
TRNOPT,FULL
nropt,full,,on
!NLGEOM,ON
PRED,ON
TIME,DT0
TIMINT,OFF
KBC,1
ESEL,ALL
NSLE,S
IC,ALL,TEMP,25
ALLS
!OUTRES,NSOL,ALL
SOLVE
TIMINT,ON
TINTP,0.005,,,1,0.5,0.2
TREF,25
lnsrch,on
cnvtol,heat,,0.03
ncnv,2
outpr,nsol,all
outres,nsol,all
CSYS,1
ESEL,S,TYPE,,2
!nsel,r,loc,x,r11+0.0002,r2
!nsel,r,loc,y,C1*D,C1*D+C1
!nsel,r,loc,z,0,l1
!ESLN,r,1
ekill,all
esel,s,live
eplot
allsel,all
!outres,all,all
ESEL,all
*DO,M,0,LNUM,1
!mj=MOD(M,2)
!*if,mj,eq,1
*afun,deg
!M=0
M1=0.001+0.001*M*2
!*DO,D,0,LNUM1-1,1
*DO,D,0,31,1
CSYS,1
ESEL,S,TYPE,,2
nsel,r,loc,x,r11,r2
nsel,r,loc,y,C1*D,C1*D+C1
nsel,r,loc,z,M1-0.001,M1+0.001
esln,s,0
ealive,all
esel,s,live
eplot
allsel
!outres,all,all
CSYS,0
x=r2*cos(C1*D)
y=r2*sin(C1*D)
z=M1
*DO,f,NEMIN,NEMAX,1
*IF,ESEL(f),EQ,1,THEN
XSY=CENTRX(f)
YSY=CENTRY(f)
ZSY=CENTRZ(f)
RR=ABS(SQRT((XSY-x)*(XSY-x)+(YSY-y)*(YSY-y)+(ZSY-z)*(ZSY-z)))
*IF,RR,LE,RADIUS1,THEN
QR=QMAX*EXP(-k*RR*RR)
SFE,f,NMFACE(f),HFLUX,,QR
*ENDIF
*ENDIF
*ENDDO
TIMINT,ON
AUTOTS,ON
KBC,1
DELTIM,0.05,0.01,0.2
TM1=dt2+(M-0)*TINC+(D-0)*TINC1
TIME,TM1
allsel,all
!outres,all,all
SOLVE
ESEL,S,live
nsel,r,loc,z,M1-0.001,M1+0.001
esln,s,0
*DO,f,1,6
SFEDELE,ALL,f,HFLUX
*ENDDO
TM=DT1+(M-0)*TINC+(D-0)*TINC1
TIME,TM
esel,s,live
eplot
alls
ASEL,S,TYPE,,4
SFA,ALL,1,CONV,100,25
ALLSEL,ALL
!outres,all,all
SOLVE
*ENDDO
*ENDDO
finish
![ANSYS热分析温度差别较大]()
载荷步45.png
![ANSYS热分析温度差别较大-1]()
在和不46.png |