Numerical simulation method of pulsed laser paint removal and use thereof

Abstract

The present disclosure provides a numerical simulation method of pulsed laser paint removal and a use thereof. This method establishes a three-dimensional (3D) temperature field model by ANSYS software to perform a numerical simulation of nanosecond pulsed laser paint removal. A high-speed moving pulsed laser is loaded on a surface of the model in a form of heat flux, and a coordinate system is moved to realize loading on different paths. A special surface mesh screening method is used to realize loading on any surface, and it ensures that laser energy distribution on a material surface is in line with reality. In addition, an element birth/death technology is combined to remove an element that exceeds a threshold, so as to intuitively present the surface morphology after laser paint removal. The present disclosure can realize the prediction of the contour of a paint layer ablated by a pulsed laser.

Claims

1. A numerical simulation method of pulsed laser paint removal, comprising the following steps: step 1: establishing a two-layer three-dimensional (3D) solid model of a paint layer and a substrate, and conducting meshing, wherein in the 3D solid model, a selected element type is an 8-node hexahedral thermal element; step 2: setting an initial condition of a temperature field of the model to an ambient temperature, and setting a load option to transient analysis and a load mode to step load; step 3: establishing a coordinate system with a center of a single-pulse spot on an upper surface of the paint layer as an origin, taking an axial incident direction of a laser as a z-axis, taking the origin of the coordinate system as a pulsed laser loading position, and defining the upper surface of the paint layer as a z=0 plane; step 4: selecting nodes having a perpendicular distance not exceeding a radius of the spot from the z-axis in the 3D solid model to form a laser irradiation coverage zone, and selecting nodes with less than 8 surviving elements attached, among the nodes in the laser irradiation coverage zone, to form a laser irradiated surface; step 5: loading a pulsed laser onto a surface of the model in a form of heat flux, the elements attached to the nodes at the laser irradiated surface being laser irradiated elements; and calculating an energy load q received by a node on an upper surface of a laser irradiated element whose center has a distance of h from the z-axis: $\begin{array}{cc}q\left(h,t\right)=\frac{12.Math.\mathrm{aP}}{f\tau .Math.\pi d^2}.Math.\exp \left(-\frac{12.Math.h^2}{d^2}\right).Math.{\left(\frac{t}{\tau }\right)}^7\exp \left[7\left(1-\frac{t}{\tau }\right)\right]& \left(1\right)\end{array}$ wherein, a is an absorption rate of laser energy by a material, P is a laser output power, f is a laser repetition frequency, d is a spot diameter, τ is a laser pulse width, and t is a time, 0<t<τ; step 6: performing iterative solution in a form of time integration according to a heat conservation law of a heat transfer theory, to calculate an instantaneous temperature of each of the nodes; and step 7: defining an instantaneous temperature of the surviving elements reaching a threshold temperature or above as a condition for an element birth/death operation, and killing a surviving element whose instantaneous temperature exceeds the threshold temperature, to obtain morphology after the removal of the paint layer by a single pulsed laser.

2. The numerical simulation method of pulsed laser paint removal according to claim 1, wherein the threshold temperature is a vaporization temperature of a material of the paint layer.

3. The numerical simulation method of pulsed laser paint removal according to claim 1, wherein the heat conservation law of the heat transfer theory is that a temperature change of the material subjected to laser irradiation satisfies: $\begin{array}{cc}\frac{\partial T}{\partial t}=k\left(\frac{{\partial }^{2}T}{\partial x^2}+\frac{{\partial }^{2}T}{\partial y^2}+\frac{{\partial }^{2}T}{\partial z^2}\right)& \left(2\right)\end{array}$ wherein, Tis a temperature, t is a time, k is a thermal diffusivity, and x, y and z obey a Cartesian coordinate system; and a boundary condition of an irradiated surface is:
−λ.Math.∂T/∂z=q(h,t)  (3) wherein, λ is a thermal conductivity, and another boundary condition is adiabatic surface.

