CN110940941B - Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect - Google Patents
Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect Download PDFInfo
- Publication number
- CN110940941B CN110940941B CN201911140997.XA CN201911140997A CN110940941B CN 110940941 B CN110940941 B CN 110940941B CN 201911140997 A CN201911140997 A CN 201911140997A CN 110940941 B CN110940941 B CN 110940941B
- Authority
- CN
- China
- Prior art keywords
- laser
- magnetic field
- sensing
- mode
- longitudinal
- 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.)
- Expired - Fee Related
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000000694 effects Effects 0.000 title claims abstract description 14
- 230000003287 optical effect Effects 0.000 claims abstract description 45
- 238000007781 pre-processing Methods 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000000691 measurement method Methods 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 39
- 239000013307 optical fiber Substances 0.000 claims description 37
- 239000000835 fiber Substances 0.000 claims description 23
- 230000006698 induction Effects 0.000 claims description 6
- 230000001902 propagating effect Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000002657 fibrous material Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000005355 Hall effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
The divisional application relates to the technical field of optical measurement, in particular to a magnetic field sensing measurement device and a magnetic field sensing measurement method based on a multi-longitudinal-mode self-mixing effect, wherein the measurement device comprises a multi-longitudinal-mode laser, a sensing unit, a vibration target, a sliding device, a light splitting element, a photoelectric detector, a signal preprocessing unit and a signal processing unit, and the measurement method comprises the following steps: the vibration target takes place to vibrate, many longitudinal mode laser outgoing laser incides on the vibration target behind the sensing element, then feedback back and form from the mixing signal in many longitudinal mode laser resonant cavities, above-mentioned in-process sensing element changes and causes from the mixing signal wave form change, make the vibration target take place to move a little through adjusting slider, form the self-mixing signal under different laser exocoel length, gather the self-mixing signal under the different exocoel length with photoelectric detector, reuse signal preprocessing unit and signal processing unit to handle, can obtain the change of sensing element, the present case measurement cost is low, the light path is simple, measurement accuracy is high.
Description
The application is divisional application with application number 201810327444.4, application date 2018, 4 and 12 months, and invention name "sensing measurement device and method based on multi-longitudinal mode self-mixing effect".
Technical Field
The invention relates to the technical field of optical measurement, in particular to a magnetic field sensing measurement device and method based on a multi-longitudinal-mode self-mixing effect.
Background
The optical measurement method is a main method in the technical field of measurement and measurement, and is mature to be applied to measurement occasions such as temperature measurement, voltage measurement, magnetic field measurement, strain measurement, liquid concentration measurement and the like at present due to the advantages of non-contact measurement, high measurement sensitivity, high measurement precision and the like.
In the field of magnetic field measurement technology, the conventional magnetic field sensor is generally realized by hall effect, Faraday magneto-optical effect, giant magneto-inductive effect, magnetic saturation effect and the like, but the methods generally have the problems of large volume, high cost, narrow measurement frequency band, small dynamic range and the like of a measurement system. With the development of optical sensing technology, optical magnetic sensors are receiving attention from researchers. The optical magnetic sensor mainly comprises a fiber grating magnetic field sensor, a Sagnac magnetic field sensor, a Michelson magnetic field sensor, a Mach-Zehnder magnetic field sensor and a Fabry-Perot magnetic field sensor. The fiber bragg grating magnetic field sensor and the Sagnac magnetic field sensor both need to be connected to a spectrometer to observe sensor output spectra under different magnetic field strengths, and are high in cost and easy to be influenced by the environment; the Michelson magnetic field sensor and the Mach-Zehnder magnetic field sensor acquire magnetic field intensity by collecting interference signals between a sensing arm and a reference arm, but signal light and reference light are in different optical paths, are greatly influenced by the environment, have complex structures and are difficult to debug; the fabry-perot type magnetic field sensor senses the magnetic field intensity by using the interference effect of light in the air cavity, but the air cavity is easily interfered by the environment and the optical path is limited to a certain extent, which is not beneficial to the measurement of the magnetic field intensity with high sensitivity.
Disclosure of Invention
Aiming at the problems existing in the prior art when the magnetic field is measured by using the optical sensing technology, the invention provides the sensing measuring device based on the multi-longitudinal-mode self-mixing effect, which can realize the sensing measurement of the magnetic field.
In order to realize the technical purpose of measuring the magnetic field, the technical scheme of the invention is as follows:
a magnetic field sensing measuring device based on a multi-longitudinal-mode self-mixing effect comprises a multi-longitudinal-mode laser with a tail fiber, a sensing unit, a vibration target, a sliding device, a light splitting element, a photoelectric detector, a signal preprocessing unit and a signal processing unit;
the vibration target can vibrate, and a reflection structure is attached to a vibration surface of the vibration target;
the sensing unit comprises a hysteresis telescopic material and a sensing optical fiber, the hysteresis telescopic material is arranged in a magnetic field to be measured, and the sensing optical fiber is fixed on the hysteresis telescopic material;
the multi-longitudinal-mode laser is used for emitting laser, a tail fiber of the multi-longitudinal-mode laser is connected with one end of a sensing optical fiber, the laser emitted from the other end of the sensing optical fiber is incident on a vibration surface of a vibration target, and is reflected by a reflecting structure and then fed back to the resonant cavity of the multi-longitudinal-mode laser along a primary circuit to form a laser self-mixing signal;
the bottom of the vibration target is fixed on the sliding device, and the vibration target can move along the laser emitting direction by adjusting the sliding device;
the light splitting element is a coupler and is used for splitting the laser self-mixing signal onto a photoelectric detector;
the photoelectric detector is used for converting the received laser signal into an electric signal and then sending the electric signal to the signal preprocessing unit;
the signal preprocessing unit is used for preprocessing the received electric signals, and the preprocessing comprises shaping, amplifying and filtering;
and the signal processing unit is used for analyzing and processing the preprocessed electric signals to obtain the magnetic field intensity of the magnetic field to be detected where the hysteresis telescopic material is located.
The magnetic field measuring method based on the measuring device comprises the following steps: the vibration target vibrates, the multi-longitudinal-mode laser emits laser to the vibration target, the emitted laser is reflected by the reflection structure and then fed back to the resonant cavity of the multi-longitudinal-mode laser to form a laser self-mixing signal, the magnetic field to be measured changes in the process, hysteresis materials change, and laser self-mixing signal waveforms change, the vibration target moves slightly along the direction of the light path of the emitted laser by adjusting the sliding device to change the distance between the vibration target and the multi-longitudinal-mode laser, so that the required laser self-mixing signals under different laser external cavity lengths are formed, the laser self-mixing signals under different laser external cavity lengths are collected by the photoelectric detector, then the signal preprocessing unit is used for preprocessing the laser self-mixing signals, and finally the signal processing unit is used for analyzing the preprocessed laser self-mixing signals, the magnetic field intensity of the magnetic field to be measured where the hysteresis telescopic material is located can be obtained, and the specific measurement and analysis method is as follows:
for laser self-mixing signals of a multi-longitudinal-mode laser, different longitudinal modes of the laser only interfere with the self-mode, the finally formed laser self-mixing signals are laser self-mixing signal intensity superposition formed by the respective longitudinal modes, and according to a related interference mixing theory model, under the condition of not considering speckle influence, the multi-longitudinal-mode laser self-mixing signal intensity is obtained:
beta in the formula (1) is the total number of oscillation starting modes in the multi-longitudinal-mode laser, j represents the jth longitudinal mode in the laser, I0Is the initial light intensity,. DELTA.IjAmplitude of variation of light intensity of j-mode laser, phitjFor the phase of the j-mode laser, phi, during one round trip of the outer cavitytj(t) is the real-time phase, k, of the round trip of the j-mode laser in the outer cavity0jWave number, op, of j mode in vacuumt(t) is the total optical path of the real-time external cavity of the laser, c.c. represents the complex conjugate of the formula, and the refractive index change caused by different longitudinal modes in the same material can be ignored in the calculation;
when the phase of the sensing unit changes, the external cavity total phase relationship is as follows:
phi in the formula (2)0jInitial phase of j-mode laser back and forth one revolution of the external cavity, delta phisjFor sensing unit phase changes caused by magnetic field changes, delta phicjFor compensating phase change, delta phi in measuring magnetic fieldsj=-δφcj,op0Is the initial optical path of the external cavity of the laser, delta opsFor changes in the optical path length of the sensing unit caused by changes in the magnetic field, delta opcTo compensate for optical path, ncIs the refractive index of air in the external cavity, and has a value of 1, nsFor sensing the refractive index of the optical fibre, LsIs a laser inTotal geometric length, L, of the actual path of transmission in the sensing fibercTo compensate for the length;
in the formula (3) < omega >0Is the angular frequency of the laser, c is the speed of light in vacuum, ngIs the refractive index of the laser resonant cavity medium group, L0Is the laser resonant cavity length;
substituting formula (3) into formula (1) to obtain:
if the waveforms of the laser self-mixing signals of different modes are not separated, the waveforms of the modes need to keep the same phase or the phase delay is integral multiple of 2 pi:
φtj=k0jopt=2mk0jngL0=mφgjformula (5)
Namely:
opt=2mngL0formula (6)
In the formula (5), m is the external cavity mode order of the laser and is a positive integer phigjThe phase of the laser is round trip in the resonant cavity of the laser, so the laser has a series of special position points, the superposed laser self-mixing signal does not generate waveform separation, and as can be known from the formula (5), when the intensity of the magnetic field to be measured changes, the phase of the light during transmission of the sensing optical fiber changes, resulting in phi of each modetjThe value m is not an integer, the waveform of the superposed laser self-mixing signal is separated, at the moment, the position of an external feedback object is changed by adjusting a sliding device to compensate the phase change, the waveform of the superposed laser self-mixing signal is changed into a complete waveform again, and then the position of the external feedback object is measured to obtain the compensation phase change delta phicjFurther obtaining the phase change delta phi of the sensing unit caused by the magnetic field changesjHere, strain of the hysteresis material causes a phase change δ φ in the sensor cellsjThe relationship of (a) is shown as follows:
in the formula (7), Ls0For the total initial geometric length of the actual path of the laser light propagating in the sensing fiber, ns0For the initial refractive index in the sensing fiber, nsIn order to sense the refractive index in the optical fiber,is the strain coefficient of the sensing optical fiber, v is the laser output frequency, a is the radius of the sensing optical fiber,for refractive index changes due to changes in the radius of the sensing fiber, this value is ignored in single mode fibers, and thus the strain induced phase change can be expressed as:
δφsj=k0jns0ξLs0epsilon formula (8)
The strain change of the magnetic induction material caused by the change of the magnetic field to be measured can be expressed as epsilon ═ CH2Where C is the expansion coefficient of the hysteresis expansion material and H is the magnetic field strength, the phase delay caused by the change in the magnetic field can be expressed as:
δφsj=k0jns0ξCH2Ls0formula (9)
Using compensating optical path in combination with sensing fiber material initial refractive index ns0The total initial geometric length L of the actual path of the laser transmitted in the sensing optical fibers0And calculating the expansion coefficient C of the magnetic hysteresis expansion material and the strain coefficient xi of the sensing optical fiber to obtain the magnetic field intensity of the magnetic field to be measured where the magnetic induction material is located.
From the above description, it can be seen that the present invention has the following advantages:
1. the sensing unit of the measuring device is a passive optical sensor, and power supply is not needed;
2. the test device has smaller volume and lower cost;
3. non-contact real-time high-precision measurement can be realized;
4. the optical path of the measuring device is a single optical path, the measuring device is small in environmental interference, simple in structure and convenient to adjust the optical path;
5. the measurement sensitivity and resolution can be adjusted by designing parameters of the sensing unit and selecting different external cavity measurement tools.
Drawings
FIG. 1 is a schematic structural diagram of an embodiment of the present invention;
fig. 2 is a diagram illustrating simulation results according to an embodiment of the present invention.
Detailed Description
Embodiment 1 of the present invention will be described in detail with reference to fig. 1 and 2, but the present invention is not limited to the claims.
As shown in fig. 1, a magnetic field sensing measuring device based on a multi-longitudinal-mode self-mixing effect includes a multi-longitudinal-mode laser 1 including a tail fiber, a sensing unit 2, a vibrating target 3, a sliding device 4, a light splitting element 5, a photodetector 6, a signal preprocessing unit 7 and a signal processing unit 8; the vibration target 3 can vibrate, and a reflection structure is attached to a vibration surface of the vibration target; the sensing unit 2 comprises a hysteresis telescopic material 23 and a sensing optical fiber 21, the hysteresis telescopic material 23 is arranged in a magnetic field to be measured, and the sensing optical fiber 21 is fixed on the hysteresis telescopic material 23; the multi-longitudinal-mode laser 1 is used for emitting laser, a tail fiber of the multi-longitudinal-mode laser is connected with one end of a sensing optical fiber 21, the laser emitted from the other end of the sensing optical fiber 21 is incident on a vibration surface of a vibration target 3, and is reflected by a reflecting structure and then fed back to a resonant cavity of the multi-longitudinal-mode laser 1 along a primary circuit to form a laser self-mixing signal; the bottom of the vibration target 3 is fixed on a sliding device 4, and the vibration target can move along the laser emitting direction by adjusting the sliding device 4; the light splitting element 5 adopts a coupler and is used for splitting the laser self-mixing signal onto the photoelectric detector 6; the photoelectric detector 6 is used for converting the received laser signal into an electric signal and then sending the electric signal to the signal preprocessing unit 7; the signal preprocessing unit 7 is used for preprocessing the received electric signals, and the preprocessing at least comprises shaping, amplifying and filtering; the signal processing unit 8 is configured to analyze and process the preprocessed electrical signals to obtain a magnetic field strength of a magnetic field to be measured where the hysteresis telescopic material is located.
In the above apparatus:
1. the sliding device 4 comprises a sliding rail 41 and a sliding block 42 arranged on the sliding rail, and the bottom of the vibration target is fixed on the sliding block 42; the slide rail 41 and the emergent laser are positioned on the same straight line;
2. the reflecting structure can be a reflecting plane mirror, and can also be a material with scattering property or reflecting property, such as a reflecting film;
3. the vibration target 3 may be a speaker 32 driven by a signal generator 31 or a piezoelectric ceramic, and two speakers 32 respectively shown by a solid line and a broken line in fig. 1 represent positions before and after the speaker slides along the sliding device;
4. the signal processing unit 8 may be a computer, an oscilloscope or a spectrometer.
The magnetic field measuring method based on the measuring device comprises the following steps: the vibration target vibrates, the multi-longitudinal-mode laser emits laser to the vibration target, the emitted laser is reflected by the reflection structure and then fed back to the resonant cavity of the multi-longitudinal-mode laser to form a laser self-mixing signal, the magnetic field to be measured changes in the process, hysteresis materials change, and laser self-mixing signal waveforms change, the vibration target moves slightly along the direction of the light path of the emitted laser by adjusting the sliding device to change the distance between the vibration target and the multi-longitudinal-mode laser, so that the required laser self-mixing signals under different laser external cavity lengths are formed, the laser self-mixing signals under different laser external cavity lengths are collected by the photoelectric detector, then the signal preprocessing unit is used for preprocessing the laser self-mixing signals, and finally the signal processing unit is used for analyzing the preprocessed laser self-mixing signals, the magnetic field intensity of the magnetic field to be measured where the hysteresis telescopic material is located can be obtained, and the specific measurement and analysis method is as follows:
for laser self-mixing signals of a multi-longitudinal-mode laser, different longitudinal modes of the laser only interfere with the self-mode, the finally formed laser self-mixing signals are laser self-mixing signal intensity superposition formed by the respective longitudinal modes, and according to a related interference mixing theory model, under the condition of not considering speckle influence, the multi-longitudinal-mode laser self-mixing signal intensity is obtained:
beta in the formula (1) is the total number of oscillation starting modes in the multi-longitudinal-mode laser, j represents the jth longitudinal mode in the laser, I0Is the initial light intensity,. DELTA.IjAmplitude of variation of light intensity of j-mode laser, phitjFor the phase of the j-mode laser, phi, during one round trip of the outer cavitytj(t) is the real-time phase, k, of the round trip of the j-mode laser in the outer cavity0jWave number, op, of j mode in vacuumt(t) is the real-time external cavity total optical path c.c. of the laser, which represents the complex conjugate of the formula, and the refractive index change caused by different longitudinal modes in the same material can be ignored in the calculation;
when the phase of the sensing unit changes, the external cavity total phase relationship is as follows:
phi in the formula (2)0jInitial phase of j-mode laser back and forth one revolution of the external cavity, delta phisjFor sensing unit phase changes caused by magnetic field changes, delta phicjFor compensating phase change, delta phi in measuring magnetic fieldsj=-δφcj,op0Is the initial optical path of the external cavity of the laser, delta opsFor changes in the optical path length of the sensing unit caused by changes in the magnetic field, delta opcTo compensate for optical path, ncIs the refractive index of air in the external cavity, and has a value of 1, nsFor sensing the refractive index of the optical fibre, LsFor the total geometrical length, L, of the actual path of the laser light propagating in the sensing fibercTo compensate for the length;
in the formula (3) < omega >0Is the angular frequency of the laser, c is the speed of light in vacuum, ngIs the refractive index of the laser resonant cavity medium group, L0Is the laser resonant cavity length;
substituting formula (3) into formula (1) to obtain:
if the waveforms of the laser self-mixing signals of different modes are not separated, the waveforms of the modes need to keep the same phase or the phase delay is integral multiple of 2 pi:
φtj=k0jopt=2mk0jngL0=mφgjformula (5)
Namely:
opt=2mngL0formula (6)
In the formula (5), m is the external cavity mode order of the laser and is a positive integer phigjThe phase of the laser is round trip in the resonant cavity of the laser, so the laser has a series of special position points, the superposed laser self-mixing signal does not generate waveform separation, and as can be known from the formula (5), when the intensity of the magnetic field to be measured changes, the phase of the light during transmission of the sensing optical fiber changes, resulting in phi of each modetjThe value m is not an integer, the waveform of the superposed laser self-mixing signal is separated, at the moment, the position of an external feedback object is changed by adjusting a sliding device to compensate the phase change, the waveform of the superposed laser self-mixing signal is changed into a complete waveform again, and then the position of the external feedback object is measured to obtain the compensation phase change delta phicjFurther obtaining the phase change delta phi of the sensing unit caused by the magnetic field changesjHere, strain of the hysteresis material causes a phase change δ φ in the sensor cellsjThe relationship of (a) is shown as follows:
in the formula (7), Ls0For the total initial geometric length of the actual path of the laser light propagating in the sensing fiber, ns0For the initial refractive index in the sensing fiber, nsIn order to sense the refractive index in the optical fiber,is the strain coefficient of the sensing optical fiber, v is the laser output frequency, a is the radius of the sensing optical fiber,for refractive index changes due to changes in the radius of the sensing fiber, this value is ignored in single mode fibers, and thus the strain induced phase change can be expressed as:
δφsj=k0jns0ξLs0epsilon formula (8)
The strain change of the magnetic induction material caused by the change of the magnetic field to be measured can be expressed as epsilon ═ CH2Where C is the expansion coefficient of the hysteresis expansion material and H is the magnetic field strength, the phase delay caused by the change in the magnetic field can be expressed as:
δφsj=k0jns0ξCH2Ls0formula (9)
Using compensating optical path in combination with sensing fiber material initial refractive index ns0The total initial geometric length L of the actual path of the laser transmitted in the sensing optical fibers0And calculating the expansion coefficient C of the magnetic hysteresis expansion material and the strain coefficient xi of the sensing optical fiber to obtain the magnetic field intensity of the magnetic field to be measured where the magnetic induction material is located.
Based on the above technical solution, an experimental apparatus is established, the experimental apparatus is a dual-mode LD laser, and simulation software is used for analog simulation, for simplicity, we only consider the intensity superposition waveform of the dual-mode LD laser self-mixing signal with the same amplitude, and specific simulation parameters are as follows:op0=14574.00mm,δopc=0mm,Ls0=10m,n1=1.45,ng=3.5,L0=300um,C=6.9×10-15A-2m-2the magnetic field strength is increased by 30k A/m.
Simulation As shown in FIG. 2, it can be seen from FIG. 2 that when the magnetic field strength is 0, the initial optical path of the external cavity of the laser is 14574.00mm, which is ngL0M is 13880, and the laser self-mixing signal waveform is not separated. When the magnetic field intensity of the sensing unit is increased by 30k A/m, the optical path of the sensing unit is slightly changed due to the change of the magnetic field intensity, the waveform of the overlapped laser self-mixing signal is separated, the length of the fine-tuning compensation external cavity is 0.727mm, and the optical path of the laser external cavity becomes n againgL0The wave form of the laser self-mixing signal after superposition disappears separately, and finally the magnetic field change of the sensing unit is obtained through measuring the compensation optical path, thereby realizing the magnetic field measurement.
The phase change sensitivity S of the magnetic field sensor can be further obtained by the formula (7)mLcAnd adjacent magnetic field intensity difference Delta Hm。SmLcAnd Δ HmThe refractive index of the sensing unit optical fiber material, the length of the sensing unit optical fiber, the expansion coefficient of the hysteresis expansion material and the strain coefficient of the optical fiber are determined together. Wherein the external cavity variation sensitivity SmLcMeans the length change of the compensation external cavity caused by the change of unit magnetic field intensity and the magnetic field intensity difference Delta H of adjacent levelmRefers to the magnetic field intensity H2Position of external cavity equiphase point (m +1 level) and magnetic field intensity H caused by (after change)1And (before change) the magnetic field intensity difference value of the adjacent level corresponding to the position (m level) of the external cavity equiphase point. In general, in the magnetic field measurement process, if the difference of the magnetic field intensity measured in two consecutive measurement intervals is larger than the adjacent magnetic field intensity difference Δ HmThe number of cycles of the waveform change of the self-mixing signal, i.e. the change of the value m, in two consecutive measurement intervals is recorded, and the length of the compensation external cavity is adjusted to restore the waveform of the laser self-mixing signal to the position where the waveform of the signal corresponding to the mth level coincides.
Expressions (10) and (11) are external cavity variation sensitivity SmLcAnd adjacent levelMagnetic field intensity difference Δ HmExpression:
when the device of this embodiment is adopted to measure the magnetic field intensity, the following advantages are provided:
1. the sensing unit of the measuring device is a passive optical sensor, and power supply is not needed;
2. the test device has smaller volume and lower cost;
3. non-contact real-time high-precision measurement can be realized;
4. the optical path of the measuring device is a single optical path, the measuring device is small in environmental interference, simple in structure and convenient to adjust the optical path;
5. the sensitivity and resolution of magnetic field intensity measurement can be adjusted by designing parameters of the sensing unit and selecting different external cavity measuring tools. In order to improve the performance of the measuring device in the embodiment, the following improvements can be made to the device:
1. an optical attenuator 9 is added to an optical path between the spectroscopic element 5 and the vibration target 3, and the intensity of the optical feedback light is adjusted by the optical attenuator 9.
2. The multi-longitudinal-mode laser 1 adopts a semiconductor laser, and integrates a photodiode into the semiconductor laser by utilizing the characteristics of the semiconductor laser to realize the function of a photoelectric detector, so that the optical path of the whole device is simplified, and a light splitting element and the photoelectric detector are removed;
3. in order to improve the collimation performance of the emitted laser, the other end of the sensing optical fiber is connected with a collimator 10, and the laser is ensured to be emitted to a vibration target in parallel through the collimator 10.
In summary, the invention has the following advantages:
1. the sensing unit of the measuring device is a passive optical sensor, and power supply is not needed;
2. the test device has smaller volume and lower cost;
3. non-contact real-time high-precision measurement can be realized;
4. the optical path of the measuring device is a single optical path, the measuring device is small in environmental interference, simple in structure and convenient to adjust the optical path;
5. the measurement sensitivity and resolution can be adjusted by designing parameters of the sensing unit and selecting different external cavity measurement tools.
It should be understood that the detailed description of the invention is merely illustrative of the invention and is not intended to limit the invention to the specific embodiments described. It will be appreciated by those skilled in the art that the present invention may be modified or substituted equally as well to achieve the same technical result; as long as the use requirements are met, the method is within the protection scope of the invention.
Claims (2)
1. The utility model provides a magnetic field sensing measuring device based on many longitudinal modes are from mixing effect which characterized in that: the device comprises a multi-longitudinal-mode laser with tail fibers, a sensing unit, a vibration target, a sliding device, a light splitting element, a photoelectric detector, a signal preprocessing unit and a signal processing unit;
the vibration target can vibrate, and a reflection structure is attached to a vibration surface of the vibration target;
the sensing unit comprises a hysteresis telescopic material and a sensing optical fiber, the hysteresis telescopic material is arranged in a magnetic field to be measured, and the sensing optical fiber is fixed on the hysteresis telescopic material;
the multi-longitudinal-mode laser is used for emitting laser, a tail fiber of the multi-longitudinal-mode laser is connected with one end of a sensing optical fiber, the laser emitted from the other end of the sensing optical fiber is incident on a vibration surface of a vibration target, and is reflected by a reflecting structure and then fed back to the resonant cavity of the multi-longitudinal-mode laser along a primary circuit to form a laser self-mixing signal;
the bottom of the vibration target is fixed on the sliding device, and the vibration target can move along the laser emitting direction by adjusting the sliding device;
the light splitting element is a coupler and is used for splitting the laser self-mixing signal onto a photoelectric detector;
the photoelectric detector is used for converting the received laser signal into an electric signal and then sending the electric signal to the signal preprocessing unit;
the signal preprocessing unit is used for preprocessing the received electric signals, and the preprocessing comprises shaping, amplifying and filtering;
and the signal processing unit is used for analyzing and processing the preprocessed electric signals to obtain the magnetic field intensity of the magnetic field to be detected where the hysteresis telescopic material is located.
2. The magnetic field measurement method of the magnetic field sensing measurement device based on the multi-longitudinal-mode self-mixing effect according to claim 1, characterized in that: the vibration target vibrates, the multi-longitudinal-mode laser emits laser to the vibration target, the emitted laser is reflected by the reflection structure and then fed back to the resonant cavity of the multi-longitudinal-mode laser to form a laser self-mixing signal, the magnetic field to be measured changes in the process, hysteresis materials change, and laser self-mixing signal waveforms change, the vibration target moves slightly along the direction of the light path of the emitted laser by adjusting the sliding device to change the distance between the vibration target and the multi-longitudinal-mode laser, so that the required laser self-mixing signals under different laser external cavity lengths are formed, the laser self-mixing signals under different laser external cavity lengths are collected by the photoelectric detector, then the signal preprocessing unit is used for preprocessing the laser self-mixing signals, and finally the signal processing unit is used for analyzing the preprocessed laser self-mixing signals, the magnetic field intensity of the magnetic field to be measured where the hysteresis telescopic material is located can be obtained, and the specific measurement and analysis method is as follows:
for laser self-mixing signals of a multi-longitudinal-mode laser, different longitudinal modes of the laser only interfere with the self-mode, the finally formed laser self-mixing signals are laser self-mixing signal intensity superposition formed by the respective longitudinal modes, and according to a related interference mixing theory model, under the condition of not considering speckle influence, the multi-longitudinal-mode laser self-mixing signal intensity is obtained:
beta in the formula (1) is the total number of oscillation starting modes in the multi-longitudinal-mode laser, j represents the jth longitudinal mode in the laser, I0Is the initial light intensity,. DELTA.IjAmplitude of variation of light intensity of j-mode laser, phitjFor the phase of the j-mode laser, phi, during one round trip of the outer cavitytj(t) is the real-time phase, k, of the round trip of the j-mode laser in the outer cavity0jWave number, op, of j mode in vacuumt(t) is the total optical path of the real-time external cavity of the laser, c.c. represents the complex conjugate of the formula, and the refractive index change caused by different longitudinal modes in the same material can be ignored in the calculation;
when the phase of the sensing unit changes, the external cavity total phase relationship is as follows:
phi in the formula (2)0jInitial phase of j-mode laser back and forth one revolution of the external cavity, delta phisjFor sensing unit phase changes caused by magnetic field changes, delta phicjFor compensating phase change, delta phi in measuring magnetic fieldsj=-δφcj,op0Is the initial optical path of the external cavity of the laser, delta opsFor changes in the optical path length of the sensing unit caused by changes in the magnetic field, delta opcTo compensate for optical path, ncIs the refractive index of air in the external cavity, and has a value of 1, nsFor sensing the refractive index of the optical fibre, LsFor the total geometrical length, L, of the actual path of the laser light propagating in the sensing fibercTo compensate for the length;
in the formula (3) < omega >0Is the angular frequency of the laser, c is the speed of light in vacuum, ngIs the refractive index of the laser resonant cavity medium group, L0Is a laser resonant cavityLength;
substituting formula (3) into formula (1) to obtain:
if the waveforms of the laser self-mixing signals of different modes are not separated, the waveforms of the modes need to keep the same phase or the phase delay is integral multiple of 2 pi:
φtj=k0jopt=2mk0jngL0=mφgjformula (5)
Namely:
opt=2mngL0formula (6)
In the formula (5), m is the external cavity mode order of the laser and is a positive integer phigjThe phase of the laser is round trip in the resonant cavity of the laser, so the laser has a series of special position points, the superposed laser self-mixing signal does not generate waveform separation, and as can be known from the formula (5), when the intensity of the magnetic field to be measured changes, the phase of the light during transmission of the sensing optical fiber changes, resulting in phi of each modetjThe value m is not an integer, the waveform of the superposed laser self-mixing signal is separated, at the moment, the position of an external feedback object is changed by adjusting a sliding device to compensate the phase change, the waveform of the superposed laser self-mixing signal is changed into a complete waveform again, and then the position of the external feedback object is measured to obtain the compensation phase change delta phicjFurther obtaining the phase change delta phi of the sensing unit caused by the magnetic field changesjHere, strain of the hysteresis material causes a phase change δ φ in the sensor cellsjThe relationship of (a) is shown as follows:
in the formula (7), Ls0For the total initial geometric length of the actual path of the laser light propagating in the sensing fiber, ns0For the initial refractive index in the sensing fiber, nsIn order to sense the refractive index in the optical fiber,is the strain coefficient of the sensing optical fiber, v is the laser output frequency, a is the radius of the sensing optical fiber,for refractive index changes due to changes in the radius of the sensing fiber, this value is ignored in single mode fibers, and thus the strain induced phase change can be expressed as:
δφsj=k0jns0ξLs0epsilon formula (8)
The strain change of the magnetic induction material caused by the change of the magnetic field to be measured can be expressed as epsilon ═ CH2Where C is the expansion coefficient of the hysteresis expansion material and H is the magnetic field strength, the phase delay caused by the change in the magnetic field can be expressed as:
δφsj=k0jns0ξCH2Ls0formula (9)
Using compensating optical path in combination with sensing fiber material initial refractive index ns0The total initial geometric length L of the actual path of the laser transmitted in the sensing optical fibers0And calculating the expansion coefficient C of the magnetic hysteresis expansion material and the strain coefficient xi of the sensing optical fiber to obtain the magnetic field intensity of the magnetic field to be measured where the magnetic induction material is located.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911140997.XA CN110940941B (en) | 2018-04-12 | 2018-04-12 | Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810327444.4A CN108775974B (en) | 2018-04-12 | 2018-04-12 | Temperature sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201911140997.XA CN110940941B (en) | 2018-04-12 | 2018-04-12 | Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810327444.4A Division CN108775974B (en) | 2018-04-12 | 2018-04-12 | Temperature sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110940941A CN110940941A (en) | 2020-03-31 |
CN110940941B true CN110940941B (en) | 2022-03-25 |
Family
ID=64033741
Family Applications (5)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911140972.XA Expired - Fee Related CN110806397B (en) | 2018-04-12 | 2018-04-12 | Liquid concentration sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201911140997.XA Expired - Fee Related CN110940941B (en) | 2018-04-12 | 2018-04-12 | Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201911140996.5A Expired - Fee Related CN110850144B (en) | 2018-04-12 | 2018-04-12 | Voltage sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201810327444.4A Expired - Fee Related CN108775974B (en) | 2018-04-12 | 2018-04-12 | Temperature sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201911141654.5A Active CN110806274B (en) | 2018-04-12 | 2018-04-12 | Strain sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911140972.XA Expired - Fee Related CN110806397B (en) | 2018-04-12 | 2018-04-12 | Liquid concentration sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Family Applications After (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911140996.5A Expired - Fee Related CN110850144B (en) | 2018-04-12 | 2018-04-12 | Voltage sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201810327444.4A Expired - Fee Related CN108775974B (en) | 2018-04-12 | 2018-04-12 | Temperature sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
CN201911141654.5A Active CN110806274B (en) | 2018-04-12 | 2018-04-12 | Strain sensing measurement device and method based on multi-longitudinal-mode self-mixing effect |
Country Status (1)
Country | Link |
---|---|
CN (5) | CN110806397B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109497947B (en) * | 2018-11-19 | 2024-05-17 | 天津大学 | Portable OCT measuring device |
CN109724648B (en) * | 2019-01-21 | 2020-09-25 | 安徽大学 | Device and method for synchronously measuring temperature and strain based on orthogonal polarization dual-wavelength laser multi-longitudinal-mode self-mixing effect |
CN113687279B (en) * | 2021-08-11 | 2022-10-11 | 北京航天控制仪器研究所 | Reflection type optical path adjustable CPT atomic magnetometer probe |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104677296A (en) * | 2015-03-09 | 2015-06-03 | 北京交通大学 | System for measurement of displacement through self-mixing interference fusion of beat waves and single waves of fiber laser |
CN105716704A (en) * | 2016-04-20 | 2016-06-29 | 安徽大学 | Micro-cavity chip type laser self-mixing vibration, displacement and speed sensing method and system |
Family Cites Families (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57133531A (en) * | 1981-02-12 | 1982-08-18 | Agency Of Ind Science & Technol | Optical information processor |
US4827483A (en) * | 1985-08-12 | 1989-05-02 | Hitachi, Ltd. | Semiconductor laser device and method of fabricating the same |
US5838439A (en) * | 1997-03-14 | 1998-11-17 | Zang; De Yu | Heterodyned self-mixing laser diode vibrometer |
CN100489752C (en) * | 2004-02-09 | 2009-05-20 | 皇家飞利浦电子股份有限公司 | Optical input device based on Doppler frequency shift and laser self-mixing |
CN100386930C (en) * | 2005-03-04 | 2008-05-07 | 清华大学 | HeNe laser device capable of outputting two beams of vertically polarized light with stabilized intensity |
WO2006126152A1 (en) * | 2005-05-24 | 2006-11-30 | Koninklijke Philips Electronics N.V. | Glucose sensor |
CN101305337B (en) * | 2005-08-30 | 2012-09-05 | 皇家飞利浦电子股份有限公司 | Method of measuring relative movement in two dimensions of an object and an optical input device |
CN100451536C (en) * | 2007-01-19 | 2009-01-14 | 清华大学 | Quasi-common path type feedback interferometer of laser in microchip |
CN100538397C (en) * | 2007-03-16 | 2009-09-09 | 清华大学 | Double-refraction external cavity displacement measuring system |
JP2011527750A (en) * | 2008-07-10 | 2011-11-04 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Gas detector |
EP2166328A1 (en) * | 2008-09-22 | 2010-03-24 | Universita' degli studi di Bari | System for optical fiber strain measure |
US8550958B2 (en) * | 2009-03-31 | 2013-10-08 | GM Global Technology Operations LLC | Shift control method for a multi-mode hybrid transmission |
CN101539454A (en) * | 2009-04-01 | 2009-09-23 | 南京师范大学 | Semiconductor laser self-mixing interference vibration meter |
KR101701711B1 (en) * | 2009-08-10 | 2017-02-03 | 코닌클리케 필립스 엔.브이. | Vertical cavity surface emitting laser with active carrier confinement |
CN102427573A (en) * | 2011-12-09 | 2012-04-25 | 安徽大学 | Real-time voice signal receiving device based on self-mixing interference measuring method |
EP2654143A1 (en) * | 2012-04-18 | 2013-10-23 | Nederlandse Organisatie voor toegepast -natuurwetenschappelijk onderzoek TNO | Frequency tunable laser system |
CN103018747A (en) * | 2012-11-26 | 2013-04-03 | 中国人民解放军第四军医大学 | Laser self-mixing distance measuring system based on uncollimated feedback external cavity |
CN103308147A (en) * | 2013-06-07 | 2013-09-18 | 中国人民解放军第四军医大学 | Vibration measurement method and system based on single-frequency laser feedback |
CN103337776B (en) * | 2013-06-09 | 2015-07-15 | 安徽大学 | All-optical fiber type self-mixing distance measuring system of laser |
CN103344184B (en) * | 2013-06-09 | 2015-11-25 | 安徽大学 | Based on the wavelength-division of the mixing certainly multiplexed multi-channel displacement sensing system of linear cavity multi-wavelength optical fiber laser |
CN105829844B (en) * | 2013-08-22 | 2018-05-01 | 昆士兰大学 | For being imaged the Optical Maser System with material analysis |
CN103439268B (en) * | 2013-09-05 | 2016-05-04 | 安徽大学 | High-sensitivity diaphragm type photoacoustic spectrum sensor based on self-mixing |
CN103575313B (en) * | 2013-11-21 | 2016-08-31 | 黑龙江大学 | Many longitudinal modes annular chamber laser sensor frequency division multiplexing device based on beat frequency technology |
CN104729424B (en) * | 2015-03-09 | 2017-10-20 | 南京信息工程大学 | Confocal dot laser microscope and its scan method based on self-mixed interference |
US10327069B2 (en) * | 2015-07-26 | 2019-06-18 | Vocalzoom Systems Ltd. | Laser microphone utilizing speckles noise reduction |
CN107850528B (en) * | 2015-07-30 | 2020-08-11 | 通快光电器件有限公司 | Laser sensor for multi-parameter detection |
CN105486224A (en) * | 2015-11-25 | 2016-04-13 | 上海珍岛信息技术有限公司 | Laser self-mixing interference measuring system |
CN105547197B (en) * | 2015-12-10 | 2017-12-19 | 中国科学院合肥物质科学研究院 | Measurement angle and the method and device of vibration while based on laser self-mixing interference |
CN205505978U (en) * | 2016-03-31 | 2016-08-24 | 天津大学 | Small -size optical measuring head based on laser is interfered from mixing |
CN105823755B (en) * | 2016-04-07 | 2019-02-15 | 南京先进激光技术研究院 | It is a kind of that sensor-based system is absorbed from mixed gas based on tunable semiconductor laser |
CN109818245B (en) * | 2016-04-20 | 2020-03-31 | 安徽大学 | Micro-cavity chip type laser self-mixing distance sensing system |
DE102016210925A1 (en) * | 2016-06-20 | 2017-12-21 | Robert Bosch Gmbh | System and method for interior sensing in a vehicle |
CN106441144A (en) * | 2016-06-24 | 2017-02-22 | 天津大学 | Three-dimensional morphology measurement method for inhibiting speckles on the basis of wide-spectrum laser |
CN106444068A (en) * | 2016-06-24 | 2017-02-22 | 天津大学 | Laser speckle inhibition method |
CN106443551B (en) * | 2016-10-12 | 2019-09-10 | 长飞光纤光缆股份有限公司 | A kind of all-fiber current transformator monitoring system |
CN106802165A (en) * | 2017-03-02 | 2017-06-06 | 阜阳师范学院 | Speed and distance synchronous measuring method and device based on laser self-mixing interference |
CN107045161B (en) * | 2017-03-28 | 2023-11-28 | 南京信息职业技术学院 | Multi-longitudinal-mode fiber laser sensing multiplexing system |
CN106997051B (en) * | 2017-06-05 | 2023-04-04 | 安徽大学 | Laser vector wind measurement method and wind measurement radar based on polarization effect and self-mixing effect |
CN110132180B (en) * | 2017-10-12 | 2021-01-05 | 安徽大学 | Mirror-surface type laser self-mixing micro-angle measuring system and method with any included angle |
CN110806306B (en) * | 2018-04-12 | 2021-03-30 | 安徽大学 | Device and method for measuring temperature change of cavity of resonant cavity of multi-longitudinal-mode laser |
-
2018
- 2018-04-12 CN CN201911140972.XA patent/CN110806397B/en not_active Expired - Fee Related
- 2018-04-12 CN CN201911140997.XA patent/CN110940941B/en not_active Expired - Fee Related
- 2018-04-12 CN CN201911140996.5A patent/CN110850144B/en not_active Expired - Fee Related
- 2018-04-12 CN CN201810327444.4A patent/CN108775974B/en not_active Expired - Fee Related
- 2018-04-12 CN CN201911141654.5A patent/CN110806274B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104677296A (en) * | 2015-03-09 | 2015-06-03 | 北京交通大学 | System for measurement of displacement through self-mixing interference fusion of beat waves and single waves of fiber laser |
CN105716704A (en) * | 2016-04-20 | 2016-06-29 | 安徽大学 | Micro-cavity chip type laser self-mixing vibration, displacement and speed sensing method and system |
Also Published As
Publication number | Publication date |
---|---|
CN110940941A (en) | 2020-03-31 |
CN110806397A (en) | 2020-02-18 |
CN108775974B (en) | 2020-03-31 |
CN110806274A (en) | 2020-02-18 |
CN110850144A (en) | 2020-02-28 |
CN110806397B (en) | 2022-01-11 |
CN110806274B (en) | 2020-12-15 |
CN110850144B (en) | 2022-02-18 |
CN108775974A (en) | 2018-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110806306B (en) | Device and method for measuring temperature change of cavity of resonant cavity of multi-longitudinal-mode laser | |
CN110940941B (en) | Magnetic field sensing measurement device and method based on multi-longitudinal-mode self-mixing effect | |
CN109724648B (en) | Device and method for synchronously measuring temperature and strain based on orthogonal polarization dual-wavelength laser multi-longitudinal-mode self-mixing effect | |
CN110631514B (en) | Pentagonal prism type angle sensing measurement device and method based on multi-longitudinal mode self-mixing effect | |
CN103364775A (en) | Optical frequency comb calibration-based dual-color laser scanning absolute distance measuring device and method | |
CN101825432A (en) | Dual-wavelength optical-fiber interference large-range high-resolution displacement measurement system | |
CN111693255B (en) | Device and method for measuring frequency drift of laser light source | |
CN102176022B (en) | Method for measuring young modulus by adopting multi-beam laser heterodyne distance measurement device | |
Sun et al. | Multimode interference-based fiber-optic ultrasonic sensor for non-contact displacement measurement | |
JP2001330669A (en) | Laser diode-type distance and displacement meter with double external resonator | |
CN100520281C (en) | Adjustable frequency difference two-frequency laser back-feed nano measuring ruler for frequency stabilizing | |
CN111579816B (en) | Acceleration measuring instrument based on photoelectric oscillator | |
CN103226205B (en) | Optical fiber sensing measurement method of laser plasma shock wave mechanical effect | |
Huang et al. | Ultra-high-resolution phase demodulation based miniature fiber-optic accelerometer at low and medium frequencies | |
CN101738167A (en) | Resonant cavity frequency stabilization-based absolute distance measurement system and implementing method thereof | |
Di Sante et al. | A novel fiber optic sensor for multiple and simultaneous measurement of vibration velocity | |
Sathitanon et al. | A fiber optic interferometric sensor for dynamic measurement | |
CN108709717B (en) | Device and method for measuring resonant cavity FSR of multi-longitudinal-mode laser by using large-amplitude laser self-mixing vibration signal | |
Suleiman et al. | Interrogation of fiber Bragg grating dynamic strain sensors by self-mixing interferometry | |
Wu et al. | A novel optical accelerometer based on slant-ended fiber | |
Pullteap et al. | Modified fringe-counting technique applied to a dual-cavity fiber Fabry-Pérot vibrometer | |
CN108037143B (en) | Method and device for measuring refractive index of gas | |
Ivaschescu | Small sinusoidal vibrations amplitude measurements with the Michelson interferometer | |
CN112098039B (en) | System and method for measuring pulsation density of hypersonic flow field | |
CN100427905C (en) | All-optical fiber spectral measurement method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220325 |
|
CF01 | Termination of patent right due to non-payment of annual fee |