Moving Heat Source Problems
Moving Heat Source Problems
(OP)
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 --> APDL
/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,1
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!




