CN111781001A - Bridge damping ratio identification method based on axle coupling - Google Patents

Bridge damping ratio identification method based on axle coupling Download PDF

Info

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
mode
acceleration response
axle coupling
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
Application number
CN202010678822.0A
Other languages
Chinese (zh)
Other versions
CN111781001B (en
Inventor
阳洋
沈小俊
唐钰昇
袁爱鹏
周武召
李卫东
徐润琴
李铁军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Traffic Engineering Quality Inspection Co ltd
Chongqing Transportation Planning And Technology Development Center Chongqing Transportation Engineering Cost Station
Chongqing University
Original Assignee
Chongqing Traffic Engineering Quality Inspection Co ltd
Chongqing Transportation Planning And Technology Development Center Chongqing Transportation Engineering Cost Station
Chongqing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing Traffic Engineering Quality Inspection Co ltd, Chongqing Transportation Planning And Technology Development Center Chongqing Transportation Engineering Cost Station, Chongqing University filed Critical Chongqing Traffic Engineering Quality Inspection Co ltd
Priority to CN202010678822.0A priority Critical patent/CN111781001B/en
Publication of CN111781001A publication Critical patent/CN111781001A/en
Application granted granted Critical
Publication of CN111781001B publication Critical patent/CN111781001B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/007Subject matter not provided for in other groups of this subclass by applying a load, e.g. for resistance or wear testing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical 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

The invention belongs to the technical field of building structure damage identification, and particularly relates to a bridge damping ratio identification method based on axle coupling, which comprises the following steps: the method comprises the following steps that firstly, acceleration response signals of each point of a bridge are collected by a moving test vehicle according to a certain sampling frequency; step two, filtering the acceleration response signal of the step one through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge; step three, assuming a damping ratio and the original frequency omega of the bridge1Knowing, an exponential function can then be obtained
Figure DDA0002585036780000011
Then, the filtered acceleration response of the second step is divided by the exponential function, and the filtered acceleration response of the undamped bridge can be obtained; step four, passing shortAnd (3) acquiring the mode of the signal processed in the step three by using a time-frequency domain decomposition method, and judging whether the assumed damping ratio is the real damping ratio by judging whether the maximum value of the mode is in the middle point of the mode. Compared with the prior art, the method has the advantages that the identification result of the bridge damping ratio is better, and the operation is convenient.

Description

Bridge damping ratio identification method based on axle coupling
Technical Field
The invention belongs to the technical field of building structure damage identification, and particularly relates to a bridge damping ratio identification method based on axle coupling.
Background
The bridge is used as an important component of national traffic engineering and plays a vital role in the development of socioeconomic development of China. If the bridge fails accidentally for some reason, property loss and even serious casualties can be caused. Therefore, it is very important to know the damage characteristics of the bridge from the perspective of structural safety, and the bridge health monitoring method at the present stage is to directly arrange the sensors on the bridge to be measured to obtain bridge parameters, which is called as a direct measurement method.
Based on the defects of direct testing, Yang's is equal to the mobile intelligent bridge detection method proposed in 2004, that is, a test vehicle runs on a bridge, and based on signal processing acquired by sensors on the test vehicle, bridge parameters are identified until an indirect measurement idea of identifying damage is to be carried out. After the students go forward, the method is applied to parameter identification of bridge frequency, mode, damping and the like, and finally bridge damage is identified. Yang et al extracted the most fundamental parameter, frequency, of the bridge based on the study of axle coupling. Empirical mode decomposition (EDM) is introduced into the team, Inherent Modal Functions (IMFs) are generated, finally, frequency is successfully identified through Fast Fourier Transform (FFT), then, frequency identification effects are improved through introducing singular spectrum analysis band-pass filtering and the like, meanwhile, based on the adverse effects of roughness on the identification effects and the characteristics of a mobile intelligent vehicle measurement method, yang et al propose that two test vehicles pass through the same road surface twice, and then signal subtraction is carried out to reduce the influence of the roughness of the road surface, in addition, the research team proposes that bridge frequency is extracted through axle contact point signals, and the result shows that the contact point responses can better extract the bridge frequency. The researcher then started to attempt to identify another parameter, the mode, of the bridge. In 2014, a Yang professor team successfully extracts a bridge mode by constructing instantaneous amplitude through a vehicle body response through Hilbert Transform (HT), in the same year, Obbrien and Malekjafarian build the bridge mode step by step from local to whole segmentation by combining concepts of singular value decomposition and short-time frequency domain decomposition based on finite element simulation of an axle coupling system, Li, Au and the like improve a mode identification method provided by Yang and the like, and excitation is applied to the vehicle body to increase the bridge response and improve the anti-noise property, so that the mode identification effect is improved. Bridge damping-another important parameter in bridge health monitoring, many studies indicate that damping may be more sensitive to bridge damage. Obrien and coworkers propose a method for detecting bridge damping changes using a detection vehicle-trailer vehicle system. This method requires extensive and repetitive simulations, and lacks a clear theoretical basis. Gonzalez et al propose a six-step algorithm for identifying bridge damping by vehicle body response of a two-degree-of-freedom detection vehicle, and accurately identify the bridge damping by adopting an iterative method. The method comprises the steps of bin setting and the like, obtaining contact point signals through signals measured by an acceleration sensor and a displacement sensor which are installed on a two-axis mobile test vehicle, and identifying the damping ratio of a bridge, wherein the contact point signals are lack of experimental verification.
The above describes the identification process of each parameter of the bridge, but the identification method of the damping ratio is not mature yet.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides the bridge damping ratio identification method based on the axle coupling, the identification result is better, and the operation is convenient.
In order to solve the technical problems, the invention adopts the following technical scheme:
the bridge damping ratio identification method based on axle coupling comprises the following steps:
the method comprises the following steps that firstly, acceleration response signals of each point of a bridge are collected by a moving test vehicle according to a certain sampling frequency;
step two, filtering the acceleration response signal of the step one through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge;
step three, assuming a damping ratio and the original frequency omega of the bridge1Knowing, an exponential function can then be obtained
Figure BDA0002585036760000023
Then, the filtered acceleration response of the second step is divided by the exponential function, and the filtered acceleration response of the undamped bridge can be obtained;
and step four, acquiring the mode of the signal processed in the step three through a short-time frequency domain decomposition method, judging whether the assumed damping ratio is the real damping ratio by judging whether the maximum value of the mode is in the mode midpoint, if so, taking the assumed damping ratio as the real damping ratio of the bridge, otherwise, re-assuming a damping ratio until the maximum value of the mode is in the mode midpoint, and finally identifying the damping ratio of the bridge.
Further, the upper limit value of the band-pass filter in the step two is
Figure BDA0002585036760000021
Further, step two is a band-pass filterHas a lower limit value of
Figure BDA0002585036760000022
Further, the sampling frequency of the test vehicle in the step one is 100 Hz.
Further, the test vehicle in the step one is a single degree of freedom vehicle.
Further, the mass of the test car in the step one is 1000 kg.
Further, the speed of the test vehicle in the first step is constantly 1 m/s.
Compared with the prior art, the method has the advantages that the identification result of the bridge damping ratio is better, and the operation is convenient.
Drawings
FIG. 1 is a flow chart of a bridge damping ratio identification method based on axle coupling according to the present invention;
FIG. 2 is a simplified model diagram of a dual-axle system based on the axle coupling bridge damping ratio identification method of the present invention;
FIG. 3 is a schematic diagram of the numbering of units and nodes in a bridge model based on the bridge damping ratio identification method of axle coupling according to the present invention;
FIG. 4 is a time domain signal diagram acquired by a test vehicle sensor based on the axle coupling bridge damping ratio identification method of the present invention;
FIG. 5 is a time domain signal diagram of the axle coupling-based bridge damping ratio identification method after filtering the collected signals;
FIG. 6 is a filtered time domain signal plot of the bridge damping ratio identification method based on axle coupling divided by the attenuation signal according to the present invention;
FIG. 7 is a fitting mode diagram of the bridge damping ratio identification method based on axle coupling according to the present invention;
FIG. 8 is a fitting mode diagram after symmetric processing according to the bridge damping ratio identification method based on axle coupling of the present invention;
FIG. 9 is a graph of the identification results of different assumed damping ratios under the A-level roughness of the bridge damping ratio identification method based on axle coupling of the present invention;
FIG. 10 is a graph of the identification results of different assumed damping ratios under B-level roughness for the axle coupling-based bridge damping ratio identification method of the present invention;
FIG. 11 is a graph of the damping ratio recognition result at 50db for the axle coupling based bridge damping ratio recognition method of the present invention;
FIG. 12 is a graph of the damping ratio identification result at 40db for the axle coupling based bridge damping ratio identification method of the present invention;
FIG. 13 is a graph of the damping ratio identification result at 30db for the axle coupling-based bridge damping ratio identification method of the present invention.
Detailed Description
In order that those skilled in the art can better understand the present invention, the following technical solutions are further described with reference to the accompanying drawings and examples.
As shown in FIG. 1, the bridge damping ratio identification method based on axle coupling of the present invention comprises the following steps:
the method comprises the following steps that firstly, acceleration response signals of each point of a bridge are collected by a moving test vehicle according to a certain sampling frequency;
step two, filtering the acceleration response signal of the step one through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge;
step three, assuming a damping ratio and the original frequency omega of the bridge1Knowing, an exponential function can then be obtained
Figure BDA0002585036760000033
Then, the filtered acceleration response of the second step is divided by the exponential function, and the filtered acceleration response of the undamped bridge can be obtained;
and step four, acquiring the mode of the signal processed in the step three through a short-time frequency domain decomposition method, judging whether the assumed damping ratio is the real damping ratio by judging whether the maximum value of the mode is in the mode midpoint, if so, taking the assumed damping ratio as the real damping ratio of the bridge, otherwise, re-assuming a damping ratio until the maximum value of the mode is in the mode midpoint, and finally identifying the damping ratio of the bridge.
As a preferred embodimentIn the second step, the upper limit value of the band-pass filter is
Figure BDA0002585036760000031
Preferably, the lower limit value of the band-pass filter in the second step is
Figure BDA0002585036760000032
Preferably, the sampling frequency of the test vehicle in the first step is 100 Hz.
As the preferred scheme, the test vehicle in the step one is a single-degree-of-freedom vehicle.
Preferably, the mass of the test car in the step one is 1000 kg.
Preferably, the speed of the test vehicle in the first step is constantly 1 m/s.
In fig. 2, a test car and a towing vehicle having a distance of △ L are driven at a constant speed v over the bridge deck, the test car and the towing vehicle being simplified to a moving spring mass m supported thereonv,2And mv,1A damper with a damping coefficient cv2, cv1 and a spring with a stiffness kv2, kv 1. The bridge is simply supported by length L, mass m per unit length, and bending stiffness EI (n order damping ratio), where the bending stiffness EI includes the function of non-structural members, such as sidewalks in actual bridges of bridge railings and bridge decks and mass m per unit length, can be obtained from actual bridge design data or estimated from cross-sectional area steps. The bridge is assumed to be stationary before the test car arrives.
The equation of motion for this axle coupling system can be written as:
Figure BDA0002585036760000041
Figure BDA0002585036760000042
Figure BDA0002585036760000043
wherein u (x, t) represents the vertical displacement of the bridge structure from the left support point x, uv1(t),uv2(t) time t and the test carriage and tractor, respectively, from their static equilibrium positionsa=t+ΔL/v=t+tsThe static displacement measured is started, it being noted that time t is the time after the vehicle enters the bridge. The dotted characters represent the derivative with respect to time t and coordinate x. Interaction force f of contact points of test vehicle and tractor and platformc1(t),fc2(t) can be written as:
Figure BDA0002585036760000044
Figure BDA0002585036760000045
g represents the acceleration of gravity
The bridge displacement can be expressed by generalized coordinates qn (t) and modal sin (n pi x/L) of the simply supported beam as:
Figure BDA0002585036760000046
assuming that the vehicle mass mv1, mv2 is much smaller than the mass of the deck plate, mv1<<m*L and mv2<<m*L, this assumption is easy to implement for an actual bridge. By substituting equation (6) into equation (3), multiplying by sin (n π x/L), and integrating from 0 to L, then according to the quadrature condition of the sine function, the nth mode balance equation of the structure can be written as:
Figure BDA0002585036760000047
ωnis the nth order angular frequency of the bridge deck
Figure BDA0002585036760000048
For the zero initial condition, the generalized coordinate q of the bridge can be obtained from equation (7)n(t) is:
Figure BDA0002585036760000051
substituting equation (9) into equation (6) yields a bridge displacement of:
Figure BDA0002585036760000052
the displacement u of the test vehicle can also be obtained by substituting the equation (10) into the equation (1)v1(t)。
In practice, the response component associated with the nth modal frequency of the bridge may be separated from the response of the test car. In this study, the upper and lower limits were used
Figure BDA0002585036760000061
And
Figure BDA0002585036760000062
the band pass filter of (1). 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).
Substituting the last term of equation (10) into equation (1) yields a vehicle displacement associated with the nth modal shape of the deck (11) as:
Figure BDA0002585036760000063
the damping ratio of the test vehicle in the formula is ξv1=cv1/(2mv1ωv1) Furthermore, in addition to
Figure BDA0002585036760000064
In the modified direct stiffness method, only a single vibration mode of the deck slab is required for damage identification. Since it is not guaranteed that the exact high frequency vibration mode is obtained in the field, the frequencies and modes in the following discussion refer to the first vibration mode of the bridge deck, unless otherwise stated. The response component R1 of the test car associated with the first vibration mode of the deck can also be written as:
Figure BDA0002585036760000065
the coefficients a1 to a2 may be determined by comparing the terms in equations (9) and (10),
Figure BDA0002585036760000066
Figure BDA0002585036760000067
the corresponding acceleration response component of the first vibration mode of the bridge deck may also be obtained by:
Figure BDA0002585036760000071
observation formula (16): when ξ1When the content is equal to 0, the content,
Figure BDA0002585036760000072
the equation is the acceleration response component of the first vibration mode of the undamped bridge deck. The relationship between them is
Figure BDA0002585036760000073
Namely, the ratio of the acceleration response obtained by the damped bridge to the acceleration response of the undamped bridge is an exponential function. By utilizing the relation, a method for identifying the damping ratio can be provided.
To verify the correctness of the above theory, the feasibility of the method was verified under single-vehicle no roughness. Numerical simulation is performed using finite elements as follows,to verify the reliability of the method. The numerical simulation of this time determines the length of bridge 30 meters and the mass m of unit length*19116kg, and a section area A of 7.965m2The section inertia moment Ix is 2.959m4, the bridge elastic modulus E is 2.9 × 1010N/m2, the mass of the test vehicle is 1000kg, the vehicle speed is constantly 1m/s, the sampling frequency is 100Hz, and the real damping ratio of the bridge is assumed to be 0.01.
When a short-time frequency domain decomposition method is used for extracting a mode, a bridge is divided into 10 units, namely E1-E10, as shown in FIG. 3, the rest numbers are unit node numbers (j is 1,2, … and 11), it is noted that the nodes of the mode obtained by the method are 2-10, and the nodes 1 and 11 are just the positions of an incoming bridge and an outgoing bridge, and the position vibration is very weak and can be almost ignored due to the constraint of a support, so the method does not consider the node mode values of the edge units.
FIG. 4 is a graph of acceleration time domain signals acquired by a test vehicle driving on a bridge at a sampling frequency of 100Hz (in this case, the signals contain vehicle frequency and bridge frequency information), and then filtered by a band-pass filter to acquire acceleration signals of a first-order frequency of the bridge. Referring to FIG. 5, assuming a damping ratio, and assuming a true damping ratio of 0.01, the damping coefficient is
Figure BDA0002585036760000075
The filtered signal is divided by the attenuation coefficient to obtain a signal diagram in fig. 6, and a short-time frequency domain decomposition method is adopted to obtain a modal signal in fig. 7, so that the observation of the signal in fig. 6 shows that compared with the non-bridge damping, the front half section is basically consistent, the backward signal result is worse, which also causes the distortion phenomenon to occur at the tail end of the signal in fig. 7, because the attenuation effect is stronger when the bridge is damped to the rear half section of the bridge, the bridge frequency signal is possibly attenuated and almost exhausted, so the reduction effect after the division is not good, at this time, the modal in fig. 8 is obtained by adopting the method of partial signal symmetry, the inventor obtains the result in fig. 8 through a large number of numerical simulations, when the assumed damping ratio is 0.01 (at this time, the true damping ratio is exactly the true damping ratio), the maximum value of the modal is centered, when the assumed damping ratio is 0.0092 (at this time, is smaller than the true damping ratio), the maximum value of the modal is deviated to the left, maximum of this modeThe value is deviated to the right, namely the deviation to the right occurs, when the damping ratio is assumed to be between 0.0092 and 0.0107, the phenomenon of deviation to the right from the left is not obvious at the moment, and the phenomenon can be taken as the real damping ratio, and the maximum error of the method is that
Figure BDA0002585036760000074
Within the error range. The recognition result is good under the condition that the bicycle is not interfered by external signals. However, various external influences still exist in the actual situation.
The above is a conclusion that the inventor can identify the applicable range of the damping ratio based on the condition that the real damping ratio is 0.01, and the inventor has obtained through a large number of numerical simulations, when the real damping ratio is less than or equal to 0.02, the method can effectively identify the damping ratio, and the reason that the real damping ratio is greater than 0.02 and cannot be identified is as described above, when the real damping ratio is greater, the bridge signal is attenuated faster, that is, the mode identification effect at the rear end of the signal is poor. In conclusion, the method is suitable for small damping identification which is less than or equal to 0.02.
Influence on roughness
The actual bridge pavement has unevenness, so the roughness of the invention adopts ISO 8608(1995) standard [36 ]]The proposed function density function is modeled. Function of density of its function Gd(n) is as follows:
Figure BDA0002585036760000081
wherein n is a spatial frequency per unit length, w is a constant of 2, and n0Is 0.1cycle/m, Gd(n0) The function value of the displacement functional density is determined by the roughness grade of the road surface. The square root of the geometric mean of the function values provided by the ISO 8608(1995), i.e. the function G of the roughness-displacement functional density of the pavement at each leveld(n0) The values are respectively:
a level: gd(n0)=4×10-6m3(ii) a B stage: gd(n0)=8×10-6m3(ii) a C level: gd(n0)=16×10-6m3
The displacement amplitude value d of the roughness of the road surface under each level of roughness can be expressed as:
Figure BDA0002585036760000082
where Δ n is the sampling interval of the spatial frequency
Then, cosine functions with different spatial frequencies are superposed to simulate the road roughness r (x), which can be expressed as:
Figure BDA0002585036760000083
in the formula, ns,iSpatial frequency of roughness, di、θiThe amplitude of the roughness and the random phase angle, respectively.
In order to study the identification condition of the bridge damping ratio under the influence of roughness, the numerical simulation adopts a method that two test vehicles pass through the same road surface twice and then the signal subtraction is carried out to reduce the influence of the road surface roughness, which is proposed by Yangmen et al. The mass of the large test vehicle body is 2000kg, the rigidity is 20000N/m, the mass of the small test vehicle body is 1000kg, and the rigidity is 10000N/m. And the real damping ratio of the bridge is assumed to be 0.01, and then a series of damping ratios are assumed and the method is applied to identify the real damping ratio. The recognition results are shown in fig. 9-10.
It can be seen from fig. 9 and 10 that at a roughness level of A, B, when the damping ratio is assumed to be exactly the true damping ratio, or in the vicinity of the true damping ratio, the mode maximum is at the mode midpoint, and when the damping ratio is assumed to deviate from the true damping ratio, the mode maximum is deviated to a different degree from the mode midpoint, i.e., the maximum is not at the mode midpoint. As shown in fig. 9, under B-class roughness, when the damping ratio is assumed to be in the range of 0.0095 to 0.0103, the mode maximum is at the mode midpoint, when the damping ratio is assumed to be less than 0.0095, the left-hand yaw phenomenon occurs, and when the damping ratio is assumed to be greater than 0.0103, the right-hand yaw phenomenon occurs, which is the phenomenon that the mode maximum is at the mode midpointThe method can help how to assume the damping ratio so as to identify the real damping ratio more effectively, and the maximum error is
Figure BDA0002585036760000091
And the error range is met.
The identification result of the damping ratio is not ideal under C-class roughness, because the roughness is too large, the effect of reducing the roughness by signal subtraction is not ideal, so that roughness elimination is not ideal, and moreover, too large roughness also causes that back-end signal acquisition is not ideal (which is also the reason of partial symmetrical processing in the introduction of the method), so that an ideal modal signal is not obtained, and the damping ratio is not ideal.
Concerning the influence of noise
In practical application, signals acquired by an acceleration sensor can not be interfered by noise, in order to explore the effectiveness of the damping ratio identification method under the noise interference, the noise immunity of the method is researched by a method of adding Gaussian white noise to acceleration signals acquired in numerical simulation, and the signal-to-noise ratio is used as an index, and is defined as follows:
Figure BDA0002585036760000092
in the formula: n is the number of data points, yiAcceleration response, σ, of test vehicle containing noise at time iiThe SNR is the signal-to-noise ratio (dB) which is the noise value at the i-th time, and the larger the SNR is, the smaller the noise influence is, the lower the interference degree of the signal is, and the smaller the SNR is, the larger the noise influence is, the larger the interference degree of the signal is.
The present invention assumes a true damping ratio of 0.01 and then assumes a damping ratio of 0.01 plus noise interference, thereby analyzing the error of the recognition result. The simulation is set to research the influence of noise under the A-level roughness, numerical simulation is carried out on the noise of each level, and the modal shape identification is obtained by taking the value of the numerical simulation.
The simulation draws a damping ratio identification model under 3 different noise levelsThe state results are shown in fig. 11, 12 and 13: when the signal-to-noise ratio is 50db, when the assumed damping ratio is 0.091-0.0108, the mode maximum value is at the mode midpoint, when the assumed damping ratio is less than or equal to 0.0091, the mode maximum value is not at the mode midpoint, namely, the left deviation phenomenon occurs, when the assumed damping ratio is more than 0.0108, the right deviation phenomenon occurs, and the maximum error at the moment is
Figure BDA0002585036760000093
And the error range is met. And as the signal-to-noise ratio is smaller and the noise interference is larger, the damping ratio can be well identified by the method, and the identification precision is more accurate. But when the signal-to-noise ratio is less than or equal to 20db, the noise interference is already large, and the identification result is poor, which indicates that the damping ratio identified by the method has certain noise resistance.
By combining the numerical simulation results, under the influence of roughness and noise, under the A, B-level roughness, when the signal-to-noise ratio is greater than or equal to 30db, the damping ratio is better identified by the method.
Finally, the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all of them should be covered in the claims of the present invention.

Claims (7)

1. The bridge damping ratio identification method based on axle coupling is characterized by comprising the following steps of:
the method comprises the following steps that firstly, acceleration response signals of each point of a bridge are collected by a moving test vehicle according to a certain sampling frequency;
step two, filtering the acceleration response signal of the step one through a band-pass filter to obtain the acceleration response of the first-order frequency of the bridge;
step three, assuming a damping ratio and the original frequency omega of the bridge1Knowing, an exponential function can then be obtained
Figure FDA0002585036750000013
Then, the filtered acceleration response of the second step is divided by the exponential function, and the filtered acceleration response of the undamped bridge can be obtained;
and step four, acquiring the mode of the signal processed in the step three through a short-time frequency domain decomposition method, judging whether the assumed damping ratio is the real damping ratio by judging whether the maximum value of the mode is in the mode midpoint, if so, taking the assumed damping ratio as the real damping ratio of the bridge, otherwise, re-assuming a damping ratio until the maximum value of the mode is in the mode midpoint, and finally identifying the damping ratio of the bridge.
2. The axle coupling-based bridge damping ratio identification method of claim 1, wherein: in step two, the upper limit value of the band-pass filter is
Figure FDA0002585036750000011
3. The axle coupling-based bridge damping ratio identification method of claim 2, wherein: in step two, the lower limit value of the band-pass filter is
Figure FDA0002585036750000012
4. The axle coupling-based bridge damping ratio identification method of claim 3, wherein: and in the step one, the sampling frequency of the test vehicle is 100 Hz.
5. The axle coupling-based bridge damping ratio identification method of claim 4, wherein: the test vehicle in the step one is a single degree of freedom vehicle.
6. The axle coupling-based bridge damping ratio identification method of claim 5, wherein: the mass of the test car in the step one is 1000 kg.
7. The axle coupling-based bridge damping ratio identification method of claim 6, wherein: in the first step, the speed of the test vehicle is constantly 1 m/s.
CN202010678822.0A 2020-07-15 2020-07-15 Bridge damping ratio identification method based on axle coupling Active CN111781001B (en)

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 (2)

* Cited by examiner, † Cited by third party
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
CN115436473A (en) * 2022-09-05 2022-12-06 重庆大学 Method for identifying structural damage based on timing synchronization theory

Citations (12)

* Cited by examiner, † Cited by third party
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 东南大学 Bridge based on move vehicle is quickly tested and parameter identification method
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 重庆大学 beam bridge damage identification method based on average value of acceleration ratio and GPSA (general purpose analysis System) algorithm

Patent Citations (12)

* Cited by examiner, † Cited by third party
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 东南大学 Bridge based on move vehicle is quickly tested and parameter identification method
CN110568069A (en) * 2019-09-11 2019-12-13 重庆大学 beam bridge damage identification method based on average value of acceleration ratio and GPSA (general purpose analysis System) algorithm

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
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 (4)

* Cited by examiner, † Cited by third party
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
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
CN106524967B (en) A kind of automobile core wheel actual travel displacement measurement and extracting method
CN111781001B (en) Bridge damping ratio identification method based on axle coupling
Piombo et al. Modelling and identification of the dynamic response of a supported bridge
CN109855823A (en) A method of Damage Identification of Bridge Structure is carried out using test carriage
Zhao et al. IRI estimation by the frequency domain analysis of vehicle dynamic responses
CN105000033B (en) A kind of check and evaluation system of track geometry irregularities
CN116481747A (en) Bridge modal shape driving identification method considering bridge damping influence
CN104032629A (en) On-line vertical track long wave irregularity monitoring method and system
Yang et al. Rigid-mass vehicle model for identification of bridge frequencies concerning pitching effect
Li et al. Indirect damage detection for bridges using sensing and temporarily parked vehicles
CN115524086A (en) Statistical moment-curvature beam type bridge damage identification method based on axle coupling vibration
Yang et al. Closely spaced modes of bridges estimated by a hybrid time–frequency method using a multi-sensor scanning vehicle: Theory and practice
CN116698316A (en) Bridge vibration mode identification method based on vehicle vibration signals
CN115032894A (en) High-speed train suspension fault diagnosis method based on T-S fuzzy data driving ToMFIR
Gong et al. Modal vibration decomposition method and its application on multi-mode vibration control of high-speed railway car bodies
Kong et al. Vehicle Response-Based Bridge Modal Identification Using Different Time–Frequency Analysis Methods
CN114264727A (en) Track-bridge system damage identification method based on dynamic response of operation train
CN112484839A (en) Bridge movement detection method and system based on two-axis vehicle response
CN113447283A (en) Failure detection method and device for under-train shock absorber
Liu et al. Vibration displacement measurement based on three axes accelerometer
CN116678576A (en) Bridge vibration mode correction system based on vehicle vibration signals
Liu et al. Smartphone-based bridge frequency identification using vehicle contact-point response
Shi et al. Theory for computing vehicle-rail contact responses from a multi-DOF test vehicle and detecting track modulus and rail damages
CN106679659A (en) Signal denoising method based on parameter-adjustable nonlinear track differentiator
CN114923649A (en) Vehicle body composite vibration identification method based on test and simulation technology

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