I would like to simulate a simple additive manufacturing process via welding.
I am using the code below and I have a problem with ekill and alive commands together with the moving heat source.
I can have welding done in WELD1, but it is not possible to continue welding in WELD2.
WELD contains both WELD1 and WELD2.
Can anybody give some advice?
CMSEL,S,WELD
EKILL,ALL
!FIRST WELD FIRST SECOND
CMSEL,S,WELD1
EALIVE,WELD1
CMSEL,ALL
*GET,EMAX,ELEM,,NUM,MAX
*GET,EMIN,ELEM,,NUM,MIN
ALLSEL
TIME_WELD=6
DT=0.1
A=0.005
B=0.005
C1=0.005
C2=0.015
TAU=0
FF=0.5
FR=1.5
Q=1500
VEL=5e-3
NPT=TIME_WELD/DT
NROPT, FULL
*DO,i,1,NPT,1
WTIME=(i/10)
TIME,WTIME+0*TIME_WELD
HCENTER=VEL*WTIME
*DO,jj,EMIN,EMAX,1
*GET, X, elem, jj, cent, X
*GET, Y, elem, jj, cent, Y
*GET, Z, elem, jj, cent, Z
CSI=Y+VEL*(TAU-WTIME)
*IF,Y,GT,HCENTER,THEN
C=C1
F=FF
*ELSE
C=C2
F=FR
*ENDIF
PART1=(6*(3**0.5)*F*Q)/(A*B*C*3.14*(3.14**0.5))
PART2=(exp(-3*(Z/A)**2))*(exp(-3*(X/B)**2))*(exp(-3*(CSI/C)**2))
QF=PART1*PART2
BFE,jj,HGEN,,QF
*ENDDO
SOLVE
*ENDDO
CMSEL,ALL
BFEDELE,ALL,ALL
ALLSEL
!SECOND WELD SECOND SECOND
CMSEL,S,WELD2
EALIVE,WELD2
*DO,i,1,NPT,1
WTIME=(i/10)
TIME,WTIME+1*TIME_WELD
HCENTER=VEL*WTIME
*DO,jj,EMIN,EMAX,1
*GET, X, elem, jj, cent, X
*GET, Y, elem, jj, cent, Y
*GET, Z, elem, jj, cent, Z
CSI=Y+VEL*(TAU-WTIME)
*IF,Y,GT,HCENTER,THEN
C=C1
F=FF
*ELSE
C=C2
F=FR
*ENDIF
PART1=(6*(3**0.5)*F*Q)/(A*B*C*3.14*(3.14**0.5))
PART2=(exp(-3*(Z/A)**2))*(exp(-3*(X/B)**2))*(exp(-3*(CSI/C)**2))
QF=PART1*PART2
BFE,jj,HGEN,,QF
*ENDDO
SOLVE
*ENDDO
CMSEL,ALL
BFEDELE,ALL,ALL
ALLSEL