WO2018233359A1 - Method for recognizing spatial distribution and statistical feature of random dynamic load - Google Patents
Method for recognizing spatial distribution and statistical feature of random dynamic load Download PDFInfo
- Publication number
- WO2018233359A1 WO2018233359A1 PCT/CN2018/083275 CN2018083275W WO2018233359A1 WO 2018233359 A1 WO2018233359 A1 WO 2018233359A1 CN 2018083275 W CN2018083275 W CN 2018083275W WO 2018233359 A1 WO2018233359 A1 WO 2018233359A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- random
- modal
- dynamic load
- random dynamic
- space
- Prior art date
Links
Images
Classifications
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Definitions
- the invention relates to a method for identifying spatial distribution of random dynamic loads and statistical features, and belongs to the technical field of structural dynamics inverse problems.
- Dynamic load information on the engineering structure is the basis for structural design and safety assessment. Dynamic load acquisition methods are broadly divided into direct measurement and indirect recognition. In many cases, dynamic loads are difficult to obtain by direct measurement. Generally, the dynamic response of the structure is measured, and the known dynamic load information is identified when the structural system is known.
- the traditional dynamic load identification method uses the structural dynamic response data of a single actual measurement to identify the excitation information that causes the secondary dynamic response, and is a deterministic dynamic load identification method.
- the existing deterministic dynamic load identification method is used to obtain information such as concentrated dynamic load, moving load and distributed dynamic load on the engineering structure. It is worth noting that the distributed dynamic load identification problem is equivalent to identifying an infinite number of concentrated dynamic loads, which is more difficult. Generally, the distributed dynamic load identification problem needs to be reduced in dimension.
- the dynamic loads acting on the actual engineering structure are not only distributed on the structure, but also random.
- the dynamic response will also appear “randomness”; therefore, the structural dynamic response of a single measured measurement can only be one of the samples of the structural random dynamic response information, and the certainty is utilized based on a certain response sample.
- the dynamic load information obtained by the dynamic load identification method can only partially reflect the random dynamic load excitation; in addition, the dynamic response error contained in a single measurement is also used as part of the “true response” in the deterministic dynamic load identification, causing the load. Identify the deviation of the results.
- the traditional deterministic distributed dynamic load identification method and the identification method suitable for centralized random dynamic load can not be used. It is necessary to develop a new method for distributed random dynamic load identification. There are measured structures and dynamic response samples to identify the spatial distribution and statistical characteristics of random dynamic loads.
- the object of the present invention is to provide a method for identifying the spatial distribution of random dynamic loads and the identification of statistical features, and solving the problem of statistical characteristics and distribution characteristics of random dynamic loads on the structure of the dynamic response sample identification structure in the time domain, for serving in distributed randomization
- the engineering structure design and safety assessment under dynamic load environment provides accurate and reliable dynamic load information.
- a method for identifying a spatial distribution of random dynamic loads and a statistical feature characterized in that the method comprises the following steps:
- the structural random vibration response is developed by using the mode shape to obtain the dynamic response of the structure in the modal space;
- the rth measurement acquisition displacement vector sample W r at the position of the structure (x 1 , x 2 , ... x n ) is expressed as:
- w r (x j , t) represents the value of the structural displacement obtained at the rth measurement at time t at x j
- N is the number of measurements, and may also be considered as the response sample size
- the modal shape function is used to calculate the modal displacement vector corresponding to the rth measurement in the modal space:
- the random dynamic load distribution is independent of its random characteristics and time history, and its distribution function and random time history are all developed by using the mode shape.
- the specific steps are as follows: the distributed random dynamic load f(x, t, ⁇ ) is expressed.
- the product of its distribution function T(x) and the random time history P(t, ⁇ ) is as follows:
- the distribution function T(x) is expanded using the modal shape:
- the random dynamic load f i (t, ⁇ ) in the i-th modal space can be expressed as:
- L represents the length of the beam
- S51 The spatial distribution of the random dynamic load can be calculated by:
- t 1 is any time
- the invention has the following advantages:
- the existing random dynamic load identification technology can only identify the random concentrated dynamic load on the structure by the measured structure dynamic response sample.
- the existing distributed random dynamic load identification method cannot be applied to the identification of non-stationary random dynamic load.
- the distributed random dynamic load time domain identification technology provided by the invention can utilize the measured structure dynamic response sample at the limited measurement point to identify the spatial distribution of the random dynamic load and the variation of the statistical characteristics with time, and has certain advancement;
- the problem of identifying the distributed random dynamic load is transformed into the estimation problem of the dynamic process of the dynamic load in the modal space, which greatly reduces the dimension and difficulty of the load identification problem, and has easy operation and computational efficiency. High features.
- Figure 1 is a logic flow diagram of the method of the present invention.
- Figure 2 is a schematic diagram of a simply supported beam under distributed random dynamic loads.
- Figure 3 is a schematic diagram showing the results of spatial distribution identification of random dynamic loads.
- Fig. 4(a) is a schematic diagram showing the results of time-varying mean value identification of random dynamic loads.
- Fig. 4(b) is a schematic diagram showing the results of time-varying variance identification of random dynamic loads.
- Figure 4(c) is a schematic diagram showing the results of the random dynamic load correlation function.
- Embodiment The distributed random dynamic load acting on the simply supported beam structure as shown in Fig. 2 is identified by the method of the present invention.
- the trapezoidal distributed random dynamic load distribution function to be identified is:
- the stochastic dynamic load component F(t, ⁇ ) of distributed random dynamic load is divided into two parts: deterministic dynamic load and random dynamic load.
- the spatial distribution and statistical characteristics of the random dynamic load are identified by the measured dynamic random response sample by using the technique of the present invention, and specifically includes the following steps:
- S1 Perform modal test on the simply supported beam, and obtain the first five natural frequencies of the structure are 3.9 Hz, 15.6 Hz, 35.1 Hz, 62.5 Hz and 97.6 Hz, respectively, and obtain the mode shapes corresponding to the natural frequencies of each order;
- the rth measurement acquisition displacement response sample vector W r at each position on the beam structure is expressed as:
- w r (x j , t) represents the value of the structural displacement obtained at the rth measurement at time t at x j
- N is the number of measurements, which is 5000 times in this embodiment, that is, the total number of samples is 5000.
- the modal shape function is used to calculate the modal displacement vector corresponding to the rth measurement in the modal space:
- q i,r (t) is the modal displacement of the structural displacement obtained in the rth measurement in the i-th modal space
- the upper right corner + sign indicates the generalized inverse.
- S31 The random dynamic load distribution is independent of its random characteristics and time history, and its distribution function and random time history are all developed by using the mode shape. Specific steps are as follows:
- the random dynamic load f i (t, ⁇ ) in the i-th modal space can be expressed as:
- L represents the length of the beam.
- ⁇ A is the linear density of the beam.
- the identification method in the present invention can accurately identify the statistical characteristics of the random dynamic load with time by using the response samples at the limited measurement points, and is suitable for the case of non-stationary random dynamic loads, and has certain advancement.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
Description
Claims (5)
- 一种随机动载荷空间分布及统计特征的识别方法,其特征在于,该方法包括如下步骤:A method for identifying a spatial distribution of random dynamic loads and a statistical feature, characterized in that the method comprises the following steps:S1.开展模态试验,获取结构的模态参数,包括固有频率和模态振型;S1. Carry out a modal test to obtain the modal parameters of the structure, including the natural frequency and the mode shape;S2.将结构随机振动响应利用模态振型展开,获取结构在模态空间的动响应;S2. The structural random vibration response is developed by using the mode shape to obtain the dynamic response of the structure in the modal space;S3.利用模态振型展开,将随空间分布的随机动载荷投影到模态空间;S3. Using a mode shape expansion, projecting a random dynamic load distributed with space into a modal space;S4.在模态空间内,由随机动响应样本反演随机动载荷样本;S4. Inverting the random dynamic load sample from the random motion response sample in the modal space;S5.由模态空间内随机动载荷样本及振型函数求解结构上随机动载荷的空间分布和统计特征。S5. Solving the spatial distribution and statistical characteristics of random dynamic loads on a structure from random dynamic load samples and mode functions in modal space.
- 根据权利要求1所述的随机动载荷空间分布及统计特征的识别方法,其特征在于,步骤S2中所述的将结构随机振动响应利用模态振型展开,获取结构在模态空间的动响应,具体步骤为:The method for identifying a spatial distribution of random dynamic loads and a statistical feature according to claim 1, wherein the random vibration response of the structure is expanded by using a mode shape to obtain a dynamic response of the structure in a modal space. The specific steps are:S21:利用多次重复测量方式获取随机振动响应的样本集合:S21: Acquiring a sample set of random vibration responses by using multiple repeated measurements:在结构上(x 1,x 2,…x n)位置处第r次测量获取位移响应样本向量W r表示为: The rth measurement acquisition displacement vector sample W r at the position of the structure (x 1 , x 2 , ... x n ) is expressed as:W r={w r(x 1,t) w r(x 2,t) … w r(x n,t)} T,r=1,…,N (1), W r ={w r (x 1 ,t) w r (x 2 ,t) ... w r (x n ,t)} T , r=1,...,N (1),其中w r(x j,t)表示第r次测量获得的结构位移在x j处t时刻的值,N为测量的次数,也可认为响应样本容量; Where w r (x j , t) represents the value of the structural displacement obtained at the rth measurement at time t at x j , and N is the number of measurements, and may also be considered as the response sample size;S22:针对单个样本,即单次实测结构振动响应,利用模态振型展开获得结构振动在模态空间内响应:S22: For a single sample, that is, a single measured structural vibration response, the modal vibration mode is used to obtain the structural vibration response in the modal space:利用模态振型函数计算模态空间内第r次测量对应的模态位移向量:The modal shape function is used to calculate the modal displacement vector corresponding to the rth measurement in the modal space:其中q i,r(t)为第r次测量获得的结构位移在第i阶模态空间中的模态位移, 表示第i阶模态振型函数在x j处的值,右上角+号表示广义逆。 Where q i,r (t) is the modal displacement of the structural displacement obtained in the rth measurement in the i-th modal space, Indicates the value of the i-th mode mode shape function at x j , and the upper right corner + sign indicates the generalized inverse.
- 根据权利要求1所述的随机动载荷空间分布及统计特征的识别方法,其特征在于,步骤S3中所述的利用模态振型展开,将随空间分布的随机动载荷投影到模态空间,具体包括以下步骤:The method for identifying a spatial distribution of random dynamic loads and a statistical feature according to claim 1, wherein the modal vibration mode expansion is performed in step S3, and a random dynamic load distributed with space is projected into the modal space. Specifically, the following steps are included:S31:令随机动载荷分布独立于其随机特性和时间历程,将其分布函数和随机时间历程均利用模态振型展开,具体步骤如下:将分布随机动载荷f(x,t,θ)表示为其分布函数T(x)与随机 时间历程P(t,θ)的乘积,如下式所示:S31: The random dynamic load distribution is independent of its random characteristics and time history, and its distribution function and random time history are all developed by using the mode shape. The specific steps are as follows: the distributed random dynamic load f(x, t, θ) is expressed. The product of its distribution function T(x) and the random time history P(t, θ) is as follows:f(x,t,θ)=T(x)·P(t,θ) (3),f(x,t,θ)=T(x)·P(t,θ) (3),将分布函数T(x)利用模态振型展开:The distribution function T(x) is expanded using the modal shape:将式(4)代入式(3),可得:Substituting equation (4) into equation (3), you can get:其中a k(t,θ)=d kP(t,θ)。 Where a k (t, θ) = d k P(t, θ).第i阶模态空间内的随机动载荷f i(t,θ)可以表示为: The random dynamic load f i (t, θ) in the i-th modal space can be expressed as:其中L表示梁的长度;Where L represents the length of the beam;S32:利用模态振型函数的正交性,建立模态空间内随机动响应与随机动载荷之间的关系。具体步骤如下:S32: Using the orthogonality of the mode shape function, establishing a relationship between the random dynamic response and the random dynamic load in the modal space. Specific steps are as follows:结构在模态空间内的动力学方程为:The dynamic equations of the structure in the modal space are:其中q i(t,θ), 和 分别为模态空间内的结构随机位移,速度和加速度,ζ i和m i分别为第i阶模态阻尼比和模态质量。将式(6)代入式(7),利用结构模态振型的正交性条件,可以获得模态空间内随机动响应与随机动载荷之间的关系式,如下: Where q i (t, θ), with The random displacement, velocity and acceleration of the structure in the modal space, ζ i and m i are the i-th modal damping ratio and modal mass, respectively. Substituting equation (6) into equation (7), using the orthogonality condition of the structural mode shape, the relationship between the random dynamic response and the random dynamic load in the modal space can be obtained as follows:
- 根据权利要求1所述的随机动载荷空间分布及统计特征的识别方法,其特征在于,步骤S4中所述的在模态空间内,由随机动响应求解随机动载荷,具体包括以下步骤:The method for identifying a spatial distribution of random dynamic loads and a statistical feature according to claim 1, wherein in the modal space, the random dynamic response is solved by the random dynamic response in the modal space, which specifically includes the following steps:S41:根据结构位移,速度和加速度之间的导数关系,由第r次测量获得的结构模态位移q i,r(t),对时间求导获取对应的模态速度和模态加速度;若测量获得的是结构加速度信号,同样可采用积分的方式获取速度和位移; S41: According to the structural displacement, the derivative relationship between the velocity and the acceleration, the structural modal displacement q i,r (t) obtained by the rth measurement, and the time derivative to obtain the corresponding modal velocity and modal acceleration; The measurement obtains the structural acceleration signal, and the speed and displacement can also be obtained by integrating;S42:由式(8)中的模态空间内随机位移q i(t,θ)的样本q i,r(t)及其导数,求解得到模态空间内随机动载荷a i(t,θ)的样本a i,r(t)。 S42: Calculating the random dynamic load a i (t, θ in the modal space by the sample q i,r (t) and its derivative of the random displacement q i (t, θ) in the modal space in the equation (8) The sample a i,r (t).
- 根据权利要求1所述的随机动载荷空间分布及统计特征的识别方法,其特征在于,步骤S5中所述的由模态空间内随机动载荷样本及振型函数求解结构上随机动载荷的空间分布和统计特征,具体包括以下步骤:The method for identifying a random dynamic load spatial distribution and a statistical feature according to claim 1, wherein the random dynamic load sample and the mode function in the modal space are used to solve the space of the random dynamic load on the structure described in step S5. Distribution and statistical characteristics, including the following steps:S51:随机动载荷的空间分布可以由下式计算:S51: The spatial distribution of the random dynamic load can be calculated by:其中t 1为任意时刻; Where t 1 is any time;S52:随机动载荷的均值μ f(x,t)可以由下式计算: S52: The mean value of the random dynamic load μ f (x, t) can be calculated by:S53:随机动载荷的方差Var f(x,t)可以由下式计算: S53: The variance of the random dynamic load Var f (x, t) can be calculated by:
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710476795.7 | 2017-06-21 | ||
CN201710476795.7A CN107423487B (en) | 2017-06-21 | 2017-06-21 | A kind of recognition methods of Random dynamic loads spatial distribution and statistical nature |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2018233359A1 true WO2018233359A1 (en) | 2018-12-27 |
Family
ID=60426599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2018/083275 WO2018233359A1 (en) | 2017-06-21 | 2018-04-17 | Method for recognizing spatial distribution and statistical feature of random dynamic load |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN107423487B (en) |
WO (1) | WO2018233359A1 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107423487B (en) * | 2017-06-21 | 2018-06-12 | 东南大学 | A kind of recognition methods of Random dynamic loads spatial distribution and statistical nature |
CN108038315B (en) * | 2017-12-15 | 2019-04-09 | 东南大学 | A kind of Random dynamic loads recognition methods based on spectrum stochastic finite meta-model |
CN108491578B (en) * | 2018-02-13 | 2019-03-12 | 东南大学 | A kind of Random dynamic loads recognition methods based on perturbative matrix |
CN110083886A (en) * | 2019-04-08 | 2019-08-02 | 西北工业大学 | A kind of helicopter vibration load simulating analysis based on finite element |
CN113392547B (en) * | 2021-05-31 | 2024-01-05 | 东南大学 | Frequency domain distribution random dynamic load identification method considering spatial correlation |
CN113392545B (en) * | 2021-05-31 | 2024-01-30 | 东南大学 | Time domain distribution random dynamic load identification method considering spatial correlation |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129520A (en) * | 2011-03-16 | 2011-07-20 | 哈尔滨工程大学 | Self-adaptive time domain identification method of dynamic load |
US8577628B2 (en) * | 2009-04-10 | 2013-11-05 | University Of South Carolina | System and method for modal identification using smart mobile sensors |
CN104090977A (en) * | 2014-07-22 | 2014-10-08 | 东南大学 | Random recognition method for bridge floor moving vehicle loads |
CN104123463A (en) * | 2014-07-22 | 2014-10-29 | 东南大学 | Time domain identification method of random dynamic loads |
CN107423487A (en) * | 2017-06-21 | 2017-12-01 | 东南大学 | A kind of recognition methods of Random dynamic loads spatial distribution and statistical nature |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201104413D0 (en) * | 2011-03-16 | 2011-04-27 | Airbus Operations Ltd | Stresses induced by random loading |
-
2017
- 2017-06-21 CN CN201710476795.7A patent/CN107423487B/en active Active
-
2018
- 2018-04-17 WO PCT/CN2018/083275 patent/WO2018233359A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8577628B2 (en) * | 2009-04-10 | 2013-11-05 | University Of South Carolina | System and method for modal identification using smart mobile sensors |
CN102129520A (en) * | 2011-03-16 | 2011-07-20 | 哈尔滨工程大学 | Self-adaptive time domain identification method of dynamic load |
CN104090977A (en) * | 2014-07-22 | 2014-10-08 | 东南大学 | Random recognition method for bridge floor moving vehicle loads |
CN104123463A (en) * | 2014-07-22 | 2014-10-29 | 东南大学 | Time domain identification method of random dynamic loads |
CN107423487A (en) * | 2017-06-21 | 2017-12-01 | 东南大学 | A kind of recognition methods of Random dynamic loads spatial distribution and statistical nature |
Also Published As
Publication number | Publication date |
---|---|
CN107423487A (en) | 2017-12-01 |
CN107423487B (en) | 2018-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018233359A1 (en) | Method for recognizing spatial distribution and statistical feature of random dynamic load | |
US10228278B2 (en) | Determining a health condition of a structure | |
WO2018233360A1 (en) | Distributed random dynamic load recognition method based on kl expansion | |
WO2019011025A1 (en) | Finite element based method for identifying distributed random dynamic load | |
Neu et al. | Fully automated operational modal analysis using multi-stage clustering | |
US20160109323A1 (en) | Method for analyzing measured signal in resonance fatigue test and apparatus using the same | |
Kim et al. | Automated peak picking using region‐based convolutional neural network for operational modal analysis | |
Yao et al. | Blind modal identification using limited sensors through modified sparse component analysis by time‐frequency method | |
Hu et al. | Local damage detection of membranes based on Bayesian operational modal analysis and three-dimensional digital image correlation | |
Zhi et al. | Identification of wind loads on super-tall buildings by Kalman filter | |
Le et al. | Computationally efficient stochastic approach for the fragility analysis of vertical structures subjected to thunderstorm downburst winds | |
Wang et al. | Non-contact identification of rotating blade vibration | |
Khoshnoudian et al. | A new damage index using FRF data, 2D-PCA method and pattern recognition techniques | |
Wang et al. | Substructural identification of jack-up platform in time and frequency domains | |
Chen et al. | Gradient-based point tracking method and its application in the modal test of a solar array model | |
Rinn et al. | Stochastic method for in-situ damage analysis | |
CN110411387B (en) | Frequency-based beam type structure additional mass preliminary positioning method | |
Allada et al. | An unsupervised learning algorithm for computer vision-based blind modal parameters identification of output-only structures from video measurements | |
Asadollahi et al. | Statistical analysis of modal properties of a cable-stayed bridge through long-term structural health monitoring with wireless smart sensor networks | |
Dan et al. | Damping estimation by 2D dr_APES and its application to a real cable-stayed bridge | |
Hsiung et al. | Analyzing the dynamic response of rotating blades in small-scale wind turbines | |
Hoell et al. | Structural damage detection in wind turbine blades based on time series representations of dynamic responses | |
Chu et al. | Bayesian spectral density approach for identification and uncertainty quantification of bridge section's flutter derivatives operated in turbulent flow | |
Li et al. | Optimal sensor placement for structural health monitoring based on KL divergence | |
Ge et al. | Structural damage detection approach using input-output monitoring system and transfer learning |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 18819855 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 18819855 Country of ref document: EP Kind code of ref document: A1 |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 08.09.2020) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 18819855 Country of ref document: EP Kind code of ref document: A1 |