Bridge dynamic deflection and dynamic strain data fusion method

Abstract

The invention discloses a bridge dynamic deflection and dynamic strain data fusion method, and belongs to a bridge structure health monitoring technology. According to the method, deflection measuring points and strain measuring points are deployed at the bottom of a bridge according to rectangular grids established in the longitudinal and transverse directions; and then taking the span trend of the bridge as a longitudinal direction, taking a bridge cross section perpendicular to the span trend of the bridge as a position reference of a data fusion point, and determining the data fusion point based on the bridge cross section to realize fusion and mutual conversion of the strain data and the deflection data. Through data fusion, the dynamic deflection sudden change and the dynamic strain hysteresis effect are effectively associated, the structure damage feature omission ratio is remarkably reduced, and meanwhile early warning of micro damage is supported.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of bridge health monitoring, and relates to a bridge dynamic deflection and dynamic strain data fusion method. BACKGROUND

[0002] In the field of bridge health monitoring, dynamic deflection and dynamic strain are core parameters for evaluating structural dynamic response. The existing technology generally adopts a separate collection and analysis mode: dynamic deflection data is obtained by a laser displacement meter or acceleration integration method, and dynamic strain data is independently collected by a resistance strain gauge and the like. However, due to the difference in sampling frequency (up to 10-100 times) between the two types of data, the time reference deviates. The traditional interpolation correction method produces a phase error of 0.5%-2% under non-stationary vibration, which weakens the dynamic response correlation. At the same time, the existing analysis adopts a single parameter threshold warning, which fails to reveal the mechanical coupling mechanism between dynamic deflection and dynamic strain, resulting in 15%-30% of structural damage feature missing; and the traditional fusion method (such as Kalman filtering) is limited by the linear assumption, and the error increases by more than 40% in nonlinear vibration with an amplitude greater than 5 mm.

[0003] In the current solution, a high-precision synchronous collection system can alleviate the data asynchrony, but the hardware cost increases by 3-5 times and the sampling rate is limited (≤200 Hz), and the physical model driven method fails more than 60% of the model due to bridge material degradation. The above defects show that the existing technology cannot realize the deep correlation of multi-physical field data and the accurate fusion of nonlinear dynamic characteristics, which restricts the reliability and real-time performance of bridge safety evaluation.

[0004] In view of the core defects in the existing bridge dynamic deflection and dynamic strain monitoring technology, the technical problems to be solved include: 1) how to eliminate the time domain reference deviation of multi-source heterogeneous data caused by the difference in sampling frequency (10-100 times) without relying on high-cost synchronous hardware; 2) how to break through the limitation of single parameter analysis, establish a nonlinear mechanical coupling relationship model between dynamic deflection and dynamic strain, and enhance the dynamic response correlation; 3) how to build a new fusion framework that fuses structural mechanics constitutive and data-driven, overcome the defect that the error of the traditional linear fusion algorithm (such as Kalman filtering) increases sharply in nonlinear vibration, and avoid the model failure of the physical model driven method due to material degradation.

[0005] The application adopts a matrix method to reveal the mechanical coupling mechanism between dynamic deflection and dynamic strain, and accurately establishes a corresponding calculation model. SUMMARY

[0006] The application aims to provide a bridge dynamic deflection and dynamic strain data fusion method to solve the problems listed in the background and realize the accurate spatiotemporal matching of multi-physical field data, the deep correlation of nonlinear dynamic characteristics and the real-time and reliable evaluation of bridge dynamic response.

[0007] Technical scheme: A bridge dynamic deflection and dynamic strain data fusion method, deflection measuring points and strain measuring points are arranged on the rectangular grid formed in the longitudinal and transverse directions at the bottom of the bridge; the span direction of the bridge is taken as the longitudinal direction, and the bridge cross section perpendicular to the span direction of the bridge is taken as the position reference of the data fusion point, thereby determining the deflection data and strain data of different position points, and the selection and calculation of the data fusion point include:

[0008] (1) If there are strain measuring points and deflection measuring points on the bridge cross section, the strain measuring points on the bridge cross section are taken as the data fusion point, and the deflection data contained in the longitudinal section of the bridge is directly fused;

