CN109781863A - Structure erosion two-stage detection method and its system based on fast vibration test - Google Patents
Structure erosion two-stage detection method and its system based on fast vibration test Download PDFInfo
- Publication number
- CN109781863A CN109781863A CN201910037785.2A CN201910037785A CN109781863A CN 109781863 A CN109781863 A CN 109781863A CN 201910037785 A CN201910037785 A CN 201910037785A CN 109781863 A CN109781863 A CN 109781863A
- Authority
- CN
- China
- Prior art keywords
- strain
- corrosion
- gauge length
- flexibility
- long gauge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a kind of structure erosion two-stage detection methods and its system based on fast vibration test, the strain flexibility of structure is judged according to long gauge length FBG sensor strain-responsive and impact force time-histories, corrosion location is positioned using the strain Flexibility Difference of corrosion front-end geometry, realizes the corrosion orientation problem of first stage；After determining corrosion location, structure unknown parameter to be identified is reduced to the rigidity of erosion unit by the rigidity of all units, establish the sensitivity matrix based on long gauge length strain flexibility to shaped structural section rigidity, according to sensitivity equation, structure erosion element stiffness is solved, etching extent is judged according to stiffness variation；Quick Test Vehicle, data collection system and the data analysis system of mating configuration percussion mechanism form the two-stage policy structure erosion detection system based on fast vibration simultaneously, it is chiefly used in the current safe state assessment and long-term behaviour assessment of structure, measuring accuracy height, Noise Resistance Ability are by force, it can be achieved that the long term monitoring of structure erosion.
Description
Fields
The invention belongs to transport structure health monitoring technique field, it is related to realizing the quick positioning of pavement structure corrosion, determines
Amount and security state evaluation technology particularly relate to a kind of structure erosion two-stage detection side based on fast vibration test
Method and its system.
Background technique
It is well known that the countries in the world infrastructure scale of construction is huge, critically important effect is played in national economy life,
And infrastructure, due to the effect that environmental attack, daily military service load even overload, causes structural within the very long military service period
Degeneration can gradually occur and the invasion of the extreme natural calamities such as earthquake, typhoon, therefore, more and more people may be met at any time
The Performance Evaluation of bridge is begun to focus on, optimizes maintenance management expense to ensure structure safety and realize.
In all kinds of infrastructure service lifes, the corrosion of reinforcing bar has always been considered as being to cause structural damage and deterioration most
One of the main reasons, the reinforced concrete structure being exposed to especially for those in corrosive environment.For many years, people endeavour
In developing reliable, efficient corrosion monitoring instrument, some scholars propose it is a kind of using Hall-effect magnetic to reinforcing bar
The innovative technology of concrete structure progress Corrosion Evaluation；Some scholars propose low frequency (LF) RFID sensor and select transient state special
Sign characterization etch state, these instruments are very sensitive to structure erosion, but if not covering corrosion area, cannot detect reinforcing bar
Corrosion.And optical fiber is small in size, light-weight, can be attached to structure in the case where not changing original structure mechanical behavior
In surface or embedded structure, thus the corrosion monitoring based on strain sensing technology is product emerging in recent years, by Fibre Optical Sensor
Device is wrapped in around the reinforcing bar after corrosion area polishes, and bar diameter caused by measurement corrosion deposits expands to characterize extent of corrosion,
The features such as strain transducer based on fiber grating (FBG) is high with its precision, sensing capability is stable, good reliability is applicable in the most.
However the problem of following challenge, seriously hinders the further development of the corrosion monitoring process based on strain: 1) passing
The limitation of sensor: traditional point strain transducer is not suitable for large scale civil engineering, because these strain gauge gauge lengths are shorter
Otherwise (about 1~2 centimetre) is unable to monitor corrosion, and reinforcement corrosion unless such sensor fixed point is mounted on corrosion area
Residue can have an impact its fixed point.2) limitation of method: usually large-scale, complex form civil engineering structure needs to consider big
The unknown parameter of amount, this will lead to, and equation solution analysis convergence rate is slow, nonuniqueness is poor, and therefore, many is based on vibration-testing
Damnification recognition method there was only the ability for positioning corrosion, but the severity of not energetic corrosion.Therefore, if both of these problems
It is all resolved, then strain measurement is effectively detected in structure erosion will greatly improve with the application in safety evaluation.
Summary of the invention
The present invention provides a kind of based on fast exactly for the prior art in deficiency present on corrosion positioning, quantitative problem
The two-stage policy structure erosion detection method and its system of speed vibration, according to long gauge length FBG sensor strain-responsive and impact
Power time-histories judges the strain flexibility of structure, orients corrosion location using the strain Flexibility Difference of corrosion front-end geometry, realizes the
The corrosion orientation problem in one stage；After determining corrosion location, structure unknown parameter to be identified is contracted by the rigidity of all units
It is kept to the rigidity of erosion unit, the sensitivity matrix based on long gauge length strain flexibility to shaped structural section rigidity is established, according to sensitive
Function and structure Bulk stiffness matrix are spent, structural unit rigidity is solved, etching extent is judged according to stiffness variation, is chiefly used in structure
Current safe state assessment and long-term behaviour assessment, measuring accuracy height, Noise Resistance Ability are by force, it can be achieved that the long-term prison of structure erosion
It surveys.
To achieve the goals above, the technical solution adopted by the present invention is that: the two-stage policy structure based on fast vibration
Corrosion detective method, comprising the following steps:
S1 is responded based on the impact force time-histories under impact vibration and the long gauge length FBG sensor for being fixed on body structure surface, is known
It Chu not structural strain flexibility；
S2, according to strain Flexibility DifferenceDamage index is calculated, is determined to corrode
Position, whereinThe long gauge length strain flexibility of respectively intact structure and the long gauge length strain of corrosion structure are soft
Degree；
S3 establishes original finite element model, obtains Bulk stiffness matrix [K], is positioned according to the corrosion that step S2 is obtained
As a result, establishing long gauge length flexibility matrix to the sensitivity matrix [S] of corrosion location section rigidity；
S4 obtains sensitivity analysis function according to the sensitivity matrix that step S3 is obtained Wherein Δ k is the variable quantity of element stiffness,For the variable quantity for straining flexibility, Nu is parameter number to be identified
Amount, M are strain columns flexible；
S5 brings sensitivity function [S] structure Bulk stiffness matrix [K] into, solves structure through the oriented corruption of step S2
The rigidity for losing unit, judges etching extent according to stiffness variation.
It as an improvement of the present invention, further include step S6, according to non-when geodesic structure is reinforced concrete structure
Linear stage Rigidity Calculation formulaSolve current structure reinforcing bars residual area to be measured, judgement
Etching extent.
As an improvement of the present invention, the step S1 further comprises:
S11 records impact force and long gauge length strain-responsive under impact vibration, judges long gauge length strain FRFs；
S12 identifies that elementary structure parameter, the structural parameters include at least structure frequency by multiple mould indicator function method
ωr, damping ratio ξr, vibration shape zoom factor QrAnd structural strain Mode Shape
S13, using improved conjugation beam method from the strain mode vibration shapeThe displacement model vibration shape is calculated
Strain flexibility is calculated using following formula in S14:
Wherein, symbol * indicates complex conjugate；
It is improved as another kind of the invention, in the step S2, takes the diagonal entry of flexibility matrix difference by normalizing
After change processing, take following indexes as damage criterion:
It improves as another kind of the invention, in the step S12, singular value decomposition is strained applied to long gauge length
FRFs is obtained:
Wherein,WithRespectively left singular vector and right singular vector, by frequency information and resistance
The singular value matrix of Buddhist nun's information composition；ULFor left singular vector, the Mode Shape information comprising all elements left sibling；URFor institute
There is the Mode Shape information of the right node of element；ηm=hm/Lm,hmTo be mounted between the long gauge length sensor of beam bottom portion central axes
Distance；LmFor long gauge length element length.
It is improved as another of the invention, vibration shape zoom factor Q in the step S12rAre as follows:
Wherein,To reinforce
The frequency response function of r rank mode afterwards.
As a further improvement of the present invention, long gauge length strains flexibility for i-th of unit section in the step S3
Rigidity kiSensitivity are as follows:
Wherein,Flexibility matrix m row q column element is strained for long gauge length；The vibration shape is moved for r component level；It is
M-th of element in the long gauge length strain vibration shape of r rank；ωrFor the r rank natural frequency of vibration；
Strain flexibility is extended to one class number of Taylor, can be obtained as follows:
To achieve the goals above, the present invention also the technical solution adopted is that: the two-stage policy knot based on fast vibration
Structure corrosion detective system, Quick Test Vehicle, data collection system and data analysis system including configuring percussion mechanism；
The Quick Test Vehicle of the configuration percussion mechanism, occurs impact vibration to geodesic structure for moving quickly through in vehicle
When, generate impact force；
The data collection system, when for acquiring and store percussion mechanism impact bridge the impact force time-histories that generates and
The Dynamic Response Information of measured structure long gauge length sensor under impact load；
The data analysis system includes structural strain pliability parameter identification module, corrosion locating module and the quantitative mould of corrosion
Block, structural strain degree parameter identification module structure frequency, damping ratio, vibration shape zoom factor, contingency model vibration for identification
Type and displacement model modal parameters, and by structural strain flexibility is calculated；The locating module is according to strain Flexibility Difference meter
Damage index is calculated, to carry out corrosion positioning；It is described to corrode quantitative module for bringing sensitivity analysis function into structure entirety
Stiffness matrix solves structural unit rigidity, to judge etching extent according to stiffness variation.
Compared with prior art, the provided two-stage policy structure erosion detection side based on fast vibration of the invention
Method and its system, the beneficial effect is that:
(1) traditional corrosion detecting method, it is necessary to which corrosion sensor is placed on corrosion by known corrosion location in advance
The single detection etching extent in position, is difficult to find unknown corrosion, and two-stage of the invention corrodes positioning and quantitative method, can be with
The identification of corrosion positioning, quantitative two levels is only realized by a kind of response of sensor, and has the testing time few, test
Precision is high and Noise Resistance Ability is strong feature, it can be achieved that reinforcement corrosion long term monitoring, be widely used in practical bridge
The good prospect that can be assessed.
(2) Modal Flexibility is a kind of damage evaluation index more more sensitive than Mode Shape and intrinsic frequency.Strain flexibility
Structural unit is defined as to the strain-responsive of unit input power, it is directly related with the rigidity of structure, it is more suitable for safety appraisement of structure.
(3) based on impact vibration impact force time-histories and the response time-histories of long gauge length FBG sensor as input, output information
Can identification structure erosion more accurately, the impact vibration of " Bian Yidong, Bian Jili " based on mobile impact vehicle detects survey side
Method can realize the quick positioning of corrosion and quantify.
(4) long gauge length optical fibre sensor realizes the macro microscopic measurement of structure by the mean strain within the scope of monitoring gauge length
(such as 1~2 meter).Compared with traditional foil gauge can only carry out " part " measurement, long gauge length optical fibre grating sensor can be caught
The structure composition information of entire overlay area is grasped, to detect unpredictalbe damage.
(5) one stage damage reason locations lead to being greatly lowered for unknown parameter in sensitivity equation, to be more advantageous to two
Stage corrodes quantitative realization, and influence factor is less, and precision is higher.
(6) two-stage corrodes positioning and quantitative method, is the novel corrosion monitoring process based on shock-testing and Fibre Optical Sensor,
Structural strain flexibility is identified based on FBG Fibre Optical Sensor, and can be gradually according to the sensitivity analysis of strain Flexibility Difference and strain flexibility
Realize that corrosion positions and quantitatively, establishes solid theoretical basis, ensure that the good operation of system.
Detailed description of the invention
Fig. 1 is that the present invention is based on the work flow diagrams of the two-stage policy structure erosion detection system of fast vibration；
Fig. 2 is the girder steel Corrosion monitoring laboratory model construction figure of the embodiment of the present invention 3；
Fig. 3 is 3 structural strain flexibility identification process datagram of the embodiment of the present invention, wherein
Fig. 3 (a) is impact force time course data figure；
Fig. 3 (b) is that long gauge length strain facies answers datagram；
Fig. 3 (c) is strain FRFs datagram；
Fig. 3 (d) is structural strain Mode Shape datagram；
Fig. 3 (e) is displacement modes Data of Mode figure；
Fig. 3 (f) is strain flexibility matrix datagram；
Fig. 3 (g) is strain prediction comparison diagram；
Fig. 4 is the corrosion orientated damage schematic diagram of three kinds of cases in embodiment 3；
Fig. 4 (a) is the corrosion orientated damage schematic diagram of case 1 in embodiment 3；
Fig. 4 (b) is the corrosion orientated damage schematic diagram of case 2 in embodiment 3；
Fig. 4 (c) is the corrosion orientated damage schematic diagram of case 3 in embodiment 3；
Fig. 5 is the corrosion quantitative data figure of three kinds of cases in embodiment 3.
Specific embodiment
Below with reference to drawings and examples, the present invention is described in detail.
Embodiment 1
Two-stage policy structure erosion detection method based on fast vibration, comprising the following steps:
S1 is responded based on the impact force time-histories under impact vibration and the long gauge length FBG sensor for being fixed on body structure surface, is known
It Chu not structural strain flexibility；
S11 records impact force and long gauge length strain-responsive under impact vibration, judges long gauge length strain FRFs；
S12 identifies elementary structure parameter by multiple mould indicator function method, and singular value decomposition (SVD) is applied to length and is answered
Become FRFs, obtain:
WhereinandRespectively left singular vector and right singular vector, be by frequency information and
Damp the singular value matrix of information composition；ULFor left singular vector, the Mode Shape information comprising all elements left sibling；URFor
The Mode Shape information of the right node of all elements；ηm=hm/Lm,hmFor be mounted on beam bottom portion central axes long gauge length sensor it
Between distance, LmFor long gauge length element length, long gauge length strain FRF can be write as following format again:
In formulaFor the strain mode vibration shape, For γrConjugation, ωrFor structure frequency, ξrFor damping ratio, L is that mode participates in matrix,
Wherein Lr=Qr·φr,drvMatrix, φ are participated in for rth moder,drvFor driving point modal vector, QrFor vibration shape zoom factor；
S13, using improved conjugation beam method from the strain mode vibration shapeThe displacement model vibration shape is calculated
As can be seen that (1) formula is identical as (2) formula form, therefore, structure frequency, damping ratio can be identified from (1) formula
And the vibration shape, but vibration shape zoom factor QrIt is also uncertain at present, it can be solved from the least-squares estimation of (3) formula:
In formulaAfter reinforcing
R rank mode frequency response function；
Strain flexibility is calculated using following formula in S14:
In formulaIndicate to be based on improving conjugate beam method from strain modeThe displacement modes vibration shape of obtained rth mode,
Symbol * indicates complex conjugate, it is notable that formula (4) requires no knowledge about the quality of structure, is shaken using the strain mode identified
TypeThe displacement modes vibration shapeVibration shape zoom factor QrAnd it is counted by the parameter that structure frequency and damping information form
It calculates.The vibration shape zoom factor of r rank mode and the relationship of modal mass are Qr=1/2j ωrMr, as modal mass MrWhen equal to 1,
Available mass normalisation Mode Shape:The Mode Shape not scaled in formula (4) is replaced
For mass normalisation Mode Shape and the corresponding strain mode vibration shape, straining flexibility will export in another form:
S2, corrosion would generally cause the variation of the rigidity of structure, these variations will lead to the variation of strain flexibility.Therefore, complete
The difference of strain flexibility and corrosion structure the strain flexibility of whole structure can corrode in a stage identifies corrosion damage position in positioning
It sets, according to following formula, calculates damage index:
WhereinFor the strain flexibility of complete structure and damaged structure；
In order to amplify corrosion front and back flexibility difference to be easier identification of damage, the diagonal entry of flexibility matrix difference is taken to pass through
As damage criterion, which passes through standardization and obtains following index fetching number after crossing normalized:
S3 establishes original finite element model, obtains Bulk stiffness matrix [K], is positioned according to the corrosion that step S2 is obtained
As a result, establishing long gauge length flexibility matrix to the sensitivity matrix [S] of corrosion location section rigidity；
During in the two-stage, corrosion is quantitative, derive based on the corrosion quantitative approach for straining damage sensitivity equation flexible,
By formula (5) it is found that straining flexibility to the derivative of i-th cell rigidity are as follows:
The equilibrium equation of not damaged structural vibration equation are as follows:
By formula (9) to element stiffness kiDerivation obtains formula (10):
Because mass matrix [M] is independently of stiffness matrix ki, thusEquation both sides all multiplied byBased on kiR rank intrinsic frequency sensitivity coefficient it is available:
Since the Mode Shape of structure is mutually indepedent, complete Mode Shape space can be formed, so r rank mould
The sensitivity coefficient of the state vibration shape can be expressed as the combining form of Mode Shape:
α in formulaLFor L rank Mode Shape weight coefficient.Formula (12) are substituted into formula (10), and are multiplied by simultaneously on both sidesSince the vibration shape is orthogonal, soIt brings into
In formula (10):
Solve αs:
Work as s=r, when αs=0, the long gauge length strain mode vibration shapeIt may be expressed as:
Therefore, the long gauge length strain mode vibration shape of r rank is to kiDerivative are as follows:
Formula (11-12), formula (14), formula (16) substitution formula (8) can be obtained to the sensitivity coefficient of long gauge length strain mode flexibility:
Strain flexibility can be extended to one class number of Taylor, as follows:
S4 obtains sensitivity analysis function according to step S3:
[S] is the sensitivity matrix uniquely determined by original configuration in formula；Δ k is the variable quantity of element stiffness；
For the variable quantity for straining flexibility；Nu is number of parameters to be identified；M is strain columns flexible.
S5 brings sensitivity function into structure Bulk stiffness matrix [K], the variable quantity including straining flexibilityEtc. correlations
Parameter, can solve structural unit rigidity, judge etching extent according to stiffness variation.
Embodiment 2
The present embodiment difference from example 1 is that: for when geodesic structure be steel construction when, that is, it is rigid that unit can be used
Degree variation directly solves the structural unit rigidity obtained by step S5, directly judges etching extent as quantitative index is corroded；
It further include step S6, according to nonlinear phase Rigidity Calculation formula when geodesic structure is reinforced concrete structureCurrent structure reinforcing bars area and corroded area to be measured are solved, judges etching extent.
Embodiment 3
Two-stage policy structure erosion detection system based on fast vibration, as shown in Figure 1, the present embodiment is used for system
For bridge experiment, Quick Test Vehicle, data collection system and data analysis system including configuring percussion mechanism；
The Quick Test Vehicle of the configuration percussion mechanism, for producing when vehicle moves quickly through bridge generation impact vibration
Raw impact force, the key position of the bridge are fixed with long gauge length sensor, the sensing theory partition structure based on area distribution
Long gauge length sensor is fixed on the structure floor central axes of key position by key position, and structure is drawn according to long gauge length position
It is divided into several long gauge length units.It should be noted that all the sensors are all answered by the generation of suitable power in advance during the installation process
Become, to obtain the accurate signal in low-level strain.Before long gauge length FBG sensor is installed, diamond dresser is used
Bridge concrete surface is polished flat, the coverings such as cement are removed, is then applied after cleaning concrete surface using acetone
Epoxy primer is added to paste long gauge length sensor comprehensively, to ensure persistently to fix with concrete surface dduring test, with impact
The mobile impact vehicle of device at the uniform velocity passes through bridge to be detected, is impacted in scheduled structure position, rushed in such a manner
Hitting vehicle, " Bian Yidong, Bian Jili " complete the shock-testing of whole bridge beam.
The data collection system carries out real-time data acquisition record, for acquiring and storing percussion mechanism impact bridge
When the Dynamic Response Information of impact force time-histories and the measured structure long gauge length sensor under impact load that generates, the present embodiment
The record of impact force time-histories is used for using NI PXIe-1082 data collection system, MOI (FBG) demodulator is demodulated for lightwave signal, and
Carry out long gauge length strain acquirement record in real time.
The data analysis system includes structural strain pliability parameter identification module, corrosion locating module and the quantitative mould of corrosion
The identification of (1) Structure Flexibility Matrix: block using the impact force of record and long gauge length strain-responsive estimation strain FRF, then passes through
CMIF method identifies basic friction angle parameter: including formation zoom factor, long gauge length strain mode etc., and utilizing improved conjugation
The displacement model vibration shape is calculated from strain mode vibration shape in beam method；(2) one stages corrosion positioning: obtain strain flexibility it
Afterwards, equation can be passed throughIt calculates damage index and carries out corrosion positioning；(3) two-stage corrodes
It is quantitative: after corrosion sites determine, quantitative objective is concentrated on into erosion unit, according to strain sensitivity function flexible:It carries out sensitivity analysis and solves structural section element stiffness, realize rotten
The determination of erosion amount.
As shown in Fig. 2, the present embodiment has studied the validity of this method using Simply-supported Steel Beams model, it is 5.76 meters of beam length, total
It is divided into 12 units, using Q235 steel, elasticity modulus 206Gpa, Unit Weight 76.97kN/m3.Assuming that at the bottom of beam
Portion is installed by the long gauge length optical fibre sensor that 12 length are 0.48m, the long gauge length strain-responsive of measurement structure, using FInite Element
Package assembly stiffness matrix and mass matrix, using Reyleigh' damping matrix.In addition to basic structure, three kinds of damage modes are considered: 1)
The single degree of impairment of 1: the 5 element stiffness of case loss 5%；2) the single damage of 2: the 5 element stiffness of case loss 10%
Situation；3) 3: the 5 element stiffness of case loss 10%, the multiple damages situation of the 9th element stiffness loss 15%.
Structure basic model and damage model are simulated in SAP2000 software, it is wide by reducing steel beam flange
Degree is simulated corrosion damage.Apply impact force respectively as excitation in the 5th node and the 8th node, while simulating
5% white noise is added in data indicates observation noise.
It strains flexibility identification: by taking complete structure as an example, illustrating the process for straining flexibility identification.Typical impact force and length
Gauge length strain-responsive as shown in Fig. 3 (a) and (b), is estimated to grow according to the long gauge length strain-responsive during impact test respectively
Gauge length strains FRFs, after determining long gauge length strain FRF matrix, carries out singular value decomposition by following formula,
Shown in the drawing such as Fig. 3 (c) of singular value as frequency function, curve number corresponds to impact test input in Fig. 3 (c)
Number, wherein structural modal is indicated with the peak value for the spectral line that circle is mark.As figure shows, the peak value singular value on CMIF figure can
It can be the pole location of structure, the frequency of first three mode is respectively 15.463Hz, 61.197Hz, 135.37Hz.Equally, it determines
Damping ratio is respectively 0.51%, 1.09% and 0.56%.FromIn can also identify structural strain Mode Shape such as Fig. 3 (d)
It is shown；The displacement modes vibration shape is calculated from the strain mode vibration shape using improved conjugate beam method, as shown in Fig. 3 (e).
By following formula, vibration shape zoom factor is calculated:
In formulaAfter reinforcing
R rank mode frequency response function, once structure frequency, damping ratio, vibration shape zoom factor, the contingency model vibration shape and displacement model
Recognition of vibration comes out, and can calculate structural strain flexibility, it may be assumed that
In formulaIndicate to be based on improving conjugate beam method from strain modeThe displacement modes vibration shape of obtained rth mode,
Symbol * indicates complex conjugate, shown in obtained strain flexibility matrix such as Fig. 3 (f), because the Simply-supported Steel Beams have 12 units, 13
A node, 12 to the 13 of dimension.
The accuracy for straining flexibility is estimated to verify, the static(al) strain that corresponding slow test measures is drawn and compares figure,
Slow test applies the power of 100N size respectively on node 3,5,8,10 respectively, as shown in Fig. 3 (g), at each occurrence,
Prediction strain and actual measurement strain have preferable consistency, and similar processing is carried out to corrosion structure, corresponding structure is obtained and answers
Become flexibility, as shown in figure 4,.
The corrosion positioning of one stage: difference matrix is calculated according to the structural strain flexibility matrix of corrosion front and back, and is extracted pair
Damage reason location of the linea angulata unit for corrosion, as a result as shown in Fig. 4 (a-c), Fig. 4 is the corrosion orientated damage signal of three kinds of cases
Figure.It can be seen from the figure that the position of structural damage can clearly be navigated to, therefore, strain Flexibility Difference is a kind of fine
Corrosion positioning index, be suitable for it is single damage and multiple damages positioning.
Two-stage corrosion is quantitative: once it is determined that the position of corrosion damage, corrosion is quantitative can only concentrate on corrosion area,
To greatly reduce the quantity of unknown parameter.Then structural unit rigidity is calculated according to sensitivity analysis, wherein with case
5 rigidity of unit of example 1,2,5 rigidity of unit of case 3, the target component that 9 rigidity of situation is lesion quantification, according to the vibration of identification
Type zoom factor determines that the proportionality coefficient of first three rank mode is respectively 0.0003,0.0012,0.0022, it is hereby achieved that quality
The vibration shape is normalized, the sensitivity coefficient of the corresponding frequency of unit under test rigidity, displacement modes and long gauge length strain mode is calculated,
Then Stiffness degradation is quantified by following formula,
Experimental result is as shown in table 1 below, and is compared with theoretical value:
1 corrosion damage quantitative result of table
Using the sensitivity function based on strain flexibility, corrosion damage can be effectively quantified, error is acceptable
In range.Fig. 5 is that the corrosion of three kinds of cases of the present embodiment is quantitatively schemed, and is more intuitively described as a result, discre value and theoretical value are kissed
It closes preferable.
The basic principles, main features and advantages of the present invention have been shown and described above.The technology of the industry
Personnel only illustrate the present invention it should be appreciated that the present invention is not limited by examples detailed above described in examples detailed above and specification
Principle, various changes and improvements may be made to the invention without departing from the spirit and scope of the present invention, these variation and
Improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention is by appended claims and its is equal
Object defines.
Claims (8)
1. the structure erosion two-stage detection method based on fast vibration test, it is characterised in that the following steps are included:
S1 is responded based on the impact force time-histories under impact vibration and the long gauge length FBG sensor for being fixed on body structure surface, is identified
Structural strain flexibility；
S2, according to strain Flexibility DifferenceDamage index is calculated, so that corrosion positioning is carried out, whereinThe long gauge length strain flexibility of respectively intact structure and the long gauge length of corrosion structure strain flexibility；
S3 establishes original finite element model, obtains Bulk stiffness matrix [K], the corrosion positioning result obtained according to step S2,
Long gauge length flexibility matrix is established to the sensitivity matrix [S] of corrosion location section rigidity；
S4 obtains sensitivity analysis equation according to the sensitivity matrix that step S3 is obtainedWherein Δ k is the variable quantity of element stiffness,For strain flexibility
Variable quantity, Nu are number of parameters to be identified, and M is strain columns flexible；
S5 brings sensitivity matrix [S] in sensitivity analysis equation into, solves structure through the oriented erosion unit of step S2
Rigidity judges etching extent according to stiffness variation.
2. the structure erosion two-stage detection method according to claim 1 based on fast vibration test, it is characterised in that:
It further include step S6, according to nonlinear phase Rigidity Calculation formula when geodesic structure is reinforced concrete structureCurrent structure reinforcing bars residual area to be measured is solved, judges etching extent.
3. the structure erosion two-stage detection method according to claim 1 or 2 based on fast vibration test, feature exist
In: the step S1 further comprises:
S11 records impact force and long gauge length strain-responsive under impact vibration, judges long gauge length strain FRFs；
S12 identifies that elementary structure parameter, the structural parameters include at least structure frequency ω by multiple mould indicator function methodr、
Damping ratio ξr, vibration shape zoom factor QrAnd structural strain Mode Shape
S13, using improved conjugation beam method from the strain mode vibration shapeThe displacement model vibration shape is calculated
Strain flexibility is calculated using following formula in S14:
Wherein, symbol * indicates complex conjugate；
4. the structure erosion two-stage detection method according to claim 3 based on fast vibration test, it is characterised in that:
In the step S2, takes the diagonal entry of flexibility matrix difference after normalized, takes following indexes as damage criterion:
5. the structure erosion two-stage detection method according to claim 3 based on fast vibration test, it is characterised in that:
In the step S12, singular value decomposition is applied to long gauge length and strains FRFs, is obtained:
Wherein,WithRespectively left singular vector and right singular vector are believed by frequency information and damping
Cease the singular value matrix of composition；ULFor left singular vector, the Mode Shape information comprising all elements left sibling；URFor all members
The Mode Shape information of the right node of element；ηm=hm/Lm,hmFor be mounted between the long gauge length sensor of beam bottom portion central axes away from
From；LmFor long gauge length element length.
6. the structure erosion two-stage detection method according to claim 4 or 5 based on fast vibration test, feature exist
In: vibration shape zoom factor Q in the step S12rAre as follows:
Wherein,It is strengthened
The frequency response function of r rank mode.
7. the structure erosion two-stage detection method according to claim 3 based on fast vibration test, it is characterised in that:
Long gauge length strain flexibility is for i-th of unit section rigidity k in the step S3iSensitivity are as follows:
Wherein,Flexibility matrix m row q column element is strained for long gauge length；The vibration shape is moved for r component level；For r rank
M-th of element in the long gauge length strain vibration shape；ωrFor the r rank natural frequency of vibration；
Strain flexibility is extended to one class number of Taylor, can be obtained as follows:
8. the structure erosion two-stage detection system based on fast vibration test, it is characterised in that: including configuration percussion mechanism
Quick Test Vehicle, data collection system and data analysis system；
The Quick Test Vehicle of the configuration percussion mechanism is produced for moving quickly through when impact vibration occurs for geodesic structure in vehicle
Raw impact force；
The data collection system, for acquiring and storing the impact force time-histories generated when percussion mechanism impacts bridge and be tested
The Dynamic Response Information of structure long gauge length sensor under impact load；
The data analysis system includes structural strain pliability parameter identification module, corrosion locating module and corrosion quantitatively module,
The structural strain degree parameter identification module for identification structure frequency, damping ratio, vibration shape zoom factor, the contingency model vibration shape and
Displacement model modal parameters, and by structural strain flexibility is calculated；The locating module calculates damage according to strain Flexibility Difference
Hurt index, to carry out corrosion positioning；It is described to corrode quantitative module for bringing sensitivity analysis function into structure overall stiffness
Matrix solves structural unit rigidity, to judge etching extent according to stiffness variation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910037785.2A CN109781863B (en) | 2019-01-16 | 2019-01-16 | Structure corrosion two-stage detection method and system based on rapid vibration test |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910037785.2A CN109781863B (en) | 2019-01-16 | 2019-01-16 | Structure corrosion two-stage detection method and system based on rapid vibration test |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109781863A true CN109781863A (en) | 2019-05-21 |
| CN109781863B CN109781863B (en) | 2021-07-23 |
Family
ID=66500664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910037785.2A Active CN109781863B (en) | 2019-01-16 | 2019-01-16 | Structure corrosion two-stage detection method and system based on rapid vibration test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109781863B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110516297A (en) * | 2019-07-19 | 2019-11-29 | 东南大学 | Multi-level finite element modeling model modification method based on long gauge length time-frequency domain information |
| CN111579647A (en) * | 2020-07-06 | 2020-08-25 | 中南大学 | Concrete member corrosion degree detection method and system based on analytic hierarchy process |
| CN113639941A (en) * | 2021-06-30 | 2021-11-12 | 山东建筑大学 | Method for testing bridge strain flexibility matrix |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04235346A (en) * | 1991-01-11 | 1992-08-24 | Mitsubishi Heavy Ind Ltd | Method for detecting grain-boundary segregation property of metal material |
| CN104198144A (en) * | 2014-09-12 | 2014-12-10 | 东南大学 | Middle and small bridge fast detecting method based on long-scale-distance optical fiber strain sensor |
| CN104931373A (en) * | 2015-06-12 | 2015-09-23 | 中国科学院金属研究所 | Corrosion fatigue crack propagation testing device |
| CN105953996B (en) * | 2016-06-30 | 2017-04-19 | 东南大学 | Bridge detection and evaluation method and equipment based on impact vibration |
| CN107449676A (en) * | 2017-08-17 | 2017-12-08 | 中广核研究院有限公司 | Crackle aperture displacement measurement apparatus |
-
2019
- 2019-01-16 CN CN201910037785.2A patent/CN109781863B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH04235346A (en) * | 1991-01-11 | 1992-08-24 | Mitsubishi Heavy Ind Ltd | Method for detecting grain-boundary segregation property of metal material |
| CN104198144A (en) * | 2014-09-12 | 2014-12-10 | 东南大学 | Middle and small bridge fast detecting method based on long-scale-distance optical fiber strain sensor |
| CN104931373A (en) * | 2015-06-12 | 2015-09-23 | 中国科学院金属研究所 | Corrosion fatigue crack propagation testing device |
| CN105953996B (en) * | 2016-06-30 | 2017-04-19 | 东南大学 | Bridge detection and evaluation method and equipment based on impact vibration |
| CN107449676A (en) * | 2017-08-17 | 2017-12-08 | 中广核研究院有限公司 | Crackle aperture displacement measurement apparatus |
Non-Patent Citations (2)
| Title |
|---|
| JIAN ZHANG,ET AL: "Strain flexibility identification of bridges from long-gauge strain measurements", 《MECHANICAL SYSTEMS AND SIGNAL PROCESSING》 * |
| 郭双林: "基于冲击振动的位移柔度和应变柔度识别方法", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110516297A (en) * | 2019-07-19 | 2019-11-29 | 东南大学 | Multi-level finite element modeling model modification method based on long gauge length time-frequency domain information |
| CN111579647A (en) * | 2020-07-06 | 2020-08-25 | 中南大学 | Concrete member corrosion degree detection method and system based on analytic hierarchy process |
| CN113639941A (en) * | 2021-06-30 | 2021-11-12 | 山东建筑大学 | Method for testing bridge strain flexibility matrix |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109781863B (en) | 2021-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104198144B (en) | Middle and small bridge fast detecting method based on long-scale-distance optical fiber strain sensor | |
| CN107389285B (en) | A kind of quick test and evaluation method of bridge changed based on temperature | |
| WO2021159569A1 (en) | Method for engineering structural fiexibility identification considering the impact of noise | |
| CN102034021A (en) | Integral and local information fusing method of structure health diagnosis | |
| CN104517036A (en) | Simply-supported piece damage identification method based on strain statistical moment | |
| CN109781863A (en) | Structure erosion two-stage detection method and its system based on fast vibration test | |
| Fayyadh et al. | Combined modal parameters-based index for damage identification in a beamlike structure: theoretical development and verification | |
| Ni et al. | Dynamic property evaluation of a long-span cable-stayed bridge (Sutong bridge) by a Bayesian method | |
| CN108280294B (en) | A kind of cable arch structure damage combined recognising method based on modal parameter | |
| Hong et al. | Condition assessment of reinforced concrete beams using dynamic data measured with distributed long-gage macro-strain sensors | |
| CN109299539A (en) | A method of utilizing curvilinear corner index evaluation cantilever beam rigidity weak link | |
| Wu et al. | Two-level damage detection strategy based on modal parameters from distributed dynamic macro-strain measurements | |
| CN101539493A (en) | Symmetric signal method for structural damage diagnosis | |
| CN110926735A (en) | Bridge structure rapid diagnosis method based on multidimensional dynamic parameters | |
| CN110470738A (en) | Structural Damage Identification based on vibratory response difference than function | |
| CN102252792B (en) | Absolute axial force test method of bar | |
| CN107664489A (en) | A kind of bridge strain and the measuring method of deformation | |
| CN106383003A (en) | Cable structure cable force measurement method and system based on flexibility identification | |
| Jiang et al. | Structural stiffness identification of traditional mortise-tenon joints based on statistical process control chart | |
| Xiang et al. | Vibration-based structural damage identification using P-CNN and time-frequency hybrid index under the conditions of uncertainties and incomplete measurements | |
| CN113639941A (en) | Method for testing bridge strain flexibility matrix | |
| Guan | Vibration-based structural health monitoring of highway bridges | |
| KR102435166B1 (en) | A method of determining the measurement location and location of a structure using a genetic algorithm and AI technology for precise strain-displacement prediction | |
| CN115713020A (en) | Rapid test and evaluation method for bearing rigidity of simply supported beam bridge based on local vibration mode | |
| Riyanti et al. | An Automatic Load Detector Design to Determine the Strength of Pedestrian Bridges Using Load Cell Sensor Based on Arduino |
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 |