CN107228772A - A kind of shield tunnel construction method for estimating damage - Google Patents
A kind of shield tunnel construction method for estimating damage Download PDFInfo
- Publication number
- CN107228772A CN107228772A CN201710434506.7A CN201710434506A CN107228772A CN 107228772 A CN107228772 A CN 107228772A CN 201710434506 A CN201710434506 A CN 201710434506A CN 107228772 A CN107228772 A CN 107228772A
- Authority
- CN
- China
- Prior art keywords
- damage
- mrow
- shield tunnel
- frequency band
- tunnel construction
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/004—Testing the effects of speed or acceleration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/007—Subject matter not provided for in other groups of this subclass by applying a load, e.g. for resistance or wear testing
Abstract
The invention discloses a kind of shield tunnel construction method for estimating damage, it is characterized in that, including, not damaged operating mode is set up according to the parameter of shield tunnel respectively and has the shield tunnel construction model of damage regime, apply traveling load on the mold, obtain the not damaged of measuring point and have damage acceleration responsive signal;To not damaged acceleration responsive and there is damage acceleration responsive signal to carry out WAVELET PACKET DECOMPOSITION respectively, obtain the wavelet-packet energy of each frequency band;According to the wavelet-packet energy of each frequency band, by the rate of change deviation of energy ratio, the damage position and/or degree of injury of shield tunnel construction are determined.The method of technical solution of the present invention, using the damnification recognition method based on data, it is not necessary to the vibration shape information of structure, non-destructive tests assessment can be carried out exactly using the acceleration signal before and after structural damage.
Description
Technical field
The invention belongs to civil engineering works structure health monitoring and damage identification technique field, and in particular to a kind of shield tunnel
Structural Damage Assessment method.
Background technology
Tunnel structure makes the daily life to everybody once going wrong as the pith in urban track traffic
Into tremendous influence and loss.In practice, the Service Environment of subterranean tunnel is complicated, with the growth of active time, tunnel structure
To inevitably occur local damage, and if potential safety hazard must be left with safeguarding without detection, cause severe accident to occur.
In fact, on the one hand during operation, the reciprocating motion influence of the dynamic load such as Under Dynamic Load of vehicle, it is also possible to produce structure
Raw different degrees of fatigue damage.On the other hand, building structure generally existing is inevitably constructed initial imperfection, although these
Defect is very small, but accumulation and the continuous effect of internal and external environment with the service time, and initial imperfection can gradually develop expansion,
The health status of tunnel structure can progressively deteriorate with the growth of service time.Therefore, periodically damaged for tunnel structure
It is very necessary that triage, which is surveyed to prevent accident,.
In recent years, monitoring structural health conditions (Structural Health Monitoring, abbreviation SHM) and diagnostic techniques exist
Considerable progress is achieved in the application of civil engineering structure (including skyscraper, large bridge and dam engineering etc.), still
This technology is applied in Tunnel Engineering almost just start starting.With both at home and abroad to the safe operation of urban track traffic
Increasingly pay attention to, the health monitoring systems of urban track traffic underground structure (predominantly shield tunnel construction) increasingly turn into research
Focus, such as:The application study of technology of wireless sensing network and MEMS intelligent sensor nodes technology in subway health monitoring systems
Deng obtaining extensive concern.And the research in terms of shield tunnel non-destructive tests, only it is based on numerical simulation on a small quantity at present
Model has been carried out.For in theory, effective shield tunnel construction damnification recognition method is carried out using magnanimity Monitoring Data
With important theory and engineering application value.
In the prior art, CN201310502959 discloses a kind of Damage Assessment Method based on vibration transmissibility function
Method, it refers to the WAVELET PACKET DECOMPOSITION energy variation amount that vibration transmissibility function before and after structural damage is replicated as damage characteristic
Mark, Damage Assessment Method is carried out using the pattern recognition function of SVMs.The vibration shape information of this method height dependency structure,
But the actually acquisition relative difficulty of vibration shape information;Although and this method is solving small sample, the knowledge of non-linear and high dimensional pattern
There is certain advantage in other problem, but when accurate damage reason location and degree of injury confirmation is carried out, with office
It is sex-limited.In addition,《Structural damage alarming technique study based on Wavelet Packet Energy Spectrum》In one text, although inquired into profit in detail
Energy ratio and the energy ratio change obtained with WAVELET PACKET DECOMPOSITION carries out the processes of non-destructive tests, but this method still have it is larger
Error, further it is accurate carry out non-destructive tests when error it is larger.
The content of the invention
For the disadvantages described above or Improvement requirement of prior art, the invention provides a kind of shield tunnel construction lesion assessment
Method.The method of technical solution of the present invention, using the damnification recognition method based on data, it is not necessary to the vibration shape information of structure, only
Damage reason location can be carried out exactly using the acceleration signal before and after structural damage and degree of injury is assessed.In addition, of the invention
The method of technical scheme, can also accurately be recognized to damage.
To achieve the above object, according to one aspect of the present invention, there is provided a kind of shield tunnel construction lesion assessment side
Method, it is characterised in that including,
S1 sets up the shield tunnel construction model of not damaged operating mode according to the parameter of shield tunnel, applies lotus on the mold
Carry, obtain the not damaged acceleration responsive signal of corresponding measuring point;
S2 sets up the shield tunnel construction model with damage regime according to the parameter of shield tunnel, applies on the mold
The load, obtain corresponding measuring point has damage acceleration responsive signal;
S3 is to the not damaged acceleration responsive and has damage acceleration responsive signal to carry out WAVELET PACKET DECOMPOSITION respectively, obtains
The wavelet-packet energy of each frequency band;
S4 calculates the corresponding not damaged acceleration responsive signal of the measuring point and has the small echo of damage acceleration responsive signal
Bag energy feature frequency band rate of change, is appointed as signatures for damage detection, according to the change of the signatures for damage detection by the rate of change
Change, determine the damage position and/or degree of injury of shield tunnel construction.
In technical solution of the present invention, its key is to set up not damaged operating mode and has the structural model conduct of damage regime
Control, by comparing the acceleration responsive signal under not damaged operating mode and having the acceleration responsive signal in the case of damage regime,
To confirm position and the degree damaged in structure.Specifically, one sensor is set every a segment distance in shield tunnel,
Test obtains the acceleration signal under different operating modes respectively.The signal measured for each sensor, by under not damaged operating mode
The acceleration responsive signal measured is compared the acceleration responsive signal measured in the case of damage regime with benchmark as benchmark,
To determine damage position and degree of injury.In actual conditions, damage regime (e.g., damage after wall by soil erosion, overloading load, rigidity
Wound etc.) it is typically what is be readily apparent that, but this damage can cause the influence of which kind of degree to be to be difficult to predict to shield tunnel
's;The damage of damage regime and the structure of shield tunnel also not necessarily synchronous appearance.
The damage of shield tunnel passes through fine wavelet-packet energy not by acceleration signal Direct Recognition
After decomposition, obtain the wavelet-packet energy of each frequency band, so using wavelet pack energy feature frequency band rate of change from multiple angles come
Identification of damage situation.The energy size that sensor is measured can not be represented directly to be used in degree of impairment, technical solution of the present invention
There is the energy ratio of acceleration responsive in the case of the acceleration responsive measured under degree of impairment and not damaged to assess the position of damage,
The rate of change of the energy ratio of acceleration responsive is commented in the case of the acceleration responsive and not damaged that are measured using having under degree of impairment
The size of assessment of loss wound.This appraisal procedure can effectively reduce error, improve the degree of accuracy of assessment result.
As the optimal technical scheme of the present invention, the wavelet-packet energy of each frequency band is defined as described in step S3
In formula, i is the WAVELET PACKET DECOMPOSITION number of plies, SijThe subsignal of j-th of frequency band, x in the frequency band obtained for i-th layer of decompositionjk
For the value of k-th of discrete point of j-th of frequency band subsignal, N is total points of j-th of frequency band subsignal.
As the optimal technical scheme of the present invention, signatures for damage detection is defined as
In formula, i is the WAVELET PACKET DECOMPOSITION number of plies, ERVjChange for the energy ratio of j-th of frequency band,It is ERVjAverage value.
It is pointed out that including some footmarks in each character of above formula, footmark i is the number of plies of WAVELET PACKET DECOMPOSITION, its generation
The certain sequence number of table;Similarly, footmark j is band number, and footmark k is the sequence number value of discrete point.Footmark i, j, k are respectively provided with necessarily
Span, when above-mentioned character footmark have determine numerical value when, what the amount representated by character was also to determine.
As the optimal technical scheme of the present invention, distinguishing indexes can determine the position of structural damage, and relatively
Determine the degree of structural damage.
Under theoretical case, because the numerical value (ERVD values) in the energy ratio deviation for occurring damage measuring point is remote substantially abnormal, because
This, can verify that the damage criterion can be recognized clearly and orientated damage, and can protrude injury region, reduce the possibility of erroneous judgement.
For the damage index figure under different-stiffness degree of impairment, its shape is similar, but amplitude is different, with the increasing of loss of rigidity,
The numerical value of damage index also accordingly increases therewith.In technical solution of the present invention, the energy ratio deviation can for identification of damage position
To reach very high precision, it can only relatively be determined for the identification of degree of injury according to the amplitude of energy ratio deviation.
As the optimal technical scheme of the present invention, acceleration responsive signal preferably passes through ABAQUS finite element software meters
Obtain.
It is used as the optimal technical scheme of the present invention, the preferred 100Hz of acceleration responsive signal sampling frequencies.
As the optimal technical scheme of the present invention, wavelet packet functions Decomposition order is preferably 8 layers.
In ABAQUS finite element software environment, it can effectively simulate the not damaged of shield tunnel and have damaged structure mould
Type.The preferred 100Hz of acceleration signal sample frequency, preferably 8 layers of wavelet packet functions Decomposition order.
In general, by the contemplated above technical scheme of the present invention compared with prior art, with following beneficial effect
Really:
1) method of technical solution of the present invention, it is not necessary to the vibration shape information of structure, only by the small of structural dynamic response signal
Ripple bag energy spectrum rate of change is calculated, and structural dynamic response is easier to obtain than vibration shape information in practice, and accuracy rate is also more
It is high;
2) method of technical solution of the present invention, goes for single injury, a variety of damage regimes, by WAVELET PACKET DECOMPOSITION,
Noise can be effectively reduced, the damage position of structure is accurately determined and assesses its degree of injury, and to the knot of shield tunnel
Structure carries out damage alarming, is of very high actual application value;
3) method of technical solution of the present invention, using wavelet-packet energy decomposition method, compared to using energy ratio or energy
Amount is turned to signatures for damage detection than change, and energy ratio deviation has higher measurement accuracy, more enough that small error is amplified into convex
Reveal and, contribute to the accurate identification and early warning of damage.
Brief description of the drawings
Fig. 1 is the damage criterion construction flow chart of the embodiment of the present invention;
Fig. 2 is the structural finite element model schematic diagram of the shield tunnel of the embodiment of the present invention;
Fig. 3 is the partial enlarged drawing of the structural finite element model schematic diagram of the shield tunnel of the embodiment of the present invention;
Fig. 4 is damage reason location design sketch of the shield tunnel construction of the embodiment of the present invention under conditions of operating mode 2;
Fig. 5 is damage reason location design sketch of the shield tunnel construction of the embodiment of the present invention under conditions of operating mode 3;
Fig. 6 is damage reason location design sketch of the shield tunnel construction of the embodiment of the present invention under conditions of operating mode 4;
Fig. 7 is damage reason location design sketch of the shield tunnel construction of the embodiment of the present invention under conditions of operating mode 5.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.As long as in addition, technical characteristic involved in each embodiment of invention described below
Not constituting conflict each other can just be mutually combined.With reference to embodiment, the present invention is described in more detail.
As shown in figure 1, be a kind of damage criterion construction flow chart of embodiment of technical solution of the present invention, including following master
Want step.
The first step:The FEM model for not damaging shield tunnel construction is set up, applies running train load, is accordingly surveyed
The acceleration responsive of point.For shield tunnel construction, damage is the change in a kind of structure, therefore it should have a ginseng
The object examined, for a kind of normal, health shield tunnel construction, the change in this structure is exactly to damage.Cause
In the embodiment of this technical solution of the present invention, there is provided not damaging the FEM model of shield tunnel as benchmark, damage needs
By being compared with the benchmark, so that it is determined that the situation of damage.In practice, in shield tunnel have load, such as subway, because
This is set up after the FEM model of undamaged shield tunnel construction, it is necessary to simulate imposed load, the technology of the present invention side thereon
In the embodiment of case, preferably apply running train load.More fully understand in order to which the structure to shield tunnel has, one in engineering
As can a series of sensor is disposed in tunnel according to certain rule, by Sensor monitoring load through changing the time when letter
Number.In the embodiment of technical solution of the present invention, measuring point preferably is set to according to actual sensing station, so as to obtain measuring point
Acceleration responsive signal.
Second step:Damage regime shield tunnel construction model is set up, phase is obtained under same running train load action
Answer the acceleration responsive of measuring point.Compared in the first step as benchmark not damaged shield tunnel construction FEM model, also
Need to set up a shield tunnel construction model with damage regime.Damage regime includes but is not limited to after loss of rigidity, wall
Soil erosion etc., damage regime can be the combination of single damage regime or multiple damage regimes.Accordingly, should
Should also apply on damage model with identical load in the first step, point position is also consistent.Pass through this control variable
Method, in the case where other conditions are consistent, only change damage regime, to ensure that the result measured only sustains damage work
The influence of condition.
3rd step:WAVELET PACKET DECOMPOSITION is carried out respectively to obtained not damaged and the acceleration signal for having damage, obtains each
The wavelet-packet energy of frequency band.From original acceleration signal, preliminary information can only obtain, it is difficult to accurately identify shield tunnel
Damage position and degree of injury.In the embodiment of technical solution of the present invention, it is preferred to use the mode of WAVELET PACKET DECOMPOSITION, to the first step
Detect that the acceleration signal obtained has carried out further analyzing and processing with second step.This method analysis divides time-frequency plane
More careful, the resolution ratio dyadic wavelet of its HFS to signal will height.Moreover, it is on the basis of Wavelet Analysis Theory
On, introduce the concept of optimal base selection.That is, frequency band is passed through after multi-level division, according to the spy of analyzed signal
Levy, adaptively choose optimal basic function, be allowed to match with signal, to improve the analysis ability of signal, drastically increase
The accuracy of signal analysis result.Specifically, vibration signal S is carried out after i layers of WAVELET PACKET DECOMPOSITION, i-th layer has j=2iIndividual one
The frequency band (i.e. frequency band) of fixed width degree, the WAVELET PACKET DECOMPOSITION coefficient of each frequency band of extraction, which is reconstructed, can obtain the son letter of each frequency band
Number Sij.It is possible to further obtain the energy of single frequency band using the subsignal.
4th step:The wavelet pack energy feature frequency band rate of change of each subinterval acceleration signal before and after being damaged with shield tunnel
It is used as signatures for damage detection;According to the change of signatures for damage detection amplitude, the damage position and damage of shield tunnel construction are determined
Hinder degree.Specifically, according to the wavelet pack energy feature frequency band obtained in the 3rd step, first according to a series of energy feature
Frequency band confirms the ratio between the average energy value of whole frequency bands of a WAVELET PACKET DECOMPOSITION layer, the energy and the average energy value of each frequency band,
The energy ratio of each energy feature frequency band is defined as, the position to determine structural damage., can be further according to the energy ratio
The energy ratio deviation for calculating each feature band, deviation is bigger, and degree of injury is bigger, so as to estimate damage journey
The size of degree.
Why Wavelet Packet Energy Spectrum is used to characterize structural dynamic characteristic, and its theoretical foundation is as follows:It is subject to when structure
The timings of F mono- are encouraged, the change of the structural dynamic characteristic as caused by structural damage can cause structural dynamic response S in different analysis rulers
Spend subsignal S on (frequency band)ijChange, i.e. redistribution of the structural dynamic response S energy in each frequency band.Relevant small echo
Bag Energy Decomposition is used for the detailed theoretical foundation for confirming damage, may refer to《Structural damage alarming based on Wavelet Packet Energy Spectrum
Technique study》One text, wherein the computational methods of wavelet-packet energy decomposition are explained in detail.
In one preferred embodiment of technical solution of the present invention, signal energy is defined as follows:
Accordingly, the ENERGY E of j-th of frequency bandjCan be by subsignal SijObtain:
In formula, i is the WAVELET PACKET DECOMPOSITION number of plies, SijThe subsignal of j-th of frequency band, x in the frequency band obtained for i-th layer of decompositionjk
For the value of k-th of discrete point of j-th of frequency band subsignal, N is total points of j-th of frequency band subsignal.Technical solution of the present invention
Embodiment in, as needed, we can carry out the WAVELET PACKET DECOMPOSITION of certain number of plies, i-th after decomposition to acceleration signal
Layer has 2iThe frequency band (frequency band) of individual one fixed width, each frequency band can be obtained by extracting the WAVELET PACKET DECOMPOSITION coefficient reconstruct of each frequency band
Subsignal Sij.The signal of each frequency band can obtain several subsignal points, x by decomposingjkAs j-th frequency band subsignal
K-th of discrete point value, N is total points of j-th frequency band subsignal.
In above-mentioned 4th step, in order to characterize the On Damage State in different stress stages, the reality of technical solution of the present invention
Apply example and define structural damage alarming index on the basis of each frequency band energy is than change.Define energy ratio deviation ERVD works
Know index for damage:
Wherein ERVjIt is the energy ratio change of jth frequency band,It is ERVjAverage value.Under theoretical case, due in hair
The numerical value (ERVD values) of the energy ratio deviation of raw damage measuring point is much larger than other measuring points, therefore, it can verify the damage criterion energy
Enough clearly identification and orientated damage, and injury region can be protruded, reduce the possibility of erroneous judgement.For under different-stiffness degree of impairment
Damage index figure, its shape is similar, but amplitude is different, with the increasing of loss of rigidity, the numerical value of damage index also it is corresponding therewith
Increase.
In the embodiment of technical solution of the present invention, energy ratio change is calculated via following process to be obtained, and energy is defined first
Than for:
According to theory analysis, with the increase of WAVELET PACKET DECOMPOSITION level, observation noise can correspondingly weaken, and observation noise
Energy be distributed on each frequency band than more uniform, the Wavelet Packet Energy Spectrum using the ratio between energy and the gross energy of each frequency band construction has
There is the wavelet packet of j-th of frequency band of preferred definition above formula in susceptibility to damage and noise robustness, therefore technical solution of the present invention
Energy ratio is used as damage location identification index.One WAVELET PACKET DECOMPOSITION layer can confirm that according to a series of energy feature frequency band
The average energy value of whole frequency bands, whole frequency bands that energy ratio is confirmed as to the energy and decomposition layer where the frequency band of a certain frequency band are equal
Value is compared, and can so be effectively reduced influence of the noise to result of calculation.
In the embodiment of technical solution of the present invention, preferably by the energy ratio I of each feature bandjNormalize and carry out by maximum
Adjustment, calculates energy ratio IjChange (i.e. Energy Ratio Variation, abbreviation ERV), i.e.,:
In formula, ERVjFor the energy ratio change of j-th of feature band;IdjFor structure under faulted condition j-th frequency band
Energy ratio,Represent j-th frequency band energy of the not damaged structure under multiple excitations than average value.Technical solution of the present invention
Embodiment in, compare be a certain frequency band energy ratio and the frequency band under not damaged structure, multiple excitations energy ratio it is flat
Average.The method avoids the influence that energy ratio error is brought for judged result so that is used as the energy of benchmark
Than with higher confidence level, so that further such that the judged result for damage has higher accuracy rate.
Adjustment is normalized by maximum in the energy ratio of each feature band, passes through the change of energy ratio, it can be determined that
The degree of impairment of structure.In the embodiment of technical solution of the present invention, according to the energy ratio, each spy can be further calculated
The energy ratio deviation of frequency band is levied, deviation is bigger, and degree of injury is bigger, so as to be estimated that degree of injury with reference to other specification
Size.In the embodiment of technical solution of the present invention, what energy ratio change was compared is that energy ratio and the frequency band of a certain frequency band exists
The average value of energy ratio under not damaged structure, multiple excitations.This processing method can ensure the not damaged state as benchmark
Energy ratio it is relatively accurate, the damage of structure can be made it is more accurate judge.
Signatures for damage detection is turned to compared to becoming using energy ratio or energy ratio, energy ratio deviation has higher survey
Accuracy of measurement, it is more enough to show especially out by small error amplification, contribute to the accurate identification and early warning of damage.
Below by way of the three-dimensional shield tunnel FEM model of foundation as the numerical example, the implementation to technical solution of the present invention
Example describes in detail.
As shown in Fig. 2 the shield tunnel construction of the embodiment analog study of technical solution of the present invention is by shield tunnel, road
Bed, soil body etc. are constituted.It is preferred that setting up the three-dimensional finite element model of railway roadbed-section of jurisdiction-country rock by ABAQUS finite element softwares.Fig. 3
It is the partial enlarged drawing of FEM model schematic diagram.Using tunnel center line position as absolute altitude, the preferred 16m of buried depth, tunnel radius is excellent
Select 2.925m, the preferred 2.525m of lining cutting internal diameter, the preferred 0.4m of section of jurisdiction thickness.The computer capacity of simulation is as follows, it is preferable that in tunnel
300m is taken on length direction, the vertical depth that calculates takes 316m, and horizontal direction takes 360m.The embodiment of technical solution of the present invention is preferred
The vibration of vertical direction is studied, it is unconfined free boundary the soil body to be applied at the top of the fixed constraint of horizontal direction, i.e. soil layer
Condition, the deformation in left and right sides bound level direction, bottom is fixed constraint.Compiled using computer language such as Fortran language
The user's favorite subroutine (such as DLOAD) for writing the offer of ABAQUS secondary developing platforms applies traveling load, and load spreads through wheel load to convert
It is applied directly to afterwards on simplified tunnel railway roadbed.Preferably, the long 300m of tunnel structure, is divided into 1 unit per 1m, is evenly arranged
31 measuring points, measuring point interval 10m.
Damage reason location and degree of injury analysis and assessment are carried out with reference to actual conditions.
Damage reason location is carried out using the method for the present invention to assess with degree of injury, each operating mode is described as follows:
Operating mode one, is designated as UDG, and any damage does not occur for tunnel;
Operating mode two, is designated as DMG1, sets tunnel duct piece element stiffness to decline at unit 100 (respective sensor numbering 11) place
3% damage;
Operating mode three, is designated as DMG2, sets tunnel duct piece element stiffness to decline at unit 100 (respective sensor numbering 11) place
6% damage;
Operating mode four, is designated as DMG3, in unit 100 (respective sensor numbering 11) and (the respective sensor numbering of unit 200
21) place sets the damage of tunnel duct piece element stiffness decline 3% simultaneously;
Operating mode five, is designated as DMG4, in unit 100 (respective sensor numbering 11) and (the respective sensor numbering of unit 180
19) place sets the damage of tunnel duct piece element stiffness decline 3% simultaneously.
Under above-mentioned six kinds of operating condition of test, apply same traveling load excitation on tunnel railway roadbed, obtain sensor volume
Response acceleration at number.
Step 1:Under above-mentioned five kinds of operating modes, the acceleration responsive time-histories of 31 measuring points, i.e. dynamic response are gathered respectively and is believed
Number, sample frequency 100Hz.
Step 2:It is 8 to choose wavelet function DB18, the WAVELET PACKET DECOMPOSITION number of plies, and the acceleration of each measuring point before and after damage is believed
Number carry out WAVELET PACKET DECOMPOSITION.
Step 3:The wavelet-packet energy of each measuring point acceleration signal in step 2 is tried to achieve by formula (2), wavelet packet is constituted
Energy spectrum vector.
Step 4:The energy ratio deviation of damage regime, as damage criterion are tried to achieve by formula (3)-(5).
Step 5:Using Sensor numbering as abscissa, step 4 gained damage criterion is ordinate, draws wavelet packet energy
Amount determines to damage the position of unit than deviation block diagram, according to the size of amplitude.As Figure 4-Figure 7.
Damage reason location is carried out using the method for the present invention to assess with degree of injury.Damage reason location result as Figure 4-Figure 7,
Damage criterion mutation value clearly demonstrates the position of damage in figure, and can be seen that by amplitude difference the difference of degree of injury.
From five kinds of damage regimes of above-described embodiment, it is proposed by the present invention can be correct based on wavelet-packet energy spectrometry
Efficiently identify out damage position and degree of injury.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, it is not used to
The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the invention etc., it all should include
Within protection scope of the present invention.
Claims (7)
1. a kind of shield tunnel construction method for estimating damage, it is characterised in that including,
S1 sets up the shield tunnel construction model of not damaged operating mode according to the parameter of shield tunnel, is further applied load on the mold,
Obtain the not damaged acceleration responsive signal of measuring point;
S2 sets up the shield tunnel construction model with damage regime according to the parameter of shield tunnel, applies described on the mold
Load, obtain measuring point has damage acceleration responsive signal;
S3 is to the not damaged acceleration responsive and has damage acceleration responsive signal to carry out WAVELET PACKET DECOMPOSITION respectively, obtains each
The wavelet-packet energy of frequency band;
S4 calculates the not damaged acceleration responsive signal of the measuring point and has the wavelet-packet energy of damage acceleration responsive signal special
The energy ratio deviation of frequency band is levied, the energy ratio deviation is appointed as signatures for damage detection, according to the signatures for damage detection, really
Determine the degree of impairment of shield tunnel construction.
2. a kind of shield tunnel construction method for estimating damage according to claim 1, wherein, each frequency described in step S3
The wavelet-packet energy of band is defined as
<mrow>
<msub>
<mi>E</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>|</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<msup>
<mn>2</mn>
<mi>i</mi>
</msup>
<mo>;</mo>
</mrow>
In formula, i is the WAVELET PACKET DECOMPOSITION number of plies, SijThe subsignal of j-th of frequency band, x in the frequency band obtained for i-th layer of decompositionjkFor jth
The value of k-th of discrete point of individual frequency band subsignal, N is total points of j-th of frequency band subsignal.
3. a kind of shield tunnel construction method for estimating damage according to claim 1 or 2, wherein, damaged described in step S4
Hinder distinguishing indexes to be defined as
<mrow>
<mi>E</mi>
<mi>R</mi>
<mi>V</mi>
<mi>D</mi>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msup>
<mn>2</mn>
<mi>i</mi>
</msup>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>ERV</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mover>
<mrow>
<mi>E</mi>
<mi>R</mi>
<mi>V</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>;</mo>
</mrow>
In formula, i is the WAVELET PACKET DECOMPOSITION number of plies, ERVjChange for the energy ratio of j-th of frequency band,It is ERVjAverage value.
4. a kind of shield tunnel construction method for estimating damage according to any one of claims 1 to 3, wherein, the identification
Index can determine the position of structural damage, and the relative degree for determining structural damage.
5. a kind of shield tunnel construction method for estimating damage according to any one of Claims 1 to 4, wherein, it is described to accelerate
Degree response signal is preferably calculated by ABAQUS finite element softwares to be obtained.
6. a kind of shield tunnel construction method for estimating damage according to any one of Claims 1 to 5, wherein, it is described to accelerate
Spend the preferred 100Hz of response signal sample frequency.
7. a kind of shield tunnel construction method for estimating damage according to any one of claim 1~6, wherein, the small echo
The bag function decomposition number of plies is preferably 8 layers.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710434506.7A CN107228772A (en) | 2017-06-09 | 2017-06-09 | A kind of shield tunnel construction method for estimating damage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710434506.7A CN107228772A (en) | 2017-06-09 | 2017-06-09 | A kind of shield tunnel construction method for estimating damage |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107228772A true CN107228772A (en) | 2017-10-03 |
Family
ID=59935942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710434506.7A Pending CN107228772A (en) | 2017-06-09 | 2017-06-09 | A kind of shield tunnel construction method for estimating damage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107228772A (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109815940A (en) * | 2019-03-05 | 2019-05-28 | 韦灼彬 | Wavelet-packet energy spectrometry damnification recognition method |
CN109840386A (en) * | 2019-03-05 | 2019-06-04 | 韦灼彬 | Damnification recognition method based on factorial analysis |
CN109871645A (en) * | 2019-03-07 | 2019-06-11 | 汕头大学 | A kind of Structural Damage Identification based on mutative scale Symbolic time series analysis |
CN109903525A (en) * | 2019-01-23 | 2019-06-18 | 昆明理工大学 | A kind of Rock And Soil unstability disaster monitoring method for early warning based on external exciting and Wavelet Packet Energy Spectrum |
CN109975136A (en) * | 2019-04-09 | 2019-07-05 | 西南交通大学 | A kind of steel-frame structure damnification recognition method based on wavelet packet analysis |
CN110017929A (en) * | 2019-03-12 | 2019-07-16 | 华中科技大学 | Vessel bump bridge load and damage method for synchronously recognizing based on minor structure sensitivity analysis |
CN111597617A (en) * | 2020-05-18 | 2020-08-28 | 中铁第六勘察设计院集团有限公司 | Method for predicting fatigue life of shield tunnel |
CN112529842A (en) * | 2020-11-23 | 2021-03-19 | 天津大学 | Multi-excitation fusion plate structure damage identification method based on wavelet packet energy |
CN113627048A (en) * | 2021-07-12 | 2021-11-09 | 暨南大学 | Method for quickly identifying structural damage based on local transfer rate function and pattern matching |
CN114755122A (en) * | 2022-04-19 | 2022-07-15 | 西南交通大学 | Testing device and testing method for subway tunnel structure full life cycle health monitoring |
CN117191952A (en) * | 2023-11-06 | 2023-12-08 | 中冶建筑研究总院有限公司 | Fatigue damage identification and life prediction method based on acoustic emission signal wavelet packet decomposition frequency band energy spectrum |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1809999A1 (en) * | 2004-09-16 | 2007-07-25 | Sekos, Inc. | Method for detecting the onset of structural collapse of burning structures |
CN103472141A (en) * | 2013-09-09 | 2013-12-25 | 北京工业大学 | Signal demodulation method for recognizing vibrant noise modulation mechanism |
CN103983692A (en) * | 2014-05-29 | 2014-08-13 | 重庆大学 | Damage alarming method for retaining wall structure system |
CN104458173A (en) * | 2014-11-27 | 2015-03-25 | 广东电网有限责任公司中山供电局 | Steel framework structure mutational damage recognition method and system |
CN106768763A (en) * | 2017-01-09 | 2017-05-31 | 温州大学 | A kind of dash-board injury method for quick |
-
2017
- 2017-06-09 CN CN201710434506.7A patent/CN107228772A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1809999A1 (en) * | 2004-09-16 | 2007-07-25 | Sekos, Inc. | Method for detecting the onset of structural collapse of burning structures |
CN103472141A (en) * | 2013-09-09 | 2013-12-25 | 北京工业大学 | Signal demodulation method for recognizing vibrant noise modulation mechanism |
CN103983692A (en) * | 2014-05-29 | 2014-08-13 | 重庆大学 | Damage alarming method for retaining wall structure system |
CN104458173A (en) * | 2014-11-27 | 2015-03-25 | 广东电网有限责任公司中山供电局 | Steel framework structure mutational damage recognition method and system |
CN106768763A (en) * | 2017-01-09 | 2017-05-31 | 温州大学 | A kind of dash-board injury method for quick |
Non-Patent Citations (2)
Title |
---|
丁幼亮 等: "基于小波包分析的Benchmark结构损伤预警试验研究", 《工程力学》 * |
刘涛 等: "基于小波包能领普的结构损伤预警方法试验研究", 《振动与冲击》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109903525A (en) * | 2019-01-23 | 2019-06-18 | 昆明理工大学 | A kind of Rock And Soil unstability disaster monitoring method for early warning based on external exciting and Wavelet Packet Energy Spectrum |
CN109840386A (en) * | 2019-03-05 | 2019-06-04 | 韦灼彬 | Damnification recognition method based on factorial analysis |
CN109815940A (en) * | 2019-03-05 | 2019-05-28 | 韦灼彬 | Wavelet-packet energy spectrometry damnification recognition method |
CN109871645A (en) * | 2019-03-07 | 2019-06-11 | 汕头大学 | A kind of Structural Damage Identification based on mutative scale Symbolic time series analysis |
CN110017929A (en) * | 2019-03-12 | 2019-07-16 | 华中科技大学 | Vessel bump bridge load and damage method for synchronously recognizing based on minor structure sensitivity analysis |
CN109975136B (en) * | 2019-04-09 | 2021-06-25 | 西南交通大学 | Steel frame structure damage identification method based on wavelet packet analysis |
CN109975136A (en) * | 2019-04-09 | 2019-07-05 | 西南交通大学 | A kind of steel-frame structure damnification recognition method based on wavelet packet analysis |
CN111597617A (en) * | 2020-05-18 | 2020-08-28 | 中铁第六勘察设计院集团有限公司 | Method for predicting fatigue life of shield tunnel |
CN111597617B (en) * | 2020-05-18 | 2022-05-27 | 中铁第六勘察设计院集团有限公司 | Method for predicting fatigue life of shield tunnel |
CN112529842A (en) * | 2020-11-23 | 2021-03-19 | 天津大学 | Multi-excitation fusion plate structure damage identification method based on wavelet packet energy |
CN113627048A (en) * | 2021-07-12 | 2021-11-09 | 暨南大学 | Method for quickly identifying structural damage based on local transfer rate function and pattern matching |
CN113627048B (en) * | 2021-07-12 | 2022-12-06 | 暨南大学 | Method for quickly identifying structural damage based on local transfer rate function and pattern matching |
CN114755122A (en) * | 2022-04-19 | 2022-07-15 | 西南交通大学 | Testing device and testing method for subway tunnel structure full life cycle health monitoring |
CN114755122B (en) * | 2022-04-19 | 2023-09-01 | 西南交通大学 | Test device and test method for full life cycle health monitoring of subway tunnel structure |
CN117191952A (en) * | 2023-11-06 | 2023-12-08 | 中冶建筑研究总院有限公司 | Fatigue damage identification and life prediction method based on acoustic emission signal wavelet packet decomposition frequency band energy spectrum |
CN117191952B (en) * | 2023-11-06 | 2024-01-23 | 中冶建筑研究总院有限公司 | Fatigue damage identification and life prediction method based on acoustic emission signal wavelet packet decomposition frequency band energy spectrum |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107228772A (en) | A kind of shield tunnel construction method for estimating damage | |
CN103852269B (en) | Bullet train runs kinetic parameter detection method | |
CN103293014B (en) | Bridge fatigue damage state and residual life evaluating method | |
WO2017202139A1 (en) | Bridge damage identification method based on long-gauge-length strain influence envelope | |
Kildashti et al. | Drive-by-bridge inspection for damage identification in a cable-stayed bridge: Numerical investigations | |
CN106197910A (en) | A kind of Bridges Detection based on Vehicle-bridge Coupling Analysis and detecting system | |
CN104090977A (en) | Random recognition method for bridge floor moving vehicle loads | |
CN106323442A (en) | Railway health monitoring method based on distributed optical fiber vibration sensing system | |
CN108698622A (en) | Track condition evaluation method, device and program | |
CN101266190A (en) | Normal traffic flow stayd-cable bridge mode parametric measurement device and method | |
CN106198058A (en) | A kind of vertical wheel impacts power method for real-time measurement based on tire pressure monitoring | |
CN104792937A (en) | Bridge head bump detection evaluation method based on vehicle-mounted gravitational acceleration sensor | |
CN106156336A (en) | A kind of Cable-Stayed Bridge Structure state evaluation system and assessment method | |
CN104568493B (en) | The quick damnification recognition method of structure based on displacement time-histories area under vehicular load | |
CN102175768A (en) | Method and device for detecting defects and failures of high-speed rail based on vibration signals | |
CN104123463A (en) | Time domain identification method of random dynamic loads | |
CN109002673A (en) | A kind of bridge foundation washout recognition methods based on vehicle braking percussion | |
CN104133960A (en) | Improved optimal arranging method of static sensors | |
Zolghadri et al. | Field verification of simplified bridge weigh-in-motion techniques | |
CN104833535A (en) | Railway vehicle tire tread scratch detection method | |
Li et al. | Modeling dynamic responses of a cross-river road shield tunnel under stochastic vehicle loads | |
CN110926735A (en) | Bridge structure rapid diagnosis method based on multidimensional dynamic parameters | |
CN106526149A (en) | Method for predicting pavement performance on basis of transport durations and traffic volumes | |
CN102998133A (en) | Energy damage identification method based on quasi-distributed acceleration data | |
Liu et al. | Turning telecommunication fiber-optic cables into distributed acoustic sensors for vibration-based bridge health monitoring |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171003 |