24小时热门版块排行榜    

CyRhmU.jpeg
查看: 708  |  回复: 1

材料种子

新虫 (初入文坛)

[交流] 激光熔覆温度场模拟

各位,最近做ansys激光熔覆模拟时,发现做出来的温度忽高忽低,比如说第一载荷步后节点温度1600,第二个就降到600多,第三个又升到1600多,第四个又降回600多,一直反复,也不知道什么情况?有哪位遇到过吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

材料种子

新虫 (初入文坛)

!***********建立热源数组*****************
QMAX=P/(3.1415*R**2)
MAX_TIME=len/LSIZE+1
MAX_X=wid/LSIZE+1
MAX_Y=len/ LSIZE+1
!**********施加对流载荷*************
ESEL,S,LIVE
ASEL,S,AREA,,9,12,1
sfa,all,1,conv,150,20
ESEL,S,LIVE
ASEL,S,AREA,,7,7,1
sfa,all,1,conv,150,20

*do,h,1,road,1
m=mod(h,2)
*if,m,eq,1,then
!h=1
!热源数组
*DIM,FLUX1,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1 !循环时间
  *DO,I,1,MAX_X,1
  FLUX1(I,0,K)=(I-1)*LSIZE
  *ENDDO   
  *DO,J,1,MAX_Y,1
  FLUX1(0,J,K)=(J-1)*LSIZE
  *ENDDO   
   *DO,I,1,MAX_X,1
      *DO,J,1,MAX_Y,1
       XCENTER=(h-1)*weld+weld/2
       YCENTER=V*(K-1)*TINC
       DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
         *IF,DISTANCE,LE,R,THEN
         FLUX1(I,J,K)=2*A*QMAX*EXP(-2*DISTANCE**2/R**2)
         *ELSE
         FLUX1(I,J,K)=0
         *ENDIF
     *ENDDO
   *ENDDO
   FLUX1(0,0,K)=(K-1)*TINC
!激活单元
                esel,s,live
        nsel,s,loc,x,XCENTER-weld/2 ,XCENTER+weld/2
        nsel,r,loc,z,0,hth
        nsel,r,loc,y,YCENTER,YCENTER+LSIZE
        esln,s,0
        ealive,all
        allsel,all

!施加热流载荷
    ESEL,S,LIVE
    NSEL,S,LOC,Z,hth,hth!选择上边界
    SF,all,HFLUX,%FLUX1% !热流加载
    ALLSEL,ALL
    time,(h-1)*(TI+tinc)+(k)*TINC
    NSUBST,2
    SOLVE
  NSEL,S,LOC,Z,hth,hth!选择上边界
  SFDELE,ALL,HFLUX !删除上一步加载的热源
  /REPLOT
  *ENDDO

*else
*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=(h-1)*weld+weld/2
       YCENTER=len-V*(K-1)*TINC
       DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
         *IF,DISTANCE,LE,R,THEN
         FLUX2(I,J,K)=2*A*QMAX*EXP(-2*DISTANCE**2/R**2)
         *ELSE
         FLUX2(I,J,K)=0
         *ENDIF
     *ENDDO
   *ENDDO
   FLUX2(0,0,K)=(K-1)*TINC
!激活单元
              esel,s,live
        nsel,s,loc,x,XCENTER-weld/2 ,XCENTER+weld/2
        nsel,r,loc,z,0,hth
        nsel,r,loc,y,YCENTER,YCENTER+LSIZE
        esln,s,0
        ealive,all
        allsel,all

!施加热流载荷
    ESEL,S,LIVE
    NSEL,S,LOC,Z,hth,hth!选择上边界
    SF,all,HFLUX,%FLUX2% !热流加载
    ALLSEL,ALL
    time,(K)*TINC+(h-1)*(TI+tinc)
   NSUBST,2  
    SOLVE
  NSEL,S,LOC,Z,hth,hth!选择上边界
  SFDELE,ALL,HFLUX !删除上一步加载的热源
  /REPLOT
  *ENDDO
*endif
*enddo
FINISH
这是求解的部分程序,希望大佬帮我看看哪错了?
2楼2021-06-21 20:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 材料种子 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见