分享移动高斯热原模型的Fortran子程序

小兰子mst

SUBROUTINE DFLUX(FLUX,SOL,KSTEP,KINC,TIME,NOEL,NPT,COORDS,
     1 JLTYP,TEMP,PRESS,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION FLUX(2), TIME(2), COORDS(3)
      CHARACTER*80 SNAME
     
      rad=0.000040         !Laser Radius
      q1=200                  !Laser Power
      pen=0.000070        !Penetration Depth
      height=0.000400     !Total Part Height
      x=COORDS(1)          !X Coordinate  
      y=COORDS(2)          !Y Coordinate
      z=COORDS(3)          !Z Coordinate
      R=0.2                     !reflectance of powder
      pi=3.14159
      v=0.8                     !scanningspeed
      heat=2*(1-R)*q1/(pen*(rad**2)*pi)            
      depth=exp(-1*abs(height-z)/pen)
      
      IF(TIME(2) .GT. 0 .AND. TIME(2) .LE. 800e-6)THEN
          shape1=exp(-2*((x-(v*TIME(2)+0.000150))**2+(y-0.000150)**2)/rad**2)   !Gaussian Distribution
          FLUX(1)=heat*shape1*depth      ! Defining the Body Flux
      ELSEIF(TIME(2).GT.800e-6)THEN
          FLUX(1)=0
      ENDIF
      
      RETURN
      END

啦啦啦啦啦

研究研究，谢谢

汉庭庭

66666666666666666

zhangxiang001

不错不错，学习一下

LYFX

楼主你好，请问你QQ是多少，想请教你啊

liuzhen.neu

学习学习

小兰子mst

LYFX 发表于 2017-10-13 14:29
楼主你好，请问你QQ是多少，想请教你啊

我的QQ号836803895