CN102998367A - Damage identification method based on virtual derivative structure - Google Patents

Damage identification method based on virtual derivative structure Download PDF

Info

Publication number
CN102998367A
CN102998367A CN2012105610813A CN201210561081A CN102998367A CN 102998367 A CN102998367 A CN 102998367A CN 2012105610813 A CN2012105610813 A CN 2012105610813A CN 201210561081 A CN201210561081 A CN 201210561081A CN 102998367 A CN102998367 A CN 102998367A
Authority
CN
China
Prior art keywords
virtual
source structure
virtual derivative
damage
characteristic frequency
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
Application number
CN2012105610813A
Other languages
Chinese (zh)
Other versions
CN102998367B (en
Inventor
杨秋伟
杨丽君
李翠红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shaoxing
Original Assignee
University of Shaoxing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Shaoxing filed Critical University of Shaoxing
Priority to CN201210561081.3A priority Critical patent/CN102998367B/en
Publication of CN102998367A publication Critical patent/CN102998367A/en
Application granted granted Critical
Publication of CN102998367B publication Critical patent/CN102998367B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a damage identification method based on a virtual derivative structure. The damage identification method comprises the following steps of: firstly, establishing a finite element model of a to-be-tested structure (also known as a source structure) by finite element software; secondly, acquiring front several modal data of vibration of the source structure by testing of a modal testing system; thirdly, constructing a series of virtual derivative structures based on the source structure and counting corresponding feature frequency parameters of the virtual derivative structures; and finally, synthesizing feature frequencies of the source structure and all the virtual derivative structures, and obtaining a source structure damage identification result containing a damage position and a degree by a hybrid frequency sensitivity diagnostic program. According to the technical scheme, the test equipment is simple, the cost is low, the defect of the existing testing technique is effectively overcome, the accuracy of the damage identification is substantially improved, and the damage identification method is particularly suitable for being applied in the damage identification of large-sized structures.

Description

A kind of damnification recognition method based on virtual derivative strucure
Technical field
The invention belongs to the Damage Assessment Method field of civil engineering work subject, relate to a kind of damnification recognition method based on virtual derivative strucure.
Background technology
Along with developing rapidly of national economy, large quantities of large-scale civil engineering structures have been built, such as high-rise building, Longspan Bridge, large span space roof truss etc.Because the growth of tenure of use, and the impact of environmental corrosion, Hazard Loads etc., these structures will inevitably be damaged, the local damage of structure may cause the rapid destruction of structural entity, bring massive losses to people's lives and properties, must make promptly and accurately identification to the damage status in the structure, so that the consequence of averting a calamity property of maintenance and reinforcement occurs.
Mostly traditional damnification recognition method is the experimental technique of visual part, such as ultrasonic method, methods of magnetic field, temperature field method etc.Yet all there is following major defect in said method: at first, said method requires probably to know in advance damage position, so anticipation is relatively poor; Secondly said method requires detected position equipment to arrive, so workload expends greatly height; The 3rd, because there is as above defective in said method, therefore be not suitable for the damage identification of large scale civil engineering structure.
In the last few years, research method whole, real-time, that be applicable to large scale structure damage identification had become the common focus of paying close attention to of numerous field of engineering technology such as building, machinery, aviation.Wherein, the method for utilizing the variation of large scale structure mode of oscillation data to come inverting to identify its damage is present a kind of new technology.The ultimate principle of these class methods is: the Produced by Modal data are functions of structural physical parameter (such as quality, rigidity etc.), therefore the variation (structural damage) of structural physical parameter must cause the variation of Produced by Modal parameter (frequency and the vibration shape), utilizes the variation of these vibration datas just can identify conversely the damage status of structure.Yet, because the limitation of measuring technology, can only measure at present the front several lower modes that obtain structural vibration, this can not satisfy the needs of large scale structure damage identification far away, in other words, the identification equation number that can be set up by the lower mode that measure to obtain is far smaller than the number of the impairment parameter of the unknown, and this is to cause present method to damage one of the reason of the most critical of recognition failures.
In view of this, the inventor is in conjunction with being engaged in Damage Assessment Method area research work experience for many years, and the defective of above-mentioned technical field is studied for a long period of time, and this case produces thus.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of damnification recognition method based on virtual derivative strucure, and is easy and simple to handle, and testing apparatus is simple, with low cost, improved the accuracy of damage identification.
To achieve these goals, technical scheme of the present invention is as follows:
A kind of damnification recognition method based on virtual derivative strucure may further comprise the steps: at first, set up by finite element software and to treat the not finite element model under the faulted condition of geodesic structure (being called again source structure); Secondly, measure former rank mode (being characteristic frequency and the corresponding vibration shape) data that obtain the source structure vibration by the mode test macro, the rank number of mode that needs to measure considers definite according to the complexity of structure, precision and the test site condition of surveying instrument; Then, a series of virtual derivative strucure and calculate corresponding characteristic frequency parameter take source structure as base configuration; At last, the characteristic frequency of comprehensive source structure and all virtual derivative strucures by hybrid frequency sensitivity diagnostic routine, draws the source structure damage recognition result that comprises damage position and degree.
Described mode test macro comprises a data acquisition instrument, a plurality of wireless senser and a cover mode analysis software, is measured former rank characteristic frequency and the corresponding Data of Mode that obtains the source structure vibration by the mode test macro; Wherein wireless senser is arranged in the position for the treatment of to occur easily in the geodesic structure damage, the quantity of wireless senser
Figure 204369DEST_PATH_IMAGE001
Can be by formula
Figure 898655DEST_PATH_IMAGE002
Determine, wherein
Figure 81375DEST_PATH_IMAGE003
Be the sum of unit in the finite element model,
Figure 907117DEST_PATH_IMAGE004
Number for the characteristic frequency measured.
Further, described former rank modal data refers to modal data between front 4 ~ 10 rank.
Further, described virtual derivative strucure can obtain by the way of virtual arrangement quality on source structure or rigidity, the position of virtual mass or rigidity is the position at wireless senser place, can arrange separately also and can make up layout, the size of virtual mass or rigidity be taken as source structure not under the faulted condition in the finite element model mass matrix or stiffness matrix corresponding to 10% ~ 15% of sensing station place numerical value; Adopt different virtual arrangement schemes just can obtain a series of virtual derivative strucure, the sum of the virtual derivative strucure that can construct and the quantity of sensor
Figure 833485DEST_PATH_IMAGE001
Between the pass be
Figure 382278DEST_PATH_IMAGE005
After virtual derivative strucure is constructed, adopt again the fundamental frequency sensitivity formula to calculate the former rank characteristic frequency parameter that obtains virtual derivative strucure.
Further, described Mixed Sensitivity diagnostic routine may further comprise the steps: at first, according to source structure and the virtual derivative strucure finite element model under the faulted condition not, calculate respectively source structure and virtual derivative strucure under the faulted condition not before
Figure 470320DEST_PATH_IMAGE004
Rank characteristic frequency and corresponding fundamental frequency sensitivity; Then before will measuring the source structure of gained
Figure 268511DEST_PATH_IMAGE004
Before rank characteristic frequency and each the virtual derivative strucure
Figure 201832DEST_PATH_IMAGE004
The rank characteristic frequency combines, and lists the one order equation of damage front and back frequency variation, solves all unknown impairment parameter by generalizde inverse
Figure 870711DEST_PATH_IMAGE006
The impairment parameter output that to calculate gained by forms such as histograms at last can accurately identify damage position and the degree of source structure.
In the technical program, only utilize the minority sensor measurement to treat former rank mode of geodesic structure (source structure) vibration, obtain more frequency information by constructing a series of virtual derivative strucure, the frequency of last joint source structure and all derivative strucures identifies the damage status of source structure.
Adopt the technical program, obtain following beneficial effect: the first, the present invention only needs the sensor seldom counted, and front several lower modes of only measuring structural vibration can carry out, and testing apparatus is simple, with low cost; Second, the present invention obtains more frequency parameter by virtual a series of derivative strucure and is used for Damage Assessment Method, effectively overcome the deficiency of existing measuring technology, thereby increased substantially the accuracy of damage identification, be particularly suitable for being applied in the damage identification of large scale structure; The 3rd, derivative strucure of the present invention is virtual construct, so need not treating actual arrangement quality or rigidity member on the geodesic structure, only needs conventional mode testing apparatus and a small amount of sensor to carry out, therefore easy to operate easy to implement, and do not increase extra test job amount; Therefore the 4th, the speed according to power increases rapidly the virtual derivative strucure quantity that the present invention can construct along with the increase of number of sensors, and the used number of sensors of the present invention is minimum under the prerequisite that satisfies damage identification needs.
Technical scheme of the present invention the present invention is described in further detail below in conjunction with drawings and Examples in order further to explain.
Description of drawings
Fig. 1 is enforcement synoptic diagram of the present invention; Label declaration wherein: the 1st, Wireless Acceleration Sensor; The 2nd, the data acquisition instrument; The 3rd, model analysis software; 1, the 2 and 3 common mode test macros that form; The 4th, treat geodesic structure (source structure); The 5th, source structure is the finite element model under the faulted condition not; The 6th, virtual derivative strucure; The 7th, the Mixed Sensitivity diagnostic routine; The 8th, the damage recognition result.
Fig. 2 is a grid structure model to be detected;
Fig. 3 is the recognition result (unit 15 losss of rigidity 15%) of single degree of impairment;
The recognition result of a plurality of degree of impairments of Fig. 4 ( unit 7 and 15 rigidity all lose 15%).
Embodiment
To Fig. 4, further detailed description is done in enforcement of the present invention below in conjunction with accompanying drawing 1.
A kind of damnification recognition method based on virtual derivative strucure may further comprise the steps: at first, set up by finite element software and to treat the not finite element model under the faulted condition of geodesic structure (being called again source structure); Secondly, measured former rank mode (being characteristic frequency and the corresponding vibration shape) data that obtain the source structure vibration by the mode test macro, the rank number of mode of need measuring considers definite according to the precision of the complexity of structure, surveying instrument and test site condition, be taken as in the present embodiment between front 4 ~ 10 rank; Then, a series of virtual derivative strucure and calculate corresponding characteristic frequency parameter take source structure as base configuration; At last, the characteristic frequency of comprehensive source structure and all virtual derivative strucures by hybrid frequency sensitivity diagnostic routine, draws the source structure damage recognition result that comprises damage position and degree.
Described source structure is the finite element model under the faulted condition not, can modeling obtains to source structure by general finite element software (such as general commercial softwares such as ANSYS or MATLAB); Repeat no more herein.In the present embodiment, described finite element model comprises stiffness matrix and mass matrix, does not consider structural damping, and adopts lumped mass matrix.
Described mode test macro comprises a data acquisition instrument, a plurality of Wireless Acceleration Sensor and a cover mode analysis software (they all can be bought by market and obtain), is measured former rank characteristic frequency and the corresponding Data of Mode that obtains the source structure vibration by the mode test macro.Wherein wireless senser is arranged in the position (determining that according to engineering experience can be arranged in span centre and beam-ends section such as girder construction, plate structure can be arranged in four jiaos and center position) for the treatment of to occur easily in the geodesic structure damage, the quantity of wireless senser
Figure 129654DEST_PATH_IMAGE001
Can be by formula Determine, wherein
Figure 152154DEST_PATH_IMAGE003
Be the sum of unit in the finite element model, Number for the characteristic frequency measured.
Described virtual derivative strucure can obtain (being not actual arrangement quality or rigidity on source structure) by the way in virtual arrangement quality or rigidity on the source structure, the position of virtual mass or rigidity is the position at wireless senser place, can arrange separately also and can make up layout, the size of virtual mass or rigidity generally be taken as source structure not under the faulted condition in the finite element model mass matrix or stiffness matrix corresponding to 10% ~ 15% of sensing station place numerical value.Adopt different virtual arrangement schemes just can obtain a series of virtual derivative strucure, the sum of the virtual derivative strucure that can construct
Figure 839804DEST_PATH_IMAGE007
Quantity with sensor
Figure 81430DEST_PATH_IMAGE001
Between the pass be
Figure 904023DEST_PATH_IMAGE005
After virtual derivative strucure is constructed, adopt again the fundamental frequency sensitivity formula to calculate the former rank characteristic frequency parameter that obtains virtual derivative strucure.
The below comes the constructing virtual derivative strucure and calculates corresponding characteristic frequency as an example of the way of virtual arrangement quality example.If source structure not mass matrix and the stiffness matrix branch under the faulted condition is
Figure 281915DEST_PATH_IMAGE008
With
Figure 882661DEST_PATH_IMAGE009
, they can be obtained structural modeling by finite element software, wherein mass matrix
Figure 877161DEST_PATH_IMAGE008
Be diagonal matrix
Figure 221555DEST_PATH_IMAGE010
(1)
In the formula (1)
Figure 719532DEST_PATH_IMAGE011
Expression degree of freedom sum.Source structure occurs before and after the damage, and fair constant and only have loss of rigidity, therefore, the finite element model under the source structure faulted condition is: mass matrix still
Figure 225600DEST_PATH_IMAGE008
, and stiffness matrix is reduced to
Figure 972976DEST_PATH_IMAGE012
,
Figure 855482DEST_PATH_IMAGE012
Be unknown stiffness matrix, it can be expressed as
(2)
Wherein
Figure 884935DEST_PATH_IMAGE014
With
Figure 854028DEST_PATH_IMAGE015
Respectively in the finite element model The impairment parameter of individual unit and element stiffness matrix.Be without loss of generality, suppose on source structure, altogether to have arranged 2 acceleration transducers, be arranged in respectively the 1st and 2 degree of freedom place, can come the virtual arrangement quality to construct 3 virtual derivative strucures in these 2 sensor layout points so, be respectively: first virtual derivative strucure be for separately in lumped mass of the 1st degree of freedom place virtual arrangement
Figure 747214DEST_PATH_IMAGE017
, its stiffness matrix of virtual derivative strucure of gained is still identical with source structure, its mass matrix
Figure 860664DEST_PATH_IMAGE018
Be the source structure mass matrix
Figure 566320DEST_PATH_IMAGE008
With additional virtual mass matrix
Figure 790628DEST_PATH_IMAGE019
Sum, namely
Figure 117705DEST_PATH_IMAGE020
,
Figure 933214DEST_PATH_IMAGE021
(3),(4)
Second virtual derivative strucure is for separately in lumped mass of the 2nd degree of freedom place virtual arrangement
Figure 80161DEST_PATH_IMAGE022
, its stiffness matrix of virtual derivative strucure of gained is still identical with source structure, its mass matrix
Figure 108160DEST_PATH_IMAGE023
Be the source structure mass matrix
Figure 820901DEST_PATH_IMAGE008
With additional virtual mass matrix Sum, namely
Figure 175976DEST_PATH_IMAGE025
,
Figure 742087DEST_PATH_IMAGE026
(5),(6)
The 3rd virtual derivative strucure is simultaneously in 2 lumped masses of the 1st and 2 degree of freedom place virtual arrangement
Figure 309334DEST_PATH_IMAGE017
With
Figure 935488DEST_PATH_IMAGE022
, its stiffness matrix of virtual derivative strucure of gained is still identical with source structure, its mass matrix
Figure 57028DEST_PATH_IMAGE027
Be the source structure mass matrix
Figure 426829DEST_PATH_IMAGE008
With additional virtual mass matrix
Figure 864895DEST_PATH_IMAGE028
Sum, namely
Figure 661950DEST_PATH_IMAGE029
,
Figure 270785DEST_PATH_IMAGE030
(7),(8)
More than after three virtual derivative strucures construct, can calculate corresponding characteristic frequency value by the fundamental frequency sensitivity formula respectively, computing formula is as follows:
Figure 178699DEST_PATH_IMAGE031
Figure 720538DEST_PATH_IMAGE032
, wherein
Figure 688494DEST_PATH_IMAGE033
,
Figure 784626DEST_PATH_IMAGE034
(9)
In the formula (9)
Figure 496230DEST_PATH_IMAGE035
With
Figure 158156DEST_PATH_IMAGE036
For with the mode test macro to source structure measure obtain the
Figure 297013DEST_PATH_IMAGE037
Rank characteristic frequency and the vibration shape,
Figure 614862DEST_PATH_IMAGE038
Be
Figure 130157DEST_PATH_IMAGE039
Of individual virtual derivative strucure
Figure 646589DEST_PATH_IMAGE037
The rank characteristic frequency.Notice
Figure 956348DEST_PATH_IMAGE040
In except arranging that the virtual mass position is not 0 to think that other element all is 0, so do not require the measurement vibration shape that degree of freedom is complete in the calculating of formula (9)
Figure 761493DEST_PATH_IMAGE036
, and as long as part vibration shape value corresponding to the virtual mass position is arranged, and this is the position of sensor just, therefore with formula (9) when calculating used all data all be to be obtained by measurement.
Described Mixed Sensitivity diagnostic routine may further comprise the steps: at first, according to source structure and the virtual derivative strucure finite element model under the faulted condition not, calculate respectively source structure and virtual derivative strucure under the faulted condition not before
Figure 814899DEST_PATH_IMAGE004
Rank characteristic frequency and corresponding fundamental frequency sensitivity; Then before will measuring the source structure of gained
Figure 435105DEST_PATH_IMAGE004
Before rank characteristic frequency and each the virtual derivative strucure
Figure 915765DEST_PATH_IMAGE004
The rank characteristic frequency combines, and lists the one order equation of damage front and back frequency variation, solves all unknown impairment parameter by generalizde inverse
Figure 473785DEST_PATH_IMAGE006
The impairment parameter output that to calculate gained by forms such as histograms at last can accurately identify damage position and the degree of source structure.
The used formula of described Mixed Sensitivity diagnostic routine is as follows: by the source structure model matrix under the faulted condition not
Figure 65304DEST_PATH_IMAGE008
With , can obtain the characteristic frequency of source structure under the faulted condition not by finding the solution following generalized eigenvalue equation
Figure 207889DEST_PATH_IMAGE041
(10)
Wherein
Figure 987626DEST_PATH_IMAGE042
With Respectively that source structure is not under the faulted condition
Figure 728366DEST_PATH_IMAGE037
Individual characteristic frequency and the vibration shape.The
Figure 550829DEST_PATH_IMAGE037
Individual characteristic frequency is about
Figure 817862DEST_PATH_IMAGE016
The one order computing formula of individual impairment parameter is
(11)
In like manner,
Figure 216799DEST_PATH_IMAGE039
Individual virtual derivative strucure not characteristic frequency and the corresponding one order computing formula under the faulted condition is
,
Figure 964493DEST_PATH_IMAGE046
(12),(13)
By equation (10)-(13) just can calculate obtain source structure and all virtual derivative strucures under faulted condition not before
Figure 967084DEST_PATH_IMAGE004
Rank characteristic frequency and corresponding fundamental frequency sensitivity.Then before will measuring the source structure of gained
Figure 772360DEST_PATH_IMAGE004
The rank characteristic frequency and calculated by formula (9) gained each virtual derivative strucure before The rank characteristic frequency combines, and lists the one order equation of all frequency variations of damage front and back, and it is
Figure 178250DEST_PATH_IMAGE047
(14)
Wherein frequency shift is vectorial
Figure 250112DEST_PATH_IMAGE048
For
Figure 628003DEST_PATH_IMAGE049
(15)
The impairment parameter vector
Figure 228749DEST_PATH_IMAGE050
For
Figure 692091DEST_PATH_IMAGE051
(16)
The Mixed Sensitivity matrix
Figure 567643DEST_PATH_IMAGE052
For
Figure 800042DEST_PATH_IMAGE053
(17)
Can find the solution all impairment parameter by equation (14) by generalizde inverse, namely
(18)
Subscript "+" expression is to the Matrix Calculating generalized inverse in the equation (18).Whether will be exported with histogrammic form by all impairment parameter of equation (18) calculating gained at last, and can damage source structure, the position of damage and degree are made accurate identification.
Fig. 2 is a grid structure model to be detected, and this structure basic parameter is: elastic modulus E=200GP a, density p=7.8 * 10 3Kg/m 3, every rod member length L=1m, rod member cross-sectional area A=0.004m 2Simulate two kinds of degree of impairments, the first situation is single damage: suppose rod member 15 losss of rigidity 15% among Fig. 2; The second situation is a plurality of damages: suppose that rod member 7 and 15 rigidity all lose 15% among Fig. 2.
The step that adopts this case institute extracting method that grid structure shown in Figure 2 is damaged identification is as follows: at first, analyze the as can be known sum of unit in its finite element model of this structural model
Figure 522327DEST_PATH_IMAGE055
If quadravalence mode (namely before only measuring
Figure 935991DEST_PATH_IMAGE056
), so according to formula Can determine the quantity of wireless senser
Figure 231023DEST_PATH_IMAGE057
Therefore, layout level and vertical two Wireless Acceleration Sensor on the node A in Fig. 2 are arranged a horizontal acceleration sensor in Node B, altogether arrange three sensors; Secondly, utilize mode test macro (1+2+3) to measure front quadravalence characteristic frequency and the corresponding vibration shape of grid structure vibration; The 3rd, utilize finite element software to set up the not finite element model 5 under the faulted condition of source structure; The 4th, in the level of the node A of source structure and vertical both direction (being the arranged direction of two sensors of node A) virtual arrangement quality respectively, the size of virtual mass all is taken as corresponding to 10% of sensing station place finite element model quality numerical value, thereby obtain two virtual derivative strucures 6, and calculated the front quadravalence characteristic frequency of each virtual derivative strucure by formula (9); The 5th, before measuring the source structure of gained Before rank characteristic frequency and each the virtual derivative strucure
Figure 355154DEST_PATH_IMAGE004
The rank characteristic frequency combines, and inputs in the Mixed Sensitivity diagnostic routine 7, i.e. exportable damage recognition result 8 to source structure.
Adopt the damage recognition result of this case institute extracting method to list among Fig. 3 and Fig. 4 under two kinds of degree of impairments.As seen from Figure 3, for single degree of impairment, this case institute extracting method damage recognition result shows that clearly rod member 15 damages, and the impairment parameter value is 13.9%, and default 15% is very approaching.As seen from Figure 4, for a plurality of degree of impairments, this case institute extracting method damage recognition result shows that equally clearly rod member 7 and 15 damages, and the impairment parameter value is respectively 12.4% and 13.7%, and is all very approaching with default 15%.In sum, the damnification recognition method easy operating that this case is carried is implemented, and recognition result accurately and reliably.
The technical program only utilizes the minority sensor measurement to treat former rank mode of geodesic structure (source structure) vibration, obtain more frequency information by constructing a series of virtual derivative strucure, the frequency of last joint source structure and all derivative strucures identifies the damage status of source structure.In the aforementioned techniques scheme, the technology contents that the present invention does not specify is prior art, repeats no more herein.
The present invention only needs the sensor seldom counted, and front several lower modes of only measuring structural vibration can carry out, and testing apparatus is simple, with low cost; Obtain more frequency parameter by virtual a series of derivative strucure and be used for Damage Assessment Method, effectively overcome the deficiency of existing measuring technology, increased substantially the accuracy of damage identification, be particularly suitable for being applied in the damage identification of large scale structure; Derivative strucure of the present invention is virtual construct, so need not treating actual arrangement quality or rigidity member on the geodesic structure, only need conventional mode testing apparatus and a small amount of sensor to carry out, therefore easy to operate easy to implement, and do not increase extra test job amount; Therefore speed according to power increases rapidly the virtual derivative strucure quantity that the present invention can construct along with the increase of number of sensors, and the used number of sensors of the present invention is minimum under the prerequisite that satisfies damage identification needs.
The above is specific embodiments of the invention only, is not the restriction to this case design, and all equivalent variations of doing according to the design key of this case all fall into the protection domain of this case.

Claims (4)

1. damnification recognition method based on virtual derivative strucure is characterized in that: may further comprise the steps: at first, set up by finite element software and to treat the not finite element model under the faulted condition of geodesic structure (being called again source structure); Secondly, measure former rank mode (being characteristic frequency and the corresponding vibration shape) data that obtain the source structure vibration by the mode test macro, the rank number of mode that needs to measure considers definite according to the complexity of structure, precision and the test site condition of surveying instrument; Then, a series of virtual derivative strucure and calculate corresponding characteristic frequency parameter take source structure as base configuration; At last, the characteristic frequency of comprehensive source structure and all virtual derivative strucures by hybrid frequency sensitivity diagnostic routine, draws the source structure damage recognition result that comprises damage position and degree;
Described mode test macro comprises a data acquisition instrument, a plurality of wireless senser and a cover mode analysis software, is measured former rank characteristic frequency and the corresponding Data of Mode that obtains the source structure vibration by the mode test macro; Wherein wireless senser is arranged in the position for the treatment of to occur easily in the geodesic structure damage, the quantity of wireless senser
Figure 2012105610813100001DEST_PATH_IMAGE001
Can be by formula
Figure 2012105610813100001DEST_PATH_IMAGE002
Determine, wherein
Figure 2012105610813100001DEST_PATH_IMAGE003
Be the sum of unit in the finite element model,
Figure 2012105610813100001DEST_PATH_IMAGE005
Number for the characteristic frequency measured.
2. a kind of damnification recognition method based on virtual derivative strucure as claimed in claim 1, it is characterized in that: described former rank modal data refers to modal data between front 4 ~ 10 rank.
3. a kind of damnification recognition method based on virtual derivative strucure as claimed in claim 1 or 2, it is characterized in that: described virtual derivative strucure can obtain by the way of virtual arrangement quality on source structure or rigidity, the position of virtual mass or rigidity is the position at wireless senser place, can arrange separately also and can make up layout, the size of virtual mass or rigidity generally be taken as source structure not under the faulted condition in the finite element model mass matrix or stiffness matrix corresponding to 10% ~ 15% of sensing station place numerical value; Adopt different virtual arrangement schemes just can obtain a series of virtual derivative strucure, the sum of the virtual derivative strucure that can construct
Figure 2012105610813100001DEST_PATH_IMAGE007
Quantity with sensor
Figure 2012105610813100001DEST_PATH_IMAGE009
Between the pass be
Figure 2012105610813100001DEST_PATH_IMAGE011
After virtual derivative strucure is constructed, adopt again the fundamental frequency sensitivity formula to calculate the former rank characteristic frequency parameter that obtains virtual derivative strucure.
4. a kind of damnification recognition method based on virtual derivative strucure as claimed in claim 1 or 2, it is characterized in that: described Mixed Sensitivity diagnostic routine, may further comprise the steps: at first, according to source structure and the virtual derivative strucure finite element model under the faulted condition not, calculate respectively source structure and virtual derivative strucure under the faulted condition not before
Figure 35991DEST_PATH_IMAGE005
Rank characteristic frequency and corresponding fundamental frequency sensitivity; Then before will measuring the source structure of gained Before rank characteristic frequency and each the virtual derivative strucure
Figure 74671DEST_PATH_IMAGE005
The rank characteristic frequency combines, and lists the one order equation of damage front and back frequency variation, solves all unknown impairment parameter by generalizde inverse
Figure 2012105610813100001DEST_PATH_IMAGE013
The impairment parameter output that to calculate gained by forms such as histograms at last can accurately identify damage position and the degree of source structure.
CN201210561081.3A 2012-12-21 2012-12-21 Damage identification method based on virtual derivative structure Expired - Fee Related CN102998367B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210561081.3A CN102998367B (en) 2012-12-21 2012-12-21 Damage identification method based on virtual derivative structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210561081.3A CN102998367B (en) 2012-12-21 2012-12-21 Damage identification method based on virtual derivative structure

Publications (2)

Publication Number Publication Date
CN102998367A true CN102998367A (en) 2013-03-27
CN102998367B CN102998367B (en) 2015-01-21

Family

ID=47927206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210561081.3A Expired - Fee Related CN102998367B (en) 2012-12-21 2012-12-21 Damage identification method based on virtual derivative structure

Country Status (1)

Country Link
CN (1) CN102998367B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105424797A (en) * 2015-11-05 2016-03-23 北京航空航天大学 Device and method for performing modal testing on inflatable flexible film structure based on hammering excitation method
CN109187769A (en) * 2018-07-27 2019-01-11 南京航空航天大学 A method of defect of pipeline is detected based on supersonic guide-wave quantification
CN111178703A (en) * 2019-12-17 2020-05-19 国网江西省电力有限公司电力科学研究院 Elastic modulus sensitivity analysis method and system for power transmission tower-platform system
CN113297773A (en) * 2021-06-08 2021-08-24 中山大学 Uncertainty quantitative analysis method and device for structural damage identification

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120006118A1 (en) * 2004-07-15 2012-01-12 Stothers Ian Mcgregor Acoustic structural integrity monitoring system and method
CN102539537A (en) * 2011-12-30 2012-07-04 大连理工大学 Frame structure damage method based on additional virtual quality
CN102645488A (en) * 2012-04-28 2012-08-22 绍兴文理学院 Steel roof truss damage diagnosis device
CN102662001A (en) * 2012-04-28 2012-09-12 绍兴文理学院 Steel roof truss damage diagnosis method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120006118A1 (en) * 2004-07-15 2012-01-12 Stothers Ian Mcgregor Acoustic structural integrity monitoring system and method
CN102539537A (en) * 2011-12-30 2012-07-04 大连理工大学 Frame structure damage method based on additional virtual quality
CN102645488A (en) * 2012-04-28 2012-08-22 绍兴文理学院 Steel roof truss damage diagnosis device
CN102662001A (en) * 2012-04-28 2012-09-12 绍兴文理学院 Steel roof truss damage diagnosis method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
PRZEMYSTAW KOLAKOWSKI ET AL.: "The virtual distortion method—a versatile reanalysis tool for structures and systems", 《STRUCTURAL AND MULTIDISCIPLINARY OPTIMIZATION》 *
QINGXIA ZHANG ET AL.: "Simultaneous identification of moving masses and structural damage", 《STRUCTURAL AND MULTIDISCIPLINARY OPTIMIZATION》 *
侯吉林: "约束子结构模型修正方法", 《中国博士学位论文全文数据库•工程科技II辑》 *
侯吉林: "约束子结构模型修正方法", 《中国博士学位论文全文数据库•工程科技II辑》, no. 4, 15 April 2011 (2011-04-15), pages 76 - 90 *
杨秋伟: "基于振动的结构损伤识别方法研究进展", 《振动与冲击》 *
杨秋伟等: "结构损伤识别的附加质量方法", 《工程力学》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105424797A (en) * 2015-11-05 2016-03-23 北京航空航天大学 Device and method for performing modal testing on inflatable flexible film structure based on hammering excitation method
CN109187769A (en) * 2018-07-27 2019-01-11 南京航空航天大学 A method of defect of pipeline is detected based on supersonic guide-wave quantification
CN109187769B (en) * 2018-07-27 2021-02-12 南京航空航天大学 Method for quantitatively detecting pipeline defects based on ultrasonic guided waves
CN111178703A (en) * 2019-12-17 2020-05-19 国网江西省电力有限公司电力科学研究院 Elastic modulus sensitivity analysis method and system for power transmission tower-platform system
CN111178703B (en) * 2019-12-17 2022-12-30 国网江西省电力有限公司电力科学研究院 Elastic modulus sensitivity analysis method and system for power transmission tower-platform system
CN113297773A (en) * 2021-06-08 2021-08-24 中山大学 Uncertainty quantitative analysis method and device for structural damage identification
CN113297773B (en) * 2021-06-08 2023-01-24 中山大学 Uncertainty quantitative analysis method and device for structural damage identification

Also Published As

Publication number Publication date
CN102998367B (en) 2015-01-21

Similar Documents

Publication Publication Date Title
CN101915733B (en) Frequency change-based structural damage degree assessment method
CN104200005A (en) Bridge damage identification method based on neural network
CN104200265A (en) Improved bridge damage identification method based on neural network
CN108280294B (en) A kind of cable arch structure damage combined recognising method based on modal parameter
CN109839441A (en) A kind of bridge Modal Parameters Identification
Ni et al. Dynamic property evaluation of a long-span cable-stayed bridge (Sutong bridge) by a Bayesian method
Gillich et al. How to correlate vibration measurements with FEM results to locate damages in beams
Yang et al. Damage identification of frame structure based on approximate Metropolis–Hastings algorithm and probability density evolution method
CN104133960A (en) Improved optimal arranging method of static sensors
CN102998367B (en) Damage identification method based on virtual derivative structure
CN107271127A (en) Based on the operational modal parameter recognition methods and device extracted from iteration pivot
CN103913512B (en) The damage reason location system of suspension cable periodic detection
Mordini et al. The finite element model updating: a powerful tool for structural health monitoring
CN110008520B (en) Structural damage identification method based on displacement response covariance parameters and Bayesian fusion
An et al. Real-time fast damage detection of shear structures with random base excitation
Żółtowski et al. The use of modal analysis in the evaluation of welded steel structures.
Kim et al. Non-baseline damage detection based on the deviation of displacement mode shape data
Schommer et al. Health monitoring based on dynamic flexibility matrix: theoretical models versus in-situ tests
CN107246944A (en) A kind of Structural Damage Identification theoretical based on statistical moment
CN106383003A (en) Cable structure cable force measurement method and system based on flexibility identification
Akintunde et al. Singular value decomposition and unsupervised machine learning for virtual strain sensing: Application to an operational railway bridge
CN109781863A (en) Structure erosion two-stage detection method and its system based on fast vibration test
Salehi et al. A structural damage detection technique based on measured frequency response functions
Guan Vibration-based structural health monitoring of highway bridges
CN101793590B (en) Structural impact damage diagnostic method based on blackboard cooperation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150121

Termination date: 20161221

CF01 Termination of patent right due to non-payment of annual fee