CN102338680A - Method for measuring micro-impulse based on multi-beam laser heterodyne second harmonic method and torsion pendulum method - Google Patents

Abstract

The invention discloses a method for measuring micro-impulse based on a multi-beam laser heterodyne second harmonic method and a torsion pendulum method, relates to a micro-impulse measurement method, and solves the problem of low measurement accuracy caused by poor effect of acquiring a laser difference frequency signal and low operation speed of signal processing in the conventional method for measuring the micro-impulse by using the multi-beam laser heterodyne second harmonic method. The multi-beam laser heterodyne second harmonic method is applicable to the micro-impulse measurement method; the measurement of the micro-impulse is converted into the measurement of a pendulum angle of torsion pendulum, and the size of the micro-impulse is obtained by directly measuring the pendulum angle of the torsion pendulum by using the multi-beam laser heterodyne second harmonic method, so measurement accuracy is effectively improved; when a rotating angle is smaller than 5 DEG, the measured micro-impulse has a linear relation with an incident angle, and a measurement error is smaller than 0.5 percent; requirements for measurement of the impulse of a laser micro-thruster can be met; and good measurement means can be provided for evaluation of the performance of the laser micro-thruster.

Description

Measure the method for little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking

Technical field

The present invention relates to a kind of method of measuring little momentum.

Background technology

The laser micro thruster has extensive and deep application prospect at microsatellite attitude and track control field; It has than leaps high, big, the minimum momentum of momentum dynamic range is little, low in energy consumption, the energy coupling efficiency high and be easy to significant advantages such as realization, lightweight and Digital Control, has received Chinese scholars and has paid close attention to widely.And momentum is an important parameter of reflection laser micro thruster performance, and characteristics are that magnitude is little, are about 10 ^-7～10 ^-5Ns.People such as the Photonic Associates Phipps of group have proposed in 1999 with the small momentum that rocks the generation of systematic survey laser micro thruster, and carry out the test of microthrust device performance parameter with it; Domestic subsequently Chinese University of Science and Technology and equipment command technology institute have also carried out correlative study.From the present both at home and abroad result of study of report, on the one hand, the noise of measuring system can influence the precision of system, in little momentum magnitude, and systematic error even reached 50%; Simultaneously, in the power effect time, target plane departs from the focal plane, and the energy coupling efficiency reduces, and this also can influence the measurement of little momentum, and therefore conventional little momentum measuring system is difficult to satisfy measurement requirement.

Laser interferance method can effectively solve above two problems that the conventionally test system exists, and improves the measuring accuracy of system.The method that adopts two corner cubes to form variate replaces original light pointer method measurement to rock the angle of rotation, has improved the precision of system greatly; Rock Push Technology quality in 2010 and be increased to 58g, overcome the out of focus problem by original 0.2g.Result of study shows that the introducing of laser interferance method has greatly improved the performance of rocking test macro, can satisfy the test request of the small momentum of laser micro thruster.But because measuring amount is more indirectly, accidental error is bigger, so measuring accuracy can be very not high yet.

And in optical measuring method; Advantages such as the laser heterodyne measurement technology has that high room and time resolution, measuring speed are fast, precision is high, the linearity good, antijamming capability is strong, dynamic response is fast, good reproducibility and measurement range are big and enjoy Chinese scholars to pay close attention to; The laser heterodyne measurement technology has been inherited the plurality of advantages of heterodyne technology and Doppler technology, is one of present superhigh precision measuring method.This method has become one of significant technology of modern ultraprecise detection and surveying instrument, is widely used in ultra precise measurement, detection, process equipment, laser radar system etc.

But the existing method of the little momentum of multi-beam laser heterodyne measurement that adopts causes measuring accuracy lower owing to laser signal difference frequency signal collection effect arithmetic speed poor, signal Processing slowly.

Summary of the invention

The lower problem of measuring accuracy that the present invention causes owing to laser difference frequency signal collection effect arithmetic speed poor, signal Processing for the method that solves the little momentum of existing employing multi-beam laser heterodyne measurement slowly, thus a kind of method of measuring little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking is provided.

Measure the method for little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking, it is to adopt the system of measuring little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking to realize that said system comprises H ₀Solid state laser, rock system, quarter-wave plate, galvanometer, polarizing beam splitter mirror PBS, convergent lens, pulsed laser, plane standard mirror, photodetector and digital information processing system;

Wherein said H ₀Solid state laser, rock system, quarter-wave plate, galvanometer, polarizing beam splitter mirror PBS, convergent lens and plane standard mirror and be positioned at vacuum chamber, this vacuum chamber has a vacuum window, and the said system of rocking is made up of standard beam, plane mirror and working medium target; On the plane of the end of crossbeam of standard beam, be adhesive with plane mirror, be symmetrically fixed with the working medium target on the opposite side plane of this crossbeam relative with this plane mirror, the reflecting surface of said plane mirror is vertical with the swaying direction of the crossbeam of standard beam; This standard beam is in the equilibrium state of level, and the optical axis of the target surface of said working medium target and pulsed laser emitted laser bundle is perpendicular;

