CN101819056B - Instrument electronic device for checking and diagnosing flow meter and method - Google Patents

Abstract

The invention discloses an instrument electronic device for checking and diagnosing a flow meter and a method. According to the embodiment of the invention, the instrument electronic device (20) for the flow meter (5) is provided. The instrument electronic device (20) comprises an interface (201) for receiving a vibration response from the flow meter (5) and a processing system (203) communicated with the interface (201), wherein the vibration response is the response to the vibration of the flow meter (5) at a fundamental harmonic frequency. The processing system (203) is used for receiving the vibration response from the interface (201) so as to determine the frequency (omega0), the response voltage (V) and the driving current (I) of the vibration response, measure the attenuation characteristic (zeta) of the flow meter (5) and determine the rigidity parameter (K) according to the frequency (omega0), the response voltage (V), the driving current (I) and the attenuation characteristic (zeta).

Description

The instrument electronic device and the method that are used for the checking and diagnosing of flowmeter

The application is that application number is 200580051626.4, the applying date is on September 19th, 2005, denomination of invention is divided an application for the application of " instrument electronic device and the method that are used for the checking and diagnosing of flowmeter ".

Technical field

The present invention relates to a kind of instrument electronic device and method of the checking and diagnosing for flowmeter.

Background technology

The statement of problem

The vibrating conduit sensor, such as Coriolis (Coriolis) mass flowmeter or vibrating tube densimeter, the motion that generally comprises the vibrating conduit of fluent material by detection comes work.Can determine the characteristic relevant with ducted material by processing from the measuring-signal of the motion converter reception relevant with pipeline, such as mass rate, density etc.By comprising pipeline and the combination quality that is contained in material wherein, rigidity and damping characteristic affect the vibration mode of vibration Material Filling system usually.

The pipeline of vibratory flowmeter can comprise one or more stream pipes.The stream pipe is forced to the vibration at the resonance frequency place, and in this case, the resonance frequency of pipe flows the density of fluid in the pipe in proportion to.Coupled vibration between the end of the sensor measuring tube that the entrance and exit of pipe part is positioned.In flow process, because Coriolis force, vibrating tube and flow mass are coupled, and cause the phase shift of vibrating between the end of pipe.This phase shift is mass flow in proportion to directly.

Typical coriolis mass flowmeters comprises to be pipeline or other transmission system and to transmit material, the fluid in the system for example, one or more pipelines of the continuous interior lines in the mud etc.Each pipeline can be regarded as having one group of eigentone, for example comprises simple crookedly, reverses, radially and coupled mode.In typical Coriolis mass flow measurement application, when the flow of material piping, with one or more vibration modes excitation pipeline configurations, and along pipeline enclosure every the moving of some place's measuring channel.Typically provide excitation by actuator, electro-mechanical devices for example, such as the sound coil type of drivers, it disturbs pipeline with periodic form.By the time delay between the measurement translator position motion or differ and to determine mass rate.Typically adopt two this sensors (or pickoff sensor (pick-off sensor)), thereby measure the vibratory response of stream pipeline or pipeline, and typically be positioned in upstream position and the downstream position of actuator.By cable, two pickoff sensors are connected to electronic equipment.This equipment receives the signal from two pickoff sensors, and processing signals, in order to obtain mass flow measurement.

Differing between two sensor signals relates to the mass rate of the material that flows through a stream pipe or a plurality of stream pipes.The mass rate of material is the time delay between two sensor signals in proportion to, and therefore by making time delay and traffic alignment factor (FCF, Flow CalibrationFactor) multiply each other and can determine mass rate, in this case, time delay comprises and differing divided by frequency.Material behavior and the cross section property of FCF reflection stream pipe.In the prior art, before the installation flowmeter enters pipeline or other pipeline, determine FCF by calibration steps.In calibration steps, fluid is with the logical flow tube of given flow, and calculate differ and flow between characteristic.

An advantage of coriolis flowmeter is that the degree of accuracy of the mass rate of measurement is not affected by the wearing and tearing of the movable part of flowmeter.Multiply by the traffic alignment factor by differing between two points that make the stream pipe and determine flow.Only input is the sinusoidal signal from sensor, upper 2 vibration of expression stream pipe.Differ from these sinusoidal signals calculating.In vibrating flow tube, there is not movable part.Therefore, the measurement that differs with the traffic alignment factor is not affected by the wearing and tearing of movable part in the flow.

FCF can relate to the stiffness characteristics of flowermeter.If the stiffness characteristics of flowermeter changes, then FCF also changes.Change therefore impact by the degree of accuracy of the flow measurement of flowmeter generation.For example by corrosion or the variation of corroding the cross section property of the variation can cause material and stream pipe.Therefore, high expectations can detect and/or quantize any variation for the rigidity of flowermeter, thereby keeps the high-caliber degree of accuracy of flowmeter.

Summary of the invention

According to embodiments of the invention, be provided for the instrument electronic device of flowmeter.This instrument electronic device comprise for receive from the interface of the vibratory response of flowmeter and with the disposal system of interface communication.This vibratory response is included in the fundamental resonance frequency place to the response of the vibration of flowmeter.This disposal system is arranged to receive the vibratory response from interface, determines the frequency (ω of vibratory response ₀), determine response voltage (V) and the drive current (I) of vibratory response, the attenuation characteristic of measuring flow meter (ζ), and from frequency (ω ₀), response voltage (V), drive current (I) and attenuation characteristic (ζ) are determined stiffness parameters (K).

According to embodiments of the invention, provide a kind of method of the stiffness parameters (K) for determining flowmeter.The method comprises the vibratory response that receives from flowmeter.This vibratory response is included in the fundamental resonance frequency place to the response of the vibration of flowmeter.The method further comprises the frequency (ω that determines vibratory response ₀), determine response voltage (V) and the drive current (I) of vibratory response and the attenuation characteristic of measuring flow meter (ζ).The method further comprises from frequency (ω ₀), response voltage (V), drive current (I) and attenuation characteristic (ζ) are determined stiffness parameters (K).

According to embodiments of the invention, provide a kind of method of the stiffness variation (Δ K) for determining flowmeter.The method comprises the vibratory response that receives from flowmeter.This vibratory response is included in the fundamental resonance frequency place to the response of the vibration of flowmeter.The method further comprises the frequency (ω that determines vibratory response ₀), determine response voltage (V) and the drive current (I) of vibratory response and the attenuation characteristic of measuring flow meter (ζ).The method further comprises from frequency (ω ₀), response voltage (V), drive current (I) and attenuation characteristic (ζ) are determined stiffness parameters (K).The method further is included in the second time t ₂The place receives the second vibratory response from flowmeter, produces the second stiffness characteristics (K from the second vibratory response ₂), compare the second stiffness characteristics (K ₂) and stiffness parameters (K), and if the second stiffness characteristics (K ₂) different from predetermined tolerance with stiffness parameters (K), then detect stiffness variation (Δ K).

According to embodiments of the invention, be provided for the instrument electronic device of flowmeter.This instrument electronic device comprises for the interface that receives from three or more vibratory responses of flowmeter.These three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response.This instrument electronic device further comprises the disposal system with interface communication, and this disposal system is arranged to receive three or more vibratory responses from interface, produce limit-residue frequency response function from these three or more vibratory responses, and determine at least one stiffness parameters (K) from limit-residue frequency response function.

According to embodiments of the invention, provide a kind of method of the stiffness variation (Δ K) for determining flowmeter.The method comprises three or more vibratory responses of reception, and these three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response.The method further comprises from these three or more vibratory responses generation limit-residue frequency response functions, and determines at least one stiffness parameters (K) from limit-residue frequency response function.

According to embodiments of the invention, provide a kind of method of the stiffness parameters (K) for determining flowmeter.The method comprises three or more vibratory responses of reception, and these three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response.The method further comprises from these three or more vibratory responses generation limit-residue frequency response functions and from limit-residue frequency response function determines at least one stiffness parameters (K).The method further is included in the second time t ₂The place receives second group of three or more vibratory response from flowmeter, produces the second stiffness characteristics (K from these second group of three or more vibratory response ₂), compare the second stiffness characteristics (K ₂) and stiffness parameters (K), and if the second stiffness characteristics (K ₂) different from predetermined tolerance from stiffness parameters (K), then detect stiffness variation (Δ K).

Various aspects of the present invention

In aspect of instrument electronic device, measure attenuation characteristic (ζ) and comprise that further the vibratory response of permissible flow meter decays to the predetermined vibration target downwards.

Instrument electronic device on the other hand in, disposal system further is arranged through the excitation of removing flowmeter, and the vibratory response of permissible flow meter decays to the predetermined vibration target downwards and measures attenuation characteristic (ζ) when measuring attenuation characteristic.

Instrument electronic device on the other hand in, stiffness parameters (K) comprises K=(I*BL _PO* BL _DR* ω ₀)/2 ζ V.

In aspect of the method, measure attenuation characteristic (ζ) and comprise that further the vibratory response of permissible flow meter decays to the predetermined vibration target downwards.

