CN110598332A - Method for calculating axial position of focus of high-power laser cutting optical system - Google Patents
Method for calculating axial position of focus of high-power laser cutting optical system Download PDFInfo
- Publication number
- CN110598332A CN110598332A CN201910874987.2A CN201910874987A CN110598332A CN 110598332 A CN110598332 A CN 110598332A CN 201910874987 A CN201910874987 A CN 201910874987A CN 110598332 A CN110598332 A CN 110598332A
- Authority
- CN
- China
- Prior art keywords
- optical system
- cutting optical
- laser
- cutting
- axial position
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Laser Beam Processing (AREA)
Abstract
The invention discloses a method for calculating the axial position of a focal point of a high-power laser cutting optical system, which adopts multi-physical-field coupling simulation software to calculate multi-physical-field coupling processes such as a light-heat-stress field and the like in the cutting optical system under simulated laser irradiation, obtains coordinate data of each surface shape point in the system, guides the coordinate data of each surface shape point into optical design software, and obtains the axial position of the focal point of the cutting optical system under laser power irradiation through optical path calculation. The invention provides a method for calculating the axial position of a focal point of a laser cutting optical system under the action of high-power laser, and provides a basis for realizing high-quality and high-efficiency cutting processing.
Description
Technical Field
The invention relates to a method for calculating the axial position of a focus of a high-power laser cutting optical system, and belongs to the technical field of laser.
Background
In a high-power laser cutting system, a laser beam passes through a cutting optical system and then is focused and coupled to the surface of a workpiece to perform cutting processing. Since the optical lens in the cutting optical system absorbs the laser to generate thermal expansion deformation, the surface shape of the corresponding lens changes. With the continuous increase of laser power, the expansion deformation of the optical lens leads to the change of the focal length of each optical lens in the cutting optical system, and further causes the axial drift of the focal point of the cutting optical system, which affects the cutting effect.
Disclosure of Invention
The invention aims to provide a method for calculating the axial position of the focal point of a high-power laser cutting optical system.
In order to achieve the above purpose, as shown in fig. 1, the technical scheme adopted by the invention is as follows:
a method for calculating the axial position of a focus of a high-power laser cutting optical system mainly comprises two steps:
firstly, simulating the coordinates of each surface shape point in the cutting optical system by adopting multi-physical field coupling software, wherein the coordinates of each surface shape point in the cutting optical system are irradiated by high-power laser, and the specific steps are as follows:
step 1-1, setting physical fields such as light, heat and force, and setting the coupling working state of a plurality of physical fields as a stable state;
step 1-2, establishing a three-dimensional geometric model of the cutting optical system, and setting physical parameters and material parameters of the cutting optical system;
step 1-3, loading an irradiation laser light source A and setting parameters of the laser light source A;
step 1-4, setting initial conditions and boundary conditions of each optical lens in the cutting optical system;
step 1-5, establishing coupling relations among multiple physical fields such as light-heat-stress fields and the like;
step 1-6, performing free tetrahedral mesh subdivision on a three-dimensional geometric model of the cutting optical system;
1-7, performing coupling calculation of multiple physical fields such as a light-heat-stress field and the like;
step 1-8, outputting coordinate data of points on each surface of the cutting optical system under the irradiation of laser power P;
secondly, calculating the axial position of the focal point of the cutting optical system under the irradiation of the laser power P by adopting optical design software, and specifically comprising the following steps:
step 2-1, establishing a laser light source B; the laser light source B and the laser light source A have the same wavelength, light field distribution function, light beam divergence angle, output light spot diameter and axial interval L with the cutting optical system;
step 2-2, importing the coordinate data of each surface shape point of the cutting optical system obtained in the step 1-8 into a non-sequence module of optical design software;
and 2-3, calculating a light path trace to obtain the axial position of the focal point of the cutting optical system under the irradiation of the laser power P.
Has the advantages that: the invention adopts multi-physics coupling software to simulate the surface shape change of each optical lens in the laser cutting optical system under the high-power laser irradiation, and then introduces the point coordinates on each optical lens into optical design software to carry out light path calculation simulation, thereby obtaining the axial position of the focus of the cutting optical system under the high-power laser irradiation and providing a theoretical basis for the focus control of high-power laser cutting equipment.
Drawings
Fig. 1 is a flowchart of calculation of the axial position of the focal point of the laser cutting optical system.
Fig. 2 is a three-dimensional geometric model of a laser cutting optical system.
Detailed Description
Embodiment 1 is a method for calculating the axial position of a focal point of a high-power laser cutting optical system.
The invention discloses a method for calculating the axial position of a focus of a high-power laser cutting optical system, which is carried out according to the following steps by combining the attached drawings 1 and 2:
firstly, simulating the coordinates of each surface shape point in the laser cutting optical system by using COMSOL (common mode of integration) simulation high-power laser in a multi-physical field coupling software, and specifically comprising the following steps:
step 1-1, entering a software main interface, selecting 'model guide', setting optical-thermal-force and other physical fields, adding ray optics, solid heat transfer and solid mechanics physical field interfaces, and setting a multi-physical field coupling working state to be a stable state;
step 1-2, establishing a three-dimensional geometric model of the cutting optical system, and setting physical parameters of each optical lens in the cutting optical system; as shown in fig. 2, the cutting optical system is composed of a plano-convex lens F1 and a plano-convex lens F2 made of K9 glass, and the distance Δ L between the plano-convex lens F1 and the plano-convex lens F2 is 60 mm; the curvature radius of the plano-convex lens F1 and the plano-convex lens F2 is 50mm, the center thickness is 4mm, the diameter is 25.4mm, the heat conductivity coefficient is 1.38W/(m × K), and the density is 2203kg/m3Constant pressure heat capacity 703J/(kg K), Young's modulus 7.31X 1010Pa, Poisson's ratio of 0.17, and coefficient of thermal expansion of 0.55 x 10-61/K, refractive index of 1.5 and absorption coefficient of 0.05 percent;
step 1-3, loading an irradiation laser light source A, setting laser power P to be 2000W, setting the wavelength to be 1064nm, setting a light field distribution function to be a Gaussian function, setting a light beam divergence angle theta to be 13.8 degrees, setting the diameter of an output light spot to be 0.2mm, and setting an axial interval L between the output light spot and the cutting optical system to be 98.7 mm;
step 1-4, setting initial conditions and boundary conditions of each optical lens in the cutting optical system; for a ray optical interface, setting the emitting position of laser at the focus of a lens, wherein the emitting direction is vertical to the surface of the lens; for a solid heat transfer interface, the convective heat flux of the optical lens is set to 10W/(m)2K); for a solid mechanical interface, setting the edge part of a lens as fixed constraint; and coupling the linear optics, solid heat transfer and solid mechanics interfaces;
step 1-5, establishing coupling relations among multiple physical fields such as light-heat-stress fields and the like;
step 1-6, performing free tetrahedral mesh subdivision on a three-dimensional geometric model of the cutting optical system;
1-7, performing coupling calculation of multiple physical fields such as a light-heat-stress field and the like;
step 1-8, outputting coordinate data of points on each surface of the cutting optical system under the irradiation of laser power P;
and secondly, calculating and simulating the axial position of the focal point of the cutting optical system by adopting optical design software Zemax, wherein the method comprises the following specific steps:
step 2-1, establishing a laser light source B; the laser light source B and the laser light source A have the same wavelength, light field distribution function, light beam divergence angle, output light spot diameter and axial interval L with the cutting optical system;
step 2-2, importing the coordinate data of each point on the surface shape of the plano-convex lens F1 and the plano-convex lens F2 obtained in the step 1-5 into a non-sequence module of optical design software;
and 2-3, calculating a light path to obtain the axial position Z of the cutting optical system focus under the irradiation of laser power 2000W.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (3)
1. A method for calculating the axial position of a focus of a high-power laser cutting optical system is characterized by mainly comprising the following two steps:
firstly, simulating the coordinates of each surface shape point in the cutting optical system by adopting multi-physical field coupling software, wherein the coordinates of each surface shape point in the cutting optical system are irradiated by high-power laser, and the specific steps are as follows:
step 1-1, setting physical fields such as light, heat and force, and setting the coupling working state of a plurality of physical fields as a stable state;
step 1-2, establishing a three-dimensional geometric model of the cutting optical system, and setting physical parameters and material parameters of the cutting optical system;
step 1-3, loading an irradiation laser light source A and setting parameters of the laser light source A;
step 1-4, setting initial conditions and boundary conditions of each optical lens in the cutting optical system;
step 1-5, establishing coupling relations among multiple physical fields such as light-heat-stress fields and the like;
step 1-6, performing free tetrahedral mesh subdivision on a three-dimensional geometric model of the cutting optical system;
1-7, performing coupling calculation of multiple physical fields such as a light-heat-stress field and the like;
step 1-8, outputting coordinate data of points on each surface of the cutting optical system under the irradiation of laser power P;
secondly, calculating the axial position of the focal point of the cutting optical system under the irradiation of the laser power P by adopting optical design software, and specifically comprising the following steps:
step 2-1, establishing a laser light source B; the laser light source B and the laser light source A have the same wavelength, light field distribution function, light beam divergence angle, output light spot diameter and axial interval L with the cutting optical system;
step 2-2, importing the coordinate data of each surface shape point of the cutting optical system obtained in the step 1-8 into a non-sequence module of optical design software;
and 2-3, calculating a light path trace to obtain the axial position of the focal point of the cutting optical system under the irradiation of the laser power P.
2. The method for calculating the axial position of the focal point of the high-power laser cutting optical system according to claim 1, wherein the coordinate data of each facet point in the cutting optical system under the irradiation of the laser power P is obtained by using multi-physical-field coupling software.
3. The method for calculating the axial position of the focal point of the high-power laser cutting optical system according to claim 1, wherein the coordinate data of each surface point in the cutting optical system obtained by the calculation of the multi-physical-field coupling software is imported into a non-sequence module of optical design software for optical path calculation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910874987.2A CN110598332A (en) | 2019-09-19 | 2019-09-19 | Method for calculating axial position of focus of high-power laser cutting optical system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910874987.2A CN110598332A (en) | 2019-09-19 | 2019-09-19 | Method for calculating axial position of focus of high-power laser cutting optical system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110598332A true CN110598332A (en) | 2019-12-20 |
Family
ID=68860211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910874987.2A Pending CN110598332A (en) | 2019-09-19 | 2019-09-19 | Method for calculating axial position of focus of high-power laser cutting optical system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110598332A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111121621A (en) * | 2019-12-24 | 2020-05-08 | 北京理工大学 | Method for analyzing position error of main lens blocking mirror of large-aperture film-based diffraction optical system |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020090122A1 (en) * | 2000-11-03 | 2002-07-11 | Baer Thomas M. | Road map image guide for automated microdissection |
| US20050213022A1 (en) * | 2004-03-25 | 2005-09-29 | Yoshitaka Kawada | Method and apparatus for correcting a defective pixel of a liquid crystal display |
| CN105531071A (en) * | 2013-09-13 | 2016-04-27 | 通快机床两合公司 | Devices and methods for monitoring, more particularly controlling, a cutting process |
| CN106994557A (en) * | 2017-04-20 | 2017-08-01 | 武汉铱科赛科技有限公司 | A kind of dynamic controllable laser-processing system and method for focal position of laser |
| CN108319113A (en) * | 2018-01-31 | 2018-07-24 | 宁波大学 | The distortion correcting method of processing micro structure in a kind of capillary glass tube |
| CN108581243A (en) * | 2018-05-15 | 2018-09-28 | 大族激光科技产业集团股份有限公司 | Laser focal shift amount removing method |
| CN108875264A (en) * | 2018-07-06 | 2018-11-23 | 厦门大学 | A kind of method for building up of the laser source model for femtosecond laser ablation emulation |
| CN108875114A (en) * | 2017-09-07 | 2018-11-23 | 湖南大学 | A kind of focusing laser beam Characteristic parameter identification method |
| CN109416419A (en) * | 2016-04-25 | 2019-03-01 | 普雷茨特两合公司 | For the beam shaped optical system of laser cutting and the equipment with beam shaped optical system |
-
2019
- 2019-09-19 CN CN201910874987.2A patent/CN110598332A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020090122A1 (en) * | 2000-11-03 | 2002-07-11 | Baer Thomas M. | Road map image guide for automated microdissection |
| US20050213022A1 (en) * | 2004-03-25 | 2005-09-29 | Yoshitaka Kawada | Method and apparatus for correcting a defective pixel of a liquid crystal display |
| CN105531071A (en) * | 2013-09-13 | 2016-04-27 | 通快机床两合公司 | Devices and methods for monitoring, more particularly controlling, a cutting process |
| CN109416419A (en) * | 2016-04-25 | 2019-03-01 | 普雷茨特两合公司 | For the beam shaped optical system of laser cutting and the equipment with beam shaped optical system |
| CN106994557A (en) * | 2017-04-20 | 2017-08-01 | 武汉铱科赛科技有限公司 | A kind of dynamic controllable laser-processing system and method for focal position of laser |
| CN108875114A (en) * | 2017-09-07 | 2018-11-23 | 湖南大学 | A kind of focusing laser beam Characteristic parameter identification method |
| CN108319113A (en) * | 2018-01-31 | 2018-07-24 | 宁波大学 | The distortion correcting method of processing micro structure in a kind of capillary glass tube |
| CN108581243A (en) * | 2018-05-15 | 2018-09-28 | 大族激光科技产业集团股份有限公司 | Laser focal shift amount removing method |
| CN108875264A (en) * | 2018-07-06 | 2018-11-23 | 厦门大学 | A kind of method for building up of the laser source model for femtosecond laser ablation emulation |
Non-Patent Citations (4)
| Title |
|---|
| CHRISTOPHER BOUCHER: "Modeling Thermally Induced Focal Shift in High-Powered Laser Systems", 《COMSOL.COM/BLOGS/MODELING-THERMALLY-INDUCED-FOCAL-SHIFT-HIGH-POWERED-LASER-SYSTEMS》, 18 November 2014 (2014-11-18), pages 1 - 10 * |
| WILLIAM H. CARTER等: "Focal shift and concept of effective Fresnel number for a Gaussian laser beam", 《APPLIED OPTICS》, vol. 21, no. 11, 1 June 1982 (1982-06-01), pages 1989 - 1994 * |
| 佟玲等: "激光焦点控制磁力驱动的控制特性实验对比分析", 《国防科技大学学报》, vol. 40, no. 03, 28 June 2018 (2018-06-28), pages 120 - 126 * |
| 陈爔等: "3A21铝合金表面激光毛化坑点形貌演变规律", 《激光与光电子学进展》, vol. 56, no. 24, 24 June 2019 (2019-06-24), pages 173 - 180 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111121621A (en) * | 2019-12-24 | 2020-05-08 | 北京理工大学 | Method for analyzing position error of main lens blocking mirror of large-aperture film-based diffraction optical system |
| CN111121621B (en) * | 2019-12-24 | 2021-04-02 | 北京理工大学 | Method for analyzing position error of main lens blocking mirror of large-aperture film-based diffraction optical system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103801838B (en) | The wide laser galvanometer scanning fast etching method of a kind of modified line | |
| CN104191089B (en) | Three-Dimensional Dynamic focusing mark etching system and method based on Laser Output Beam | |
| EP2716397A1 (en) | Laser working head, laser working device, optical system for laser working device, laser working method, and laser focusing method | |
| CN103212786B (en) | Laser processing device and method thereof | |
| CN110977152A (en) | SLM double-laser combined machining system | |
| CN102658431B (en) | Device capable of automatically diagnosing and correcting divergence angle and beam quality of laser beam | |
| CN105140765A (en) | Resonant cavity module of fiber laser and fiber laser thereof | |
| CN110598332A (en) | Method for calculating axial position of focus of high-power laser cutting optical system | |
| CN204287551U (en) | The laser of a kind of macro-energy, narrow spaces and the coupling device of optical fiber | |
| JP2016514621A (en) | Apparatus for generating a laser beam having a linear intensity distribution | |
| CN103252575A (en) | Optical transmission method and system for laser material machining | |
| CN211564832U (en) | SLM double-laser composite processing device | |
| CN103885186A (en) | Astigmatism eliminating light beam shaping system based on prism pair and cylindrical mirror | |
| CN105710369A (en) | Device for manufacturing three-dimensional body layer by layer, control method and scanning method | |
| CN203607373U (en) | Laser annealing device | |
| CN103313817A (en) | Laser processing system | |
| CN104966986A (en) | Directive test system for assembling laser array | |
| CN203838413U (en) | Anastigmatic light beam shaping system based on prisms and cylindrical mirror | |
| CN218657326U (en) | Multi-focus laser ablation system | |
| CN105925792A (en) | Laser shock processing system | |
| CN103692090A (en) | Laser cutting device and method for LOCA (Liquid Optical Clear Adhesive) of touch panel | |
| CN114167549A (en) | Optical fiber laser beam combining device | |
| CN210427969U (en) | ZOOM cutting device based on adjustable annular light spot of aspherical mirror | |
| JP6808892B2 (en) | Combined wave optics | |
| CN111291458B (en) | Method for determining three-dimensional coordinates of focusing mirror surface profile of ECRH system antenna |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191220 |