Summary of the invention
Technical matters: the objective of the invention is when the Cable Structure bearing has sedimentation, at the health monitoring problem of cable system in the Cable Structure, provide a kind of based on to the hybrid monitoring of multiclass parameter, can monitor cable system in the Cable Structure rationally and effectively support settlement is arranged the time based on the progressive health monitoring method of the cable system of hybrid monitoring.
Technical scheme: the present invention is made up of the two large divisions.Be respectively: one, set up the method for required knowledge base of cable system health monitoring systems and parameter, and based on knowledge base (containing parameter) with based on the laddering health status appraisal procedure of cable system to the hybrid monitoring of multiclass parameter and actual measurement Cable Structure support coordinate; Two, the software and hardware part of health monitoring systems.
First of the present invention: set up the method for required knowledge base of cable system health monitoring systems and parameter, and based on knowledge base (containing parameter) with based on the laddering health status appraisal procedure of cable system to the hybrid monitoring of multiclass parameter and actual measurement Cable Structure support coordinate.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 damage vector d when at first needing to set up or set up this circulation beginning
_{o} ^{i}(i=1,2,3 ...), set up the initial mechanical calculating benchmark model A of Cable Structure
_{o}(for example finite element benchmark model, A in the present invention
_{o}Be 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} _{o}Bring 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 represents cycle index, i.e. the i time circulation.
If total N root rope in the cable system, the cable system initial damage vector that needs during the i time circulation beginning is designated as d
_{o} ^{i}(as the formula (1)), use d
_{o} ^{i}Cable Structure when representing this time circulation beginning is (with Mechanics Calculation benchmark model A
^{i}The health status of cable system expression).
${d}_{o}^{i}={\left[\begin{array}{cccccc}{d}_{o1}^{i}& {d}_{o2}^{i}& \·\·\·& {d}_{\mathrm{oj}}^{i}& \·\·\·& {d}_{\mathrm{oN}}^{i}\end{array}\right]}^{T}---\left(1\right)$
D in the formula (1)
^{i} _{Oj}(i=1,2,3, J=1,2,3 ...., when N) the i time circulation of expression begins, Mechanics Calculation benchmark model A
^{i}In the initial damage value of j root rope of cable system, d
^{i} _{Oj}Being to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, represents the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%.
Set up cable system initial damage vector during circulation beginning for the first time and (be designated as d according to formula (1)
^{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 damage vector d
^{1} _{o}If when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, perhaps can think when the structure original state is the not damaged state vectorial d
^{1} _{o}Each element numerical value get 0.
The i time (i=2,3,4,5,6 ...) the cable system initial damage vector d of needs when circulation begins
^{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 described below.
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}
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, rope force data, draw-bar pull data, Cable Structure support 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) and design drawing, as-constructed drawing, utilize mechanics method (for example finite element method) to set up A
_{o}If 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
_{o}No matter which kind of method to obtain A with
_{o}, based on A
_{o}The 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
_{o}Strain 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
_{o}Cable Structure support coordinate data form initial Cable Structure support coordinate vector U
_{o}A
_{o}And U
_{o}Be 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
^{1}Just equal A
_{o}
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 described below.
Existing Mechanics Calculation benchmark model A
^{1}With cable system initial damage vector d
^{1} _{o}After, model A
^{1}In the damage of each rope by vectorial d
^{1} _{o}Express.At A
^{1}The basis on, the damage of all ropes is changed to 0, mechanical model A
^{1}The damage that is updated to all ropes all is that 0 mechanical model (is designated as A
^{0}), mechanical model A
^{0}Be actually the mechanical model of intact Cable Structure correspondence.Might as well claim model A
^{0}Not damaged model A for Cable Structure
^{0}
Monitored 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
_{1}The M of individual appointment rope
_{1}Individual 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
_{1}Individual cable force measurement value or calculated value characterize the rope force information of structure.M
_{1}Be one and be not less than 0 integer.
The monitored strain data of structure can be by K on the structure
_{2}L individual specified point, that reach each specified point
_{2}The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
_{2}The 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
_{2}Be one and be not less than 0 integer.
The monitored angle-data of structure is by K on the structure
_{3}L individual specified point, that cross each specified point
_{3}H individual appointment straight line, each appointment straight line
_{3}Individual 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
_{3}Be one and be not less than 0 integer.
The monitored shape data of structure is by K on the structure
_{4}L individual specified point, that reach each specified point
_{4}The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
_{4}The 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
_{4}Be 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 ".
The present invention " initial value vector C of monitored amount
^{i} _{o}" (i=1,2,3 ...) initial value (referring to formula (2)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation of expression begins, C
^{i} _{o}Full name be the initial value vector of monitored amount " the i time circulation ".
${C}_{o}^{i}={\left[\begin{array}{cccccc}{C}_{o1}^{i}& {C}_{o2}^{i}& \·\·\·& {C}_{\mathrm{ok}}^{i}& \·\·\·& {C}_{\mathrm{oM}}^{i}\end{array}\right]}^{T}---\left(2\right)$
C in the formula (2)
^{i} _{Ok}(i=1,2,3, K=1,2,3 ...., M; M 〉=N; ) k monitored amount when being the i time circulation beginning, in the Cable Structure.Vector C
^{i} _{o}Be to be formed according to certain series arrangement by the monitored amount of previously defined M, this is put in order there is no 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 C of the 1st the monitored amount that circulates
^{1} _{o}" (seeing formula (2)) be made up of measured data, because according to model A
^{1}The 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 " the vectorial C of the initial value of the i time monitored amount of circulation that needs
^{i} _{o}", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) calculate before the loop ends and to obtain, concrete grammar is described below.
Second step: in Cable Structure military service process, in circulation each time, constantly (all data are formed current cable structure actual measurement support coordinate vector U to actual measurement acquisition Cable Structure support coordinate current data
^{Ti}, vectorial U
^{Ti}Definition mode and vectorial U
_{o}Identical).For simplicity, for the i time circulation, the Cable Structure support coordinate current data the when last time is upgraded current Mechanics Calculation benchmark model is designated as current cable structural bearings coordinate vector U
^{Ti} _{o}Set up and renewal A
^{Ti} _{o}Method be: in round-robin each time zero hour, the current Mechanics Calculation benchmark model A of Cable Structure
^{Ti} _{o}Just equal A
^{i}(i=1,2,3,4,5,6 ...).In 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
^{Ti}If, U
^{Ti}Equal U
^{Ti} _{o}, then do not need A
^{Ti} _{o}Upgrade; If U
^{Ti}Be not equal to U
^{Ti} _{o}, then need A
^{Ti} _{o}Upgrade, at this moment U
^{Ti}With U
_{o}Difference 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
^{Ti} _{o}Method be: at A
_{o}The basis on make that the health status of rope is cable system initial damage vector d
^{i} _{o}, more further to A
_{o}In 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 current support displacement vector V, to A
_{o}In the Cable Structure bearing apply the constraint of current support displacement after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
^{Ti} _{o}, upgrade A
^{Ti} _{o}The time, U
^{Ti} _{o}All 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 3rd step: circulation each time needs to set up " unit damage monitored quantitative change matrix " and " nominal unit damage vector ", and " unit damage monitored quantitative change matrix " that the i time circulation set up is designated as Δ C
^{i}(i=1,2,3 ...)." nominal unit damage vector " that the i time circulation set up is designated as D
^{i} _{u}Δ C in circulation each time
^{i}And D
^{i} _{u}Bring in constant renewal in, promptly upgrading current Mechanics Calculation benchmark model A
^{Ti} _{o}The time, upgrade Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector " be designated as D
^{i} _{u}
Set up and renewal Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector " be designated as D
^{i} _{u}Process as follows:
Current Mechanics Calculation benchmark model A in Cable Structure
^{Ti} _{o}The 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 to increase unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again on the basis of original damage (original damage can be 0, can not be 0 also).Calculate for convenient, when setting 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 unit damage (in subsequent step, damage numerical value that calculate, rope---be called name damage 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, the name that the formula that therefore must foundation hereinafter provides will calculate is damaged and is converted into true damage.)。The rope of appearance damage was different from the rope that appearance damages in other time calculating during a round-robin calculated each time together, and supposed that each time the unit damage value of the rope that damage is arranged can be different from the unit damage value of other ropes, with " the vectorial D of nominal unit damage
^{i} _{u}" (as the formula (3)) write down the unit damage of the supposition of all ropes in each time circulation, circulation time is designated as D for the first time
^{1} _{u}Calculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, the current calculated value of the M of appointment monitored amount in front, the current calculated value that calculates gained M monitored amount is each time formed one " the current numerical value vector of the calculating of monitored amount ", and (when hypothesis j root rope had unit damage, available formula (4) was represented the current numerical value vector of the calculating C of M monitored amount of all appointments
^{1} _{Tj}); The current numerical value vector of the calculating of the monitored amount that calculates each time deducts the initial value vector C of monitored amount
^{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 unit damage is arranged or numbering etc.) (when j root rope has unit damage, used δ C under this condition
^{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 vectorial D for formula (4)
^{1} _{u}J element D
_{Uj}Gained), the numerical value change of monitored amount vector δ C
^{1} _{j}Each element representation since when calculating supposition the unit damage (D for example of the Na Gensuo (for example j root rope) of unit damage is arranged
_{Uj}), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage D of supposition
_{Uj}Rate of change; There is N root rope that N " the numerical value change vector of monitored amount " just arranged, the numerical value change vector of each monitored amount has M (general, the individual element of M 〉=N) is formed " the unit damage monitored quantitative change matrix Δ C that M * N element arranged successively by this N " the numerical value change vector of monitored amount "
^{1}" (the capable N row of M), each vectorial δ C
^{1} _{j}(j=1,2,3 ...., N) be matrix Δ C
^{1}One row, Δ C
^{1}Definition as the formula (8).
${D}_{u}^{i}={\left[\begin{array}{cccccc}{D}_{u1}^{i}& {D}_{u2}^{i}& \·\·\·& {D}_{\mathrm{uj}}^{i}& \·\·\·& {D}_{\mathrm{uN}}^{i}\end{array}\right]}^{T}---\left(3\right)$
Nominal unit damage vector D in the formula (3)
^{i} _{u}Element D
^{i} _{Uj}(i=1,2,3, J=1,2,3 ...., N) the unit damage numerical value of the j root rope of supposition in the i time circulation of expression, vectorial D
^{i} _{u}In the numerical value of each element can be the same or different.
${C}_{\mathrm{tj}}^{i}={\left[\begin{array}{cccccc}{C}_{\mathrm{tk}1}^{i}& {C}_{\mathrm{tk}2}^{i}& \·\·\·& {C}_{\mathrm{tjk}}^{i}& \·\·\·& {C}_{\mathrm{tjM}}^{i}\end{array}\right]}^{T}---\left(4\right)$
Elements C in the formula (4)
^{i} _{Tjk}(i=1,2,3 ...; J=1,2,3 ...., N; K=1,2,3 ...., M; When the i time circulation of the expression of M 〉=N) has unit damage owing to j root rope, according to the current numerical value of calculating of the monitored amount of pairing k the appointment of coding rule.
$\mathrm{\δ}{C}_{j}^{i}=\frac{{C}_{\mathrm{tj}}^{i}-{C}_{o}^{i}}{{D}_{\mathrm{uj}}^{i}}---\left(5\right)$
The subscript i of each amount in the formula (5) (i=1,2,3 ...) the i time circulation of expression, subscript j (j=1,2,3 ...., N) expression j root rope has unit damage, D in the formula
^{i} _{Uj}Be vectorial D
^{i} _{u}In j element.Vector δ C
^{i} _{j}Definition as the formula (6), δ C
^{i} _{j}K (k=1,2,3 ...., M; The individual element δ C of M 〉=N)
^{i} _{Jk}Represent to set up matrix Δ C in the i time circulation
^{i}The time, suppose that the change amount of calculating a gained k monitored amount when j root rope has unit damage is with respect to the unit damage D that supposes
^{i} _{Uj}Rate of change, it defines as the formula (7).
${\mathrm{\δC}}_{j}^{i}={\left[\begin{array}{cccccc}{\mathrm{\δC}}_{j1}^{i}& \mathrm{\δ}{C}_{j2}^{i}& \·\·\·& {\mathrm{\δC}}_{\mathrm{jk}}^{i}& \·\·\·& \mathrm{\δ}{C}_{\mathrm{jM}}^{i}\end{array}\right]}^{T}---\left(6\right)$
$\mathrm{\δ}{C}_{\mathrm{jk}}^{i}=\frac{{C}_{\mathrm{tjk}}^{i}-{C}_{\mathrm{ok}}^{i}}{{D}_{\mathrm{uj}}^{i}}---\left(7\right)$
The definition of each amount has been previously described in the formula (7).
${\mathrm{\ΔC}}^{i}=\left[\begin{array}{cccccc}\mathrm{\δ}{C}_{1}^{i}& \mathrm{\δ}{C}_{2}^{i}& \·\·\·& \mathrm{\δ}{C}_{j}^{i}& \·\·\·& \mathrm{\δ}{C}_{N}^{i}\end{array}\right]---\left(8\right)$
Vectorial δ C in the formula (8)
^{i} _{j}(i=1,2,3 ....,, j=1,2,3 ...., N) in the i time circulation of expression, because j root rope has unit damage D
^{i} _{Uj}Cause, the relative value of all monitored amounts changes.Matrix Δ C
^{i}The coding rule of row (subscript j) and front vector d
^{i} _{o}The coding rule of subscript j of element identical.
In Cable Structure military service process, in circulation each time, constantly actual measurement obtains Cable Structure support coordinate current data, in case monitor U
^{Ti}Be not equal to U
^{Ti} _{o}, then needed to get back to second step to A
^{Ti} _{o}Upgrade, to A
^{Ti} _{o}Enter this step after upgrading again to Δ C
^{i}Upgrade.In fact each time the circulation in Δ C
^{i}Bring in constant renewal in, promptly upgrading current Mechanics Calculation benchmark model A
^{Ti} _{o}Afterwards, upgrade Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}
The 4th step: the current health status of identification cable system.Detailed process is as follows.
I (i=1,2,3 ...) in the inferior circulation, cable system " current (calculating or actual measurement) numerical value vector C of monitored amount
^{i}" " the initial value of monitored amount vector C together
^{i} _{o}", " unit damage monitored quantitative change matrix Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between linear approximate relationship, shown in (9) or formula (10).
${C}^{i}={C}_{o}^{i}+\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}---\left(9\right)$
${C}^{i}-{C}_{o}^{i}=\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}---\left(10\right)$
Current (calculating or actual measurement) numerical value vector C of monitored amount in formula (9) and the formula (10)
^{i}Definition be similar to the initial value vector C of monitored amount
^{i} _{o}Definition, see formula (11); The vectorial d of the current name damage of cable system
^{i} _{c}Definition see formula (12).
${C}^{i}={\left[\begin{array}{cccccc}{C}_{1}^{i}& {C}_{2}^{i}& \·\·\·& {C}_{k}^{i}& \·\·\·& {C}_{M}^{i}\end{array}\right]}^{T}---\left(11\right)$
Elements C in the formula (11)
^{i} _{k}(i=1,2,3 ....; K=1,2,3 ...., M; M 〉=N) be the i time circulation time Cable Structure, according to the current numerical value of the monitored amount of the pairing k of being numbered of coding rule.
${d}_{c}^{i}={\left[\begin{array}{cccccc}{d}_{c1}^{i}& {d}_{c2}^{i}& \·\·\·& {d}_{\mathrm{cj}}^{i}& \·\·\·& {d}_{\mathrm{cN}}^{i}\end{array}\right]}^{T}---\left(12\right)$
D in the formula (12)
^{i} _{Cj}(i=1,2,3 ....; J=1,2,3 ...., N) be the current nominal impairment value of cable system j root rope in the i time circulation, vectorial d
^{i} _{c}The coding rule and the matrix Δ C of subscript j of element
^{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).
${e}^{i}=\mathrm{abs}(\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}-{C}^{i}+{C}_{o}^{i})---\left(13\right)$
Abs () is the function that takes absolute value in the formula (13), and 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 C of monitored amount
^{i}" directly find the solution and obtain the vectorial d of Suo Dangqian name damage
^{i} _{c}If done like this, the damage vector d that obtains
^{i} _{c}In element in addition bigger negative value can appear, just negative damage, this obviously is irrational.Therefore obtaining rope damages vectorial d
^{i} _{c}Acceptable separating (promptly have reasonable error, but can be more accurately from cable system, determine the position and the degree of injury thereof of damaged cable) become a rational solution, available formula (14) is expressed this method.
$\mathrm{abs}(\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}-{C}^{i}+{C}_{o}^{i})\≤{g}^{i}---\left(14\right)$
Abs () is the function that takes absolute value in the formula (14), vectorial g
^{i}Description departs from the reasonable deviation of ideal linearity relation (formula (9) or formula (10)), is defined by formula (15).
${g}^{i}={\left[\begin{array}{cccccc}{g}_{1}^{i}& {g}_{2}^{i}& \·\·\·& {g}_{k}^{i}& \·\·\·& {g}_{M}^{i}\end{array}\right]}^{T}---\left(15\right)$
G in the formula (15)
^{i} _{k}(i=1,2,3 ....; K=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (9) or the formula (10) in the i time circulation has been described.Vector g
^{i}Can be according to the error vector e of formula (13) definition
^{i}Tentative calculation is selected.
Initial value vector C in monitored amount
^{i} _{o}(survey or calculate), Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}The current numerical value vector C of (calculating) and monitored amount
^{i}When (actual measurement obtains) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain the vectorial d of the current name damage of cable system
^{i} _{c}Acceptable separating, the current actual damage of cable system vector d
^{i}The element of (formula (16) is seen in definition) can calculate according to formula (17), has just obtained Suo Dangqian actual damage vector d
^{i}Thereby, can be by d
^{i}Determine to decrease the position and the degree of injury of rope, just realized the health monitoring of cable system.
${d}^{i}={\left[\begin{array}{cccccc}{d}_{1}^{i}& {d}_{2}^{i}& \·\·\·& {d}_{j}^{i}& \·\·\·& {d}_{N}^{i}\end{array}\right]}^{T}---\left(16\right)$
D in the formula (16)
^{i} _{j}(j=1,2,3, J=1,2,3 ...., N) the actual damage value of j root rope in the i time circulation of expression, formula (17), d are seen in its definition
^{i} _{j}Being to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, represents the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, vectorial d
^{i}The coding rule of element and formula (1) in vectorial d
^{i} _{o}The coding rule of element identical.
${d}_{j}^{i}=1-(1-{d}_{\mathrm{oj}}^{i})(1-{d}_{\mathrm{cj}}^{i})---\left(17\right)$
D in the formula (17)
^{i} _{Oj}(i=1,2,3,4, J=1,2,3 ...., N) be vectorial d
^{i} _{o}J element, d
^{i}c
_{j}Be vectorial d
^{i} _{c}J element.
The 5th step: judge whether to finish this (the i time) circulation, if, then finish the preceding tailing in work of 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 as follows.
In this (the i time) circulation, try to achieve the vectorial d of current name damage
^{i} _{c}After, at first, set up mark vector F according to formula (18)
^{i}, formula (19) has provided mark vector F
^{i}The definition of j element; If mark vector F
^{i}Element be 0 entirely, then in this circulation, continue health monitoring and calculating to cable system; If mark vector F
^{i}Element be not 0 entirely, then finish subsequent step after, enter next time circulation.So-called subsequent step is: at first, according to formula (20) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation
^{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
^{0}The basis on, the health status situation that makes rope is 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+1}At last, by to Mechanics Calculation benchmark model A
^{I+1}The 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 C of monitored amount of circulation
^{I-1} _{o}" (i=1,2,3,4 ...).
${F}^{i}={\left[\begin{array}{cccccc}{F}_{1}^{i}& {F}_{2}^{i}& \·\·\·& {F}_{j}^{i}& \·\·\·& {F}_{N}^{i}\end{array}\right]}^{T}---\left(18\right)$
Mark vector F in the formula (18)
^{i}Subscript i represent the i time the circulation, its element F
^{i} _{j}(j=1,2,3 ..., subscript j N) represents the damage characteristic of j root rope, can only get 0 and 1 two amount, concrete value rule is seen formula (19).
${F}_{j}^{i}=\left\{\begin{array}{ccc}0,& \mathrm{if}& {d}_{\mathrm{cj}}^{i}<{D}_{\mathrm{uj}}^{i}\\ 1,& \mathrm{if}& {d}_{\mathrm{cj}}^{i}\≥{D}_{\mathrm{uj}}^{i}\end{array}\right.---\left(19\right)$
Element F in the formula (19)
^{i} _{j}Be mark vector F
^{i}J element, D
^{i} _{Uj}Be nominal unit damage vector D
^{i} _{u}J element (seeing formula (3)), d
^{i} _{Cj}Be the vectorial d of the current name damage of cable system
^{i} _{c}J element (seeing formula (12)), they all represent the relevant information of j root rope.
${d}_{\mathrm{oj}}^{i+1}=1-(1-{d}_{\mathrm{oj}}^{i})(1-{D}_{\mathrm{uj}}^{i}{F}_{j}^{i})---\left(20\right)$
D in the formula (20)
^{i} _{Uj}Be nominal unit damage vector D
^{i} _{u}J element (seeing formula (3)), d
^{i} _{Cj}Be the vectorial d of the current name damage of cable system
^{i} _{c}J element (seeing formula (12)).
Second portion 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 section should be finished the process that first of the present invention sets, promptly finish needed among the present invention, can be with functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
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
_{1}The M of individual appointment rope
_{1}Individual 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
_{1}Individual cable force measurement value or calculated value characterize the rope force information of structure; M
_{1}Be 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
_{2}L individual specified point, that reach each specified point
_{2}The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
_{2}The variation of the tested strain of all of individual specified point; Each total M
_{2}Individual strain measurement value or calculated value characterize structural strain, M
_{2}Be K
_{2}And L
_{2}Long-pending; M
_{2}Be 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
_{3}L individual specified point, that cross each specified point
_{3}H individual appointment straight line, each appointment straight line
_{3}Individual 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}Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M
_{3}Be K
_{3}, L
_{3}And H
_{3}Long-pending; M
_{3}Be 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
_{4}L individual specified point, that reach each specified point
_{4}The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
_{4}The variation of all coordinate components of individual specified point; Each total M
_{4}Individual measurement of coordinates value or calculated value characterize planform, M
_{4}Be K
_{4}And L
_{4}Long-pending; M
_{4}Be 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
_{3}And M
_{4}Sum, definition parameter K, K is M
_{1}, K
_{2}, K
_{3}And K
_{4}Sum, 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. 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 damage vector d
^{1} _{o}If when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, perhaps can think when the structure original state is the not damaged state vectorial d
^{1} _{o}Each element numerical value get 0;
D. setting up cable system initial damage vector d
^{1} _{o}The time, directly measure the monitored amount of all appointments that calculate Cable Structure, form " the initial value vector C of monitored amount
^{1} _{o}";
E. setting up cable system initial damage vector d
^{1} _{o}Initial value vector C with monitored amount
^{1} _{o}The time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, and actual measurement obtains the initial geometric data of Cable Structure;
F. set up the initial mechanical calculating benchmark model A of Cable Structure
_{o}, set up initial Cable Structure support 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
^{1}The 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 support coordinate data, Cable Structure modal data, 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 mechanical calculating benchmark model A of Cable Structure
_{o}If 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 mechanical calculating benchmark model A of Cable Structure equally
_{o}No matter which kind of method to obtain A with
_{o}, based on A
_{o}The 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
_{o}Cable Structure support coordinate data form initial Cable Structure support coordinate vector U
_{o}A
_{o}And U
_{o}Be 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 represents 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 represents 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
^{1}Just equal A
_{o}For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A
^{Ti} _{o}", A in circulation each time
^{Ti} _{o}Can bring in constant renewal in as required, during circulation beginning each time, A
^{Ti} _{o}Equal A
^{i}Equally for sake of convenience, name " Cable Structure actual measurement support coordinate vector U
^{Ti}", in circulation each time, 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
^{Ti}, vectorial U
^{Ti}Element and vectorial U
_{o}The coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, for the i time circulation, the last time is upgraded A
^{Ti} _{o}The time Cable Structure support coordinate current data be designated as current cable structural bearings coordinate vector U
^{Ti} _{o}During circulation beginning for the first time, A
^{T1} _{o}Equal A
^{1}, U
^{T1} _{o}Equal U
_{o}
When g. circulation begins each time, make A
^{Ti} _{o}Equal A
^{i}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
^{Ti}, according to current cable structure actual measurement support coordinate vector U
^{Ti}, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
^{Ti} _{o}With current cable structural bearings coordinate vector U
^{Ti} _{o}
H. at the current Mechanics Calculation benchmark model of Cable Structure A
^{Ti} _{o}The basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector D
^{i} _{u}
I. 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}".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 the element each vector, that numbering is identical that occurs before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
J. define the vectorial d of the current name damage of cable system
^{i} _{c}With current actual damage vector d
^{i}, the element number of damage vector equals the quantity of rope, is one-to-one relationship between the element of damage vector and the rope, and the element numerical value of damage vector is represented the degree of injury 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}", " unit damage monitored quantitative change matrix Δ C
^{i}" and " the vectorial d of current name damage
^{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} _{c}Other outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage
^{i} _{c}
${C}^{i}={C}_{o}^{i}+\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}$ Formula 1
1. the current actual damage vector d that utilizes formula 2 to express
^{i}With initial damage vector d
^{i} _{o}With the vectorial d of current name damage
^{i} _{c}Element between relation, calculate current actual damage vector d
^{i}All elements.
${d}_{j}^{i}=1-(1-{d}_{\mathrm{oj}}^{i})(1-{d}_{\mathrm{cj}}^{i})$ Formula 2
J=1 in the formula 2,2,3 ..., N.
Because current actual damage vector d
^{i}Element numerical value represent the degree of injury of corresponding rope, so according to current actual damage vector d
^{i}Just can define the impaired and degree of injury of which rope, promptly realize the health monitoring of cable system in the Cable Structure; If the numerical value of a certain element of current actual 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.
M. try to achieve the vectorial d of current name damage
^{i} _{c}After, set up mark vector F according to formula 3
^{i}, formula 4 has provided mark vector f
^{i}The definition of j element;
${F}^{i}={\left[\begin{array}{cccccc}{F}_{1}^{i}& {F}_{2}^{i}& \·\·\·& {F}_{j}^{i}& \·\·\·& {F}_{N}^{i}\end{array}\right]}^{T}$ Formula 3
${F}_{j}^{i}=\left\{\begin{array}{ccc}0,& \mathrm{if}& {d}_{\mathrm{cj}}^{i}<{D}_{\mathrm{uj}}^{i}\\ 1,& \mathrm{if}& {d}_{\mathrm{cj}}^{i}\≥{D}_{\mathrm{uj}}^{i}\end{array}\right.$ Formula 4
Element F in the formula 4
^{i} _{j}Be mark vector F
^{i}J element, D
^{i} _{Uj}Be nominal unit damage vector D
^{i} _{u}J element, d
^{i} _{Cj}Be the vectorial d of the current name damage of cable system
^{i} _{c}J element, they all represent the relevant information of j root rope.J=1 in the formula 4,2,3 ..., N.
If mark vector F n.
^{i}Element be 0 entirely, then get back to step g and continue this circulation; If mark vector F
^{i}Element be not 0 entirely, then enter next step, be step o.
O. according to formula 5 calculate next time, i.e. the i+1 time required initial damage of circulation vector d
^{I+1} _{o}Each element d
^{I+1} _{Oj}
${d}_{\mathrm{oj}}^{i+1}=1-(1-{d}_{\mathrm{oj}}^{i})(1-{D}_{\mathrm{uj}}^{i}{F}_{j}^{i})$ Formula 5
D in the formula 5
^{i} _{Uj}Be nominal unit damage vector D
^{i} _{u}J element, d
^{i} _{Cj}Be the vectorial d of the current name damage of cable system
^{i} _{c}J element, F
^{i} _{j}Be mark vector F
^{i}J element.J=1 in the formula 5,2,3 ..., N.
P. at the current Mechanics Calculation benchmark model of Cable Structure A
^{Ti} _{o}The basis on, the health status that makes rope is 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}, promptly the Mechanics Calculation benchmark model is upgraded;
Q. pass through Mechanics Calculation benchmark model A
^{I+1}Calculate corresponding to model A
^{I+1}The 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} _{o}
R. 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} _{o}Equal A
^{I+1}
S. set up the required current cable structural bearings coordinate vector U that next time, promptly circulates for the i+1 time
^{Ti+1} _{o}, promptly get U
^{Ti+1} _{o}Equal U
^{Ti} _{o}
T. get back to step g, beginning is circulation next time.
In step g, according to current cable structure actual measurement support coordinate vector U
^{Ti}, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
^{Ti} _{o}With current cable structural bearings coordinate vector U
^{Ti} _{o}Concrete grammar be:
G1. actual measurement obtains current cable structure actual measurement support coordinate vector U
^{Ti}After, compare U
^{Ti}And U
^{Ti} _{o}If, U
^{Ti}Equal u
^{Ti} _{o}, then do not need A
^{Ti} _{o}Upgrade;
G2. actual measurement obtains current cable structure actual measurement support coordinate vector U
^{Ti}After, compare U
^{Ti}And U
^{Ti} _{o}If, U
^{Ti}Be not equal to U
^{Ti} _{o}, then need A
^{Ti} _{o}Upgrade, update method is: calculate U earlier
^{Ti}With U
_{o}Poor, U
^{Ti}With U
_{o}Difference be exactly that the current cable structural bearings is about setting up A
_{o}The time the current support displacement of Cable Structure bearing, represent support displacement with current support displacement vector V, between element among the current 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 current support displacement vector V, and wherein support displacement is exactly the support settlement amount at the component of gravity direction; Upgrade A
^{Ti} _{o}Method be: at A
_{o}The basis on make that the health status of rope is cable system initial damage vector d
^{i} _{o}, more further to A
_{o}In 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 current support displacement vector V, to A
_{o}In the Cable Structure bearing apply the constraint of current support displacement after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
^{Ti} _{o}, upgrade A
^{Ti} _{o}The time, U
^{Ti} _{o}All 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} _{o}The basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector D
^{i} _{u}Concrete grammar be:
H1. the i time when beginning circulation, directly h2 obtains Cable Structure unit damage monitored quantitative change matrix Δ C to the listed method of step h4 set by step
^{i}With nominal unit damage vector D
^{i} _{u}At other constantly, in step g to A
^{Ti} _{o}After upgrading, h2 is to the listed method acquisition of step h4 Cable Structure unit damage monitored quantitative change matrix Δ C set by step
^{i}With nominal unit damage vector D
^{i} _{u}If, in step g not to A
^{Ti} _{o}Upgrade, 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} _{o}The 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 to increase unit damage again on the basis of original damage, the rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, and supposition each time has the unit damage value of the rope of damage can be different from the unit damage value of other ropes, uses " nominal unit damage vector D
^{i} _{u}" write down the unit damage of the supposition of all ropes, calculate each time that all specify the current numerical value of monitored amount in the Cable Structure, the current numerical value of the monitored amount of all that calculate is formed one " the current numerical value vector of the calculating of monitored amount " each time; When hypothesis j root rope has unit damage, available C
^{i} _{Tj}" the current numerical value vector of the calculating of monitored amount " that expression is corresponding; When giving each vectorial element numbering in this step, should use same coding rule with other vector among the present invention, can guarantee any one 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 numerical value of the calculating of monitored amount vector 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 " numerical value change vector δ a C of monitored amount after the unit damage value of supposition
^{i} _{j}"; There is N root rope that N " the numerical value change vector of monitored amount " just arranged;
H4. form " the unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount "
^{i}"; " unit damage monitored quantitative change matrix Δ C
^{i}" each row corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " unit damage monitored quantitative change matrix " and the vectorial d of current name damage
^{i} _{c}With current actual damage vector d
^{i}The element coding rule identical.
Beneficial effect: system and method disclosed by the invention occurs under the situation of sedimentation at the Cable Structure bearing, having under the synchronously impaired condition of more rope monitoring and evaluation very exactly go out the health status (position and the degree of injury that comprise all damaged cables) of cable system.This is because " the current numerical value vector C of monitored amount
^{i}" " the initial value of monitored amount vector C together
^{i} _{o}", " unit damage monitored quantitative change matrix Δ C
^{i}" and " the vectorial d of current name damage
^{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; for overcoming this obstacle, the invention discloses under a kind of situation that occurs sedimentation at the Cable Structure bearing, in the minizone, approach the health monitor method of this nonlinear relationship with linear relationship.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, and the health status of the cable system that obtains according to its judgement also is reliable.Therefore, occur at the Cable Structure bearing under the situation of sedimentation, system and method disclosed by the invention is very useful to effective health monitoring of cable system.
Embodiment
When support settlement is arranged,, 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 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.
Occur at the Cable Structure bearing under the situation of sedimentation, the algorithm that the present invention adopts 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: determine type, position and the quantity of monitored amount, and numbering.Detailed process is:
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.
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
_{1}The M of individual appointment rope
_{1}Individual 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
_{1}Individual cable force measurement value or calculated value characterize the rope force information of structure.M
_{1}Be 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
_{2}L individual specified point, that reach each specified point
_{2}The strain of individual assigned direction is described, and the variation of structural strain data is exactly K
_{2}The variation of the tested strain of all of individual specified point.Each total M
_{2}Individual strain measurement value or calculated value characterize structural strain, M
_{2}Be K
_{2}And L
_{2}Long-pending.M
_{2}Be 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
_{3}L individual specified point, that cross each specified point
_{3}H individual appointment straight line, each appointment straight line
_{3}Individual 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}Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M
_{3}Be K
_{3}, L
_{3}And H
_{3}Long-pending.M
_{3}Be 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
_{4}L individual specified point, that reach each specified point
_{4}The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K
_{4}The variation of all coordinate components of individual specified point.Each total M
_{4}Individual measurement of coordinates value or calculated value characterize planform, M
_{4}Be K
_{4}And L
_{4}Long-pending.M
_{4}Be 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
_{3}And M
_{4}Sum, definition parameter K, K is M
_{1}, K
_{2}, K
_{3}And K
_{4}Sum, 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: 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 damage vector d
^{1} _{o}If when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, perhaps can think when the structure original state is the not damaged state vectorial d
^{1} _{o}Each element numerical value get 0.
The 3rd step: setting up cable system initial damage vector d
^{1} _{o}The time, directly measure the monitored amount of all appointments that calculate Cable Structure, form " the initial value vector C of monitored amount
^{1} _{o}".
The 4th step: setting up cable system initial damage vector d
^{1} _{o}Initial value vector C with monitored amount
^{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.Simultaneously, calculate the initial Suo Li and the Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) of all ropes of Cable Structure 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 just can determine according to these coordinate datas the geometric properties of Cable Structure.For 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 mechanical calculating benchmark model A of Cable Structure
_{o}, set up initial Cable Structure support 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
^{1}The 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 support coordinate data, Cable Structure modal data, 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 mechanical calculating benchmark model A of Cable Structure
_{o}If 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 mechanical calculating benchmark model A of Cable Structure equally
_{o}No matter which kind of method to obtain A with
_{o}, based on A
_{o}The 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
_{o}Cable Structure support coordinate data form initial Cable Structure support coordinate vector U
_{o}A
_{o}And U
_{o}Be 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 represents 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 represents 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
^{1}Just equal A
_{o}For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A
^{Ti} _{o}", A in circulation each time
^{Ti} _{o}Can bring in constant renewal in as required, during circulation beginning each time, A
^{Ti} _{o}Equal A
^{i}Equally for sake of convenience, name " Cable Structure actual measurement support coordinate vector U
^{Ti}", in circulation each time, 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
^{Ti}, vectorial U
^{Ti}Element and vectorial U
_{o}The coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, for the i time circulation, the last time is upgraded A
^{Ti} _{o}The time Cable Structure support coordinate current data be designated as current cable structural bearings coordinate vector U
^{Ti} _{o}During circulation beginning for the first time, A
^{T1} _{o}Equal A
^{1}, U
^{T1} _{o}Equal U
_{o}
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 subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate measurement subsystem, signal conditioner etc.), 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 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 finished functions such as monitoring that the present invention's progressive health monitoring method of the cable system of hybrid monitoring " when support settlement is arranged based on " required by task wants, record, control, storage, calculating, notice, warning (all work that can finish with computing machine in this specific 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 finish necessary evaluation work automatically.
The 7th step: the step begins circulation running thus, is designated as the i time circulation for sake of convenience, i=1 wherein, and 2,3,4,5 ....Actual measurement (for example measuring with total powerstation) 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
^{Ti}, according to current cable structure actual measurement support coordinate vector U
^{Ti}, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary
^{Ti} _{o}With current cable structural bearings coordinate vector U
^{Ti} _{o}Concrete grammar is:
Actual measurement obtains current cable structure actual measurement support coordinate vector U
^{Ti}After, compare U
^{Ti}And U
^{Ti} _{o}If, U
^{Ti}Equal U
^{Ti} _{o}, then do not need A
^{Ti} _{o}Upgrade;
Actual measurement obtains current cable structure actual measurement support coordinate vector U
^{Ti}After, compare U
^{Ti}And U
^{Ti} _{o}If, U
^{Ti}Be not equal to U
^{Ti} _{o}, then need A
^{Ti} _{o}Upgrade, update method is: calculate U earlier
^{Ti}With U
_{o}Poor, U
^{Ti}With U
_{o}Difference be exactly that the current cable structural bearings is about setting up A
_{o}The time the current support displacement of Cable Structure bearing, represent support displacement with current support displacement vector V, between element among the current 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 current support displacement vector V, and wherein support displacement is exactly the support settlement amount at the component of gravity direction; Upgrade A
^{Ti} _{o}Method be: at A
_{o}The basis on make that the health status of rope is cable system initial damage vector d
^{i} _{o}, more further to A
_{o}In 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 current support displacement vector V, to A
_{o}In the Cable Structure bearing apply the constraint of current support displacement after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
^{Ti} _{o}, upgrade A
^{Ti} _{o}The time, U
^{Ti} _{o}All 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} _{o}The basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector D
^{i} _{u}Concrete grammar is:
A. the i time when beginning circulation or in the 7th step to A
^{Ti} _{o}After upgrading, directly b obtains Cable Structure unit damage monitored quantitative change matrix Δ C to the listed method of steps d set by step
^{i}With nominal unit damage vector D
^{i} _{u}At other constantly, if in step g not to A
^{Ti} _{o}Upgrade, 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} _{o}The 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 to increase unit damage again on the basis of original damage, the rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating, and supposition each time has the unit damage value of the rope of damage can be different from the unit damage value of other ropes, uses " nominal unit damage vector D
^{i} _{u}" write down the unit damage of the supposition of all ropes, calculate each time that all specify the current numerical value of monitored amount in the Cable Structure, the current numerical value of the monitored amount of all that calculate is formed one " the current numerical value vector of the calculating of monitored amount " each time; When hypothesis j root rope has unit damage, available C
^{i} _{Tj}" the current numerical value vector of the calculating of monitored amount " that expression is corresponding; When giving each vectorial element numbering in this step, should use same coding rule with other vector among the present invention, can guarantee any one 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 numerical value of the calculating of monitored amount vector 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 " numerical value change vector δ a C of monitored amount after the unit damage value of supposition
^{i} _{j}"; There is N root rope that N " the numerical value change vector δ C of monitored amount just arranged
^{i} _{j}" (j=1,2,3 ..., N).
D. form " the unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount "
^{i}"; " unit damage monitored quantitative change matrix Δ C
^{i}" each row corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " unit damage monitored quantitative change matrix " and the vectorial d of current name damage
^{i} _{c}With current actual damage vector d
^{i}The element coding rule identical.
Reach in this step when giving each vectorial element numbering thereafter, should use same coding rule with other vector among the present invention, can guarantee any one 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 linear relationship error vector e
^{i}With vectorial g
^{i}Utilize data (" the initial value vector C of monitored amount of front
^{i} _{o}", " unit damage monitored quantitative change matrix Δ C
^{i}"), when the 8th step calculated each time, promptly in calculating each time, have only in the hypothesis cable system increase unit damage again on the basis of rope in original damage in, calculate each time and form a vectorial d of damage
^{i} _{t}, damage vectorial d
^{i} _{t}Element number equal the quantity of rope, vectorial d
^{i} _{t}All elements in have only the numerical value of an element to get to calculate each time in hypothesis increase the unit damage value of the rope of 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 unit damage, be identical with the element of the same numbering of other vectors with the corresponding relation of this rope; With C
^{i} _{Tj}, C
^{i} _{o}, Δ C
^{i}, d
^{i} _{t}Bring formula (13) into, obtain a linear relationship error vector e
^{i}, calculate a linear relationship error vector e each time
^{i}There is N root rope that N calculating is just arranged, N linear relationship error vector e just arranged
^{i}, with this N linear relationship error vector e
^{i}Obtain a vector after the addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e
^{i}Vector g
^{i}Equal final error vector e
^{i}With vectorial g
^{i}Be kept on the hard disc of computer of operation health monitoring systems software, use for health monitoring systems software.
The tenth step: will " the initial value vector C of monitored amount
^{i} _{o}" and " unit damage monitored quantitative change matrix Δ C
^{i}" etc. parameter be kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.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}".
The 11 step: according to " the current numerical value vector C of monitored amount
^{i}" " the initial value of monitored amount vector C together
^{i} _{o}", " unit damage monitored quantitative change matrix Δ C
^{i}" and " the vectorial d of current name damage
^{i} _{c}" between the linear approximate relationship (formula (9)) that exists, calculate the vectorial d of the current name damage of cable system according to multi-objective optimization algorithm
^{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 the vectorial d of current name damage with the goal programming method only
^{i} _{c}Process, 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 (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (21) and the formula (22), γ in the formula (21)
^{i}Be a real number, R is a real number field, and area of space Ω has limited vectorial d
^{i} _{c}Span (the present embodiment requirements vector d of each element
^{i} _{c}Each element be not less than 0, be not more than 1).The meaning of formula (21) is to seek the real number γ of an absolute value minimum
^{i}, make formula (22) be met.G (d in the formula (22)
^{i} _{c}) by formula (23) definition, weighing vector W in the formula (22)
^{i}With γ
^{i}Product representation formula (22) in G (d
^{i} _{c}) and vectorial g
^{i}Between the deviation that allows, g
^{i}Definition referring to formula (15), its value will the 8th the step calculate.Vector W during actual computation
^{i}Can with vectorial 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 the vectorial d of current name damage
^{i} _{c}
minimize γ
^{i} (21)
γ
^{i}∈R，
${d}_{c}^{i}\∈\mathrm{\Ω}$
$G\left({d}_{c}^{i}\right)-{W}^{i}{\mathrm{\γ}}^{i}\≤{g}^{i}---\left(22\right)$
$G\left({d}_{c}^{i}\right)=\mathrm{abs}(\mathrm{\Δ}{C}^{i}\·{d}_{c}^{i}-{C}^{i}+{C}_{o}^{i})---\left(23\right)$
Try to achieve the vectorial d of current name damage
^{i} _{c}After, can be according to the vectorial d of the current actual damage that formula (17) obtain
^{i}Each element, current actual damage vector d
^{i}Have reasonable error exactly but can be more exactly from all ropes, determine the position of damaged cable and separating of degree of injury thereof.If the current actual damage vector d that solves
^{i}The numerical value of a certain element be 0, represent that the pairing rope of this element is intact, 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 12 step:, promptly try to achieve the vectorial d of current name damage in the i time circulation in this circulation
^{i} _{c}After, at first, set up mark vector F according to formula (18), formula (19)
^{i}If mark vector F
^{i}Element be 0 entirely, then got back to for the 7th step and continue this circulation; If mark vector F
^{i}Element be not 0 entirely, then enter next step, i.e. the 13 step.
The 13 step: according to the initial damage vector d that formula (20) calculates next time, i.e. the i+1 time circulation is required
^{I+1} _{o}Each element d
^{I+1} _{Oj}
The 14 step: at the current Mechanics Calculation benchmark model of Cable Structure A
^{Ti} _{o}The basis on, the health status that makes rope is the vectorial d that previous step calculates
^{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 15 step: by to Mechanics Calculation benchmark model A
^{I+1}Calculate corresponding to model A
^{I+1}The numerical value of all monitored amounts of structure, these numerical value are formed next time, required vectorial C promptly circulates for the i+1 time
^{I+1} _{o}, i.e. the initial value vector of monitored amount.
The 16 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} _{o}Equal A
^{I+1}
The 17 step: set up next time, required current cable structural bearings coordinate vector U promptly circulates for the i+1 time
^{Ti+1} _{o}, promptly get U
^{Ti+1} _{o}Equal U
^{Ti} _{o}
The 18 step: got back to for the 7th step, beginning is circulation next time.