Summary of the invention
Technical matters: the objective of the invention is when the Cable Structure bearing has sedimentation,, disclose a kind of health monitor method when support settlement is arranged based on the cable system of hybrid monitoring at the health monitoring problem of cable system in the Cable Structure.
Technical scheme: the present invention is made up of three parts.Be respectively the method for setting up required knowledge base of cable system health monitoring systems and parameter, based on knowledge base (containing parameter) with based on to the hybrid monitoring of multiclass parameter and the cable system health status appraisal procedure of actual measurement Cable Structure support coordinate, the software and hardware part of health monitoring systems.
First of the present invention: foundation is used for the knowledge base of cable system health monitoring and the method for parameter.Specific as follows:
1. set up the initial mechanical calculating benchmark model A of Cable Structure
o(for example finite element benchmark model) and current Mechanics Calculation benchmark model A
t oThe method of (for example finite element benchmark model).A in the present invention
oBe constant.A
t oBring in constant renewal in.Set up A
o, set up and upgrade A
t oMethod as follows.
Set up A
oThe time, measured data according to the Cable Structure in the Cable Structure completion (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support 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) and design drawing, as-constructed drawing, 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 the Cable Structure shape data, the rope force data, the draw-bar pull data, Cable Structure support coordinate data, measured datas such as Cable Structure modal data, to cable-stayed bridge, suspension bridge and the bridge type data of Yan Shiqiao, the 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), design drawing according to these data and Cable Structure, as-constructed drawing utilizes 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 support coordinate data form initial Cable Structure support coordinate vector U
o
In Cable Structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data (all data composition current cable structure actual measurement support coordinate vector U
t, vectorial U
tDefinition mode identical with vectorial U).Cable Structure support coordinate current data when for simplicity, being upgraded current Mechanics Calculation benchmark model the last time is designated as current cable structural bearings coordinate vector U
t oSet up and renewal A
t oMethod be: at initial time, A
t oJust equal A
o, U
t oJust equal U
oIn Cable Structure military service process, constantly actual measurement obtains Cable Structure support coordinate data and obtains current cable structure actual measurement support coordinate vector U
tIf, U
tEqual U
t o, then do not need A
t oUpgrade; If U
tBe not equal to U
t o, then need A
t oUpgrade, at this moment U
tWith U
oDifference be exactly the Cable Structure bearing about initial position (corresponding to A
o) support displacement (representing support displacement with support displacement vector V, is exactly support settlement in the displacement of gravity direction), upgrade A
t oMethod be to A
oIn the Cable Structure bearing apply the current Mechanics Calculation benchmark model A that obtains upgrading after the displacement constraint (its numerical value take from support displacement vector V)
t o, upgrade A
t oThe time, U
t oAll elements numerical value is also used U
tAll elements numerical value replaces, and has promptly upgraded U
t o, so just obtained correctly corresponding to A
t oU
t oMonitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate are described below respectively:
If total N root rope in the cable system, the monitored rope force data of structure is by M on the structure
1The M of individual appointment rope
1Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope.Each total M
1Individual cable force measurement value or calculated value characterize the rope force information of structure.M
1Be one and be not less than 0 integer.
The monitored strain data of structure can be by K on the structure
2L individual specified point, that reach each specified point
2The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
2The variation of the tested strain of all of individual specified point.Each total M
2(M
2=K
2* L
2) individual strain measurement value or calculated value characterize structural strain.M
2Be one and be not less than 0 integer.
The monitored angle-data of structure is by K on the structure
3L individual specified point, that cross each specified point
3H individual appointment straight line, each appointment straight line
3Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M
3(M
3=K
3* L
3* H
3) individual angle coordinate component measurement value or calculated value characterize the angle information of structure.M
3Be one and be not less than 0 integer.
The monitored shape data of structure is by K on the structure
4L individual specified point, that reach each specified point
4The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
4The variation of all coordinate components of individual specified point.Each total M
4(M
4=K
4* L
4) individual measurement of coordinates value or calculated value characterize planform.M
4Be one and be not less than 0 integer.
Comprehensive above-mentioned monitored amount, total has M (M=M
1+ M
2+ M
3+ M
4) individual monitored amount, definition parameter K (K=M
1+ K
2+ K
3+ K
4), K and M must not be less than the quantity N of rope.Because M monitored amount is dissimilar, so the present invention is called " method of discerning the support cable that needs adjustment Suo Li based on hybrid monitoring ".
For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " monitored amount ".
Among the present invention with monitored amount initial vector C
oThe vector (seeing formula (1)) that the initial value of all monitored amounts of expression Cable Structure is formed.Requirement is obtaining A
oThe time obtain C
oBecause of subject to the foregoing, the monitored amount of calculating gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial monitored amount reliably, in the narration of back, will represent this calculated value and measured value with prosign.
C
o=[C
O1C
O2C
OjC
OM]
T(1) C in the formula (1)
Oj(j=1,2,3 ...., M; M 〉=N) is the original bulk of j monitored amount in the Cable Structure, this component according to coding rule corresponding to specific j monitored amount.T represents the transposition (back together) of vector.
The vector of forming by the currency of all monitored amounts in the Cable Structure with the current numerical value vector of monitored amount C among the present invention (formula (2) is seen in definition).
C=[C
1C
2C
jC
M]
T(2) C in the formula (2)
j(j=1,2,3 ...., M; M 〉=N) is the currency of j monitored amount in the Cable Structure, this component C
jAccording to coding rule and C
OjCorresponding to same " monitored amount ".
2. set up and upgrade the method for Cable Structure unit damage monitored quantitative change matrix Δ C.
Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current Mechanics Calculation benchmark model A
t oThe time, upgrade Cable Structure unit damage monitored quantitative change matrix Δ C.Concrete grammar is as follows:
Current Mechanics Calculation benchmark model A in Cable Structure
t oThe basis on carry out several times and calculate, equal the quantity of all support cables on the calculation times numerical value.Calculating each time in the hypothesis cable system has only a support cable that unit damage D is arranged
u(for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), the rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, calculate the current calculated value that all utilizes mechanics method (for example finite element method) to calculate all monitored amounts of Cable Structure each time, (when hypothesis i root rope had unit damage, available formula (3) was represented monitored amount calculation current vector C to monitored amount calculation current vector of current calculated value composition of the monitored amount of all that calculate each time
t i); Calculate monitored amount calculation current vector each time and deduct monitored amount initial vector, the gained vector is exactly that the monitored quantitative change vector of (is mark with the position of support cable that unit damage is arranged or numbering etc.) (when i root rope has unit damage, is used δ C under this condition
iRepresent monitored quantitative change vector, formula (4) is seen in definition, formula (4) deducts formula (1) gained for formula (3)), each element representation of monitored quantitative change vector supposition owing to calculating has the change amount of the pairing monitored amount of this element that the unit damage of the Na Gensuo of unit damage causes; There is N root rope that N monitored quantitative change vector just arranged, because M monitored amount arranged, so each monitored quantitative change vector has M element, being made up of successively this N monitored quantitative change vector has the unit damage monitored of M * N element quantitative change matrix Δ C, and the definition of Δ C as the formula (5).
Elements C in the formula (3)
Tj i(i=1,2,3 ...., N; J=1,2,3 ...., M; M 〉=when N) expression has unit damage owing to i root rope, according to the current calculated amount of the pairing j of coding rule monitored amount.
Δ C in the formula (5)
J, i(i=1,2,3 ...., N; J=1,2,3 ...., M; The expression of M 〉=N) only since i root rope have that unit damage causes, according to the variation (algebraic value) of the current numerical value of calculating of the individual monitored amount of the pairing j of coding rule.Monitored quantitative changeization vector δ C
iBe actually the row among the matrix Δ C, that is to say that formula (5) also can write an accepted way of doing sth (6).
ΔC=[δC
1?δC
2…δC
i…δC
N] (6)
3. the current numerical value vector of monitored amount C (calculating or actual measurement) is with monitored amount initial vector C
o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D
uAnd the linear approximate relationship between the vectorial d of current damage, shown in (7) or formula (8).The definition of the vectorial d of the current damage of cable system is referring to formula (9).
d=[d
1?d
2…d
i…d
N]
T (9)
D in the formula (9)
i(i=1,2,3 ...., N) be the current damage of i root rope (or pull bar) in the cable system; d
iBeing to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, represents to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%.
Represented that rope thoroughly lost load-bearing capacity at 100% o'clock if establish rope damage and be, when actual damage is not too big, (for example be not more than 30% damage) so, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (7) or formula (8) is less with the error of actual conditions.Error with linear relationship shown in the linear relationship error vector e expression (7) of formula (10) definition or the formula (8).
Abs () is the function that takes absolute value in the formula (10), and each element of the vector of trying to achieve in the bracket is taken absolute value.
Second portion of the present invention: based on the cable system health status appraisal procedure of knowledge base (containing parameter) and the monitored amount of actual measurement.
Because there are certain error in formula (7) or the represented linear relationship of formula (8), therefore can not be simply directly find the solution and obtain the vectorial d of current damage according to formula (7) or formula (8) and the vectorial C of the current numerical value of the monitored amount of actual measurement.If done like this, the element among the vectorial d of current damage that obtains even bigger negative value can occur, just negative damage, this obviously is irrational.Therefore acceptable the separating (promptly have reasonable error, but can determine the position and the degree of injury thereof of damaged cable more accurately from cable system) that obtains the vectorial d of current damage becomes a rational solution, and available formula (11) is expressed this method.
Abs () is the function that takes absolute value in the formula (11), and vectorial g describes the reasonable deviation that departs from ideal linearity relation (formula (7) or formula (8)), is defined by formula (12).
G=[g
1g
2G
jG
M]
T(12) g in the formula (12)
j(j=1,2,3 ...., M) maximum allowable offset of the ideal linearity relation that departs from shown in formula (7) or the formula (8) has been described.Vector g can be selected according to the error vector e tentative calculation of formula (10) definition.
At monitored amount initial vector C
o, unit damage monitored quantitative change matrix Δ C, the current numerical value of the monitored amount of actual measurement vector C and unit damage D
uWhen (setting before calculating Δ C, is scalar) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (11), obtain acceptable the separating of the vectorial d of current damage, thereby determine the position and the degree of injury of damaged cable.
Third part of the present invention: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (comprising monitored amount monitoring system, Cable Structure support coordinate monitoring system), signal picker and computing machine etc.Require to monitor in real time or quasi real time each monitored amount, require to monitor in real time or quasi real time each Cable Structure support coordinate.
Software should the following function of tool: the data in real time that software section at first transmits according to monitoring system or quasi real time analyze and obtain current cable structure actual measurement support coordinate vector U
t, the current numerical value of monitored amount vector C, the Mechanics Calculation benchmark model A of the Cable Structure of reading pre-stored then
o, initial Cable Structure support coordinate vector U
o, current Mechanics Calculation benchmark model A
t o, current cable structural bearings coordinate vector U
t o, cable system unit damage monitored quantitative change matrix Δ C, monitored amount initial vector C
oWith unit damage value D
u, relatively the current cable structure is surveyed support coordinate vector U
tWith current cable structural bearings coordinate vector U
t o, work as U
tAnd U
t oWhen identical, find the solution formula (11) according to suitable algorithm (for example multi-objective optimization algorithm), obtain the noninferior solution of the vectorial d of current damage of cable system, just have reasonable error but can be more exactly from cable system, determine the position of damaged cable and separating of degree of injury thereof; Work as U
tAnd U
t oWhen inequality, earlier to current Mechanics Calculation benchmark model A
t oWith current cable structural bearings coordinate vector U
t oUpgrade, again at new A
t oThe basis on upgrade Δ C according to aforementioned " setting up and upgrade the method for Cable Structure unit damage monitored quantitative change matrix Δ C ", find the solution formula (11) according to suitable algorithm (for example multi-objective optimization algorithm) equally, obtain the noninferior solution of the vectorial d of current damage of cable system, just have reasonable error but can be more exactly from cable system, determine the position of damaged cable and separating of degree of injury thereof.
The inventive method specifically comprises:
A. establish total N root rope, at first determine 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;
The support cable with monitored Suo Li of appointment when b. determining hybrid monitoring is established total N root rope in the cable system, and the monitored rope force data of structure is by M on the structure
1The M of individual appointment rope
1Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope; Each total M
1Individual cable force measurement value or calculated value characterize the rope force information of structure; M
1Be one and be not less than 0 integer; The measured point with monitored strain of appointment when determining hybrid monitoring, the monitored strain data of structure can be by K on the structure
2L individual specified point, that reach each specified point
2The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
2The variation of the tested strain of all of individual specified point; Each total M
2Individual strain measurement value or calculated value characterize structural strain, M
2Be K
2And L
2Long-pending; M
2Be to be not less than 0 integer; The measured point with monitored angle of appointment when determining hybrid monitoring, the monitored angle-data of structure is by K on the structure
3L individual specified point, that cross each specified point
3H individual appointment straight line, each appointment straight line
3Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total M
3Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M
3Be K
3, L
3And H
3Long-pending; M
3Be one and be not less than 0 integer; When determining hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by K on the structure
4L individual specified point, that reach each specified point
4The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
4The variation of all coordinate components of individual specified point; Each total M
4Individual measurement of coordinates value or calculated value characterize planform, M
4Be K
4And K
4Long-pending; M
4Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of total, M is M
1, M
2, M
3And M
4Sum, definition parameter K, K is M
1, K
2, K
3And K
4Sum, K and M must not be less than the quantity N of rope; Because M monitored amount is dissimilar, so the present invention is called " health monitor method based on the cable system of hybrid monitoring when support settlement is arranged "; For simplicity, in the present invention that this step is listed " all monitored parameters of structure during hybrid monitoring " abbreviate " monitored amount " as;
C. directly measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial vector C
oObtain monitored amount initial vector C in actual measurement
oThe time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure support coordinate data of structure, and initial Cable Structure support coordinate data are formed initial Cable Structure support coordinate vector U
o
D. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data and the initial Cable Structure support coordinate vector U of rope
oSet up the initial mechanical calculating benchmark model A of Cable Structure
oAnd set up for the first time the current Mechanics Calculation benchmark model A of Cable Structure
t o, the measured data of Cable Structure comprises the initial geometric data of the initial rope force data of all ropes of Cable Structure, initial Cable Structure support coordinate data and Cable Structure at least; Set up for the first time the current Mechanics Calculation benchmark model A of Cable Structure
t oThe time, the current Mechanics Calculation benchmark model A of Cable Structure
t oJust equal the initial mechanical calculating benchmark model A of Cable Structure
oCurrent Mechanics Calculation benchmark model A corresponding to Cable Structure
t oCable Structure support coordinate data form current cable structural bearings coordinate vector U
t o, set up the current Mechanics Calculation benchmark model A of Cable Structure the first time
t oThe time, U
t oJust equal U
o
E. go on foot the k circulation in step from entering here by e; In structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, and all Cable Structure support coordinate current datas are formed current cable structure actual measurement support coordinate vector U
t
F. according to current cable structure actual measurement support coordinate vector U
t, upgrade current Mechanics Calculation benchmark model A where necessary
t oWith current cable structural bearings coordinate vector U
t o
G. at current Mechanics Calculation benchmark model A
t oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D
u
H. 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. define the vectorial d of the current damage of cable system, the element number of the vectorial d of the current damage of cable system equals the quantity of rope, between the element of the vectorial d of the current damage of cable system and the rope is one-to-one relationship, and the element numerical value of the vectorial d of the current damage of cable system is represented the degree of injury or the health status of corresponding rope;
J. the current numerical value vector C of the monitored amount of foundation is with monitored amount initial vector C
o, Cable Structure unit damage monitored quantitative change matrix Δ C, unit damage scalar D
uAnd the linear approximate relationship that exists between the vectorial d of the current damage of cable system to be asked, this linear approximate relationship can be expressed as formula 1, and other amount in the formula 1 except that d is known, finds the solution formula 1 and just can calculate the vectorial d of the current damage of cable system; Because the element numerical value of the vectorial d of current damage is represented the degree of injury of corresponding rope, so define the impaired and degree of injury of which rope according to current damage vector, has promptly realized the health monitoring of cable system in the Cable Structure; If the numerical value of a certain element of current cable damage vector is 0, represent that the pairing rope of this element is intact, do not damage; 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;
Formula 1
K. get back to the e step, begin to have e to go on foot the k circulation next time in step.
In step f, according to current cable structure actual measurement support coordinate vector U
t, upgrade current Mechanics Calculation benchmark model A where necessary
t oWith current cable structural bearings coordinate vector U
t oConcrete grammar be:
F1. actual measurement obtains current cable structure actual measurement support coordinate vector U in step e
tAfter, compare U
tWith in U
t oIf, U
tEqual U
t o, A then
t oAnd U
t oRemain unchanged;
F2. actual measurement obtains current cable structure actual measurement support coordinate vector U in step e
tAfter, compare U
tAnd U
t oIf, U
tBe not equal to U
t o, then need A
t oAnd U
t oUpgrade, update method is: calculate U earlier
tWith U
oPoor, U
tWith U
oDifference be exactly the current support displacement of Cable Structure bearing about initial position, V represents support displacement with the support displacement vector, between element among the support displacement vector V and the support displacement component is one-to-one relationship, the numerical value of an element is corresponding to the displacement of an assigned direction of an appointment bearing among the support displacement vector V, and wherein support displacement is exactly the support settlement amount at the component of gravity direction; Upgrade A
t oMethod be to A
oIn the Cable Structure bearing apply the constraint of current support displacement, the numerical value of current support displacement constraint is just taken from the numerical value of corresponding element among the support displacement vector V, to A
oIn the Cable Structure bearing apply the current Mechanics Calculation benchmark model A that obtains upgrading after the support displacement constraint
t o, upgrade A
t oThe time, U
t oAll elements numerical value is also with the e U in step
tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U
t o, so just obtained correctly corresponding to A
t oU
t o
In step g, at current Mechanics Calculation benchmark model A
t oThe basis on, obtain Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D by the several times Mechanics Calculation
uConcrete grammar be:
G1. Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current Mechanics Calculation benchmark model A
t oWith current cable structural bearings coordinate vector U
t oThe time, must upgrade Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D simultaneously
u
G2. at the current Mechanics Calculation benchmark model A of Cable Structure
t oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have N root rope that N calculating is just arranged, calculate each time in the hypothesis cable system and have only a rope that unit damage scalar D is arranged
uThe rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, calculate the current calculated value of all monitored amounts in the Cable Structure each time, the current calculated value of the monitored amount of all that calculate is formed a monitored amount calculation current vector each time;
G3. the monitored amount calculation current vector that calculates each time deducts monitored amount initial vector and obtains a monitored quantitative change vector; There is N root rope that N monitored quantitative change vector just arranged;
G4. form the Cable Structure unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N monitored quantitative change vector; Each row of Cable Structure unit damage monitored quantitative change matrix Δ C are corresponding to a monitored quantitative change vector.
Beneficial effect: method disclosed by the invention occurs at the Cable Structure bearing under the situation of sedimentation, many ropes synchronous when impaired very exactly monitoring and evaluation go out the health status (position and the degree of injury that comprise all damaged cables) of cable system, effective health monitoring of the cable system of system and method disclosed by the invention when having support settlement is very useful.
Embodiment
When support settlement is arranged,, the invention discloses a kind of system and method for the health status of each root rope in the cable system of can monitoring rationally and effectively in the Cable Structure at the health monitoring of the cable system of Cable Structure.The following describes of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
The present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status of the cable system of Cable Structure.During concrete enforcement, the following step is a kind of in the various steps that can take.
The first step: establish total N root rope, at first determine 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.
The support cable with monitored Suo Li of appointment when determining hybrid monitoring is established total N root rope in the cable system, and the monitored rope force data of structure is by M on the structure
1The M of individual appointment rope
1Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope.Each total M
1Individual cable force measurement value or calculated value characterize the rope force information of structure.M
1Be one and be not less than 0 integer.When reality was selected the rope of monitored Suo Li, the rope that can select those Suo Li to be easy to measure was monitored rope.
The measured point with monitored strain of appointment when determining hybrid monitoring, the monitored strain data of structure can be by K on the structure
2L individual specified point, that reach each specified point
2The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
2The variation of the tested strain of all of individual specified point.Each total M
2Individual strain measurement value or calculated value characterize structural strain, M
2Be K
2And L
2Long-pending.M
2Be one and be not less than 0 integer.A near point each fixed endpoint that can be exactly each root rope (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge) with the measured point of monitored strain, this point generally should not be a stress concentration point, avoiding occurring excessive strain measurement value, the fixed endpoint of the rope of the monitored Suo Li of appointment or in its vicinity when these points generally should all be not hybrid monitoring yet.
The measured point with monitored angle of appointment when determining hybrid monitoring, the monitored angle-data of structure is by K on the structure
3L individual specified point, that cross each specified point
3H individual appointment straight line, each appointment straight line
3Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M
3Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M
3Be K
3, L
3And H
3Long-pending.M
3Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) or a near point it of each root rope, and the point of monitored angle-data generally should all not be chosen as " fixed endpoint of the rope of the monitored Suo Li of appointment or point in its vicinity in the hybrid monitoring " and " point of the monitored strain of appointment or point in its vicinity in the hybrid monitoring "; Can only measure an angle coordinate of specifying straight line at each specified point, for example measure the body structure surface normal of specified point or the tangent line angle coordinate with respect to the acceleration of gravity direction, in fact be exactly measurement of dip angle here.
When determining hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by K on the structure
4L individual specified point, that reach each specified point
4The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
4The variation of all coordinate components of individual specified point.Each total M
4Individual measurement of coordinates value or calculated value characterize planform, M
4Be K
4And L
4Long-pending.M
4Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (for example being the stiff end of drag-line on bridge of cable-stayed bridge) of each root rope; Here Xuan Ding monitored point should all not selected " fixed endpoint of the rope of the monitored Suo Li of appointment or point in its vicinity in the hybrid monitoring ", " point of the monitored strain of appointment or point in its vicinity in the hybrid monitoring " and " point of the monitored angle-data of appointment or point in its vicinity in the hybrid monitoring " for use.
Comprehensive above-mentioned monitored amount, total be total M monitored amount with regard to hybrid monitoring, and M is M
1, M
2, M
3And M
4Sum, definition parameter K, K is M
1, K
2, K
3And K
4Sum, K and M must not be less than the quantity N of rope.Because M monitored amount is dissimilar, so the present invention is called " method of discerning the support cable that needs adjustment Suo Li based on hybrid monitoring ".For simplicity, in the present invention that this step is listed " all monitored parameters of structure during hybrid monitoring " abbreviate " monitored amount " as.
Second step: directly measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial vector C
oObtain monitored amount initial vector C in actual measurement
oThe time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data (is exactly its initial bridge type data for cable-stayed bridge) and the initial Cable Structure support coordinate data of structure, and initial Cable Structure support coordinate data are formed initial Cable Structure support coordinate vector U
o
The 3rd step: according to the Non-Destructive Testing data and the initial Cable Structure support coordinate vector U of the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure (data such as initial Suo Li, structural modal data that comprise structure original geometric form data, strain data, all ropes, to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, strain data, rope force data, bridge), rope
oSet up the initial mechanical calculating benchmark model A of Cable Structure
oAnd set up for the first time the current Mechanics Calculation benchmark model A of Cable Structure
t oSet up for the first time the current Mechanics Calculation benchmark model A of Cable Structure
t oThe time, the current Mechanics Calculation benchmark model A of Cable Structure
t oInitial mechanical calculating benchmark model A with Cable Structure
oIdentical; Current Mechanics Calculation benchmark model A corresponding to Cable Structure
t oCable Structure support coordinate data form current cable structural bearings coordinate vector U
t oSet up for the first time the current Mechanics Calculation benchmark model A of Cable Structure
t oThe time, U
t oJust equal U
oBased on initial mechanical calculating benchmark model A
oThe computational data that calculates structure must be very near its measured data, and error generally must not be greater than 5%.
The 4th step: in structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, and all data are formed current cable structure actual measurement support coordinate vector U
t
The 5th step: according to current cable structure actual measurement support coordinate vector U
t, upgrade current Mechanics Calculation benchmark model A where necessary
t oWith current cable structural bearings coordinate vector U
t oSurvey support coordinate vector U in the current cable structure that the actual measurement of the 4th step obtains
tAfter, compare U
tAnd U
t oIf, U
tEqual U
t o, A then
t oAnd U
t oRemain unchanged; If U
tBe not equal to U
t o, then need A
t oAnd U
t oUpgrade, update method is: calculate U earlier
tWith U
oPoor, U
tWith U
oDifference be exactly the current support displacement of Cable Structure bearing about initial position, V represents support displacement with the support displacement vector, between element among the support displacement vector V and the support displacement component is one-to-one relationship, the numerical value of an element is corresponding to the displacement of an assigned direction of an appointment bearing among the support displacement vector V, and wherein support displacement is exactly the support settlement amount at the component of gravity direction; Upgrade A
t oMethod be to A
oIn the Cable Structure bearing apply the constraint of current support displacement, the numerical value of current support displacement constraint is just taken from the numerical value of corresponding element among the support displacement vector V, to A
oIn the Cable Structure bearing apply the current Mechanics Calculation benchmark model A that obtains upgrading after the support displacement constraint
t o, upgrade A
t oThe time, U
t oAll elements numerical value is also with the 4th U that goes on foot
tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U
t o, so just obtained correctly corresponding to A
t oU
t o
The 6th step: at current Mechanics Calculation benchmark model A
t oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D
uConcrete grammar is: Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current Mechanics Calculation benchmark model A
t oWith current cable structural bearings coordinate vector U
t oThe time, must upgrade Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D simultaneously
uCurrent Mechanics Calculation benchmark model A in Cable Structure
t oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have N root rope that N calculating is just arranged, calculate each time in the hypothesis cable system and have only a rope that unit damage D is arranged
u(for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), the rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, calculate the current calculated value of all monitored amounts in the Cable Structure each time, the current calculated value of the monitored amount of all that calculate is formed a monitored amount calculation current vector C each time; Calculating monitored amount calculation current vector each time deducts monitored amount initial vector and obtains a monitored quantitative change vector; There is N root rope that N monitored quantitative change vector just arranged; Form the unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N monitored quantitative change vector; Each row of unit damage monitored quantitative change matrix are corresponding to a monitored quantitative change vector.
The 7th step: set up linear relationship error vector e and vectorial g.Utilize data (the monitored amount initial vector C of front
o, unit damage monitored quantitative change matrix Δ C), when the 6th step calculated each time, promptly in calculating the hypothesis cable system each time, have only a rope that unit damage D is arranged
uThe rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, calculate each time and all utilize mechanics method (for example adopting finite element method) to calculate that all specify the current numerical value of monitored amount in the cable system in the Cable Structure, when calculate forming a monitored amount calculation current vector C each time, calculate to form a vectorial d of damage each time, this damages in all elements of vectorial d has only the numerical value of an element to get D
u, the numerical value of other element gets 0, damages that numerical value is D among the vectorial d
uThe unit damage degree D of element unique damaged cable when calculating corresponding to this time
uWith C, C
o, Δ C, D
u, d brings formula (10) into, obtains a linear relationship error vector e, calculates a linear relationship error vector e each time; Have N root rope that N calculating is just arranged, N linear relationship error vector e just arranged, will obtain a vector after this N the linear relationship error vector e addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e.Vector g equals final error vector e.
The 8th 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 subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate measurement subsystem, signal conditioner etc.), Cable Structure support coordinate monitoring system (for example measuring with total powerstation), signal (data) collector, the computing machine and the panalarm of communicating by letter.The support coordinate of each monitored amount, each Cable Structure 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; When monitoring rope when damage is arranged, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 9th step: with monitored amount initial vector C
o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D
uParameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
The tenth step: the establishment and the cable system health monitoring systems software of Cable Structure when installation and operation has support settlement on computers, this software will be finished functions such as monitoring that the present invention's " cable system health monitoring of Cable Structure when support settlement is arranged " required by task wants, record, control, storage, calculating, notice, warning (all work that can finish with computing machine in this specific implementation method)
The 11 step: actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount ";
The 12 step: the current numerical value vector C of the monitored amount of foundation is with monitored amount initial vector C
o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D
uAnd the linear approximate relationship (formula (7)) that exists between the vectorial d of the current damage of cable system (forming) by all Suo Dangqian amount of damage, calculate the noninferior solution of the vectorial d of the current damage of cable system according to multi-objective optimization algorithm, just have reasonable error but can from all ropes, determine the position of damaged cable and separating of degree of injury thereof more exactly.
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, this implementation step is that example provides the process of finding the solution the vectorial d of current damage with the goal programming method only, and the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.
According to the goal programming method, formula (7) can transform the multi-objective optimization question shown in an accepted way of doing sth (13) and the formula (14), γ is a real number in the formula (13), R is a real number field, area of space Ω has limited the span (each element of present embodiment requirements vector d is not less than 0, is not more than 1) of each element of vectorial d.The meaning of formula (13) is to seek the real number γ of a minimum, makes formula (14) be met.G (d) is defined by formula (15) in the formula (14), the deviation that allows between middle G (d) of the product representation formula (14) of weighing vector W and γ and the vectorial g in the formula (14), and the definition of g is referring to formula (12), and its value calculates in the 7th step.Vector W can be identical with vectorial g during actual computation.The concrete programming of goal programming method realizes having had universal program directly to adopt.Use goal programming method just can be in the hope of the vectorial d of the current damage of cable system.
minimize γ
(13)
γ∈R,d∈Ω
G(d)-Wγ≤g (14)
The element number of the vectorial d of the current damage of cable system equals the quantity of rope, is one-to-one relationship between the element of the vectorial d of the current damage of cable system and the rope, and the element numerical value of the vectorial d of the current damage of cable system is represented the degree of injury or the health status of corresponding rope.If the numerical value of a certain element of the vectorial d of the current damage of cable system that solves is 0, represent that the pairing rope of this element is intact, do not damage; 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 13 step: the computing machine in the health monitoring systems regularly generates cable system health condition form automatically or by the personnel operation health monitoring systems.
The 14 step: under specified requirements, the automatic operation communication panalarm of the computing machine in the health monitoring systems to monitor staff, owner and (or) personnel of appointment report to the police.
The 15 step: got back to for the 4th step, beginning goes on foot the 15 circulation that goes on foot by the 4th.