CN111649908A - Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction - Google Patents
Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction Download PDFInfo
- Publication number
- CN111649908A CN111649908A CN202010546766.5A CN202010546766A CN111649908A CN 111649908 A CN111649908 A CN 111649908A CN 202010546766 A CN202010546766 A CN 202010546766A CN 111649908 A CN111649908 A CN 111649908A
- Authority
- CN
- China
- Prior art keywords
- balance
- force measurement
- dynamic characteristic
- model
- characteristic compensation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
- G01M9/062—Wind tunnel balances; Holding devices combined with measuring arrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention relates to a heavenly translation dynamic characteristic compensation method and a device based on wavelet reconstruction, wherein the method comprises the following steps: installing a model, a balance and a support rod to obtain natural frequency; respectively installing unidirectional accelerometers in the axial direction, the normal direction and the lateral direction of the model; simulating an impact process under a windless condition, acquiring force measurement signals output by a balance and a unidirectional accelerometer, and respectively calculating dynamic characteristic compensation coefficients in three directions; performing a wind tunnel force measurement test, impacting the model and acquiring force measurement data output by the balance and the unidirectional accelerometer; performing wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer in three directions, performing high-frequency filtering by taking the inherent frequency as a filtering cut-off frequency, and performing dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient; and calculating aerodynamic force of the wind tunnel force measurement test. The invention can compensate the dynamic characteristic of the balance signal, improves the reliability of the force measurement data and has strong universality.
Description
Technical Field
The invention relates to the technical field of wind tunnel force measurement tests, in particular to a astronomical dynamic characteristic compensation method and device based on wavelet reconstruction.
Background
Balances for measuring strain forces, such as six-component balances, rod balances, ring balances, box balances, etc., are important devices in wind tunnel force measurements. When the wind tunnel test is carried out, the balance is arranged in the model, when the model is impacted, the balance can vibrate along with the model, and the model can be recovered to a stable state after a period of time, and in the process, the force measurement data of the balance is interfered by the vibration, so that the precision of the test data of the test is seriously influenced.
Therefore, in order to overcome the defects, a technical scheme capable of correcting the influence of model oscillation caused by impact in the wind tunnel force measurement test on the balance force measurement needs to be provided.
Disclosure of Invention
Technical problem to be solved
The invention aims to solve the technical problem that the influence of model oscillation caused by impact on measurement is not considered in the prior art, so that the reliability of the force measurement result of a balance is low.
(II) technical scheme
In order to solve the technical problem, the invention provides a compensation method for the astronomical motion state characteristics based on wavelet reconstruction, which comprises the following steps:
s1, installing the model, the balance and the support rod, and acquiring the integral natural frequency of the model, the balance and the support rod;
s2, respectively installing unidirectional accelerometers in the axial direction, the normal direction and the lateral direction of the model, wherein the frequency response of the unidirectional accelerometers is greater than the natural frequency;
s3, under the windless condition, simulating the impact process of the model in a wind tunnel force measurement test, acquiring force measurement signals output by the balance and the unidirectional accelerometers, and respectively calculating dynamic characteristic compensation coefficients in three directions;
s4, performing a wind tunnel force measurement test, impacting the model, and acquiring force measurement data output by the balance and each one-way accelerometer;
s5, performing wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer in three directions, performing high-frequency filtering by taking the inherent frequency as a filtering cut-off frequency, and performing dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient to obtain a compensated balance signal;
and S6, calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
Preferably, in step S3, when the dynamic characteristic compensation coefficients in the three directions are calculated, the following formula is used for each direction:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance, N representing the total number of points of the force measurement signal participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point of the force measuring signal output by the balance, YiAnd the voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer is represented.
Preferably, in step S4, when the force measurement data output by the balance and the unidirectional accelerometers are obtained, the original data output by the balance and the unidirectional accelerometers in three directions are recorded, an effective data segment and zero point data in the original data are determined, and zero point data are subtracted from the effective data segment, where the zero point data are obtained from a zero point state of stable damped oscillation after the model is impacted.
Preferably, in step S5, when performing wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer in three directions, the decomposition is performed by using db10 wavelet basis.
Preferably, in the step S5, when the db10 wavelet basis is adopted for decomposition, the decomposition is performed 5-10 times.
Preferably, in step S5, when the dynamic characteristic compensation is performed on the force measurement data output by the balance according to the corresponding dynamic characteristic compensation coefficient, the compensated balance signal S in any direction is obtainedend＝sWT-ty×sWT-JSDWherein s isWTRepresenting the reconstructed signal, s, of the balance after high-frequency filteringWT-JSDAnd representing a reconstructed signal of the unidirectional accelerometer obtained after high-frequency filtering.
Preferably, in the step S1, when the natural frequency of the model, the balance and the whole strut is obtained, a manual tapping mode is adopted, and the method includes the following steps:
knocking the center of mass of the model, the balance and the support rod by a force hammer, and recording the force measurement data output by the balance; and determining the integral natural frequency of the model, the balance and the support rod by performing frequency spectrum analysis on the force measurement data output by the balance.
Preferably, in the step S3, when an impact process received by the model in the wind tunnel force measurement test is simulated, the center of mass of the model, the balance and the support rod is knocked by the hammer, and a direction in which the hammer knocks the model is consistent with a direction in which the hammer impacts the model in the wind tunnel force measurement test in the step S4.
The invention also provides a compensation device for the characteristic of the astronomical motion state based on the wavelet reconstruction, which comprises the following components: the system comprises three unidirectional accelerometers, a coefficient unit, a signal unit and a calculation unit; wherein the content of the first and second substances,
the unidirectional accelerometers are respectively arranged in the axial direction, the normal direction and the lateral direction of the model, and the frequency response of the unidirectional accelerometers is greater than the integral natural frequency of the model, the balance and the support rod;
the coefficient unit is connected with the balance and each one-way accelerometer and is used for acquiring force measurement signals output by the balance and each one-way accelerometer in the process of simulating the impact on the model in a wind tunnel force measurement test under a windless condition and respectively calculating dynamic characteristic compensation coefficients in three directions;
the signal unit is connected with the coefficient unit, the balance and each one-way accelerometer and is used for acquiring force measurement data output by the balance and each one-way accelerometer during a wind tunnel force measurement test, performing wavelet decomposition on the force measurement data output by the balance and each one-way accelerometer in three directions, performing high-frequency filtering by taking inherent frequency as filtering cut-off frequency, and performing dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient to acquire a compensated balance signal;
and the calculation unit is connected with the signal unit and is used for calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
Preferably, when the coefficient unit calculates the dynamic characteristic compensation coefficients in three directions, the following formula is adopted for each direction to calculate:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance, N representing the total number of points of the force measurement signal participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point of the force measuring signal output by the balance, YiAnd the voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer is represented.
(III) advantageous effects
The technical scheme of the invention has the following advantages: the invention provides a compensation method and a device for astronomical dynamic characteristics based on wavelet reconstruction. The technical scheme provided by the invention can effectively improve the reliability of the force measurement result of the balance, is easy to realize, has strong universality and is suitable for dynamic characteristic compensation of various balances.
Drawings
Fig. 1 is a schematic diagram of steps of a astronomical dynamic characteristic compensation method based on wavelet reconstruction in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
As shown in fig. 1, a skyward dynamic characteristic compensation method based on wavelet reconstruction provided by an embodiment of the present invention includes the following steps:
s1, installing the model, the balance and the support rod used in the wind tunnel force measurement test, and acquiring the integral natural frequency f of the model, the balance and the support rod0。
f0The natural frequency obtained by considering the model, the balance and the supporting rod as a whole can be obtained in a manner that the natural frequency is obtained by referring to the prior art, and in order to obtain more accurate natural frequency, the knocking position needs to be as close to the mass center of the whole as possible. Preferably, in step S1, when acquiring the natural frequency of the model, the balance and the whole strut, the method of manual tapping is adopted, which includes the following steps:
knocking the center of mass of the model, the balance and the support rod by a force hammer, and recording the force measurement data output by the balance; determining the natural frequency f of the model, the balance and the support rod as a whole by carrying out frequency spectrum analysis on the force measurement data output by the balance0. The balance is preferably a six-component balance, although other component balances, as well as bar balances, ring balances, cassette balances or other shaped balances, may also be used.
S2, respectively installing unidirectional accelerometers in the axial direction, the normal direction and the lateral direction of the model, wherein the frequency response of the unidirectional accelerometers is greater than the integral natural frequency f of the model, the balance and the support rod obtained in the step S10。
In the step, the installation of the unidirectional accelerometer does not need to distinguish the positive direction and the negative direction, and the direction adjustment is realized by the positive sign and the negative sign of the calculated dynamic characteristic compensation coefficient in the subsequent dynamic characteristic compensation process.
And S3, simulating the impact process of the model in the wind tunnel force measurement test under the windless condition, acquiring the force measurement signals output by the balance and the unidirectional accelerometers, and respectively calculating the dynamic characteristic compensation coefficients in three directions.
The step S3 is performed without opening the wind tunnel, and is intended to obtain the dynamic characteristic compensation coefficients corresponding to the axial direction, the normal direction, and the lateral direction by a dynamic calibration method. Three directions, namely axial direction, normal direction and lateral direction, and the corresponding three balance output parameters are VN, VZ and VA which are common parameters of the balance used for measuring the strain force.
And S4, performing a wind tunnel force measurement test, impacting the model, and acquiring force measurement data output by the balance and each unidirectional accelerometer.
In the wind tunnel force test, the model is impacted and oscillates, and in step S4, data of the complete process from the stable state of the model before impact to the stable state of the model after impact is acquired, so as to further analyze the impact process.
S5, respectively carrying out wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer obtained in the step S4 in three directions, and carrying out wavelet decomposition on the force measurement data with the natural frequency f0And carrying out high-frequency filtering on the filtering cut-off frequency to realize wavelet reconstruction, and carrying out dynamic characteristic compensation on the force measurement data output by the balance according to the corresponding dynamic characteristic compensation coefficient to obtain a compensated balance signal.
In step S5, dynamic characteristic compensation is performed according to the dynamic characteristic compensation coefficients obtained in step S3, the force measurement data output by the balance obtained in step S4, and the force measurement data output by the unidirectional accelerometer, in the axial direction, the normal direction, and the lateral direction, respectively, to finally obtain compensated balance signals corresponding to the three directions.
And S6, calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
In step S6, the aerodynamic force (coefficient) is calculated from the compensated balance signals (i.e., the balance output parameters VN, VZ, VA with the dynamic characteristics compensated). How to calculate the aerodynamic force (coefficient) of the wind tunnel test according to the balance output parameters VN, VZ, VA is the prior art, and is not described herein again.
Preferably, in step S3, when simulating the impact process on the model in the wind tunnel force measurement test under the windless condition, a manual tapping mode is adopted, and the center of mass of the model, the balance and the support rod is tapped by the force hammer, and the direction of the force hammer tapping the model is consistent with the direction of the impact model in the wind tunnel force measurement test in step S4. The force of the force hammer can not be consistent with the impact force of the impact model in the wind tunnel force measurement test.
Preferably, in step S3, when the dynamic characteristic compensation coefficients in the three directions are calculated, the following formula is used for each direction:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance (corresponding to the direction), N representing the total point number of the force measuring signals participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point, Y, of the force measuring signal output by the balance (corresponding to the direction)iThe voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer (corresponding to the direction) is represented.
Further, in order to ensure the reliability of the calculation result, the force measurement signals output by the balance and the unidirectional accelerometer which respectively participate in the calculation in the three directions are taken from the same time period, and the number N of the force measurement signals is the same.
When the model is subjected to the impact process in the wind tunnel force measurement test in the step S3, after the impact, the model is not subjected to an external force, the force measurement signals output by the balance and the unidirectional accelerometers exhibit damping attenuation curve attenuation, and particularly, in order to more accurately obtain the dynamic characteristic compensation coefficient, the force measurement signals (output by the balance and the unidirectional accelerometers) involved in the calculation in the step S3 should be taken from the relatively smooth section in the corresponding damping attenuation curve.
Preferably, in step S4, when acquiring force measurement data output by the balance and the unidirectional accelerometers, respectively recording raw data output by the balance and the unidirectional accelerometers in three directions, determining an effective data segment and zero point data in the raw data, and subtracting the zero point data from the effective data segment, wherein the zero point data is obtained from a zero point state of the stable damped oscillation after (returning to) the model is impacted.
For either direction, let sorgFor valid data sections, s, in the raw data output by the balanceJSD-orgIs a valid data segment, s, in the raw data output by the unidirectional accelerometerzeroZero point data s of stable damping oscillation of model after impact in original data output by balanceJSD-zeroRespectively deducting zero point data of the model after being impacted from the original data of the balance and the unidirectional accelerometer to obtain force measurement data s and s output by the balance and the unidirectional accelerometerJSDI.e. s ═ sorg-szeroAnd sJSD＝sJSD-org-sJSD-zero。
Preferably, in step S5, when the force measurement data output by the balance and the unidirectional accelerometer are wavelet decomposed in three directions, the decomposition is performed by using db10 wavelet basis. Further, when the db10 wavelet basis is used for decomposition, the decomposition is performed 5-10 times, wherein the decomposition is preferably performed 10 times. The decomposition times affect the accuracy of the reconstructed signal obtained after the high-frequency filtering, and the larger the decomposition times, the finer the decomposition, and the smaller the decomposition times, which may result in the low signal spectrum resolution.
Preferably, in step S5, when performing high-frequency filtering with the natural frequency as the filtering cutoff frequency, the components obtained by performing wavelet decomposition and the residual components are respectively subjected to FFT conversion, combining the natural frequency f0And determining a filtering component so as to obtain a reconstruction signal of the balance and a reconstruction signal of the unidirectional accelerometer.
Taking the db10 wavelet basis as an example of 10-time decomposition, for any direction, the force measurement data s output by the balance is decomposed for 10 times to obtain 10 components, and the relationship from high frequency to low frequency is D1～D10Force measurement data s output by the unidirectional accelerometerJSDAfter 10 decompositions, 10 components are obtained, corresponding to D from high frequency to low frequencyJSD-1～DJSD-10Together with the residual component (assuming that the residual component of the force measurement data s output by the balance is A)10And force measurement data s output by the unidirectional accelerometerJSDThe remaining component of (A)JSD-10) Respectively, for a total of 11 components. FFT conversion is respectively carried out on the respective 11 components, and the filter cut-off frequency is determined according to the natural frequency f0Determining, combining FFT transforms and f0Determining the filter component to obtain the reconstruction signals of the balance (corresponding to the direction) and the one-way accelerometer as s respectivelyWTAnd sWT-JSDNamely:
wherein the component D of the force-measuring data s output by the balanceMAnd force measurement data s output by the unidirectional accelerometerJSDComponent D ofJSD-HAre not greater than the natural frequency f0。
In step S5, when the dynamic characteristic compensation is performed on the force measurement data output by the balance according to the corresponding dynamic characteristic compensation coefficient, for any direction, the compensated balance signal is equal to the product of the dynamic characteristic compensation coefficient subtracted from the reconstructed signal of the balance obtained after the high-frequency filtering and the reconstructed signal of the unidirectional accelerometer, that is, Send＝sWT-ty×sWT-JSD。
In some preferred embodiments, the present invention further provides a compensation apparatus for compensating for a characteristic of a celestial motion based on wavelet reconstruction, the apparatus including: three one-way accelerometers, a coefficient unit, a signal unit and a calculation unit, wherein:
the unidirectional accelerometer is respectively arranged in the axial direction, the normal direction and the lateral direction of the model, and the frequency response of the unidirectional accelerometer is greater than the integral natural frequency of the model, the balance and the support rod.
And the coefficient unit is connected with the balance and each unidirectional accelerometer and is used for acquiring the balance and force measurement signals output by each unidirectional accelerometer in the process of simulating the impact on a model in a wind tunnel force measurement test under a windless condition and respectively calculating the dynamic characteristic compensation coefficients in three directions.
The signal unit is connected with the coefficient unit, the balance and the unidirectional accelerometers, and is used for acquiring the force measurement data output by the balance and the unidirectional accelerometers when a wind tunnel force measurement test is carried out, respectively carrying out wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometers in three directions, carrying out high-frequency filtering by taking the inherent frequency as a filtering cut-off frequency, and carrying out dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient to obtain a compensated balance signal;
and the calculation unit is connected with the signal unit and is used for calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
Preferably, when the coefficient unit calculates the dynamic characteristic compensation coefficients in three directions, the following formula is used for each direction:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance, N representing the total number of points of the force measurement signal participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point of the force measuring signal output by the balance, YiAnd the voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer is represented.
Preferably, when the signal unit acquires force measurement data output by the balance and the unidirectional accelerometers, the signal unit records original data output by the balance and the unidirectional accelerometers in three directions respectively, determines an effective data segment and zero data in the original data, and deducts zero data from the effective data segment, wherein the zero data is obtained from a zero state of stable damped oscillation after the model is impacted.
Preferably, when the signal unit carries out wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer in three directions, the decomposition is carried out by adopting a db10 wavelet base. Further, when the db10 wavelet basis is adopted for decomposition, the decomposition is carried out 5-10 times, and preferably 10 times.
Preferably, when the signal unit performs dynamic characteristic compensation on the force measurement data output by the balance according to the corresponding dynamic characteristic compensation coefficient, the compensated balance signal s in any direction is obtainedend＝sWT-ty×sWT-JSDWherein s isWTRepresenting the reconstructed signal, s, of the balance after high-frequency filteringWT-JSDAnd representing a reconstructed signal of the unidirectional accelerometer obtained after high-frequency filtering.
In summary, the invention provides a compensation method and a compensation device for the astronomical dynamic characteristics based on wavelet reconstruction. Meanwhile, the invention utilizes the characteristic that wavelet reconstruction has multi-resolution, carries out high-frequency filtering through wavelet decomposition, has high calculation speed, can better observe the local characteristics of signals and is beneficial to reducing the interference of high-frequency signals.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A skyward dynamic characteristic compensation method based on wavelet reconstruction is characterized by comprising the following steps:
s1, installing the model, the balance and the support rod, and acquiring the integral natural frequency of the model, the balance and the support rod;
s2, respectively installing unidirectional accelerometers in the axial direction, the normal direction and the lateral direction of the model, wherein the frequency response of the unidirectional accelerometers is greater than the natural frequency;
s3, under the windless condition, simulating the impact process of the model in a wind tunnel force measurement test, acquiring force measurement signals output by the balance and the unidirectional accelerometers, and respectively calculating dynamic characteristic compensation coefficients in three directions;
s4, performing a wind tunnel force measurement test, impacting the model, and acquiring force measurement data output by the balance and each one-way accelerometer;
s5, performing wavelet decomposition on the force measurement data output by the balance and the unidirectional accelerometer in three directions, performing high-frequency filtering by taking the inherent frequency as a filtering cut-off frequency, and performing dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient to obtain a compensated balance signal;
and S6, calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
2. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 1, wherein: in step S3, when the dynamic characteristic compensation coefficients in the three directions are calculated, the following formula is used for each direction:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance, N representing the total number of points of the force measurement signal participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point of the force measuring signal output by the balance, YiAnd the voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer is represented.
3. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 2, wherein: in step S4, when the force measurement data output by the balance and the unidirectional accelerometers are obtained, the original data output by the balance and the unidirectional accelerometers in three directions are recorded, an effective data segment and zero point data in the original data are determined, and zero point data is subtracted from the effective data segment, where the zero point data is obtained from a zero point state of stable damped oscillation after the model is impacted.
4. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 2, wherein: in the step S5, when the force measurement data output by the balance and the unidirectional accelerometer are respectively subjected to wavelet decomposition in three directions, a db10 wavelet basis is used for decomposition.
5. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 4, wherein: in the step S5, decomposition is performed 5-10 times when db10 wavelet basis is adopted for decomposition.
6. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 4, wherein: in step S5, when the dynamic characteristic compensation is performed on the force measurement data output by the balance according to the corresponding dynamic characteristic compensation coefficient, the compensated balance signal S in any direction is obtainedend＝sWT-ty×sWT-JSDWherein s isWTRepresenting the reconstructed signal, s, of the balance after high-frequency filteringWT-JSDAnd representing a reconstructed signal of the unidirectional accelerometer obtained after high-frequency filtering.
7. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 1, wherein: in the step S1, when the natural frequency of the model, the balance, and the whole support rod is obtained, a manual tapping mode is adopted, which includes the following steps:
knocking the center of mass of the model, the balance and the support rod by a force hammer, and recording the force measurement data output by the balance; and determining the integral natural frequency of the model, the balance and the support rod by performing frequency spectrum analysis on the force measurement data output by the balance.
8. The wavelet reconstruction-based astronomical motion dynamic characteristic compensation method according to claim 1, wherein: in the step S3, when an impact process on the model in the wind tunnel force measurement test is simulated, the center of mass of the model, the balance and the support rod is knocked by the force hammer, and the direction in which the force hammer knocks the model is consistent with the direction in which the model is knocked when the wind tunnel force measurement test is performed in the step S4.
9. A compensation device for the characteristic of a astronomical motion state based on wavelet reconstruction is characterized by comprising: the system comprises three unidirectional accelerometers, a coefficient unit, a signal unit and a calculation unit; wherein the content of the first and second substances,
the unidirectional accelerometers are respectively arranged in the axial direction, the normal direction and the lateral direction of the model, and the frequency response of the unidirectional accelerometers is greater than the integral natural frequency of the model, the balance and the support rod;
the coefficient unit is connected with the balance and each one-way accelerometer and is used for acquiring force measurement signals output by the balance and each one-way accelerometer in the process of simulating the impact on the model in a wind tunnel force measurement test under a windless condition and respectively calculating dynamic characteristic compensation coefficients in three directions;
the signal unit is connected with the coefficient unit, the balance and each one-way accelerometer and is used for acquiring force measurement data output by the balance and each one-way accelerometer during a wind tunnel force measurement test, performing wavelet decomposition on the force measurement data output by the balance and each one-way accelerometer in three directions, performing high-frequency filtering by taking inherent frequency as filtering cut-off frequency, and performing dynamic characteristic compensation on the force measurement data output by the balance according to a corresponding dynamic characteristic compensation coefficient to acquire a compensated balance signal;
and the calculation unit is connected with the signal unit and is used for calculating the aerodynamic force of the wind tunnel force measurement test according to the compensated balance signal and the aerodynamic force formula.
10. The wavelet reconstruction-based astronomical dynamic characteristic compensation apparatus according to claim 9, wherein: when the coefficient unit calculates the dynamic characteristic compensation coefficients in three directions, the following formula is adopted for calculation in each direction:
wherein the content of the first and second substances,tyrepresenting the compensation coefficient of the dynamic characteristic of the balance, N representing the total number of points of the force measurement signal participating in the calculation, i ∈ [1, N]，DiIndicating the voltage value at the ith point of the force measuring signal output by the balance, YiAnd the voltage value of the ith point of the force measurement signal output by the unidirectional accelerometer is represented.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010546766.5A CN111649908B (en) | 2020-06-16 | 2020-06-16 | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010546766.5A CN111649908B (en) | 2020-06-16 | 2020-06-16 | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111649908A true CN111649908A (en) | 2020-09-11 |
| CN111649908B CN111649908B (en) | 2022-03-01 |
Family
ID=72347216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010546766.5A Active CN111649908B (en) | 2020-06-16 | 2020-06-16 | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111649908B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112345199A (en) * | 2020-10-29 | 2021-02-09 | 中国空气动力研究与发展中心高速空气动力研究所 | Temporary-impulse high-speed wind tunnel attack angle sensor vibration influence correction method |
| CN112525471A (en) * | 2020-12-07 | 2021-03-19 | 中国空气动力研究与发展中心高速空气动力研究所 | Rolling damping compensation method for wind tunnel free rolling dynamic experiment |
| CN112697387A (en) * | 2020-12-23 | 2021-04-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method for analyzing validity of measurement data of film resistance thermometer in wind tunnel aerodynamic heat test |
| CN112729764A (en) * | 2020-12-14 | 2021-04-30 | 华南理工大学 | High-frequency base balance power calibration method and system for underdetermined condition |
| CN117909659A (en) * | 2024-03-20 | 2024-04-19 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method and device for processing aerodynamic force signals of balance in wind tunnel dynamic test of variant aircraft |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103940577A (en) * | 2014-04-21 | 2014-07-23 | 西南交通大学 | Wind tunnel test balance assessment method based on acceleration signal energy |
| CN108195555A (en) * | 2018-01-16 | 2018-06-22 | 中国空气动力研究与发展中心超高速空气动力研究所 | Optical fibre balance aerodynamics force measurement system and measuring method |
| CN108344556A (en) * | 2018-04-28 | 2018-07-31 | 中国空气动力研究与发展中心超高速空气动力研究所 | A kind of hypersonic wind tunnel balance dynamic characteristic verifying attachment |
| CN109359340A (en) * | 2018-09-18 | 2019-02-19 | 中国科学院力学研究所 | The measurement method and device of six component aerodynamic force of high-speed rain movable model type |
| CN110132527A (en) * | 2019-06-24 | 2019-08-16 | 中国空气动力研究与发展中心高速空气动力研究所 | A kind of model vibration monitoring method in the wind tunnel test based on balance signal |
-
2020
- 2020-06-16 CN CN202010546766.5A patent/CN111649908B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103940577A (en) * | 2014-04-21 | 2014-07-23 | 西南交通大学 | Wind tunnel test balance assessment method based on acceleration signal energy |
| CN108195555A (en) * | 2018-01-16 | 2018-06-22 | 中国空气动力研究与发展中心超高速空气动力研究所 | Optical fibre balance aerodynamics force measurement system and measuring method |
| CN108344556A (en) * | 2018-04-28 | 2018-07-31 | 中国空气动力研究与发展中心超高速空气动力研究所 | A kind of hypersonic wind tunnel balance dynamic characteristic verifying attachment |
| CN109359340A (en) * | 2018-09-18 | 2019-02-19 | 中国科学院力学研究所 | The measurement method and device of six component aerodynamic force of high-speed rain movable model type |
| CN110132527A (en) * | 2019-06-24 | 2019-08-16 | 中国空气动力研究与发展中心高速空气动力研究所 | A kind of model vibration monitoring method in the wind tunnel test based on balance signal |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112345199A (en) * | 2020-10-29 | 2021-02-09 | 中国空气动力研究与发展中心高速空气动力研究所 | Temporary-impulse high-speed wind tunnel attack angle sensor vibration influence correction method |
| CN112345199B (en) * | 2020-10-29 | 2022-07-22 | 中国空气动力研究与发展中心高速空气动力研究所 | Method for correcting impact of vibration of attack angle sensor of temporary-impulse high-speed wind tunnel |
| CN112525471A (en) * | 2020-12-07 | 2021-03-19 | 中国空气动力研究与发展中心高速空气动力研究所 | Rolling damping compensation method for wind tunnel free rolling dynamic experiment |
| CN112525471B (en) * | 2020-12-07 | 2022-03-25 | 中国空气动力研究与发展中心高速空气动力研究所 | Rolling damping compensation method for wind tunnel free rolling dynamic experiment |
| CN112729764A (en) * | 2020-12-14 | 2021-04-30 | 华南理工大学 | High-frequency base balance power calibration method and system for underdetermined condition |
| CN112729764B (en) * | 2020-12-14 | 2022-03-29 | 华南理工大学 | High-frequency base balance power calibration method and system for underdetermined condition |
| CN112697387A (en) * | 2020-12-23 | 2021-04-23 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method for analyzing validity of measurement data of film resistance thermometer in wind tunnel aerodynamic heat test |
| CN117909659A (en) * | 2024-03-20 | 2024-04-19 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method and device for processing aerodynamic force signals of balance in wind tunnel dynamic test of variant aircraft |
| CN117909659B (en) * | 2024-03-20 | 2024-05-14 | 中国空气动力研究与发展中心超高速空气动力研究所 | Method and device for processing aerodynamic force signals of balance in wind tunnel dynamic test of variant aircraft |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111649908B (en) | 2022-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111649908B (en) | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction | |
| US11372129B2 (en) | Post-compensation method for motion errors of rotating accelerometer gravity gradiometer | |
| CN111157208A (en) | Satellite micro-vibration isolation simulation measurement system and method | |
| CN209247819U (en) | It is a kind of for calibrating the impact test apparatus of acceleration transducer | |
| CN114964673B (en) | Structural frequency response function correction method for frequency spectrum leakage error | |
| CN106404914B (en) | The structural damage of a kind of pair of Wooden Pagoda at Yingxian County and the test method of safe condition | |
| CN112345199B (en) | Method for correcting impact of vibration of attack angle sensor of temporary-impulse high-speed wind tunnel | |
| CN113740929B (en) | Ground testing device and method for space inertial sensor | |
| CN112881743B (en) | Speed and acceleration calculation method and measurement device based on regularization algorithm | |
| CN105509867B (en) | It is a kind of to measure the micro- angular oscillation method of satellite structure using optical fibre gyro | |
| CN109633205B (en) | Temperature compensation method for quartz resonance accelerometer | |
| CN108573084B (en) | Environmental vibration test method and system | |
| CN109060241B (en) | Sensor equivalent inertia force correction method for pulse force calibration | |
| CN114577318B (en) | Vehicle-mounted weighing module and sensing method thereof | |
| JPH09145530A (en) | Fatigue evaluation apparatus for piping system | |
| CN115628801A (en) | Laser measuring device based on DataSocket communication | |
| CN109359340A (en) | The measurement method and device of six component aerodynamic force of high-speed rain movable model type | |
| CN109871658A (en) | The multi-pose optimal estimation method measured for guided missile warhead rotary inertia and the product of inertia | |
| CN103712794B (en) | Gear train assembly impact shock energy slippages rate detection system and method | |
| CN116839780B (en) | Weak electromagnetic force measuring device and method for wireless power transmission coupler | |
| CN218444386U (en) | Passive dynamic vibration absorber natural frequency testing and debugging device | |
| RU2743778C1 (en) | Method for determining the non-stationary force and a device for its implementation | |
| CN110569568A (en) | Method for eliminating suspension rigidity and sensor quality influence in free modal test | |
| Bruns et al. | Traceable calibration for the dynamic | |
| JP3335765B2 (en) | Vibration characteristic analyzer |
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 |