Laser Heating Simulation in ANSYS

[@Daredevil4798v]

How would you know the measure of the material removed from the contour plot

[@rezaamiri4520]

Hello, any one has the amended code ? please kindly help! mine is not working

[@kushagrakirtiman3292]

can you explain what the last plot mean? why the temperature was shown on left edge of the geometry?

[@pinalrana]

This is the axis-symmetric model, so the temperature is shown on the left edge.

[@bgborden6852]

@@pinalrana Can you explain what changes should I make if I want to put the laser beam in the middle of the workpiece

[@vincentfeng9802]

When i try to upload the code, it says line 1: parse error

[@thenoobers4628]

Hey Pinal, do you mind sharing the code file again? The link in the description seems to be broken

[@pinalrana]

Updated

[@jessefonarcus]

Im getting this error when i interpret 'function returning float returns double'

[@pinalrana]

Replace float with double

[@kushagrakirtiman3292]

Can you plz explain the temperature contour plot you are getting is top view or front view?

[@pinalrana]

Front view

[@bgborden6852]

@@pinalrana Thank you, so it showing the depth values of laser ablation?

[@pinalrana]

@@bgborden6852 Yes

[@pinalrana]

@@bgborden6852 Yes

[@bgborden6852]

Do you know why I am getting max temp as 5000K while you are only getting around 444K ?

[@pinalrana]

You may have used different laser intensity or material properties.

[@kushagrakirtiman3292]

I am getting temperature of 5000K with float I = 7.95455*pow(10,11), which is used in your code. Can you tell how are you getting max temp of just 400K? I have used the same material as you have used

[@pinalrana]

Actually, you are getting a temperature more than 5000 K. By default In Ansys the maximum temperature limit is 5000 K. However you can change the default value to some higher number. In my case, I have used a different source file. You are using the updated version of the source file, which shows more accurate results.

[@kushagrakirtiman3292]

@@pinalrana Do you know how to change the limit and is not drilling as temperature is quite is high more than melting point of the material?

[@pinalrana]

@@kushagrakirtiman3292 Workbench > Setup > Solution > Solution control > Limits > Maximum static temperature

[@kushagrakirtiman3292]

@@pinalrana thank you so much, Cant we say this simulation as drilling since max temp is more than the melting point

[@pinalrana]

@@kushagrakirtiman3292 Its depends on your assumptions and mechanism of material removal.

[@bgborden6852]

Is this not multiple pulses because your time step is 100 and the pulse on time is in ms? and if not can you suggest the changes that we can make in order to simulate two or three pulses ?

[@pinalrana]

You need to include time term in the code and make appropriate changes based on pulse on and pulse off time.

[@alinaderi333]

Hi Pinal, Thanks for your interesting tutorial. I have an error while adding the code. (Error, line 32: function returning float returns double). Can you please check what is the problem?

[@alinaderi333]

The name of the file is Laser_heating_1.c

[@pinalrana]

You have to write float instead of double.

[@alinaderi333]

@@pinalrana Dear Pinal, the error is removed. Thanks.

[@alinaderi333]

Excuse me, Did you use different parameters for the laser system in this code?

[@pinalrana]

@@alinaderi333 possible, I an not sure

[@anandasinha9098]

sir...not able to access the code file

[@pinalrana]

Sorry for the later response. I have updated the link.

[@anandasinha9098]

@@pinalrana thank u so much sir.... thank u :-)

[@kushagrakirtiman3292]

Q = (I/6.28)*(1-R)*a*exp(a*(y1))*exp(-1*(y*y)/(r0*r0)); Can plz tell from where exp(-1*(y*y)/(r0*r0)) coz in the formula it is exp(-2(r/w)^2) and why are u dividing the whole term by 6.28 moreover the need to add absorption and reflectivity

[@pinalrana]

Various types of profiles can be used for simulation such as 1/e, 1/e2, FWHM, etc. I have used a 1/e profile.

[@pinalrana]

a is absorption coefficient and R is reflectivity

[@kushagrakirtiman3292]

@@pinalrana Thank you so much, can you please tell why 6.28 is used to divide this?

[@pinalrana]

@@kushagrakirtiman3292 Its obtained from the integration of Beer-Lambert law.

[@jessefonarcus]

Sir you uploaded docx file.. it is asking for .c file. What should i do?

[@pinalrana]

Conver it into .c file.

[@jessefonarcus]

@@pinalrana i cant convert it...can you please help me?

[@kushagrakirtiman3292]

how to simulate multiple pulses ?

[@pinalrana]

It's computationally intensive. Still, it is possible to do in Ansys simply by editing code.

[@kushagrakirtiman3292]

@@pinalrana Thanks a lot, Can you tell me how to beam at the canter of the workpiece as in the video u have put it on the corner

[@kushagrakirtiman3292]

@@pinalrana Can you guide what changes should I make

[@rajivsinha4878]

Excellent simulation brother, I am new to C language so can you tell what the following line in your code does? C_CENTROID(x,c,t) y=x[0] y1=x[1] y2=x[2]

[@pinalrana]

Its defined co-ordinates in X, Y, and Z

[@rajivsinha4878]

@@pinalrana it would be great help if u explain that in a bit more detail. Thanks in advance

[@pinalrana]

@@rajivsinha4878 y=x[0] represents co-ordinates in X direction, y1=x[1] in Y, and y2=x[2] in Z direction.

[@rajivsinha4878]

@@pinalrana thanks a lot, So what does this line do? Coz I cannot understand like how the gaussian distribution was applied on the top face?

[@divyayadav5834]

Thanks for this video....can you please share any contact details....i needed a small help...for my project

[@pinalrana]

You can connect me on Linkedin: www.linkedin.com/in/pinal-rana-207413104/