CN101241042A - Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method - Google Patents

Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method Download PDF

Info

Publication number
CN101241042A
CN101241042A CNA200810030820XA CN200810030820A CN101241042A CN 101241042 A CN101241042 A CN 101241042A CN A200810030820X A CNA200810030820X A CN A200810030820XA CN 200810030820 A CN200810030820 A CN 200810030820A CN 101241042 A CN101241042 A CN 101241042A
Authority
CN
China
Prior art keywords
optical element
optical
foozle
beam quality
error
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.)
Granted
Application number
CNA200810030820XA
Other languages
Chinese (zh)
Other versions
CN101241042B (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN200810030820XA priority Critical patent/CN101241042B/en
Publication of CN101241042A publication Critical patent/CN101241042A/en
Application granted granted Critical
Publication of CN101241042B publication Critical patent/CN101241042B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Optical Devices Or Fibers (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

The invention discloses a quantitative analysis method of optical manufacturing error and beam quality relation of powerful laser system, comprising following steps: 1) adopting the optical detecting device to obtain manufacturing error of the detected optical element; 2) calculating error grads distribution according to the manufacturing error of the optical element; 3) calculating transfer function of the error grads to beam phase modulation; 4) calculating transmission influence of the optical system to the laser beam according to structural parameter of the optical system; 5) calculating beam quality index of target according to quantitative relation mode between the manufacturing error and beam quality of the optical element; 6) judging qualification of the optical element according to request of the powerful laser system to the beam quality for providing data analytical result to instruct modified process, which is significant to improve working performance of the powerful laser system, to design and control the manufacturing error of the optical element under confirmative processing condition.

Description

