24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1420  |  回复: 0
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 ZJ970617 的 5 个金币

ZJ970617

新虫 (初入文坛)

[交流] Abaqus高斯热源子程序

小弟最近在研究激光仿真,对旋转高斯体热源子程序方面有些小问题,希望能和大家多多交流。


   SUBROUTINE DFLUX(FLUX,SOL,JSTEP,JINC,TIME,NOEL,NPT,COORDS,JLTYP,
     1                 TEMP,PRESS,SNAME)
C
c      use GlobalVar
      INCLUDE 'ABA_PARAM.INC'
C
      parameter(one=1.d0)
      DIMENSION COORDS(3),FLUX(2),TIME(2)
      CHARACTER*80 SNAME
C     高斯旋转体热源参数,H为热源高度,R0为热源作用有效半径
      Q=2000
      effi=0.8
      P=Q*effi
      R0=0.004
      H=0.005
      v=0.01
      pi=3.1415
      
        x=COORDS(1)
        y=COORDS(2)
        z=COORDS(3)
      
C     初始热源中心坐标
        x0=0
        y0=0   
      z0=0
C     热源形状及热源公式,d当前时刻焊接斑点中心跟焊接初始位置的距离
      d=v*TIME(2)
      Cs=3/R0**2
      K=3*Cs*P/(pi*H*(1-1/(EXP(3)))*EXP(-3*Cs/LOG(H/Z)*((X-d)**2+Y**2)
      
      
C     JLTYP=1,体热源
      JLTYP=1
        FLUX(1)=3*Cs*P/(pi*H*(1-1/(EXP(3)))*EXP(-3*Cs/LOG(H/Z)*((X-d)**2+Y**2)
      RETURN
      end subroutine DFLUX
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ZJ970617 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见