The method on the other hand in, measures attenuation characteristic (ζ) and further comprise the excitation of removing flowmeter, and the vibratory response of permissible flow meter decays to the predetermined vibration target downwards in the measurement attenuation characteristic.

The method on the other hand in, stiffness parameters (K) comprises K=(I*BL _PO* BL _DR* ω ₀)/2 ζ V.

The method on the other hand in, produce the second stiffness characteristics (K from the second vibratory response ₂) comprise from second frequency, the second response voltage, the second drive current and the second damping characteristic produce the second stiffness characteristics (K ₂).

The method on the other hand in, the method further comprises if the second stiffness parameters (K ₂) different from the predetermined stiffness tolerance with stiffness parameters (K), then detect stiffness variation (Δ K).

The method on the other hand in, the method further comprises from K ₂With K relatively quantize stiffness variation (Δ K).

In an embodiment of instrument electronic device, disposal system further is arranged to determines damping parameter (C) from limit-residue frequency response function.

In another embodiment of instrument electronic device, disposal system further is arranged to determines mass parameter (M) from limit-residue frequency response function.

In another embodiment of instrument electronic device, disposal system further is arranged to calculates limit (λ), left residue (R from limit-residue frequency response function _L) and right residue (R _R).

In another embodiment of instrument electronic device, these three or more vibratory responses comprise at least one tone (tone) that is higher than fundamental frequency response and at least one tone that is lower than fundamental frequency response.

In another embodiment of instrument electronic device, these three or more vibratory responses comprise at least two tones that are higher than fundamental frequency response and at least two tones that are lower than fundamental frequency response.

In another embodiment of instrument electronic device, limit-residue frequency response function comprises first order pole-residue frequency response function.

In another embodiment of instrument electronic device, limit-residue frequency response function comprises first order pole-residue frequency response function, and this first order pole-residue frequency response function comprises

$H\left(\mathrm{\ω}\right)=R/(\mathrm{j\ω}-\mathrm{\λ})+\stackrel{\‾}{R}/(\mathrm{j\ω}-\stackrel{\‾}{\mathrm{\λ}}).$

In another embodiment of instrument electronic device, limit-residue frequency response function comprises first order pole-residue frequency response function, and this first order pole-residue frequency response function comprises And wherein according to equation M=1/2jR ω _d, K=(ω _n) ²M and C=2 ζ ω _nM determines stiffness parameters (K), damping parameter (C) and mass parameter (M).

In another embodiment of instrument electronic device, limit-residue frequency response function comprises second order limit-residue frequency response function.

In another embodiment of instrument electronic device, limit-residue frequency response function comprises second order limit-residue frequency response function, and this second order limit-residue frequency response function comprises

$\stackrel{.}{H}\left(\mathrm{\ω}\right)=\frac{\stackrel{.}{X}\left(\mathrm{\ω}\right)}{F\left(\mathrm{\ω}\right)}=\frac{\mathrm{j\ω}}{-M{\mathrm{\ω}}^2+\mathrm{jC\ω}+K}.$

In another embodiment of instrument electronic device, limit-residue frequency response function comprises second order limit-residue frequency response function, and this second order limit-residue frequency response function comprises

And foundation wherein

Determine stiffness parameters (K), according to M=K/ (ω _n) ²Determine mass parameter (M), and foundation

Determine damping parameter (C).

In an embodiment of the method, described determine to comprise from limit-residue frequency response function further determine damping parameter (C).

In another embodiment of the method, described determine to comprise from limit-residue frequency response function further determine mass parameter (M).

In another embodiment of the method, calculate limit (λ), left residue (R described definite further comprising from limit-residue frequency response function _L) and right residue (R _R).

In another embodiment of the method, these three or more vibratory responses comprise at least one tone that is higher than fundamental frequency response and are lower than at least one tone of fundamental frequency response.

In another embodiment of the method, these three or more vibratory responses comprise at least two tones that are higher than fundamental frequency response and at least two tones that are lower than fundamental frequency response.

In another embodiment of the method, limit-residue frequency response function comprises first order pole-residue frequency response function.

In another embodiment of the method, limit-residue frequency response function comprises first order pole-residue frequency response function, and this first order pole-residue frequency response function comprises

$H\left(\mathrm{\ω}\right)=R/(\mathrm{j\ω}-\mathrm{\λ})+\stackrel{\‾}{R}/(\mathrm{j\ω}-\stackrel{\‾}{\mathrm{\λ}}).$

In another embodiment of the method, limit-residue frequency response function comprises first order pole-residue frequency response function, and this first order pole-residue frequency response function comprises

And wherein according to equation M=1/2jR ω _d, K=(ω _n) ²M and C=2 ζ ω _nM determines stiffness parameters (K), damping parameter (C) and mass parameter (M).

In another embodiment of the method, limit-residue frequency response function comprises second order limit-residue frequency response function.

In another embodiment of the method, limit-residue frequency response function comprises second order limit-residue frequency response function, and this second order limit-residue frequency response function comprises

$\stackrel{.}{H}\left(\mathrm{\ω}\right)=\frac{\stackrel{.}{X}\left(\mathrm{\ω}\right)}{F\left(\mathrm{\ω}\right)}=\frac{\mathrm{j\ω}}{-M{\mathrm{\ω}}^2+\mathrm{jC\ω}+K}.$

In another embodiment of the method, limit-residue frequency response function comprises second order limit-residue frequency response function, and this second order limit-residue frequency response function comprises And foundation wherein

Determine stiffness parameters (K), according to M=K/ (ω _n) ²Determine mass parameter (M), and foundation

Determine damping parameter (C).

In another embodiment of the method, the method further comprises if the second stiffness characteristics (K ₂) different from the predetermined stiffness tolerance from stiffness parameters (K), then detect stiffness variation (Δ K).

In another embodiment of the method, the method further comprises from K and K ₂relatively quantize stiffness variation (Δ K).

Description of drawings

Identical reference number represents element identical on whole accompanying drawings.

Fig. 1 shows the flowmeter that comprises instrument assembly and instrument electronic device;

Fig. 2 shows the instrument electronic device according to embodiments of the invention;

Fig. 3 is according to the process flow diagram of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter;

Fig. 4 is according to the process flow diagram of embodiments of the invention for the method for the stiffness variation (Δ K) of determining flowmeter

Fig. 5 shows the instrument electronic device according to another embodiment of the present invention;

Fig. 6 is according to the process flow diagram of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter;

Fig. 7 shows the embodiment that limit (λ) according to embodiments of the invention and residue (R) are found the solution;

Fig. 8 is the M that illustrates according to embodiments of the invention, the calcspar of the calculating of C and K systematic parameter;

Fig. 9 shows the rigidity estimating system based on FRF according to the integral body of embodiments of the invention;

Figure 10 is according to the process flow diagram of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter;

Figure 11 shows according to embodiments of the invention from the M of equation (29) to second order limit-residue response, the embodiment that C and K find the solution;

Figure 12 shows the rigidity estimating system based on FRF according to the integral body of embodiments of the invention.

Embodiment

Fig. 1-12 and following declarative description concrete example so that how instruction those skilled in the art obtain and utilize optimal mode of the present invention.In order to instruct inventive principle, simplified and omitted some traditional aspects.It will be appreciated by those skilled in the art that the distortion from these examples that falls into scope of the present invention.It will be appreciated by those skilled in the art that to be combined in many ways feature described below, in order to form a plurality of distortion of the present invention.As a result, the present invention is not limited to object lesson described below, but only limits by claim and their equivalent.

Fig. 1 shows the coriolis flowmeter 5 that comprises instrument assembly 10 and instrument electronic device 20.Mass rate and the density of instrument assembly 10 response rapidoprints.Via lead-in wire 100, instrument electronic device 20 is connected to instrument assembly 10, in order to provide density, mass rate and temperature information by path 26, and out of Memory not related to the present invention.Described a kind of coriolis flowmeter structure, yet it will be obvious to those skilled in the art that the vibrating tube densimeter of the additional measurement capability that provides by coriolis mass flowmeters can be provided not in the present invention.

Instrument assembly 10 comprise a pair of manifold 150 and 150 ', have flange neck 110 and 110 ' flange (flange) 103 and 103 ', the stream pipe 130 and 130 of pair of parallel ', driving mechanism 180, temperature sensor 190, and a pair of speed pickup 170L and 170R.Two basic straight entrance legs (leg) 131 and 131 of stream pipe 130 and 130 ' have ' and outlet leg 134 and 134 ', it is at stream pipe assembling block 120 and 120 ' locate to restrain toward each other.Stream pipe 130 and 130 ' crooked along their two symmetric position places of length and to run through their whole length substantially parallel.Brace 140 and 140 ' be used for limiting axle W and W ', each stream circumference of cannon bone is around this shaft vibration.

Stream pipe 130 and 130 ' side leg 131,131 ' and 134,134 ' regularly be attached to stream pipe assembling block 120 and 120 ', and these pieces be attached to regularly again manifold 150 and 150 '.This provides the continuous closed material path by Coriolis instrument assembly 10.