The quantitative analysis method of strong laser system optical foozle and beam quality relation
Technical field
The invention belongs to the optical assessment method, relate to the analytical approach of quantitative relationship between the optical element foozle that is used for strong laser system and the beam quality.
Background technology
In strong laser system, by optical transmitting system laser beam is outputed on the focal plane of far field, thereby reduce the angle of divergence, and laser beam is focused on, to realize high cumulative requirement by Focused Optical system.There is face shape error in optical element in manufacture process, according to the Small-scale Self-focusing theory, the optical element foozle will form the small sample perturbations to light field amplitude and phase place, after transmission, because non-linear growth, the disturbance of amplitude and phase place all will form the modulation to light intensity; Optical element intermediate frequency, high frequency error be as the noise source of beam intensity and phase perturbation, also is to cause one of main root that focal spot secondary lobe and Nonlinear Self-Focusing destroy.The parameter that is usually used in estimating beam quality also has M except the Si Telieer ratio 2The factor, far-field divergence angle and encircled power compare etc.; Compare with other index, encircled power is than more being applicable to energy transport, coupled mode application scenario.Quantitative relationship at optical element foozle and beam quality is analyzed, for the serviceability that improves strong laser system, the foozle of optical element is significant under the determinacy processing conditionss such as appropriate design and the digital control processing of control small abrasive nose, magnetorheological polishing, ion beam polishing.
The U.S. is in the process of development laser ICF engineering " national portfire (NIF) ", the LLNL laboratory has proposed to estimate the new method of optical element foozle: adopt error gradient GRMS to characterize the low-frequency range of face shape error, but it does not provide the quantitative formula with the beam quality index.Photoelectric Technology Inst., Chinese Academy of Sciences starts with from the frequency domain distribution of optical element foozle, utilize power spectral density function PSD, some diffraction function to analyze encircled power, with the energy loss situation in definite corresponding band error range, and then the qualification of evaluation optical element; But this method is based on frequency-domain analysis, and what provide is a kind of globality evaluation conclusion, can't determine the distribution character of optical element error gradient and far field light intensity, thereby can not provide guidance for next step carries out optics correction processing.
In strong laser system, what beam quality indexs such as encircled power ratio, light distribution were played main influence is the face shape error gradient of optical element.At present than between also do not have clear and definite, general quantitative analysis method about error gradient and encircled power both at home and abroad, this area research still belongs to blank.
We once proposed the influence relation (referring to be published in " optical precision engineering " paper of 2007 year 9th phase " optics face shape error to the influence of encircled power ratio ") of optical component surface shape error to far field encircled power ratio, its emphasis is to utilize Gaussian phase error simulated optical element mismachining tolerance, analyze its influence, and compare to determine the scope of application of encircled power than computing formula by theoretical analysis and simulation result to far field encircled power ratio.
The weak point of this evaluation method is:
One: only analyze at optical transmitting system, analysis result not necessarily is applicable to other optical system;
Its two: only Theoretical Calculation and simulation result are compared, and do not have concrete experimental result to verify;
Its three: be not applied to during concrete example estimates.
Summary of the invention
The technical problem to be solved in the present invention: overcome the deficiency of existing assessment technique and method, propose the quantitative analysis method of a kind of strong laser system optical foozle and beam quality relation.This method is for the serviceability that improves strong laser system, and the foozle of optical element is significant under appropriate design, the control determinacy processing conditions.
Realization the object of the invention technical scheme is finished by following steps:
1) adopt optical detection apparatus to obtain the foozle data of the measured optical unit;
2), calculate its error gradient and distribute according to the foozle data of optical element;
3) according to the optical element foozle laser beam is carried out the characteristics of phase modulation (PM), by formula:
H s ( u , v ) = exp [ - k 2 σ Δ 2 ( u 2 + v 2 ) ]
Error of calculation gradient is to the transport function of light beam phase modulation (PM), and in the formula, k=2 π/λ is a wave number, and λ is an optical maser wavelength, σ ΔBe the root-mean-square value of optical element foozle gradient, (u v) is the volume coordinate of optical element;
4) according to the structural parameters of optical system, by formula:
H 0 ( u , v ) = T ( u , v ) · π A 0 2 ω 0 2 2 · exp [ - L 12 ( u 2 + v 2 ) 8 L 0 ω 0 2 L 2 2 ]
The calculating optical system is to the transmission influence of laser beam, and in the formula, (u v) is a transmission factor to T, A 0, ω 0Be respectively the amplitude and the waist width of incident Gaussian beam, L 1, L 2Be transmission range, L 0=k 2ω 0 4+ 4L 1 2, L 12=4L 0L 2 2+ 8L 0L 1L 2+ L 0 2For different optical systems, and corresponding transmission factor T (u, v) different.
5) by the quantitative relationship model between optical element foozle and the beam quality:
PIB ( r ) = ∫ 0 r ρ · I ( ρ ) · dρ ∫ 0 + ∞ ρ · I ( ρ ) · dρ
Calculate the beam quality index on the target surface, described quantitative relationship model is according to the light distribution of laser beam on optical element surface:
Figure S200810030820XD00034
Make ρ 2=x 2+ y 22L 2 2(f x 2+ f y 2) after obtain;
6) according to the requirement of strong laser system to beam quality, judge the qualification of optical element, provide the data analysis result for instructing to revise to process.
The encircled power of correspondence is more qualified than being considered as greater than 84% when focal spot radius r=w (w is the laser beam beamwidth).When not meeting the demands, determine the bigger zone of optical element surface error gradient distribution amplitude, provide the data analysis result for instructing to revise to process.
The inventive method compared with prior art, advantage is:
1, the inventive method is based on the Gradient distribution of optical element foozle, directly set up the relational model between itself and the beam quality index, realize the numerical solution of beam quality indexs such as far-field intensity distribution, encircled power under the determinacy optics processing conditions compare, and unqualified in the specific light channel structure of this class of optical transmitting system;
2, method is verified analysis result by experiment;
3, can analyze the optical element foozle respectively to the transmission influence of the modulation effects of laser beam, optical system structure parameter by step 3), step 4) to laser beam.
The inventive method is for the serviceability that improves strong laser system, and the foozle of optical element is significant under appropriate design, the control determinacy processing conditions.
Description of drawings
Fig. 1 is the strong laser system light channel structure synoptic diagram of measuring optical element foozle to laser beam impact;
Fig. 2 is the optical element foozle distribution plan that utilizes the Wyko wavefront interferometer to measure;
Fig. 3 is the analogous diagram of optical element foozle to the far-field intensity distribution influence;
Fig. 4 is after incoming laser beam is influenced by the optical element foozle, the far-field intensity distribution figure that adopts the light intensity analyser to measure;
Fig. 5 is the influence curve comparison diagram of optical element foozle gradient to the encircled power ratio.
Embodiment
Embodiment:
By the inventive method, the face shape error of analyzing Φ 100mm devitrified glass level crossing after the I.B.M. under the laser beam effect to the quantitative effect of beam quality.
The applied strong laser system light channel structure of the inventive method comprises laser instrument 1, diaphragm 2, optical element 3, long focus lens 4 and light intensity analyser 5 to be tested as shown in Figure 1.Analyze the quantitative relationship of optics foozle and beam quality by the inventive method.
1. adopt the Wyko wavefront interferometer to detect the face shape error data of Φ 100mm devitrified glass level crossing after I.B.M., and utilize Wyko wavefront interferometer Survey Software to carry the influence that analytic function is removed constant term and inclination item, measurement result as shown in Figure 2.As can be seen from the figure the face shape error at irradiation area center (promptly being caused by temperature distortion) maximum illustrates and is heated at most, and is consistent with the Gaussian incident beam.
Also can adopt optical detection apparatuss such as Hartmann's array, shearing interferometer, also can be equipped with corresponding catoptron in case of necessity, remove the face shape error constant term in measurement or the analytic process and the influence of the item that tilts from standard apparatus or LASER Light Source etc.
2. according to the measurement data of optical element foozle, calculate its Gradient distribution.
By the Wyko wavefront interferometer obtain optical element the foozle data W (x, y) after, obtain the error gradient W of each measurement point X, Y two directions respectively X(x i, y j) and W Y(x i, y j).Because error gradient is very responsive to noise,, need to adopt high-precision gradient calculation formula in order to guarantee computational accuracy.We adopt Richard's extrapolation formula to find the solution:
W X ( x i , y j ) = 4 3 W ( x i + 1 , y j ) - W ( x i - 1 , y j ) 2 h x - 1 6 W ( x i + 2 , y j ) - W ( x i - 2 , y j ) 2 h x
W Y ( x i , y j ) = 4 3 W ( x i , y j + 1 ) - W ( x i , y j - 1 ) 2 h y - 1 6 W ( x i , y j + 2 ) - W ( x i , y j - 2 ) 2 h y
In the formula, (x y) is the foozle of optical element, h to W x, h yRepresent the sampling interval of margin of error strong point respectively in X, Y direction.Sampling interval is more little, and the error of calculation of following formula is more little.H in an embodiment x=h y=109.37 μ m, following formula can obtain sufficiently high computational accuracy.
Based on the aforementioned calculation result, the Gradient distribution of optical element foozle is:
Δ ( x , y ) = W X 2 ( x , y ) + W Y 2 ( x , y )
3. according to the optical element foozle laser beam is carried out the characteristics of phase modulation (PM), calculate its modulation transfer function.
Under determinacy optics processing conditionss such as the digital control processing of computer control small abrasive nose, magnetorheological polishing, ion beam polishing, optical element foozle gradient to the transport function that light beam carries out phase modulation (PM) is:
H s(u,v)=E{exp[i2kW(x-u,y-v)-i2kW(x,y)]}
The employing first approximation is expressed as:
H s ( u , v ) ≅ E { exp [ i 2 ku ∂ W ( x , y ) ∂ x + i 2 kv ∂ W ( x , y ) ∂ y ] }
According to Probability Statistics Theory, can get:
H s ( u , v ) = exp [ - 2 k 2 ( u 2 σ Δx 2 + v 2 σ Δy 2 + 2 ruv σ Δx σ Δy ) ]
Generally speaking, the root-mean-square value approximately equal and the correlativity of X, Y deflection error Gradient distribution are very little, and σ is arranged this moment Δ x≈ σ Δ y, r ≈ 0, i.e. σ Δ 2≈ 2 σ Δ x 2≈ 2 σ Δ y 2Thereby, have:
H s ( u , v ) = exp [ - k 2 σ Δ 2 ( u 2 + v 2 ) ]
In the formula, k=2 π/λ is a wave number, and λ is an optical maser wavelength, σ ΔBe the gradient root-mean-square value of optical element foozle, (u v) is the volume coordinate of optical element.
4. according to the structural parameters of optical system, the calculating optical system is to the transmission influence of laser beam.
In strong laser system, incident beam commonly used is a Gaussian beam.At strong laser system shown in Figure 1, owing to before face shape error is modulated laser beam, do not exist other optical element that laser beam is produced the transmission influence, so transmission factor T (u, v)=1, thereby optical system to the transmission influence of Gaussian beam is:
H 0 ( u , v ) = ∫ - ∞ + ∞ ∫ - ∞ + ∞ P D ( x , y ) P D * ( x - u , y - v ) exp [ ik 2 L 2 ( 2 xu + 2 yv - u 2 - v 2 ) ] dxdy
In the formula, P D(x y) is subjected to light intensity distribution of amplitudes behind the design of Optical System error effect for laser beam.
P D ( x , y ) = A 0 k ω 0 2 k ω 0 2 + i 2 L 1 exp ( ik L 1 - k ( x 2 + y 2 ) k ω 0 2 + i 2 L 1 )
In the formula, A 0, ω 0Be respectively the amplitude and the waist width of incident Gaussian beam, L 1, L 2Be the Laser Transmission distance.
With P D(x, y) expression formula substitution H 0(u, v) in and carry out Fourier integral and obtain:
H 0 ( u , v ) = π A 0 2 ω 0 2 2 · exp [ - L 12 ( u 2 + v 2 ) 8 L 0 ω 0 2 L 2 2 ]
In the formula, L 0=k 2ω 0 4+ 4L 1 2, L 12=4L 0L 2 2+ 8L 0L 1L 2+ L 0 2
5. calculate the beam quality index on the target surface.
After the transmission influence of laser beam process optical system and the modulation effects of optics foozle, the light distribution on target surface is:
Figure S200810030820XD00074
Light channel structure at Fig. 1 calculates:
I ( f x , f y ) = 4 π 2 A 0 2 L 0 ω 0 4 λ 2 ( L 12 + 8 k 2 L 0 ω 0 2 L 2 2 σ Δ 2 ) · exp { - 8 L 2 2 ω 0 2 L 0 π 2 ( f x 2 + f y 2 ) L 12 + 8 k 2 L 0 ω 0 2 L 2 2 σ Δ 2 }
In the formula, f x=x/ λ L 2, f y=y/ λ L 2
For circular bore optical element, make ρ 2=x 2+ y 22L 2 2(f x 2+ f y 2), then the encircled power in the radius r than PIB is:
PIB ( r ) = ∫ 0 r ρ · I ( ρ ) · dρ ∫ 0 + ∞ ρ · I ( ρ ) · dρ
The substitution of I (ρ) expression formula found the solution:
PIB ( r ) = 1 - exp { - 2 ω 0 2 L 0 k 2 r 2 L 12 + 8 k 2 L 0 ω 0 2 L 2 2 σ Δ 2 }
According to the above-mentioned relation model, can draw the quantitative relationship between the optical element foozle and beam quality index under the determinacy processing conditions easily.The far-field intensity distribution figure of Fig. 3 for coming out at the optical element foozle emulation that Fig. 2 surveyed finds that its distribution is approximated to Gaussian, and be closely similar with incident beam.
6. according to the requirement of strong laser system, judge the qualification of optical element, provide the data analysis result for instructing to revise to process to beam quality.
In strong laser system shown in Figure 1, the encircled power of correspondence is more qualified than being considered as greater than 84% when focal spot radius r=w (w is the laser beam beamwidth).When not meeting the demands, determine the bigger zone of error gradient distribution amplitude, provide the data analysis result for instructing to revise to process.
In the present embodiment, the encircled power of calculating ratio is 86.47%, greater than 84%, so meet the demands.
By the target practice experiment of indoor laser system, analyze the influence relation of optical element foozle gradient, experimental result such as Fig. 4, shown in Figure 5 to light distribution on the target surface, encircled power ratio.The incoming laser beam of measuring with PRIMES light intensity analyser is subjected to the far-field intensity distribution situation after the optical element foozle modulation effects as seen from Figure 4, and Fig. 5 represents to adopt resulting three encircled powers of model formation Theoretical Calculation, numerical simulation and experimental technique than curve comparison diagram respectively.We find three curve basically identicals, and this has verified the correctness of the inventive method.
According to the evaluation procedure that the inventive method is set up, can not practice shooting to test and just learn accurately that the optical element foozle waits the quantitative effect of beam quality index to encircled power ratio by the outfield; And then, judge the qualification of optical element according to the requirement of strong laser system to beam quality, provide the data analysis result for instructing to revise to process.

