×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Are you an
Engineering professional?
Join Eng-Tips Forums!
• Talk With Other Members
• Be Notified Of Responses
• Keyword Search
Favorite Forums
• Automated Signatures
• Best Of All, It's Free!

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

#### Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

# Ansys thermal analysis |Simulation | Moving heat soure

## Ansys thermal analysis |Simulation | Moving heat soure

(OP)
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.
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

#### Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

#### Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

#### Resources

Low-Volume Rapid Injection Molding With 3D Printed Molds
Learn methods and guidelines for using stereolithography (SLA) 3D printed molds in the injection molding process to lower costs and lead time. Discover how this hybrid manufacturing process enables on-demand mold fabrication to quickly produce small batches of thermoplastic parts. Download Now
Examine how the principles of DfAM upend many of the long-standing rules around manufacturability - allowing engineers and designers to place a partâ€™s function at the center of their design considerations. Download Now
Taking Control of Engineering Documents
This ebook covers tips for creating and managing workflows, security best practices and protection of intellectual property, Cloud vs. on-premise software solutions, CAD file management, compliance, and more. Download Now

Close Box

# Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

• Talk To Other Members
• Notification Of Responses To Questions
• Favorite Forums One Click Access
• Keyword Search Of All Posts, And More...

Register now while it's still free!