| 查看: 710 | 回复: 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 |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有14人回复
-
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
ANSYS四维TABLE建立移动三维热源，程序总是出错，求原因解释！
已经有16人回复
- 求助高斯热源中对流边界的设置
已经有4人回复
- 请教一下，如何利用Ansys在一个长方形中创建50个椭圆。
已经有12人回复
- ansys焊接热应力分析温度场是对的，残余应力分布是错的
已经有11人回复
- 最近学ansys遇到的一些难题！！！
已经有5人回复
- 【讨论】ansys如何加载高斯分布的热流密度？
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100