H ₀The solid state laser emission of lasering beam is to the front surface of polarizing beam splitter mirror PBS; After the quarter-wave plate transmission, be transmitted into the plane of incidence of galvanometer through the folded light beam of this polarizing beam splitter mirror PBS; After the quarter-wave plate transmission, be emitted to polarizing beam splitter mirror once more through the folded light beam after the vibration mirror reflected; After this polarizing beam splitter mirror transmission, be incident to the plane of incidence of pasting the plane mirror on the standard beam, the folded light beam of this plane mirror is with incident angle θ ₀Oblique incidence to plane standard mirror, the reflected light of this plane standard mirror focus on the photosurface of photodetector through another vacuum window of this vacuum chamber after the convergent lens transmission, and photodetector output electric signal is given digital information processing system; Said digital information processing system is used for obtaining the suffered little momentum of crossbeam of standard beam according to the signal that receives continuously;

The method of measuring little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking is realized by following steps:

At first, adopt pulsed laser to send pulse laser excitation working medium target, make this working medium target produce plasma jet, the effect of regurgitating of the plasma spraying that is produced is rotated the crossbeam of standard beam;

Simultaneously, open H ₀The driving power of solid state laser and galvanometer, galvanometer is done uniformly accelerated motion under the driving power effect, and the optical frequency that difference incides the galvanometer front surface is constantly modulated;

Then, the digital information processing system signal that the continuous acquisition photodetector sends in rocking system's swing process, and all signals of continuous acquisition are handled, obtain the suffered little momentum of crossbeam of standard beam;

The suffered little momentum of the crossbeam of said standard beam is according to the crossbeam angle of oscillation θ ' acquisition of standard beam:

The relational expression of momentum and rotational angle is:

$I=\frac{2\mathrm{J\ω}}{D}\·{\mathrm{\θ}}^{\′}=\frac{4\mathrm{\πJ}}{\mathrm{DT}}\·{\mathrm{\θ}}^{\′}$

In the formula, k=4 π J/DT, wherein, J is the moment of inertia of the system of rocking, and T rocks the damping period of system for this, and D is a crossbeam length, and θ ' is the pivot angle of standard beam;

Make k=4 π J/DT, then said little momentum is:

I＝k·θ′

The pivot angle θ ' of said standard beam (8) is:

${\mathrm{\θ}}^{\′}=\frac{{\mathrm{\θ}}_0}{2}$

The pivot angle θ ' of said standard beam is the signal of the photodetector that obtains constantly according to difference, obtains through the multi-beam laser heterodyne method of quadratic harmonics, and detailed process is:

Incident field according to plane standard mirror:

E(t)＝E ₀exp(iω ₀t)

And the vibration equation of Doppler's galvanometer:

x(t)＝a(t ²/2)

Rate equation with Doppler's galvanometer:

v(t)＝at

Acquisition becomes through the catoptrical frequency of Doppler's galvanometer:

ω＝ω ₀(1+at/c)

ω in the formula ₀Be the laser angular frequency, a is a vibration acceleration, E ₀Be constant, c is the light velocity;

Then t-l/c arrives the reflection light field of plane standard mirror front surface constantly and is:

$E_1\left(t\right)={\mathrm{\α}}_1{E}_0\exp \left\{i\right[{\mathrm{\ω}}_0(1+\frac{a(t-l/c)}{c})t+{\mathrm{\ω}}_0\frac{\frac{a{(t-l/c)}^2}{2}}{c}\left]\right\}$

Wherein, l is the light path of galvanometer front surface to plane standard mirror front surface;

And through the light of the front surface transmission of plane standard mirror in difference constantly by plane standard mirror after surperficial m-1 secondary reflection, and then obtain m-1 bundle transmitted light, the light field of said m-1 bundle transmitted light is respectively:

$E_2\left(t\right)={\mathrm{\α}}_2{E}_0\exp \left\{i\right[{\mathrm{\ω}}_0(1+a\frac{t-\frac{1}{c}-\frac{2\mathrm{nd}\cos \mathrm{\θ}}{c}}{c})t+{\mathrm{\ω}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2\mathrm{nd}\cos \mathrm{\θ}}{c})}^2}{c}+2\mathrm{nd}\cos \mathrm{\θ})}{c}\left]\right\}$

$\begin{array}{c}.\\ .\\ .\end{array}$

$\begin{array}{c}.\\ .\\ .\end{array}$

$E_m\left(t\right)={\mathrm{\α}}_m{E}_0\exp \left\{i\right[{\mathrm{\ω}}_0(1+a\frac{t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\θ}}{c}}{c})]t$

$+{\mathrm{\ω}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\θ}}{c})}^2}{2}+2(m-1)\mathrm{nd}\cos \mathrm{\θ})}{c}\left]\right\}$

Wherein, m representes the number of folded light beam, and n is the refractive index of plane standard mirror, α ₁=r, α ₂=β β ' r ' ..., α _m=β β ' r ' ^(2m-3)R is light is injected plane standard mirror front surface from surrounding medium a reflectivity; To be light inject light from surrounding medium to β injects the transmissivity of plane standard mirror front surface from surrounding medium, and r ' is the reflectivity for the back surface of plane standard mirror, the transmissivity when β ' be the front surface of reflected light ejaculation plane standard mirror on back surface of plane standard mirror; D is the thickness of plane standard mirror, and θ is the refraction angle of light when inciding plane standard mirror front surface;

