CN111515548B - Method for optimizing laser processing scanning track of micro-curvature radius antenna - Google Patents

Method for optimizing laser processing scanning track of micro-curvature radius antenna Download PDF

Info

Publication number
CN111515548B
CN111515548B CN202010273752.0A CN202010273752A CN111515548B CN 111515548 B CN111515548 B CN 111515548B CN 202010273752 A CN202010273752 A CN 202010273752A CN 111515548 B CN111515548 B CN 111515548B
Authority
CN
China
Prior art keywords
curve
laser
scanning track
ablation depth
laser processing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010273752.0A
Other languages
Chinese (zh)
Other versions
CN111515548A (en
Inventor
马建伟
贺广智
贾振元
韩东旭
姜文文
刘巍
司立坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202010273752.0A priority Critical patent/CN111515548B/en
Publication of CN111515548A publication Critical patent/CN111515548A/en
Application granted granted Critical
Publication of CN111515548B publication Critical patent/CN111515548B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00—Working by laser beam, e.g. welding, cutting or boring
    • B23K26/36—Removing material
    • B23K26/362—Laser etching
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00—Working by laser beam, e.g. welding, cutting or boring
    • B23K26/08—Devices involving relative movement between laser beam and workpiece
    • B23K26/082—Scanning systems, i.e. devices involving movement of the laser beam relative to the laser head
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
    • B23K—SOLDERING OR UNSOLDERING; WELDING; CLADDING OR PLATING BY SOLDERING OR WELDING; CUTTING BY APPLYING HEAT LOCALLY, e.g. FLAME CUTTING; WORKING BY LASER BEAM
    • B23K26/00—Working by laser beam, e.g. welding, cutting or boring
    • B23K26/70—Auxiliary operations or equipment
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Data Mining & Analysis (AREA)
  • Plasma & Fusion (AREA)
  • Mechanical Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Laser Beam Processing (AREA)

Abstract

The invention discloses a method for optimizing a laser processing scanning track of a micro-curvature radius antenna, belongs to the field of special processing, and relates to the method for optimizing the laser processing scanning track of the micro-curvature radius antenna, which considers the problem of asymmetrical erosion groove profile. The method is based on a nanosecond laser processing complex pattern part surface energy dynamic distribution model, the position of the maximum value of the surface energy density of the part is solved, and the offset of the position of the maximum ablation depth is calculated. And (4) calculating an actual maximum ablation depth curve according to an ideal antenna strip line design curve Frenet frame equation, and solving the optimized laser scanning track. And solving a complex scanning track optimization curve based on an osculating circle discrete approximation strategy according to the circular pattern laser scanning track optimization result, and realizing the laser processing scanning track optimization of the antenna with the micro-curvature radius. The method is effective and reliable, can be applied to the optimization of the laser processing scanning track of the high-performance antenna of the aircraft with the micro-curvature radius characteristic, and has important practical application significance for improving the service performance of the antenna.

Description

Method for optimizing laser processing scanning track of micro-curvature radius antenna
Technical Field
The invention belongs to the field of special processing, and relates to a micro-curvature radius antenna laser processing scanning track optimization method considering the problem of asymmetrical erosion groove profile.
Background
The high-performance antenna is a key device of high-end equipment of an aircraft. The antenna takes a polyimide shell coated with copper on the surface as a blank, a symmetrical micro-groove with a Gaussian cross section is formed by ablating the copper-coated layer through nanosecond multi-pulse laser, and a copper antenna strip line with high precision and complex geometric characteristics is formed on the surface layer of the micro-groove. In order to meet the requirement for improving the performance of an aircraft, the antenna is often designed to have a local micro-curvature radius strip line structure, and the accurate forming of the local micro-curvature radius strip line is the key for ensuring the electrical performance of the high-performance antenna, and gradually becomes a hotspot for researching the laser processing of the antenna. When the laser etching straight groove is machined, the light spot overlapping rule depends on the motion parameters of a laser light source and a machine tool, the light spots are uniformly distributed along a scanning straight line track, the accumulated amount of the surface energy of the part is symmetrical about the center of the light spot in a laser radiation area, and the straight groove has the characteristic of Gaussian symmetry of the profile of the cross section. However, during laser processing of a local micro-curvature radius strip line, laser spots are projected to the surface of a part along a complex scanning track, the highest position of the overlapping rate of the spots deviates from the center of the spots to the center of curvature of the scanning track, namely the center of the laser spots is not the position of the maximum energy accumulation of laser, so that the symmetry of the laser energy accumulation on the surface of the part during laser processing is lost, the position deviation of the maximum ablation depth of a laser erosion groove is induced, the curve of the maximum ablation depth deviates from the design curve of an ideal antenna strip line, and the performance of the antenna is directly influenced. Therefore, the method considers the asymmetric profile problem of the micro-groove of the micro-curvature radius antenna in laser processing, and the optimization of the laser processing scanning track is a key link for realizing the nanosecond multi-pulse laser precision processing of the micro-curvature radius antenna, and has important practical application significance for improving the processing quality of the antenna and the service performance of an aircraft.
Prior art document 1, "a laser processing correction method", of patent publication No. CN101508055, to Zhannian army et al, processes a part of a predetermined size using a laser processing apparatus, draws a spline curve according to a deviation value of a processing result, corrects a processing parameter, and re-processes, and repeats the above processes until a laser processing precision satisfies a technical requirement of the part. The method is weak in theoretical basis, the test result is limited to the specific processing working condition of the material, and the method is lack of guiding significance for the accurate laser processing of the characteristic pattern with small curvature radius. Technical literature 2 "a proproach to minimize build errors in direct metal laser sintering", y.ning et al, IEEE Transactions on Automation Science and Engineering, 2006, 3 (1): 73-80 which takes into account the effect of the geometry of the part on the metal laser sintering (DMLS) accuracy, optimizes the laser scanning speed by compensating for the effect of the geometry on DMLS, and improves the dimensional accuracy of the part. The method explores the influence of the geometric shape of the part on the processing precision, however, for the laser processing rule of the part with complex geometric characteristics, the proposed scanning speed optimization method cannot fundamentally solve the problem that the laser processing quality of the micro-curvature radius pattern cannot meet the requirement.
Disclosure of Invention
Aiming at the limitations and defects of the prior art, the invention provides a method for optimizing the laser processing scanning track of a micro-curvature radius antenna. The method considers the relevance of the geometric characteristics of the laser erosion groove profile and the curvature radius of a scanning track in the nanosecond multi-pulse laser processing process, solves the position of the maximum value of the surface energy density of the part based on a nanosecond multi-pulse laser processing complex pattern part surface energy dynamic distribution model, and calculates the offset of the position of the maximum ablation depth; then, calculating an actual maximum ablation depth curve by relying on Frenet frame equations at all positions of an ideal antenna strip line design curve, and solving an optimized laser processing scanning track through a second-order partial differential equation set; and finally, according to the optimization result of the laser scanning track of the circular pattern, solving a complex scanning track optimization curve based on an osculating circle discrete approximation strategy, and realizing the optimization of the laser processing scanning track of the antenna with the micro-curvature radius. The method effectively and reliably solves the problem that the maximum ablation depth position generated when the nanosecond multi-pulse laser is used for processing the micro-curvature radius antenna deviates from the ideal antenna strip line design position, can be applied to the optimization of the laser processing scanning track of the high-performance antenna of the aircraft with the micro-curvature radius characteristic, and has important practical application significance for improving the processing quality of the antenna and the service performance of the aircraft.
The technical scheme adopted by the invention is a micro-curvature radius antenna laser processing scanning track optimization method which is characterized in that the method considers the correlation between the geometrical characteristics of a laser erosion groove profile and the curvature radius of a scanning track in the nanosecond multi-pulse laser processing process, solves the position of the maximum energy density of the surface of a part based on a nanosecond multi-pulse laser processing complex pattern part surface energy dynamic distribution model, and calculates the offset of the position of the maximum ablation depth; then, calculating an actual maximum ablation depth curve according to Frenet frame equations at all positions of an ideal antenna strip line design curve, and solving an optimized laser processing scanning track through a second-order partial differential equation set; and finally, according to the optimization result of the laser scanning track of the circular pattern, solving a complex scanning track optimization curve based on an osculating circle discrete approximation strategy, and realizing the optimization of the laser processing scanning track of the antenna with the micro-curvature radius.
The method comprises the following specific steps:
step 1 solving the offset of the maximum ablation depth position
The distribution formula of the nanosecond multi-pulse laser single-pulse energy density on the beam waist plane is as follows:
Figure BDA0002444049940000031
wherein, w0Is the radius of the girdling, F0And x and y are coordinates of points on the beam waist plane. When nanosecond multi-pulse laser is processed along a circular arc track with the radius of R, a dynamic distribution model of the surface energy of the target material is as follows:
Figure BDA0002444049940000032
wherein f is the laser repetition frequency, v is the laser scanning speed, k is the relative position of the facula, and the length R from any point in the laser scanning area on the surface of the target material to the center of the circular arc and the radius R of the circular arc satisfy:
r＝R+kw (3)
as can be seen from the monotonicity law of equation (2), equation (2) has a unique maximum value point in the whole scanning area of the laser, that is:
Figure BDA0002444049940000041
wherein i is R and w0The ratio of (A) to (B):
Figure BDA0002444049940000042
as can be seen from the formula (2), FkIs not always 0, and must be guaranteed if formula (4) has a solution
Figure BDA0002444049940000043
The theoretical solution of equation (4) is therefore:
Figure BDA0002444049940000044
furthermore, in the antenna design process, in order to ensure the processability of the antenna surface, the antenna strip line geometric characteristics are generally constrained as follows:
i≥2 (7)
as can be seen from equations (6) and (7),
Figure BDA0002444049940000045
so that the relative position k of the maximum of the energy density is within the whole laser scanning areadComprises the following steps:
Figure BDA0002444049940000046
thus, the offset l of the location of maximum ablation depthdCan be expressed as:
Figure BDA0002444049940000047
step 2: laser processing scanning track after optimization based on Frenet frame solution
Assuming that the ideal antenna strip line design curve is r (r)(s), the Frenet frame equation { r(s): α(s), β(s) } is expressed as:
Figure BDA0002444049940000051
Figure BDA0002444049940000052
Figure BDA0002444049940000053
wherein s is the curve arc length, α(s) is the curve tangent vector, and β(s) is the curve normal vector. The curvature κ(s) of the curve can be obtained from equation (12). t is characteristic of arc length s, then s is expressed as:
s＝s(t) (13)
Figure BDA0002444049940000054
thus, equations (10-12) can be further expressed as:
Figure BDA0002444049940000055
Figure BDA0002444049940000056
Figure BDA0002444049940000057
when the laser is processed along the curve r (t), the actual maximum ablation depth curve r is obtained by considering the maximum ablation depth position offset in the step 1d(t) deviating the normal vector direction β (t) from the ideal curve r (t) at each position along the curve r (t) by a distance:
Figure BDA0002444049940000058
thus, the actual maximum ablation depth curve is:
rd(t)＝r(t)+ld(t)·β(t) (19)
therefore, the laser machining scan trajectory must be optimized to ensure that the actual laser maximum ablation depth position is closer to the ideal position.
Let the optimized laser processing scanning track be ro＝ro(t) curve curvature κo(t) of (d). Similarly, when the laser ablates the target along the track, the maximum ablation depth curve also follows the normal vector direction βo(t) deviating from the scanning track of laser processing by the following position deviation:
Figure BDA0002444049940000061
the actual maximum ablation depth curve of the scan trajectory optimization post-processing should be:
r(t)＝ro(t)+lo(t)·βo(t) (21)
therefore, the optimized laser processing scanning track r can be solved by using a second-order partial differential equation systemo(t)：
Figure BDA0002444049940000062
And step 3: method for solving complex scanning track optimization curve based on osculating circle discrete approximation strategy
Step 2 provides a method for solving the optimized laser processing scanning track by using a second-order partial differential equation set, however, if an ideal antenna strip line design curve r is r (t) and an analytic expression is complex or has no analytic expression, the solving process of the equation set is complicated and difficult, and the method is difficult to apply to the optimization problem of the laser processing scanning track of the antenna with actual complex geometric characteristics.
Therefore, the idea of discrete numerical approximation is adopted, the complex curve is divided into discrete circular arcs for processing, the problem of solving the second-order partial differential equation system of the complex curve is solved based on the motion of the pulse laser on the vector circle scanning track and the position deviation rule of the maximum ablation depth curve, and the complex scanning track optimization method based on the osculating circle discrete approximation strategy is established.
For a circular pattern of radius R, R (θ), its Frenet frame equation { R (θ): α (θ), β (θ) } can be expressed as:
α(θ)＝e1(θ) (23)
β(θ)＝e(θ) (24)
Figure BDA0002444049940000071
wherein e is1(θ) is a unit tangent vector, e (θ) is a unit normal vector, and κ (θ) is a curvature.
Substituting the formulas (23-25) into the formulas (15-19), and scanning the laser along the ideal circle curve to obtain the actual maximum ablation depth circle rd(θ)＝RdComprises the following steps:
Figure BDA0002444049940000072
substituting the formulas (23-25) into the formula (22) to obtain the optimized laser scanning circle ro(θ)＝RoTo solve the system of equations:
Figure BDA0002444049940000073
solving the system of equations yields:
Figure BDA0002444049940000074
therefore, when the laser is scanned along the curve r (t), the laser spot is dispersed into the point P based on the pulse laser spot distribution point on the curve r (t)1To PnThe osculating circle of the curve at each discrete point is o.c.i. Calculating each discrete point P according to equation (26)iThe actual maximum ablation depth position is offset distance along the normal vector direction of the osculating circle, and the actual maximum ablation depth position corresponding to each discrete point isPi', from Pi' fitting the actual maximum ablation depth curve, which should theoretically match the actual maximum ablation depth curve r during laser machiningd(t) complete anastomosis.
The optimized laser processing scanning track calculation method is similar to the actual maximum ablation depth curve calculation method, and pulse laser spot distribution points on the curve r (t) are dispersed into points P based on the pulse laser spot distribution points1To PnThe osculating circle of the curve at each discrete point is o.c.i. Calculating each discrete point P according to equation (28)iPosition P of discrete point of optimized laser processing scanning track corresponding to direction of normal vector of osculating circlei", from Pi' laser processing scanning track r after fitting optimizationo(t) of (d). Theoretically, when the laser light follows the optimized scanning track ro(t) during processing, the maximum ablation depth curve coincides with the ideal antenna strip design curve r (t).
The invention has the following remarkable effects and benefits: the method is a micro-curvature radius antenna laser processing scanning track optimization method considering the problem of asymmetrical profile of an etched groove, clarifies the position deviation rule of a laser processing maximum ablation depth curve, and realizes the solution of the actual maximum ablation depth curve position when the micro-curvature radius antenna is processed by laser; designing a curve Frenet frame equation based on an ideal antenna strip line, calculating an actual maximum ablation depth curve, and solving an optimized laser processing scanning track through a second-order partial differential equation set; a complex scanning track optimization method based on an osculating circle discrete approximation strategy is established, and nanosecond multi-pulse laser precision machining of a micro-curvature radius antenna is realized. The method can effectively and reliably solve the problem that the maximum ablation depth position generated when the nanosecond multi-pulse laser is used for processing the micro-curvature radius antenna deviates from the ideal antenna strip line design position, can be applied to the optimization of the laser processing scanning track of the high-performance antenna of the aircraft with the micro-curvature radius characteristic, and has important practical application significance for improving the processing quality of the antenna and the service performance of the aircraft.
Drawings
FIG. 1-flowchart of the overall process.
FIG. 2-geometric model of the rose trefoil lineAnd Frenet frame schematic diagram. Wherein r (theta) is a rose trefoil line, e1Is a unit tangent vector, e is a unit normal vector, P1To P7Seven theoretical points are respectively arranged on the single leaf at the interval of 6 degrees.
FIG. 3-results of nanosecond pulsed laser ablation of the three-leaf rose lines, D1To D7Respectively is a measuring benchmark at a position 200 mu m away from the normal vector direction of a theoretical point on a single leaf.
FIG. 4 is a schematic diagram of measurement of offset of maximum ablation depth of nanosecond pulsed laser with three leaves and roses.
FIG. 5 shows the results of the actual maximum ablation depth offset and the optimized scan trajectory maximum ablation depth offset test at different positions and processing parameters. Wherein, figure a) is P3The measurement results are shown in b) as P4The measurement results are shown in FIG. c) as P5And measuring the result. In each figure: line L1-a calculated theoretical maximum ablation depth offset; line L2Scanning speeds of 2m/min, 3.5m/min and 5m/min, a laser repetition frequency of 30kHz, and a laser energy density of 4kJ/m2Average maximum ablation depth offset under conditions, line L3-optimizing the average maximum ablation depth offset for the same parameter after scanning the track; line L4Laser energy density of 2kJ/m2、3kJ/m2And 4kJ/m2Average maximum ablation depth offset at a laser repetition frequency of 30kHz and a scanning speed of 3m/min, line L5-optimizing the average maximum ablation depth offset for the same parameter after scanning the track; line L6Repetition frequencies of 20kHz, 30kHz and 40kHz and laser energy density of 4kJ/m2Average maximum ablation depth offset at a scanning speed of 3m/min, line L7-optimizing the average maximum ablation depth offset for the same parameter after scanning the track.
Detailed Description
The detailed description of the embodiments of the invention is provided with reference to the accompanying drawings.
In the process of ablating the antenna with the micro-curvature radius characteristic by the nanosecond-level multi-pulse laser, the maximum ablation depth curve deviates from an ideal curve due to the asymmetrical ablation micro-groove profile caused by uneven energy distribution of the nanosecond-level pulse laser, and the performance of the high-performance antenna is directly influenced. Aiming at the limitations and defects of the prior art, the invention provides a micro-curvature radius antenna laser processing scanning track optimization method considering the problem of asymmetrical erosion groove profile, and the whole flow is shown in figure 1.
Example with waist radius w0Nanosecond multi-pulse laser with the wavelength of 532nm and the ablation diameter r of 20 mu m on the copper target0The specific solving process of the method is explained in detail by using MATLAB software and a verification test as an example of a 500-mu m trilobe rose line.
The first step solves for the offset of the maximum ablation depth curve position:
from equation (1), the waist radius w0The energy distribution formula of the nanosecond multi-pulse laser with the diameter of 20 mu m on the beam waist plane is as follows:
Figure BDA0002444049940000101
the polar equation and the curvature calculation formula of the three-leaf rose line are as follows:
r(θ)＝r0×cos(3θ) (29)
Figure BDA0002444049940000102
the geometric model of the three-leaf rose line and the Frenet frame are shown in the attached figure 2, wherein the value range of theta is set to be [0 degrees, 180 degrees ].
According to the formula (2), when the nanosecond laser processes the target along the arc track with the radius of R, the dynamic distribution model of the surface energy of the target is as follows:
Figure BDA0002444049940000103
calculating the maximum ablation depth position deviation l according to the formulas (3) to (9)dObtaining the position of the curve of maximum ablation depthOffset amount:
Figure BDA0002444049940000111
and secondly, solving the optimized laser processing scanning track based on a Frenet standard frame. Obtaining the laser processing scanning track r after the solution and optimization according to the formula (22)oThe second order differential equation set of (θ) and its solution:
Figure BDA0002444049940000112
Figure BDA0002444049940000113
and thirdly, solving a complex scanning track optimization curve based on an osculating circle discrete approximation strategy.
Fitting the actual maximum ablation depth curve r according to equation (26) and equation (28)d(theta) and optimized laser machining scan trajectory ro(θ)。
The change of the curvature and the curvature radius of the three-leaf rose line is periodic, and the experiment verifies that only one leaf in the three-leaf rose line is selected. Seven points are taken at equal angular intervals on the leaf, see fig. 2. The radius of curvature, the maximum ablation depth offset, and the compensation for the optimized scan trajectory at these points were calculated separately and the results are shown in table 1.
TABLE 1 theoretical calculation results of laser processing of three-leaf rose lines
Figure BDA0002444049940000114
Figure BDA0002444049940000121
FIG. 3 shows the results of laser ablation with three-leaf rose lines for nanosecond pulses to measure the maximum ablation depth shift of laser ablation at each locationDistance and optimization effect, D1To D7Respectively is a measuring benchmark at a position 200 mu m away from the normal vector direction of a theoretical point on a single leaf. In order to verify the effectiveness of the laser processing scanning track optimization result, a nanosecond pulse ablation result of the trilobed rose line is measured by using a super-depth-of-field three-dimensional microscope system VHX-600E, as shown in the attached figure 4.
Selecting the point P with the maximum ablation depth offset3、P4And P5And evaluating the optimization effect of the laser scanning track. P3、P4And P5The position offset of the theoretical maximum ablation depth curve of the point under different processing parameters and the offset comparison result of the actual maximum ablation depth curve compared with the ideal curve after the scanning track of the optimized laser processing are respectively shown in the attached figures 5a), b) and c). The test results show that the laser processing is carried out along the ideal three-leaf rose line, P3、P4And P5Are 0.47, 2.18, and 0.51 μm, respectively; correspondingly, the machining is carried out along the optimized laser scanning track under the same machining parameter setting, P3、P4And P5The average offsets of the maximum ablation depth position of (a) from the ideal trilobed line are 0.10 μm, 0.51 μm and 0.13 μm, respectively. Point P for optimized laser machining results compared to non-optimized laser machining results3、P4And P5The average maximum ablation depth offset is reduced by 78.72%, 76.61%, and 74.51%, respectively.
In conclusion, the theoretical prediction result is well consistent with the test result, and the method can effectively and reliably solve the problem that the maximum ablation depth position generated when the nanosecond multi-pulse laser processing micro-curvature radius antenna deviates from the ideal antenna strip line design position, can be applied to the laser processing scanning track optimization of the high-performance antenna of the high-speed aircraft with the micro-curvature radius characteristic, is beneficial to realizing the precision processing of the high-performance antenna of the high-speed aircraft, and has important practical application significance for improving the processing quality of the antenna and improving the service performance of the high-speed aircraft.

