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 PDF

Info

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
Application number
CN201710282612.8A
Other languages
Chinese (zh)
Other versions
CN107092756B (en
Inventor
顾玥
乔洋
朱庆华
吴建铭
王坤东
陈桦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN201710282612.8A priority Critical patent/CN107092756B/en
Publication of CN107092756A publication Critical patent/CN107092756A/en
Application granted granted Critical
Publication of CN107092756B publication Critical patent/CN107092756B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design 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

A kind of angular-rate sensor modeling method based on MHD effect
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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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
CN201710282612.8A 2017-04-26 2017-04-26 Angular velocity sensor modeling method based on magnetohydrodynamics effect Active CN107092756B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (7)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
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
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
Löfdahl et al. Characteristics of a hot-wire microsensor for time-dependent wall shear stress measurements
CN105300269B (en) A kind of wireless accurate strain gauge means and a kind of wireless accurate strain measurement method
CN106840155A (en) MHD angular-rate sensors and high accuracy gyroscope instrument combination measurement method
CN102288815B (en) Temperature compensator for giant magneto-resistance effect current sensor
CN105965380B (en) Electric eddy current measurement device for wafer surface metallic film polishing process
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
CN104156535B (en) A kind of induction type acceleration transducer emulation mode based on circuit equivalent model
CN201096557Y (en) Sub micro Newton grade force measuring apparatus
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
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