Total light field that photodetector receives is:

E(t)＝E ₁(t)+E ₂(t)+…+E _m(t)

Then the photocurrent of photodetector output is expressed as:

$I^{\&prime;}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{1}{Z}\underset{\mathrm{<figure-callout\; id="1"\; label="S"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">S</figure-callout>}}{\&Integral;\&Integral;}\frac{1}{2}[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]{[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]}^{*}\mathrm{ds}$

Wherein, e is an electron charge, and Z is the intrinsic impedance of detector surface medium, and η is a quantum efficiency, and S is the area of detector photosurface, and h is a Planck's constant, and v is a laser frequency;

Above-mentioned photocurrent obtains second harmonic signal behind low-pass filter (17) electric current of intermediate frequency is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;<figure-callout\; id="2"\; label="e"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">e</figure-callout>}}{2\mathrm{hv}}\frac{1}{Z}\underset{s}{\&Integral;\&Integral;}\underset{p=1}{\stackrel{\&infin;}{\mathrm{\&Sigma;}}}\underset{j=p+2}{\overset{\&infin;}{\mathrm{\&Sigma;}}}(E_p\left(t\right){E_j}^{*}\left(t\right)+{E_p}^{*}\left(t\right)E_j\left(t\right))\mathrm{ds}$

With the formula substitution following formula of all light fields, the arrangement result is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p{E_0}^2\cos (\frac{8{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c}-\frac{4l{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^3}-\frac{8p{\mathrm{\&omega;}}_0{\mathrm{an}}^2{\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">d</figure-callout>}}^2{\mathrm{<figure-callout\; id="2"\; label="cos"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">cos</figure-callout>}}^2\mathrm{\&theta;}}{c^3})$

Ignore l/c ³Event after can be reduced to:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}{E_0}^2\cos (\frac{8a{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c})\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p$

In the formula, p and j are natural number;

And then the frequency that obtains interference signal is:

f＝8andcosθω ₀/(2πc ²)＝4andω ₀/(πc ²)＝Kcosθ

Can know according to refraction law:

$\cos \mathrm{\&theta;}=\cos \left[\arcsin \left(\frac{\sin {\mathrm{\&theta;}}_0}{n}\right)\right]$

Then scale-up factor is:

K＝4andω ₀/(πc ²)

With the light source angle frequencies omega ₀, the thickness d of standard flat mirror, refractive index n, a is relevant for the galvanometer acceleration.Therefore be easy to obtain the relational expression of momentum I and refraction angle θ, said momentum I is:

$I=\frac{k\&CenterDot;\arcsin \left\{n\sin \right[\arccos \left(\cos \mathrm{\&theta;}\right)\left]\right\}}{2}.$

Beneficial effect: the present invention adopts the multi-beam laser heterodyne method of quadratic harmonics to be applied in little momentum measuring method, and laser difference frequency signal collection effect is better, and the arithmetic speed of signal Processing is very fast, and the precision of measurement is higher.

Description of drawings

Fig. 1 of the present inventionly measures the structural representation of the system of little momentum based on the multi-beam laser heterodyne second harmonic of the method for rocking, and mark 5 is the skylight; Fig. 2 is a multi-beam laser principle of interference synoptic diagram.The different laser incident angle of Fig. 3 θ ₀Under the situation, multi-beam laser heterodyne second harmonic is measured the corresponding multi-beam laser heterodyne second harmonic signal Fourier transform frequency spectrum of minute angle; A is 5.976mrad, and B is 6.723mrad, and C is 7.470, and mrad, D are 8.217mrad, and E is 8.964mrad, and F is 9.711mrad, and G is 10.458mrad, and H is 11.205mrad.

Embodiment

Embodiment one, combination Fig. 1 and Fig. 2 explain this embodiment; Measure the method for little momentum based on the multi-beam laser heterodyne second harmonic of the method for rocking; It is that the system that employing is measured little momentum based on the multi-beam laser heterodyne second harmonic of the method for rocking realizes that said system comprises H ₀ Solid state laser 10, rock system, quarter-wave plate 12, galvanometer 13, polarizing beam splitter mirror PBS11, convergent lens 15, pulsed laser 6, plane standard mirror 14, photodetector 2 and digital information processing system 1;

Wherein said H ₀ Solid state laser 10, rock system, quarter-wave plate 12, galvanometer 13, polarizing beam splitter mirror PBS11, convergent lens 15 and plane standard mirror 14 and be positioned at vacuum chamber 4; This vacuum chamber 4 has a vacuum window 3, and the said system of rocking is made up of standard beam 8, plane mirror 9 and working medium target 7; On the plane of the end of crossbeam of standard beam 8, be adhesive with plane mirror 9; Be symmetrically fixed with working medium target 7 on the opposite side plane of this crossbeam relative with this plane mirror 9, the reflecting surface of said plane mirror 9 is vertical with the swaying direction of the crossbeam of standard beam 8; This standard beam 8 is in the equilibrium state of level, and the optical axis of the target surface of said working medium target 7 and pulsed laser emitted laser bundle is perpendicular;

H ₀The solid state laser emission of lasering beam is to the front surface of polarizing beam splitter mirror PBS11; After quarter-wave plate 12 transmissions, be transmitted into the plane of incidence of galvanometer 13 through the folded light beam of this polarizing beam splitter mirror PBS11; Folded light beam after galvanometer 13 reflections is emitted to polarizing beam splitter mirror once more after quarter-wave plate 12 transmissions; After this polarizing beam splitter mirror transmission, be incident to the plane of incidence of the plane mirror of pasting on standard beam 89; The folded light beam of this plane mirror 9 with incident angle θ 0 oblique incidence to the plane standard mirror 14; The reflected light of this plane standard mirror 14 focuses on the photosurface of photodetector 2 through another vacuum window 3 of this vacuum chamber 4 after convergent lens 15 transmissions, and photodetector 2 output electric signal are given digital information processing system 1; Said digital information processing system 1 is used for obtaining the suffered little momentum of crossbeam of standard beam 8 according to the signal that receives continuously;

The method of measuring little momentum based on the multi-beam laser heterodyne second harmonic of the method for rocking is:

At first, adopt pulsed laser 6 to send pulse laser excitation working medium target 7, make this working medium target 7 produce plasma jet, the effect of regurgitating of the plasma spraying that is produced is rotated the crossbeam of standard beam 8;

Simultaneously, open H ₀The driving power of solid state laser 10 and galvanometer 13, galvanometer 13 is done uniformly accelerated motion under the driving power effect, and the optical frequency that difference incides galvanometer 13 front surfaces is constantly modulated;

Then, digital information processing system 1 signal that continuous acquisition photodetector 2 sends in rocking system's swing process, and all signals of continuous acquisition are handled, obtain the suffered little momentum of crossbeam of standard beam 8;

The suffered little momentum of the crossbeam of said standard beam 8 is according to the crossbeam angle of oscillation θ ' acquisition of standard beam 8:

The relational expression of momentum and rotational angle is:

$I=\frac{2\mathrm{J\&omega;}}{D}\&CenterDot;{\mathrm{\&theta;}}^{\&prime;}=\frac{4\mathrm{\&pi;J}}{\mathrm{DT}}\&CenterDot;{\mathrm{\&theta;}}^{\&prime;}---\left(1\right)$

In the formula, k=4 π J/DT, wherein, J is the moment of inertia of the system of rocking, and T rocks the damping period of system for this, and D is a crossbeam length, and θ ' is the pivot angle of standard beam 8;

Then said little momentum is:

I＝k·θ′ (2)

The pivot angle θ ' of said standard beam 8 is:

${\mathrm{\&theta;}}^{\&prime;}=\frac{{\mathrm{\&theta;}}_0}{2}$

θ ₀It is the incident angle of plane standard mirror 14;

Can find out from formula (2), under the small angle approximation condition,, just can obtain the size of little momentum as long as known the pivot angle of standard beam.

The pivot angle θ ' of said standard beam 8 is signals of the photodetector that obtains constantly according to difference, obtains through the multi-beam laser heterodyne method of quadratic harmonics, and detailed process is:

As shown in Figure 2; Because light beam can constantly reflect and reflect between the front and rear surfaces of plane standard mirror 14; And this reflection and refraction for reflected light and transmitted light at infinity or the interference on the lens focal plane contribution is all arranged; So when interference is discussed, must consider repeatedly reflection and refraction effect, multi-beam laser promptly should be discussed interfere.

But; Because the optical mixing that transmit plane standard mirror front surface of laser after the reflected light of plane standard mirror front surface and back surface reflection k time and k+1 time; The amplitude of two difference frequency signals that produce differs 2～3 one magnitude; Through after the Fourier transform, in order to collect laser difference frequency signal and the arithmetic speed that improves signal Processing preferably, so the E of our the k secondary reflection only considering to be detected here _kE behind the surperficial k+2 secondary reflection of light and back _K+2The humorous frequency difference of the secondary that optical mixing produced.

When laser with incident angle θ ₀During oblique incidence, the incident field of plane standard mirror 14:

E(t)＝E ₀exp(iω ₀t) (3)

And the vibration equation of Doppler's galvanometer 13:

x(t)＝a(t ²/2) (4)

Rate equation with Doppler's galvanometer 13:

v(t)＝at (5)

Acquisition becomes through the catoptrical frequency of Doppler's galvanometer 13:

ω＝ω ₀(1+at/c) (6)

ω in the formula ₀Be the laser angular frequency, a is a vibration acceleration, E ₀Be constant, c is the light velocity;

Then t-l/c arrives the reflection light field of plane standard mirror front surface constantly and is:

$E_1\left(t\right)={\mathrm{\&alpha;}}_1{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+\frac{a(t-l/c)}{c})t+{\mathrm{\&omega;}}_0\frac{\frac{a{(t-l/c)}^2}{2}}{c}\left]\right\}---\left(7\right)$

Wherein l is the light path of galvanometer front surface to plane standard mirror front surface;

And through the light of the front surface transmission of plane standard mirror 14 in difference constantly by plane standard mirror (14) after surperficial m-1 secondary reflection, and then obtain m-1 bundle transmitted light, the light field of said m-1 bundle transmitted light is respectively:

$E_2\left(t\right)={\mathrm{\&alpha;}}_2{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+a\frac{t-\frac{1}{c}-\frac{2\mathrm{nd}\cos \mathrm{\&theta;}}{c}}{c})t+{\mathrm{\&omega;}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2\mathrm{nd}\cos \mathrm{\&theta;}}{c})}^2}{c}+2\mathrm{nd}\cos \mathrm{\&theta;})}{c}\left]\right\}$

