- 积分
- 0
- 注册时间
- 2017-7-9
- 仿真币
-
- 最后登录
- 1970-1-1
|
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 |
|