4. The numerical simulation method of pulsed laser paint removal according to claim 1, wherein the iterative solution is:
{T.sub.n+1}={T.sub.n}+(1−θ)t{{dot over (T)}.sub.n}+θt{{dot over (T)}.sub.n+1}  (4) wherein, {T.sub.n} is a temperature vector at a current time, {T.sub.n+1} is a temperature vector at a next time, θ is an Euler parameter, t is a time step, and {{dot over (T)}.sub.n} and {{dot over (T)}.sub.n+1} are change rates of the temperature vectors at the current and next times, respectively; and according to Eq. (4), a temperature change curve of the material of the paint layer at different depths in a heat conduction process is derived.

5. The numerical simulation method of pulsed laser paint removal according to claim 4, wherein according to the temperature change curve of the material of the paint layer at different depths in the heat conduction process, a time for heat transferred to the paint layer to reach a vaporization temperature at a maximum depth is determined as a time to perform an element birth/death operation.

6. The numerical simulation method of pulsed laser paint removal according to claim 1, wherein the coordinate system is moved by programming with APDL language to control the loading position; a movement path of the origin of the coordinate system is taken as a movement path for loading the pulsed laser; and steps 3 to 7 are repeated for each movement, until an entire laser paint removal process is completed, to obtain simulated morphology after the removal of the paint layer by the pulsed laser loaded on a specific path.

7. A use of the numerical simulation method of pulsed laser paint removal according to claim 6, wherein the use is designed to determine a parameter, a processing trajectory and a scanning rate of a pulsed laser to remove a paint layer, according to surface morphology simulated by loading pulsed lasers with different parameters on different paths at different scanning rates.

8. The use of the numerical simulation method of pulsed laser paint removal according to claim 7, wherein according to a zone of a single-pulsed laser-ablated pit in which a paint layer is effectively removed, a selection window of the scanning rate is reduced, and an overlap rate of adjacent ablated pits is controlled to be 0-50%.

9. The use of the numerical simulation method of pulsed laser paint removal according to claim 7, wherein single-direction scanning is conducted track by track on a movement path for loading the pulsed laser, and adjacent scanning tracks are arranged alternately.

10. The use of the numerical simulation method of pulsed laser paint removal according to claim 7, wherein a scanning interval between adjacent scanning tracks is determined based on a zone where the paint layer is effectively removed according to the simulated morphology obtained by single-track laser irradiation on a surface of the paint layer, and an overlap rate of adjacent ablated trenches is controlled to be 0-50%.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flowchart of a numerical simulation method of pulsed laser paint removal according to the present disclosure;

(2) FIG. 2A and FIG. 2B show a node attached with elements, where FIG. 2A shows an internal node with 8 surviving elements attached; and FIG. 2B shows a node with 1 to 7 surviving elements attached; and the figure only lists some possible situations, not all of them;

(3) FIG. 3 shows a temperature change curve of a paint layer with time at different depths under a laser power of 100 W, where 170° C. in the figure is a vaporization temperature of the paint layer;

(4) (a) and (b) of FIG. 4 respectively show a section simulated by single-point single-pulsed laser ablation and single-point double-pulsed laser ablation, where a color indicates a residual temperature;

(5) FIG. 5 shows a loading position that goes downward as a material is removed;

(6) FIG. 6 shows a simulation result of a single trench after laser paint removal; and

(7) FIG. 7 shows a simulation result of laser paint removal with different scanning intervals, where different colors indicate different depths.

REFERENCE NUMERALS

(8) 1. node; and 2. element.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(9) The present disclosure will be further described in detail below with reference to the drawings and examples, but the protection scope of the present disclosure is not limited thereto.

(10) The present disclosure provides a numerical simulation method of pulsed laser paint removal, which performs a numerical simulation of nanosecond pulsed laser paint removal based on an ANSYS temperature field, as shown in FIG. 1. A high-speed moving pulsed laser is loaded on a surface of a model in a form of heat flux, and a coordinate system is moved to realize loading on different paths. A special surface mesh screening method is used to realize loading on any surface in circumstances such as material removal. In addition, an element birth/death technology is combined to remove an element that exceeds a threshold, so as to predict the surface morphology after laser paint removal.

