Learning Forum

[SteBir]

Hi,

I use an APDL command snippet to achieve a double ellipsoid according to Goldak equation, but the Temperature rise at the end of the welding line is strange and artificial.

The APDL snippet is:

 

 

 

CMSEL,ALL

*GET,EMAX,ELEM,,NUM,MAX

*GET,EMIN,ELEM,,NUM,MIN

ALLSEL

 

TIME_WELD=80

DT=1

A=0.003

B=0.003

C1=0.004

C2=0.016

TAU=0

FF=0.6

FR=1.4

Q=1000

VEL=1.73e-3

 

NPT=TIME_WELD/DT

NROPT, FULL

 

*DO,i,1,NPT,1

WTIME=(i)

TIME,WTIME

 

HCENTER=VEL*WTIME

 

*DO,jj,EMIN,EMAX,1

 

X=CENTRX(jj)

Y=CENTRY(jj)

Z=CENTRZ(jj)

 

CSI=-Z+VEL*(TAU-WTIME)

 

 

*IF,Z,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*(X/A)**2))*(exp(-3*(Y/**2))*(exp(-3*(CSI/C)**2))

 

QF=PART1*PART2

 

BFE,jj,HGEN,,QF

 

*ENDDO

 

SOLVE

 

*ENDDO

CMSEL,ALL

BFEDELE,ALL,ALL

ALLSEL

This is when the Laser spot leaves the plate.

Has anyone a suggestion what went wrong?

Best regards

[Shekhar]

Have you performed any experiments and compared the thermal profile with the simulated result?

If not, then try it. You will find it the same as you got in the simulated result.

[SteBir]

Thanks for your answer.

You're right, I made some experimental test to validate the model and it's like you said.

Best regards

Ste

[ChristianDaniel]

Hello SteBir,nRight now I am trying to implement a Goldak Heat Source aswell. Just yesterday I stumbled upon a youtube video where a very similar code gets demonstrated (

&ab_channel=S.B.). Could you find any validation for the code? Is it suitable to simulate the welding process?nI would be really grateful for any tips you can share with me.nBest Regards,nDanieln