Nonlinear system dynamic parameter identification method based on high-order dynamic modal decomposition

Abstract

The invention relates to the field of nonlinear system dynamic parameter identification, in particular to a nonlinear system dynamic parameter identification method based on high-order dynamic modal decomposition. The scheme comprises the following steps: obtaining a time sequence of a nonlinear system through test or numerical simulation; determining the optimal time delay amount of the time sequence according to a mutual information method; reconstructing an original time sequence into a high-order time sequence matrix through a phase space according to the Tarkens embedding theorem; obtaining a similar matrix of the high-order time sequence system matrix; obtaining a high-order dynamic mode of the nonlinear system; judging whether the high-order dynamic modal characteristic spectrum is located in or close to a unit circle in a complex plane, and if not, returning to the step of reconstructing the original time sequence into the high-order time sequence matrix through the phase space; and if yes, obtaining modal frequency and damping ratio parameters of the nonlinear system. By introducing a mutual information method and a phase-space reconstruction theory, the analysis process in the prior art is optimized, and the method is suitable for nonlinear system dynamic parameter identification.

Description

Technical Field

[0001] The present invention relates to the field of dynamic parameter identification of nonlinear systems, and particularly to a method for identifying dynamic parameters of nonlinear systems based on high-order dynamic mode decomposition. Background Art

[0002] In the field of nonlinear systems, nonlinear systems are often accompanied by complex random evolution phenomena that are difficult to understand. If time series can be used as samples and the essential characteristics, related action mechanisms, and spatio-temporal evolution laws of these complex phenomena can be intuitively displayed through a certain method or technical analysis, it is of great significance for studying the transient changes and dynamic behaviors of nonlinear systems, and can promote the analysis and calculation efficiency of nonlinear system dynamics. Therefore, a modal decomposition method based on feature extraction technology - the flow field reduced-order model has been proposed by researchers. Its essence is a data-driven technology, which decomposes a high-dimensional unsteady system into the superposition of dynamic modes or coherent structures on a low-dimensional coordinate system, so as to describe the spatio-temporal evolution of nonlinear systems in a low-dimensional space. The currently more commonly used flow field reduced-order method is the proper orthogonal decomposition method.

[0003] The proper orthogonal decomposition technology decomposes a nonlinear system into several spatially orthogonal modes, which is a process of converting a multi-dimensional isotropic random field into a set of uncorrelated one-dimensional spatial patterns. These patterns consist of a series of continuous random processes and are sorted according to the energy (i.e., eigenvalues) of each mode, so as to select the main modes of the random field.

[0004] Although the application scope of the proper orthogonal decomposition technology is very wide, due to the existence of the covariance matrix, the proper orthogonal decomposition analysis is only limited to the second-order characteristics of variables, and it cannot directly identify single-frequency dynamic coherent structures to explain the time evolution characteristics of random systems.

[0005] In addition, a nonlinear system is a random, complex high-dimensional dynamic system. The acquisition of the dynamic information of the system often directly depends on the time series of multiple variables. However, if the data is not fully extended to the dimension of the original dynamic system, direct modal analysis may not be able to discover the fuzzy dynamic characteristics hidden in the original data. Summary of the Invention

[0006] The purpose of the present invention is to overcome the shortcomings of the prior art and provide a method for identifying dynamic parameters of nonlinear systems based on high-order dynamic mode decomposition, which can well identify modal parameters such as the natural frequency and damping ratio of dynamic structures, and can separate nonlinear coupling signals to obtain single-frequency signals, having a good decoupling effect.

[0007] The present invention adopts the following technical solutions to achieve the above purpose. The present invention provides a method for identifying dynamic parameters of nonlinear systems based on high-order dynamic mode decomposition, including: S1. Obtain the time series of the nonlinear system through experiments or numerical simulations; S2. Determine the optimal time delay of the time series according to the mutual information method; S3. According to the Takens embedding theorem, reconstruct the original time series into a high-order time series matrix through phase space; S4. Obtain the similarity matrix of the high-order time series system matrix; S5. Obtain the high-order dynamic modes of the nonlinear system; S6. Judge whether the characteristic spectrum of the high-order dynamic mode is located on or close to the unit circle in the complex plane. If so, go to step S7; otherwise, return to step S3; S7. Obtain the modal frequency and damping ratio parameters of the nonlinear system.