(11) Taking nanosecond pulsed laser paint removal with a spot diameter of 50 μm and a laser energy of 100 W as an example, the simulation process of this method is described in detail below. The simulation process includes the following steps:

(12) Step 1: A two-layer three-dimensional (3D) solid model of a paint layer and a substrate is established, where an upper layer is defined by epoxy paint, and a lower layer is defined by 304 stainless steel. In the simulation, a density is set to 1,062 and 7,750 kg/m.sup.3, specific heat is set to 2,153 and 461 J.Math.(kg.Math.° C.).sup.−1, and a thermal conductivity is set to 0.3 and 16 W.Math.(m.Math.° C.).sup.−1, respectively. Meshing is conducted. In the 3D solid model, an element type is selected as an 8-node hexahedral thermal element solid 70. The element has 8 nodes, and each node has a single degree of freedom (DOF). It can realize the transfer of uniform heat flow, and can be used for 3D static or transient thermal analysis.

(13) Step 2: An initial condition of the temperature field model is set to an ambient temperature of 20° C., and a load option is set to transient analysis, where each load step includes 20 load substeps, and a load mode is set to step load.

(14) Step 3: A coordinate system is established with a center of a single-pulse spot on an upper surface of the paint layer as an origin, an axial incident direction of a laser is taken as a z-axis, the origin of the coordinate system is taken as a pulsed laser loading position, and the upper surface of the paint layer is defined as a z=0 plane.

(15) Step 4: Nodes having a perpendicular distance not exceeding a radius (25 μm) of the spot from the z-axis in the 3D solid model are selected to form a laser irradiation coverage zone, and nodes with less than 8 surviving elements attached are selected among the nodes in the laser irradiation coverage zone to form a laser irradiated surface. An internal node of the model is generally attached with 8 elements, as shown in FIG. 2A. A load located on the upper surface of the paint layer and a surface of an ablated pit is generally attached with less than 8 surviving elements, as shown in FIG. 2B. The figure only lists some possible situations, not all of them.

(16) Step 5: The pulsed laser is loaded to a surface of the model in a form of heat flux. Elements attached to the nodes on the laser radiation energy surface are laser irradiated elements. This method realizes loading on any surface in circumstances such as material removal, and ensures that the laser energy distribution on the material surface is in line with reality. An energy load q received by a node on an upper surface of a laser irradiated element whose center has a distance of h from the z-axis is calculated:

(17) $\begin{array}{cc}q\left(h,t\right)=\frac{12.Math.\mathrm{aP}}{f\tau .Math.\pi d^2}.Math.\exp \left(-\frac{12.Math.h^2}{d^2}\right).Math.{\left(\frac{t}{\tau }\right)}^7\exp \left[7\left(1-\frac{t}{\tau }\right)\right]& \left(1\right)\end{array}$

(18) where, a is an absorption rate of laser energy by a material, which is 0.3; P is a laser output power, which is 100 W; f is a laser repetition frequency, which is 100 kHz; d is a spot diameter, which is 50 μm; r is a laser pulse width, which is 100 ns; and t is a time, 0<t<τ.

(19) Step 6: Iterative solution is conducted on the model in a form of integration according to a heat conservation law of a heat transfer theory, to calculate an instantaneous temperature of each of the nodes.

(20) The heat conservation law of the heat transfer theory is that a temperature change of the material subjected to laser irradiation satisfies:

(21) $\begin{array}{cc}\frac{\partial T}{\partial t}=k\left(\frac{{\partial }^{2}T}{\partial x^2}+\frac{{\partial }^{2}T}{\partial y^2}+\frac{{\partial }^{2}T}{\partial z^2}\right)& \left(2\right)\end{array}$

(22) where, T is a temperature, t is a time, k is a thermal diffusivity, and x, y and z obey a Cartesian coordinate system; and a boundary condition of an irradiated surface is:
—λ∂T/∂z=q(h,t)  (3)

(23) where, λ is a thermal conductivity, and another boundary condition is adiabatic surface.

