Simulation | Heat source model | Welding

[Sra1ksk]

Am working on simulation of the moving heat source of double ellisopidal heat source model (In ansys ACT extension Gaussian moving heat source is only available)..To find out the temperature distribution along the weldzone. Giving convection and moving heat source (UDF) as input. In the code the mathemtical equation of double ellipsoidal is given. But am not getting the temperature distribution reference to the double ellipsoidal graph.
Error with temperature profile in the ouput. Can any one help me. please.
Code is attached. Can you please check the code..

CMSEL,ALL

*GET,EMAX,ELEM,,NUM,MAX

*GET,EMIN,ELEM,,NUM,MIN

ALLSEL

AUTOTS,ON

TIME_WELD=100

DT=1 ! Time step

A=0.003

B=0.004

C1=0.004

C2=0.016 !A, B AND C ARE JUST THE HEAT SOURCE DIMENSIONS

TAU=0 !TIME DELAY

FF=0.6 !FRACTION OF HEAT ON THE FRONT PART OF THE HEAT SOURCE

FR=1.4 !FRACTION OF HEAT ON THE REAR PART OF THE HEAT SOURCE

Q=3840 !WELDING ENERGY J/S (W)

VEL=1.73e-3 !WELDING SPEED

!!EE=2.71 !I HAD TO DEFINE THE EULER NUMBER AS 2.71 INSTEAD OF USING exp() DUE TO UNSTABILITIES ON MY RESULTS

NPT=TIME_WELD/DT

NROPT, FULL

*DO,i,1,NPT,1 !TIME INCREMENT

WTIME=(i/2)

TIME,WTIME

HCENTER=VEL*WTIME !HEAT SOURCE CENTER

*DO,jj,EMIN,EMAX,1 !SELECTS THE CENTROID OF EVERY ELEMENT

X=CENTRX(jj)

Y=CENTRY(jj)

Z=CENTRZ(jj)

CSI=Z+(VEL*(TAU-WTIME)) !MOVING COORDINATE (IT MOVES THROUGH Z AS TIME PASSES

*IF,Z,GT,HCENTER,THEN !DEFINES WHETHER IT'S THE FRONT PART OF THE HEAT SOURCE OR THE REAR PART

C=C1

F=FF

*ELSE

C=C2

F=FR

*ENDIF

PART1=(6*(3**0.5)*Q*F)/(A*B*C*3.14*(3.14**0.5))

PART2=(Exp(-3*(X/A)**2))*(Exp(-3*(Y/B)**2))*(Exp(-3*((CSI/C)**2)))

QF=PART1*PART2 !HEAT SOURCE EQUATION

BFE,jj,HGEN,,QF

*ENDDO

SOLVE

*ENDDO

CMSEL,ALL

BFEDELE,ALL,ALL

ALLSEL