CN114964673B - Structural frequency response function correction method for frequency spectrum leakage error - Google Patents
Structural frequency response function correction method for frequency spectrum leakage error Download PDFInfo
- Publication number
- CN114964673B CN114964673B CN202210376662.3A CN202210376662A CN114964673B CN 114964673 B CN114964673 B CN 114964673B CN 202210376662 A CN202210376662 A CN 202210376662A CN 114964673 B CN114964673 B CN 114964673B
- Authority
- CN
- China
- Prior art keywords
- frequency response
- displacement
- response function
- frequency
- fitted
- 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.)
- Active
Links
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
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
- G01M7/025—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/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Discrete Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
The invention belongs to the field of modal identification of structural data analysis, and provides a structural frequency response function correction method for a frequency spectrum leakage error. And establishing a theoretical displacement frequency response function to be fitted by adopting a rational-fractal polynomial method, and performing curve fitting on the discrete actually-measured displacement frequency response data by adopting a least square method to finally obtain the fitted displacement frequency response function. And performing partial fractional expansion on the fitted displacement frequency response function, and then extracting an imaginary part of each order of modal residue to construct a modified structure displacement frequency response function. The method can directly eliminate the influence of the frequency spectrum leakage error on the structure displacement frequency response function, namely, the real displacement frequency response function of the structure is accurately obtained, the calculation process is clear and simple, and iteration is not needed.
Description
Technical Field
The invention belongs to the field of modal identification of structural data analysis, and relates to a frequency response function correction method of an engineering structure.
Background
The modal analysis method is one of important means for solving the dynamic characteristic design of the modern complex structure, and modal parameters reflect the dynamic characteristics of the structure and can be used for vibration control, model correction, structural state evaluation and the like, so that the modal parameter identification of the structure by using the measured data is very important.
The frequency response function reflects the transfer capability of the system to different input signals, and is a non-parametric model for describing the characteristics of the dynamic system. For any system, the stability of the system can be directly analyzed through a frequency response function, the system is comprehensively designed and corrected, a powerful tool is provided for application control theory analysis and solving of complex problems in engineering practice, and important positions are occupied in aspects of modal analysis, parameter identification, model correction, fault diagnosis and the like.
At present, frequency response function estimation methods are mature and have a plurality of types, and the commonly used frequency response estimation methods mainly comprise three methods, namely an H1 method, an H2 method and an Hv method. On this basis, the spectral leakage error is generally reduced by applying a windowing function to the time domain data, such as a rectangular window, a hamming window, an exponential window, and the like. Schoukens further compared the spectral leakage principle of rectangular windows and hamming windows, and indicated that hamming windows actually convert leakage errors of rectangular windows into interpolation errors. Antoni proposed a method of weighted overlap-average of data, demonstrating that leakage errors and random errors can be reduced by the overlap ratio of adjacent data blocks. When the modal parameters of the engineering structure are identified by actually measuring the frequency response function through a frequency domain method, the estimation precision of the frequency response function directly leads to the quality of the modal parameter identification effect. However, the above frequency response estimation method can only reduce the influence of the spectrum leakage error, thereby causing unstable modal parameter identification. Therefore, how to correct the structural frequency response function for the spectrum leakage error is crucial, which also has reference value for the practical application of many frequency domain identification methods.
Disclosure of Invention
The invention aims to provide a displacement frequency response function correction method for a proportional damping structure, and solves the problem that a fitted displacement frequency response function is inaccurate due to frequency spectrum leakage errors.
The technical scheme of the invention is as follows:
a structure frequency response function correction method aiming at a frequency spectrum leakage error comprises the steps of firstly respectively carrying out fast Fourier transform on a limited recording length displacement response and an excitation force signal of a structure under pulse excitation, and obtaining discrete measured displacement frequency response data through a frequency response calculation method. And establishing a theoretical displacement frequency response function to be fitted by adopting a rational fraction polynomial method, and performing curve fitting on the discrete actually-measured displacement frequency response data by adopting a least square method to finally obtain the fitted displacement frequency response function. Performing partial fractional expansion on the fitted displacement frequency response function, then extracting the imaginary part of each order of modal residue, calculating the corresponding spectrum leakage error, and finally constructing a corrected structure displacement frequency response function; the method comprises the following steps:
(1) Applying single-point pulse excitation to the structure, acquiring excitation force response by using a force sensor, acquiring displacement response at a measuring point of the structure by using a displacement sensor, and respectively transforming the excitation force response and the displacement response to a frequency domain by adopting fast Fourier transform to obtain a corresponding excitation force frequency spectrum and a corresponding displacement response frequency spectrum;
(2) Obtaining discrete actually-measured displacement frequency response data by utilizing the ratio of the displacement response frequency spectrum to the excitation force frequency spectrum at each frequency point;
(3) Adopting a rational numerator polynomial method to establish a theoretical displacement frequency response function to be fitted as follows:
wherein, N represents the number of peak values in the amplitude diagram of the actually measured displacement frequency response data; a. The k And B k Undetermined coefficients in numerator and denominator of a theoretical displacement frequency response function are respectively represented, k =0,1, \8230, 2N-1; j represents an imaginary unit and satisfies j 2 =-1;
(4) Performing curve fitting on actually measured displacement frequency response data by a least square method to obtain undetermined coefficient A k And B k Obtaining a fitted displacement frequency response function;
(5) And (3) carrying out partial fraction expansion on the fitted displacement frequency response function:
wherein, mu i Residue representing the ith mode, symbol "+" representing conjugate, s i A pole representing the ith order mode, i =1,2, \ 8230;, N;
(6) The frequency and damping ratio are found by:
w i =|s i |
wherein, w i And ζ i The ith order frequency and the damping ratio of the structure are respectively, | and Re (-) respectively represent the amplitude and the real part of the extracted data;
(7) The corrected residue is found by:
wherein T is the finite recording length of the displacement response,damped frequencies for order i, imag (·) represents the imaginary part of the extracted data;
(8) Establishing a corrected structure displacement frequency response function:
the invention has the beneficial effects that: the method can directly eliminate the influence of the frequency spectrum leakage error on the structure displacement frequency response function, namely, the real displacement frequency response function of the structure is accurately obtained, the calculation process is clear and simple, and iteration is not needed.
Detailed Description
The following further illustrates embodiments of the present invention in conjunction with the technical solutions.
Taking a 3-layer frame structure as an example, the mass matrix and stiffness matrix are as follows:
the damping matrix adopts a Rayleigh damping matrix, namely C =0.2M. And applying pulse excitation to the layer 1 frame, wherein the response signal is the displacement of the layer 2 frame. The specific implementation mode of the method is as follows:
(1) Obtaining a finite record length displacement response x of a structure T (t)=[x(t 1 ),x(t 2 ),…,x(t l )]And a pulsed excitation force signal f (t) = [ f (t) = 1 ),f(t 2 ),…,f(t l )]Respectively transforming the time domain displacement response and the force signal to a frequency domain by adopting fast Fourier transform, wherein the expression is X T (ω)=[X T (ω 1 ),X T (ω 2 ),…,X T (ω l )]And F (ω) = [ F (ω) 1 ),F(ω 2 ),…,F(ω l )]Wherein l represents the number of data points of the displacement response or excitation force, ω i Represents the ith frequency point;
(2) Obtaining discrete measured displacement frequency response data by using the ratio of the displacement response frequency spectrum to the excitation force frequency spectrum, wherein the expression is
H T (ω)=X T (ω)./F(ω)=[H T (ω 1 ),H T (ω 2 ),…,H T (ω l )]
Wherein, the symbol "/" represents a dot division, that is, each column of elements in the vector is divided separately;
(3) By adopting a rational fraction polynomial method, a theoretical displacement frequency response function to be fitted is established as
Wherein, the number of peak values in the actually measured displacement frequency response data amplitude diagram expressed by N is 3 k (k =0,1, \ 8230;, 2N-1) and B k (k =0,1, \8230;, 2N) represents undetermined coefficients in numerator and denominator of theoretical shift frequency response function, respectively, j represents an imaginary unit and satisfies j 2 =-1;
(4) Performing curve fitting on actually measured frequency response data by a least square method to obtain undetermined coefficient A k And B k Obtaining a fitted displacement frequency response function;
(5) Partial fractional expansion is carried out on the fitted displacement frequency response function
Wherein, mu i Residue representing the ith mode of the frequency response function, symbol "+" representing the conjugate, s i A pole representing an ith order mode;
(6) The frequency and damping ratio are found by:
ω i =|s i |
wherein, ω is i And ζ i I-th order frequency and damping ratio of the structure, respectively, | and Re (·) represent the amplitude and real part of the extracted data, respectively.
(7) The corrected residue number [ mu' = j [ -5.3020, -1.7656, -2.3914 ] was obtained by the following equation]×10 -2 :
Wherein T is the finite recording length of the displacement response,damped frequencies for order i, imag (·) represents the imaginary part of the extracted data;
(8) Establishing a modified structure displacement frequency response function
Claims (1)
1. A structural frequency response function correction method for a frequency spectrum leakage error is characterized in that firstly, fast Fourier transform is respectively carried out on finite record length displacement response and excitation force response of a structure under pulse excitation, and discrete actual measurement displacement frequency response data are obtained through a frequency response calculation method; establishing a theoretical displacement frequency response function to be fitted by adopting a rational fraction polynomial method, and performing curve fitting on discrete actually-measured displacement frequency response data by adopting a least square method to finally obtain a fitted displacement frequency response function; performing partial fractional expansion on the fitted displacement frequency response function, then extracting the imaginary part of each order of modal residue and calculating the corresponding spectrum leakage error, and finally constructing a corrected structure displacement frequency response function; the method comprises the following steps:
(1) Applying single-point pulse excitation to the structure, acquiring excitation force response by using a force sensor, acquiring displacement response at a measuring point of the structure by using a displacement sensor, and respectively transforming the excitation force response and the displacement response to a frequency domain by adopting fast Fourier transform to obtain a corresponding excitation force frequency spectrum and a corresponding displacement response frequency spectrum;
(2) Obtaining discrete actually-measured displacement frequency response data by utilizing the ratio of the displacement response frequency spectrum to the excitation force frequency spectrum at each frequency point;
(3) Adopting a rational numerator polynomial method to establish a theoretical displacement frequency response function to be fitted as follows:
wherein, N represents the number of peak values in the amplitude diagram of the actually measured displacement frequency response data; a. The k And B k Undetermined coefficients in numerator and denominator of a theoretical displacement frequency response function are respectively represented, k =0,1, \8230, 2N-1; j represents an imaginary unit and satisfies j 2 =-1;
(4) Performing curve fitting on the actually measured displacement frequency response data by a least square method to obtain an undetermined coefficient A k And B k Obtaining a fitted displacement frequency response function;
(5) And (3) carrying out partial fractional expansion on the fitted displacement frequency response function:
wherein, mu i Residue representing the ith mode, symbol "+" representing conjugate, s i A pole representing the ith order mode, i =1,2, \ 8230;, N;
(6) The frequency and damping ratio are found by:
w i =|s i |
wherein, w i And ζ i The ith order frequency and the damping ratio of the structure are respectively, | and Re (-) respectively represent the amplitude and the real part of the extracted data;
(7) The corrected residue is found by the following equation:
wherein T is the finite recording length of the displacement response,damped frequencies for order i, imag (·) representing the imaginary part of the extracted data;
(8) Establishing a modified structure displacement frequency response function:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210376662.3A CN114964673B (en) | 2022-04-12 | 2022-04-12 | Structural frequency response function correction method for frequency spectrum leakage error |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210376662.3A CN114964673B (en) | 2022-04-12 | 2022-04-12 | Structural frequency response function correction method for frequency spectrum leakage error |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114964673A CN114964673A (en) | 2022-08-30 |
CN114964673B true CN114964673B (en) | 2023-02-14 |
Family
ID=82976881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210376662.3A Active CN114964673B (en) | 2022-04-12 | 2022-04-12 | Structural frequency response function correction method for frequency spectrum leakage error |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114964673B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116299147B (en) * | 2023-03-13 | 2023-11-28 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105953996B (en) * | 2016-06-30 | 2017-04-19 | 东南大学 | Bridge detection and evaluation method and equipment based on impact vibration |
CN106960068B (en) * | 2016-09-30 | 2020-07-17 | 中国人民解放军海军工程大学 | Rapid modal damping ratio calculation method based on pulse excitation response frequency spectrum |
CN106844935B (en) * | 2017-01-18 | 2020-04-24 | 大连理工大学 | Large-damping engineering structure modal parameter identification method |
CN109060284B (en) * | 2018-08-07 | 2020-07-10 | 广东工业大学 | Test mode analysis method based on DIC technology |
CN113050596A (en) * | 2021-03-12 | 2021-06-29 | 北京强度环境研究所 | Method for accurately acquiring modal parameters of air rudder under random excitation |
-
2022
- 2022-04-12 CN CN202210376662.3A patent/CN114964673B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN114964673A (en) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Real-time dynamic displacement monitoring with double integration of acceleration based on recursive least squares method | |
Eichstädt et al. | Deconvolution filters for the analysis of dynamic measurement processes: a tutorial | |
CN114964673B (en) | Structural frequency response function correction method for frequency spectrum leakage error | |
CN101403774B (en) | Harmonic wave analysis method based on non-synchronous sampling | |
CN102759448B (en) | Gearbox fault detection method based on flexible time-domain averaging | |
CN113804404B (en) | Light source sweep frequency nonlinear correction method for optical frequency domain polarization crosstalk measurement | |
CN111649908B (en) | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction | |
CN104880172B (en) | Measuring flatness of road surface method and device based on Kalman filtering | |
CN108896274B (en) | Distributed optical fiber strain demodulation method based on subset window length optimization algorithm | |
CN109902408B (en) | Load identification method based on numerical operation and improved regularization algorithm | |
CN113820003A (en) | Acceleration real-time reconstruction dynamic displacement method suitable for bridge vibration monitoring | |
CN105043433A (en) | Neural-network-based rapid compensation method for photoelectric encoder | |
CN105938508A (en) | Method for accurately calculating frequency and amplitude of vibration or pressure fluctuation signal | |
CN110263762B (en) | Output-based semi-submersible ocean platform energy transfer path analysis method | |
CN104748704A (en) | Thin-walled structure ultrasonic resonance thickness measurement frequency spectrum analysis interpolation correction method | |
Matania et al. | Algorithms for spectrum background estimation of non-stationary signals | |
CN106872773A (en) | A kind of the multiple-pulse Precision Method of Freuqency Measurement and device of single carrier frequency pulse signal | |
CN111125613A (en) | Method for improving noise-resistant capacity of Duffing chaotic oscillator for detecting weak resonant signal | |
CN114280366A (en) | Sinusoidal signal frequency estimation method based on improved frequency interpolation algorithm | |
CN113156206A (en) | Time-frequency combined noise-containing signal parameter estimation new algorithm | |
CN107732940B (en) | Power system stabilizer parameter optimization test method based on ADPSS | |
CN109632071B (en) | Method and device for generating underwater acoustic environment noise data based on time-frequency characteristics | |
CN108801296B (en) | Sensor frequency response function calculation method based on error model iterative compensation | |
CN108132399B (en) | Simplified interpolation method for improving electric energy quality analysis precision of digital substation | |
CN112965455B (en) | Device and method for testing dynamic characteristics of actuator |
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 |