(24) The iterative solution is:
{T.sub.n+1}={T.sub.n}+(1−θ)t{{dot over (T)}.sub.n}+θt{{dot over (T)}.sub.n+1}  (4)

(25) where, {T.sub.n} is a temperature vector at a current time, {T.sub.n+1} is a temperature vector at a next time, θ is an Euler parameter, t is a time step, and {{dot over (T)}.sub.n} and {{dot over (T)}.sub.n+1} are change rates of the temperature vectors at the current and next times, respectively; and according to Eq. (4), a temperature change curve of the material of the paint layer at different depths in a heat conduction process is derived, as shown in FIG. 3.

(26) Step 7: An instantaneous temperature of the surviving elements reaching a threshold temperature or above is defined as a condition for an element birth/death operation, and a surviving element whose instantaneous temperature exceeds the threshold temperature is killed to obtain morphology after the removal of the paint layer by a single pulsed laser. According to the temperature change curve of the material of the paint layer at different depths in the heat conduction process, a time for heat transferred to the paint layer to reach a vaporization temperature at a maximum depth is determined as a time to perform an element birth/death operation.

(27) FIG. 3 shows a temperature change of the paint layer at different depths under single pulse ablation. A maximum depth of the paint layer to reach the vaporization temperature is about 30 μm, and the time is about 1.1 μs. Therefore, it is determined that the element birth/death operation is conducted at about 1.1 us after loading. By killing the element exceeding a boiling point of epoxy paint, that is, 170° C., the morphology after the removal of the paint layer by the single pulsed laser is obtained, as shown in (a) of FIG. 4.

(28) Without moving the coordinate system, Steps 3 to 7 are repeated to achieve single-point double-pulsed laser ablation. The second loading is shown in FIG. 5, where the laser is loaded on a surface of a pit after material removal instead of the original paint layer. The two pulse loads are solved, and a simulated section of the paint layer ablated by single-point double-pulsed laser ablation is shown in (b) of FIG. 4.

(29) The coordinate system is moved by programming with APDL language to control the loading position. A movement path of the origin of the coordinate system is taken as a movement path for loading the pulsed laser. Steps 3 to 7 are repeated for each movement, until an entire laser paint removal process is completed, to obtain simulated morphology after the removal of the paint layer by the pulsed laser loaded on a specific path.

(30) The movement of the coordinate system is controlled to control the movement of the spot, such that the spot moves along a linear path at 5,000 mm/s. Steps 3 to 7 are repeated until all the loads are solved, and a simulation result of laser paint removal of a single trench is obtained. The surface morphology after paint removal is shown in FIG. 6.

(31) The numerical simulation method of pulsed laser paint removal in the present disclosure can determine a parameter, a processing trajectory and a scanning rate of a pulsed laser to remove a paint layer, according to surface morphology simulated by loading pulsed lasers with different parameters on different paths at different scanning rates.

(32) Single-direction scanning is conducted track by track on a movement path for loading the pulsed laser, and adjacent scanning tracks are arranged alternately. According to a zone of a single-pulsed laser-ablated pit in which a paint layer is effectively removed, a selection window of the scanning rate is reduced, and an overlap rate of adjacent ablated pits is controlled to be 0-50%. A scanning interval between adjacent scanning tracks is determined based on the zone where the paint layer is effectively removed according to the simulated morphology obtained by single-track laser irradiation on a surface of the paint layer, and the overlap rate of adjacent ablated trenches is controlled to be 0-50%. The coordinate system is controlled to move at 5,000 mm/s for constant scanning, and the multiple scanning tracks are arranged to complete area scanning simulation of laser paint removal. Different scanning intervals can be controlled. A simulation result of laser paint removal with different scanning intervals is shown in FIG. 7.

(33) The simulation results of the surface morphology obtained by constant scanning with different scanning intervals are compared. It is intuitively found that optimal laser paint removal is achieved by single-direction scanning with a scanning interval of 70 μm.

(34) The above examples are preferred implementations of the present disclosure, but the present disclosure is not limited to the above implementations. Any obvious improvement, substitution or modification made by those skilled in the art without departing from the essence of the present disclosure should fall within the protection scope of the present disclosure.