[0009] (2) If there are only strain measuring points on the bridge cross section, the strain measuring points on the bridge cross section are selected as the data fusion point, the deflection data is calculated by cubic spline interpolation with the boundary condition that the longitudinal strain of the bridge at the support is approximately 0, and then the fusion of the deflection data and the strain data is performed;

[0010] (3) If there are only deflection measuring points on the bridge cross section, the deflection measuring point position on the section is selected as the data fusion point, the curve expression is constructed by cubic spline interpolation of the measured data of several deflection measuring points on the same beam bottom, and the strain data of the data fusion point is obtained based on the quadratic derivation of the curve expression, and then the fusion of the deflection data and the strain data is performed.

[0011] In the method, for more than one deflection measuring point on the bridge cross section, when the deflection of the data fusion point on the non-deflection monitoring longitudinal section is obtained, a spatial correction term is added to the interpolation deflection expression to calculate the deflection data of the data fusion point;

[0012] The spatial correction term is determined by cubic spline interpolation method according to the deflection data of the deflection measuring points included in the matrix network, if the four measuring points corresponding to the matrix network are represented as A( , , ), B( , , ), C( , , ), D( , , ), then the spatial correction term is calculated by the following formula :

[0013] .

[0014] Further, in the method, for more than one strain measuring point existing on the bridge cross section, strain data of one of the strain measuring points is selected for data fusion.

[0015] The data fusion mode is to construct a fusion unit based on a rectangular network and a determined data fusion point, with adjacent data fusion points in the longitudinal bridge direction as two ends of the data fusion unit; then the fusion unit is expressed as a unit geometric matrix, the deflection corresponds to a bridge displacement column vector, and the strain corresponds to a bridge deformation column vector, and the positioning vectors are array integrated to form an overall geometric pseudo-matrix and a deformation pseudo-matrix for data fusion and storage.

[0016] Among them, for high-frequency data with a sampling frequency greater than 1Hz, the two ends of the data fusion unit are adjacent data fusion points in the longitudinal bridge direction, and a coordinate system is established on the two-dimensional plane at the top of the bridge, the deflection difference at the two ends of the constructed data fusion unit and the unit length are used to calculate the rotation angle, and the rotation angle and the unit length are used to calculate the strain, and the calculation formula is expressed as:

[0017]

[0018] In the formula, is the rotation angle, denotes the longitudinal bridge direction displacement of the end of the data fusion unit, denotes the longitudinal bridge direction displacement of the head of the data fusion unit, and is the representation of deflection data, denotes the elongation deformation of the bridge along the axis, which is the change representation of the strain force, and the calculation of the data fusion point is:

[0019] .

[0020] Further, the method ignores the transverse bridge direction displacement, does not establish a data fusion unit between adjacent data fusion points in the transverse bridge direction, only considers the longitudinal bridge direction dynamic strain data, ignores the torsional deformation, and considers that the change of the deflection data in the beam bottom area is continuous.

[0021] Further, the method requires that there be strain measuring points at both ends of the bridge longitudinal section where the data fusion point is located, and that strain sensors be deployed.

[0022] Beneficial effects: compared with the prior art, the method provided by the application does not rely on a complex machine learning model, but performs partition interval probability statistical modeling based on measured data characteristics, directly reconstructs a continuous time sequence through dynamic sampling, has advantages such as simplified simulation steps, high calculation efficiency, and clear physical meaning, and the specific advantages are as follows:

[0023] 1. The present invention optimizes the deployment of deflection and strain measurement points, and then reconstructs the time base of dynamic deflection-dynamic strain monitoring data based on the consideration of the time inversion algorithm. It is suitable for synchronization of low, medium and high frequency heterogeneous data, overcomes the strict requirements of traditional methods on data frequency matching, and is a universal algorithm architecture.

[0024] 2. The present invention can effectively correlate sudden changes in dynamic deflection with the hysteresis effect of dynamic strain, significantly reducing the missed detection rate of structural damage characteristics while supporting early warning of micro-damage.

[0025] 3. By adjusting the number of data points, the computing efficiency and accuracy can be flexibly balanced, and manual experience parameter intervention is supported to achieve controllable accuracy adjustment.