Claims (1)

1, the quantitative analysis method of a kind of strong laser system optical foozle and beam quality relation is characterized in that finishing by following steps:
1) adopt optical detection apparatus to obtain the foozle data of the measured optical unit;
2), calculate its error gradient and distribute according to the foozle data of optical element;
3) according to the optical element foozle laser beam is carried out the characteristics of phase modulation (PM), by formula
H s ( u , v ) = exp [ - k 2 σ Δ 2 ( u 2 + v 2 ) ]
Error of calculation gradient is to the transport function of light beam phase modulation (PM), and in the formula, k=2 π/λ is a wave number, and λ is an optical maser wavelength, σ ΔBe the root-mean-square value of optical element foozle gradient, (u v) is the volume coordinate of optical element;
4) according to the structural parameters of optical system, by formula:
H 0 ( u , v ) = T ( u , v ) · π A 0 2 ω 0 2 2 · exp [ - L 12 ( u 2 + v 2 ) 8 L 0 ω 0 2 L 2 2 ]
The calculating optical system is to the transmission influence of laser beam, and in the formula, (u v) is a transmission factor to T, A 0, ω 0Be respectively the amplitude and the waist width of incident Gaussian beam, L 1, L 2Be transmission range, L 0=k 2ω 0 4+ 4L 1 2, L 12=4L 0L 2 2+ 8L 0L 1L 2+ L 0 2For different optical systems, and corresponding transmission factor T (u, v) different;
5) by the quantitative relationship model between optical element foozle and the beam quality:
PIB ( r ) = ∫ 0 r ρ · I ( ρ ) · dρ ∫ 0 + ∞ ρ · I ( ρ ) · dρ
Calculate the beam quality index on the target surface, described quantitative relationship model is according to the light distribution of laser beam on optical element surface:
Figure S200810030820XC00014
Make ρ 2=x 2+ y 2After obtain;
6) according to the requirement of strong laser system to beam quality, judge the qualification of optical element, provide the data analysis result for instructing to revise to process.
CN200810030820XA 2008-03-14 2008-03-14 Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method Expired - Fee Related CN101241042B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810030820XA CN101241042B (en) 2008-03-14 2008-03-14 Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810030820XA CN101241042B (en) 2008-03-14 2008-03-14 Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method