$\begin{array}{c}.\\ .\\ .\end{array}---\left(8\right)$

$\begin{array}{c}.\\ .\\ .\end{array}$

$E_m\left(t\right)={\mathrm{\&alpha;}}_m{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+a\frac{t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\&theta;}}{c}}{c})]t$

$+{\mathrm{\&omega;}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\&theta;}}{c})}^2}{2}+2(m-1)\mathrm{nd}\cos \mathrm{\&theta;})}{c}\left]\right\}$

Wherein, m representes the number of folded light beam, and n is the refractive index of plane standard mirror, α ₁=r, α ₂=β β ' r ' ..., α _m=β β ' r ' ^(2m-3)R is light is injected plane standard mirror 14 front surfaces from surrounding medium a reflectivity; To be light inject light from surrounding medium to β injects the transmissivity of plane standard mirror 14 front surfaces from surrounding medium, and r ' is the reflectivity for the back surface (plane mirror) of plane standard mirror 14, the transmissivity when β ' be the front surface of reflected light ejaculation plane standard mirror 14 on back surface of plane standard mirror 14; D is the thickness of plane standard mirror 14, and θ is the refraction angle of light when inciding plane standard mirror 14 front surfaces;

Total light field that photodetector receives is:

E(t)＝E ₁(t)+E ₂(t)+…+E _m(t) (9)