[0026] 4. Using software-layer data fusion instead of hardware-layer synchronization solution can reduce the overall cost of the monitoring system, while supporting real-time processing of distributed sensor networks, and has a wide range of applicable scenarios. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] Figure 1 This is a flow chart of an embodiment of the present invention;

[0028] Figure 2 It is an array diagram of the deflection and strain measurement points at the bottom of the bridge;

[0029] Figure 3 is the displacement deformation unit diagram;

[0030] Figure 4 is a set of measured deflection data obtained in the embodiment;

[0031] Figure 5 is a set of strain measurement data obtained in the embodiment;

[0032] Figure 6 is the original strain-time monitoring data in the embodiment. DETAILED DESCRIPTION

[0033] The implementation and application of the present invention will be further introduced below with reference to the accompanying drawings.

[0034] First, the implementation process of a bridge dynamic deflection and dynamic strain data fusion method provided by the present invention includes the following steps:

[0035] S1. Optimize the deployment structure of strain sensors and deflection sensors on bridges.

[0036] First, the span direction of the bridge is determined as the longitudinal direction of the bridge, and the direction perpendicular to the longitudinal direction of the bridge is the transverse direction. Then, based on the bottom of the bridge, the strain measurement points and deflection measurement points are deployed according to the network matrix. The sensors are all at the bottom of the beam, and the beam height is only calculated. Considering the same, ignoring the longitudinal displacement of the bridge, only taking the longitudinal dynamic strain data, ignoring the torsional deformation, and considering that the deflection data in the beam bottom area is continuous.

[0037] In order to increase practicability, considering that a single piece of beam often has only one deflection sensor in the same cross section, and the transverse connection condition between multiple pieces of beams in the actual application scene is difficult to determine, the embodiment does not establish a unit between the adjacent data fusion points in the transverse direction of the bridge.

[0038] S2, determination of data of different types of fusion points.

[0039] The bridge dynamic deflection and dynamic strain fusion to be solved by the application is established in a coordinate system with the bridge bottom as a two-dimensional plane, corresponding to the longitudinal direction of the bridge and the transverse direction of the bridge. Then, the bridge cross section is taken as a reference for the position of the data fusion point, and the deployment of the sensor is considered to be at the bridge bottom position, so that there may be a strain measuring point or a deflection measuring point at the tangent of the bridge cross section and the bridge bottom. The specific conditions include:

[0040] 1) In the longitudinal position, the principle of minimizing interpolation is adopted, if the longitudinal strain and deflection measuring point are in the same section, the strain monitoring position of the section is selected, and interpolation is not needed.

[0041] 2) The strain monitoring position should be selected, and the deflection is obtained by cubic spline interpolation with the natural boundary condition that the longitudinal strain of the beam at the support is approximately 0.

[0042] 3) For the selected deflection monitoring position, the strain is obtained by the cubic spline interpolation deflection curve expression , and then the second derivative of the expression is obtained , that is , is a function of the height of the beam.

[0043] 4) For the positions without deflection data and strain data, generally, the data at the position has no practical significance and will also cause a large amount of calculation, and should be discarded. However, for the application, the calculation and derivation can be performed in combination with the above 2) and 3), and it is required that strain sensors should be distributed on both sides of the bridge cross section.

[0044] After fusion, each data fusion point has deflection and strain data.

[0045] S3, strain data and deflection data fusion based on the bridge longitudinal section.

[0046] In combination with the description of the application, the application is not limited to the details of the application described above, and various modifications and changes can be made without departing from the spirit and scope of the application. Figure 2 Figure 2 ​The strain measuring point and deflection measuring point arrangement structure of the bottom of the two bridges are included, on the bottom of the upper bridge, the bridge longitudinal strain measuring point and deflection measuring point are located in the same bridge longitudinal section, then are directly fused, on the bottom of the lower bridge, only the deflection measuring point is arranged on the bridge longitudinal section, and the data fusion point needs to be inferred from the deflection data.

[0047] In step S3, if there are multiple deflection measuring points in the transverse direction of the same longitudinal section, when obtaining the deflection of the data fusion point on the non-deflection monitoring longitudinal section, a spatial correction term needs to be added to the interpolation deflection expression The spatial correction term determines the dynamic deflection and dynamic strain data of the data fusion point by using the interpolation method of the surrounding sensors:

[0048] Suppose that the four corner points of the rectangular grid where the target point is located are measured (A, B, C, D)

[0049] The four measuring points are represented as A( , , ), B( , , ), C( , , ), D( , , ), then the spatial correction term is calculated by the following formula :

[0050] .

[0051] In the actual scene, the end of the box girder of the bridge is generally not provided with measuring points, but considering the replacement of new and old sensors, there can be multiple sensors on the bridge cross section, and thus for the present application, the following considerations exist:

[0052] In other embodiments, for high-frequency dynamic deflection data, frequent interpolation operation has too large computational load, therefore for the case of dense deflection measuring points, the following method can be used to improve the interpolation operation time interval according to the measured data, and the strain is solved by using the matrix method within the time interval, specifically:

[0053] Each displacement degree of freedom can establish an equation, the total number of node degrees of freedom is N, each unit has three unknown deformation parameters, and the total number of unknown deformation parameters M of the entire structure is 3 times the number of units.

[0054] In the global coordinate system, the deflection data in the data fusion unit is represented by the column displacement vector , e represents the data fusion unit number, and the longitudinal bridge displacement is approximately considered as 0, Vertical displacement of the bridge, For rotation, only in the vertical plane containing the unit axis, which can be integrated from the strain curve, there is the following expression:

[0055]

[0056] and respectively represent the column displacement vectors at the two ends of the data fusion unit, corresponding to the longitudinal displacement of the two ends of the bridge, and respectively correspond to the longitudinal displacement of the first and last ends of the data fusion unit. and respectively correspond to the first and last ends of the data fusion unit.

[0057] The corresponding deformation vector in the data fusion unit is represented as , , is the elongation deformation along the axis, which can be derived from the deflection interpolation curve, is the rotation of the first end of the data fusion unit relative to the new axis, is the rotation of the last end relative to the new axis. Let be 0 or derived from the deflection interpolation curve, and the total rotation is obtained by solving the equation set. The strain of the data fusion point is calculated as shown in Figure 3 .

[0058] Each longitudinal and transverse bridge adjacent data fusion point is considered as the two ends of the unit, then the vertical displacement and axial isometric deformation of the two ends of the unit are known, and the remaining elements are calculated according to the matrix, is the displacement vector of the unit rod end, is the unit geometry matrix, ; the unit geometry matrix equation in the global coordinate system is .

[0059] Deflection integration: the calculated single element displacement column vector is integrated in array (the k end matrix of the unit displacement vector is collected, (k=1,2, e is the unit number), and the integration rule is as follows: , (k, e) is the unit end positioning vector, k is the local code, and e is the global code).

[0060] Deflection integration: the calculated single element deformation column vector is integrated in array, and the integration rule is as follows: , (e) is the positioning vector, e is the global code, The whole deformation pseudo-matrix has the same column number as the cell number, so that the full-bridge strain is stored in a concentrated manner.

[0061] Combining Figure 2 and Figure 4 , the longitudinal position is determined, and the principle of minimizing interpolation is followed. If the longitudinal strain and deflection measuring points are in the same section, the strain monitoring position of the section is selected without interpolation, such as the upper half e = 1, 3, 5 and k = 2 section. Then the data of the measured points in Figure 4 is stored in the space corresponding to the data fusion point, as shown in Figure 5 .

[0062] For the case of lacking measured data, Table 1 shows the deflection measurement at a certain time. The deflection is obtained by cubic spline interpolation with the natural boundary condition that the longitudinal strain of the beam at the support is approximately 0, as shown in Figure 6 , such as the upper half e = 1, 3, 5 and k = 1 and e = 6 and k = 2 section.

[0063] Table 1. Deflection measurement Beam span position (cm) 0 4625 9250 13875 18500 Deflection (cm) 0 -4.14 -11.67 -5.01 0

[0064] As shown in Table 2, the strain is obtained by twice differentiating the cubic spline interpolation deflection expression, such as the lower half e = 1, 3, 5 and k = 2 section.

[0065] Table 2. Strain obtained by twice differentiation Beam span position (cm) 0 4625 9250 13875 18500 Deflection (με) 0 42 51 46 0

