CN114964673A - 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
- CN114964673A CN114964673A CN202210376662.3A CN202210376662A CN114964673A CN 114964673 A CN114964673 A CN 114964673A CN 202210376662 A CN202210376662 A CN 202210376662A CN 114964673 A CN114964673 A CN 114964673A
- 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.)
- Granted
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)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Computing Systems (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 include 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 proposes a method for weighted overlap-and-average data, which proves that the leakage error and random error can be reduced by the overlap ratio of adjacent data blocks. When the modal parameter of the engineering structure is identified by actually measuring the frequency response function by 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 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 Respectively representing undetermined coefficients in a numerator and a denominator of a theoretical displacement frequency response function, wherein k is 0,1, … and 2N-1; j represents an imaginary number sheetBit and satisfy 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 mode, i ═ 1,2, …, 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 beneficial results of the invention are as follows: 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 is 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 pulse 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 )]Where 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 dot division, i.e. 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, A k (k-0, 1, …,2N-1) and B k (k is 0,1, …,2N) represents the undetermined coefficients in the numerator and denominator of the theoretical shift response function, j represents the 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 the ith 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 is 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 aiming at a frequency spectrum leakage error is characterized in that firstly, the finite record length displacement response and the excitation force response of a structure under pulse excitation are respectively subjected to fast Fourier transform, and discrete measured 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 fractional 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:
n represents the number of peak values in an amplitude diagram of the actually measured displacement frequency response data; a. the k And B k Respectively representing undetermined coefficients in a numerator and a denominator of a theoretical displacement frequency response function, wherein k is 0,1, … and 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 mode, i ═ 1,2, …, 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 notation of the displacement responseThe length of the recording is determined,
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 true CN114964673A (en) | 2022-08-30 |
| CN114964673B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116299147A (en) * | 2023-03-13 | 2023-06-23 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
| CN118094197A (en) * | 2024-04-23 | 2024-05-28 | 南京航空航天大学 | Structural modal parameter identification method based on optical fiber strain sensing and fitting method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105953996A (en) * | 2016-06-30 | 2016-09-21 | 东南大学 | Bridge detection and evaluation method and equipment based on impact vibration |
| CN106844935A (en) * | 2017-01-18 | 2017-06-13 | 大连理工大学 | A kind of big damping engineering structure Modal Parameters Identification |
| CN106960068A (en) * | 2016-09-30 | 2017-07-18 | 中国人民解放军海军工程大学 | A kind of damping ratios quick calculation method based on pulse excitation response spectrum |
| CN109060284A (en) * | 2018-08-07 | 2018-12-21 | 广东工业大学 | A kind of Experimental modal 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
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105953996A (en) * | 2016-06-30 | 2016-09-21 | 东南大学 | Bridge detection and evaluation method and equipment based on impact vibration |
| CN106960068A (en) * | 2016-09-30 | 2017-07-18 | 中国人民解放军海军工程大学 | A kind of damping ratios quick calculation method based on pulse excitation response spectrum |
| CN106844935A (en) * | 2017-01-18 | 2017-06-13 | 大连理工大学 | A kind of big damping engineering structure Modal Parameters Identification |
| CN109060284A (en) * | 2018-08-07 | 2018-12-21 | 广东工业大学 | A kind of Experimental modal 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 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116299147A (en) * | 2023-03-13 | 2023-06-23 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
| CN116299147B (en) * | 2023-03-13 | 2023-11-28 | 中国科学院声学研究所 | One-dimensional structure internal sound source positioning method based on acoustic coherence technology |
| CN118094197A (en) * | 2024-04-23 | 2024-05-28 | 南京航空航天大学 | Structural modal parameter identification method based on optical fiber strain sensing and fitting method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114964673B (en) | 2023-02-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114964673B (en) | Structural frequency response function correction method for frequency spectrum leakage error | |
| CN107590317B (en) | Generator dynamic estimation method considering model parameter uncertainty | |
| Eichstädt et al. | Deconvolution filters for the analysis of dynamic measurement processes: a tutorial | |
| CN107608935B (en) | Time rearrangement compression transformation-based impact signal time-frequency analysis and reconstruction method | |
| CN106487359B (en) | The building method of Morphologic filters based on self-adapting multi-dimension AVG-Hat transformation | |
| CN102759448B (en) | Gearbox fault detection method based on flexible time-domain averaging | |
| CN111649908B (en) | Heaven-horizontal dynamic characteristic compensation method and device based on wavelet reconstruction | |
| CN109902408B (en) | Load identification method based on numerical operation and improved regularization algorithm | |
| CN113804404B (en) | Light source sweep frequency nonlinear correction method for optical frequency domain polarization crosstalk measurement | |
| CN112362343A (en) | Distributed fault feature extraction method for gearbox under variable rotating speed based on frequency modulation dictionary | |
| CN113820003A (en) | Acceleration real-time reconstruction dynamic displacement method suitable for bridge vibration monitoring | |
| CN110263762B (en) | Output-based semi-submersible ocean platform energy transfer path analysis method | |
| CN105938508A (en) | Method for accurately calculating frequency and amplitude of vibration or pressure fluctuation signal | |
| CN106326530A (en) | Time-varying structure modal parameter identification method based on right matrix fraction model | |
| JP4509837B2 (en) | Early earthquake specifications estimation method and system | |
| 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 | |
| CN117147022B (en) | Force sensor nonlinear compensation method and system | |
| CN107561934B (en) | Output-only structure modal parameter identification method based on multi-reference time domain transfer rate | |
| CN111125613A (en) | Method for improving noise-resistant capacity of Duffing chaotic oscillator for detecting weak resonant signal | |
| CN111353415B (en) | Method for detecting harmonic component in impulse response | |
| CN111062073B (en) | Maritime structure Laplace domain dynamic response transient separation method | |
| CN108801296B (en) | Sensor frequency response function calculation method based on error model iterative compensation | |
| CN107732940B (en) | Power system stabilizer parameter optimization test method based on ADPSS | |
| CN108132399B (en) | Simplified interpolation method for improving electric energy quality analysis precision of digital substation |
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 |