[0008] Furthermore, step S2 specifically includes: Let the random variable X represent the wind pressure coefficient sequence K(t), and Y represent the lagged wind pressure coefficient sequence K(t + τ ), τ where represents the time delay. The specific process of the mutual information method is as follows: Among them, H ( X ) and H ( Y ) are the marginal entropies of X and Y respectively, H ( X , Y ) is the joint entropy of X and Y , MI ( X , Y ) and X ( Y ) are the mutual information coefficients of p and x respectively, p ( y ) and X ( Y ) are the marginal distributions of p and x respectively, y ( X , Y ) is the joint distribution of x, y and X and Y ,

[0009] Furthermore, step S3 specifically includes: According to the Takens embedding theorem, the time seriesK ( N ) is constructed as m a multi-dimensional high-order time series matrix Q : where q j represents the matrix column vector at the j th moment, and the embedding dimension m is determined by whether the modulus of the high-order dynamic mode eigenvalue spectrum is equal to or close to 1.

[0010] Furthermore, step S4 specifically includes: The time-shifted sequence matrices are respectively Q 1 ={ q 1 , q 2 ,..., q j-1} and Q 2 ={ q 2 , q 3 ,..., q j}. Based on the assumption of the Koopman operator theory, Q 1 and Q 2 have a linear mapping A , that is: , A represents the system matrix; By performing singular value decomposition on the matrix Q 1 , the orthogonal subspace of the similarity transformation can be obtained: where U and V are unitary matrices, Σ is a diagonal matrix, H represents the conjugate transpose; The similarity matrix à is obtained through the minimization problem of the Frobenius norm: The matrix à is A a similarity transformation of

[0011] Furthermore, step S5 specifically includes: No. j Higher order dynamic modes φ j is a complex mode, expressed as: ; The logarithmic mapping of the eigenvalues ​​is defined as the eigenspectrum s j : ; in, λ j , w j They are respectively j eigenvalues ​​and eigenvectors, Δ t The time interval representing the time series of a nonlinear system.

[0012] Furthermore, step S7 specifically includes: The solution of the structural dynamics characteristic spectrum of the damped nonlinear system is: In the formula ω j , ζ j are the natural frequencies and the corresponding modal damping ratios, respectively; Transform the solution of the structural dynamics eigenvalues ​​of the damped nonlinear system to obtain the modal frequencies f j and modal damping ratio ζ j : In the formula, || represents the modulus of the complex number, and Re() represents the real part of the complex number.

[0013] The beneficial effects of the present invention are: The present invention combines the mutual information method and the Takens embedding theorem with the dynamic modal decomposition method to form a high-order dynamic modal decomposition method, which reconstructs the phase space of the nonlinear system and expands the data dimension. It can mine the fuzzy dynamic characteristics hidden in the data set and make the dynamic modal decomposition mode obtained by decomposition more neutral and stable.

[0014] The present invention can well identify modal parameters such as the natural frequency and damping ratio of the dynamic structure through a high-order dynamic modal decomposition method, and can separate nonlinear coupled signals to obtain single-frequency signals, thus having a good decoupling effect.