When connection has hole 102,102 ' flange 103 and 103 ' time, via inlet end 104 and endpiece 104 ', enter the production line (not shown), the measured rapidoprint of this production line delivery, by the hole 101 in the flange 103, material enters the end 104 of flowmeter, is directed to the stream pipe assembling block 120 with surface 121 by manifold 150.In manifold 150 inside, logical flow tube 130 and 130 ', material is separated and route (route).Based on current stream pipe 130 and 130 ', rapidoprint manifold 150 ' in reconfigured with independent stream, and be routed to thereafter by have bolt hole 102 ' flange 103 ' be connected to production line (not shown) endpiece 104 '.

Stream pipe 130 and 130 ' selected and be assembled to suitably stream pipe assembling block 120 and 120 ', thereby have respectively essentially identical mass distribution, around moment of inertia and the Young modulus of bending axis W-W and W '-W '.These bending axis by brace 140 and 140 '.Because the Young modulus of stream pipe is along with temperature change, and should change calculating affect flow and density, thus resistance temperature detector (RTD) 190 be assembled to stream pipe 130 ', in order to measure continuously the temperature that flows pipe.The temperature domination that the temperature of stream pipe and the voltage that manifests for the given electric current leap RTD by the there are thus led to the material of flow tube.Cross over that the voltage that depends on temperature that RTD manifests is used for compensation by instrument electronic device 20 in known method because the stream pipe 130 and 130 that any variation of stream pipe temperature causes ' the variation of elastic modulus.By going between 195, this RTD is connected to instrument electronic device 20.

Around they each bending axis W-W and the relative direction of W '-W ' on and be called under the first out-phase beam mode of flowmeter, by driver 180 drive two stream pipes 130 and 130 '.This driving mechanism 180 can comprise any one of a plurality of known layouts, such as be assembled to stream pipe 130 ' magnet, and the relative coil that is assembled to stream pipe 130, and in order to vibrate two stream pipes, alternating current is by described layout.Via lead-in wire 185, apply suitable driving signal to driving mechanism 180 by instrument electronic device 20.

Instrument electronic device 20 receives the RTD temperature signals at lead-in wire 195, and left and right rate signal manifests at go between 165L and 165R respectively.Instrument electronic device 20 is created in the driving signal that manifests on the lead-in wire 185, so that driving element 180 and vibrating tube 130 and 130 '.Instrument electronic device 20 is processed left and right rate signal and RTD signal, in order to calculate mass rate and density by instrument assembly 10.On path 26, by instrument electronic device 20 this information is applied to use device 29 with out of Memory.

Fig. 2 shows the instrument electronic device 20 according to embodiments of the invention.Instrument electronic device 20 can comprise interface 201 and disposal system 203.Instrument electronic device 20 receives such as the vibratory response 210 from instrument assembly 10.Instrument electronic device 20 is processed vibratory response 210, thereby the flow characteristics of the stream material of instrument assembly 10 is flow through in acquisition.In addition, in foundation instrument electronic device 20 of the present invention, vibratory response 210 is also processed, thereby determines the stiffness parameters (K) of instrument assembly 10.In addition, instrument electronic device 20 can be along with two or more this vibratory responses of change process of time, thus the stiffness variation in the measuring instrument assembly 10 (Δ K).Can carry out rigidity under mobile or non-current condition determines.The non-current advantage of determining in the vibratory response that obtains, to provide the noise level that reduces.

As previously discussed, the cross section property of traffic alignment factor (FCF) reflection material behavior and stream pipe.By the multiply each other mass rate of the stream material of determining to flow through flowmeter of the time delay (or differing/frequency) that will measure and FCF.This FCF can relate to the stiffness characteristics of instrument assembly.If the stiffness characteristics of instrument assembly changes, FCF also changes so.Therefore the change of the rigidity of flowmeter will affect the degree of accuracy of the flow measurement that produces by flowmeter.

The present invention is significant, determines because it allows instrument electronic device 20 to carry out rigidity in the field, and does not carry out actual stream calibration testing.It allows not have the rigidity of calibration test platform or other concrete equipment or concrete fluid to determine.This expects, is expensive, difficulty and time consuming because carry out the stream calibration in the field.Yet better and easier calibration inspection is expected, because the rigidity of instrument assembly 10 can change along with the time.This change may be because such as the corrosion of stream pipe, the erosion of stream pipe and the deadline factor of instrument assembly 10 is caused.

Utilize mathematical model can explain the present invention.By open loop, the second order driving model can represent the vibratory response of flowmeter, comprising:

$M\stackrel{.}{x}+C\stackrel{.}{x}+\mathrm{Kx}=f---\left(1\right)$

Wherein, f is the power that is applied to this system, and M is the quality of system, and C is damping characteristic, and K is the stiffness characteristics of system.Item K comprises K=M (ω ₀) ², and a C comprises C=M2 ζ ω ₀, wherein ζ comprises lag characteristic, and ω ₀=2 π f ₀, f wherein ₀Be take Hz as unit instrument assembly 10 intrinsic/resonance frequency.In addition, x is the physical displacement distance of vibration,

Be the speed of stream pipe displacement, and x is acceleration.This is commonly referred to as the MCK model.This formula can be rearranged into following form:

$M[{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">s</figure-callout>}}^2+2\mathrm{\&zeta;}{\mathrm{\&omega;}}_{0}s+{\mathrm{\&omega;}}_0^2]x=f---\left(2\right)$

Equation (2) can be further processed into the transport function form.With the transport function form, make firmly top offset item, comprising:

$\frac{x}{f}=\frac{s}{M[{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">s</figure-callout>}}^2+2\mathrm{\&zeta;}{\mathrm{\&omega;}}_{0}s+{\mathrm{\&omega;}}_0^2]}---\left(3\right)$

Known magnetic equation can be used for reduced equation (3).Two equations applicatory are:

$V={\mathrm{BL}}_{\mathrm{PO}}*\stackrel{.}{x}---\left(4\right)$

And

f＝BL _DR*I (5)

The sensor voltage VEMF of equation (4) (at pickoff sensor 170L or 170R place) equals pick-up sensitivity factor BL _POMultiply by the pickup velocity of motion

For each pickoff sensor, common known or measurement pick-up sensitivity factor BL _POThe power (f) that driver 180 by equation (5) produces equals driver sensitivity factor BL _DRMultiply by the drive current (I) that is supplied to driver 180.Usually known or measure the driver sensitivity factor BL of driver 180 _DRFactor BL _POAnd BL _DRAll be the function of temperature, and can be proofreaied and correct by temperature survey.

By the transport function with magnetic equation (4) and (5) substitution equation (3), the result is:

$\frac{V}{I}=\frac{{\mathrm{BL}}_{\mathrm{PO}}*{\mathrm{BL}}_{\mathrm{DR}}*s}{M[{\mathrm{<figure-callout\; id="2"\; label="s"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">s</figure-callout>}}^2+2\mathrm{\&zeta;}{\mathrm{\&omega;}}_{0}s+{\mathrm{\&omega;}}_0^2]}---\left(6\right)$

If instrument assembly 10 is driven into open loop at resonance, namely at resonance/natural frequency ω ₀(the ω wherein of place ₀=2 π f ₀), then equation (6) can be write as again:

${\left(\frac{V}{I}\right)}_{{\mathrm{\&omega;}}_0}=\frac{{\mathrm{BL}}_{\mathrm{PO}}*{\mathrm{BL}}_{\mathrm{DR}}*{\mathrm{\&omega;}}_0}{2\mathrm{\&zeta;}\left[M{\mathrm{\&omega;}}_0^2\right]}---\left(7\right)$

By substituting rigidity, equation (7) is simplified as:

${\left(\frac{V}{I}\right)}_{{\mathrm{\&omega;}}_0}=\frac{{\mathrm{BL}}_{\mathrm{PO}}*{\mathrm{BL}}_{\mathrm{DR}}*{\mathrm{\&omega;}}_0}{2\mathrm{\&zeta;K}}---\left(8\right)$

Here, stiffness parameters (K) can be separated, thereby obtain:

$K=\frac{I*{\mathrm{BL}}_{\mathrm{PO}}*{\mathrm{BL}}_{\mathrm{DR}}*{\mathrm{\&omega;}}_0}{2\mathrm{\&zeta;V}}---\left(9\right)$

As a result, by measuring/quantification attenuation characteristic (ζ) and driving voltage (V) and drive current (I), can determine stiffness parameters (K).Can be determined from the response voltage that picks up (V) by vibratory response and drive current (I).Discuss below the step of definite stiffness parameters (K) in more detail in conjunction with Fig. 3.

In use, can follow the tracks of in time stiffness parameters (K).For example, statistical technique can be used for definite any variation (namely stiffness parameters (K)) in time.The statistics variations of stiffness parameters (K) can represent that FCT changes for concrete flowmeter.

The invention provides a kind of stiffness parameters (K) that does not rely on the calibration density value that is stored or recovers.This and prior art form contrast, and wherein known stream material is used for the factory calibrated operation, in order to obtain can be used for all the in the future density criterion of calibration operation.The invention provides a kind of only from the independent stiffness parameters (K) that obtains of the vibratory response of flowmeter.The invention provides a kind of rigidity detection/calibration steps that does not need the factory calibrated step.

Via the lead-in wire 100 of Fig. 1, the reception vibratory response 210 of interface 201 from speed pickup 170L and 170R.Interface 201 can carry out any must or the signal conditioning of expectation, such as the format of any mode, amplification, buffering etc.Alternately, some or all that can the executive signal condition in disposal system 203.In addition, interface 201 can allow the communication between instrument electronic device 20 and the external device (ED).Interface 201 can carry out any in electricity, light or the radio communication.

Interface 201 is coupled with the Aristogrid (not shown) in one embodiment, and wherein sensor signal comprises analog sensor signal.This Aristogrid sampling and digitized simulation vibratory response, and produce digital vibratory response 210.

The operation of disposal system 203 management instrument electron devices 20, and processing is from the flow measurement of instrument assembly 10.Disposal system 203 is carried out one or more handling procedures, and therefore processes flow measurement, thereby produces one or more properties of flow.

Disposal system 203 can comprise multi-purpose computer, microprocessing systems, logical circuit or some other treating apparatus general or customization.Can be in a plurality for the treatment of apparatus distributed processing system (DPS) 203.Disposal system 203 can comprise the monoblock type of any mode or electronic storage medium independently, such as storage system 204.

Storage system 204 can be stored flow meter parameter and data, software program, constant and variable.In one embodiment, storage system 204 comprises the program of carrying out by disposal system 203, such as the rigidity program 230 of the stiffness parameters (K) of determining flowmeter 5.

Rigidity program 230 can be arranged 203 one-tenth vibratory responses that receive from flowmeter of disposal system in one embodiment, and this vibratory response is included in the fundamental resonance frequency place to the response of the vibration of flowmeter, determines the frequency (ω of vibratory response ₀), determine response voltage (V) and the drive current (I) of vibratory response, the attenuation characteristic of measuring flow meter (ζ), and from frequency (ω ₀), response voltage (V), drive current (I) and attenuation characteristic (ζ) are determined stiffness parameters (K) (seeing Fig. 3 and relevant discussion).

Rigidity program 230 can arrange that 203 one-tenth of disposal systems receive vibratory response in one embodiment, determines frequency, determines response voltage (V) and drive current (I), measures attenuation characteristic (ζ) and definite stiffness parameters (K).Rigidity program 230 arranges that further 203 one-tenth of disposal systems are at the second time t in this embodiment ₂The place receives the vibratory response from flowmeter, the second vibratory response is repeated to determine and measuring process, thereby produce the second stiffness characteristics (K ₂), compare the second stiffness characteristics (K ₂) and stiffness parameters (K), and if the second stiffness characteristics (K ₂) different from tolerance 224 from stiffness parameters (K), then detect stiffness variation (Δ K) (seeing Fig. 4 and relevant discussion).

In one embodiment, storage system 204 storages are used for the variable of operate flow meter 5.Storage system 204 storage of variables in one embodiment, such as vibratory response 210, for example it can receive from speed/pickoff sensor 170L and 170R.

In one embodiment, storage system 204 storage of constant, coefficient and operating variable.For example, storage system 204 the second stiffness characteristics 221 that can store definite stiffness characteristics 220 and produce at subsequently time point place.Storage system 204 can be stored working value, such as the frequency 212 of vibratory response 210, and the voltage 213 of vibratory response 210 and the drive current 214 of vibratory response 210.Storage system 204 can further be stored the Vibration Targets 226 of flowmeter 5 and the attenuation characteristic 215 of measuring.In addition, storage system 204 can storage of constant, and threshold value or scope are such as tolerance 224.In addition, storage system 204 can be stored the in time data of cycle accumulation, such as stiffness variation 228.

Fig. 3 is according to the process flow diagram 300 of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter.In step 301, receive vibratory response from flowmeter.This vibratory response is that flowmeter is to the response of the vibration at fundamental resonance frequency place.This vibration can be continuous or intermittence.The stream material can flow through instrument assembly 10, perhaps can be stable.

In step 302, determine the frequency of vibratory response.By any method, step or hardware can be determined frequencies omega from vibratory response ₀

In step 303, determine voltage (V or the V of vibratory response _EMF), and drive current (I).Can obtain this voltage and drive current from vibratory response not processed or that be conditioned.

In step 304, the damping characteristic of measuring flow meter.Vibratory response by the permissible flow meter decays to Vibration Targets downwards and measures damping characteristic, measures simultaneously attenuation characteristic.Can carry out in several ways this attenuation.Can reduce drive signal amplitude, in fact driver 180 can carry out the braking (in suitable flowmeter) of instrument assembly 10, and perhaps driver 180 can be by hoisting power only until till reaching target.In one embodiment, Vibration Targets comprises the level that reduces in the driving set-point.For example, if it is current at the 3.4mV/Hz place to drive the set-point, so for damping measurement, drives the set-point and can be reduced to lower value, such as 2.5mV/Hz.By this way, instrument electronic device 20 can allow instrument assembly 10 slide simply until vibratory response is mated till this new driving target basically.

In step 305, from frequency, voltage, drive current and attenuation characteristic (ζ) are determined stiffness parameters (K).Can determine stiffness parameters (K) according to top equation (9).Except determining and follow the tracks of the rigidity (K) that the method is also determined and tracking damping parameter (C) and mass parameter (M).

Can iteration, periodically or at random carry out the method 300.Can locate to carry out the method in predetermined boundary mark (landmark), such as locating at predetermined hour that operates, when the stream changes in material, etc.

Fig. 4 is according to the process flow diagram 400 of embodiments of the invention for the method for the stiffness variation (Δ K) of determining flowmeter.In step 401, receive vibratory response from flowmeter, as previously discussed.

In step 402, determine the frequency of vibratory response, as previously discussed.

In step 403, determine voltage and the drive current of vibratory response, as previously discussed.

In step 404, the attenuation characteristic of measuring flow meter (ζ), as previously discussed.

In step 405, from frequency, voltage, drive current and attenuation characteristic (ζ) are determined stiffness parameters (K), as previously discussed.

In step 406, at the second time t ₂The place receives the second vibratory response.At time t ₂The place is from generation of vibration second vibratory response of instrument assembly 10.

In step 407, produce the second stiffness characteristics K from the second vibratory response ₂For example, utilize step 401 can produce the second stiffness characteristics K to 405 ₂

In step 408, the second stiffness characteristics K ₂(K) compares with stiffness parameters.This comparison that relatively is included in the stiffness characteristics that the different time place obtains is in order to detect stiffness variation (Δ K).

In step 409, determine K ₂And any stiffness variation between the K (Δ K).This stiffness variation determines to adopt statistics or the mathematical method for any mode of the marked change of determining rigidity.This stiffness variation (Δ K) can be stored, to be used in the future purposes and/or to be transferred into remote location.In addition, this stiffness variation (Δ K) can trigger the alert if in the instrument electronic device 20.Stiffness variation in one embodiment (Δ K) is at first made comparisons with tolerance 224.If stiffness variation (Δ K) surpasses tolerance 224, then determine error condition.Except determining and follow the tracks of the rigidity (K) that the method is also determined and tracking damping parameter (C) and mass parameter (M).

Can iteration, periodically or at random carry out the method 400.Can carry out the method at predetermined boundary mark place, such as locating at predetermined hour that operates, when the stream changes in material, etc.

Fig. 5 shows the instrument electronic device 20 according to another embodiment of the present invention.Instrument electronic device 20 can comprise interface 201 in this embodiment, disposal system 203 and storage system 204, as previously discussed.Instrument electronic device 20 is such as receiving three or more vibratory responses 505 from instrument assembly 10.This instrument electronic device 20 is processed these three or more vibratory responses 505, thereby the properties of flow of the stream material of instrument assembly 10 is flow through in acquisition.In addition, these three or more vibratory responses 505 also can be processed to determine the stiffness parameters (K) of instrument assembly 10.Instrument electronic device 20 can further be determined damping parameter (C) and mass parameter (M) from these three or more vibratory responses 505.These instrument component parameters can be for detection of the variation of instrument assembly 10, as previously discussed.

Storage system 204 can storage processing program, such as rigidity program 506.Storage system 204 can the storing received data, such as vibratory response 505.Storage system 204 can be stored value in advance programming or user's input, such as rigidity tolerance 516, and damping tolerance 517 and quality tolerance 518.Storage system 204 can be stored working value, such as limit (pole) (λ) 508 and residue (residue) (R) 509.Storage system 204 can be stored definite final value, such as rigidity (K) 510, and damping (C) 511 and quality (M) 512.This storage system 204 can be stored the fiducial value of in time cycle generation and operation, such as the second rigidity (K ₂) 521, the second quality (M ₂) 522, stiffness variation (Δ K) 530, damping change (Δ C) 531, and mass change (Δ M) 532.Stiffness variation (Δ K) 530 can comprise the variation such as the stiffness parameters (K) of the instrument assembly 10 of measuring in time.Stiffness variation (Δ K) 530 can be used to detect and determine the physical change of instrument assembly 10 in time, such as corrosion and etching effect.In addition, can measure in time the mass parameter (M) 512 with trace meter assembly 10, and it is stored in the mass change (Δ M) 532, and can measure in time damping parameter (C) 511, and it is stored in the damping change (Δ C) 531.Mass change (Δ M) 532 can represent in the instrument assembly 10 existence of the increase of stream material, and damping change (Δ C) 531 can represent to flow the variation of pipe, comprises material degeneration, corrodes and corrosion, breaks etc.

In operation, instrument electronic device 20 receives three or more vibratory responses 505, and utilizes rigidity program 506 to process vibratory response 505.In one embodiment, these three or more vibratory responses 505 comprise five vibratory responses 505, as will be discussed below.Instrument electronic device 20 is determined limit (λ) 508 and residue (R) 509 from vibratory response 505.Limit (λ) 508 and residue (R) 509 can comprise first order pole and residue, perhaps can comprise second order limit and residue.Instrument electronic device 20 is determined stiffness parameters (K) 510, damping parameter (C) 511 and mass parameter (M) 512 from limit (λ) 508 and residue (R) 509.Instrument electronic device 20 can further be determined the second rigidity (K ₂) 520, can determine stiffness variation (Δ K) 530, from stiffness parameters (K) 510 and the second rigidity (K ₂) 520 can determine stiffness variation (Δ K) 530, and can compare stiffness variation (Δ K) 530 and rigidity tolerance 516.If stiffness variation (Δ K) 530 surpasses rigidity tolerance 516, error log and/or Error processing program that then instrument electronic device 20 can any mode of initialization.Equally, instrument electronic device 20 can further be followed the tracks of damping and mass parameter in time, and can determine and record the second damping (C ₂) the 521 and second quality (M ₂), and the damping change that obtains (Δ C) 531 and mass change (Δ M) 532.Damping change (Δ C) 531 and mass change (Δ M) 532 can similarly be compared with damping tolerance 517 and quality tolerance 518.

Utilize mathematical model can describe the present invention.By open loop, the second order driving model can represent the vibratory response of flowmeter, comprising:

$M\stackrel{.}{x}+C\stackrel{.}{x}+\mathrm{Kx}=f\left(t\right)---\left(10\right)$

Wherein f is the power that is applied to system, and M is the mass parameter of system, and C is damping parameter, and K is stiffness parameters.Item K comprises K=M (ω ₀) 2, and a C comprises C=M2 ζ ω ₀, ω wherein ₀=2 π f ₀, and f ₀It is the resonance frequency of the instrument assembly 10 take Hz as unit.Item ζ comprises the attenuation characteristic measurement that obtains from vibratory response, as previously discussed.In addition, x is the physical displacement distance of vibration, The speed of stream pipe displacement, and

It is acceleration.This is commonly referred to as the MCK model.This formula can be re-arranged into following form:

$({\mathrm{<figure-callout\; id="2"\; label="ms"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">ms</figure-callout>}}^2+\mathrm{cs}+k)X\left(s\right)=F\left(s\right)+(\mathrm{ms}+c)x\left(0\right)+m\stackrel{.}{x}\left(0\right)---\left(11\right)$

Equation (11) can be further processed into the transport function form, ignores simultaneously starting condition.The result is:

Further processing can become first order pole-residue frequency response function form by transformation equation (12), comprising:

$H\left(\mathrm{\&omega;}\right)=\frac{R}{(\mathrm{j\&omega;}-\mathrm{\&lambda;})}+\frac{\stackrel{\&OverBar;}{R}}{(\mathrm{j\&omega;}-\stackrel{\&OverBar;}{\mathrm{\&lambda;}})}---\left(13\right)$

Wherein λ is limit, and R is residue, and (j) comprises-1 square root, and ω is circular excitation frequency (take the radian per second as unit).

This systematic parameter comprise by limit limit intrinsic/resonance frequency (ω _n), damped natural frequency (ω _d) and attenuation characteristic (ζ).

ω _n＝|λ| (14)

ω _d＝imag(λ) (15)

$\mathrm{\&zeta;}=\frac{\mathrm{real}\left(\mathrm{\&lambda;}\right)}{{\mathrm{\&omega;}}_n}---\left(16\right)$

Can obtain stiffness parameters (K), damping parameter (C) and the mass parameter (M) of system from limit and residue.

$M=\frac{1}{2\mathrm{jR}{\mathrm{\&omega;}}_d}---\left(17\right)$

$K={\mathrm{\&omega;}}_n^{2}M---\left(18\right)$

C＝2ζω _nM (19)

Therefore, can calculated rigidity parameter (K), damping parameter (C) and mass parameter (M) according to the good estimation of limit (λ) and residue (R).

Can estimate limit and residue from the frequency response function of measuring.Utilize the direct or iterative calculation method of a certain mode can estimate limit (λ) and residue (R).

Mainly by first formation of equation (13), the complex conjugate item is only shared " residue " part little, approximate constant of response near the driving frequency response.As a result, equation (13) can be simplified as:

$H\left(\mathrm{\&omega;}\right)=\frac{R}{(\mathrm{j\&omega;}-\mathrm{\&lambda;})}---\left(20\right)$

In equation (20), H (ω) is the frequency response function (FRF) of measuring, and it obtains from these three or more vibratory responses.In this was derived, H was made of divided by the power input displacement output.Yet, picking up in the typical situation at the voice coil loudspeaker voice coil of coriolis flowmeter, the FRF of measurement is (namely

) according to speed divided by power.Therefore, equation (20) can be transformed into following form:

$\stackrel{.}{H}\left(\mathrm{\&omega;}\right)=H\left(\mathrm{\&omega;}\right)\&CenterDot;\mathrm{j\&omega;}=\frac{\mathrm{j\&omega;R}}{(\mathrm{j\&omega;}-\mathrm{\&lambda;})}---\left(21\right)$

Equation (21) can further be rearranged the form that limit (λ) and residue (R) are found the solution easily that is paired in.

$\stackrel{.}{H}\mathrm{j\&omega;}-\stackrel{.}{H}\mathrm{\&lambda;}=\mathrm{j\&omega;R}$

$\stackrel{.}{H}=R+\frac{\stackrel{.}{H}}{\mathrm{j\&omega;}}\mathrm{\&lambda;}$

$\left[\begin{array}{cc}1& \frac{\stackrel{.}{H}}{\mathrm{j\&omega;}}\end{array}\right]\left\{\begin{array}{c}R\\ \mathrm{\&lambda;}\end{array}\right\}=\stackrel{.}{H}---\left(22\right)$

Equation (22) forms excessively definite (over-determined) system of equation.Can calculate solving equation (22), so that from speed/power FRF

Determine limit (λ) and residue (R).Item H, R and λ are plural numbers.

In one embodiment, being forced to frequency (forcing frequency) ω is 5 tones.5 tones comprise driving frequency and surpass 2 tones of driving frequency and be lower than 2 tones of driving frequency in this embodiment.These tones can separate 0.5-2Hz with fundamental frequency.Yet, be forced to frequencies omega and can comprise more multi-tone or still less tone, such as driving frequency and on and under 1 tone.Yet, the degree of accuracy of 5 tone impact results and obtain the good compromise of this result between the required processing time.

It may be noted that and in preferred FRF measures, concrete driving frequency and vibratory response are measured two FRF.Pick up (RPO) from driver to the right side and obtain a FRF measurement, and pick up FRF of (LPO) acquisition from driver a to left side and measure.The method is known as single input, many outputs (SIMO).In difference new feature of the present invention, the SIMO technology is used for estimating better limit (λ) and residue (R).Before, two FRF are respectively applied to provide two independent limits (λ) and residue (R) estimation.Can recognize that two FRF share public limit (λ), but independent residue (R _L) and (R _R), two measurements can be advantageously combined that more limit and the residue of robust are determined in order to obtain.

$\left[\begin{array}{ccc}1& 0& \frac{{\stackrel{.}{H}}_{\mathrm{LPO}}}{\mathrm{j\&omega;}}\\ 0& 1& \frac{{\stackrel{.}{H}}_{\mathrm{RPO}}}{\mathrm{j\&omega;}}\end{array}\right]\left\{\begin{array}{c}{R}_L\\ R_R\\ \mathrm{\&lambda;}\end{array}\right\}=\stackrel{.}{H}---\left(23\right)$

Can any amount of mode solving equation (23).In one embodiment, by the recurrent least square method solving equation.In another embodiment, by pseudoinverse technology solving equation.In another embodiment since all measure simultaneously available, Q-R decomposition technique that therefore can Application standard.This Q-R decomposition technique comes into question in modern control theory (Modern Control Theory) (William Brogan, copyright 1991, Prentice Hall, pp.222-224,168-172).

In use, can follow the tracks of in time stiffness parameters (K) and damping parameter (C) and mass parameter (M).For example, statistical technique can be used for determining stiffness parameters (K) any variation (namely stiffness variation (Δ K)) in time.The statistics variations of stiffness parameters (K) can represent to change for the FCF of concrete flowmeter.

The invention provides a kind of stiffness parameters (K) that does not rely on storage or recover the calibration density value.This is opposite with prior art, wherein utilizes known stream material in the factory calibrated operation, in order to obtain can be used for all the in the future density criterion of calibration operation.The invention provides and a kind ofly only obtain stiffness parameters (K) from the vibratory response of flowmeter.The invention provides a kind of rigidity detection/calibration steps that does not need the factory calibrated step.

Fig. 6 is according to the process flow diagram 600 of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter.In step 601, accept three or more vibratory responses.Can receive these three or more vibratory responses from flowmeter.These three or more vibratory responses can comprise basic fundamental frequency response and two or more non-fundamental frequency response.In one embodiment, receive a tone that surpasses fundamental frequency response, and receive a tone that is lower than fundamental frequency response.In another embodiment, receive two or more tones that surpass fundamental frequency response, and receive two or more tones that are lower than fundamental frequency response.

In one embodiment, tone on the fundamental frequency response and under interval equidistantly basically.Alternately, tone interval equidistantly not.

In step 602, produce first order pole-residue frequency response from these three or more vibratory responses.This first order pole-residue frequency response has the form that provides in the equation (23).

In step 603, determine mass parameter (M) from first order pole-residue frequency response.Determine mass parameter (M) by first order pole (λ) and the single order residue (R) of determining vibratory response.Then, determine natural frequency ω from first order pole (λ) and single order residue (R) _n, damped natural frequency ω _d, and attenuation characteristic (ζ).Subsequently, damped natural frequency ω _d, residue (R) and imaginary term (j) are inserted in the equation (17), in order to obtain mass parameter (M).

In step 604, find the solution definite stiffness parameters (K) from equation (18).This is found the solution and adopts natural frequency ω _nAnd be inserted into equation (18) from the mass parameter (M) that step 603 is determined, in order to obtain stiffness parameters (K).

In step 605, find the solution definite damping parameter (C) from equation (19).This is found the solution and adopts attenuation characteristic (ζ), natural frequency ω _n, and the mass parameter (M) of determining.

Fig. 7 shows the embodiment that limit (λ) according to embodiments of the invention and residue (R) are found the solution.This embodiment is followed equation (23).The FRF input is on the left side of figure.These FRF input is five frequencies (four frequency test signal and driving frequency) in this embodiment, calculates the FRF coefficient at this frequency place.FRF_L and FRF_R input is to pick up multiple FRF coefficient at the driver that those frequency places are calculated, corresponding in the equation (23) With These FRF coefficients enter the B input of QR solver piece 701.Take item by item as the basis, be formed for the A matrix of QR solver piece 701 from the FRF coefficient divided by j ω, and this A matrix comprises 1 and 0 row, in order to meet equation (23).This matrix is shaped as suitable [10 * 3] complex dimension again, and enters the A input of QR solver piece 701.The x vector output of QR solver piece 701 comprises left and right residue R _LAnd R _RWith limit λ.In order to process, these outputs spread out of from QR piece 701, in order to produce systematic parameter.

Fig. 8 is the M that illustrates according to embodiments of the invention, the calcspar of the calculating of C and K systematic parameter.This embodiment is determined M, C and K systematic parameter from limit and the residue estimation of each equation (14-16) and equation (17-19).These residues are pure imaginary numbers to real normal modal model.Yet, owing to the noise in the measurement data with owing to models fitting numerical precision problem, will usually have a certain real part.Therefore, use the absolute value of residue, it has each equation (17) except the similar effect of j.Utilize limit and calculus of residues quality and the rigidity of each equation (17-18).It may be noted that to have " left side " and " right side " quality and rigidity the quality rigidity of namely calculating from the FRF of LPO/ driver and RPO/ driver.Because the asymmetry of coil and magnet and structure self, from right to left, quality can be different with the rigidity estimation.The variation of difference or difference ratio represent the non-uniform change of quality or rigidity and can be used for providing about the additional diagnostic information for the integrality of the variation of FCF or flow.

Two other outputs of calculating from systematic parameter are ratio of damping, Z (zeta) or ζ, and natural frequency ω _nThis embodiment provides more to cross and determines or the better global parameter group of estimating.

ω _nEstimation obtained good quality check for the closed loop drive system.If drive real work at the resonance place, so for the natural frequency estimation, driving frequency will meet (agreeto) several millihertzs hereby in.If this difference greater than several millihertzs hereby, then can arrange warning label, the expression drive system is worked irrelevantly, and perhaps current rigidity estimation is shady.

Fig. 9 shows the rigidity estimating system based on FRF according to the integral body of embodiments of the invention.Exist to seven different inputs of rigidity estimation subsystem, it represents (five left sides up, and two at rightmost) by the pentagon as signal source." RawDrive " and " RawPOs " input is the original reading of pick-up voltage and drive current.For example by extracting (decimation), then these signals, and are fed and enter FRF coefficient estimate subsystem to 2kHz by down sample.This " CmdmA " input is the instruction current that obtains from the output of the digital drive system of correspondence." StiffnessEnable (rigidity enables) " estimation is the logic input, allows the digital drive system to control when the FCF checking algorithm is effective." freq " input is driving frequency, as estimating by the digital drive system.It is the input to measuring signal generator subsystem and Rigidity Calculation subsystem.

FRF Rigidity Calculation piece 902 output system parameter estimation M and K left and the right side and ζ and FreqEst.These are the Main Diagnosis output for the FCF verification.This figure also illustrates frequency differential alarming block 903 and the frequency differential error block 904 of implementing driving quality check discussed above by the natural frequency of relatively driving frequency and estimation.

Measuring FRF needs current measurement usually, needs additional modulus (A/D) converter.Yet this embodiment utilizes the command current of calibration, has avoided the needs for additional A/D converter.CL input selection piece 906 and CL output calibration piece 907 are carried out calibration algorithm.This calibration steps utilizes " test signal FRF " piece 901, in order to calculate at a state place of steering logic, actual (RawDrive) electric current is to the frequency response function of command current (CmdmA).In FCF check logic state procedure, by calculate and proofread and correct the FRF between original POs and the instruction current for the raw data of command current FRF coefficient, in order to provide the FRFs that processes for further.

FRF rigidity estimating algorithm is in chart central authorities' left side place outputs " test signal " output of figure.This test signal output is added into the excitation that drives four test frequency places that order before being included in output immediately.In the time can carrying out the FCF verification, these test signals are added into digital drive signals.

This logic is such: when the FCF verification is turn-offed, digital drive signals is just by switch or other device, in this case, by interpolation filter, from it basic speed (generally being 4kHz) to suitable output speed (generally being 8kHz) to this digital drive signals of up-sampling.In the time that the FCF verification can be carried out, be added into digital drive signals from 2 to 4kHz by the test signal to up-sampling.Then this driving signal is made of closed loop driving frequency signal and 4 test tones, and then it is all by to up-sampling filter.

The expectation of FCF checking routine is transparent to drive system.In one embodiment, from picking up the removal test signal, so that good frequency and amplitude evaluation are guaranteed in driving to closed loop.This utilization is tuned to one group of notch filter of the precise frequency of test signal and finishes.

In another embodiment, limit-residue method can adopt second order limit-residue frequency response function, thereby realizes better result.Second order limit-residue provides more real match than first order pole-residue method to real number.Compromise is larger numerical value complicacy and the processing time of increase.

The MCK embodiment of rigidity estimation starts from simple second order model, shown in following equation (24).Owing to the picking up of flowmeter survey speed, be not the position, equation is by micronized, and then estimated at concrete frequencies omega place.

$H\left(s\right)=\frac{X\left(s\right)}{F\left(s\right)}=\frac{1}{{\mathrm{<figure-callout\; id="2"\; label="Ms"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">Ms</figure-callout>}}^2+\mathrm{Cs}+K}$

$\stackrel{.}{H}\left(s\right)=\frac{\stackrel{.}{X}\left(s\right)}{F\left(s\right)}=\frac{s}{{\mathrm{<figure-callout\; id="2"\; label="Ms"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">Ms</figure-callout>}}^2+\mathrm{Cs}+K}$

$\stackrel{.}{H}\left(\mathrm{\&omega;}\right)=\frac{\stackrel{.}{X}\left(\mathrm{\&omega;}\right)}{F\left(\mathrm{\&omega;}\right)}=\frac{\mathrm{j\&omega;}}{-M{\mathrm{\&omega;}}^2+\mathrm{jC\&omega;}+K}---\left(24\right)$

Since target be measurement from drive current (or power) and pick-up voltage (or speed) to M, C and K find the solution, therefore expediently, rewrite equation (24) is so that the separation unknown quantity.This produces equation (25).

$K-M{\mathrm{\&omega;}}^2+\mathrm{jC\&omega;}=\frac{\mathrm{j\&omega;}}{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)}---\left(25\right)$

At this some place, equation can be separated into real part and imaginary part.

$K-M{\mathrm{\&omega;}}^2=\mathrm{Re}\left\{\frac{\mathrm{j\&omega;}}{\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)}\right\}$

$\mathrm{C\&omega;}=\mathrm{Im}\left\{\frac{\mathrm{j\&omega;}}{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)}\right\}---\left(26\right)$

Launch

Equation (26) can be write as again:

$K-M{\mathrm{\&omega;}}^2=\frac{\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{{\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}$

$\mathrm{C\&omega;}=\frac{\mathrm{\&omega;Re}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{{\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}---\left(27\right)$

The second equation is simple, Algebraic Algorithm now.For the first of further reduced equation, adopt the resonant drive frequency of measuring.Because

Therefore can set up:

$K-K/{\mathrm{\&omega;}}_{\mathrm{<figure-callout\; id="2"\; label="n"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">n</figure-callout>}}^2=\frac{\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{{\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}$

$K\left(\frac{{\mathrm{\&omega;}}_n^2-{\mathrm{\&omega;}}^2}{{\mathrm{\&omega;}}_n^2}\right)=\frac{\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{{\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}$

$K=\frac{{\mathrm{\&omega;}}_{\mathrm{<figure-callout\; id="2"\; label="n"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">n</figure-callout>}}^2\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{({\mathrm{\&omega;}}_n^2-{\mathrm{\&omega;}}^2){\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}---\left(28\right)$

As long as ω ≠ ω _nFind the solution from this for K and to return M, in equation (29), provide M, three solutions of C and K.

$K=\frac{{\mathrm{\&omega;}}_{\mathrm{<figure-callout\; id="2"\; label="n"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">n</figure-callout>}}^2\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{({\mathrm{\&omega;}}_n^2-{\mathrm{\&omega;}}^2){\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}$

$M=\frac{K}{{\mathrm{\&omega;}}_{\mathrm{<figure-callout\; id="2"\; label="n"\; filenames="HSA00000109999600071.png,HSA00000109999600091.png"\; state="\{\{state\}\}">n</figure-callout>}}^2}$

$C=\frac{\mathrm{Re}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{{\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}---\left(29\right)$

It may be noted that given resonance frequency omega _n, a concrete frequencies omega ₁The driver at place picks up FRF is enough to solving equation, and definite parameter M, C and K.This is useful especially; When obtaining FRF at a plurality of frequencies place, be each estimation average of each coefficient simply to the least square fitting of data.This is than the simpler good processing mode of the pseudoinverse that typically will have to be performed.It may be noted that ω ≠ ω _nRestriction got rid of the use of finding the solution middle resonant drive FRF to K or M.This is not surprising especially, because only determine the height of resonance place peak value by damping.A latent defect of the method is: must not rely on each other from the parameter of left and right picks up data estimation.This and limit-residue method forms contrast, in this case, leftly and right picks up to estimate that identical limit obtains some advantages by limiting, and regardless of their difference on amplitude.

Figure 10 is according to the process flow diagram 1000 of embodiments of the invention for the method for the stiffness parameters (K) of determining flowmeter.In step 1001, receive three or more vibratory responses, as previously discussed.

In step 1002, produce second order limit-residue frequency response from these three or more vibratory responses.This second order limit-residue frequency response has the form that provides in equation (24).

In step 1003, find the solution definite stiffness parameters (K) from equation (29).This is found the solution and adopts natural frequency ω _n, one or more periodicity pitch ω, the imaginary part of FRF is (namely

Imaginary part), and the amplitude of FRF is (namely

Absolute value).

In step 1004, determine mass parameter (M) from second order limit-residue frequency response.From equation (29) find the solution definite mass parameter (M), and utilize stiffness parameters (K) and natural frequency ω _nObtain this mass parameter (M).

In step 1005, determine damping parameter (C) from second order limit-residue frequency response.From equation (29) find the solution definite damping parameter (C), and utilize these one or more periodicity pitch ω, the real part of FRF (namely

Real part), and the amplitude of FRF is (namely

Absolute value) obtain this damping parameter (C).

Figure 11 shows according to embodiments of the invention from the M of equation (29) for second order limit-residue response, the embodiment that C and K find the solution.Input is revealed as the oval input port at the place, the left side of figure.These are the driving frequency ω that measures _ drive, and it is used for equation (29) as ω _n, five frequencies are calculated FRF coefficients (four frequency test signals and driving frequency are represented by ω _ test) at these five frequency places, and the driver that calculates at those frequency places pick up multiple FRF coefficient (

Or Hdot).It abandons (discard) driving frequency FRF by the selector switch piece, because can not be used for as find the solution at M and the K of front description.K finds the solution and is calculated as:

$K=\frac{\mathrm{\&omega;Im}\left\{\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right\}}{(1-{\mathrm{\&omega;}}^2/{\mathrm{\&omega;}}_n^2){\left|\stackrel{.}{H}\left(\mathrm{\&omega;}\right)\right|}^2}---\left(30\right)$

It is the equivalent form of value of the method for solving that provides in the equation (29).For finding the solution of C be with equation (29) in derivation find the solution identical form, and directly calculate M from the method for solving for K.It may be noted that the operation of averaging is applied to each coefficient estimate.This mean value result in the method for solving is the least square fitting to the input data.At last, provide M, C and K estimation, attenuation characteristic (ζ or Z) is calculated as:

$\mathrm{\&zeta;}=\frac{{\mathrm{C\&omega;}}_n}{2K}---\left(31\right)$

Attenuation characteristic (ζ) is considered to the parameter more useful than damping parameters C.Therefore, mass M, stiffness K and attenuation characteristic (ζ) are the output of measuring.

Figure 12 shows the rigidity estimating system based on FRF according to the integral body of embodiments of the invention.Exist to six of system different inputs, by the pentagon as signal source represent (three left sides up, and three below the right)." RawDrive " and " RawPOs " input is from the original reading that picks up with drive current.By withdrawal device (Decimator) piece 1201 down samples these to 2kHz, and then these are presented and enter FRF coefficient estimate subsystem." DriveDemod " input is sine and the cosine signal from the driving frequency of digital drive system acquisition.These signals are by the sinusoidal curve combination that produces with the test frequency place, and are fed as the basis that is used for demodulation and enter FRF coefficient estimate subsystem." StiffnessEnable " estimation is the logic input, allows the digital drive system to control when the rigidity estimating algorithm is effective." freq " input is driving frequency, as estimating by the digital drive system.It is the input to test signal generation block 1204 and Rigidity Calculation piece 1206." Temp " input is the temperature reading from flowmeter that is imported in the temperature correction piece 1207.FRF rigidity estimating algorithm output system parameter estimation, and " test signal " output of locating in the leftmost side of figure.This test signal output comprises the excitation at four test frequency places that are added into drive command.

These input and output are formed to the interface whole (bulk) of digital drive.To the output of drive assembly, test signal is added into immediately drives order.In order to make this FCF checking routine transparent to drive system, must be from picking up the removal test signal.This can utilize one group of notch filter of the precise frequency that is tuned to test signal to finish in one embodiment.

The test signal FRF piece 1208 of Figure 11 is carried out demodulation.It is demodulated at each place of five incoming frequencies (frequency test signal of four generations and driving frequency) to pick up and drive signal.After utilizing sinusoidal and the cosine basis carries out complex demodulation, the real number of each signal and imaginary part are extracted to lower frequency downwards, and are low pass filtered to 0.4Hz.These signals must be unpolluted in this zone, because any spectrum component among the 0.4Hz of test signal will be not suppressed, and will manifest with output.The FRF that then is used for this frequency place of estimation with the plural coefficient of drive current that picks up that is used for each frequency place.To a plurality of sample mean power spectrum, and output is estimated than low rate FRF.

Can adopt according to instrument electronic device of the present invention and method according to arbitrary embodiment, thereby several advantages are provided, if desired.The invention provides the stiffness parameters (K) of the stream pipe stiffness that basically relates to flowmeter.The invention provides do not rely on for generation of storage or recover the stiffness parameters (K) of calibration value.Only the invention provides the stiffness parameters (K) that obtains from the vibratory response of flowmeter.Similarly, the invention provides mass parameter (M) and damping parameter (C) from vibratory response.

The invention provides does not need the rigidity of factory calibrated step detection/calibration steps.The present invention can carry out rigidity/FCF calibration steps in the field.The present invention can locate to carry out rigidity/FCF calibration steps at any time.The present invention can carry out the rigidity that do not need calibration testing ring and/or known flow material/FCF calibration steps.The present invention can carry out the in time rigidity of the stiffness variation of definite flowmeter/FCF calibration steps.

Claims (41)

1. the instrument electronic device (20) that is used for flowmeter (5), this instrument electronic device (20) comprises that for the interface (201) that receives from three or more vibratory responses of flowmeter (5) wherein these three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response; And the disposal system (203) of communicating by letter with interface (201), wherein this instrument electronic device (20) further comprises:

Disposal system (203) is arranged to receive these three or more the vibratory responses from interface (201), produce single orders or higher order pole-residue frequency response function more from these three or more vibratory responses, and from this single order or more higher order pole-residue frequency response function determine at least stiffness parameters K.

2. the instrument electronic device of claim 1 (20), wherein disposal system (203) further be arranged to from this single order or more higher order pole-residue frequency response function determine damping parameter C.

3. the instrument electronic device of claim 1 (20), wherein disposal system (203) further be arranged to from this single order or more higher order pole-residue frequency response function determine mass parameter M.

4. the instrument electronic device of claim 1 (20), wherein disposal system (203) further be arranged to from this single order or more higher order pole-residue frequency response function calculate limit λ, left residue R _LWith right residue R _R

5. the instrument electronic device of claim 1 (20), wherein these three or more vibratory responses comprise at least one tone that is higher than fundamental frequency response and are lower than at least one tone of fundamental frequency response.

6. the instrument electronic device of claim 1 (20), wherein these three or more vibratory responses comprise at least two tones that are higher than fundamental frequency response and at least two tones that are lower than fundamental frequency response.

7. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function.

8. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

Wherein the R item comprises residue, and the λ item comprises that limit and ω item comprise circular excitation frequency.

9. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

And wherein according to equation M=1/2jR ω _d, K=(ω _n) ²M and C=2 ζ ω _nM determines stiffness parameters K, damping parameter C and mass parameter M, and wherein the M item comprises mass parameter, the ζ item comprises attenuation characteristic, ω _nItem comprises natural frequency and ω _dItem comprises damped natural frequency.

10. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function.

11. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises

Wherein the M item comprises mass parameter, and the C item comprises damping parameter, and the K item comprises stiffness parameters, and the ω item comprises circular excitation frequency.

12. the instrument electronic device of claim 1 (20), wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises

And foundation wherein $K=\left({\left({\mathrm{\&omega;}}_n\right)}^2\mathrm{\&omega;Im}\right[\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\left]\right)/({\left({\mathrm{\&omega;}}_n\right)}^2-{\mathrm{\&omega;}}^2){\left|\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\right|}^2)$ Determine stiffness parameters K, according to M=K/ (ω _n) ²Determine mass parameter M, and foundation

Determine damping parameter C, wherein the ω item comprises circular excitation frequency.

13. a method that is used for the stiffness parameters K of definite flowmeter, the method comprises:

Receive three or more vibratory responses, wherein these three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response;

Produce single orders or higher order pole-residue frequency response function more from these three or more vibratory responses; And

From this single order or more higher order pole-residue frequency response function determine at least stiffness parameters K.

14. the method for claim 13, further comprise from this single order or more higher order pole-residue frequency response function determine damping parameter C.

15. the method for claim 13, further comprise from this single order or more higher order pole-residue frequency response function determine mass parameter M.

16. the method for claim 13, further comprise from this single order or more higher order pole-residue frequency response function calculate limit λ, left residue R _LWith right residue R _R

17. the method for claim 13, wherein these three or more vibratory responses comprise at least one tone that is higher than fundamental frequency response and at least one tone that is lower than fundamental frequency response.

18. the method for claim 13, wherein these three or more vibratory responses comprise at least two tones that are higher than fundamental frequency response and at least two tones that are lower than fundamental frequency response.

19. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function.

20. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

Wherein the R item comprises residue, and the λ item comprises that limit and ω item comprise circular excitation frequency.

21. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

And wherein according to equation M=1/2jR ω _d, K=(ω _n) ²M and C=2 ζ ω _nM determines stiffness parameters K, damping parameter C and mass parameter M, and wherein the M item comprises mass parameter, the ζ item comprises attenuation characteristic, ω _nItem comprises natural frequency and ω _dItem comprises damped natural frequency.

22. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function.

23. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises

Wherein the M item comprises mass parameter, and the C item comprises damping parameter, and the K item comprises stiffness parameters, and the ω item comprises circular excitation frequency.

24. the method for claim 13, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises $\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)=\frac{\stackrel{\&CenterDot;}{X}\left(\mathrm{\&omega;}\right)}{F\left(\mathrm{\&omega;}\right)}=\frac{\mathrm{j\&omega;}}{-M{\mathrm{\&omega;}}^2+\mathrm{jC\&omega;}+K},$ And foundation wherein $K=({\left({\mathrm{\&omega;}}_n\right)}^2\mathrm{\&omega;}$ $\mathrm{Im}\left[\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\right])/({\left({\mathrm{\&omega;}}_n\right)}^2-{\mathrm{\&omega;}}^n){\left|\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\right|}^2)$ Determine stiffness parameters K, according to M=K/ (ω _n) ²Determine mass parameter M, and foundation

Determine damping parameter C, wherein the ω item comprises circular excitation frequency.

25. the method for claim 13 further comprises:

At the second time t ₂The place receives second group of three or more vibratory response from flowmeter;

Produce the second stiffness characteristics K from these second group of three or more vibratory response ₂

Compare the second stiffness characteristics K ₂With stiffness parameters K; And

If the second stiffness characteristics K ₂Different from predetermined tolerance from stiffness parameters K are then detected stiffness variation Δ K.

26. the method for claim 25 further comprises if the second stiffness characteristics K ₂Different from the predetermined stiffness tolerance from stiffness parameters K are then detected stiffness variation Δ K.

27. the method for claim 25 further comprises from stiffness parameters K and the second stiffness characteristics K ₂relatively quantize stiffness variation Δ K.

28. a method that is used for the stiffness variation Δ K of definite flowmeter, the method comprises:

Receive three or more vibratory responses, wherein these three or more vibratory responses comprise basic fundamental frequency response and two or more non-fundamental frequency response;

Produce single orders or higher order pole-residue frequency response function more from these three or more vibratory responses;

From this single order or more higher order pole-residue frequency response function determine at least stiffness parameters K;

At the second time t ₂The place receives second group of three or more vibratory response from flowmeter;

Produce the second stiffness characteristics K from these second group of three or more vibratory response ₂

Compare the second stiffness characteristics K ₂With stiffness parameters K; And

If the second stiffness characteristics K ₂Different from predetermined tolerance from stiffness parameters K are then detected stiffness variation Δ K.

29. the method for claim 28 further comprises if the second stiffness characteristics K ₂Different from the predetermined stiffness tolerance from stiffness parameters K are then detected stiffness variation Δ K.

30. the method for claim 28 further comprises from stiffness parameters K and the second stiffness characteristics K ₂relatively quantize stiffness variation Δ K.

31. the method for claim 28, wherein said determine to comprise from this single order or more higher order pole-residue frequency response function further determine damping parameter C.

32. the method for claim 28, wherein said determine to comprise from this single order or more higher order pole-residue frequency response function further determine mass parameter M.

33. the method for claim 28, wherein said determine further to comprise from this single order or more higher order pole-residue frequency response function calculate limit λ, left residue R _LWith right residue R _R

34. the method for claim 28, wherein these three or more vibratory responses comprise at least one tone that is higher than fundamental frequency response and at least one tone that is lower than fundamental frequency response.

35. the method for claim 28, wherein these three or more vibratory responses comprise at least two tones that are higher than fundamental frequency response and at least two tones that are lower than fundamental frequency response.

36. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function.

37. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

Wherein the R item comprises residue, and the λ item comprises that limit and ω item comprise circular excitation frequency.

38. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise first order pole-residue frequency response function, this first order pole-residue frequency response function comprises

And wherein according to equation M=1/2jR ω _d, K=(ω _n) ²M and C=2 ζ ω _nM determines stiffness parameters K, damping parameter C and mass parameter M, and wherein the M item comprises mass parameter, the ζ item comprises attenuation characteristic, ω _nItem comprises natural frequency and ω _dItem comprises damped natural frequency.

39. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function.

40. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises

Wherein the M item comprises mass parameter, and the C item comprises damping parameter, and the K item comprises stiffness parameters, and the ω item comprises circular excitation frequency.

41. the method for claim 28, wherein this single order or more higher order pole-residue frequency response function comprise second order limit-residue frequency response function, this second order limit-residue frequency response function comprises $\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)=\frac{\stackrel{\&CenterDot;}{X}\left(\mathrm{\&omega;}\right)}{F\left(\mathrm{\&omega;}\right)}=\frac{\mathrm{j\&omega;}}{-M{\mathrm{\&omega;}}^2+\mathrm{jC\&omega;}+K},$ And foundation wherein $K=\left({\left({\mathrm{\&omega;}}_n\right)}^2\mathrm{\&omega;Im}\right[\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\left]\right)/({\left({\mathrm{\&omega;}}_n\right)}^2-{\mathrm{\&omega;}}^2){\left|\stackrel{\&CenterDot;}{H}\left(\mathrm{\&omega;}\right)\right|}^2)$ Determine stiffness parameters K, according to M=K/ (ω _n) ²Determine mass parameter M, and foundation

Determine damping parameter C, wherein the ω item comprises circular excitation frequency.

Images