24小时热门版块排行榜

返回列表

虫之恋

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3.5
帖子: 4
在线: 1小时
虫号: 2758678
注册: 2013-10-28
专业: 金属结构材料

[求助] 我写了一个焊接移动热源命令，总是出错，求大神帮看看已有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

回复此楼