CN115728383B - Bridge structure damage positioning method, device, computer equipment and medium - Google Patents
Bridge structure damage positioning method, device, computer equipment and medium Download PDFInfo
- Publication number
- CN115728383B CN115728383B CN202211258909.8A CN202211258909A CN115728383B CN 115728383 B CN115728383 B CN 115728383B CN 202211258909 A CN202211258909 A CN 202211258909A CN 115728383 B CN115728383 B CN 115728383B
- Authority
- CN
- China
- Prior art keywords
- bridge
- principal component
- analyzed
- matrix
- dynamic vibration
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 147
- 238000012545 processing Methods 0.000 claims abstract description 57
- 238000006073 displacement reaction Methods 0.000 claims abstract description 30
- 238000000926 separation method Methods 0.000 claims abstract description 30
- 238000000513 principal component analysis Methods 0.000 claims abstract description 25
- 238000004590 computer program Methods 0.000 claims abstract description 23
- 230000004044 response Effects 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims description 14
- 239000013598 vector Substances 0.000 claims description 11
- 230000035772 mutation Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 9
- 238000005070 sampling Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000009434 installation Methods 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
The application relates to the technical field of structural damage detection, and provides a bridge structural damage positioning method, a bridge structural damage positioning device, computer equipment, a storage medium and a computer program product. The application can improve the efficiency and accuracy of the damage positioning of the bridge structure, has strong applicability and reduces the calculated amount of the damage positioning of the bridge structure. The method comprises the following steps: combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed, performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix, performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed, obtaining instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information, and determining the structural damage position of the bridge to be analyzed according to a time sequence corresponding to the instantaneous frequency.
Description
Technical Field
The present application relates to the field of structural damage detection technology, and in particular, to a method, an apparatus, a computer device, a storage medium, and a computer program product for locating damage to a bridge structure.
Background
With the construction of the bridge, the bridge provides important support for economic development. In order to avoid bridge accidents, ensure the operation safety of the bridge, timely find and repair early damage of the bridge structure is an important measure for reducing the occurrence of bridge accidents, and structural damage detection is an important measure for timely finding the damage of the bridge structure.
The traditional technology is used for positioning the damage of the bridge structure by a local damage detection method, but the efficiency of positioning the damage of the bridge structure by the technology is lower.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, an apparatus, a computer device, a computer readable storage medium and a computer program product for locating damage to a bridge structure.
In a first aspect, the application provides a method for locating damage to a bridge structure. The method comprises the following steps:
combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed;
performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix;
performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed;
obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information;
and determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
In one embodiment, after performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix, the method further includes:
determining a target principal component signal according to the principal component contribution rate corresponding to the principal component score matrix;
performing signal separation processing on target principal component signals in a principal component score matrix to obtain dynamic vibration information corresponding to a bridge to be analyzed, wherein the method comprises the following steps:
and carrying out signal separation processing on the target principal component signal through a window filtering model to obtain dynamic vibration information corresponding to the bridge to be analyzed.
In one embodiment, the target principal component signal is a first order principal component signal in a principal component score matrix;
performing signal separation processing on the target principal component signal through a window filtering model to obtain dynamic vibration information corresponding to the bridge to be analyzed, wherein the method comprises the following steps:
performing signal separation processing on the first-order main component signals through a window filtering model to obtain first-order structural vibration mode signals corresponding to the bridge to be analyzed;
and determining dynamic vibration information of the first-order modal frequency corresponding to the bridge to be analyzed according to the difference value between the first-order main component signal and the first-order structural vibration mode signal.
In one embodiment, performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix, including:
determining a covariance matrix corresponding to the array matrix;
determining a feature vector matrix corresponding to the covariance matrix;
and obtaining a principal component score matrix according to the eigenvector matrix and the array matrix.
In one embodiment, determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency includes:
obtaining a damage characteristic quantity index curve according to a time sequence corresponding to the instantaneous frequency;
determining time information of the preset object passing through the structural damage position according to the mutation position in the damage characteristic quantity index curve;
and determining the position of the structural damage according to the time information and the speed information corresponding to the preset object.
In one embodiment, obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information includes:
and performing Hilbert transform processing on the dynamic vibration information to obtain the instantaneous frequency of the bridge to be analyzed.
In a second aspect, the application further provides a device for positioning damage of the bridge structure. The device comprises:
the array matrix obtaining module is used for combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed;
the principal component score matrix obtaining module is used for carrying out principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix;
the dynamic vibration information obtaining module is used for carrying out signal separation processing on the target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed;
the instantaneous frequency obtaining module is used for obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information;
and the structural damage position determining module is used for determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
In a third aspect, the present application also provides a computer device. The computer device comprises a memory storing a computer program and a processor which when executing the computer program performs the steps of:
combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed; performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix; performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed; obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information; and determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
In a fourth aspect, the present application also provides a computer-readable storage medium. The computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed; performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix; performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed; obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information; and determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
In a fifth aspect, the present application also provides a computer program product. The computer program product comprises a computer program which, when executed by a processor, implements the steps of:
combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed; performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix; performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed; obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information; and determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
The bridge structure damage positioning method, the bridge structure damage positioning device, the computer equipment, the storage medium and the computer program product combine the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed, perform principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix, perform signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed, obtain instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information, and determine the structure damage position of the bridge to be analyzed according to a time sequence corresponding to the instantaneous frequency. According to the scheme, displacement sensors are arranged on a bridge to be analyzed or at the bottom of the bridge to be analyzed at equal intervals, the installation direction is perpendicular to the bridge deck direction, vertical displacement response information, sent by each displacement sensor, when a measured vehicle load passes through the bridge at a constant speed is received, all the vertical displacement response information is combined to obtain an array matrix corresponding to the bridge to be analyzed, main component analysis processing is carried out on the array matrix to obtain a corresponding main component score matrix, signal separation processing is carried out on target main component signals in the main component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed, instantaneous frequency of the bridge to be analyzed is obtained according to the dynamic vibration information, a damage characteristic index curve is drawn according to a time sequence corresponding to the instantaneous frequency, the time when the vehicle passes through the damage position is determined according to the mutation position of the damage characteristic index curve, and the structural damage position of the bridge to be analyzed is determined according to the time when the vehicle speed is multiplied by the vehicle passes through the damage position, so that the efficiency and accuracy of bridge structural damage positioning are improved, applicability is high, and the calculated amount of structural damage positioning of the bridge is reduced.
Drawings
FIG. 1 is a flow chart of a method for locating damage to a bridge structure in one embodiment;
FIG. 2 is a flow chart illustrating the steps for determining dynamic vibration information in one embodiment;
FIG. 3 is a flowchart illustrating a step of obtaining a corresponding principal component score matrix according to one embodiment;
FIG. 4 is a flow chart illustrating a step of determining a structural damage location of a bridge to be analyzed in one embodiment;
FIG. 5 is a block diagram of a bridge structure damage location device in one embodiment;
fig. 6 is an internal structural diagram of a computer device in one embodiment.
Detailed Description
The present application will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present application more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the application.
In one embodiment, as shown in fig. 1, a method for positioning damage to a bridge structure is provided, where the method is applied to a terminal or a server for illustration, and includes the following steps:
and step S101, combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed.
In this step, the bridge to be analyzed may refer to a bridge whose position of the loss of the structure to be determined; each preset position point can be arranged on the bridge or at the bottom of the bridge, and each preset position point can be arranged at equal intervals; the vertical displacement response information can be collected by displacement sensors, and the installation direction of each displacement sensor can be perpendicular to the direction of the bridge deck (bridge).
Specifically, displacement sensors are arranged on a bridge to be analyzed or at the bottom of the bridge to be analyzed at equal intervals, the installation direction is perpendicular to the direction of a bridge deck, vertical displacement response information when a vehicle load passes through the bridge at a constant speed is acquired by the displacement sensors, and the vertical displacement response information is combined to obtain a displacement sensor array matrix W corresponding to the bridge to be analyzed m×n ，W m×n Can be expressed as
Where n=1, 2, … …, N is the signal sample length; m=1, 2, … …, M is the number of measurement points.
Step S102, performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix.
In this step, the principal component analysis processing may refer to principal component analysis computation (PCA).
Specifically, a matrix W is arranged in a matrix array m×n Principal component analysis calculation was performed to obtain PCA (W m×n )＝[V m×m ，Y m×n ，Λ m×m ]Wherein V is m×m As a characteristic vector matrix, Y m×n For a corresponding principal component score matrix, the principal component score matrix may be represented as Y m×n ＝[Y (1) ，Y (2) ……Y (n) ]Wherein Y is (1) 、Y (2) 、Y (n) Representing column vectors, Λ m×m For the eigenvalue matrix corresponding to V, Λ is a diagonal matrix of m×m, and eigenvalues are arranged from large to small.
And step S103, performing signal separation processing on the target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed.
In this step, the target principal component signal may be the first column data time series signal Y in the principal component scoring matrix Y 1(t) First column data time series signal Y 1(t) May be a first order principal component time series signal; the dynamic vibration information may be bridge first order modal frequency dynamic vibration information.
Specifically, signal separation processing is carried out on target principal component signals in the principal component score matrix, and dynamic vibration information corresponding to the bridge to be analyzed is obtained.
And step S104, obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information.
In this step, the instantaneous frequency may be a first order modal frequency of the bridge.
Specifically, according to the dynamic vibration information, the instantaneous frequency of the bridge to be analyzed is obtained.
Step S105, determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
Specifically, according to a time sequence corresponding to the instantaneous frequency, a damage characteristic quantity index curve is drawn, the time when a vehicle (vehicle load) passes through the damage position is determined according to the abrupt change position of the damage characteristic quantity index curve, the structural damage position of the bridge to be analyzed is determined according to the time when the vehicle speed is multiplied by the vehicle passes through the damage position, wherein the time when the vehicle load enters the bridge can be set to be the time of the time sequence t=0.
In the bridge structure damage positioning method, the vertical displacement response information of each preset position point of the bridge to be analyzed is combined to obtain the array matrix corresponding to the bridge to be analyzed, the main component analysis processing is carried out on the array matrix to obtain the corresponding main component score matrix, the signal separation processing is carried out on the target main component signals in the main component score matrix to obtain the dynamic vibration information corresponding to the bridge to be analyzed, the instantaneous frequency of the bridge to be analyzed is obtained according to the dynamic vibration information, and the structure damage position of the bridge to be analyzed is determined according to the time sequence corresponding to the instantaneous frequency. According to the scheme, displacement sensors are arranged on a bridge to be analyzed or at the bottom of the bridge to be analyzed at equal intervals, the installation direction is perpendicular to the bridge deck direction, vertical displacement response information, sent by each displacement sensor, when a measured vehicle load passes through the bridge at a constant speed is received, all the vertical displacement response information is combined to obtain an array matrix corresponding to the bridge to be analyzed, main component analysis processing is carried out on the array matrix to obtain a corresponding main component score matrix, signal separation processing is carried out on target main component signals in the main component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed, instantaneous frequency of the bridge to be analyzed is obtained according to the dynamic vibration information, a damage characteristic index curve is drawn according to a time sequence corresponding to the instantaneous frequency, the time when the vehicle passes through the damage position is determined according to the mutation position of the damage characteristic index curve, and the structural damage position of the bridge to be analyzed is determined according to the time when the vehicle speed is multiplied by the vehicle passes through the damage position, so that the efficiency and accuracy of bridge structural damage positioning are improved, applicability is high, and the calculated amount of structural damage positioning of the bridge is reduced.
In one embodiment, the method may further determine the target principal component signal by: determining a target principal component signal according to the principal component contribution rate corresponding to the principal component score matrix; the step S103 of performing signal separation processing on the target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed specifically includes: and carrying out signal separation processing on the target principal component signal through a window filtering model to obtain dynamic vibration information corresponding to the bridge to be analyzed.
In this embodiment, the window filtering model may be a window filtering function MAF (a moving window is added to the moving average filter).
Specifically, after principal component analysis processing is performed on the array matrix to obtain a corresponding principal component score matrix, to determine the principal component order containing most of the original data information, the principal component cumulative contribution rate CCR needs to be calculated to determine, which is defined as
Wherein CCR (CCR) i Representing the cumulative contribution rate of the first-order principal component, when CCR i Not less than 90%, the first i-order principal component can be considered to contain most of the original data information, thereby extracting core data information and removing the target of redundant information. All principal components can form a principal component score matrix, the principal component contribution rate is used as an index to measure the degree that all order modes in the vibration system participate in the vibration in an energy form, the principal component contribution rate of the vibration system is sequentially reduced from a low order to a high order, particularly the first order principal component contribution rate is highest, the vibration system is mainly in the first order vibration form and the first order principal component contains most vibration information, the degree that the first order vibration participates in the vibration in the energy form is the largest, and the energy that the modes of other orders participate in the vibration is sequentially reduced from the low order to the high order, so that a target principal component signal is determined to be the first order principal component signal in the principal component score matrix. And carrying out signal separation processing on the target principal component signal through a window filtering model to obtain dynamic vibration information corresponding to the bridge to be analyzed.
According to the technical scheme, the target principal component signals are determined, and the window filtering model is utilized to conduct signal separation processing on the target principal component signals, so that dynamic vibration information corresponding to the bridge to be analyzed is obtained, the target principal component signals can be accurately determined, the accuracy of signal separation processing is improved, more accurate dynamic vibration information is facilitated to be obtained, and the accuracy of damage positioning of the bridge structure is facilitated to be improved subsequently.
In one embodiment, as shown in fig. 2, the method may further determine dynamic vibration information of a first-order modal frequency corresponding to the bridge to be analyzed by the following steps, which specifically includes: step S201, performing signal separation processing on the first-order main component signals through a window filtering model to obtain first-order structural vibration mode signals corresponding to the bridge to be analyzed; step S202, determining dynamic vibration information of first-order modal frequencies corresponding to the bridge to be analyzed according to the difference between the first-order principal component signals and the first-order structural vibration mode signals.
In this embodiment, the target principal component signal is the first order principal component signal in the principal component score matrix.
Specifically, the moving load (vehicle load) is applied to the bridge provided with the displacement sensor, the bridge displacement signal is decomposed by using the Principal Component Analysis (PCA) derivation method, the full-bridge modal shape is extracted, the physical explanation of the bridge responding to the principal component Y under the moving load is firstly deduced, the result shows that each component of the principal component Y consists of the modal shape of the bridge and the dynamic component expressed by the vibration mode frequency, and the principal component Y serving as a time sequence in the result of PCA is finally expressed as
Wherein,
where c is the damping coefficient, v is the moving load speed, l is the length of the simply supported bridge (which may be the length of a span simply supported beam), ω n Is the nth order structural frequency, m is the number of sensors, k is the sampling number, ρ is the mass per unit length of the bridge, f is the external force, q m (t) represents the m-th order generalized coordinates, and since the actual bridge damping coefficient c is a high order small amount compared with other parameters of the formula and can be ignored, the principal component Y can be finally simplified into:
the principal component also contains a component representing the modeIs a trigonometric function with the angular frequency and a trigonometric function with the omega as the angular frequency representing the dynamic characteristic component, the physical significance of the trigonometric function is respectively the structural vibration mode and the dynamic vibration information under the modal frequency of the order, meanwhile, the filtering function specifically filters the dynamic information represented by the trigonometric function with the omega as the angular frequency, and the dynamic information is filtered to be identified>The signal of Hilbert transform is dynamic vibration information at the mode frequency, so that a structural vibration mode signal is obtained through a window filter function MAF (a moving window is added on the basis of a moving average filter), the structural vibration mode signal is subtracted from a principal component signal to obtain dynamic vibration information at the mode frequency, and the principal component score matrix first column data (first order principal component signal) Y 1 (t) after MAF smoothing, obtaining a first order structure mode signal +.>Where t is the time point and k is the size of the moving window, +.>Wherein f s For signal sampling frequency, f 1 For the fundamental frequency of the detected signal, the first-order structural vibration mode signal is removed from the first-order main component signal, so that dynamic vibration information c (t) and +.>
According to the technical scheme, the dynamic vibration information of the first-order modal frequency corresponding to the bridge to be analyzed is determined according to the difference value between the first-order main component signal and the first-order structure vibration mode signal, so that the accuracy and the efficiency of determining the dynamic vibration information are improved, and the efficiency and the accuracy of positioning the damage of the bridge structure are improved.
In one embodiment, as shown in fig. 3, the principal component analysis processing is performed on the array matrix in the step S102, and the obtaining of the corresponding principal component score matrix specifically includes: step S301, determining a covariance matrix corresponding to the array matrix; step S302, determining a feature vector matrix corresponding to the covariance matrix; step S303, obtaining a principal component score matrix according to the eigenvector matrix and the array matrix.
Specifically, in a first step, a matrix W (i.e., W m×n ) Is set to the covariance matrix C of (C),
second, calculating eigenvectors and eigenvalues corresponding to the covariance matrix C, wherein V represents the eigenvector matrix, each column vector of V is the eigenvector corresponding to the covariance matrix C, Λ represents the eigenvalue matrix as a diagonal matrix, and eigenvalue λ n Arranged diagonally from large to small, the eigenvalue λ (λ i ) Can be used for calculating the principal component cumulative contribution rate CCR (CCR) i ) The method comprises the steps of carrying out a first treatment on the surface of the In the third step, the principal component Y may be represented as Y (n) ＝WV (n) Wherein Y is (n) Representing the nth order principal component, reflecting the projection of the original data onto the nth dimension vector, V (n) Representing the nth order feature vector.
According to the technical scheme, the principal component score matrix is obtained through the feature vector matrix and the array matrix, and the accuracy and the efficiency of obtaining the principal component score matrix can be improved, so that the efficiency and the accuracy of the damage positioning of the bridge structure can be improved.
In one embodiment, as shown in fig. 4, the determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency in the step S105 specifically includes: step S401, obtaining a damage characteristic quantity index curve according to a time sequence corresponding to the instantaneous frequency; step S402, determining time information of the preset object passing through the structural damage position according to the mutation position in the damage characteristic quantity index curve; step S403, determining the structural damage position according to the time information and the speed information corresponding to the preset object.
In this embodiment, the time information that the preset object passes through the structural damage position may be a time when the vehicle load (moving load) passes through the structural damage position; the speed information corresponding to the preset object may be a running speed of a vehicle load (moving load); the damage characteristic amount index curve may be a curve with time on the abscissa (i.e., time corresponding to each instantaneous frequency) and instantaneous frequency on the ordinate.
Specifically, a damage characteristic quantity index omega curve is drawn according to an instantaneous frequency omega (n) time sequence, the time when the vehicle passes through the damage position is determined according to the abrupt change position of the damage characteristic quantity index omega curve, and the structural damage position of the bridge is determined by multiplying the vehicle speed by the time when the vehicle passes through the damage position, namely converting a time axis into a spatial position axis.
For example, the speed of the moving load trolley is set to be constant v, the length of the bridge is assumed to be 6m, the moving speed of the trolley is assumed to be 0.5m/s, the sampling frequency of the sensor is 500hz, the signal sampling length=500×6/0.5=6000, and the position damage index of 2400 th sampling point is assumed to be suddenly changed, which corresponds to the damage of the bridge at the position of 0.5×2400/500=2.4=0.4L.
According to the technical scheme, the structural damage position is determined according to the time information and the speed information corresponding to the preset object, so that the efficiency and the accuracy of positioning the damage of the bridge structure are improved.
In one embodiment, the obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information in the step S104 specifically includes: and performing Hilbert transform processing on the dynamic vibration information to obtain the instantaneous frequency of the bridge to be analyzed.
In this embodiment, the Hilbert transform is the Hilbert transform.
In particular, the dynamic vibration information c (t) of the bridge under the first-order modal frequency is taken to carry out Hilbert transformation,
wherein PV represents the Cauchy principal value, thereby constructing an analytic signal z (t) as
z(t)＝c(t)+jH[(c(t))]＝a(t)e jΦ(t) ，
Wherein a (t) is an amplitude function, phi (t) is a phase function,
differentiating the instantaneous phase to obtain the instantaneous frequency
And a time series ω (t) of instantaneous frequencies can be obtained.
According to the technical scheme, the transient frequency of the bridge to be analyzed is obtained by carrying out Hilbert transform processing on the dynamic vibration information, so that the accuracy and the efficiency of obtaining the transient frequency of the bridge to be analyzed are improved, and the accuracy and the efficiency of the damage positioning of the bridge structure are improved.
It should be understood that, although the steps in the flowcharts related to the embodiments described above are sequentially shown as indicated by arrows, these steps are not necessarily sequentially performed in the order indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in the flowcharts described in the above embodiments may include a plurality of steps or a plurality of stages, which are not necessarily performed at the same time, but may be performed at different times, and the order of the steps or stages is not necessarily performed sequentially, but may be performed alternately or alternately with at least some of the other steps or stages.
Based on the same inventive concept, the embodiment of the application also provides a bridge structure damage positioning device for realizing the above related bridge structure damage positioning method. The implementation scheme of the device for solving the problem is similar to that described in the above method, so the specific limitation in the embodiments of the device for positioning damage to a bridge structure provided below may be referred to as limitation of the method for positioning damage to a bridge structure, which is not described herein.
In one embodiment, as shown in fig. 5, a bridge construction damage positioning apparatus is provided, and the apparatus 500 may include:
the array matrix obtaining module 501 is configured to combine vertical displacement response information of each preset position point of a bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed;
the principal component score matrix obtaining module 502 is configured to perform principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix;
a dynamic vibration information obtaining module 503, configured to perform signal separation processing on the target principal component signal in the principal component score matrix, to obtain dynamic vibration information corresponding to the bridge to be analyzed;
the instantaneous frequency obtaining module 504 is configured to obtain, according to the dynamic vibration information, an instantaneous frequency of the bridge to be analyzed;
and the structural damage position determining module 505 is configured to determine the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
In one embodiment, the apparatus 500 further comprises: the target principal component signal determining module is used for determining the target principal component signal according to the principal component contribution rate corresponding to the principal component score matrix; the dynamic vibration information obtaining module 503 is further configured to perform signal separation processing on the target principal component signal through a window filtering model, so as to obtain dynamic vibration information corresponding to the bridge to be analyzed.
In one embodiment, the target principal component signal is a first order principal component signal in the principal component score matrix; the dynamic vibration information obtaining module 503 is further configured to perform signal separation processing on the first-order principal component signal through the window filtering model, so as to obtain a first-order structural vibration mode signal corresponding to the bridge to be analyzed; and determining dynamic vibration information of the first-order modal frequency corresponding to the bridge to be analyzed according to the difference value between the first-order principal component signal and the first-order structural vibration mode signal.
In one embodiment, the principal component score matrix obtaining module 502 is further configured to determine a covariance matrix corresponding to the array matrix; determining a feature vector matrix corresponding to the covariance matrix; and obtaining the principal component score matrix according to the eigenvector matrix and the array matrix.
In one embodiment, the structural damage position determining module 505 is further configured to obtain a damage characteristic quantity index curve according to the time sequence corresponding to the instantaneous frequency; determining time information of the preset object passing through the structural damage position according to the mutation position in the damage characteristic quantity index curve; and determining the damage position of the structure according to the time information and the speed information corresponding to the preset object.
In one embodiment, the instantaneous frequency obtaining module 504 is further configured to perform hilbert transform processing on the dynamic vibration information to obtain an instantaneous frequency of the bridge to be analyzed.
All or part of each module in the bridge structure damage positioning device can be realized by software, hardware and a combination thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, which may be a terminal, and the internal structure of which may be as shown in fig. 6. The computer device includes a processor, a memory, a communication interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless mode can be realized through WIFI, a mobile cellular network, NFC (near field communication) or other technologies. The computer equipment also comprises an input/output interface, wherein the input/output interface is a connecting circuit for exchanging information between the processor and the external equipment, and the input/output interface is connected with the processor through a bus and is called as an I/O interface for short. The computer program when executed by the processor is used for realizing a bridge structure damage positioning method. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, can also be keys, a track ball or a touch pad arranged on the shell of the computer equipment, and can also be an external keyboard, a touch pad or a mouse and the like.
It will be appreciated by those skilled in the art that the structure shown in FIG. 6 is merely a block diagram of some of the structures associated with the present inventive arrangements and is not limiting of the computer device to which the present inventive arrangements may be applied, and that a particular computer device may include more or fewer components than shown, or may combine some of the components, or have a different arrangement of components.
In an embodiment, there is also provided a computer device comprising a memory and a processor, the memory having stored therein a computer program, the processor implementing the steps of the method embodiments described above when the computer program is executed.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored which, when executed by a processor, carries out the steps of the method embodiments described above.
In an embodiment, a computer program product is provided, comprising a computer program which, when executed by a processor, implements the steps of the method embodiments described above.
The user information (including but not limited to user equipment information, user personal information, etc.) and the data (including but not limited to data for analysis, stored data, presented data, etc.) related to the present application are information and data authorized by the user or sufficiently authorized by each party.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, database, or other medium used in embodiments provided herein may include at least one of non-volatile and volatile memory. The nonvolatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical Memory, high density embedded nonvolatile Memory, resistive random access Memory (ReRAM), magnetic random access Memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric Memory (Ferroelectric Random Access Memory, FRAM), phase change Memory (Phase Change Memory, PCM), graphene Memory, and the like. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, and the like. By way of illustration, and not limitation, RAM can be in the form of a variety of forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), and the like. The databases referred to in the embodiments provided herein may include at least one of a relational database and a non-relational database. The non-relational database may include, but is not limited to, a blockchain-based distributed database, and the like. The processor referred to in the embodiments provided in the present application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic unit, a data processing logic unit based on quantum computing, or the like, but is not limited thereto.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The foregoing examples illustrate only a few embodiments of the application and are described in detail herein without thereby limiting the scope of the application. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the application, which are all within the scope of the application. Accordingly, the scope of the application should be assessed as that of the appended claims.
Claims (10)
1. A method for locating damage to a bridge structure, the method comprising:
combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain a matrix corresponding to the bridge to be analyzed;
performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix;
performing signal separation processing on target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed;
obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information;
and determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
2. The method of claim 1, wherein after performing principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix, the method further comprises:
determining the target principal component signal according to the principal component contribution rate corresponding to the principal component score matrix;
the step of performing signal separation processing on the target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed, includes:
and carrying out signal separation processing on the target principal component signal through a window filtering model to obtain dynamic vibration information corresponding to the bridge to be analyzed.
3. The method of claim 2, wherein the target principal component signal is a first order principal component signal in the principal component scoring matrix;
the method for obtaining the dynamic vibration information corresponding to the bridge to be analyzed comprises the following steps of:
performing signal separation processing on the first-order main component signal through the window filtering model to obtain a first-order structural vibration mode signal corresponding to the bridge to be analyzed;
and determining dynamic vibration information of the first-order modal frequency corresponding to the bridge to be analyzed according to the difference value between the first-order principal component signal and the first-order structural vibration mode signal.
4. The method of claim 1, wherein performing principal component analysis on the array matrix to obtain a corresponding principal component score matrix comprises:
determining a covariance matrix corresponding to the array matrix;
determining a feature vector matrix corresponding to the covariance matrix;
and obtaining the principal component score matrix according to the eigenvector matrix and the array matrix.
5. The method according to claim 1, wherein determining the structural damage location of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency comprises:
obtaining a damage characteristic quantity index curve according to the time sequence corresponding to the instantaneous frequency;
determining time information of the preset object passing through the structural damage position according to the mutation position in the damage characteristic quantity index curve;
and determining the damage position of the structure according to the time information and the speed information corresponding to the preset object.
6. The method according to claim 1, wherein the obtaining the instantaneous frequency of the bridge to be analyzed from the dynamic vibration information comprises:
and performing Hilbert transform processing on the dynamic vibration information to obtain the instantaneous frequency of the bridge to be analyzed.
7. A bridge construction damage locating device, the device comprising:
the array matrix obtaining module is used for combining the vertical displacement response information of each preset position point of the bridge to be analyzed to obtain an array matrix corresponding to the bridge to be analyzed;
the principal component score matrix obtaining module is used for carrying out principal component analysis processing on the array matrix to obtain a corresponding principal component score matrix;
the dynamic vibration information obtaining module is used for carrying out signal separation processing on the target principal component signals in the principal component score matrix to obtain dynamic vibration information corresponding to the bridge to be analyzed;
the instantaneous frequency obtaining module is used for obtaining the instantaneous frequency of the bridge to be analyzed according to the dynamic vibration information;
and the structural damage position determining module is used for determining the structural damage position of the bridge to be analyzed according to the time sequence corresponding to the instantaneous frequency.
8. The apparatus of claim 7, wherein the principal component score matrix obtaining module is further configured to determine a covariance matrix corresponding to the array matrix; determining a feature vector matrix corresponding to the covariance matrix; and obtaining the principal component score matrix according to the eigenvector matrix and the array matrix.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any of claims 1 to 6 when the computer program is executed.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211258909.8A CN115728383B (en) | 2022-10-14 | 2022-10-14 | Bridge structure damage positioning method, device, computer equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211258909.8A CN115728383B (en) | 2022-10-14 | 2022-10-14 | Bridge structure damage positioning method, device, computer equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115728383A CN115728383A (en) | 2023-03-03 |
| CN115728383B true CN115728383B (en) | 2023-12-05 |
Family
ID=85293553
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211258909.8A Active CN115728383B (en) | 2022-10-14 | 2022-10-14 | Bridge structure damage positioning method, device, computer equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115728383B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107733541A (en) * | 2017-11-29 | 2018-02-23 | 广东工业大学 | Method, apparatus, equipment and the computer-readable recording medium of frequency spectrum perception |
| CN109406076A (en) * | 2018-11-19 | 2019-03-01 | 暨南大学 | A method of beam bridge structure damage reason location is carried out using the mobile principal component of displacement sensor array output |
| CN110619106A (en) * | 2019-07-25 | 2019-12-27 | 中际物联科技（深圳）有限公司 | Bridge damage positioning method and quantitative method thereof |
| CN110657882A (en) * | 2019-09-23 | 2020-01-07 | 暨南大学 | Bridge real-time safety state monitoring method utilizing single-measuring-point response |
-
2022
- 2022-10-14 CN CN202211258909.8A patent/CN115728383B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107733541A (en) * | 2017-11-29 | 2018-02-23 | 广东工业大学 | Method, apparatus, equipment and the computer-readable recording medium of frequency spectrum perception |
| CN109406076A (en) * | 2018-11-19 | 2019-03-01 | 暨南大学 | A method of beam bridge structure damage reason location is carried out using the mobile principal component of displacement sensor array output |
| CN110619106A (en) * | 2019-07-25 | 2019-12-27 | 中际物联科技（深圳）有限公司 | Bridge damage positioning method and quantitative method thereof |
| CN110657882A (en) * | 2019-09-23 | 2020-01-07 | 暨南大学 | Bridge real-time safety state monitoring method utilizing single-measuring-point response |
Non-Patent Citations (1)
| Title |
|---|
| 基于移动主成分法的桥梁损伤识别研究;赵展;林逸洲;钟晓林;张鹏;;广州建筑(第02期);30-33 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115728383A (en) | 2023-03-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2019232667A1 (en) | Sensor arrangement method for reducing uncertainty in structural mode identification | |
| Cumbo et al. | Advanced optimal sensor placement for Kalman-based multiple-input estimation | |
| Dang et al. | Discrete Hopfield neural network based indoor Wi-Fi localization using CSI | |
| CN103018729A (en) | Method for calculating radar scattering cross section of metal cylindrical calibration body | |
| CN115728383B (en) | Bridge structure damage positioning method, device, computer equipment and medium | |
| Li et al. | Nonlinear model identification from multiple data sets using an orthogonal forward search algorithm | |
| CN114692529B (en) | CFD high-dimensional response uncertainty quantification method and device, and computer equipment | |
| CN112182739A (en) | Aircraft structure non-probability credible reliability topological optimization design method | |
| CN115810012B (en) | Transmission tower inclination detection method, device, equipment and storage medium | |
| CN110378041B (en) | Power grid state estimation method, system, equipment and computer storage medium | |
| CN115420418B (en) | Air pressure measuring method and device, electronic equipment and readable storage medium | |
| Zhang et al. | Damage identification through finite element model updating using wavelet damage function | |
| CN111366972B (en) | Goaf power transmission line monitoring method and device | |
| CN114358443B (en) | Average material removal rate prediction method, device, electronic device and storage medium | |
| CN116611369B (en) | Interpolation method and device based on smoothness magnitude and candidate template point number | |
| CN116663343A (en) | Sensor measuring point layout method and device based on sparse reconstruction theory | |
| CN110380413B (en) | PMU (phasor measurement Unit) arranging method, system, equipment and computer medium in power grid | |
| Liu et al. | Measurements Selection for Bias Reduction in Structural Damage Identification | |
| CN114065447A (en) | Power transmission tower monitoring sensor arrangement method and device, terminal equipment and storage medium | |
| CN116702571A (en) | Numerical simulation method and device based on multiple smoothness measurement factors | |
| CN117741646A (en) | Method, device, computer equipment and storage medium for determining speed of target object | |
| CN116796280A (en) | Data processing method, apparatus, device, storage medium and computer program product | |
| CN117473312A (en) | Bearing state prediction method, bearing state prediction device, computer equipment and storage medium | |
| CN116777073A (en) | Feedback correction high-resolution wind energy resource development amount and distribution determination method | |
| CN115309281A (en) | Method and system for determining position and computer readable storage medium |
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 |