[Unkinunki]

Hi,
I am currently working with a transient thermal analysis for a welding process. For this model, I use Goldak’s double ellipsoid heat source, which is defined using APDL commands. Geometry, mesh and timesteps are determined in Ansys Mechanical.
In order to achieve a finger type penetration, I have to combine two heat sources. The problem is that I have not been able to combine them simultaneously. As a first trial, I have been trying to simulate two heat sources at different coordinate points. However, only one of them is applied. Does anyone know how to solve this problem?

I have attached my APDL codes and some figures to give a better understanding.

Thanks in advance.

[kipa]

cut body ，one body add one source

[Juan Daniel]

Hi,
¿Do you still need help about this topic?

[shekhars.ip.17@nitj.ac.in]

Hi All,
I have been trying to model the heat source and succeded for a single bead only. For 2nd bead, the modification in the code has not worked yet.
Can you give me some lead?

[Juan Daniel]

Hi, can you share your code, I like to check some errors or to give you some tips.

[shekhars.ip.17@nitj.ac.in]

Here is my code for two moving heat source one after one:-

C*** CONSTANTS DEFINITION

VOLT = 175 ! VOLTAGE (VOLTS)
CURR = 12.5 ! CURRENT (AMPERES)
EFF = 0.90 ! EFFICIENCY

VEL = 4E-3 ! WELDING SPEED (m/s)

Q = EFF*VOLT*CURR ! HEAT INPUT (W)

A = 5E-3 ! ELLIPSOID DIMENSIONS (m)
B = 2.5E-3
C1 = 8E-3
C2 = 10E-3

FF = 0.6 ! HEAT FRACTIONS
FR = 1.4

TAU = (C1/VEL) ! LAG FACTOR (s)
TTIME = 65 ! TOTAL TIME (s)
!!TIME_WELD = 56 ! WELDING TIME (s)
DT = 1 ! TIME INCREMENT (s)

PI = 3.14
EE = 2.718

C*** HEAT SOURCE

CMSEL,S,WELD_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

!!NPT = TIME_WELD/DT ! NUMBER OF LOADSTEPS 1

*DO,ii,1,28,1

WTIME = (ii/2)
TIME,WTIME

*DO,jj,EMIN,EMAX,1

CSYS,20
*GET,X,ELEM,jj,CENT,X
*GET,Y,ELEM,jj,CENT,Y
*GET,Z,ELEM,jj,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART1 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART2 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF1 = PART1*PART2

BFE,jj,HGEN,,QF1

*ENDDO

SOLVE

*ENDDO

CMSEL,S,Substrate_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

*DO,uu,29,60,1

WTIME = (uu/2)
TIME,WTIME

*DO,kk,EMIN,EMAX,1

CSYS,30
*GET,X,ELEM,kk,CENT,X
*GET,Y,ELEM,kk,CENT,Y
*GET,Z,ELEM,kk,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART3 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART4 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF2 = PART3*PART4

BFE,kk,HGEN,,QF2

*ENDDO

SOLVE

*ENDDO

CMSEL,ALL
BFEDELE,ALL,ALL
ALLSEL

TIME,TTIME
DELTIM,DT*51.jpg

2.jpg

[Juan Daniel]

My apologies for the delay, Y had many troubles using command BF because APDL identified the nodes but the program warned that they were not selected. I tried many ways to fix this warning message and I did it with de CMSEL Command.

	SF,NSURF,CONV,H,TAR
		CSYS,1000
		CMSEL,S,NSPOT
		BF,ALL,HGEN,%BEAM%
		!SF,NSPOT,HFLUX,%BEAM%
		ALLSEL,ALL				!SELECT ALL ENTITIES

If your trouble persist, send me warning message from ANSYS APDL.

Originally Posted by shekhars.ip.17@nitj.ac.in

Here is my code for two moving heat source one after one:-

C*** CONSTANTS DEFINITION

VOLT = 175 ! VOLTAGE (VOLTS)
CURR = 12.5 ! CURRENT (AMPERES)
EFF = 0.90 ! EFFICIENCY

VEL = 4E-3 ! WELDING SPEED (m/s)

Q = EFF*VOLT*CURR ! HEAT INPUT (W)

A = 5E-3 ! ELLIPSOID DIMENSIONS (m)
B = 2.5E-3
C1 = 8E-3
C2 = 10E-3

FF = 0.6 ! HEAT FRACTIONS
FR = 1.4

TAU = (C1/VEL) ! LAG FACTOR (s)
TTIME = 65 ! TOTAL TIME (s)
!!TIME_WELD = 56 ! WELDING TIME (s)
DT = 1 ! TIME INCREMENT (s)

