24小时热门版块排行榜

返回列表

帅帅小城

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 146
帖子: 2
在线: 2.1小时
虫号: 2832131
注册: 2013-11-27
性别: GG
专业: 金属材料的合金相、相变及

[交流] 请大神帮我检查一下ansys高斯表面移动热源名命令流，不知为什么经常出现错误。

/PREP7
ET,1,SHELL57
R,1,0.0001
TOFFST,273
MPTEMP,1,0
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,800
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1800
MPTEMP,8,2000
MPTEMP,9,2800
MPDATA,KXX,1,,15
MPDATA,KXX,1,,17
MPDATA,KXX,1,,21
MPDATA,KXX,1,,25
MPDATA,KXX,1,,28
MPDATA,KXX,1,,32
MPDATA,KXX,1,,43
MPDATA,KXX,1,,70
MPDATA,KXX,1,,84
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,800
MPTEMP,5,1000
MPTEMP,6,1350
MPTEMP,7,1400
MPTEMP,8,1500
MPTEMP,9,2000
MPTEMP,10,3000
MPDATA,C,1,,462
MPDATA,C,1,,512
MPDATA,C,1,,571
MPDATA,C,1,,630
MPDATA,C,1,,680
MPDATA,C,1,,760
MPDATA,C,1,,1600
MPDATA,C,1,,1600
MPDATA,C,1,,760
MPDATA,C,1,,680
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,500
MPTEMP,3,1200
MPTEMP,4,1500
MPTEMP,5,2000
MPTEMP,6,3000
MPDATA,DENS,1,,7960
MPDATA,DENS,1,,7760
MPDATA,DENS,1,,7430
MPDATA,DENS,1,,7320
MPDATA,DENS,1,,6800
MPDATA,DENS,1,,6000
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,HF,1,,18                   !表面传热系数
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,200
MPTEMP,3,500
MPTEMP,4,800
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,3000
MPDATA,EMIS,1,,0.2                !辐射率
MPDATA,EMIS,1,,0.45
MPDATA,EMIS,1,,0.54
MPDATA,EMIS,1,,0.58
MPDATA,EMIS,1,,0.59
MPDATA,EMIS,1,,0.6
MPDATA,EMIS,1,,0.62
RECTNG,0,0.045,0,0.0008
RECTNG,0,0.045,0.0008,0.0018
RECTNG,0,0.045,0.0018,0.005
RECTNG,0,0.045,0.005,0.045
aglue,all
numcmp,all
LESIZE,12,,,20
LESIZE,13,,,20
LESIZE,7,,,20
MSHKEY,1
MSHAPE,0,2D
AMESH,4
LESIZE,10,,,2
LESIZE,11,,,2
MSHKEY,1
MSHAPE,0,2D
AMESH,3
LREFINE,5,,,1,1,OFF
LESIZE,8,,,2
LESIZE,9,,,2
MSHKEY,1
MSHAPE,0,2D
AMESH,2
LREFINE,3,,,1,1,OFF
LESIZE,2,,,4,
LESIZE,4,,,4,
MSHKEY,1
MSHAPE,0,2D
AMESH,1
ALLSEL,ALL
NUMCMP,ALL
ALLSEL,ALL
/REPLOT
FINISH
/SOLU
ANTYPE,TRANS
TRNOPT,FULL
TUNIF,293
nlgeom,off
nropt,full
pi=3.1415926
I=3.3
u=24
efficient=0.9
r=0.003
Qm=3*i*u*efficient/(pi*r**2)
kk=0.85
Qmax=Qm*kk
lz=0.0001
lx=0.045
ly=0.045
lsize=0.0025
v=0.0062
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
SFA,2,,CONV,68,293
SFA,3,,CONV,68,293
SFA,4,,CONV,68,293
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=1e05*Qmax*exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
OUTRES,BASIC,ALL
KBC,1
TSRES,ERASE
TIMINT,ON
LNSRCH,ON
TM1=0.005+(K-1)*TINC
TIME,TM1
SFA,1,,HFLUX,%FLUX2%
AUTOTS,-1
NSUBST,4
/GST,ON
/REP,FAST
/PSF,HFLUX,,1
/REPLOT
SOLVE
SFADELE,1,,HFLUX
TM2=K*TINC
TIME,TM2
SFA,1,,CONV,68,293
AUTOTS,-1
NSUBST,4
/GST,ON
SOLVE
SFADELE,1,,CONV
SAVE
*ENDDO
SAVE

回复此楼