Publications (2)

Publication Number Publication Date
CN101241042A true CN101241042A (en) 2008-08-13
CN101241042B CN101241042B (en) 2010-09-08

Family

ID=39932742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810030820XA Expired - Fee Related CN101241042B (en) 2008-03-14 2008-03-14 Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method

Country Status (1)

Country Link
CN (1) CN101241042B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107167301A (en) * 2017-07-11 2017-09-15 中国人民解放军国防科学技术大学 The method for evaluating laser beam quality Improvement
CN109060317A (en) * 2018-09-07 2018-12-21 西安工业大学 The characterisitic parameter pilot system and its course of work of long-distance propagation of laser beam

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107167301A (en) * 2017-07-11 2017-09-15 中国人民解放军国防科学技术大学 The method for evaluating laser beam quality Improvement
CN109060317A (en) * 2018-09-07 2018-12-21 西安工业大学 The characterisitic parameter pilot system and its course of work of long-distance propagation of laser beam

Also Published As

Publication number Publication date
CN101241042B (en) 2010-09-08

Similar Documents

Publication Publication Date Title
US9852518B2 (en) Method and system for calculating laser beam spot size
CN103543125B (en) All-optical gas detection method and device based on Michelson interference principle
CN105890878B (en) Measure the measurement apparatus and method of speculum damage threshold in real time using femtosecond laser
CN101980056B (en) Method and device for determining atmospheric turbulence parameter based on M<2> factor and light scintillation index
CN105137415A (en) Device and method for laser rangefinder receiving field-of-view calibration and optical axis parallelism measurement
CN104279978A (en) Three-dimensional figure detecting device and measuring method
CN105387933B (en) A kind of broadband Brewster window regulating device and method
CN101477047B (en) Nonlinear absorption measuring method based on lens geometric optical imaging
CN103712781A (en) Device and method for measuring multi-incidence-angle polarization interference in birefringence optical wedge optical axis direction
CN100535626C (en) Method for measuring focus and equivalent f coefficient using optical grating type wave-front curvature sensing unit
CN107607195B (en) A kind of beam quality measurement method obtained in real time based on complex amplitude
CN101241042B (en) Strong laser system optical manufacture error and light beam quality relationship quantitative analysis method
CN102680116B (en) Wave front aberration detection method and detection device
CN103884436A (en) Light beam phase on-line measuring device and method
CN103017664B (en) Method and system for calibrating laser beam analyzer
CN104634542A (en) Large-aperture optical element secondary exposure phase measuring device and measuring method
CN105403534A (en) Method for measuring transient state optical nonlinearity of material
CN102645408A (en) Phase object Z-scan-based pump-probe method
CN203163700U (en) Apparatus for measuring partially coherent Gaussian light beam wavefront phase radius
CN103542803B (en) Based on the synchronous phase shift interference device of Darman raster
CN102628713B (en) Curvature wavefront sensor based on digital micromirror device
CN106770335B (en) A kind of position phase defect detecting system and method based on reflection type point diffraction interferometer
CN104198053A (en) Wavefront detection method based on sub-wavelength grating array wavefront sensor
CN103063162B (en) A kind of method of measure portion coherent Gaussian beam Wave-front phase radius
CN209357082U (en) A kind of image correlativity measuring device based on spatial light modulator

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100908

Termination date: 20150314

EXPY Termination of patent right or utility model