CN111781001A - Identification method of bridge damping ratio based on vehicle-bridge coupling - Google Patents
Identification method of bridge damping ratio based on vehicle-bridge coupling Download PDFInfo
- Publication number
- CN111781001A CN111781001A CN202010678822.0A CN202010678822A CN111781001A CN 111781001 A CN111781001 A CN 111781001A CN 202010678822 A CN202010678822 A CN 202010678822A CN 111781001 A CN111781001 A CN 111781001A
- Authority
- CN
- China
- Prior art keywords
- bridge
- damping ratio
- vehicle
- acceleration response
- identification method
- 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
- 238000013016 damping Methods 0.000 title claims abstract description 117
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000008878 coupling Effects 0.000 title claims abstract description 28
- 238000010168 coupling process Methods 0.000 title claims abstract description 28
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 28
- 238000012360 testing method Methods 0.000 claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 33
- 230000001133 acceleration Effects 0.000 claims abstract description 31
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 12
- 238000006073 displacement reaction Methods 0.000 description 10
- 238000004088 simulation Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 3
- 230000003746 surface roughness Effects 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/007—Subject matter not provided for in other groups of this subclass by applying a load, e.g. for resistance or wear testing
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域technical field
本发明属于建筑结构损伤识别技术领域,具体涉及基于车桥耦合的桥梁阻尼比识别方法。The invention belongs to the technical field of building structure damage identification, in particular to a bridge damping ratio identification method based on vehicle-bridge coupling.
背景技术Background technique
桥梁作为国民交通工程的重要组成部分,对我国的社会经济发展起到至关重要的作用。如果桥梁由于某种原因意外失效,可能会造成财产损失甚至重大人员伤亡。因此从结构安全的角度了解桥梁的损伤特性便显得尤为重要,现阶段的桥梁健康监测方法是将传感器直接布置在待测桥梁上以获取桥梁参数,称之为直接测量法,但该方法存在传感器在持续测量中易老化损坏,同时耗时耗力、阻碍交通、成本高等大量不利因素。As an important part of national traffic engineering, bridges play a vital role in the social and economic development of our country. If the bridge fails unexpectedly for some reason, it may cause property damage and even heavy casualties. Therefore, it is particularly important to understand the damage characteristics of bridges from the perspective of structural safety. The current bridge health monitoring method is to directly arrange sensors on the bridge to be tested to obtain bridge parameters, which is called direct measurement method, but this method has sensors It is prone to aging and damage during continuous measurement, and at the same time, it is time-consuming and labor-intensive, hindering traffic, and costing a lot of disadvantages.
基于直接测试的不足,杨永斌等于2004年提出的移动智能桥梁检测法,即测试车辆在桥梁上运行,基于测试车辆上传感器采集的信号处理,识别桥梁参数直至识别损伤的间接量测思想应运而出。之后各国学者前赴后继,将该方法应用于桥梁频率、模态、和阻尼等参数识别,并最终识别桥梁损伤。Yang等人基于车桥耦合的研究基础提取桥梁最基本的参数—频率。该团队引入经验模态分解(EDM),生成固有模态函数(IMFs),最后对其进行快速傅里叶变换(FFT)成功识别频率、然后通过引入奇异谱分析带通滤波等改进频率识别效果,同时基于粗糙度对识别效果的不利影响和移动智能车测量法的特征,yang等人提出了利用两辆测试车两次通过同一路面,再进行信号相减而降低路面粗糙度的影响,此外该研究团队提出利用车桥接触点信号提取桥梁频率,结果表明接触点响应能够更好地提取桥频。此后研究者开始尝试识别桥梁的另外一个参数—模态。2014年yang教授团队通过车体响应经过希尔伯特变换(HT)构造瞬时振幅从而成功提取桥梁模态,同年Obrien和Malekjafarian基于车桥耦合系统有限元模拟,结合奇异值分解和短时频域分解的概念从局部到整体分段、分步构造了桥梁模态,Li和Au等人改进了Yang等提出的模态识别方法,通过在车体上施加激励以增大桥梁响应,提高抗噪性,从而改进模态识别效果。桥梁阻尼--桥梁健康监测中的另一个重要参数,许多研究指出,阻尼可能对桥梁损伤更敏感。Obrien和coworkers提出了一种利用检测车-拖车车辆系统检测桥梁阻尼变化的方法。该方法需要大量重复的模拟,同时缺乏清晰的理论基础。Gonzalez等提出了双自由度检测车的车体响应识别桥梁阻尼的六步算法,采用迭代法对桥梁阻尼进行了较为准确的识别。张斌等通过安装于两轴移动试验车辆上的加速度传感器和位移传感器所测信号获取接触点信号,用于识别桥梁阻尼比,但缺乏实验验证。Based on the shortcomings of direct testing, Yang Yongbin et al. proposed a mobile intelligent bridge detection method in 2004, that is, the test vehicle runs on the bridge, based on the signal processing collected by the sensors on the test vehicle, the indirect measurement idea of identifying bridge parameters and identifying damage came out as the times require . After that, scholars from all over the world have applied this method to the identification of parameters such as bridge frequency, mode, and damping, and finally identified bridge damage. Yang et al. extracted the most basic parameter of the bridge—frequency based on the research basis of vehicle-bridge coupling. The team introduced empirical mode decomposition (EDM) to generate intrinsic mode functions (IMFs), and finally performed fast Fourier transform (FFT) on it to successfully identify frequencies, and then improved the frequency identification effect by introducing singular spectrum analysis bandpass filtering, etc. At the same time, based on the adverse effect of roughness on the recognition effect and the characteristics of the mobile smart vehicle measurement method, Yang et al. proposed to use two test vehicles to pass the same road twice, and then perform signal subtraction to reduce the impact of road roughness. In addition, The research team proposed to use the vehicle-bridge contact point signal to extract the bridge frequency, and the results show that the contact point response can better extract the bridge frequency. Since then, researchers have begun to try to identify another parameter of the bridge - the mode. In 2014, Professor Yang's team successfully extracted the bridge mode by constructing the instantaneous amplitude through the Hilbert transform (HT) of the vehicle body response. The concept of decomposition constructs the bridge modes from local to whole segment and step by step. Li and Au et al. improved the mode identification method proposed by Yang et al. By applying excitation on the vehicle body to increase the bridge response and improve anti-noise so as to improve the mode recognition effect. Bridge damping - another important parameter in bridge health monitoring, many studies have pointed out that damping may be more sensitive to bridge damage. Obrien and coworkers proposed a method to detect changes in bridge damping using a vehicle-trailer vehicle system. This method requires a large number of repeated simulations and lacks a clear theoretical basis. Gonzalez et al. proposed a six-step algorithm for identifying bridge damping from the body response of a two-degree-of-freedom detection vehicle, and used an iterative method to identify bridge damping more accurately. Zhang Bin et al. obtained the contact point signal through the signals measured by the acceleration sensor and the displacement sensor installed on the two-axis mobile test vehicle, which were used to identify the bridge damping ratio, but lacked experimental verification.
以上内容讲述了桥梁各参数的识别过程,然而阻尼比的识别方法尚不成熟。The above content describes the identification process of bridge parameters, but the identification method of damping ratio is not yet mature.
发明内容SUMMARY OF THE INVENTION
针对现有技术中存在的上述不足之处,本发明提供了基于车桥耦合的桥梁阻尼比识别方法,识别结果较好,便于操作。In view of the above deficiencies in the prior art, the present invention provides a bridge damping ratio identification method based on vehicle-bridge coupling, which has better identification results and is easy to operate.
为了解决上述技术问题,本发明采用了如下技术方案:In order to solve the above-mentioned technical problems, the present invention adopts the following technical solutions:
基于车桥耦合的桥梁阻尼比识别方法,包括以下步骤:The identification method of bridge damping ratio based on vehicle-bridge coupling includes the following steps:
步骤一、利用移动的测试车按一定采样频率采集桥梁各点的加速度响应信号;Step 1. Use the mobile test vehicle to collect the acceleration response signal of each point of the bridge at a certain sampling frequency;
步骤二、通过带通滤波器对步骤一的加速度响应信号进行滤波处理,获取桥梁一阶频率的加速度响应;Step 2: Filter the acceleration response signal of Step 1 through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge;
步骤三、假定一阻尼比,桥梁原频率ω1已知,则可获得指数函数的数值,然后将步骤二滤波后的加速度响应除以该指数函数,即可获得无阻尼桥梁滤波后的加速度响应;Step 3. Assuming a damping ratio and the original bridge frequency ω 1 is known, the exponential function can be obtained , and then divide the filtered acceleration response in
步骤四、通过短时频域分解法获取步骤三处理后信号的模态,通过辨别模态的最大值是否处在模态中点从而判断该假定阻尼比是否为真实阻尼比,如若模态的最大值处在模态中点,则可将该假定阻尼比作为桥梁真实阻尼比,否则重新假定一阻尼比,直至模态的最大值处在模态中点,最终识别桥梁阻尼比。
进一步,步骤二中带通滤波器的上限值为 Further, the upper limit of the band-pass filter in
进一步,步骤二中带通滤波器的下限值为 Further, the lower limit of the bandpass filter in
进一步,步骤一中测试车的采样频率为100Hz。Further, the sampling frequency of the test vehicle in step 1 is 100 Hz.
进一步,步骤一中的测试车为单自由度车。Further, the test vehicle in step 1 is a single-degree-of-freedom vehicle.
进一步,步骤一中测试车的质量为1000kg。Further, the mass of the test vehicle in step 1 is 1000kg.
进一步,步骤一中测试车的车速恒为1m/s。Further, the speed of the test vehicle in step 1 is always 1 m/s.
本发明与现有技术相比,桥梁阻尼比的识别结果较好,便于操作。Compared with the prior art, the invention has better identification results of the bridge damping ratio and is easy to operate.
附图说明Description of drawings
图1为本发明基于车桥耦合的桥梁阻尼比识别方法的流程图;Fig. 1 is the flow chart of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention;
图2为本发明基于车桥耦合的桥梁阻尼比识别方法的双车桥系统简化模型图;2 is a simplified model diagram of a double-axle system based on a vehicle-axle coupling-based bridge damping ratio identification method of the present invention;
图3为本发明基于车桥耦合的桥梁阻尼比识别方法的桥梁模型中单元、节点编号示意图;3 is a schematic diagram of the number of elements and nodes in the bridge model of the bridge damping ratio identification method based on vehicle-bridge coupling according to the present invention;
图4为本发明基于车桥耦合的桥梁阻尼比识别方法的测试车传感器采集时域信号图;FIG. 4 is a time-domain signal diagram of a test vehicle sensor collection based on a vehicle-bridge coupling-based bridge damping ratio identification method of the present invention;
图5为本发明基于车桥耦合的桥梁阻尼比识别方法的对采集信号滤波后的时域信号图;Fig. 5 is the time domain signal diagram after filtering the collected signal of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention;
图6为本发明基于车桥耦合的桥梁阻尼比识别方法的除以衰减信号后的滤波时域信号图;6 is a filtered time-domain signal diagram of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention divided by the attenuation signal;
图7为本发明基于车桥耦合的桥梁阻尼比识别方法的拟合模态图;Fig. 7 is the fitting modal diagram of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention;
图8为本发明基于车桥耦合的桥梁阻尼比识别方法的做对称处理后的拟合模态图;8 is a fitting modal diagram after symmetrical processing of the bridge damping ratio identification method based on vehicle-bridge coupling according to the present invention;
图9为本发明基于车桥耦合的桥梁阻尼比识别方法在A级粗糙度下不同假定阻尼比的识别结果图;Fig. 9 is the identification result diagram of different assumed damping ratios under the A-level roughness of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention;
图10为本发明基于车桥耦合的桥梁阻尼比识别方法在B级粗糙度下不同假定阻尼比的识别结果图;10 is a diagram showing the identification results of different assumed damping ratios under the B-level roughness of the bridge damping ratio identification method based on vehicle-bridge coupling of the present invention;
图11为本发明基于车桥耦合的桥梁阻尼比识别方法在50db下阻尼比识别结果图;11 is a diagram showing the identification result of the damping ratio of the bridge damping ratio identification method based on vehicle-bridge coupling under 50db;
图12为本发明基于车桥耦合的桥梁阻尼比识别方法在40db下阻尼比识别结果图;12 is a diagram showing the identification result of the damping ratio of the bridge damping ratio identification method based on vehicle-bridge coupling under 40db;
图13为本发明基于车桥耦合的桥梁阻尼比识别方法在30db下阻尼比识别结果图。FIG. 13 is a diagram showing the identification result of the damping ratio of the bridge damping ratio identification method based on vehicle-bridge coupling under 30db.
具体实施方式Detailed ways
为了使本领域的技术人员可以更好地理解本发明,下面结合附图和实施例对本发明技术方案进一步说明。In order to enable those skilled in the art to better understand the present invention, the technical solutions of the present invention are further described below with reference to the accompanying drawings and embodiments.
如图1所示,本发明的基于车桥耦合的桥梁阻尼比识别方法,包括以下步骤:As shown in FIG. 1 , the method for identifying bridge damping ratio based on vehicle-bridge coupling of the present invention includes the following steps:
步骤一、利用移动的测试车按一定采样频率采集桥梁各点的加速度响应信号;Step 1. Use the mobile test vehicle to collect the acceleration response signal of each point of the bridge at a certain sampling frequency;
步骤二、通过带通滤波器对步骤一的加速度响应信号进行滤波处理,获取桥梁一阶频率的加速度响应;Step 2: Filter the acceleration response signal of Step 1 through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge;
步骤三、假定一阻尼比,桥梁原频率ω1已知,则可获得指数函数的数值,然后将步骤二滤波后的加速度响应除以该指数函数,即可获得无阻尼桥梁滤波后的加速度响应;Step 3. Assuming a damping ratio and the original bridge frequency ω 1 is known, the exponential function can be obtained , and then divide the filtered acceleration response in
步骤四、通过短时频域分解法获取步骤三处理后信号的模态,通过辨别模态的最大值是否处在模态中点从而判断该假定阻尼比是否为真实阻尼比,如若模态的最大值处在模态中点,则可将该假定阻尼比作为桥梁真实阻尼比,否则重新假定一阻尼比,直至模态的最大值处在模态中点,最终识别桥梁阻尼比。
作为优选方案,步骤二中带通滤波器的上限值为 As a preferred solution, the upper limit of the band-pass filter in
作为优选方案,步骤二中带通滤波器的下限值为 As a preferred solution, the lower limit of the band-pass filter in
作为优选方案,步骤一中测试车的采样频率为100Hz。As a preferred solution, the sampling frequency of the test vehicle in step 1 is 100 Hz.
作为优选方案,步骤一中的测试车为单自由度车。As a preferred solution, the test vehicle in step 1 is a single-degree-of-freedom vehicle.
作为优选方案,步骤一中测试车的质量为1000kg。As a preferred solution, the mass of the test vehicle in step 1 is 1000kg.
作为优选方案,步骤一中测试车的车速恒为1m/s。As a preferred solution, the speed of the test vehicle in step 1 is always 1m/s.
在图2中,距离为△L的测试车和牵引车以恒定速度v在桥面上行驶。测试车和牵引车简化为支撑在其上的运动弹簧质量mv,2和mv,1,阻尼系数为cv2、cv1的缓冲器和刚度为kv2、kv1的弹簧。桥梁为长度L,每单位长度的质量m*和弯曲刚度EI(n阶阻尼比)的简单支撑,此处的弯曲刚度EI包括非结构性构件的作用,例如桥梁栏杆和桥梁桥面板实际桥梁中的人行道和单位长度m*的质量可以从实际桥梁设计数据中获得,也可以通过横截面积步估算。假定桥梁在测试车到达之前处于静止状态。In Figure 2, the test vehicle and the tractor with a distance of ΔL travel on the bridge at a constant speed v. The test vehicle and the tractor are simplified to the moving spring masses m v,2 and m v,1 supported on them, the buffers with damping coefficients cv2, cv1 and springs with stiffness kv2, kv1. The bridge is a simple support of length L, mass per unit length m* and bending stiffness EI (nth order damping ratio), where bending stiffness EI includes the effect of non-structural members such as bridge railings and bridge decks in real bridges The sidewalk and mass per unit length m* can be obtained from actual bridge design data or estimated from the cross-sectional area step. The bridge is assumed to be stationary until the test vehicle arrives.
该车桥耦合系统的运动方程可写为:The equation of motion of the vehicle-axle coupling system can be written as:
其中u(x,t)表示桥梁结构距左支撑点x处的垂直位移,uv1(t),uv2(t)分别是测试车和牵引车从其静态平衡位置和时间ta=t+ΔL/v=t+ts开始测量的静态位移,应注意,时间t是指车进入桥后的时间。带点的字符表示相对于时间t和坐标x的导数。测试车和牵引车与平台接触点的相互作用力fc1(t),fc2(t)可分别写为:where u(x,t) represents the vertical displacement of the bridge structure from the left support point x, u v1 (t), u v2 (t) are the test vehicle and the tractor from their static equilibrium position and time t a = t + ΔL/v=t+t s The static displacement measured at the beginning, it should be noted that the time t refers to the time after the vehicle enters the bridge. The dotted character represents the derivative with respect to time t and coordinate x. The interaction forces f c1 (t) and f c2 (t) of the contact point between the test vehicle and the tractor and the platform can be written as:
g表示重力加速度g is the acceleration of gravity
桥梁位移可以用简支梁的广义坐标qn(t)和模态sin(nπx/L)表示为:The bridge displacement can be expressed as:
假设车辆质量mv1,mv2比桥面板的质量小得多,mv1<<m*L和mv2<<m*L,对于实际的桥梁,此假设很容易实现。通过将等式(6)代入等式(3),乘以sin(nπx/L),并从0到L积分,然后根据正弦函数的正交条件,可以将结构的第n个模态平衡方程写为:Assuming that the vehicle mass mv1, mv2 is much smaller than the mass of the bridge deck, m v1 << m * L and m v2 << m * L, this assumption is easy to realize for a practical bridge. By substituting equation (6) into equation (3), multiplying by sin(nπx/L), and integrating from 0 to L, the nth modal equilibrium equation of the structure can then be calculated from the quadrature condition of the sine function written as:
ωn is为桥面板的第n阶角频率ω n is the nth angular frequency of the bridge deck
对于零初始条件,可以从等式(7)获得桥的广义坐标qn(t)为:For zero initial conditions, the generalized coordinates qn( t ) of the bridge can be obtained from equation (7) as:
将公式(9)代入公式(6),得出桥位移为:Substituting formula (9) into formula (6), the bridge displacement is obtained as:
将方程式(10)代入方程式(1)也可得出测试车的位移uv1(t)。Substituting equation (10) into equation (1) also yields the displacement u v1 (t) of the test vehicle.
实际上,可以将与桥梁第n个模态频率相关的响应分量与测试车的响应分开。在这项研究中采用了上下限值为和的带通滤波器。产生的信号是来自桥梁结构第n个振动模式的瞬态响应,该瞬态响应与方程式(10)的最后一项直接相关。In fact, the response components related to the nth modal frequency of the bridge can be separated from the response of the test vehicle. In this study, the upper and lower limits were used as and bandpass filter. The resulting signal is the transient response from the nth vibration mode of the bridge structure, which is directly related to the last term of equation (10).
将方程式(10)的最后一项代入方程式(1),可获得与桥面板的第n模态形状有关(11)的车辆位移为:Substituting the last term of equation (10) into equation (1), the vehicle displacement related to the shape of the nth mode of the bridge deck (11) can be obtained as:
上述公式中测试车阻尼比为ξv1=cv1/(2mv1ωv1),此外In the above formula, the damping ratio of the test vehicle is ξ v1 =c v1 /(2m v1 ω v1 ), in addition
在改进的直接刚度法中,仅需将桥面板的单一振动模式用于损伤识别。由于不能保证在现场获取精确的高频率振动模态,因此除非另有说明,否则以下讨论中的频率和模态是指桥面板的第一振动模态。与桥面的第一振动模态相关的测试车的响应分量R1也可以写成:In the modified direct stiffness method, only a single vibration mode of the bridge deck needs to be used for damage identification. Since accurate high frequency vibration modes cannot be guaranteed in the field, unless otherwise stated, frequencies and modes in the following discussion refer to the first vibration mode of the bridge deck. The response component R1 of the test vehicle related to the first vibration mode of the bridge deck can also be written as:
可以通过将等式(9)和(10)中的项进行比较来确定系数A1至A2,The coefficients A1 to A2 can be determined by comparing the terms in equations (9) and (10),
桥面板的第一振动模态的相应加速度响应分量也可以通过以下方式获得:The corresponding acceleration response components of the first vibration mode of the bridge deck can also be obtained by:
观察式(16):当ξ1=0时,Observation formula (16): when ξ 1 =0,
该式即为无阻尼桥面板第一振动模态的的加速度响应分量。则它们之间的关系为This formula is the acceleration response component of the first vibration mode of the undamped bridge deck. Then the relationship between them is
即有阻尼桥梁所获得的加速度响应与无阻尼桥梁的加速度响应之比即为一个指数函数。利用这一关系即可提出识别阻尼比方法。That is, the ratio of the acceleration response obtained by the damped bridge to the acceleration response of the undamped bridge is an exponential function. Using this relationship, a method for identifying the damping ratio can be proposed.
为了验证上述理论的正确性,在单车无粗糙度下验证该方法的可行性。下面利用有限元进行数值模拟,来验证该方法的可靠性。此次数值模拟拟定桥梁桥长30米,单位长度质量m*=19116kg,截面面积A=7.965m2,截面惯性矩Ix=2.959m4,桥梁弹性模量E=2.9×1010N/m2,测试车质量1000kg,车速恒为1m/s,采样频率为100Hz,假定桥梁真实阻尼比为0.01。In order to verify the correctness of the above theory, the feasibility of the method is verified under the condition of no roughness on a bicycle. Numerical simulation is carried out using finite element to verify the reliability of the method. In this numerical simulation, the bridge length is 30 meters, the mass per unit length m * =19116kg, the cross-sectional area A= 7.965m2 , the cross-sectional inertia moment Ix=2.959m4, the bridge elastic modulus E=2.9×1010N/m2, and the mass of the test vehicle 1000kg, the vehicle speed is constant 1m/s, the sampling frequency is 100Hz, and the real damping ratio of the bridge is assumed to be 0.01.
在利用短时频域分解法提取模态时,将桥梁分为10个单元,分别为E1~E10,如图3所示,其余数字为单元节点编号(j=1,2,…,11),需要注意的是,本发明得到模态的节点为2~10,节点1和11正好为入桥和出桥的位置,由于支座的约束,该位置振动很弱,几乎可以忽略,故本发明不考虑边单元的节点模态值。When using the short-time frequency domain decomposition method to extract the modes, the bridge is divided into 10 units, which are E1~E10, as shown in Figure 3, and the remaining numbers are the unit node numbers (j=1,2,...,11) , it should be noted that the nodes of the mode obtained by the present invention are 2 to 10, and
图4是行驶在桥梁上的测试车以100Hz的采样频率获取的加速度时域信号图(此时信号含有车频和桥频信息),然后通过带通滤波器滤波,获取桥梁一阶频率的加速度信号。如图5,此时假定一阻尼比,假定真实阻尼比为0.01,则衰减系数为常数,利用上述滤波后信号除以该衰减系数得到图6信号图,同时采用短时频域分解法获取图7模态信号,观察图6信号可知,与无桥梁阻尼相比,前半段基本吻合,越往后,信号结果越差,这也导致图7信号在尾端出现失真现象,这是因为桥梁阻尼至桥梁后半段衰减作用更强,桥频信号可能衰减殆尽,所以相除后还原效果不好,此时采用部分信号对称的做法获得图8的模态,发明人通过大量数值模拟得出图8结果,当假定阻尼比为0.01时(此时正好为真实阻尼比),模态最大值居中,当假定的阻尼比为0.0092时(此时小于真实阻尼比),该模态最大值偏左,即出现左偏,当假定的阻尼比为0.0107时(此时大于真实阻尼比),该模态最大值偏右,即出现右偏,而当假定阻尼比处于0.0092-0.0107之间时,此时左偏右偏现象不明显,均可作为真实阻尼比,此时该方法的最大误差为处于误差范围内。可以得知在单车无外界信号干扰下的识别结果很好。但实际情况下外界各类影响还是存在的。Figure 4 is the acceleration time-domain signal diagram obtained by the test vehicle driving on the bridge at the sampling frequency of 100Hz (the signal contains the vehicle frequency and bridge frequency information at this time), and then filtered by the band-pass filter to obtain the acceleration of the first-order frequency of the bridge Signal. As shown in Figure 5, a damping ratio is assumed at this time, and the actual damping ratio is assumed to be 0.01, then the attenuation coefficient is a constant, the signal diagram in Figure 6 is obtained by dividing the above filtered signal by the attenuation coefficient, and the short-time frequency domain decomposition method is used to obtain the modal signal in Figure 7. Observing the signal in Figure 6, it can be seen that compared with no bridge damping, the first half is basically match, the further back, the worse the signal result, which also leads to the distortion of the signal in Figure 7 at the tail end, this is because the bridge damping to the second half of the bridge has a stronger attenuation effect, and the bridge frequency signal may be attenuated completely, so divide The post-reduction effect is not good. At this time, the mode of Fig. 8 is obtained by the method of partial signal symmetry. The inventor obtained the result of Fig. 8 through a large number of numerical simulations. When the damping ratio is assumed to be 0.01 (this is exactly the real damping ratio), The modal maximum value is centered. When the assumed damping ratio is 0.0092 (at this time, it is smaller than the real damping ratio), the modal maximum value is left, that is, left deviation occurs. When the assumed damping ratio is 0.0107 (at this time, it is greater than the real damping ratio), the modal maximum value is to the right, that is, right deviation occurs, and when the assumed damping ratio is between 0.0092-0.0107, the phenomenon of left deviation to right deviation is not obvious at this time, which can be used as the real damping ratio. At this time, the method The maximum error is within the margin of error. It can be seen that the recognition result of the bicycle without external signal interference is very good. But in reality, various external influences still exist.
以上是发明人基于真实阻尼比为0.01的情况下所得出的结论,为得到该方法所能识别阻尼比的适用范围,发明人经过大量数值模拟得出,在真实阻尼比小于或等于0.02时,该方法都能有效识别阻尼比,而真实阻尼比大于0.02无法识别的原因如上述所述,当真实阻尼比越大,桥梁信号衰减的越快,即在信号后端的模态识别效果越不佳。综上该方法适用于小于或等于0.02的小阻尼识别。The above is the conclusion drawn by the inventor based on the fact that the real damping ratio is 0.01. In order to obtain the applicable range of the damping ratio that can be identified by this method, the inventor has obtained through a large number of numerical simulations that when the real damping ratio is less than or equal to 0.02, This method can effectively identify the damping ratio, but the reason why the real damping ratio is greater than 0.02 cannot be identified is as described above. When the real damping ratio is larger, the bridge signal attenuates faster, that is, the mode identification effect at the back end of the signal is worse. . In conclusion, this method is suitable for the identification of small damping less than or equal to 0.02.
关于粗糙度的影响About the effect of roughness
实际桥梁路面是存在不平整度的,故本发明粗糙度采用ISO 8608(1995)标准[36]所建议的功能密度函数来模拟。其功能密度函数Gd(n)如下所示:The actual bridge pavement has unevenness, so the roughness of the present invention is simulated by the functional density function suggested by the ISO 8608 (1995) standard [36]. Its functional density function G d (n) is as follows:
式中,n为每单位长度的空间频率、w为常数2、n0为0.1cycle/m、Gd(n0)为位移功能密度函数值,由路面粗糙度等级确定。ISO 8608(1995)标准提供的函数值的几何平均值的平方根,即各级路面粗糙度位移功能密度函数Gd(n0)取值分别为:In the formula, n is the spatial frequency per unit length, w is a constant 2, n 0 is 0.1 cycle/m, and G d (n 0 ) is the displacement function density function value, which is determined by the road surface roughness grade. The square root of the geometric mean value of the function value provided by the ISO 8608 (1995) standard, that is, the functional density function G d (n 0 ) of the road surface roughness displacement at all levels is as follows:
A级:Gd(n0)=4×10-6m3;B级:Gd(n0)=8×10-6m3;C级:Gd(n0)=16×10-6m3 Class A: G d (n 0 )=4×10 −6 m 3 ; Class B: G d (n 0 )=8×10 −6 m 3 ; Class C: G d (n 0 )=16×10 − 6 m 3
各级粗糙度下的路面粗糙度位移振幅值d可表示为:The displacement amplitude value d of the pavement roughness at all levels of roughness can be expressed as:
式中,Δn为空间频率的采样间隔where Δn is the sampling interval of the spatial frequency
接着以不同空间频率的余弦函数叠加来模拟路面粗糙度r(x),可表示为:Then, the cosine functions of different spatial frequencies are superimposed to simulate the road surface roughness r(x), which can be expressed as:
式中,ns,i为粗糙度的空间频率,di、θi分别为粗糙度的幅值与随机相位角。where ns ,i is the spatial frequency of roughness, d i and θ i are the amplitude and random phase angle of roughness, respectively.
为了研究在有粗糙度的影响下桥梁阻尼比的识别情况,此次数值模拟采用杨永斌等人提出的利用两辆测试车两次通过同一路面,再进行信号相减而降低路面粗糙度的影响的方法。其中大测试车体质量2000kg,刚度为20000N/m,小测试车质量1000kg,刚度为10000N/m.并分别研究A级,B级,C级粗糙度下阻尼识别精度的确定。并假定桥梁真实阻尼比为0.01.然后再假定一系列阻尼比并运用上述方法来识别真实阻尼比。识别结果如图9-10所示。In order to study the identification of bridge damping ratio under the influence of roughness, this numerical simulation adopts the method proposed by Yang Yongbin et al. to use two test vehicles to pass through the same road twice, and then perform signal subtraction to reduce the influence of road roughness. method. Among them, the mass of the large test car is 2000kg and the stiffness is 20000N/m, and the mass of the small test car is 1000kg and the stiffness is 10000N/m. The determination of the damping identification accuracy under the roughness of grade A, grade B and grade C is studied respectively. And assume that the real damping ratio of the bridge is 0.01. Then assume a series of damping ratios and use the above method to identify the real damping ratio. The recognition result is shown in Figure 9-10.
由图9、图10可以看出在A、B级粗糙度下,当假定阻尼比正好为真实阻尼比时,或者在真实阻尼比附近,模态最大值都处在模态中点,而当假定阻尼比偏离真实阻尼比较大时,都会出现不同程度的左偏或者右偏,即最大值不在模态中点。正如图9所示,在B级粗糙度下,当假定阻尼比处于0.0095-0.0103时,此时模态最大值都处在模态中点,而当假定阻尼比小于0.0095时,即会出现左偏现象,当假定阻尼比大于0.0103时,即会出现右偏现象,这一现象同时可以帮助我们如何假定阻尼比,从而更有效率的去识别真实阻尼比,而此时的最大误差为符合误差范围。It can be seen from Figure 9 and Figure 10 that under the roughness of grades A and B, when the damping ratio is assumed to be exactly the real damping ratio, or near the real damping ratio, the modal maximum value is at the modal midpoint, and when It is assumed that when the damping ratio deviates from the real damping relatively large, there will be different degrees of left deviation or right deviation, that is, the maximum value is not at the midpoint of the mode. As shown in Figure 9, under the B-level roughness, when the damping ratio is assumed to be 0.0095-0.0103, the modal maximum value is at the modal midpoint, and when the damping ratio is assumed to be less than 0.0095, there will be left When the damping ratio is assumed to be greater than 0.0103, there will be a right-bias phenomenon. This phenomenon can also help us how to assume the damping ratio, so as to identify the real damping ratio more efficiently, and the maximum error at this time is within the margin of error.
在C级粗糙度下阻尼比识别结果不理想,这是因为此时粗糙度过大,利用信号相减降低粗糙度的效果不理想,所以粗糙度消除不理想,此外粗糙度过大同样导致后端信号获取不理想(这也是在方法介绍中作部分对称处理的原因),从而得不到理想的模态信号,阻尼比因此识别不理想。The identification result of the damping ratio is not ideal under the C-level roughness. This is because the roughness is too large at this time, and the effect of reducing the roughness by signal subtraction is not ideal, so the roughness elimination is not ideal. In addition, the roughness is too large. The acquisition of the terminal signal is not ideal (this is also the reason for the partial symmetry processing in the method introduction), so that the ideal modal signal cannot be obtained, and the damping ratio is therefore not ideal for identification.
关于噪音的影响About the effect of noise
在实际应用中,加速度传感器采集的信号不可避免会受到噪音的干扰,为了探究在噪音干扰下该阻尼比识别方法的有效性,本发明通过在数值模拟中获得的加速度信号添加高斯白噪声的方法来研究该方法的抗噪性,并用信噪比作为指标,信噪比的定义如下:In practical applications, the signal collected by the acceleration sensor will inevitably be interfered by noise. In order to explore the effectiveness of the damping ratio identification method under noise interference, the present invention adds Gaussian white noise to the acceleration signal obtained in the numerical simulation. To study the noise resistance of the method, and use the signal-to-noise ratio as an indicator, the signal-to-noise ratio is defined as follows:
式中:N为数据点个数,yi为第i时刻含有噪声的测试车加速度响应,σi为第i时刻的噪声值,SNR为信噪比,单位为dB,其值越大,表示噪音影响越小,信号被干扰程度低,其值越小,表示噪音影响越大,信号被干扰程度越大。In the formula: N is the number of data points, yi is the acceleration response of the test vehicle with noise at the i-th moment, σ i is the noise value at the i-th moment, and SNR is the signal-to-noise ratio, in dB. The smaller the impact of noise, the lower the degree of signal interference, the smaller the value, the greater the impact of noise, the greater the degree of signal interference.
本发明假定真实阻尼比为0.01,然后假设0.01的阻尼比加之噪音的干扰,从而分析识别结果的误差。此次模拟设置为在A级粗糙度下研究噪音的影响,并对每个水平的噪声进行数值模拟,取其值来获取模态振型识别。The present invention assumes that the real damping ratio is 0.01, and then assumes the damping ratio of 0.01 plus the interference of noise, so as to analyze the error of the identification result. The simulation was set up to study the effect of noise at a roughness class A, and numerically simulate the noise at each level, taking its value to obtain modal mode shape identification.
此次模拟绘制了3种不同噪音水平下阻尼比识别模态结果,由图11、图12和图13可知:当信噪比为50db时,当假定的阻尼比处于0.091-0.0108时,此时模态最大值都处在模态中点,而当假定阻尼比小于或等于0.0091时,此时模态最大值都不处在模态中点,即会出现左偏现象,当假定阻尼比大于0.0108时,即会出现右偏现象,而此时的最大误差为符合误差范围。而随着信噪比越小,噪音干扰越大,该方法同样能够很好地识别阻尼比,而且识别精度更加准确。但是信噪比小于或等于20db时,此时噪音的干扰已经较大,此时识别的结果较差,这也就说明该方法识别阻尼比有一定的抗噪性。This simulation plots the modal results of damping ratio identification under three different noise levels. It can be seen from Figure 11, Figure 12 and Figure 13 that: when the signal-to-noise ratio is 50db, when the assumed damping ratio is 0.091-0.0108, then The modal maximum value is all at the modal midpoint, and when the assumed damping ratio is less than or equal to 0.0091, the modal maximum value is not at the modal midpoint, that is, a left-bias phenomenon occurs. When the assumed damping ratio is greater than When 0.0108, there will be a right-bias phenomenon, and the maximum error at this time is within the margin of error. With the smaller the signal-to-noise ratio, the greater the noise interference, the method can also identify the damping ratio well, and the identification accuracy is more accurate. However, when the signal-to-noise ratio is less than or equal to 20db, the interference of noise is already large at this time, and the recognition result is poor at this time, which means that the method has a certain anti-noise performance for the recognition damping ratio.
综合上述数值模拟结果,在有粗糙度和噪音的影响下,在A、B级粗糙度下,在信噪比大于等于30db时,阻尼比的识别结果较好,所以该方法能够识别阻尼比。Based on the above numerical simulation results, under the influence of roughness and noise, under the roughness of grades A and B, when the signal-to-noise ratio is greater than or equal to 30db, the identification result of the damping ratio is better, so this method can identify the damping ratio.
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the preferred embodiments, those of ordinary skill in the art should understand that the technical solutions of the present invention can be Modifications or equivalent substitutions without departing from the spirit and scope of the technical solutions of the present invention should be included in the scope of the claims of the present invention.
Claims (7)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010678822.0A CN111781001B (en) | 2020-07-15 | 2020-07-15 | Bridge damping ratio identification method based on axle coupling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010678822.0A CN111781001B (en) | 2020-07-15 | 2020-07-15 | Bridge damping ratio identification method based on axle coupling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111781001A true CN111781001A (en) | 2020-10-16 |
CN111781001B CN111781001B (en) | 2022-07-19 |
Family
ID=72767209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010678822.0A Active CN111781001B (en) | 2020-07-15 | 2020-07-15 | Bridge damping ratio identification method based on axle coupling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111781001B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326787A (en) * | 2020-10-20 | 2021-02-05 | 中国电建集团重庆工程有限公司 | Beam bridge identification method based on multipoint rapid static acquisition of exclusive test car |
CN115326333A (en) * | 2022-08-17 | 2022-11-11 | 绍兴文理学院 | A method and system for identifying modal damping of building structures |
CN115436473A (en) * | 2022-09-05 | 2022-12-06 | 重庆大学 | Method for identifying structural damage based on timing synchronization theory |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101587007A (en) * | 2009-06-19 | 2009-11-25 | 同济大学 | Output-only wavelet analytical method for recognizing flexible bridge structure kinetic parameter |
CN102621225A (en) * | 2012-03-30 | 2012-08-01 | 东南大学 | Method for testing damping characteristic parameter of road surface and bridge deck pavement material |
US20150111062A1 (en) * | 2013-10-22 | 2015-04-23 | Mo-How Herman Shen | High strain damping method including a face-centered cubic ferromagnetic damping coating, and components having same |
CN106202817A (en) * | 2016-07-26 | 2016-12-07 | 南京航空航天大学 | One compares obtaining value method based on field measurement large cooling column comprehensive equivalent damping |
CN106802222A (en) * | 2017-01-13 | 2017-06-06 | 重庆大学 | A kind of bridge damnification diagnostic method based on Vehicle-Bridge Coupling System |
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 |
CN107727338A (en) * | 2017-06-01 | 2018-02-23 | 重庆大学 | A kind of bridge damnification diagnostic method based on Vehicle-Bridge Coupling System |
KR101864665B1 (en) * | 2017-03-30 | 2018-06-05 | 동아대학교 산학협력단 | Method for estimating damping ratio of a beam through electromechanical signatures of piezoelectric transducer mounted on the beam |
CN109813511A (en) * | 2019-01-31 | 2019-05-28 | 东南大学 | Rapid test and parameter identification method of bridge based on moving vehicle |
CN109855823A (en) * | 2019-01-25 | 2019-06-07 | 重庆大学 | A method of Damage Identification of Bridge Structure is carried out using test carriage |
CN110568069A (en) * | 2019-09-11 | 2019-12-13 | 重庆大学 | Damage identification method for girder bridges based on average value of acceleration ratio and GPSA algorithm |
-
2020
- 2020-07-15 CN CN202010678822.0A patent/CN111781001B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101587007A (en) * | 2009-06-19 | 2009-11-25 | 同济大学 | Output-only wavelet analytical method for recognizing flexible bridge structure kinetic parameter |
CN102621225A (en) * | 2012-03-30 | 2012-08-01 | 东南大学 | Method for testing damping characteristic parameter of road surface and bridge deck pavement material |
US20150111062A1 (en) * | 2013-10-22 | 2015-04-23 | Mo-How Herman Shen | High strain damping method including a face-centered cubic ferromagnetic damping coating, and components having same |
CN106202817A (en) * | 2016-07-26 | 2016-12-07 | 南京航空航天大学 | One compares obtaining value method based on field measurement large cooling column comprehensive equivalent damping |
CN106960068A (en) * | 2016-09-30 | 2017-07-18 | 中国人民解放军海军工程大学 | A kind of damping ratios quick calculation method based on pulse excitation response spectrum |
CN106802222A (en) * | 2017-01-13 | 2017-06-06 | 重庆大学 | A kind of bridge damnification diagnostic method based on Vehicle-Bridge Coupling System |
CN106844935A (en) * | 2017-01-18 | 2017-06-13 | 大连理工大学 | A kind of big damping engineering structure Modal Parameters Identification |
KR101864665B1 (en) * | 2017-03-30 | 2018-06-05 | 동아대학교 산학협력단 | Method for estimating damping ratio of a beam through electromechanical signatures of piezoelectric transducer mounted on the beam |
CN107727338A (en) * | 2017-06-01 | 2018-02-23 | 重庆大学 | A kind of bridge damnification diagnostic method based on Vehicle-Bridge Coupling System |
CN109855823A (en) * | 2019-01-25 | 2019-06-07 | 重庆大学 | A method of Damage Identification of Bridge Structure is carried out using test carriage |
CN109813511A (en) * | 2019-01-31 | 2019-05-28 | 东南大学 | Rapid test and parameter identification method of bridge based on moving vehicle |
CN110568069A (en) * | 2019-09-11 | 2019-12-13 | 重庆大学 | Damage identification method for girder bridges based on average value of acceleration ratio and GPSA algorithm |
Non-Patent Citations (4)
Title |
---|
YANG YANG,ET AL: "Feasibility Study of Tractor-Test Vehicle Technique for Practical Structural Condition Assessment of Beam-Like Bridge Deck", 《REMOTE SENSING》 * |
利远翔,雷晓燕,张斌: "高速列车一无砟轨道一桥梁耦合系统振动特性分析", 《华东交通大学学报》 * |
张大志: "《力学基础实验》", 31 January 2018 * |
黄方林等: "识别结构模态阻尼比的一种新方法", 《土木工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112326787A (en) * | 2020-10-20 | 2021-02-05 | 中国电建集团重庆工程有限公司 | Beam bridge identification method based on multipoint rapid static acquisition of exclusive test car |
CN112326787B (en) * | 2020-10-20 | 2024-05-14 | 中国电建集团重庆工程有限公司 | Beam bridge identification method based on dedicated test vehicle multipoint rapid static acquisition |
CN115326333A (en) * | 2022-08-17 | 2022-11-11 | 绍兴文理学院 | A method and system for identifying modal damping of building structures |
CN115326333B (en) * | 2022-08-17 | 2024-11-26 | 绍兴文理学院 | A method and system for identifying modal damping of building structures |
CN115436473A (en) * | 2022-09-05 | 2022-12-06 | 重庆大学 | Method for identifying structural damage based on timing synchronization theory |
CN115436473B (en) * | 2022-09-05 | 2024-07-16 | 重庆大学 | Method for identifying structural damage based on time sequence synchronization theory |
Also Published As
Publication number | Publication date |
---|---|
CN111781001B (en) | 2022-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111781001B (en) | Bridge damping ratio identification method based on axle coupling | |
CN109855823B (en) | Method for identifying damage of bridge structure by using test vehicle | |
CN105139086B (en) | Amplitude Estimation Method of Track Height Irregularity Based on Optimal Confidence Rule Reasoning | |
CN106524967B (en) | A kind of automobile core wheel actual travel displacement measurement and extracting method | |
CN109357822A (en) | A rapid test and evaluation method for bridges based on time-varying dynamic characteristics of vehicle-bridge coupling system | |
CN116481747B (en) | Bridge modal shape driving identification method considering bridge damping influence | |
CN103674578B (en) | Bullet train runs the detection method of dynamic performance state | |
CN106451498B (en) | A kind of low-frequency oscillation modal identification method based on improvement Generalized Morphological | |
CN115524086A (en) | Statistical moment-curvature beam type bridge damage identification method based on axle coupling vibration | |
CN104390795B (en) | Method and device for analyzing vibration modes of rail vehicles | |
CN110184885A (en) | A method of surface evenness is tested based on smart phone | |
CN116698316A (en) | Recognition Method of Bridge Vibration Shape Based on Vehicle Vibration Signal | |
CN115200700A (en) | Modal parameter identification method based on Welch method and covariance random subspace method | |
CN106482920A (en) | A kind of generation method of the vibration environment specification spectrum for motor train unit converter | |
CN116842775A (en) | A driving identification method for bridge modal shapes that eliminates the influence of road surface roughness | |
CN115032894A (en) | A high-speed train suspension fault diagnosis method based on T-S fuzzy data-driven ToMFIR | |
CN114264727A (en) | A method for damage identification of track-bridge system based on dynamic response of operating trains | |
CN110399683B (en) | Bridge impact coefficient extraction method based on frequency domain amplitude spectrum similarity filtering technology | |
Yang et al. | Bridge damage identification based on synchronous statistical moment theory of vehicle–bridge interaction | |
CN116558745A (en) | Bridge frequency extraction method under combined action of mobile and non-mobile vehicles | |
Kong et al. | Vehicle response-based bridge modal identification using different time–frequency analysis methods | |
CN116678576A (en) | Bridge vibration mode correction system based on vehicle vibration signals | |
Bouferrouk | Methods of calculating aerodynamic force on a vehicle subject to turbulent crosswinds | |
Liu et al. | Smartphone-based bridge frequency identification using vehicle contact-point response | |
CN116972798A (en) | Bridge pavement unevenness recognition method based on limited vehicle response |
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 |