[0015] The non - linear system reconstructed by the high - order dynamic mode decomposition of the present invention can better describe and fit the local characteristics of the original non - linear system. This is because the high - order dynamic mode decomposition directly reconstructs the non - linear system instead of reconstructing the energy field. The low - frequency modes decomposed by the dynamic mode decomposition method contain most of the pulsation energy, explaining the dominant frequency of the non - linear system. Therefore, the high - order dynamic mode decomposition method has more advantages in revealing the dynamic mechanism and characteristics of the random field. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] Figure 1 is a flowchart of a method for identifying dynamic parameters of a non - linear system based on high - order dynamic mode decomposition provided by the present invention; Figure 2 is a schematic diagram of the measuring point arrangement provided by the embodiment of the present invention; Figure 3 is a schematic diagram of the wave height response provided by the present invention; Figure 4 is a schematic diagram of the power spectrum provided by the present invention; Figure 5 is a graph of the mutual information coefficient varying with time delay provided by the present invention; Figure 6 is a schematic diagram of the high - order dynamic mode characteristic spectrum after phase - space reconstruction provided by the present invention; Figure 7 is a curve graph of the 1st - order dynamic mode coefficient and the number of snapshots provided by the present invention; Figure 8 is a curve graph of the power spectrum of the 1st - order mode coefficient and frequency provided by the present invention; Figure 9 is a curve graph of the mode shape and the degree - of - freedom serial number provided by the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0017] To make the objectives, 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 accompanying drawings in the embodiments of the present invention.

[0018] The present invention provides a method for identifying dynamic parameters of a non - linear system based on high - order dynamic mode decomposition, as Figure 1 shown, including: S1. Obtain the time series of the non - linear system through experiments or numerical simulations; S2. Determine the optimal time delay of the time series according to the mutual information method; Let the random variable X represent the wind pressure coefficient sequence K(t), and Y represent the lagged wind pressure coefficient sequence K(t + τ ), τ representing the time delay, then the specific process of the mutual information method is as follows: Among them, H ( X )、 H ( Y ) are respectively X 、 Y 's marginal entropy, H ( X , Y ) is X 、 Y 's joint entropy, MI ( X , Y ) are respectively X 、 Y 's mutual information coefficient, p ( x )、 p ( y ) are respectively X 、 Y 's marginal distribution, p ( x , y ) is X 、 Y 's joint distribution, x, y is the specific value of the random variable X 、 Y .

[0019] S3. According to the Takens embedding theorem, the original time series is reconstructed into a high-order time series matrix through phase space; According to the Takens embedding theorem, the time series K ( N ) is constructed into a m -dimensional high-order time series matrix Q : In the formula, q j represents the matrix column vector at the j th moment, and the embedding dimension m is determined by whether the modulus of the high-order dynamic mode characteristic spectrum is equal to or close to 1.

[0020] S4. Obtain the similarity matrix of the high-order time series system matrix; The time-shifted sequence matrices are respectively Q 1 ={ q 1 , q 2 ,..., q j-1} and Q 2 ={ q2 , q 3 ,... q j}, based on the assumptions of Koopman operator theory, Q 1 and Q 2 there exists a linear mapping A , that is: , A represents the system matrix; By performing singular value decomposition on the matrix Q 1 , the orthogonal subspace of the similarity transformation can be obtained: where U and V are unitary matrices, Σ is a diagonal matrix, H represents the conjugate transpose; The similarity matrix à is obtained through the minimization problem of the Frobenius norm: The matrix à is A a similarity transformation of

[0021] S5. Obtain the high-order dynamic modes of the nonlinear system; Step S5 specifically includes: The j th high-order dynamic mode φ j is a complex mode, expressed as: ; To facilitate the analysis of high-order dynamic modes, the logarithmic mapping of the eigenvalues of the damped nonlinear system is defined as the characteristic spectrum s j : ; where, λ j , w j are the j th eigenvalue and eigenvector of Ã, respectively.

[0022] S6. Determine whether the characteristic spectrum of the high-order dynamic mode is located on or close to the unit circle in the complex plane. If so, go to step S7; otherwise, return to step S3; S7. Obtain the modal frequency and damping ratio parameters of the nonlinear system.