Claims (1)

1. A micro-curvature radius antenna laser processing scanning track optimization method is characterized in that relevance of laser erosion groove profile geometric characteristics and scanning track curvature radius in a nanosecond multi-pulse laser processing process is considered, based on a nanosecond multi-pulse laser processing complex pattern part surface energy dynamic distribution model, a part surface energy density maximum value position is solved, and maximum ablation depth position offset is calculated; then, calculating an actual maximum ablation depth curve according to Frenet frame equations at all positions of an ideal antenna strip line design curve, and solving an optimized laser processing scanning track through a second-order partial differential equation system; finally, according to the optimization result of the laser scanning track of the circular pattern, a complex scanning track optimization curve based on an osculating circle discrete approximation strategy is solved, and the optimization of the laser processing scanning track of the antenna with the micro-curvature radius is realized; the method comprises the following specific steps:
step 1: solving for the offset of the maximum ablation depth position
The distribution formula of the nanosecond multi-pulse laser single-pulse energy density on the beam waist plane is as follows:
Figure FDA0002988237290000011
wherein, w0Is the radius of the girdling, F0The laser energy density is shown, and x and y are coordinates of points on a beam waist plane; when nanosecond multi-pulse laser is processed along a circular arc track with the radius of R, a dynamic distribution model of the surface energy of the target material is as follows:
Figure FDA0002988237290000012
wherein f is the laser repetition frequency, v is the laser scanning speed, k is the relative position of the facula, and the length R from any point in the laser scanning area on the surface of the target material to the center of the circular arc and the radius R of the circular arc satisfy:
r＝R+kw0 (3)
as can be seen from the monotonicity law of equation (2), equation (2) has a unique maximum value point in the whole scanning area of the laser, that is:
Figure FDA0002988237290000021
wherein i is R and w0The ratio of (A) to (B):
Figure FDA0002988237290000022
as can be seen from the formula (2), FkIs not always 0, and must be guaranteed if formula (4) has a solution
Figure FDA0002988237290000023
The theoretical solution of equation (4) is therefore:
Figure FDA0002988237290000024
furthermore, in the antenna design process, in order to ensure the processability of the antenna surface, the antenna strip line geometric characteristics are generally constrained as follows:
i≥2 (7)
from equations (6) and (7):
Figure FDA0002988237290000025
so that the relative position k of the maximum of the energy density is within the whole laser scanning areadComprises the following steps:
Figure FDA0002988237290000026
thus, the offset l of the location of maximum ablation depthdExpressed as:
Figure FDA0002988237290000027
step 2: laser processing scanning track after optimization based on Frenet frame solution
Assuming that the ideal antenna strip line design curve is r (r)(s), the Frenet frame equation { r(s): α(s), β(s) } is expressed as:
Figure FDA0002988237290000031
Figure FDA0002988237290000032
Figure FDA0002988237290000033
wherein s is the curve arc length, alpha(s) is the curve tangent vector, and beta(s) is the curve normal vector; the curve curvature κ(s) is obtained from equation (12); t is characteristic of arc length s, then s is expressed as:
s＝s(t) (13)
Figure FDA0002988237290000034
accordingly, equations (10) - (12) are further expressed as:
Figure FDA0002988237290000035
Figure FDA0002988237290000036
Figure FDA0002988237290000037
when the laser is usedWhen processing along the curve r (t), the actual maximum ablation depth curve r is obtained by considering the maximum ablation depth position offset obtained in step 1d(t) deviating the normal vector direction β (t) from the ideal curve r (t) at each position along the curve r (t) by a distance:
Figure FDA0002988237290000038
thus, the actual maximum ablation depth curve is:
rd(t)＝r(t)+ld(t)·β(t) (19)
therefore, the laser processing scanning track must be optimized to ensure that the actual laser maximum ablation depth position is closer to the ideal position;
let the optimized laser processing scanning track be ro＝ro(t) curve curvature κo(t); similarly, when the laser ablates the target along the track, the maximum ablation depth curve also follows the normal vector direction βo(t) deviating from the scanning track of laser processing by the following position deviation:
Figure FDA0002988237290000041
the actual maximum ablation depth curve of the scan trajectory optimization post-processing should be:
r(t)＝ro(t)+lo(t)·βo(t) (21)
therefore, the optimized laser processing scanning track r is solved by a second-order partial differential equation systemo(t)：
Figure FDA0002988237290000042
And step 3: method for solving complex scanning track optimization curve based on osculating circle discrete approximation strategy
Step 2, a method for solving the optimized laser processing scanning track by using a second-order partial differential equation set is provided, however, if an ideal antenna strip line design curve r is r (t) and an analytic expression is complex or has no analytic expression, the solving process of the equation set is complicated and difficult, and the method is difficult to apply to the optimization problem of the laser processing scanning track of the antenna with actual complex geometric characteristics;
therefore, the idea of discrete numerical approximation is adopted, the complex curve is divided into discrete circular arcs for processing, the problem of solving a second-order partial differential equation system of the complex curve is solved based on the motion of pulse laser on a vector circle scanning track and the position deviation rule of a maximum ablation depth curve, and a complex scanning track optimization method based on a osculating circle discrete approximation strategy is established;
for a circular pattern of radius R, R (θ), the Frenet frame equation { R (θ): α (θ), β (θ) } is expressed as:
α(θ)＝e1(θ) (23)
β(θ)＝e(θ) (24)
Figure FDA0002988237290000051
wherein e is1(theta) is a unit tangent vector, e (theta) is a unit normal vector, and kappa (theta) is a curvature;
substituting equations (23) - (25) into equations (15) - (19), the laser scans the actual maximum ablation depth circle r obtained during processing along the ideal circle curved(θ)＝RdComprises the following steps:
Figure FDA0002988237290000052
substituting equations (23) - (25) into equation (22) to obtain the optimized laser scanning circle ro(θ)＝RoTo solve the system of equations:
Figure FDA0002988237290000053
solving the system of equations yields:
Figure FDA0002988237290000054
therefore, when the laser is scanned along the curve r (t), the laser spot is dispersed into the point P based on the pulse laser spot distribution point on the curve r (t)1To PnThe osculating circle of the curve at each discrete point is o.c.i(ii) a Calculating each discrete point P according to equation (26)iThe actual maximum ablation depth position is the offset distance along the tangential circle normal vector direction, and the actual maximum ablation depth position corresponding to each discrete point is Pi', from Pi' fitting the actual maximum ablation depth curve, which should theoretically match the actual maximum ablation depth curve r during laser machiningd(t) complete anastomosis;
the optimized laser processing scanning track calculation method is similar to the actual maximum ablation depth curve calculation method, and pulse laser spot distribution points on the curve r (t) are dispersed into points P based on the pulse laser spot distribution points1To PnThe osculating circle of the curve at each discrete point is o.c.i(ii) a Calculating each discrete point P according to equation (28)iPosition P of discrete point of optimized laser processing scanning track corresponding to direction of normal vector of osculating circlei", from Pi' laser processing scanning track r after fitting optimizationo(t); theoretically, when the laser light follows the optimized scanning track ro(t) during processing, the maximum ablation depth curve coincides with the ideal antenna strip design curve r (t).
CN202010273752.0A 2020-04-09 2020-04-09 Method for optimizing laser processing scanning track of micro-curvature radius antenna Active CN111515548B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010273752.0A CN111515548B (en) 2020-04-09 2020-04-09 Method for optimizing laser processing scanning track of micro-curvature radius antenna

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010273752.0A CN111515548B (en) 2020-04-09 2020-04-09 Method for optimizing laser processing scanning track of micro-curvature radius antenna

