| 查看: 1053 | 回复: 2 | |||
[求助]
我写了一个焊接移动热源命令，总是出错，求大神帮看看 已有1人参与
|
|
命令流如下 /CLEAR,START /filname,ex-1 /COM,ANSYS RELEASE 12.0 ELLEN_LI 04/03/2014 /PREP7 ET,1,SOLID70 MPTEMP,1,100,273,298,373,438,500,600 MPDATA,KXX,1,1,0.02,0.3,0.8,0.9,1,1.2,1.8 MPDATA,C,1,1,400,1000,1186,1832.5,2093,2200,2100 MPDATA,DENS,1,1,420,420,450,600,900,900,950 WIDTH_BASE=0.1 HEIGHT_BASE=0.02 LENGTH=0.1 LAYER=1 WIDTH_SINTERING=0.05 HEIGHT_SINTERING=0.005*LAYER DT=0.00002 TEMP=298 INIT_TEMP=363 COFF_CONV=30 VELOCITY=2.0 STEPDIS_X=0.001 STEPDIS_Y=0.002 LASERPOWER=5000 RADIUS=0.005 AREA=3.14159265*(RADIUS**2) FACTOR=0.7 Q=LASERPOWER*FACTOR QM=3/AREA*Q STEPTIME_X=STEPDIS_X/VELOCITY STEPNUM_LENGTH=LENGTH/STEPDIS_X-10 STEPNUM_WIDTH=2*WIDTH_SINTERING/STEPDIS_Y-5 TOTALNUM_STEP=STEPNUM_LENGTH*STEPNUM_WIDTH TIME_LENGTH=STEPNUM_LENGTH*STEPTIME_X TIME_LAYER=TOTALNUM_STEP*STEPTIME_X BLOCK,0,LENGTH,-WIDTH_BASE,-WIDTH_SINTERING,0,-HEIGHT_BASE BLOCK,0,LENGTH,-WIDTH_SINTERING,WIDTH_SINTERING,0,-HEIGHT_BASE BLOCK,0,LENGTH,WIDTH_SINTERING,WIDTH_BASE,0,-HEIGHT_BASE BLOCK,0,LENGTH,-WIDTH_SINTERING,WIDTH_SINTERING,0,HEIGHT_SINTERING VGLUE,ALL NUMCMP,ALL LPLOT,ALL /PNUM,LINE,1 /REPLOT SAVE,ALL LSEL,S,LOC,X,DT,LENGTH-DT LESIZE,ALL,,,LENGTH/STEPDIS_X,,,,,0, ALLSEL,ALL LSEL,S,LOC,Y,-WIDTH_SINTERING+DT,WIDTH_SINTERING*2/STEPDIS_X,,,,,0, LESIZE,ALL,,,WIDTH_SINTERING*2/STEPDIS_X,,,,,0, ALLSEL,ALL LESIZE,17,,,20,1,,,,0 LESIZE,18,,,20,1,,,,0 LESIZE,14,,,20,1,,,,0 LESIZE,23,,,20,1,,,,0 LESIZE,13,,,20,1,,,,0 LESIZE,7,,,20,1,,,,0 LESIZE,21,,,20,1,,,,0 LESIZE,2,,,20,1,,,,0 LESIZE,5,,,20,1,,,,0 LESIZE,4,,,20,1,,,,0 LESIZE,25,,,20,1,,,,0 LESIZE,28,,,20,1,,,,0 LESIZE,24,,,20,1,,,,0 LESIZE,22,,,20,1,,,,0 LESIZE,27,,,20,1,,,,0 LESIZE,26,,,20,1,,,,0 ALLSEL,ALL LESIZE,19,,,4,1,,,,0 LESIZE,20,,,4,1,,,,0 LESIZE,15,,,4,1,,,,0 LESIZE,16,,,4,1,,,,0 LESIZE,11,,,4,1,,,,0 LESIZE,12,,,4,1,,,,0 LESIZE,10,,,4,1,,,,0 LESIZE,9,,,4,1,,,,0 ALLSEL,ALL LESIZE,33,,,1,1,,,,0 LESIZE,34,,,1,1,,,,0 LESIZE,35,,,1,1,,,,0 LESIZE,36,,,1,1,,,,0 ALLSEL,ALL LESIZE,6,,,10,,,,,0 LESIZE,8,,,10,,,,,0 LESIZE,31,,,10,,,,,0 LESIZE,32,,,10,,,,,0 LESIZE,29,,,10,,,,,0 LESIZE,30,,,10,,,,,0 LESIZE,1,,,10,,,,,0 LESIZE,3,,,10,,,,,0 VSEL,ALL TYPE,1 MAT,1 MSHAPE,0,3D MSHKEY,1 VMESH,ALL NSEL,S,LOC,Z,-HEIGHT_BASE,0 IC,ALL,TEMP.INIT_TEMP ALLSEL,ALL NSEL,S,LOC,Z,-HEIGHT_BASE,0 NSEL,R,LOC,Y,WIDTH_BASE SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL NSEL,S,LOC,Z,-HEIGHT_BASE,0 NSEL,R,LOC,Y,-WIDTH_BASE SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL NSEL,S,LOC,Z,-HEIGHT_BASE,0 NSEL,R,LOC,X,0 SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL NSEL,S,LOC,Z,-HEIGHT_BASE,0 NSEL,R,LOC,X,LENGTH SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL NSEL,S,LOC,Z,0 NSEL,R,LOC,Y,WIDTH_SINTERING,WIDTH_BASE SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL NSEL,S,LOC,Z,0 NSEL,R,LOC,Y,-WIDTH_BASE,-WIDTH_SINTERING SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL FINISH /SOLU ANTYPE,4 TRNOPT,FULL LUMPM,0 DELTIM,0.001,0.0002,0.002 CNVTOL,HEAT,,0.01,1,0.000001 OUTRES,NSOL,LAST NSEL,S,LOC,Z,O,HEIGHT_SINTERING ESLN,S,1,ALL EKILL,ALL ALLSEL,ALL ESEL,S,LIVE K=1!*DO,K,1,LAYER,1 NSEL,S,LOC,Z,(K-1)*HEIGHT_SINTERING/LAYER,K*HEIGHT_SINTERING/LAYER ESLN,S,1 EALIVE,ALL ALLSEL,ALL ESEL,S,LIVE EPLOT NSEL,S,LOC,Z,HEIGHT_SINTERING/LAYER*(K-1),HEIGHT_SINTERING/LAYER*K IC,ALL,TEMP,INIT_TEMP ALLSEL,ALL *DO,J,1,STEPNUM_WIDTH,2 *DO,I,1,STEPNUM_LENGTH,1 *SET,XCENTER,(I+1)*STEPDIS_X *SET,YCENTER,J*STEPDIS_Y-WIDTH_SINTERING *SET,ZCENTER,HEIGHT_SINTERING/LAYER*K LOCAL,11,1,XCENTER,YCENTER,ZCENTER NSEL,S,LOC,X,,RADIUS NSEL,R,LOC,Y,,360 NSEL,R,LOC,Z,HEIGHT_SINTERING/LAYER*(K-1),HEIGHT_SINTERING/LAYER*K CM,NODEF,NODE LOCAL,12,0,XCENTER,YCENTER,ZCENTER X=CENTRX(I) Y=CENTRY(I) HF=QM*exp(-3*((X-XCENTER)**2+(Y-YCENTER)**2)/RADIUS**2) SF,NODEF,HFLUX,HF ALLSEL,ALL CSYS,0 NSEL,S,LOC,X,0-DT,XCENTER-RADIUS,0.001 NSEL,A,LOC,X,0-DT,XCENTER+RADIUS,LENGTH, NSEL,A,LOC,Y,0-DT,-WIDTH_SINTERING,YCENTER-RADIUS, NSEL,A,LOC,Y,0-DT,YCENTER+RADIUS,WIDTH_SINTERING, NSEL,R,LOC,Z,HEIGHT_SINTERING/LAYER*K, SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL TIME,I*STEPTIME_X+(J-1)*TIME_LENGTH+(K-1)*TIME_LAYER OUTRES,,ALL SOLVE SFDELE,ALL,HFLUX CMDELE,NODEF SFDELE,ALL,CONV ALLSEL,ALL *IF,I,EQ,STEPNUM_LENGTH,THEN *DO,L,STEPNUM_LENGTH,1,-1 *SET,XCENTER,(L+1)*STEPDIS_X *SET,YCENTER,(J+1)*STEPDIS_Y-WIDTH_SINTERING *SET,ZCENTER,HEIGHT_SINTERING/LAYER*K LOCAL,11,XCENTER,YCENTER,ZCENTER NSEL,S,LOC,X,RADIUS NSEL,R,LOC,Y,360 NSEL,R,LOC,Z,HEIGHT_SINTERING/LAYER*(K-1),HEIGHT_SINTERING/LAYER*K CM,NODEF,NODE LOCAL,12,O,XCENTER,YCENTER,ZCENTER X=CENTRX(I) Y=CENTRY(I) HF=QM*exp(-3*((X-XCENTER)**2+(Y-YCENTER)**2)/RADIUS**2) SF,NODEF,HFLUX,HF ALLSEL,ALL NSEL,S,LOC,Z,HEIGHT_SINTERING/LAYER*K NSEL,U,LOC,Y,YCENTER-RADIUS,YCENTER+RADIUS NSEL,U,LOC,X,XCENTER-RADIUS,XCENTER+RADIUS SF,ALL,CONV,COFF_CONV,TEMP ALLSEL,ALL CSYS,0 TIME,(STEPNUM_LENGTH-L+1)*STEPTIME_X+J*TIME_LENGTH+(K-1)*TIME_LAYER OUTRES,,ALL SOLVE SFDELE,ALL,HFLUX CMDELE,NODEF SFDELE,ALL,CONV *ENDDO *ENDIF *ENDDO *ENDDO SAVE,ALL |
回复此楼» 收录本帖的淘帖专辑推荐
 焊接有限元分析 |
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
2楼2014-04-14 23:30:25
11050209
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 938.2
- 散金: 100
- 帖子: 321
- 在线: 68.3小时
- 虫号: 3896498
- 注册: 2015-05-28
- 专业: 金属材料的凝固与结晶学
3楼2016-05-14 13:53:05