Then the photocurrent of photodetector output can be expressed as:

$I^{\&prime;}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{1}{Z}\underset{\mathrm{<figure-callout\; id="1"\; label="S"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">S</figure-callout>}}{\&Integral;\&Integral;}\frac{1}{2}[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]{[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]}^{*}---\left(10\right)$

Wherein, e is an electron charge, and Z is the intrinsic impedance of detector surface medium, and η is a quantum efficiency, and S is the area of detector photosurface, and h is a Planck's constant, and v is a laser frequency;

Only consider E this moment _kAnd E _K+2The second harmonic difference frequency signal that optical mixing produced; Dc terms is through can filtering behind the low-pass filter, therefore, only considers here to exchange; This exchanges item and is commonly referred to electric current of intermediate frequency, and the electric current of intermediate frequency that behind low-pass filter 17, obtains second harmonic signal through photocurrent is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{2\mathrm{hv}}\frac{1}{Z}\underset{s}{\&Integral;\&Integral;}\underset{p=1}{\stackrel{\&infin;}{\mathrm{\&Sigma;}}}\underset{j=p+2}{\overset{\&infin;}{\mathrm{\&Sigma;}}}(E_p\left(t\right){E_j}^{*}\left(t\right)+{E_p}^{*}\left(t\right)E_j\left(t\right))\mathrm{ds}---\left(11\right)$

With the formula substitution following formula of all light fields, the arrangement result is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p{E_0}^2\cos (\frac{8{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c}-\frac{4l{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^3}-\frac{8p{\mathrm{\&omega;}}_0{\mathrm{an}}^2{\mathrm{<figure-callout\; id="2"\; label="d"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">d</figure-callout>}}^2{\mathrm{<figure-callout\; id="2"\; label="cos"\; filenames="HDA0000065395410000011.png"\; state="\{\{state\}\}">cos</figure-callout>}}^2\mathrm{\&theta;}}{c^3})---\left(12\right)$

Ignore l/c ³Event after can be reduced to:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}{E_0}^2\cos (\frac{8a{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c})\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p---\left(13\right)$

In the formula, p and j are natural number;

Can see that through (13) formula the information of standard mirror refraction angle, plane θ is all arranged in intermediate frequency item difference on the frequency that multiple beam heterodyne second harmonic mensuration obtains and the phase differential.Be primarily aimed at intermediate frequency item intermediate frequency rate variance and analyze, because adopt Fourier transform to be easy to realize frequency measurement.At this moment, according to (9) formula, can be designated as the frequency of interference signal:

f＝8andcosθω ₀/(2πc ²)＝4andω ₀/(πc ²)＝Kcosθ (14)

Can know according to refraction law:

$\cos \mathrm{\&theta;}=\cos \left[\arcsin \left(\frac{\sin {\mathrm{\&theta;}}_0}{n}\right)\right]---\left(15\right)$

Then scale-up factor is:

K＝4andω ₀/(πc ²) (16)

Therefore be easy to obtain the relational expression of momentum I and refraction angle θ, said momentum I is:

$I=\frac{k\&CenterDot;\arcsin \left\{n\sin \right[\arccos \left(\cos \mathrm{\&theta;}\right)\left]\right\}}{2}---\left(17\right)$

Rock the suffered little momentum of crossbeam of pivot angle and then realization measurement standard beam 8 through direct measurement.

The difference that embodiment two, this embodiment and embodiment one described multi-beam laser heterodyne second harmonic based on the method for rocking are measured the method for little momentum is; Multi-beam laser heterodyne second harmonic based on the method for rocking is measured in the system of little momentum; Digital information processing system 1 is made up of wave filter 17, prime amplifier 16, modulus converter A/D and digital signal processing controller DSP; Wave filter 17 carries out the signal of photodetector 2 outputs to send to prime amplifier 16 after the filtering; This prime amplifier 16 sends to modulus converter A/D after the signal that receives is amplified; Modulus converter A/D becomes digital signal to send to digital signal processing controller DSP the analog signal conversion that receives; Be solidified with fft algorithm among this digital signal processing controller DSP, digital signal processing controller DSP is used for the signal that receives is continuously handled, and obtains the suffered little momentum of crossbeam of standard beam 8 after the demodulation.

Below through concrete emulation experiment, verify effect of the present invention: measure the variation of adding standard beam front and back system cycle, calibrate the moment of inertia of system, the experimental result of demarcation is as shown in table 1, can obtain the size of k value according to calibration result.

Table 1 rocks the calibration result of parameter

In order to verify the feasibility of this method; Based on system shown in Figure 1, the pressure in the vacuum chamber is adjusted to 10.0Pa, utilize MATLAB to simulate working medium and be PVC (PVC)+2%C; Thickness is 180 μ m; The initial current of the pulse laser that pulsed laser sends is 5A, and the pulsewidth of said pulse laser is 50ms, and pulse laser and working medium interact and produces little momentum.

Get H _oSolid state laser wavelength X 2050nm; Standard beam 8 long D=15cm; Refractive index n=1.493983 of plane standard mirror 14, its thickness d are 2cm; The photosurface aperture of detector is R=1mm.The sensitivity of said photosurface is 1A/W.Get Doppler's galvanometer acceleration a=2 * 10 ³M/s ²

The different laser incident angle θ that utilized MATLAB emulation ₀Under the situation; It is as shown in Figure 3 that multi-beam laser heterodyne second harmonic is measured the corresponding multi-beam laser heterodyne signal Fourier transform frequency spectrum of minute angle; As can be seen from Figure 3, along with the increase of incident angle, the relative position of frequency spectrum moves promptly along with incident angle θ to the low frequency direction ₀The increase frequency reduce.Reason is: under the constant situation of scale-up factor K, because frequency f and incident angle θ ₀Relation is:

f＝Kcosθ＝Kcos[arcsin(sinθ ₀/n)] (18)

Incident angle θ ₀With the frequency relation of being inversely proportional to, as incident angle θ ₀Cos θ reduces thereupon during increase.Therefore, along with incident angle θ ₀The relative position of increase frequency spectrum move to the low frequency direction, Fig. 3 has verified the correctness of front theoretical analysis well.Also verified the feasibility of technical scheme according to the invention.

Need to prove that because heterodyne detection is a kind of detection mode of nearly diffraction limit, detection sensitivity is high, so the signal to noise ratio (S/N ratio) of the heterodyne second harmonic signal of Fig. 3 is very high.

Utilize above-mentioned multi-beam laser heterodyne second harmonic mensuration, continuous analog eight groups of data, obtained the simulation result of the little momentum of testing sample under the different incidence angles situation, as shown in table 2.

Table 2 different incidence angles feelings θ ₀Under the condition, the actual value I of little momentum _ActualWith simulation value I _i

Time	1	2	3	4	5	6	7	8
									Actual incident angle θ 0(mrad)	5.976	6.723	7.470	8.217	8.964	9.711	10.458	11.205
I Actual(×10 -6N·s)	22.183	24.956	27.729	30.502	33.275	36.048	38.820	41.593
									I i(×10 -6N·s)	22.299	24.832	27.867	30.392	33.415	35.930	38.943	41.449

Need to prove: utilize the emulation experiment data of table 2, can calculate the averaging analog value of little momentum according to (2) formula, the maximum relative error that finally can obtain the analogue value is less than 0.5%, and the precision that can find out this method is very high.Simultaneously; The analysis data it can also be seen that; Under the situation of small angle approximation, systematic error that environment brings and reading error are negligible in emulation, and the error in the emulation experiment mainly comes from trueness error and the round-off error in the computation process after the Fast Fourier Transform (FFT) (FFT).

Utilize method of the present invention to carry out the numerical simulation of pulse laser and PVC working medium target coupling little momentum that produces in the above-mentioned emulation experiment.The result shows that this method range of linearity is big and resolution is high, and the advantage of this angle-measuring method is to rotation sensitive; Insensitive to translation; Therefore test macro also has than strong anti-interference ability vibration, and particularly low-frequency vibration can return to working state of system within several seconds; Not only reduced measuring error, also reduced the requirement of measuring equipment with bad border of experiment.Simultaneously, when rotational angle less (less than 5 °), momentum and the incident angle surveyed are linear, and measuring error can satisfy laser micro thruster momentum Testing requirement less than 0.5%, for the performance of assessing the laser micro thruster provides good measurement means.

Claims (2)

1. measure the method for little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking, it is that the system that employing is measured little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking realizes that said system comprises H ₀Solid state laser (10), rock system, quarter-wave plate (12), galvanometer (13), polarizing beam splitter mirror PBS (11), convergent lens (15), pulsed laser (6), plane standard mirror (14), photodetector (2) and digital information processing system (1);

Wherein said H ₀Solid state laser (10), rock system, quarter-wave plate (12), galvanometer (13), polarizing beam splitter mirror PBS (11), convergent lens (15) and plane standard mirror (14) and be positioned at vacuum chamber (4); This vacuum chamber (4) has a vacuum window (3), and the said system of rocking is made up of standard beam (8), plane mirror (9) and working medium target (7); On the plane of the end of crossbeam of standard beam (8), be adhesive with plane mirror (9); Be symmetrically fixed with working medium target (7) on the opposite side plane of this crossbeam relative with this plane mirror (9), the reflecting surface of said plane mirror (9) is vertical with the swaying direction of the crossbeam of standard beam (8); This standard beam (8) is in the equilibrium state of level, and the optical axis of the target surface of said working medium target (7) and pulsed laser emitted laser bundle is perpendicular;

H ₀The solid state laser emission of lasering beam is to the front surface of polarizing beam splitter mirror PBS (11); After quarter-wave plate (12) transmission, be transmitted into the plane of incidence of galvanometer (13) through the folded light beam of this polarizing beam splitter mirror PBS (11); Folded light beam after galvanometer (13) reflection is emitted to polarizing beam splitter mirror once more after quarter-wave plate (12) transmission; After this polarizing beam splitter mirror transmission, be incident to the plane of incidence of the plane mirror of pasting on standard beam (8) (9), the folded light beam of this plane mirror (9) is with incident angle θ ₀Oblique incidence to plane standard mirror (14); The reflected light of this plane standard mirror (14) is after convergent lens (15) transmission; Another vacuum window (3) through this vacuum chamber (4) focuses on the photosurface of photodetector (2), and photodetector (2) output electric signal is given digital information processing system (1); Said digital information processing system (1) is used for obtaining the suffered little momentum of crossbeam of standard beam (8) according to the signal that receives continuously;

It is characterized in that: the method for measuring little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking is realized by following steps:

At first, adopt pulsed laser (6) to send pulse laser excitation working medium target (7), make this working medium target (7) produce plasma jet, the effect of regurgitating of the plasma spraying that is produced is rotated the crossbeam of standard beam (8);

Simultaneously, open H ₀The driving power of solid state laser (10) and galvanometer (13), galvanometer (13) is done uniformly accelerated motion under the driving power effect, and the optical frequency that difference incides galvanometer (13) front surface is constantly modulated;

Then, digital information processing system (1) signal that continuous acquisition photodetector (2) sends in rocking system's swing process, and all signals of continuous acquisition are handled, obtain the suffered little momentum I of crossbeam of standard beam (8);

The suffered little momentum of the crossbeam of said standard beam (8) is according to the crossbeam angle of oscillation θ ' acquisition of standard beam (8):

The relational expression of momentum and rotational angle is:

$I=\frac{2\mathrm{J\&omega;}}{D}\&CenterDot;{\mathrm{\&theta;}}^{\&prime;}=\frac{4\mathrm{\&pi;J}}{\mathrm{DT}}\&CenterDot;{\mathrm{\&theta;}}^{\&prime;}$

In the formula, k=4 π J/DT, wherein, J is the moment of inertia of the system of rocking, and T rocks the damping period of system for this, and D is a crossbeam length, and θ ' is the pivot angle of standard beam (8);

Then said little momentum is:

I＝k·θ′

The pivot angle θ ' of said standard beam (8) is:

${\mathrm{\&theta;}}^{\&prime;}=\frac{{\mathrm{\&theta;}}_0}{2}$

The pivot angle θ ' of said standard beam (8) is the signal of the photodetector that obtains constantly according to difference, obtains through the multi-beam laser heterodyne method of quadratic harmonics, and detailed process is:

Incident field according to plane standard mirror (14):

E(t)＝E ₀exp(iω ₀t)

And the vibration equation of Doppler's galvanometer (13):

x(t)＝a(t ²/2)

Rate equation with Doppler's galvanometer (13):

v(t)＝at

Acquisition becomes through the catoptrical frequency of Doppler's galvanometer (13):

ω＝ω ₀(1+at/c)

ω in the formula ₀Be the laser angular frequency, a is a vibration acceleration, E ₀Be constant, c is the light velocity;

Then t-l/c arrives the reflection light field of plane standard mirror (14) front surface constantly and is:

$E_1\left(t\right)={\mathrm{\&alpha;}}_1{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+\frac{a(t-l/c)}{c})t+{\mathrm{\&omega;}}_0\frac{\frac{a{(t-l/c)}^2}{2}}{c}\left]\right\}$

Wherein l is the light path of galvanometer front surface to plane standard mirror front surface;

And through the light of the front surface transmission of plane standard mirror (14) in difference constantly by plane standard mirror (14) after surperficial m-1 secondary reflection, and then obtain m-1 bundle transmitted light, the light field of said m-1 bundle transmitted light is respectively:

$E_2\left(t\right)={\mathrm{\&alpha;}}_2{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+a\frac{t-\frac{1}{c}-\frac{2\mathrm{nd}\cos \mathrm{\&theta;}}{c}}{c})t+{\mathrm{\&omega;}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2\mathrm{nd}\cos \mathrm{\&theta;}}{c})}^2}{c}+2\mathrm{nd}\cos \mathrm{\&theta;})}{c}\left]\right\}$