PI = 3.14
EE = 2.718

C*** HEAT SOURCE

CMSEL,S,WELD_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

!!NPT = TIME_WELD/DT ! NUMBER OF LOADSTEPS 1

*DO,ii,1,28,1

WTIME = (ii/2)
TIME,WTIME

*DO,jj,EMIN,EMAX,1

CSYS,20
*GET,X,ELEM,jj,CENT,X
*GET,Y,ELEM,jj,CENT,Y
*GET,Z,ELEM,jj,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART1 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART2 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF1 = PART1*PART2

BFE,jj,HGEN,,QF1

*ENDDO

SOLVE

*ENDDO

CMSEL,S,Substrate_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

*DO,uu,29,60,1

WTIME = (uu/2)
TIME,WTIME

*DO,kk,EMIN,EMAX,1

CSYS,30
*GET,X,ELEM,kk,CENT,X
*GET,Y,ELEM,kk,CENT,Y
*GET,Z,ELEM,kk,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART3 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART4 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF2 = PART3*PART4

BFE,kk,HGEN,,QF2

*ENDDO

SOLVE

*ENDDO

CMSEL,ALL
BFEDELE,ALL,ALL
ALLSEL

TIME,TTIME
DELTIM,DT*5Attachment 77950

Attachment 77951

[dharmik_85i]

Originally Posted by shekhars.ip.17@nitj.ac.in

Here is my code for two moving heat source one after one:-

C*** CONSTANTS DEFINITION

VOLT = 175 ! VOLTAGE (VOLTS)
CURR = 12.5 ! CURRENT (AMPERES)
EFF = 0.90 ! EFFICIENCY

VEL = 4E-3 ! WELDING SPEED (m/s)

Q = EFF*VOLT*CURR ! HEAT INPUT (W)

A = 5E-3 ! ELLIPSOID DIMENSIONS (m)
B = 2.5E-3
C1 = 8E-3
C2 = 10E-3

FF = 0.6 ! HEAT FRACTIONS
FR = 1.4

TAU = (C1/VEL) ! LAG FACTOR (s)
TTIME = 65 ! TOTAL TIME (s)
!!TIME_WELD = 56 ! WELDING TIME (s)
DT = 1 ! TIME INCREMENT (s)

PI = 3.14
EE = 2.718

C*** HEAT SOURCE

CMSEL,S,WELD_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

!!NPT = TIME_WELD/DT ! NUMBER OF LOADSTEPS 1

*DO,ii,1,28,1

WTIME = (ii/2)
TIME,WTIME

*DO,jj,EMIN,EMAX,1

CSYS,20
*GET,X,ELEM,jj,CENT,X
*GET,Y,ELEM,jj,CENT,Y
*GET,Z,ELEM,jj,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART1 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART2 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF1 = PART1*PART2

BFE,jj,HGEN,,QF1

*ENDDO

SOLVE

*ENDDO

CMSEL,S,Substrate_1 ! SELECT THE WELD COMPONENT

*GET,EMAX,ELEM,,NUM,MAX ! GET THE MAXIMUM ELEMENT NUMBER
*GET,EMIN,ELEM,,NUM,MIN ! GET THE MINIMUM ELEMENT NUMBER

ALLSEL

*DO,uu,29,60,1

WTIME = (uu/2)
TIME,WTIME

*DO,kk,EMIN,EMAX,1

CSYS,30
*GET,X,ELEM,kk,CENT,X
*GET,Y,ELEM,kk,CENT,Y
*GET,Z,ELEM,kk,CENT,Z

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

*IF,(TAU-WTIME),GT,0,THEN

C = C1
F = FF

*ELSE

C = C2
F = FR

*ENDIF

PART3 = 6*(3**0.5)*F*Q/(A*B*C*PI*(PI**0.5))
PART4 = (EE**(-3*(X/A)**2))*(EE**(-3*(Z/B)**2))*(EE**(-3*(CSI/C)**2))
QF2 = PART3*PART4

BFE,kk,HGEN,,QF2

*ENDDO

SOLVE

*ENDDO

CMSEL,ALL
BFEDELE,ALL,ALL
ALLSEL

TIME,TTIME
DELTIM,DT*5

Attachment 77950

Attachment 77951

I think for that you have to named selection that bead component and change location value for in do loop for I,j,k may be it works

[skv6]

hi everyone

I need a heat source for pipe joint welding please help me any for APDL GOLDAK heat source code

thank you

[skv6]

I need a heat source for pipe joint welding please help me any for APDL GOLDAK heat source code

[shekhars.ip.17@nitj.ac.in]

Share your attempt with us. then we can help