bernardirafael
Mechanical
Hello all! I am currently trying to make a moving heat source on ANSYS Workbench v14.5, and since I am a newbie I am using a script I've found on a Mastery's Degree thesis regarding welding. It uses the same heat source that I want to use, called double ellipsoid.
The problem is that it's not working, so I'd like to ask your help to solve this problem. The script is shown below:
I am the one who commented those lines based on what I've figured out so far. It sometimes plots random heat spots on the plate, and some other times it plots a weird heat source located nearly [-0.005,0,0] (meters). The temperature grows near that spot ([-0.005,0,0]) as some random points of heat. So far I couldn't makea moving heat source. Can someone help me on this? It would be really appreciated! Thanks in advance!
The problem is that it's not working, so I'd like to ask your help to solve this problem. The script is shown below:
Code:
/PREP7
ET,1,SOLID70
/SOL
CMSEL,ALL
*GET,EMAX,ELEM,,NUM,MAX
*GET,EMIN,ELEM,,NUM,MIN
ALLSEL
TIME_WELD=200
DT=1
A=0.002 !Dimention A of the ellipsoid
B=0.003 !Dimention B of the ellipsoid
C1=0.004 !Dimension C1 of the ellipsoid 1
C2=0.016 !Dimension C2 of the ellipsoid 2
TAU=-1 !Delay to position the heat source on the plate
FF=0.4 !Proportion Factor (front)
FR=1.6 !Proportion Factor (rear)
Q=2400 !Energy
VEL=0.0025 !Welding velocity
NPT=TIME_WELD/DT
*DO,ii,1,NPT,1
WTIME=(ii/2)
TIME,WTIME
HCENTER=VEL*WTIME
*DO,jj,EMIN,EMAX,1
X=CENTRZ(jj)
y=CENTRZ(jj)
Z=CENTRZ(jj)
CSI=Z+VEL*(TAU-WTIME) !Moving coordinate system equation, making the global equation a function of [x,y,z,t].
! The Heat source must move only through the Z coordinate.
*IF,Z,GT,HCENTER,THEN
!It should make a verification to determine the ellipsoid on the front of the heat source and the other ellipsoid
!on the rear of the heat source
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/B)**2))*(exp(-3*((CSI/C)**2)))
QF=PART1*PART2 !Global equation, a function of [x,y,z,t]
BFE,ii,HGEN,,QF
*ENDDO
SOLVE
*ENDDO
TIME,100
DELTIM,0.5,0.01,0.5,1I am the one who commented those lines based on what I've figured out so far. It sometimes plots random heat spots on the plate, and some other times it plots a weird heat source located nearly [-0.005,0,0] (meters). The temperature grows near that spot ([-0.005,0,0]) as some random points of heat. So far I couldn't makea moving heat source. Can someone help me on this? It would be really appreciated! Thanks in advance!