Publications (2)

Publication Number Publication Date
CN111515548A CN111515548A (en) 2020-08-11
CN111515548B true CN111515548B (en) 2021-05-18

Family

ID=71911707

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010273752.0A Active CN111515548B (en) 2020-04-09 2020-04-09 Method for optimizing laser processing scanning track of micro-curvature radius antenna

Country Status (1)

Country Link
CN (1) CN111515548B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113703399B (en) * 2021-08-26 2022-09-16 广东拓斯达科技股份有限公司 Motion trajectory planning method, device, equipment and storage medium
CN116551215B (en) * 2023-07-06 2023-09-05 北京新科以仁科技发展有限公司 Laser scanning control method, device, equipment and storage medium of laser

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106001933A (en) * 2016-05-10 2016-10-12 湖南大捷智能装备有限公司 Optimization method for laser cutting trimming line
WO2016172470A1 (en) * 2015-04-24 2016-10-27 Testrake Aviation, Inc. Optical device and method of making same
CN108955697A (en) * 2018-04-11 2018-12-07 哈尔滨工程大学 A kind of remote sensing satellite posture planing method towards multi-curvature dynamic imaging target
CN109358568A (en) * 2018-12-17 2019-02-19 大连理工大学 Curved surface subregion machining locus topology design method based on vector field
CN110919193A (en) * 2019-12-31 2020-03-27 上海工程技术大学 Nanosecond laser processing method for ceramic surface hole based on processing track optimization

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10744539B2 (en) * 2017-10-27 2020-08-18 The Boeing Company Optimized-coverage selective laser ablation systems and methods

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016172470A1 (en) * 2015-04-24 2016-10-27 Testrake Aviation, Inc. Optical device and method of making same
CN106001933A (en) * 2016-05-10 2016-10-12 湖南大捷智能装备有限公司 Optimization method for laser cutting trimming line
CN108955697A (en) * 2018-04-11 2018-12-07 哈尔滨工程大学 A kind of remote sensing satellite posture planing method towards multi-curvature dynamic imaging target
CN109358568A (en) * 2018-12-17 2019-02-19 大连理工大学 Curved surface subregion machining locus topology design method based on vector field
CN110919193A (en) * 2019-12-31 2020-03-27 上海工程技术大学 Nanosecond laser processing method for ceramic surface hole based on processing track optimization