[0023] The solution for calculating the structural dynamic characteristic spectrum of the damped nonlinear system is as follows: In the formula ω j and ζ j are the natural frequency and the corresponding modal damping ratio respectively; The solution of the structural dynamic characteristic values of the damped nonlinear system is deformed to obtain the modal frequency f j and the modal damping ratio ζ j : In the formula, || represents the modulus of a complex number, and Re() represents the real part of a complex number.

[0024] The following combines specific implementation cases and appendices Figures 2 to 9 to elaborate on the present invention in detail.

[0025] (1) Experiment In the embodiment of the present invention, a shaking table test of a tuned liquid damper with a grid is adopted for high-order dynamic modal method analysis. A tuned liquid damper is a dynamic vibration damper commonly used in the control of super high-rise buildings. The natural frequency, damping ratio, and vibration mode of the tuned liquid damper system are important dynamic parameters for the design of the tuned liquid damper.

[0026] Since the liquid sloshing of the tuned liquid damper will exhibit obvious nonlinear characteristics, the performance parameters are more obtained by the shaking table test method. In the embodiment of the present invention, a shaking table test of the tuned liquid damper based on colored noise excitation is adopted. Figure 2 shows the schematic diagram of the measuring point layout of the tuned liquid damper model provided by the present invention. The clear dimension L×B×H (length×width×liquid depth) of the tuned liquid damper model is 2.1m×0.64m×0.44m. The model scale ratio, frequency scale ratio, and damping ratio scale ratio are 1:10, √10:1, and 1:1 respectively. Two grids are arranged, and the consistency ratio S =0.55. Four digital wave gauges are used in the experiment to measure the wave height change of the liquid in the tuned liquid damper along the long side direction, the sampling frequency is 100Hz, and an acceleration sensor is used to measure the acceleration of the shaking table surface, and the sampling frequency is 25Hz. According to the wave theory, the frequency and modal damping ratio of the 1st order mode of the liquid sloshing of the rectangular tuned liquid damper with a grid are 0.46Hz and 0.0477 respectively.

[0027] Figure 3 And Figure 4 respectively represent the wave height response and its power spectrum of the system wave height displacement gauge. Figure 4 In frepresents the frequency, and the vertical axis S(f) represents the power spectrum. It can be seen that the time history curve and power spectrum of the liquid surface wave height response show obvious non-linear characteristics compared with the results of the linear system, and there are coupling phenomena of different degrees between different signals. The first-order frequency is 0.452 Hz, which is in good agreement with the theoretical value, indicating the mutual applicability of the test model and the theoretical formula.

[0028] (2)Construct a high-order time series matrix Through the mutual information analysis of the test data, the delay time for the phase space reconstruction of the data is obtained. Figure 5 It represents the schematic diagram of the time delay when the mutual information coefficient decays from 1 to zero for the first time. The sampling times corresponding to the mutual information coefficient decreasing to 0.05 are 9.02 s respectively. At this time, it is considered that the data are independent of each other and used as the delay time for reconstructing the high-order time series matrix.

[0029] (3)High-order dynamic modal analysis The high-order time series matrix is analyzed through steps S4 - S5, as Figure 6 shown. The characteristic spectrum is located on or close to the unit circle in the complex plane, indicating that the eigenvalues of the system matrix are stable or neutrally stable.

[0030] (4)The first-order modal frequency and damping ratio of the tuned liquid damper Through the analysis of step S7, as Figures 7 - 9 , the first-order frequency and damping ratio are obtained as 0.452 Hz and 0.0487 respectively, which are in good agreement with the results of the theoretical formula. It can be seen that the high-order dynamic modal method has good applicability for identifying the dynamic modal parameters of the non-linear tuned liquid damper system. In addition, the high-order dynamic modal coefficient curve is a single-frequency decay curve. Combining with the power spectrum curve of the time mode shows that the high-order dynamic modal decomposition technology can well separate the coupled signals and play a decoupling role.