[0066] For the lower half e = 1, 3, 5 and k = 1 and e = 6 and k = 2 section, both the deflection and strain measured data are lacking, and the data of Table 1 and Table 2 are obtained by cubic spline interpolation deflection expression and its twice differentiation.

[0067] For the lower half e = 2, the measured data at a certain time is shown in Table 3.

[0068] Table 3. Strain obtained by twice differentiation Length l (cm) k = 1 deflection v1 (cm) k = 2 deflection v2 (cm) 2312 -3.6 -6.8

[0069] , .

Claims

1. A bridge dynamic deflection and dynamic strain data fusion method, characterized in that: Deflection and strain measurement points are deployed in a rectangular grid at the bottom of the bridge in the longitudinal and transverse directions. The span of the bridge is used as the longitudinal direction, and the cross section of the bridge perpendicular to the span of the bridge is used as the reference for the location of the data fusion point. The deflection and strain data at different locations are determined from this. The selection and calculation of the data fusion points include: (1) If there are strain measurement points and deflection measurement points on the bridge cross section, the strain measurement points on the bridge cross section are used as data fusion points to directly fuse the deflection data contained in the bridge cross section; (2) If there are only strain measurement points on the bridge cross section, the strain measurement points on the bridge cross section are selected as data fusion points. The deflection data are calculated by cubic spline interpolation under the boundary condition that the longitudinal strain of the bridge at the support is approximately 0. Then, the deflection data and strain data are fused. (3) If there is only a deflection measurement point on the cross section of the bridge, the position of the deflection measurement point on the cross section is selected as the data fusion point. The curve expression is constructed by cubic spline interpolation of the measured data of several deflection measurement points on the same beam bottom, and the strain data of the data fusion point is obtained by performing a quadratic derivative based on the curve expression. Then, the deflection data and strain data are fused.

2. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 is characterized in that: In this method, when there are more than one deflection measurement points on the bridge cross section, when obtaining the deflection corresponding to the data fusion point on the non-deflection monitoring cross section, the deflection data of the data fusion point is calculated by adding a spatial correction term through the interpolation deflection expression; The spatial correction term is determined by using the cubic spline interpolation method based on the deflection data of the deflection measuring points included in the matrix network. If the four measuring points corresponding to the matrix network are represented as A( , , ), B( , , ), C( , , ),D( , , ), the spatial correction term is calculated by the following formula : 。 3. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 is characterized in that: In this method, if there are more than one strain measuring points on the bridge cross section, the strain data of one of the strain measuring points is selected for data fusion.

4. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 is characterized in that: The data fusion method is to construct a fusion unit based on a rectangular network and determined data fusion points, with the adjacent data fusion points in the longitudinal bridge direction as the two ends of the data fusion unit; then the fusion unit is represented as a unit geometry matrix, the deflection corresponds to the bridge displacement column vector, the strain corresponds to the bridge deformation column vector, and the array is integrated according to the positioning vector to form an overall geometric pseudo-matrix and a deformation pseudo-matrix for data fusion and storage.

5. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 or 4, characterized in that: For high-frequency data with a sampling frequency greater than 1 Hz, the adjacent data fusion points in the longitudinal direction of the bridge are used as the two ends of the data fusion unit, and the coordinate system is established with the top of the bridge as the two-dimensional plane. The deflection difference at both ends of the constructed data fusion unit and the unit length are used to calculate the angle, and the rotation angle and unit length are used to calculate the strain. The calculation formula is expressed as follows: ， Where, For corners, represents the longitudinal bridge displacement of the end of the data fusion unit, Indicates the longitudinal bridge displacement of the head end of the data fusion unit, and is a representation of the deflection data, Indicates the elongation deformation of the bridge along the axis, which is a representation of the change in strain force, and then obtains the data fusion point The calculation is: 。 6. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 or 4, characterized in that: The method ignores transverse displacement, does not establish data fusion units between adjacent data fusion points in the transverse direction, only considers longitudinal dynamic strain data, ignores torsional deformation, and assumes that the change of deflection data in the beam bottom area is continuous.

7. The bridge dynamic deflection and dynamic strain data fusion method according to claim 1 or 4, characterized in that: This method requires that strain measurement points exist before and after the longitudinal position of the bridge cross section where the data fusion point is located, and strain sensors are deployed.

Images