Also Published As

Publication number Publication date
CN111515548A (en) 2020-08-11

Similar Documents

Publication Publication Date Title
CN111515548B (en) Method for optimizing laser processing scanning track of micro-curvature radius antenna
CN102189335B (en) For the manufacture of laser processing device and the method for rotation symmetric tool
CN108491352B (en) Ablation depth solving method based on laser energy dynamic distribution model
CN110497092B (en) Laser processing method of low side wall taper angle blind groove
CN102151984A (en) Laser machining method and device applicable for complicated curved surface
JP2008502485A (en) Continuous butt welding method using plasma and laser, and metal pipe manufacturing method using the same
CN106676519A (en) Method for performing laser cladding on inclined basal body
US20210331276A1 (en) A roller laser texturing processing equipment and its processing method
CN109593919A (en) Bearing surface laser-quenching apparatus and method based on the scanning of Distributed Three-dimensional light beam
CN110340485A (en) A kind of molten product method of oriented energy five axis of deposition of cantilever design
CN111575702B (en) Laser cladding method and system
CN109518180A (en) A kind of device and method of adaptive laser deposition reparation
RU2477679C2 (en) Method of repairing metal plate worn-out end force part
US11878369B1 (en) Laser scanning welding method for lap joints based on linear energy density regulation
CN211939504U (en) Spiral bevel gear femtosecond laser processing system
CN113031516A (en) Plane compensation trajectory optimization method and device considering diameter of milling cutter
CN107824963A (en) The method for improving high reflecting metal surface laser absorption rate
CN112247365B (en) Nanosecond pulse laser tilt machining ablation profile prediction method
CN110977174A (en) Pulse laser high-speed same-point interval multiple processing system and processing method
CN115890024A (en) Method for processing rectangular inclined hole by ultrafast laser
US20220241898A1 (en) Method for manufacturing machine parts, such as, but not limited to compressor, expander or vacuum pump parts and machine part manufactured by said method
CN111085779B (en) Laser modulation welding method for outer part containing process step and thin plate substrate
CN112695157A (en) Gear laser quenching method and device capable of obtaining continuous uniform hardened layer
CN107052550A (en) A kind of galvanized steel plain sheet welding method
CN116393815B (en) Laser scanning welding method considering welding inclination angle

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant