24小时热门版块排行榜    

查看: 657  |  回复: 0

帅帅小城

新虫 (初入文坛)

[交流] 请大神帮我检查一下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
回复此楼
成功需要不断的努力
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 帅帅小城 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见