| 查看: 1012 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 焊接有限元分析 |
» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有9人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有20人回复
-
2026国自然放榜时间
已经有9人回复
-
这样的filecode谁见过
已经有22人回复
-
据悉今年马上要出结果了
已经有6人回复
-
面上项目filecode邪修
已经有7人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
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
20