CN102323089A - Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement - Google Patents
Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement Download PDFInfo
- Publication number
- CN102323089A CN102323089A CN201110143101A CN201110143101A CN102323089A CN 102323089 A CN102323089 A CN 102323089A CN 201110143101 A CN201110143101 A CN 201110143101A CN 201110143101 A CN201110143101 A CN 201110143101A CN 102323089 A CN102323089 A CN 102323089A
- Authority
- CN
- China
- Prior art keywords
- vector
- current
- cable structure
- rope
- time
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention provides a progressive method for identifying a slack cable on the basis of angle monitoring during support generalized displacement. According to angle monitoring, taking into consideration the fact that linear relationships between the current value vector of the monitored amount and an initial value vector of the monitored amount, a virtual unit damage monitored amount variation matrix and the current nominal virtual damage vector are approximate with each other through monitoring the support generalized coordinate of the structure during identifying the loose supporting cable, the invention provides a method for approaching nonlinear relation in a sectional way by using the linear relationships in order to overcome the defects. By using the method, the virtual damaged cable can be identified; after the real damaged cable is distinguished by using a nondestructive testing method and other methods, the rest virtual damaged cables are loose supporting cables; and then the cable length needing to be adjusted can be determined according to the relationship between the looseness degree and the virtual damage degree.
Description
Technical field
(for example the bearing generalized displacement refers to that bearing is along the angular displacement around X, Y, Z axle of the displacement of the lines of X, Y, Z axle and bearing when the bearing generalized displacement; Corresponding to the bearing generalized displacement; The bearing generalized coordinate refers to that bearing is about the coordinate of X, Y, Z axle and the bearing angular coordinate about X, Y, Z axle); Rope supporting structure (particularly large-scale Cable Structure is discerned in the monitoring that the present invention is based on angle equivalent; For example large-scale cable-stayed bridge, suspension bridge) cable system (refer to all supporting ropes) in the supporting rope of need adjustment Suo Li, and provide the long adjustment amount of concrete rope, genus engineering structure security fields.
Background technology
Cable system is Cable Structure (particularly large-scale Cable Structure normally; For example large-scale cable-stayed bridge, suspension bridge) key components; Owing to reason such as lax, new construction is completed, and the Suo Li of supporting rope can change the lax variation that also can cause the supporting cable force of its supporting rope behind the structure long service usually after a period of time; These change the variation that all will cause structural internal force; Safety to structure causes harmful effect, will cause the inefficacy of structure when serious, and therefore accurately and timely discerning needs the supporting rope of adjustment Suo Li to be very important.
The health status of supporting cable system changes after (for example take place lax, damage etc.); Can cause the variation of the measurable parameter of structure; For example caused variation (the for example arbitrarily variation of the angle coordinate of the straight line of any this point of mistake in any section of body structure surface of angle coordinate of any imaginary line of the every bit of Cable Structure; The perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily); In fact the variation of structure angle has comprised the health status information of cable system; That is to say the health status that to utilize the structural angle degrees of data to judge structure; Can discern damaged cable based on angle monitor (the present invention is called monitored angle-data " monitored amount ", and the back is mentioned " monitored amount " and just is meant monitored angle-data), monitored amount is except the influence that receives the cable system health status; Also can receive the influence of Cable Structure bearing generalized displacement (usually can take place), also not have a kind of disclosed, effective health monitoring systems and method to solve this problem at present.Therefore can monitor based on monitored amount and discern the rope that needs adjustment Suo Li; Like this when the bearing generalized displacement; A method of can rational and effective setting up the relation of (specifically the characteristic parameter according to rope characterizes the rope that needs adjustment Suo Li) between the characteristic parameter of monitored amount with all ropes just must be arranged, and the recognition result of the supporting rope of the need adjustment Suo Li that sets up based on this method just can be more credible.
Summary of the invention
Technical matters:The objective of the invention is when the Cable Structure bearing has generalized displacement; To identification problem in the cable system in the Cable Structure, that need the supporting rope of adjustment Suo Li, a kind of structure health monitoring method monitoring, that can discern the supporting rope that needs adjustment Suo Li rationally and effectively based on angle equivalent is disclosed.
Reason according to the Suo Li that supports rope changes can change the three kinds of situation that be divided into the Suo Li of supporting rope: the one, and the supporting rope has received damage, for example supports rope and localized cracks and corrosion or the like have occurred; The 2nd, supporting rope and not damaged; But variation has also taken place in Suo Li; The one of the main reasons that this variation occurs is that variation has taken place the Suo Changdu (be called drift, the present invention specially refers to support the drift of that section rope between rope two supporting end points) under the supporting rope free state (this moment, Suo Zhangli claimed that also Suo Li is 0); The 3rd, supporting rope and not damaged, but the Cable Structure bearing has had rotation, also can cause the variation of structural internal force, also will cause the variation of Suo Li certainly.One of main purpose of the present invention is exactly when the bearing generalized displacement, identify drift the support cable that changes has taken place, and identify the change amount of their drift, and this change amount provides direct foundation for the Suo Li adjustment of this rope.The reason that supporting rope drift changes is not single, and for ease, the present invention is referred to as slack line with the supporting rope that drift changes.
Technical scheme:The present invention is made up of the two large divisions.Be respectively: one, set up required knowledge base of the health monitoring systems be used for discerning supporting rope cable system, that need adjustment Suo Li and parameter method, based on knowledge base (containing parameter), based on the generalized displacement of actual measurement Cable Structure bearing, based on the method for the supporting rope of need adjustment Suo Li monitoring, the identification Cable Structure of monitored amount equivalent; Two, the software and hardware part of health monitoring systems.
First of the present invention: set up required knowledge base of the health monitoring systems be used for discerning supporting rope cable system, that need adjustment Suo Li and parameter method, based on knowledge base (containing parameter), based on the generalized displacement of actual measurement Cable Structure bearing, based on the method for the supporting rope of need adjustment Suo Li monitoring, the identification Cable Structure of monitored amount equivalent.Can be successively circularly as follows, laddering carrying out, to obtain the health status assessment of cable system more accurately.
The first step: during circulation beginning each time, the cable system initial virtual damage vector in the time of at first need setting up or set up this circulation beginning
d o i (
i=1,2,3 ...Because in fact the supporting rope possibly be lax and not damage is the expression difference, claim here " virtual lesion ", the back with), set up the initial Mechanics Calculation benchmark model A of Cable Structure
o(for example finite element benchmark model, A in the present invention
oBe constant), set up the current Mechanics Calculation benchmark model A of Cable Structure
Ti o(finite element benchmark model for example, A in circulation each time
Ti oBring in constant renewal in), set up the Mechanics Calculation benchmark model A of Cable Structure
i(finite element benchmark model for example,
i=1,2,3 ...).Letter i is except the place of representing number of steps significantly, and alphabetical in the present invention i only representes cycle index, i.e. the i time circulation.
If it is total in the cable system
NThe root rope, " initial virtual damage vector is designated as the cable system that needs during the i time circulation beginning
d o i " (shown in (1)), use
d o i Cable Structure is (with Mechanics Calculation benchmark model A when representing this time circulation beginning
iThe health status of cable system expression).
In the formula (1)
d i Oj (
i=1,2,3,
; j=1,2,3 ....,
N) expression circulates for the i time when beginning, Mechanics Calculation benchmark model A
iIn cable system
jThe initial virtual impairment value of root rope,
d i Oj Be to represent at 0 o'clock
jRoot rope not damaged does not have lax, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes in the time of between 0 and 100%
jThe load-bearing capacity of root rope forfeiture corresponding proportion.In the formula (1)
TThe transposition of expression vector (back together).
Setting up cable system initial virtual damage vector during circulation beginning for the first time (is designated as according to formula (1)
d 1 o ) time, the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up cable system initial virtual damage vector
d 1 o If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vector
d 1 o Each element numerical value get 0.
The i time (
i=2,3,4,5,6 ...) when beginning circulation the cable system initial virtual damage vector that needs
d i o , be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
The Mechanics Calculation benchmark model that need set up during the i time circulation beginning or the Mechanics Calculation benchmark model of having set up are designated as A
i
Design drawing, as-constructed drawing and the measured data of the Cable Structure in being completed according to Cable Structure (comprises that the Non-Destructive Testing data etc. of rope can express measured datas such as the data of the health status of rope, Cable Structure shape data, structural angle degrees of data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge), utilize mechanics method (for example finite element method) to set up A
oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, structure is surveyed; The measured data that obtains Cable Structure (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope); According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method (for example finite element method) to set up A
oNo matter which kind of method to obtain A with
o, based on A
oThe Cable Structure computational data that calculates (to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge) must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this
oStrain computational data, Suo Li computational data, Cable Structure shape computational data and displacement computational data, Cable Structure angle-data etc. under the analog case of calculating gained, the measured data when truly taking place near institute's analog case reliably.Corresponding to A
oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector
U oA
oWith
U oBe constant, only when circulation beginning for the first time, set up.
The Mechanics Calculation benchmark model of the Cable Structure of setting up during circulation beginning for the first time is designated as A
1, A
1Just equal A
oA
1The health status of corresponding rope by
d 1 o Describe.
The i time (
i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins
i, be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Existing Mechanics Calculation benchmark model A
1With cable system initial virtual damage vector
d 1 o After, model A
1In the virtual lesion of each rope by vector
d 1 o Express.At A
1The basis on, the virtual lesion value of all ropes is changed to 0, mechanical model A
1The virtual lesion that is updated to all ropes all is that 0 mechanical model (is designated as A
0), mechanical model A
0Be actually the excellent corresponding mechanical model of lax Cable Structure that do not have.Might as well claim model A
0For the not damaged of Cable Structure does not have relaxation model A
0
" the whole monitored angle-data of structure " is by on the structure
KIndividual specified point, cross each specified point
LIndividual appointment straight line, each specifies straight line
HIndividual angle coordinate component is described, and the variation of structure angle is exactly the variation of the angle coordinate component of all appointments all specified points, all appointment straight lines.Each total
M (M=K * L * H)Individual angle coordinate component measurement value or calculated value characterize the angle information of structure.
KWith
MMust not be less than the quantity of supporting rope
NFor simplicity, in the present invention " the monitored angle-data of structure " abbreviated as " monitored amount ".
The present invention is with " the initial value vector of monitored amount
C i o " (
i=1,2,3 ...) expression the i time (
i=1,2,3,4,5,6 ...) initial value (referring to formula (2)) of the monitored amount of all appointments when circulation begins,
C i o Full name be " the initial value vector of the i time monitored amount of circulation ".
In the formula (2)
C i Ok (
i=1,2,3,
k=1,2,3,
., M; M>=N;) be in the i time when beginning circulation, the Cable Structure the
kIndividual monitored amount.Vector
C i o Define by the front
MIndividual monitored amount forms according to certain series arrangement, and this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
During circulation beginning for the first time, " the initial value vector of the 1st the monitored amount that circulates
C 1 o " (seeing formula (2)) be made up of measured data, because according to model A
1The initial value of calculating the monitored amount of gained approaches corresponding measured value reliably, in the narration of back, will represent this calculated value composition of vector and measured value composition of vector with prosign.
The i time (
i=2,3,4,5,6 ...) when beginning circulation need " initial value of the i time monitored amount of circulation is vectorial
C i o ", be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Second step: in Cable Structure military service process, in circulation each time, constantly (all data are formed current cable structure actual measurement bearing generalized coordinate vector to actual measurement acquisition Cable Structure bearing generalized coordinate current data
U Ti , vector
U Ti Definition mode with the vector
U oIdentical).For simplicity, for the i time circulation, the Cable Structure bearing generalized coordinate current data the when last time is upgraded current Mechanics Calculation benchmark model is designated as current cable structural bearings generalized coordinate vector
U Ti oSet up and renewal A
Ti oMethod be: in round-robin each time zero hour, the current Mechanics Calculation benchmark model A of Cable Structure
Ti oJust equal A
i(
i=1,2,3,4,5,6 ...).In Cable Structure military service process, constantly actual measurement obtains Cable Structure bearing generalized coordinate data and obtains current cable structure actual measurement bearing generalized coordinate vector
U Ti If,
U Ti Equal
U Ti o, then need be to A
Ti oUpgrade; If
U Ti Be not equal to
U Ti o, then need be to A
Ti oUpgrade, at this moment
U Ti With
U oDifference be exactly the Cable Structure bearing about initial position (corresponding to A
o) the bearing generalized displacement (with bearing generalized displacement vector
VThe generalized displacement of expression bearing).Upgrade A
Ti oMethod be: at A
oThe basis on make that the health status of rope is cable system initial damage vector
d i o , more further to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint, the numerical value of current bearing generalized displacement constraint is just taken from current bearing generalized displacement vector
VThe numerical value of middle corresponding element is to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
Ti o, upgrade A
Ti oAfter,
U Ti oAll elements numerical value is used
U Ti All elements numerical value replaces, and has promptly upgraded
U Ti o, so just obtained correctly corresponding to A
Ti o U Ti o
The 3rd step: circulation each time needs to set up " the monitored numerical quantity transformation matrices of virtual unit damage " and " nominal virtual unit damage vector ",
," the monitored numerical quantity transformation matrices of virtual unit damage " that the i time circulation set up is designated as
Δ C i (
i=1,2,3 ...)." nominal virtual unit damage vector " that the i time circulation set up is designated as
D i u In circulation each time
Δ C i With
D i u Bring in constant renewal in, promptly upgrading current Mechanics Calculation benchmark model A
Ti oThe time, upgrade the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i With nominal virtual unit damage vector
D i u
Set up and the monitored numerical quantity transformation matrices of renewal virtual unit damage
Δ C i With nominal virtual unit damage vector
D i u Process following:
Current Mechanics Calculation benchmark model A in Cable Structure
Ti oThe basis on carry out several times and calculate, equal the quantity of all ropes on the calculation times numerical value.Calculating each time in the hypothesis cable system has only a rope on the basis of original virtual lesion (original virtual lesion can be 0, can not be 0 also), to increase virtual unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is virtual unit damage) again.Calculate for convenient; When setting virtual unit damage in the circulation each time can all be structural health conditions during this time circulation beginning as being healthy fully, and set on this basis virtual unit damage (in subsequent step, damage numerical value that calculate, rope---be called nominal virtual lesion
d i c (
i=1,2,3 ...), all with respect to this time when beginning circulation, with the health status of rope as being healthy fully speech, therefore must according to after the formula that provides of the literary composition nominal virtual lesion that will calculate be converted into actual and virtual and damage).The rope that occurs virtual lesion during a round-robin calculates each time together is different from the rope of the virtual lesion that occurs in other time calculating; And supposition each time has the virtual unit damage value of the rope of virtual lesion can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
D i u " (shown in (3)) write down the virtual unit damage of the supposition of all ropes in each time circulation, circulation time is designated as for the first time
D 1 u , calculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, appointment in front
MThe current calculated value of individual monitored amount calculates gained each time
MThe current calculated value of individual monitored amount is formed one " the current numerical value vector of the calculating of monitored amount " (when hypothesis the
jWhen the root rope had unit damage, available formula (4) was represented all appointments
MThe current numerical value vector of the calculating of individual monitored amount
C 1 Tj ); The current numerical value vector of the calculating of the monitored amount that calculates each time deducts the initial value vector of monitored amount
C 1 o , the gained vector is exactly that " the numerical value change vector of monitored amount " of (is mark with the position of rope that virtual unit damage is arranged or numbering etc.) is (when the under this condition
jWhen the root rope has virtual unit damage, use
δ C 1 j The numerical value change vector of representing monitored amount,
δ C 1 j Definition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vector for formula (4)
D 1 u jIndividual element
D Uj Gained), the numerical value change of monitored amount vector
δ C 1 j Each element representation since when calculating supposition the Na Gensuo (for example the of virtual unit damage is arranged
jThe root rope) virtual unit damage (for example
D Uj ), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the virtual unit damage of supposition
D Uj Rate of change; Have
NThe root rope just has
NIndividual " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has
M(it is general,
M>=N) individual element, by this
NIndividual " the numerical value change vector of monitored amount " formed successively to be had
M * N" the unit damage monitored quantitative change matrix of individual element
Δ C 1 "
(MOK
NRow
), each vector
δ C 1 j (
j=1,2,3 ....,
N) be matrix
Δ C 1 One row,
Δ C 1 Definition suc as formula shown in (8).
Nominal virtual unit damage vector in the formula (3)
D i u Element
D i Uj (
i=1,2,3,
j=1,2,3 ....,
N) expression the
iSuppose in the inferior circulation
jThe virtual unit damage numerical value of root rope, vector
D i u In the numerical value of each element can be the same or different.
Element in the formula (4)
C i Tjk (
i=1,2,3,
j=1,2,3 ....,
N; k=1,2,3 ....,
M; M>=N) expression the
iInferior circulation is because the
jWhen the root rope has virtual unit damage, according to coding rule pairing
kThe current numerical value of the calculating of the monitored amount of individual appointment.
The subscript of each amount in the formula (5)
i(
i=1,2,3 ...) expression the
iInferior circulation, subscript
j(
j=1,2,3 ....,
N) expression the
jThe root rope has virtual unit damage, in the formula
D i Uj It is vector
D i u In
jIndividual element.Vector
δ C i j Definition suc as formula shown in (6),
δ C i j k(
k=1,2,3 ....,
M; M>=N) individual element
δ C i Jk Expression the
iIn the inferior circulation, set up matrix
Δ C i The time, suppose
jWhen having virtual unit damage, the root rope calculates gained the
kThe change amount of individual monitored amount is with respect to the virtual unit damage of supposition
D i Uj Rate of change, it defines suc as formula shown in (7).
(6)
(7)
The definition of each amount was narrated in front in the formula (7).
Vector in the formula (8)
δ C i j (
i=1,2,3 ....,,
j=1,2,3 ....,
N) expression the
iIn the inferior circulation, because the
jThe root rope has virtual unit damage
D i Uj Cause, the relative value of all monitored amounts changes.Matrix
Δ C i Row (subscript
j) coding rule and front vector
d i o The subscript of element
jCoding rule identical.
In Cable Structure military service process, in circulation each time, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, in case monitor
U Ti Be not equal to
U Ti o, then need get back to for second step to A
Ti oUpgrade, to A
Ti oIt is right to get into this step after upgrading again
Δ C i Upgrade.In fact in circulation each time
Δ C i Bring in constant renewal in, promptly upgrading current Mechanics Calculation benchmark model A
Ti oAfterwards, upgrade the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i
The 4th step: the current health status of identification cable system.Detailed process is following.
The
i(
i=1,2,3 ...) in the inferior circulation, cable system " current (calculating or actual measurement) numerical value vector of monitored amount
C i " " initial value of monitored amount is vectorial together
C i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between linear approximate relationship, shown in (9) or formula (10).
Current (calculating or actual measurement) numerical value vector of monitored amount in formula (9) and the formula (10)
C i Definition be similar to the initial value vector of monitored amount
C i o Definition, see formula (11); Cable system " current nominal virtual lesion vector
d i c " definition see formula (12).
Element in the formula (11)
C i k (
i=1,2,3 ....;
k=1,2,3 ....,
M; M>=N) be
iInferior circulation time Cable Structure, be numbered according to coding rule is pairing
kThe current numerical value of monitored amount.
In the formula (12)
d i Cj (
i=1,2,3 ....;
j=1,2,3 ....,
N) be
iCable system in the inferior circulation
jThe current nominal virtual lesion value of root rope, vector
d i c The subscript of element
jCoding rule and matrix
Δ C i The coding rule of row identical.
When the rope actual damage was not too big, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure was also less, and the represented a kind of like this linear relationship of formula (9) or formula (10) is less with the error of actual conditions, and error can be used error vector
e i (formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).
(13)
In the formula (13)
Abs ()Be the function that takes absolute value, each element of the vector of trying to achieve in the bracket is taken absolute value.
Because there are certain error in formula (9) or the represented linear relationship of formula (10), therefore can not be simply according to formula (9) or formula (10) and " current (actual measurement) numerical value vector of monitored amount
C i " come directly to find the solution and obtain " current nominal virtual lesion vector
d i c ".If done like this, the vector that obtains
d i c In element in addition bigger negative value can appear, just negative damage, this obviously is irrational.Therefore obtain vector
d i c Acceptable separating (promptly have reasonable error, but can confirm the position and the virtual lesion degree thereof of virtual damaged cable more exactly) become a rational solution, available formula (14) is expressed this method.
In the formula (14)
Abs ()Be the function that takes absolute value, vector
g i Description departs from ideal linearity relation (formula (9) or formula (10))
Reasonable deviation, define by formula (15).
In the formula (15)
g i k (
i=1,2,3 ....;
k=1,2,3 ....,
M) described
iThe maximum allowable offset that departs from the ideal linearity relation shown in formula (9) or the formula (10) in the inferior circulation.Vector
g i Can be according to the error vector of formula (13) definition
e i Tentative calculation is selected.
At " the initial value vector of monitored amount
C i o " (survey or calculate), " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " (calculating) and " the current numerical value vector of monitored amount
C i " when (actual measurement obtains) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain " current nominal virtual lesion vector
d i c " acceptable separating, " current actual virtual lesion vector then
d i " element of (formula (16) is seen in definition) can calculate according to formula (17), just obtained " " current actual virtual lesion vector
d i Thereby, can by
d i Confirm the position and the virtual lesion degree of virtual damaged cable,, need just to have confirmed rope and the long adjustment amount of rope thereof of adjustment Suo Li then according to position and the relax level of below the method for narration being confirmed slack line.
In the formula (16)
d i j (
i=1,2,3,
; j=1,2,3 ....,
N) expression the
iIn the inferior circulation
jThe actual virtual lesion value of root rope, formula (17) is seen in its definition,
d i j Be to represent at 0 o'clock
jRoot rope not damaged does not have lax, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes in the time of between 0 and 100%
jThe load-bearing capacity of root rope forfeiture corresponding proportion, vector
d i Coding rule and the formula (1) of element in vector
d i o The coding rule of element identical.
In the formula (17)
d i Oj (
i=1,2,3,4,
; j=1,2,3 ....,
N) be vector
d i o jIndividual element,
d i Cj It is vector
d i c jIndividual element.
Narration has obtained the actual virtual lesion vector of Suo Dangqian below
d i After, how to confirm the position and the relax level of slack line.
If it is total in the cable system
NRoot supporting rope, structure rope force data by
NThe Suo Li of root supporting rope describes.Available " initial rope force vector
F o " represent that all support the initial Suo Li (formula (18) is seen in definition) of ropes in the Cable Structure.Because the initial Suo Li that calculates gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial Suo Li reliably, in the narration of back, will represent this calculated value and measured value with prosign.
In the formula (18)
F Oi (
j=1,2,3,
., N) be in the Cable Structure
jThe initial Suo Li of root supporting rope, this element is according to the Suo Li of coding rule corresponding to appointment supporting rope.Vector
F o Be constant, irrelevant with cycle index, after when circulation beginning for the first time, confirming, just no longer change.Setting up the Mechanics Calculation benchmark model A of Cable Structure
1The time used vector
F o
Use " current cable force vector among the present invention
F i " expression the
iThe current cable power (formula (19) is seen in definition) of all supporting ropes in the Cable Structure that inferior circulation time actual measurement obtains.
In the formula (19)
F i j (
i=1,2,3,4,
; j=1,2,3,
., N) be
iIn the inferior circulation time Cable Structure
jThe current cable power of root supporting rope.
Among the present invention, under supporting rope original state (not damaged, do not have lax), and the supporting rope is when being in free state (free state refers to that Suo Li is 0, back with), and the length of supporting rope is called initial drift, with " initial drift vector
l o " represent that all support the initial drift (formula (20) is seen in definition) of ropes in the Cable Structure.
In the formula (20)
l Oj (
j=1,2,3,
., N) be in the Cable Structure
jThe initial drift of root supporting rope.Vector
l o Be constant, irrelevant with cycle index, after when circulation beginning for the first time, confirming, just no longer change.
Among the present invention, with " current drift vector
l i " expression the
iThe current drift (formula (21) is seen in definition) of all supporting ropes in the inferior circulation time Cable Structure.
In the formula (21)
l i j (
i=1,2,3,4,
; j=1,2,3,
., N) be
iIn the inferior circulation time Cable Structure
jThe current drift of root supporting rope.
Among the present invention, with " drift changes vectorial Δ
l i " (or claiming supporting Suo Dangqian relax level vector) expression the
iThe change amount (formula (22) and formula (23) are seen in definition) of the drift of all supporting ropes in the inferior circulation time Cable Structure.
Δ in the formula (22)
l i j (
i=1,2,3,4,
; j=1,2,3,
., N) be current (
iInferior circulation time) in the Cable Structure
jThe change amount of the drift of root supporting rope, formula (23), Δ are seen in its definition
l i j Be not that 0 rope is a slack line, Δ
l i j Numerical value be the slack of rope, and expression cable system the
jThe current relax level of root supporting rope, the long adjustment amount of rope of this rope when also being adjustment Suo Li.
Through slack line is carried out the relax level identification that slack line is carried out in the mechanics equivalence with damaged cable, the mechanical condition of equivalence is in the present invention:
The mechanics parameters of initial drift, geometrical property parameter and material when one, the nothing of the rope of two equivalences relaxes with not damaged is identical;
Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.
When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing slack line with the damaged cable of equivalence, Cable Structure any variation can not take place, vice versa.
Among the present invention, the
iInferior circulation time is with the
j(its current relax level is used Δ to individual supporting rope
l i j The current actual virtual lesion degree of definition) carrying out the virtual impaired supporting rope of equivalence is used
d i j Expression (
d i j Definition see formula (16) and formula (17)).Lax the
jThe current relax level Δ of individual supporting rope
l i j (Δ
l i j Definition see formula (22)) with the current actual virtual lesion degree of damaged cable of equivalence
d i j Between relation confirm by aforementioned two mechanics equivalent conditions.Δ
l i j With
d i j Between physical relationship can adopt accomplished in many ways, for example can be directly confirm (referring to formula (24)) according to aforementioned equivalent condition, also can adopt based on the Ernst equivalent elastic modulus to replace in the formula (24)
ERevise the back and confirm (referring to formula (25)), also can adopt and confirm based on other methods such as trial and error procedure of finite element method.
In formula (24) and the formula (25)
EBe the elastic modulus of this supporting rope,
ABe the cross-sectional area of this supporting rope,
F i j Be the current cable power of this supporting rope,
d i j Be the current actual virtual lesion degree of this supporting rope,
ω i Be the weight of the unit length of this supporting rope,
l i Jx It is the horizontal range of two supporting end points of this supporting rope.Item in the formula (25) in [] is the Ernst equivalent elastic modulus of this supporting rope, can just can confirm to support Suo Dangqian relax level vector Δ by formula (24) or formula (25)
l i Formula (25) is the correction to formula (24).
The 5th step: judge whether to finish this (the
iInferior) circulation, if then accomplish the tailing in work before this loop ends, for next time (promptly the
i+ 1 time,
i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following.
This (
iInferior) try to achieve " current nominal virtual lesion vector in the circulation
d i c " after, at first, set up " mark vector according to formula (26)
B i ", formula (27) has provided " mark vector
B i "
jThe definition of individual element; If " mark vector
B i " element be 0 entirely, then in this circulation, continue health monitoring and calculating to cable system; If " mark vector
B i " element be not 0 entirely, then accomplish subsequent step after, get into circulation next time.So-called subsequent step is: at first, calculate next time (promptly according to formula (28)
i+ 1 time,
i=1,2,3,4 ...) required " the initial virtual damage vector of circulation
d i+
1 o " each element
d i+
1 Oj The second, at Mechanics Calculation benchmark model A
i(
i=1,2,3,4 ...) or the not damaged model A of Cable Structure
0The basis on, make the health status of rope do
d i+
1 o The back upgrade and to obtain next time (the i+1 time,
i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation
I+1At last, through to Mechanics Calculation benchmark model A
I+1The initial value that calculates monitored amount, by its form next time (promptly the i+1 time,
i=1,2,3,4 ...) required " the initial value vector of monitored amount of circulation
C I+1 o " (
i=1,2,3,4 ...).
Mark vector in the formula (26)
B i Subscript
iExpression the
iInferior circulation, its element
B i j (
j=1,2,3 ..., subscript N)
jExpression the
jThe damage characteristic of root rope can only be got 0 and 1 two amount, and concrete value rule is seen formula (27).
Element in the formula (27)
B i j It is " mark vector
B i "
jIndividual element,
D i Uj It is " nominal virtual unit damage vector
D i u "
jIndividual element (seeing formula (3)),
d i Cj It is " current nominal virtual lesion vector
d i c "
jIndividual element (seeing formula (12)), they all represent
jThe relevant information of root rope.
In the formula (28)
D i Uj It is " nominal virtual unit damage vector
D i u "
jIndividual element (seeing formula (3)),
d i Cj It is " current nominal virtual lesion vector
d i c "
jIndividual element (seeing formula (12)).
Second portion of the present invention: the software and hardware part of health monitoring systems.Hardware components comprises monitoring system (monitoring the horizontal range of monitored amount, the generalized displacement of monitoring Cable Structure bearing, monitoring Suo Li, monitoring supporting rope two supporting end points), signal picker and computing machine etc.Requirement is monitored each monitored amount in real time or quasi real time, monitors the Suo Li of each support cable, is monitored the horizontal range that each support cable two supports end points.Software should be used following function by tool: software section should be accomplished the process that first of the present invention sets, and promptly accomplishes functions such as needed among the present invention, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.
The inventive method specifically comprises:
A. establish total N root rope, at first confirm the coding rule of rope, with rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule;
B. confirm the measured point of appointment, give all specified point numberings; Confirm the measured straight line of each measurement point, gave the measured straight line numbering of all appointments; Confirm the measured angle coordinate component of each measured straight line, give all measured angle coordinate component numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored angle-data of structure " is made up of above-mentioned all measured angle coordinate components; For simplicity, in the present invention " the monitored angle-data of structure " abbreviated as " monitored amount "; The quantity of measurement point must not be less than the quantity of rope; The quantity sum of all measured angle coordinate components must not be less than the quantity of rope;
C. the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up initial virtual damage vector
d i o , wherein i representes cycle index, back i and subscript i represent cycle index, and i=1,2,3, Circulation time for the first time
d i o Be designated as
d 1 o If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think perhaps that the structure original state is lax for not having, during the not damaged state, vector
d 1 o Each element numerical value get 0.
D. setting up initial virtual damage vector
d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount
C 1 o
E. setting up initial virtual damage vector
d 1 o Initial value vector with monitored amount
C 1 o The time, directly measure the initial Suo Li that calculates all supporting ropes, form initial rope force vector
F o Simultaneously, obtain the initial drift that all support ropes, form initial drift vector according to structural design data, completion data
l o Simultaneously, obtain the initial geometric data of Cable Structure according to structural design data, completion data or actual measurement; Simultaneously, actual measurement or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information;
F. set up the initial Mechanics Calculation benchmark model A of Cable Structure
o, set up initial Cable Structure bearing generalized coordinate vector
U o, the Mechanics Calculation benchmark model A of the Cable Structure that needs when setting up circulation beginning for the first time
1The measured data of the Cable Structure in Cable Structure completion; This measured data comprises measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data, all ropes; And the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope; According to design drawing and as-constructed drawing, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure
oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, this Cable Structure is surveyed; Obtain the measured data of Cable Structure equally; According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure equally
oNo matter which kind of method to obtain A with
o, based on A
oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A
oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector
U oA
oWith
U oBe constant, only when circulation beginning for the first time, set up; The Mechanics Calculation benchmark model of the Cable Structure of setting up during the i time circulation beginning is designated as A
i, wherein i representes cycle index; Alphabetical i is except the place of representing number of steps significantly in the application form of the present invention, and alphabetical i only representes cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A
1, A among the present invention
1Just equal A
oFor narrating conveniently name " the current Mechanics Calculation benchmark model of Cable Structure A
Ti o", A in circulation each time
Ti oCan bring in constant renewal in as required, during circulation beginning each time, A
Ti oEqual A
iIt is convenient to be similarly narration, name " Cable Structure actual measurement bearing generalized coordinate vector
U Ti ", in circulation each time, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector
U Ti , vector
U Ti Element with the vector
U oThe generalized coordinate of the equidirectional of the element representation same abutment of same position; For the purpose of narrating conveniently,, the last time is upgraded A for the i time circulation
Ti oThe time Cable Structure bearing generalized coordinate current data be designated as current cable structural bearings generalized coordinate vector
U Ti oDuring circulation beginning for the first time, A
T1 oEqual A
1,
U T1 oEqual
U oA
1The health status of corresponding rope by
d 1 o Describe; Mechanics Calculation benchmark model A
iThe health status of corresponding rope by
d i o Describe; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;
When g. circulation begins each time, make A
Ti oEqual A
iActual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector
U Ti , according to current cable structure actual measurement bearing generalized coordinate vector
U Ti , upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
Ti oWith current cable structural bearings generalized coordinate vector
U Ti o
H. at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining the monitored numerical quantity transformation matrices of Cable Structure virtual unit damage
Δ C i With nominal virtual unit damage vector
D i u
I. actual measurement obtains the current cable power of all supporting ropes of Cable Structure, forms the current cable force vector
F i Simultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount
C i ".Actual measurement calculates the space coordinates of two supporting end points of all support cables, and the space coordinates of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.When numbering to the element of the institute's directed quantity that occurred before this step and this step; Should use same coding rule; Can guarantee like this before this step and this step with each vector that occurs afterwards, number identical element, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
J. define current nominal virtual lesion vector to be asked
d i c With current actual virtual lesion vector
d i The damage vector
d i o ,
d i c With
d i Element number equal the quantity of rope, be one-to-one relationship between the element of damage vector and the rope, the element numerical value of damage vector is represented the virtual lesion degree or the health status of corresponding rope;
K. according to " the current numerical value vector of monitored amount
C i " " initial value of monitored amount is vectorial together
C i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes in the formula 1
d i c Other outer amount is known, finds the solution formula 1 and just can calculate current nominal virtual lesion vector
d i c
L. the current actual virtual lesion vector that utilizes formula 2 to express
d i Element
d i j With initial virtual damage vector
d i o Element
d i Oj With current nominal virtual lesion vector
d i c Element
d i Cj Between relation, calculate current actual virtual lesion vector
d i All elements.
In the formula 2
j=1,2,3 ..., N.
Because current actual virtual lesion vector
d i Element numerical value represent the current actual virtual lesion degree of corresponding rope, promptly actual relax level or actual damage degree, current actual virtual lesion vector
d i In numerical value be not that the corresponding supporting rope of 0 element is exactly problematic supporting rope, problematic supporting Suo Keneng is slack line, also possibly is damaged cable, its numerical response the degree of lax or damage;
M. from the problematic supporting rope that l identified the step, identify damaged cable, remaining is exactly slack line.
N. be utilized in the current actual virtual lesion vector that the l step obtains
d i Obtain the current actual virtual lesion degree of slack line, be utilized in the current cable force vector that the i step obtains
F i , be utilized in the volume coordinates that i goes on foot two supporting end points of all supporting ropes that obtain, be utilized in the initial drift vector of e step acquisition
l o Be utilized in elastic modulus, density, the initial cross sectional area data of all ropes of e step acquisition; Through with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual virtual lesion degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa.Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of supporting rope drift, has just confirmed the long adjustment amount of rope of the supporting rope that those need adjust Suo Li.The lax identification and the damage identification of support cable have so just been realized.Institute's demand power is by the current cable force vector during calculating
F i Corresponding element provides.
O. try to achieve current nominal virtual lesion vector
d i c After, set up mark vector according to formula 3
B i , formula 4 has provided mark vector
B i jThe definition of individual element;
Element in formula 3, the formula 4
B i j It is mark vector
B i jIndividual element,
D i Uj It is nominal virtual unit damage vector
D i u jIndividual element,
d i Cj It is current nominal virtual lesion vector
d i c jIndividual element, they all represent
jThe relevant information of root rope.In the formula 4
j=1,2,3 ..., N.
If mark vector p.
B i Element be 0 entirely, then get back to g step and continue this circulation; If mark vector
B i Element be not 0 entirely, then get into next step, i.e. q step.
Q. according to formula 5 calculate next time, promptly
iThe initial virtual damage vector that+1 circulation is required
d i+
1 o Each element
d i+
1 Oj
In the formula 5
D i Uj Be
iThe nominal virtual unit damage vector of inferior circulation
D i u jIndividual element,
d i Cj Be
iThe current nominal virtual lesion vector of inferior circulation
d i c jIndividual element,
B i j Be
iInferior circulation mark vector
B i jIndividual element.In the formula 5
j=1,2,3 ..., N.
R. at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on, make the health status of rope do
d i+
1 o The back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time
I+1
S. pass through Mechanics Calculation benchmark model A
I+1Calculate corresponding to model A
I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, the initial value of the required monitored amount that promptly circulates for the i+1 time is vectorial
C I+1 o
T. set up the required current Mechanics Calculation benchmark model of the Cable Structure A that next time, promptly circulates for the i+1 time
Ti+1 o, promptly get A
Ti+1 oEqual A
I+1
The current cable structural bearings generalized coordinate vector of u. set up next time, i.e. the i+1 time circulation is required
U Ti+
1 o, promptly get
U Ti+
1 oEqual
U Ti o
V. get back to step g, beginning is circulation next time.
In step g, according to current cable structure actual measurement bearing generalized coordinate vector
U Ti , upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
Ti oWith current cable structural bearings generalized coordinate vector
U Ti oConcrete grammar be:
G1. actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector
U Ti After, relatively
U Ti With
U Ti oIf,
U Ti Equal
U Ti o, then need be to A
Ti oUpgrade;
G2. actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector
U Ti After, relatively
U Ti With
U Ti oIf,
U Ti Be not equal to
U Ti o, then need be to A
Ti oUpgrade, update method is: calculate earlier
U Ti With
U oPoor,
U Ti With
U oDifference be exactly that the current cable structural bearings is about setting up A
oThe time the current bearing generalized displacement of Cable Structure bearing, with current bearing generalized displacement vector
VThe generalized displacement of expression bearing, current bearing generalized displacement vector
VIn element and bearing generalized displacement component between be one-to-one relationship, current bearing generalized displacement vector
VIn the numerical value of an element corresponding to the rotation of an assigned direction of an appointment bearing; Upgrade A
Ti oMethod be: at A
oThe basis on make that the health status of rope is cable system initial damage vector
d i o , more further to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint, the numerical value of current bearing generalized displacement constraint is just taken from current bearing generalized displacement vector
VThe numerical value of middle corresponding element is to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
Ti o, upgrade A
Ti oThe time,
U Ti oAll elements numerical value is also used
U Ti All elements numerical value replaces, and has promptly upgraded
U Ti o, so just obtained correctly corresponding to A
Ti o U Ti o
In step h, at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining the monitored quantitative change matrix of Cable Structure virtual unit damage
Δ C i With nominal virtual unit damage vector
D i u Concrete grammar be:
H1. the i time when beginning circulation, directly h2 obtains the monitored quantitative change matrix of Cable Structure virtual unit damage to the listed method of step h4 set by step
Δ C i With nominal virtual unit damage vector
D i u In the moment of non-the i time circulation beginning, in step g to A
Ti oAfter upgrading, directly h2 obtains the monitored quantitative change matrix of Cable Structure virtual unit damage to the listed method of step h4 set by step
Δ C i With nominal virtual unit damage vector
D i u In the moment of non-the i time circulation beginning, if in step g not to A
Ti oUpgrade, then directly change step I herein over to and carry out follow-up work;
H2. at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have
NThe root rope just has
NInferior calculating; Calculating each time in the hypothesis cable system has only a rope on the basis of original virtual lesion, to increase virtual unit damage again; The rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating; And supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
D i u " write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is each time formed one " the current numerical value vector of the calculating of monitored amount ".When hypothesis the
jWhen the root rope has unit damage, available
C i Tj " the current evaluation vector of monitored amount that expression is corresponding
C i Tj ".When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target;
H3. that calculates each time " the current evaluation vector of monitored amount
C i Tj " deduct " initial value of monitored amount vector
C i o " obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; Have
NThe root rope just has
NIndividual " the numerical value change vector of monitored amount ";
H4. by this
NIndividual " the numerical value change vector of monitored amount " formed successively to be had
N" the monitored numerical quantity transformation matrices of virtual unit damage of row
Δ C i "; " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " each row corresponding to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector
d i c With current actual virtual lesion vector
d i The element coding rule identical.
Beneficial effect:System and method disclosed by the invention occurs under the situation of generalized displacement at the Cable Structure bearing, having under the synchronously impaired or lax condition of more rope monitoring and evaluation very exactly go out the health status (position and relax level or the degree of injury that comprise all slack lines and damaged cable) of cable system.This is because " the current numerical value vector of monitored amount
C i " " initial value of monitored amount is vectorial together
C i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between linear relationship be similar to; be actually nonlinear relation; when particularly more the or extent of damage is big at damaged cable; the nonlinear characteristic of the relation between the above-mentioned amount is more obvious,, the invention discloses a kind of health monitor method that in the minizone, approaches this nonlinear relationship with linear relationship for overcoming this obstacle.In fact the present invention has used the method with linear relationship piecewise approximation nonlinear relationship; Big interval is divided into minizone one by one; In each minizone internal linear relation all is enough accurately; The health status of the cable system that obtains according to its judgement also is reliably, and therefore system and method disclosed by the invention is very useful to effective health monitoring of cable system.
Embodiment
When the bearing generalized displacement,, the invention discloses a kind of system and method for health status of each root rope of the cable system that can monitor Cable Structure rationally and effectively to the health monitoring of the cable system of Cable Structure.The following explanation of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
Occur at the Cable Structure bearing under the situation of generalized displacement, the present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status (relax level and the extent of damage that comprise rope) of the cable system of Cable Structure.During practical implementation, the following step is a kind of in the various steps that can take.
The first step: confirm type, position and the quantity of monitored amount, and numbering.Detailed process is:
At first confirm the coding rule of rope, all ropes are numbered by this rule.This numbering will be used to generate the vector sum matrix in subsequent step.
If total N root rope, the coding rule of at first definite rope, with rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule.
Confirm that (promptly all characterize the specified point of structure angle displacements, are provided with for the measured point of appointment
KIndividual specified point), give all specified point numberings; Confirmed that (establishing each measurement point had for the measured straight line of each measurement point
LIndividual appointment straight line), give the measured straight line numbering of all appointments; Confirm that (establishing each measured straight line has for the measured angle coordinate component of each measured straight line
HIndividual angle coordinate component), give all measured angle coordinate component numberings.Above-mentioned numbering will be used to generate the vector sum matrix equally in subsequent step." the whole monitored angle-data of structure " is by on top definite structure
KIndividual specified point, cross each specified point
LIndividual appointment straight line, each specifies straight line
HIndividual angle coordinate component is described, and the variation of structure angle is exactly the variation of the angle coordinate component of all appointments all specified points, all appointment straight lines.Each total
M (M=K * L * H)Individual angle coordinate component measurement value or calculated value characterize the angle information of structure.
KWith
MMust not be less than the quantity of supporting rope
NFor simplicity, in the present invention " the monitored angle-data of structure " is called " monitored amount ".Can only measure an angle coordinate of specifying straight line at each specified point, for example measure the angle coordinate of the body structure surface normal of specified point with respect to the acceleration of gravity direction, in fact be exactly measurement of dip angle here.
Second step: the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up initial virtual damage vector
d 1 o If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, perhaps can think when the structure original state is not damaged, no relaxed state vector
d 1 o Each element numerical value get 0.
The 3rd step: setting up initial virtual damage vector
d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms " the initial value vector of monitored amount
C 1 o "; Simultaneously, directly measure the initial Suo Li of all supporting ropes that calculate Cable Structure, form " initial rope force vector
F o "; Simultaneously, obtain the initial drift of all ropes, form " the initial drift vector of supporting rope according to structural design data, completion data
l o "; Simultaneously, actual measurement or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information.
The 4th step: setting up initial virtual damage vector
d 1 o The time, can adopt ripe measuring method to carry out cable force measurement, strain measurement, measurement of angle and volume coordinate and measure.Calculate Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) after directly measuring or measuring; The original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is to confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
Set up the initial Mechanics Calculation benchmark model A of Cable Structure
o, set up initial Cable Structure bearing generalized coordinate vector
U o, the Mechanics Calculation benchmark model A of the Cable Structure that needs when setting up circulation beginning for the first time
1The measured data of the Cable Structure in Cable Structure completion; This measured data comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data; And the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope; According to design drawing and as-constructed drawing, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure
oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, this Cable Structure is surveyed; Obtain the measured data of Cable Structure equally; According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure equally
oNo matter which kind of method to obtain A with
o, based on A
oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween generally must not be greater than 5%; Corresponding to A
oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector
U oA
oWith
U oBe constant, only when circulation beginning for the first time, set up; The Mechanics Calculation benchmark model of the Cable Structure of setting up during the i time circulation beginning is designated as A
i, wherein i representes cycle index; Alphabetical i is except the place of representing number of steps significantly in the application form of the present invention, and alphabetical i only representes cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A
1, A among the present invention
1Just equal A
oFor narrating conveniently name " the current Mechanics Calculation benchmark model of Cable Structure A
Ti o", A in circulation each time
Ti oCan bring in constant renewal in as required, during circulation beginning each time, A
Ti oEqual A
iIt is convenient to be similarly narration, name " Cable Structure actual measurement bearing generalized coordinate vector
U Ti ", in circulation each time, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector
U Ti , vector
U Ti Element with the vector
U oThe generalized coordinate of the equidirectional of the element representation same abutment of same position; For the purpose of narrating conveniently,, the last time is upgraded A for the i time circulation
Ti oThe time Cable Structure bearing generalized coordinate current data be designated as current cable structural bearings generalized coordinate vector
U Ti oDuring circulation beginning for the first time, A
T1 oEqual A
1,
U T1 oEqual
U oA
1The health status of corresponding rope by
d 1 o Describe; Mechanics Calculation benchmark model A
iThe health status of corresponding rope by
d i o Describe.
The 5th step: the hardware components of pass line structural healthy monitoring system.Hardware components comprises at least: monitored amount monitoring system (for example containing measurement of angle sensor, signal conditioner etc.), Cable Structure bearing generalized coordinate monitoring system (for example containing total powerstation, angular transducer, signal conditioner etc.), cable force monitoring system (for example containing acceleration transducer, signal conditioner etc.), respectively support horizontal range monitoring system, signal (data) collector, the computing machine and the panalarm of communicating by letter of rope two supporting end points.The horizontal range of the Suo Li of each monitored amount, each supporting rope and each root supporting rope two supporting end points all must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector with the signal that monitors; Signal is delivered to computing machine through signal picker; Computing machine then is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; Have lax or during damage when monitoring rope, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 6th step: establishment and the cable system health monitoring systems software of pass line structure on supervisory control comuter.All move this software at circulation time each time, this software is all the time in operation in other words.This software will be accomplished functions such as monitoring that the present invention's progressive method of the identification slack line of angle monitor " during bearing generalized displacement based on " required by task wants, record, control, storage, calculating, notice, warning (all work that can accomplish with computing machine in this practical implementation method); And can regularly or by the personnel operation health monitoring systems generate cable system health condition form; Can also be according to the condition of setting (for example damage reach a certain value), notice or prompting monitor staff notify specific technician to accomplish necessary evaluation work automatically.
The 7th step: the step begins circulation running thus, for narration conveniently is designated as the i time circulation, and i=1 wherein, 2,3,4,5 ...Actual measurement (for example measuring the angle coordinate of the outer normal of seating plane with angular transducer) obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector
U Ti , according to current cable structure actual measurement bearing generalized coordinate vector
U Ti , upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
Ti oWith current cable structural bearings generalized coordinate vector
U Ti oConcrete grammar is:
During circulation beginning each time, make A
Ti oEqual A
i
Actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector
U Ti After, relatively
U Ti With
U Ti oIf,
U Ti Equal
U Ti o, then need be to A
Ti oUpgrade;
Actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector
U Ti After, relatively
U Ti With
U Ti oIf,
U Ti Be not equal to
U Ti o, then need be to A
Ti oUpgrade, update method is: calculate earlier
U Ti With
U oPoor,
U Ti With
U oDifference be exactly that the current cable structural bearings is about setting up A
oThe time the current bearing generalized displacement of Cable Structure bearing, with current bearing generalized displacement vector
VThe generalized displacement of expression bearing, current bearing generalized displacement vector
VIn element and bearing generalized displacement component between be one-to-one relationship, current bearing generalized displacement vector
VIn the numerical value of an element corresponding to the rotation of an assigned direction of an appointment bearing; Upgrade A
Ti oMethod be: at A
oThe basis on make that the health status of rope is cable system initial damage vector
d i o , more further to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint, the numerical value of current bearing generalized displacement constraint is just taken from current bearing generalized displacement vector
VThe numerical value of middle corresponding element is to A
oIn the Cable Structure bearing apply current bearing generalized displacement constraint after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
Ti o, upgrade A
Ti oThe time,
U Ti oAll elements numerical value is also used
U Ti All elements numerical value replaces, and has promptly upgraded
U Ti o, so just obtained correctly corresponding to A
Ti o U Ti o
The 8th step: at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining the monitored numerical quantity transformation matrices of Cable Structure virtual unit damage
Δ C i With nominal virtual unit damage vector
D i u Concrete grammar is:
A. the i time when beginning circulation, directly b obtains the monitored quantitative change matrix of Cable Structure virtual unit damage to the listed method of steps d set by step
Δ C i With nominal virtual unit damage vector
D i u In the moment of non-the i time circulation beginning, in the 7th step to A
Ti oAfter upgrading, directly b obtains the monitored quantitative change matrix of Cable Structure virtual unit damage to the listed method of steps d set by step
Δ C i With nominal virtual unit damage vector
D i u In the moment of non-the i time circulation beginning, if in the 7th step not to A
Ti oUpgrade, then directly changing for the 9th step herein over to carries out follow-up work;
B. at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have
NThe root rope just has
NInferior calculating; Calculating each time in the hypothesis cable system has only a rope on the basis of original virtual lesion, to increase virtual unit damage again; The rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating; And supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
D i u " write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is each time formed one " the current numerical value vector of the calculating of monitored amount ".When hypothesis the
jWhen the root rope has unit damage, available
C i Tj " the current evaluation vector of monitored amount that expression is corresponding
C i Tj ".When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target;
C. that calculates each time " the current evaluation vector of monitored amount
C i Tj " deduct " initial value of monitored amount vector
C i o " obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; Have
NThe root rope just has
NIndividual " the numerical value change vector of monitored amount ";
D. by this
NIndividual " the numerical value change vector of monitored amount " formed successively to be had
N" the monitored numerical quantity transformation matrices of virtual unit damage of row
Δ C i "; " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " each row corresponding to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector
d i c With current actual virtual lesion vector
d i The element coding rule identical.
In this step, reach when giving each vectorial element numbering thereafter; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.
The 9th step: set up the linear relationship error vector
e i And vector
g i Utilize data (" the initial value vector of monitored amount of front
C i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i "), when the 8th step calculated each time, promptly in calculating each time, have only in the hypothesis cable system increase virtual unit damage again on the basis of rope at original virtual lesion in, calculate each time and form " virtual lesion a vector
d i t ", the virtual lesion vector
d i t Element number equal the quantity of rope, the virtual lesion vector
d i t All elements in have only the numerical value of an element to get to calculate each time in hypothesis increase the virtual unit damage value of the rope of virtual unit damage,
d i t The numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the rope of virtual unit damage, be identical with the element of the same numberings of other vectors with the corresponding relation of this rope; Will
C i Tj , C i o ,
Δ C i ,
d i t Bringing formula (13) into (notes, in the formula (13)
C i Use
C i Tj Bring into,
d i c Use
d i t Bring into), obtain a linear relationship error vector
e i , calculate a linear relationship error vector each time
e i Have
NThe root rope just has
NInferior calculating just has
NIndividual linear relationship error vector
e i , with this
NIndividual linear relationship error vector
e i Obtain a vector after the addition, with this vector each element divided by
NAfter the new vector that obtains be exactly final linear relationship error vector
e i Vector
g i Equal final error vector
e i With vector
g i Be kept on the hard disc of computer of operation health monitoring systems software, supply health monitoring systems software to use.
Will " initial rope force vector
F o ", " initial value of monitored amount vector
C i o ", " nominal virtual unit damage vector
D i u ", " initial drift vector
l o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and the parameters such as unit weight of the elastic modulus of all ropes, initial cross sectional area, rope be kept on the hard disc of computer of operation health monitoring systems software with the mode of data file.
The tenth step: actual measurement obtains the current cable power of all supporting ropes of Cable Structure, forms the current cable force vector
F i Simultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount
C i ".Actual measurement calculates the space coordinates of two supporting end points of all support cables, and the space coordinates of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.
The 11 step: according to " current (calculating or actual measurement) numerical value vector of monitored amount
C i " " initial value of monitored amount is vectorial together
C i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between the linear approximate relationship (seeing formula (9)) that exists, calculate the current nominal virtual lesion vector of cable system according to multi-objective optimization algorithm
d i c Noninferior solution.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted Sum Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution current nominal virtual lesion vector with the goal programming method only
d i c Process, the concrete implementation procedure of other algorithm can realize according to the requirement of its specific algorithm in a similar fashion.
According to the goal programming method, formula (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (29) and the formula (30), in the formula (29)
γ i Be a real number,
RBe real number field, area of space Ω has limited vector
d i c Span (the present embodiment requirements vector of each element
d i c Each element be not less than 0, be not more than 1).The meaning of formula (29) is to seek the minimum real number of an absolute value
γ i , make formula (30) be met.In the formula (30)
G (d i c )By formula (31) definition, weighing vector in the formula (30)
W i With
γ i Product representation formula (30) in
G (d i c )With vector
g i Between the deviation that allows,
g i Definition referring to formula (15), its value will the 8th the step calculate.Vector during actual computation
W i Can with vector
g i Identical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of current name damage vector
d i c
Try to achieve current nominal virtual lesion vector
d i c After
,Can be vectorial according to the current actual virtual lesion that formula (17) obtain
d i Each element, current actual virtual lesion vector
d i Have reasonable error exactly but can be more exactly from all ropes, confirm the position of problematic rope (be virtual damaged cable, possibly be impaired also possibly be lax) and separating of virtual lesion degree.If the current actual virtual lesion vector that solves
d i The numerical value of a certain element be 0, represent that the pairing rope of this element is intact, not damage or lax; If its numerical value is 100%, represent that then the pairing rope of this element has completely lost load-bearing capacity; If its numerical value between 0 and 100%, is then represented this rope and has been lost the load-bearing capacity of corresponding proportion.
The 12 step: because current actual virtual lesion vector
d i Element numerical value represent the virtual lesion degree of corresponding rope; So it is impaired or relaxed and possible degree of injury or relax level just to define which Suo Keneng according to current actual virtual lesion vector; But damage has taken place actually or has taken place to relax in these ropes, need differentiate.The method of differentiating is varied; Can be through removing the protective seam of supporting rope; To the visual discriminating of supporting Suo Jinhang; Perhaps carry out visual discriminating by optical imaging apparatus, also can be through lossless detection method to supporting rope impaired discriminating the whether, UT (Ultrasonic Testing) is exactly a kind of present widely used lossless detection method.Those do not find to damage and the virtual lesion degree is not that 0 supporting rope is exactly that lax rope has taken place to differentiate the back, need adjust the rope of Suo Li exactly, can be accomplished the finding the solution of relax level (being the long adjustment amount of rope) of slack line by software.
The 13 step: in this circulation, promptly the
iTry to achieve current nominal virtual lesion vector in the inferior circulation
d i c After, at first, set up mark vector according to formula (26), formula (27)
B i
The 14 step: if mark vector
B i Element be 0 entirely, then got back to for the 7th step and continue this circulation; If mark vector
B i Element be not 0 entirely, then get into next step, i.e. the 15 step.
The 15 step: according to formula (28) calculate next time, promptly the
iThe initial virtual damage vector that+1 circulation is required
d i+
1 o Each element
d i+
1 Oj
The 16 step: at the current Mechanics Calculation benchmark model of Cable Structure A
Ti oThe basis on, the health status that makes rope is the last vector that calculates of step
d i+
1 o After, obtain new Mechanics Calculation benchmark model, next time promptly the required Mechanics Calculation benchmark model A of (the i+1 time) circulation
I+1
The 17 step: through to Mechanics Calculation benchmark model A
I+1Calculate corresponding to model A
I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, required vector promptly circulates for the i+1 time
C I+1 o , promptly the initial value of monitored amount is vectorial
The 18 step: set up next time, the required current Mechanics Calculation benchmark model of Cable Structure A promptly circulates for the i+1 time
Ti+1 o, promptly get A
Ti+1 oEqual A
I+1
The 19 step: the current cable structural bearings generalized coordinate vector of set up next time, i.e. the i+1 time circulation is required
U Ti+
1 o, promptly get
U Ti+
1 oEqual
U Ti o
The 20 step: got back to for the 7th step, beginning is circulation next time.
Claims (3)
- During a bearing generalized displacement based on the progressive method of the identification slack line of angle monitor, it is characterized in that this method comprises:A. establish total N root rope, at first confirm the coding rule of rope, with rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule;B. confirm the measured point of appointment, give all specified point numberings; Confirm the measured straight line of each measurement point, gave the measured straight line numbering of all appointments; Confirm the measured angle coordinate component of each measured straight line, give all measured angle coordinate component numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored angle-data of structure " is made up of above-mentioned all measured angle coordinate components; For simplicity, " the monitored angle-data of structure " abbreviated as " monitored amount "; The quantity of measurement point must not be less than the quantity of rope; The quantity sum of all measured angle coordinate components must not be less than the quantity of rope;C. the data that Non-Destructive Testing data of utilizing rope etc. can be expressed the health status of rope are set up initial virtual and are damaged vectorial d i o, wherein i representes cycle index, back i and subscript i represent cycle index, and i=1,2,3 ...; Circulation time d for the first time i oBe designated as d 1 oIf when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think perhaps that the structure original state is lax for not having, during the not damaged state, vectorial d 1 oEach element numerical value get 0;D. damage vectorial d setting up initial virtual 1 oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount 1 oE. damage vectorial d setting up initial virtual 1 oInitial value vector C with monitored amount 1 oThe time, directly measure the initial Suo Li that calculates all supporting ropes, form initial rope force vector F oSimultaneously, obtain the initial drift that all support ropes, form initial drift vector l according to structural design data, completion data oSimultaneously, obtain the initial geometric data of Cable Structure according to structural design data, completion data or actual measurement; Simultaneously, actual measurement or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information;F. set up the initial Mechanics Calculation benchmark model A of Cable Structure o, set up initial Cable Structure bearing generalized coordinate vector U o, the Mechanics Calculation benchmark model A of the Cable Structure that needs when setting up circulation beginning for the first time 1The measured data of the Cable Structure in Cable Structure completion; This measured data comprises measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data, all ropes; And the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope; According to design drawing and as-constructed drawing, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, this Cable Structure is surveyed; Obtain the measured data of Cable Structure equally; According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial Mechanics Calculation benchmark model A of Cable Structure equally oNo matter which kind of method to obtain A with o, based on A oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U oA oAnd U oBe constant, only when circulation beginning for the first time, set up; The Mechanics Calculation benchmark model of the Cable Structure of setting up during the i time circulation beginning is designated as A i, wherein i representes cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A 1, A among the present invention 1Just equal A oFor narrating conveniently name " the current Mechanics Calculation benchmark model of Cable Structure A Ti o", A in circulation each time Ti oCan bring in constant renewal in as required, during circulation beginning each time, A Ti oEqual A iIt is convenient to be similarly narration, name " Cable Structure actual measurement bearing generalized coordinate vector U Ti", in circulation each time, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector U Ti, vectorial U TiElement and vectorial U oThe generalized coordinate of the equidirectional of the element representation same abutment of same position; For the purpose of narrating conveniently,, the last time is upgraded A for the i time circulation Ti oThe time Cable Structure bearing generalized coordinate current data be designated as current cable structural bearings generalized coordinate vector U Ti oDuring circulation beginning for the first time, A T1 oEqual A 1, U T1 oEqual U oA 1The health status of corresponding rope is by d 1 oDescribe; Mechanics Calculation benchmark model A iThe health status of corresponding rope is by d i oDescribe; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;When g. circulation begins each time, make A Ti oEqual A iActual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas are formed current cable structure actual measurement bearing generalized coordinate vector U Ti, according to current cable structure actual measurement bearing generalized coordinate vector U Ti, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary Ti oWith current cable structural bearings generalized coordinate vector U Ti oH. at the current Mechanics Calculation benchmark model of Cable Structure A Ti oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining the monitored numerical quantity transformation matrices of Cable Structure virtual unit damage Δ C iWith nominal virtual unit damage vector D i uI. actual measurement obtains the current cable power of all supporting ropes of Cable Structure, forms current cable force vector F iSimultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount i".Actual measurement calculates the space coordinates of two supporting end points of all support cables, and the space coordinates of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.When numbering to the element of the institute's directed quantity that occurred before this step and this step; Should use same coding rule; Can guarantee like this before this step and this step with each vector that occurs afterwards, number identical element, represent same monitored amount, corresponding to vectorial defined relevant information under this element;J. define current nominal virtual lesion vector d to be asked i cWith current actual virtual lesion vector d iDamage vectorial d i o, d i cAnd d iElement number equal the quantity of rope, be one-to-one relationship between the element of damage vector and the rope, the element numerical value of damage vector is represented the virtual lesion degree or the health status of corresponding rope;K. according to " the current numerical value vector C of monitored amount i" " the initial value of monitored amount vector C together i o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C i" and " current nominal virtual lesion vector d i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 i cOther outer amount is known, finds the solution formula 1 and just can calculate current nominal virtual lesion vector d i cL. the current actual virtual lesion vector d that utilizes formula 2 to express iElement d i jDamage vectorial d with initial virtual i oElement d i OjWith current nominal virtual lesion vector d i cElement d i CjBetween relation, calculate current actual virtual lesion vector d iAll elements,J=1 in the formula 2,2,3 ..., N,Because current actual virtual lesion vector d iElement numerical value represent the current actual virtual lesion degree of corresponding rope, promptly actual relax level or actual damage degree, current actual virtual lesion vector d iIn numerical value be not that the corresponding supporting rope of 0 element is exactly problematic supporting rope, problematic supporting Suo Keneng is slack line, also possibly is damaged cable, its numerical response the degree of lax or damage;M. from the problematic supporting rope that l identified the step, identify damaged cable, remaining is exactly slack line;N. be utilized in the current actual virtual lesion vector d that the l step obtains iObtain the current actual virtual lesion degree of slack line, be utilized in the current cable force vector F that the i step obtains i, all that are utilized in i step acquisition support two volume coordinates that support end points of ropes, are utilized in the vectorial l of initial drift that the e step obtains oBe utilized in elastic modulus, density, the initial cross sectional area data of all ropes of e step acquisition; Through with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual virtual lesion degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two supporting ropes in structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition; Relax level is exactly the change amount of supporting rope drift; Just confirmed the long adjustment amount of rope of the supporting rope that those need adjust Suo Li; So just realized the lax identification and the damage identification of supporting rope, institute's demand power is by current cable force vector F during calculating iCorresponding element provides;O. try to achieve current nominal virtual lesion vector d i cAfter, set up mark vector B according to formula 3 i, formula 4 has provided mark vector B iThe definition of j element;Element B in formula 3, the formula 4 i jBe mark vector B iJ element, D i UjBe nominal virtual unit damage vector D i uJ element, d i CjBe current nominal virtual lesion vector d i cJ element, they all represent the relevant information of j root rope.J=1 in the formula 4,2,3 ..., N;If mark vector B p. iElement be 0 entirely, then get back to g step and continue this circulation; If mark vector B iElement be not 0 entirely, then get into next step, i.e. q step;Q. according to formula 5 calculate next time, promptly the i+1 time the circulation required initial virtual damage vectorial d I+1 oEach element d I+1 OjD in the formula 5 i UjBe the nominal virtual unit damage vector of the i time circulation D i uJ element, d i CjBe the current nominal virtual lesion vector d of the i time circulation i cJ element, B i jBe the i time circulation mark vector B iJ element.J=1 in the formula 5,2,3 ..., N;R. at the current Mechanics Calculation benchmark model of Cable Structure A Ti oThe basis on, the health status that makes rope is d I+1 oThe back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time I+1S. pass through Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, the vectorial C of initial value of the required monitored amount that promptly circulates for the i+1 time I+1 oT. set up the required current Mechanics Calculation benchmark model of the Cable Structure A that next time, promptly circulates for the i+1 time Ti+1 o, promptly get A Ti+1 oEqual A I+1The current cable structural bearings generalized coordinate vector U of u. set up next time, i.e. the i+1 time circulation is required Ti+1 o, promptly get U Ti+1 oEqual U Ti oV. get back to step g, beginning is circulation next time.
- 2. based on the progressive method of the identification slack line of angle monitor, it is characterized in that in step g, during bearing generalized displacement according to claim 1 according to current cable structure actual measurement bearing generalized coordinate vector U Ti, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary Ti oWith current cable structural bearings generalized coordinate vector U Ti oConcrete grammar be:G1. actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector U TiAfter, compare U TiAnd U Ti oIf, U TiEqual U Ti o, then need be to A Ti oUpgrade;G2. actual measurement obtains current cable structure actual measurement bearing generalized coordinate vector U TiAfter, compare U TiAnd U Ti oIf, U TiBe not equal to U Ti o, then need be to A Ti oUpgrade, update method is: calculate U earlier TiWith U oPoor, U TiWith U oDifference be exactly that the current cable structural bearings is about setting up A oThe time the current bearing generalized displacement of Cable Structure bearing; Represent the bearing generalized displacement with current bearing generalized displacement vector V; Between element among the current bearing generalized displacement vector V and the bearing generalized displacement component is one-to-one relationship, and the numerical value of an element is corresponding to the rotation of an assigned direction of an appointment bearing among the current bearing generalized displacement vector V; Upgrade A Ti oMethod be: at A oThe basis on make that the health status of rope is cable system initial damage vector d i o, more further to A oIn the Cable Structure bearing apply current bearing generalized displacement constraint, the numerical value of current bearing generalized displacement constraint is just taken from the numerical value of corresponding element among the current bearing generalized displacement vector V, to A oIn the Cable Structure bearing apply current bearing generalized displacement constraint after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal Ti o, upgrade A Ti oThe time, U Ti oAll elements numerical value is also used U TiAll elements numerical value replaces, and has promptly upgraded U Ti o, so just obtained correctly corresponding to A Ti oU Ti o
- 3. based on the progressive method of the identification slack line of angle monitor, it is characterized in that in step h, during bearing generalized displacement according to claim 1 at the current Mechanics Calculation benchmark model of Cable Structure A Ti oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining the monitored quantitative change matrix of Cable Structure virtual unit damage Δ C iWith nominal virtual unit damage vector D i uConcrete grammar be:H1. the i time when beginning circulation, directly h2 obtains the monitored quantitative change matrix of Cable Structure virtual unit damage Δ C to the listed method of step h4 set by step iWith nominal virtual unit damage vector D i uIn the moment of non-the i time circulation beginning, in step g to A Ti oAfter upgrading, directly h2 obtains the monitored quantitative change matrix of Cable Structure virtual unit damage Δ C to the listed method of step h4 set by step iWith nominal virtual unit damage vector D i uIn the moment of non-the i time circulation beginning, if in step g not to A Ti oUpgrade, then directly change step I herein over to and carry out follow-up work;H2. at the current Mechanics Calculation benchmark model of Cable Structure A Ti oThe basis on carry out the several times Mechanics Calculation; Equal the quantity of all ropes on the calculation times numerical value; There is N root rope that N calculating is just arranged; Calculating each time in the hypothesis cable system has only a rope on the basis of original virtual lesion, to increase virtual unit damage again; The rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating, and supposes that each time the virtual unit damage value of the rope that virtual unit damage is arranged can be different from the virtual unit damage value of other ropes, with " nominal virtual unit damage vector D i u" write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is each time formed one " the current numerical value vector of the calculating of monitored amount "; When hypothesis j root rope has unit damage, available C i Tj" the current evaluation vector C of monitored amount that expression is corresponding i Tj"; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target;H3. that calculates each time " the current evaluation vector C of monitored amount i Tj" deduct " initial value of monitored amount vector C i o" obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; There is N root rope that N " the numerical value change vector of monitored amount " just arranged;H4. form " the monitored numerical quantity transformation matrices of the virtual unit damage Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount " i"; " the monitored numerical quantity transformation matrices of virtual unit damage Δ C i" each row corresponding to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector d i cWith current actual virtual lesion vector d iThe element coding rule identical.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110143101A CN102323089A (en) | 2011-05-31 | 2011-05-31 | Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110143101A CN102323089A (en) | 2011-05-31 | 2011-05-31 | Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102323089A true CN102323089A (en) | 2012-01-18 |
Family
ID=45450878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110143101A Pending CN102323089A (en) | 2011-05-31 | 2011-05-31 | Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102323089A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706642A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Slack cable identification method on basis of angle monitoring during support settlement and temperature variation |
CN102706622A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Loosened cable progressive-type identification method based on angle monitoring at moment of generalized displacement of supporting seat and temperature variation |
CN102706607A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Damaged cable approximant identification method on the basis of angle monitoring during temperature variation |
CN102706586A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Slack cable progressive identification method on basis of angle monitoring during support settlement and temperature variation |
CN104990748A (en) * | 2015-07-23 | 2015-10-21 | 东南大学 | Simplified method of progressively recognizing damage cable load generalized displacement through angle monitoring |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000213929A (en) * | 1999-01-22 | 2000-08-04 | Ishikawajima Harima Heavy Ind Co Ltd | Remote monitoring system for structure |
WO2008123510A1 (en) * | 2007-03-30 | 2008-10-16 | Kyoto University | Displacement measuring method, displacement measuring apparatus and target for displacement measurement |
CN101782946A (en) * | 2010-03-17 | 2010-07-21 | 东南大学 | Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement |
CN101788404A (en) * | 2010-03-17 | 2010-07-28 | 东南大学 | Progressive method for identifying loose support cable based on angle monitoring during support settlement |
CN101793624A (en) * | 2010-03-17 | 2010-08-04 | 东南大学 | Progressive method for distinguishing slack supporting cable based on mixed monitoring during support settlement |
-
2011
- 2011-05-31 CN CN201110143101A patent/CN102323089A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000213929A (en) * | 1999-01-22 | 2000-08-04 | Ishikawajima Harima Heavy Ind Co Ltd | Remote monitoring system for structure |
WO2008123510A1 (en) * | 2007-03-30 | 2008-10-16 | Kyoto University | Displacement measuring method, displacement measuring apparatus and target for displacement measurement |
CN101782946A (en) * | 2010-03-17 | 2010-07-21 | 东南大学 | Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement |
CN101788404A (en) * | 2010-03-17 | 2010-07-28 | 东南大学 | Progressive method for identifying loose support cable based on angle monitoring during support settlement |
CN101793624A (en) * | 2010-03-17 | 2010-08-04 | 东南大学 | Progressive method for distinguishing slack supporting cable based on mixed monitoring during support settlement |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102706642A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Slack cable identification method on basis of angle monitoring during support settlement and temperature variation |
CN102706607A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Damaged cable approximant identification method on the basis of angle monitoring during temperature variation |
CN102706586A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Slack cable progressive identification method on basis of angle monitoring during support settlement and temperature variation |
CN102706607B (en) * | 2012-05-29 | 2015-07-29 | 东南大学 | Based on the approximant recognition methods of the damaged cable of angle monitor during temperature variation |
CN102706586B (en) * | 2012-05-29 | 2015-07-29 | 东南大学 | The slack line progressive-type recognition method of angle monitor when support settlement and temperature variation |
CN102706642B (en) * | 2012-05-29 | 2015-10-07 | 东南大学 | Based on the slack line recognition methods of angle monitor when having support settlement and a temperature variation |
CN102706622A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Loosened cable progressive-type identification method based on angle monitoring at moment of generalized displacement of supporting seat and temperature variation |
CN104990748A (en) * | 2015-07-23 | 2015-10-21 | 东南大学 | Simplified method of progressively recognizing damage cable load generalized displacement through angle monitoring |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101788403B (en) | Progressive method for identifying loose support cable based on strain monitoring during support settlement | |
CN101782946B (en) | Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement | |
CN101782945B (en) | Method for identifying loose supporting ropes based on space coordinate monitoring during support settlement | |
CN101793623B (en) | Progressive method for distinguishing slack supporting cable based on cable force monitoring during support settlement | |
CN102323090A (en) | Progressive method for identifying slack cable on the basis of cable force monitoring during support generalized displacement | |
CN101788402B (en) | Method for indentifying loose supporting cable based on angle monitoring during support settlement | |
CN101782473B (en) | Progressive type healthy monitoring method of rope system based on angle monitoring during support settlement | |
CN101788401B (en) | Method for indentifying loose supporting cable based on strain monitoring during support settlement | |
CN101699248B (en) | Progressive health monitoring method for identifying loose carrying cable based on angle monitoring | |
CN101788404B (en) | Progressive method for identifying loose support cable based on angle monitoring during support settlement | |
CN102323089A (en) | Progressive method for identifying slack cable on the basis of angle monitoring during support generalized displacement | |
CN102323088A (en) | During the bearing generalized displacement based on the progressive method of the identification slack line of hybrid monitoring | |
CN101793624B (en) | Progressive method for distinguishing slack supporting cable based on mixed monitoring during support settlement | |
CN102323075A (en) | Progressive method for identifying damaged cable, slack cable and support generalized displacement on the basis of angle monitoring | |
CN101793622B (en) | Method for distinguishing slack supporting cable based on cable force monitoring during support settlement | |
CN102288433A (en) | Progressive method for distinguishing slack cable based on angle monitoring in angular displacement of supporting seat | |
CN102323092A (en) | Cable monitoring based loose cable monitoring method applied in supporting seat generalized displacement | |
CN102288428A (en) | Method for distinguishing slack cable based on hybrid monitoring in angular displacement of supporting seat | |
CN102252863B (en) | Strain monitoring-based progressive method for identifying slack cables during generalized displacement of support saddle | |
CN102288430B (en) | Progressive method for recognizing slack cable based on hybrid monitoring during angular displacement of support | |
CN102288436B (en) | Method for identifying slack supporting cables based on space coordinate monitoring during angular displacement of support | |
CN102323099A (en) | Space coordinate monitoring based loose cable progressive recognizing method applied in supporting seat generalized displacement | |
CN102297775A (en) | Progressive method for identifying loose cables based on cable stress monitoring during supporting-base angular displacement | |
CN102323093A (en) | Angle monitoring based cable system progressive health monitoring method applied in supporting seat generalized displacement | |
CN102288434B (en) | Method for identifying slack cables based on angle monitoring during angular displacement of support |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20120118 |