CN106248326A - The damage location identification method of high-rise frame structure based on natural frequency rate of change - Google Patents
The damage location identification method of high-rise frame structure based on natural frequency rate of change Download PDFInfo
- Publication number
- CN106248326A CN106248326A CN201610589651.8A CN201610589651A CN106248326A CN 106248326 A CN106248326 A CN 106248326A CN 201610589651 A CN201610589651 A CN 201610589651A CN 106248326 A CN106248326 A CN 106248326A
- Authority
- CN
- China
- Prior art keywords
- damage
- frequency
- change rate
- rate
- frequency change
- 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
- 230000006378 damage Effects 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000001133 acceleration Effects 0.000 claims abstract description 27
- 239000000284 extract Substances 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000015556 catabolic process Effects 0.000 claims description 5
- 238000006731 degradation reaction Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 4
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000009435 building construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000004643 material aging Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0066—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by exciting or detecting vibration 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
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses a kind of damage location identification method of high-rise frame structure based on natural frequency rate of change, comprise the following steps: extract different damage position index PLUnder frequency change rate kf；Set up frequency change rate kfWith damage position index PLRelational expression；Measured the acceleration responsive value of a certain layer of high-rise frame structure by acceleration transducer, and this acceleration responsive value is extracted frequency change rate kf；It is updated to frequency change rate kfWith damage position index PLRelational expression, inverse goes out damage degree index PL.Response signal before existing damnification recognition method must use structure currently to respond signal and structural damage is analyzed the position that just can recognize that damage, and the present invention has only to use structure currently to respond signal and is analyzed just can identifying the damage position of structure.The damage position of structure is formulated by the present invention with frequency change rate, utilizes the formula of this invention can calculate the damage position of structure accurately.
Description
Technical field
The present invention relates to locations of structures damnification recognition method, the damage location identification side of a kind of high-rise frame structure
Method.
Background technology
In current building construction, the application of high-rise frame structure is more and more extensive, but along with material aging, environment corruption
The effects such as erosion, load effect, artificial destruction and natural mutation effect, these effects will inevitably cause structure
Damage accumulation and degradation resistance, therefore take the means of science, from durability and the angle of safety of structure, to high-rise frame
The degree of impairment of structure carries out detecting, identify, assess extremely the most urgent.
It is currently based on the recognition methods of high-rise frame structure damage position, is all by the response signal before and after damage is entered
Row relative analysis, thus the damage position of location structure.
But, in real life, always cause the loss of data before structural damage because of various reasons, so
Structure frequency before identification of damage is caused huge interference, thus inevitably has influence on the knowledge of structural damage information
Not.
If the natural frequency rate of change using structure identifies the damage of structure, it is not necessary to the response letter before structural damage
Number, it is only necessary to provide the response signal that structure is current, it is possible to the damage position of analytical structure, so greatly reduce damage and know
Other workload.How to identify the research of the damage position of structure about natural frequency rate of change, yet there are no relevant report.
Summary of the invention
For solving the problem that existing method exists, the present invention to provide a kind of response signal essence that structure can be utilized current
Really identify the damage location identification method of the high-rise frame structure based on natural frequency rate of change of the damage position of structure.
For achieving the above object, the technical solution used in the present invention is as follows: high-rise frame based on natural frequency rate of change
The damage location identification method of structure, comprises the following steps:
A, set up frequency change rate k by emulation modefWith damage position index PLRelational expression
A1, different damage position index P of extractionLUnder frequency change rate kf
Utilize numerical value emulation method obtain different layers Stiffness degradation rate identical time structure certain layer acceleration responsive value,
This acceleration responsive value is extracted frequency change rate kf；
Frequency change rate
In formula: df is the variable quantity of consecutive sampling times inner structure frequency, dt is sampling time interval；
Damage position index
In formula: cnFor the damage layer of structure, c is total number of plies of structure；
A2, set up frequency change rate kfWith damage position index PLRelational expression
The frequency change rate k that step A1 is obtainedfUse least-square fitting approach, obtain frequency change rate kfWith damage
Positioning index PLRelational expression
G(PL,kf)=k'PL-kf+C2=0 (1)
In formula: kfFor frequency change rate, PLFor damage position index, k', C2For constant.
After B, utilization damage, the acceleration responsive value inverse of actual measurement goes out damage position index PL
B1, measured the acceleration responsive value of a certain layer of high-rise frame structure by acceleration transducer, and to this acceleration
Response value extracts frequency change rate kf；
B2, the frequency change rate k that step B1 is obtainedfSubstitution formula (1), inverse goes out damage degree index PL, calculate publicity such as
Under:
PL=(kF-C2)/k' (2)
Thus obtain damage position index P of structureL。
Further, described extraction frequency change rate kfMethod, comprise the following steps:
S1, data prediction
Time-Frequency Analysis Method is used to obtain the graphics of T/F-amplitude acceleration responsive value.Then extract each
The frequency corresponding to first peak point of amplitude in the frequency-amplitude X-Y scheme in moment, the single order of the most each moment structure is solid
There is frequency；
S2, extraction frequency change rate kf
Use least-square fitting approach, the first natural frequency in each moment obtained in step S1 is carried out Linear Quasi
Conjunction obtains the time dependent function of frequency:
F (t)=kf·t+C1 (3)
In formula: f (t) is frequency, t is the time, C1For constant；Frequency change rate k is i.e. obtained by formula (3)fValue.
Further, the Time-Frequency Analysis Method described in step S1 is Short Time Fourier Transform method in short-term.
Further, described numerical value emulation method include matlab numerical value emulation method, ansys numerical value emulation method or
Abaqus numerical value emulation method.
Further, the certain layer described in step A1 is to be provided with the place layer of acceleration transducer.
Further, described acceleration transducer is arranged on any layer.
Compared with the conventional method, the present invention has the following advantages:
1, the response signal before existing damnification recognition method must use structure currently to respond signal and structural damage enters
Row relative analysis just can recognize that the position of damage, and the present invention has only to use structure currently to respond signal and is analyzed the most permissible
Identify the damage position of structure.
2, existing damnification recognition method can only identify the damage position of structure by the difference of signal before and after damage, this
Invent and the damage position of structure is formulated with frequency change rate, utilize the formula of this invention can calculate structure accurately
Damage position.
Accompanying drawing explanation
Fig. 1 is the flow chart of the present invention.
Fig. 2 is to extract frequency change rate kfFlow chart.
Fig. 3 is frequency change rate kfWith damage position index PLCurve chart.
Detailed description of the invention
Illustrate how to use the damage position of first natural frequency rate of change identification high-rise frame structure below in conjunction with the accompanying drawings.
For high-rise frame structure, when the Stiffness degradation rate of different layers is identical, low layer damage damages relative to high level
Bigger to structural damage, so identical degree of injury, damage layer is the highest, and the damage factor of structure is the least, thus structure frequency
Rate rate of change is the biggest.By substantial amounts of test data analyzer is found that frequency change rate linearly closes with the damage position of structure
System.
Embodiment:
Practical structures is the high-rise frame structure of Floor 12, it is assumed that the Stiffness degradation rate of a certain layer of known structure is 70% He
The acceleration responsive value that the acceleration transducer of the tenth layer obtains, seeks the damage position of structure.
According to the flow process shown in Fig. 1-2, matlab numerical value emulation method is used to simulate certain Floor 12 high-rise frame structure, point
The not rigidity at three to the Floor 12 of structure reduces input El Centro seismic wave under ten operating modes such as 70% grade and obtains the tenth layer
Acceleration responsive value, to the tenth layer of acceleration responsive value under each operating mode according to step A extract frequency change rate kf, result
As shown in table 1.Draw PL-kfFigure and fitting result are as shown in Figure 3.
Thus obtain PL-kfRelational expression as follows:
PL=(kf+0.0128)/0.0095
The acceleration responsive value that the acceleration transducer of practical structures the tenth layer is obtained according to step A obtain frequency
Rate of change kf=-0.01207 substitution above formula obtains
PL=(-0.01207+0.0128)/0.0095=0.077
Damage layer cn=PL× c=0.077 × 12=0.924 ≈ 1, so, the damage position of structure is at ground floor.
Table 1PL-kfTables of data
| PL | 3/12 | 4/12 | 5/12 | 6/12 | 7/12 | 8/12 | 9/12 | 10/12 | 11/12 | 1 |
| kf | -0.0106 | -0.0097 | -0.0088 | -0.008 | -0.0071 | -0.0062 | -0.0055 | -0.0047 | -0.0041 | -0.0036 |
The present invention is not limited to the present embodiment, any equivalent concepts in the technical scope of present disclosure or change
Become, be all classified as protection scope of the present invention.
Claims (6)
1. the damage location identification method of high-rise frame structure based on natural frequency rate of change, it is characterised in that: include following
Step:
A, set up frequency change rate k by emulation modefWith damage position index PLRelational expression
A1, different damage position index P of extractionLUnder frequency change rate kf
Utilize numerical value emulation method obtain different layers Stiffness degradation rate identical time structure certain layer acceleration responsive value, to this
Acceleration responsive value extracts frequency change rate kf；
Frequency change rate
In formula: df is the variable quantity of consecutive sampling times inner structure frequency, dt is sampling time interval；
Damage position index
In formula: cnFor the damage layer of structure, c is total number of plies of structure；
A2, set up frequency change rate kfWith damage position index PLRelational expression
The frequency change rate k that step A1 is obtainedfUse least-square fitting approach, obtain frequency change rate kfWith damage position
Index PLRelational expression
G(PL,kf)=k'PL-kf+C2=0 (1)
In formula: kfFor frequency change rate, PLFor damage position index, k', C2For constant；
After B, utilization damage, the acceleration responsive value inverse of actual measurement goes out damage position index PL
B1, measured the acceleration responsive value of a certain layer of high-rise frame structure by acceleration transducer, and to this acceleration responsive
Value extracts frequency change rate kf；
B2, the frequency change rate k that step B1 is obtainedfSubstitution formula (1), inverse goes out damage degree index PL, calculate publicity as follows:
PL=(kf-C2)/k' (2)
Thus obtain damage position index P of structureL。
The damage location identification method of high-rise frame structure based on natural frequency rate of change the most according to claim 1,
It is characterized in that: described extraction frequency change rate kfMethod, comprise the following steps:
S1, data prediction
Time-Frequency Analysis Method is used to obtain the graphics of T/F-amplitude acceleration responsive value；Then each moment is extracted
Frequency-amplitude X-Y scheme in the frequency corresponding to first peak point of amplitude, the intrinsic frequency of single order of the most each moment structure
Rate；
S2, extraction frequency change rate kf
Use least-square fitting approach, the first natural frequency in each moment obtained in step S1 is carried out linear fit and obtains
To the time dependent function of frequency:
F (t)=kf·t+C1 (3)
In formula: f (t) is frequency, t is the time, C1For constant；Frequency change rate k is i.e. obtained by formula (3)fValue.
The damage location identification method of high-rise frame structure based on natural frequency rate of change the most according to claim 2,
It is characterized in that: the Time-Frequency Analysis Method described in step S1 is Short Time Fourier Transform method.
The damage location identification method of high-rise frame structure based on natural frequency rate of change the most according to claim 1,
It is characterized in that: described numerical value emulation method includes matlab numerical value emulation method, ansys numerical value emulation method or abaqus
Numerical value emulation method.
The damage location identification method of high-rise frame structure based on natural frequency rate of change the most according to claim 1,
It is characterized in that: the certain layer described in step A1 is to be provided with the place layer of acceleration transducer.
The damage location identification method of high-rise frame structure based on natural frequency rate of change the most according to claim 1,
It is characterized in that: described acceleration transducer is arranged on any layer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610589651.8A CN106248326B (en) | 2016-07-26 | 2016-07-26 | The damage location identification method of high-rise frame structure based on intrinsic frequency change rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610589651.8A CN106248326B (en) | 2016-07-26 | 2016-07-26 | The damage location identification method of high-rise frame structure based on intrinsic frequency change rate |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106248326A true CN106248326A (en) | 2016-12-21 |
| CN106248326B CN106248326B (en) | 2018-06-22 |
Family
ID=57604463
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610589651.8A Active CN106248326B (en) | 2016-07-26 | 2016-07-26 | The damage location identification method of high-rise frame structure based on intrinsic frequency change rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106248326B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110163134A (en) * | 2019-05-10 | 2019-08-23 | 大连理工大学 | A kind of structural damage area recognizing method based on weighted band-wise least square |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3733555A (en) * | 1970-08-19 | 1973-05-15 | Marshall Research & Dev Corp | Shock spectrum analyzer and synthesizer system |
| JP3358944B2 (en) * | 1996-08-08 | 2002-12-24 | 株式会社東芝 | Analysis apparatus, analysis method, frequency response analysis apparatus, and frequency response analysis method |
| CN103106305A (en) * | 2013-02-01 | 2013-05-15 | 北京工业大学 | Space grid structure model step-by-step correction method based on actual measurement mode |
| CN103116759A (en) * | 2013-02-01 | 2013-05-22 | 北京工业大学 | Method applicable to identifying damage to space grid structure |
| CN104200265A (en) * | 2014-07-28 | 2014-12-10 | 东北大学 | Improved bridge damage identification method based on neural network |
| CN104200004A (en) * | 2014-07-28 | 2014-12-10 | 东北大学 | Optimized bridge damage identification method based on neural network |
| CN104239736A (en) * | 2014-09-25 | 2014-12-24 | 广西科技大学 | Structure damage diagnosis method based on power spectrum and intelligent algorithms |
| CN104964837A (en) * | 2015-06-12 | 2015-10-07 | 广东电网有限责任公司电力科学研究院 | Method and system for monitoring rigidity damage of structure based on EMD |
| CN105547625A (en) * | 2016-01-15 | 2016-05-04 | 华北电力科学研究院有限责任公司 | Porcelain sleeve modal testing method, device and system |
-
2016
- 2016-07-26 CN CN201610589651.8A patent/CN106248326B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3733555A (en) * | 1970-08-19 | 1973-05-15 | Marshall Research & Dev Corp | Shock spectrum analyzer and synthesizer system |
| JP3358944B2 (en) * | 1996-08-08 | 2002-12-24 | 株式会社東芝 | Analysis apparatus, analysis method, frequency response analysis apparatus, and frequency response analysis method |
| CN103106305A (en) * | 2013-02-01 | 2013-05-15 | 北京工业大学 | Space grid structure model step-by-step correction method based on actual measurement mode |
| CN103116759A (en) * | 2013-02-01 | 2013-05-22 | 北京工业大学 | Method applicable to identifying damage to space grid structure |
| CN104200265A (en) * | 2014-07-28 | 2014-12-10 | 东北大学 | Improved bridge damage identification method based on neural network |
| CN104200004A (en) * | 2014-07-28 | 2014-12-10 | 东北大学 | Optimized bridge damage identification method based on neural network |
| CN104239736A (en) * | 2014-09-25 | 2014-12-24 | 广西科技大学 | Structure damage diagnosis method based on power spectrum and intelligent algorithms |
| CN104964837A (en) * | 2015-06-12 | 2015-10-07 | 广东电网有限责任公司电力科学研究院 | Method and system for monitoring rigidity damage of structure based on EMD |
| CN105547625A (en) * | 2016-01-15 | 2016-05-04 | 华北电力科学研究院有限责任公司 | Porcelain sleeve modal testing method, device and system |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110163134A (en) * | 2019-05-10 | 2019-08-23 | 大连理工大学 | A kind of structural damage area recognizing method based on weighted band-wise least square |
| CN110163134B (en) * | 2019-05-10 | 2022-03-29 | 大连理工大学 | Structural damage region identification method based on sub-band weighted least square |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106248326B (en) | 2018-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102156873B (en) | Chaos-based method for detecting and classifying early single-point faults of mechanical component | |
| CN108399277B (en) | Bridge damage identification method based on temperature and strain correlation | |
| CN104331595B (en) | The mobile principal component correlation coefficient process of bridge damnification early warning | |
| EA201370085A1 (en) | NEW SYSTEMS AND METHODS OF NON-DESTRUCTIVE AIRCRAFT CONTROL | |
| CN101996292B (en) | Method for analyzing safety property of software based on sequence clustering | |
| CN103913512B (en) | The damage reason location system of suspension cable periodic detection | |
| Jin et al. | Adaptive reference updating for vibration-based structural health monitoring under varying environmental conditions | |
| CN103262048A (en) | Operation management device, operation management method, and program | |
| CN106482639B (en) | The low velocity impact position identifying method calculated based on approximate entropy | |
| CN108956111A (en) | A kind of the abnormal state detection method and detection system of mechanical part | |
| CN104636319A (en) | Text duplicate removal method and device | |
| CN102930207A (en) | API log monitoring method and device | |
| CN103198243A (en) | Method and device for identifying emulational application programs | |
| CN104502450A (en) | Simple pipe installation stress identification method | |
| Soo Lon Wah et al. | Damage detection under temperature conditions using PCA–an application to the Z24 Bridge | |
| CN104180823A (en) | Temperature compensation method and device | |
| CN106706109B (en) | Vibration source identification method and system based on time domain two-dimensional characteristics | |
| CN106248326A (en) | The damage location identification method of high-rise frame structure based on natural frequency rate of change | |
| Zhang et al. | Automatic identification of structural modal parameters based on density peaks clustering algorithm | |
| CN110781612A (en) | Fault diagnosis method and device for ball screw, computer device and storage medium | |
| CN102778632A (en) | Double normalization recognition method for directly forecasting and recognizing transformer winding fault type | |
| CN106227941B (en) | The damage extent identification method of high-rise frame structure based on intrinsic frequency change rate | |
| CN109669875B (en) | Application crash test method and device, electronic equipment and storage medium | |
| CN106708786A (en) | Method and system for calculating problem severity of iron tower based on sensor detection | |
| CN105721459A (en) | Risk evaluation method for virtual environment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20161221 Assignee: Dalian Yude Technology Co.,Ltd. Assignor: DALIAN University Contract record no.: X2023210000014 Denomination of invention: Damage location identification method of high-rise frame structure based on natural frequency change rate Granted publication date: 20180622 License type: Common License Record date: 20230215 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20161221 Assignee: Dalian Big Data Industry Development Research Institute Co.,Ltd. Assignor: DALIAN University Contract record no.: X2023210000224 Denomination of invention: Damage location identification method for high-rise frame structures based on natural frequency variation rate Granted publication date: 20180622 License type: Common License Record date: 20231129 |
|
| OL01 | Intention to license declared | ||
| OL01 | Intention to license declared |