$\begin{array}{c}.\\ .\\ .\end{array}$

$\begin{array}{c}.\\ .\\ .\end{array}$

$E_m\left(t\right)={\mathrm{\&alpha;}}_m{E}_0\exp \left\{i\right[{\mathrm{\&omega;}}_0(1+a\frac{t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\&theta;}}{c}}{c})]t$

$+{\mathrm{\&omega;}}_0\frac{(a\frac{{(t-\frac{l}{c}-\frac{2(m-1)\mathrm{nd}\cos \mathrm{\&theta;}}{c})}^2}{2}+2(m-1)\mathrm{nd}\cos \mathrm{\&theta;})}{c}\left]\right\}$

Wherein, m representes the number of folded light beam, and n is the refractive index of plane standard mirror, α ₁=r, α ₂=β β ' r ' ..., α _m=β β ' r ' ^(2m-3)R is light is injected plane standard mirror (14) front surface from surrounding medium a reflectivity; To be light inject light from surrounding medium to β injects the transmissivity of plane standard mirror (14) front surface from surrounding medium, and r ' is the reflectivity for the back surface of plane standard mirror (14), the transmissivity when β ' be the front surface of reflected light ejaculation plane standard mirror (14) on back surface of plane standard mirror (14); D is the thickness of plane standard mirror (14), and θ is the refraction angle of light when inciding plane standard mirror (14) front surface;

