CN107092756A - A kind of angular-rate sensor modeling method based on MHD effect - Google Patents
A kind of angular-rate sensor modeling method based on MHD effect Download PDFInfo
- Publication number
- CN107092756A CN107092756A CN201710282612.8A CN201710282612A CN107092756A CN 107092756 A CN107092756 A CN 107092756A CN 201710282612 A CN201710282612 A CN 201710282612A CN 107092756 A CN107092756 A CN 107092756A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- sensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
Abstract
The invention discloses a kind of angular-rate sensor modeling method based on MHD effect, comprise the steps of：Step S1, sets up the TRANSFER MODEL of sensor probe；Step S2, sets up the transmission function of primary amplifying transformer；Step S3, sets up the transmission characteristic of rear end instrument amplifying circuit；Step S4, according to the result obtained by step S1, S2 and S3, builds the block mold of sensor.Present system gives sensor probe model, transformer model, the sensor entire process model that rear end amplifying circuit model and three are constituted.This method can quantitative analysis go out under some design parameters, the corresponding Frequency Response of sensor；Also Selection and Design can be optimized to the relevant parameter in sensor, directive function is played in actual development process to sensor by its TRANSFER MODEL.
Description
Technical field
The present invention relates to the in-orbit wideband attitude measurement field of spacecraft, and in particular to one kind is based on MHD effect
Angular-rate sensor modeling method.
Background technology
Camera optical axis measurement accuracy is not enough with bandwidth during camera imaging, will cause image geometry Quality Down, influences
Satellite is without control positioning precision.Platform stance situation of change within the measurable 0~10Hz of traditional measurement sensor, and satellite body and
The structural vibration of its annex, the operating micro-vibration of executing agency all correspond to 0.001Hz~500Hz wideband information.Load diagram
As the thousands of Hz of correspondence data message.Using 2~1000Hz high bandwidth angular-rate sensors, satellite high-precision attitude jitter is carried out
Measurement, can expand payload platform posture and determine bandwidth and improve attitude determination accuracy level.
Current high frequency angular oscillation measuring method has following several：(1) angular oscillation is obtained using the combination of multi-thread vibrating sensor
Information, this method is needed to be calculated indirectly, and number of sensors is more, and the raising of measurement accuracy is limited.(2) optical fibre gyro is utilized
Measurement angle vibration information, this mode precision is high, single expensive, is not suitable for civil applications development.(3) magnetic current bulk effect is utilized
Angular-rate sensor carries out vibration measurement, and such sensor bulk is small, light weight, and precision is high, has wide civil applications to be worth, closely
Ji Nian China has just carried out the art research work, and development time is relatively short.
The content of the invention
It is an object of the invention to provide a kind of angular-rate sensor modeling method based on MHD effect, the party
Method can obtain sensor block mold according to sensor relevant design parameter, be the estimation of sensor amplitude versus frequency characte and sensor ginseng
Number optimized Selection design lays the foundation.
To reach above-mentioned purpose, the invention provides a kind of angular-rate sensor modeling based on MHD effect
Method, is comprised the steps of：
Step S1, sets up the TRANSFER MODEL of sensor probe；
Step S2, sets up the transmission function of primary amplifying transformer；
Step S3, sets up the transmission characteristic of rear end instrument amplifying circuit；
Step S4, according to the result obtained by step S1, S2 and S3, builds the block mold of sensor.
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, the step S1 is specific
Comprise the steps of：
Step S11, the electric-field intensity inside conductor fluid is obtained according to Ohm's law：
J=σ (E+V × B) (1)
J is current strength in formula, and E is electric-field intensity, and B is magnetic induction intensity, and σ is electrical conductivity；When under rectangular coordinate system
Magnetic field condition is Bx=0, By=-B0,Bz=0, then it can obtain z-axis directional current intensity：
U in formulaiFor lower plate speed, u is the movement velocity of conductor fluid infinitesimal, and r is O points and lower plate distance；
Step S12, calculates the electromagnetic force for acting on conductor fluid infinitesimal, Fe=JyBz-JzBy=B0Jz；
Step S13, according to Hartmann's constant physical significance, calculating acts on the viscous force of conductor fluid micelle；
I.e.FuFor the viscous force of micelle, M is Hartmann's constant,Wherein h is conductor fluid
Ring thickness, η, ρ, v are respectively the resistivity, density, kinematic viscosity coefficient of conductor fluid；
Step S14, electromagnetic force and viscous force according to suffered by kinematical equation and fluid infinitesimal, obtains conductor fluid micro-
The movement velocity of group：
ρ is the density of conducting liquid in formula；
Step S15, according to the law of electromagnetic induction and the movement velocity of fluid infinitesimal, obtains probe output voltage and angular speed
Transitive relation：
L is electrical conduction current body cavity height, r in formulaRMSFor the root mean square radii of conductor fluid ring.
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, in the step S2, structure
The equivalent circuit of transformer is built, according to the transmission function of Kirchhoff's law to transformer：
L in formula1For primary inductance, L2For secondary inductance, L12For primary secondary coil mutual inductance, R1For primary line
The resistance of circle, R2For the resistance of secondary coil, C2For the parasitic capacitance of secondary coil.
The above-mentioned angular speed based on MHD effect passes speed modeling method, wherein, in the step S3,
The transmission characteristic of rear end instrument amplifying circuit is regard as constant value K in the frequency range of sensor responseU。
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, will in the step S4
Sensor probe, primary amplifying transformer transmission function corresponding with amplifying circuit with rear end instrument combines, and is sensed
Device entirety TRANSFER MODEL：
In formula：
The advantageous effect of present invention is that：There is provided a kind of angular-rate sensor based on MHD effect
Modeling method, system gives sensor probe model, transformer model, the sensing that rear end amplifying circuit model and three are constituted
Device entire process model.This method can quantitative analysis go out under some design parameters, the corresponding Frequency Response of sensor；Also can be by it
TRANSFER MODEL, Selection and Design is optimized to the relevant parameter in sensor, plays finger in actual development process to sensor
Lead effect.
Brief description of the drawings
Fig. 1 is the flow chart of the angular-rate sensor modeling method of the invention based on MHD effect；
Fig. 2 is conductor fluid flowing equivalent model schematic diagram in sensor probe；
Fig. 3 is the schematic equivalent circuit of transformer.
Embodiment
Below in conjunction with accompanying drawing, by specific embodiment, the invention will be further described, and these embodiments are merely to illustrate
The present invention, is not limiting the scope of the invention.
As shown in figure 1, the invention provides a kind of angular-rate sensor modeling method based on MHD effect,
Comprise the steps of：
Step S1, sets up the TRANSFER MODEL of sensor probe；
Step S2, sets up the transmission function of primary amplifying transformer；
Step S3, sets up the transmission characteristic of rear end instrument amplifying circuit；
Step S4, according to the result obtained by step S1, S2 and S3, builds the block mold of sensor.
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, the step S1 is specific
Comprise the steps of：
Step S11, the electric-field intensity inside conductor fluid is obtained according to Ohm's law：
J=σ (E+V × B) (1)
J is current strength in formula, and E is electric-field intensity, and B is magnetic induction intensity, and σ is electrical conductivity；When under rectangular coordinate system
Magnetic field condition is Bx=0, By=-B0,Bz=0, then it can obtain z-axis directional current intensity：
U in formulaiFor lower plate speed, u is the movement velocity of conductor fluid infinitesimal, and r is O points and lower plate distance, such as Fig. 2 institutes
Show；
Step S12, calculates the electromagnetic force for acting on conductor fluid infinitesimal, Fe=JyBz-JzBy=B0Jz；
Step S13, according to Hartmann's constant physical significance, calculating acts on the viscous force of conductor fluid micelle；
I.e.FuFor the viscous force of micelle, M is Hartmann's constant,Wherein h is conductor fluid
Ring thickness, η, ρ, v are respectively the resistivity, density, kinematic viscosity coefficient of conductor fluid；
Step S14, electromagnetic force and viscous force according to suffered by kinematical equation and fluid infinitesimal, obtains conductor fluid micro-
The movement velocity of group：
ρ is the density of conducting liquid in formula；
Step S15, according to the law of electromagnetic induction and the movement velocity of fluid infinitesimal, obtains probe output voltage and angular speed
Transitive relation：
L is electrical conduction current body cavity height, r in formulaRMSFor the root mean square radii of conductor fluid ring.
Therefore the transmission function G of sensor probe1(s) it is
In formula：
B0：Perpendicular to the magnetic field intensity of conductor fluid annulated column, 0.24T is set to；
l：Wire cutting cd length (i.e. the height of ring-shaped chamber), is set to 16.6mm；
rRMS：The root mean square radii of conductor fluid ring, is set to 9.9mm；
h：The thickness of conductor fluid ring, is set to 1.4mm；
ν：Kinematic viscosity coefficientIt is set to 7.5 × 10-8m2/s；
M：The ratio of Harmann number, measurement magnetic force and viscous force, numerical value is set to 10.74.
Transmission function can must be popped one's head in for proportional component K by formula (7)1, differential s and an inertial elementSeries connection
Combination.In addition, magnetic field intensity B0, conductor fluid height l, the root mean square radii r of conductor fluid ringRMS, Harmann number M, electrical conduction current
The kinematic viscosity coefficient ν of body numerical value determines parameter K1, the as height of amplitude-versus-frequency curve；Harmann number M, conductor fluid
Thickness h, the kinematic viscosity coefficient ν of conductor fluid determine the handing-over frequency of amplitude-versus-frequency curve
The numerical value of parameter above is substituted into formula (7), the transmission characteristic that can obtain probe segment is：
Join frequency
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, in the step S2, structure
The equivalent circuit of transformer is built, as shown in figure 3, according to the transmission function of Kirchhoff's law to transformer：
L in formula1For primary inductance, L2For secondary inductance, L12For primary secondary coil mutual inductance, R1For primary line
The resistance of circle, R2For the resistance of secondary coil, C2For the parasitic capacitance of secondary coil.
Selection and Design can be optimized to transformer, secondary coil and core material according to transformer model.
The above-mentioned angular speed based on MHD effect passes speed modeling method, wherein, in the step S3,
The transmission characteristic of rear end instrument amplifying circuit is regard as constant value K in the frequency range of sensor responseU.Amplified by adjusting rear end instrument
Circuit gain KU, angular-rate sensor is arranged to proper ratio coefficient.
The above-mentioned angular-rate sensor modeling method based on MHD effect, wherein, will in the step S4
Sensor probe, primary amplifying transformer transmission function corresponding with amplifying circuit with rear end instrument combines, and is sensed
Device entirety TRANSFER MODEL：
In formula：
In summary, present system gives sensor probe model, transformer model, rear end amplifying circuit model and
The sensor entire process model that three is constituted.This method can quantitative analysis go out under some design parameters, the corresponding frequency response of sensor
Characteristic；Also Selection and Design can be optimized to the relevant parameter in sensor by its TRANSFER MODEL, sensor is actually being ground
Directive function is played during system.
Although present disclosure is discussed in detail by above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read the above, for the present invention's
A variety of modifications and substitutions all will be apparent.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (5)
1. a kind of angular-rate sensor modeling method based on MHD effect, it is characterised in that comprise the steps of：
Step S1, sets up the TRANSFER MODEL of sensor probe；
Step S2, sets up the transmission function of primary amplifying transformer；
Step S3, sets up the transmission characteristic of rear end instrument amplifying circuit；
Step S4, according to the result obtained by step S1, S2 and S3, builds the block mold of sensor.
2. the angular-rate sensor modeling method as claimed in claim 1 based on MHD effect, it is characterised in that
The step S1 is specifically comprised the steps of：
Step S11, the electric-field intensity inside conductor fluid is obtained according to Ohm's law：
J=σ (E+V × B) (1)
J is current strength in formula, and E is electric-field intensity, and B is magnetic induction intensity, and σ is electrical conductivity；Magnetic field under rectangular coordinate system
Condition is Bx=0, By=-B0,Bz=0, then it can obtain z-axis directional current intensity：
<mrow>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&sigma;B</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mi>y</mi>
<mi>r</mi>
</mfrac>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
U in formulaiFor lower plate speed, u is the movement velocity of conductor fluid infinitesimal, and r is O points and lower plate distance；
Step S12, calculates the electromagnetic force for acting on conductor fluid infinitesimal, Fe=JyBz-JzBy=B0Jz；
Step S13, according to Hartmann's constant physical significance, calculating acts on the viscous force of conductor fluid micelle；
I.e.FuFor the viscous force of micelle, M is Hartmann's constant,Wherein h is that conductor fluid ring is thick
Degree, η, ρ, v are respectively the resistivity, density, kinematic viscosity coefficient of conductor fluid；
Step S14, electromagnetic force and viscous force according to suffered by kinematical equation and fluid infinitesimal, obtains conductor fluid micelle
Movement velocity：
<mrow>
<mi>&rho;</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>u</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>F</mi>
<mo>=</mo>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
<mo>+</mo>
<msub>
<mi>F</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>B</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>M</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mi>y</mi>
<mi>r</mi>
</mfrac>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
ρ is the density of conducting liquid in formula；
Step S15, according to the law of electromagnetic induction and the movement velocity of fluid infinitesimal, obtains the biography of probe output voltage and angular speed
Pass relation：
<mrow>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Blr</mi>
<mrow>
<mi>R</mi>
<mi>M</mi>
<mi>S</mi>
</mrow>
</msub>
<mi>s</mi>
</mrow>
<mrow>
<mi>s</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>M</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mi>v</mi>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mfrac>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
L is electrical conduction current body cavity height, r in formulaRMSFor the root mean square radii of conductor fluid ring.
3. the angular-rate sensor modeling method as claimed in claim 1 based on MHD effect, it is characterised in that
In the step S2, the equivalent circuit of transformer is built, according to the transmission function of Kirchhoff's law to transformer：
<mrow>
<mfrac>
<mrow>
<msubsup>
<mi>V</mi>
<mi>O</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>V</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>L</mi>
<mi>s</mi>
</mrow>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>CR</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
L in formula1For primary inductance, L2For secondary inductance, L12For primary secondary coil mutual inductance, R1For primary coil
Resistance, R2For the resistance of secondary coil, C2For the parasitic capacitance of secondary coil.
4. the angular speed based on MHD effect passes speed modeling method as claimed in claim 1, it is characterised in that
In the step S3, the transmission characteristic of rear end instrument amplifying circuit is regard as constant value K in the frequency range that sensor is respondedU。
5. the angular-rate sensor modeling method as claimed in claim 1 based on MHD effect, it is characterised in that
In the step S4, by sensor probe, primary amplifying transformer and rear end instrument are integrated with the corresponding transmission function of amplifying circuit
Together, the overall TRANSFER MODEL of sensor is obtained：
<mrow>
<mfrac>
<mrow>
<msub>
<mi>V</mi>
<mi>O</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&omega;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>ks</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msub>
<mi>d</mi>
<mn>4</mn>
</msub>
<msup>
<mi>s</mi>
<mn>4</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>3</mn>
</msub>
<msup>
<mi>s</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula：
<mrow>
<mi>d</mi>
<mn>2</mn>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mrow>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>vC</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>M</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<mi>d</mi>
<mn>3</mn>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mfrac>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mfrac>
<mi>v</mi>
<msup>
<mi>h</mi>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>L</mi>
<mn>12</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>M</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<mi>d</mi>
<mn>4</mn>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>L</mi>
<mn>12</mn>
<mn>2</mn>
</msubsup>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710282612.8A CN107092756B (en) | 2017-04-26 | 2017-04-26 | Angular velocity sensor modeling method based on magnetohydrodynamics effect |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710282612.8A CN107092756B (en) | 2017-04-26 | 2017-04-26 | Angular velocity sensor modeling method based on magnetohydrodynamics effect |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107092756A true CN107092756A (en) | 2017-08-25 |
| CN107092756B CN107092756B (en) | 2020-10-02 |
Family
ID=59637135
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710282612.8A Active CN107092756B (en) | 2017-04-26 | 2017-04-26 | Angular velocity sensor modeling method based on magnetohydrodynamics effect |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107092756B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108549787A (en) * | 2018-07-13 | 2018-09-18 | 哈尔滨工业大学 | A kind of rocket large amplitude liquid sloshing method for establishing model based on movement pulsation ball |
| CN109190266A (en) * | 2018-09-10 | 2019-01-11 | 石家庄铁道大学 | A kind of simplified modeling method of the rigid multibody dynamics based on ADAMS software |
| CN109696181A (en) * | 2018-12-25 | 2019-04-30 | 上海航天控制技术研究所 | The equivalent detection circuit of MHD angular oscillation sensor and its modification method of frequency bandwidth characteristics |
| CN113569375A (en) * | 2021-04-26 | 2021-10-29 | 上海卫星工程研究所 | Non-contact magnetic suspension actuator transfer characteristic modeling and ground calibration method and system |
| CN113822354A (en) * | 2021-09-17 | 2021-12-21 | 合肥工业大学 | Micro-nano probe dynamic characteristic compensation method based on Bayesian inverse calculus modeling |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0444713A1 (en) * | 1990-03-02 | 1991-09-04 | Nippon Soken, Inc. | Device for detecting angular velocity |
| US6173611B1 (en) * | 1999-04-30 | 2001-01-16 | A-Tech Corporation | MHD sensor for measuring microradian angular rates and displacements |
| CN102878989A (en) * | 2012-09-21 | 2013-01-16 | 航天东方红卫星有限公司 | Triaxial angular vibration measuring method through adopting satellite-borne linear accelerometers |
| CN203537322U (en) * | 2013-09-18 | 2014-04-09 | 深圳市绿恩环保技术有限公司 | Weak photoelectric signal amplification circuit and processing board |
| CN103941034A (en) * | 2014-04-03 | 2014-07-23 | 天津大学 | Magnetofluid composite angular velocity sensor |
| CN105486315A (en) * | 2015-11-26 | 2016-04-13 | 航天东方红卫星有限公司 | Method for adjusting remote sensing satellite absolute calibration attitude relative to moon |
| CN105486882A (en) * | 2015-11-30 | 2016-04-13 | 北京遥测技术研究所 | Detection device for MHD angular-rate sensor |
-
2017
- 2017-04-26 CN CN201710282612.8A patent/CN107092756B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0444713A1 (en) * | 1990-03-02 | 1991-09-04 | Nippon Soken, Inc. | Device for detecting angular velocity |
| US6173611B1 (en) * | 1999-04-30 | 2001-01-16 | A-Tech Corporation | MHD sensor for measuring microradian angular rates and displacements |
| CN102878989A (en) * | 2012-09-21 | 2013-01-16 | 航天东方红卫星有限公司 | Triaxial angular vibration measuring method through adopting satellite-borne linear accelerometers |
| CN203537322U (en) * | 2013-09-18 | 2014-04-09 | 深圳市绿恩环保技术有限公司 | Weak photoelectric signal amplification circuit and processing board |
| CN103941034A (en) * | 2014-04-03 | 2014-07-23 | 天津大学 | Magnetofluid composite angular velocity sensor |
| CN105486315A (en) * | 2015-11-26 | 2016-04-13 | 航天东方红卫星有限公司 | Method for adjusting remote sensing satellite absolute calibration attitude relative to moon |
| CN105486882A (en) * | 2015-11-30 | 2016-04-13 | 北京遥测技术研究所 | Detection device for MHD angular-rate sensor |
Non-Patent Citations (8)
| Title |
|---|
| MENGJIE XU,ETC: "Error analysis of theoretical model of angular velocity sensor based on magnetohydrodynamics at low frequency", 《SENSORS AND ACTUATORS A: PHISICAL》 * |
| YUE JI,ETC: "Study on magnetohydrodynamics angular rate sensor under non-uniform magnetic field", 《SENSOR REVIEW》 * |
| YUE JI,ETC: "Theoretical and Experimental Study of Radial Velocity Generation for Extending Bandwidth of Magnetohydrodynamic Angular Rate Sensor at Low Frequency", 《SENSORS 2015》 * |
| 刘少谦: "磁流变橡胶减振器的设计与研究", 《全国优秀硕士论文全文数据库 工程科技Ⅱ辑》 * |
| 刘腾骏，等: "敏捷光学卫星多模式推扫成像时的偏流角研究", 《上海航天》 * |
| 张瀚, 等: "MHD角振动传感器宽带特性设计研究", 《第十三届全国敏感元件与传感器学术会议论文集》 * |
| 谭天乐，等: "航天器微振动测试、隔离、抑制技术综述", 《上海航天》 * |
| 霍红庆，等: "MHD角速度传感器磁场设计研究", 《真空与低温》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108549787A (en) * | 2018-07-13 | 2018-09-18 | 哈尔滨工业大学 | A kind of rocket large amplitude liquid sloshing method for establishing model based on movement pulsation ball |
| CN109190266A (en) * | 2018-09-10 | 2019-01-11 | 石家庄铁道大学 | A kind of simplified modeling method of the rigid multibody dynamics based on ADAMS software |
| CN109190266B (en) * | 2018-09-10 | 2022-12-06 | 石家庄铁道大学 | Multi-rigid-body dynamics simplified modeling method based on ADAMS software |
| CN109696181A (en) * | 2018-12-25 | 2019-04-30 | 上海航天控制技术研究所 | The equivalent detection circuit of MHD angular oscillation sensor and its modification method of frequency bandwidth characteristics |
| CN109696181B (en) * | 2018-12-25 | 2020-10-20 | 上海航天控制技术研究所 | Equivalent detection circuit of MHD angular vibration sensor and correction method of frequency band characteristic of equivalent detection circuit |
| CN113569375A (en) * | 2021-04-26 | 2021-10-29 | 上海卫星工程研究所 | Non-contact magnetic suspension actuator transfer characteristic modeling and ground calibration method and system |
| CN113822354A (en) * | 2021-09-17 | 2021-12-21 | 合肥工业大学 | Micro-nano probe dynamic characteristic compensation method based on Bayesian inverse calculus modeling |
| CN113822354B (en) * | 2021-09-17 | 2022-12-06 | 合肥工业大学 | Micro-nano probe dynamic characteristic compensation method based on Bayesian inverse calculus modeling |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107092756B (en) | 2020-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107092756A (en) | A kind of angular-rate sensor modeling method based on MHD effect | |
| CN107391818B (en) | A kind of Vibrating modal parameters recognition methods based on state observer | |
| CN105093138A (en) | Magnetic field detecting sensor and magnetic field detecting apparatus using the same | |
| WO2020253795A1 (en) | Modal localization effect-based weak magnetic field measurement device and method | |
| CN109186638B (en) | Accelerometer servo circuit with controllable current scale factor and manufacturing process thereof | |
| CN103235189B (en) | A kind of micro-resistance high-precision measuring method based on double-current voltage ratio method and realize the measuring system of the method | |
| CN102252746B (en) | Big-damping-ratio-based double parameter speed and acceleration output vibration pickup | |
| CN105628976A (en) | Method, processor, and system for performance parameter calibration of MEMS acceleration sensor | |
| CN106121637A (en) | A kind of system and method for detecting accident well | |
| CN108049865B (en) | Electric field calibration method for induction logging while drilling | |
| CN106840155A (en) | MHD angular-rate sensors and high accuracy gyroscope instrument combination measurement method | |
| CN105158509B (en) | A kind of three-dimensional flow change rate sensor and measuring method based on flexoelectric effect | |
| CN103323155B (en) | Container Special pull sensor and system | |
| CN105300269B (en) | A kind of wireless accurate strain gauge means and a kind of wireless accurate strain measurement method | |
| CN105965380B (en) | eddy current measuring device for wafer surface metal film polishing process | |
| CN102288815B (en) | Temperature compensator for giant magneto-resistance effect current sensor | |
| CN104156535B (en) | A kind of induction type acceleration transducer emulation mode based on circuit equivalent model | |
| Chattopadhyay et al. | Modification of the Maxwell–Wien bridge for accurate measurement of a process variable by an inductive transducer | |
| CN207991561U (en) | A kind of device adjusted for resistance strain gage output | |
| CN109708793A (en) | Stress test system and test method | |
| CN106643637B (en) | Inspection type principal strain measurement method for vibration measurement | |
| CN205725711U (en) | Highly integrated digital integrator for squirrel cage type hollow-core current transformer | |
| CN102645642B (en) | Temperature property testing device of amorphous metal fibers | |
| CN204854720U (en) | Angular transducer sensing head structure and processing circuit thereof | |
| CN202041541U (en) | Two-parameter velocity and acceleration output vibration pickup based on high damping ratio |
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 |