[0031] The above is only the preferred embodiment of the present invention. It should be understood that the present invention is not limited to the form disclosed herein, should not be regarded as excluding other embodiments, but can be used in various other combinations, modifications and environments, and can be changed within the scope of the concept described herein through the above teachings or the technology or knowledge in related fields. And the changes and alterations made by those skilled in the art without departing from the spirit and scope of the present invention shall fall within the protection scope of the appended claims of the present invention.

Claims

1. A nonlinear system dynamic parameter identification method based on high-order dynamic modal decomposition, characterized in that: include: S1. Obtain the time series of nonlinear system through experiments or numerical simulation; S2, determining the optimal time delay of the time series according to the mutual information method; S3. According to Takkens embedding theorem, the original time series is reconstructed into a high-order time series matrix through phase space; S4, obtaining the similarity matrix of the high-order time series system matrix; S5. Obtaining high-order dynamic modes of nonlinear systems; S6, judging whether the high-order dynamic modal characteristic spectrum is located at or close to the unit circle in the complex plane, if so, proceeding to step S7, otherwise returning to step S3; S7. Obtain the modal frequency and damping ratio parameters of the nonlinear system.

2. The method for identifying dynamic parameters of a nonlinear system based on high-order dynamic modal decomposition according to claim 1, characterized in that: Step S2 specifically includes: Let the random variable X Represents the wind pressure coefficient sequence K ( t ), Y Represents the lagged wind pressure coefficient series K ( t + τ ), τ Represents the time delay, then the specific process of the mutual information method is as follows: in, H ( X ), H ( Y ) are respectively X , Y The marginal entropy of H ( X , Y )for X , Y The joint entropy of MI ( X , Y ) are respectively X , Y The mutual information coefficient of p ( x ), p ( y ) are respectively X , Y The marginal distribution of p ( x , y )for X , Y The joint distribution of x, y is a random variable X , Y The specific value of .

3. The method for identifying dynamic parameters of a nonlinear system based on high-order dynamic modal decomposition according to claim 1, characterized in that: Step S3 specifically includes: According to Takkens embedding theorem, the time series K ( N ) is constructed as m dimensional high-order time series matrix Q : In the formula q j Indicates j The matrix column vector at each moment, embedding dimension m It is determined by whether the mode of the characteristic spectrum of the high-order dynamic mode is equal to or close to 1.

4. The method for identifying dynamic parameters of a nonlinear system based on high-order dynamic modal decomposition according to claim 1, characterized in that: Step S4 specifically includes: The two time-shift sequence matrices are Q 1={ q 1, q 2, ..., q j-1 }and Q 2={ q 2, q 3, ..., q j }, based on the Koopman operator theory assumption, Q 1 and Q 2There exists a linear mapping A, namely: , A represents the system matrix; Through the matrix Q 1 Perform singular value decomposition to obtain the orthogonal subspace of similarity transformation: in U and V is a unitary matrix, Σ is a diagonal matrix, H represents conjugate transpose; The similarity matrix is ​​obtained by minimizing the Frobenius norm à : matrix à yes A Similarity transformation of .

5. The method for identifying dynamic parameters of a nonlinear system based on high-order dynamic modal decomposition according to claim 4 is characterized in that: Step S5 specifically includes: No. j Higher order dynamic modes φ j is a complex mode, expressed as: ; The logarithmic mapping of the eigenvalues ​​is defined as the eigenspectrum s j : ; in, λ j , w j They are respectively j eigenvalues ​​and eigenvectors, Δ t The time interval representing the time series of a nonlinear system.

6. The method for identifying dynamic parameters of a nonlinear system based on high-order dynamic modal decomposition according to claim 5, characterized in that: Step S7 specifically includes: The solution of the structural dynamics characteristic spectrum of the damped nonlinear system is: In the formula ω j , ζ j are the natural frequencies and the corresponding modal damping ratios, respectively; Transform the solution of the structural dynamics eigenvalues ​​of the damped nonlinear system to obtain the modal frequencies f j and modal damping ratio ζ j : In the formula, || represents the modulus of the complex number, and Re() represents the real part of the complex number.

Images