Total light field that photodetector receives is:

E(t)＝E ₁(t)+E ₂(t)+…+E _m(t)

Then the photocurrent of photodetector output is expressed as:

$I^{\&prime;}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{1}{Z}\underset{S}{\&Integral;\&Integral;}\frac{1}{2}[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]{[E_1\left(t\right)+E_2\left(t\right)+...+E_m\left(t\right)+...]}^{*}\mathrm{ds}$

Wherein, e is an electron charge, and Z is the intrinsic impedance of detector surface medium, and η is a quantum efficiency, and S is the area of detector photosurface, and h is a Planck's constant, and v is a laser frequency;

Above-mentioned photocurrent obtains second harmonic signal behind low-pass filter (17) electric current of intermediate frequency is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{2\mathrm{hv}}\frac{1}{Z}\underset{s}{\&Integral;\&Integral;}\underset{p=1}{\stackrel{\&infin;}{\mathrm{\&Sigma;}}}\underset{j=p+2}{\overset{\&infin;}{\mathrm{\&Sigma;}}}(E_p\left(t\right){E_j}^{*}\left(t\right)+{E_p}^{*}\left(t\right)E_j\left(t\right))\mathrm{ds}$

With the formula substitution following formula of all light fields, the arrangement result is:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p{E_0}^2\cos (\frac{8{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c}-\frac{4l{\mathrm{\&omega;}}_0\mathrm{and}\cos \mathrm{\&theta;}}{c^3}-\frac{8p{\mathrm{\&omega;}}_0{\mathrm{an}}^2{d}^2{\cos }^2\mathrm{\&theta;}}{c^3})$

Ignore l/c ³Event after can be reduced to:

$I_{\mathrm{if}}=\frac{\mathrm{\&eta;e}}{\mathrm{hv}}\frac{\mathrm{\&pi;}}{Z}{E_0}^2\cos (\frac{8a{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c^2}t-\frac{4{\mathrm{\&omega;}}_0\mathrm{nd}\cos \mathrm{\&theta;}}{c})\underset{p=1}{\overset{\&infin;}{\mathrm{\&Sigma;}}}{\mathrm{\&alpha;}}_{p+2}{\mathrm{\&alpha;}}_p$

In the formula, p and j are natural number;

And then the frequency that obtains interference signal is:

f＝8and?cosθω ₀/(2πc ²)＝4andω ₀/(πc ²)＝K?cosθ

Can know according to refraction law:

$\cos \mathrm{\&theta;}=\cos \left[\arcsin \left(\frac{\sin {\mathrm{\&theta;}}_0}{n}\right)\right]$

Then scale-up factor is:

K＝4andω ₀/(πc ²)

Then said momentum I is:

$I=\frac{k\&CenterDot;\arcsin \left\{n\sin \right[\arccos \left(\cos \mathrm{\&theta;}\right)\left]\right\}}{2}.$

2. the method for measuring little momentum based on the multi-beam laser heterodyne method of quadratic harmonics and the method for rocking according to claim 1; It is characterized in that measuring in the system of little momentum based on the multi-beam laser heterodyne second harmonic of the method for rocking; Digital information processing system (1) is made up of wave filter (17), prime amplifier (16), analog to digital converter (A/D) and digital signal processing controller (DSP); Wave filter (17) carries out the signal of photodetector (2) output to send to prime amplifier (16) after the filtering; This prime amplifier (16) sends to analog to digital converter (A/D) after the signal that receives is amplified; Analog to digital converter (A/D) becomes digital signal to send to digital signal processing controller (DSP) analog signal conversion that receives; Be solidified with fft algorithm in this digital signal processing controller (DSP), digital signal processing controller (DSP) is used for the signal that receives is continuously handled, and obtains the suffered little momentum of crossbeam of standard beam (8) after the demodulation.

Images