24小时热门版块排行榜    

查看: 899  |  回复: 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
回复此楼

» 收录本帖的淘帖专辑推荐

焊接有限元分析

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11050209

银虫 (正式写手)

引用回帖:
2楼: Originally posted by lvss2007 at 2014-04-14 23:30:25
提示什么错误呢,看了下,基本上都是用表数组加载热源载荷的,你这种形式的还是第一次见

我想问下数组怎么加载热源啊
3楼2016-05-14 13:53:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

lvss2007

新虫 (初入文坛)

【答案】应助回帖

提示什么错误呢,看了下,基本上都是用表数组加载热源载荷的,你这种形式的还是第一次见
2楼2014-04-14 23:30:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见