Summary of the invention
Technical matters: the invention discloses a kind of based on the identification damaged cable of space coordinate monitoring and the progressive method of support displacement, this method be based on space coordinate monitoring, adopt health monitor method progressive method, that can discern support displacement and damaged cable rationally and effectively.
Technical scheme: structures such as cable-stayed bridge, suspension bridge, truss-frame structure have a common ground; Be exactly that they have many parts that bear tensile load; Like suspension cable, main push-towing rope, hoist cable, pull bar or the like; The common ground of this class formation is to be support unit with rope, cable or the rod member that only bears tensile load, and the present invention is " Cable Structure " with such structure representation for simplicity.In the military service process of Cable Structure; The supporting system of Cable Structure (refers to all ropeway carrying-ropes, reaches all rod members that only bear tensile load that play supporting role; For simplicity, this patent is called " cable system " with whole support unit unifications of this class formation, but in fact cable system not only refers to support rope; Also comprise the rod member that only bears tensile load) can be impaired; The bearing of Cable Structure also displacement possibly occur simultaneously, and these safety that change to Cable Structure are a kind of threats, and establishing the quantity of rope and the quantity sum of support displacement component is N.For narrate convenient for the purpose of, the unified rope of being assessed of claiming of the present invention is " by evaluation object " with support displacement, gives by the evaluation object serial number, the present invention is with representing this numbering with variable j, j=1,2,3 ..., N, so we can say that N is individual by evaluation object.
The present invention is made up of the two large divisions.Be respectively: one, set up by the method for required knowledge base of evaluation object health monitoring systems and parameter, based on knowledge base (containing parameter) and the volume coordinate of surveying Cable Structure by evaluation object health status appraisal procedure; Two, the software and hardware part of health monitoring systems.
First of the present invention: foundation is used for by the method for the knowledge base of evaluation object health monitoring and parameter.Can be successively circularly as follows, laddering carrying out:
The first step: when beginning circulation each time, in the time of at first need setting up or set up this circulation beginning by evaluation object initial health vector d
_{o} ^{i}(i=1,2,3 ...), set up the initial Mechanics Calculation benchmark model A of Cable Structure
_{o}(for example finite element benchmark model, A in the present invention
_{o}Be constant), set up the Mechanics Calculation benchmark model A of Cable Structure
^{i}(finite element benchmark model for example, i=1,2,3 ...).Letter i is except the place of representing number of steps significantly, and alphabetical in the present invention i only representes cycle index, i.e. the i time circulation.
The Cable Structure that needs during the i time circulation beginning " initial health vector d
_{o} ^{i}" (shown in (1)), use d
_{o} ^{i}Cable Structure is (with Mechanics Calculation benchmark model A when representing the i time circulation beginning
^{i}The initial health of Cable Structure expression).
${d}_{o}^{i}={\left[\begin{array}{cccccccccc}{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 j of cable system by the current health status of evaluation object, if should be the rope (or pull bar) in the cable system, d so by evaluation object
_{i}Represent its current damage, d
_{i}Being to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be displacement component, a d so of a bearing by evaluation object
_{i}Represent its present bit shift value.T representes the transposition (back together) of vector in the formula (1).
Set up the initial health vector during circulation beginning for the first time and (be designated as d according to formula (1)
^{1} _{o}) time, the data and the support displacement that utilize the Non-Destructive Testing data etc. of rope can express the health status of rope are measured and are set up by evaluation object initial health vector d
^{1} _{o}If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vectorial d
^{1} _{o}In get 0 with each element numerical value of Suo Xiangguan.
The i time (i=2,3,4,5,6 ...) when beginning circulation need by evaluation object initial health vector d
^{1} _{o}, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
The Mechanics Calculation benchmark model that need set up during the i time circulation beginning or the Mechanics Calculation benchmark model of having set up are designated as A
^{i}
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
_{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 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); According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method (for example finite element method) to set up A
_{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.A
_{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}A
^{1}Corresponding by the health status of evaluation object by d
^{1} _{o}Describe.
The i time (i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins
^{i}, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Existing Mechanics Calculation benchmark model A
^{1}With by evaluation object initial health vector d
^{1} _{o}After, model A
^{1}In each by the health status of evaluation object by vectorial d
^{1} _{o}Express.At A
^{1}The basis on, all are changed to 0 by the health status numerical value of evaluation object, mechanical model A
^{1}Be updated to one all be that 0 mechanical model (is designated as A by the health status of evaluation object
^{0}), mechanical model A
^{0}Be actually the corresponding mechanical model of Cable Structure of excellent no support displacement.Might as well claim model A
^{0}For the not damaged of Cable Structure does not have support displacement model A
^{0}
" the whole monitored spatial data of structure " described by the volume coordinate specified point of K on the structure, that reach L assigned direction of each specified point, and the variation of structure space coordinate data is exactly the variation of all volume coordinate components of K specified point.(individual volume coordinate measured value of M=K * L) or calculated value characterize the structure space coordinate information to each total M.K and M must not be less than by the quantity N of evaluation object.
For simplicity, in the present invention " the monitored spatial data of structure " abbreviated as " monitored amount ".When mentioning " so-and-so matrix of monitored amount or so-and-so vector " in the back, also can be read as " volume coordinate so-and-so matrix or so-and-so vector ".
Among the present invention with monitored amount initial value vector C
^{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 the i time monitored amount of circulation ".
${C}_{o}^{i}={\left[\begin{array}{cccccccccc}{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 that M monitored amount by front definition forms according to certain series arrangement, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
During circulation beginning for the first time, " the initial value vector 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 obtain, concrete grammar is in the back literary composition narration.
Second step: circulation each time needs to set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector ", and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C
^{i}, " nominal unit damage vector " that the i time circulation set up is designated as D
^{i} _{u}, i=1,2,3 ...
The Cable Structure " unit damage monitored numerical quantity transformation matrices " that circulation is for the first time set up is designated as Δ C
^{1}Set up Δ C
^{1}Process following:
Mechanics Calculation benchmark model A in Cable Structure
^{1}The basis on carry out several times and calculate, equal N on the calculation times numerical value.Calculating hypothesis each time has only one by evaluation object unit damage to be arranged; Concrete, if should be a supporting rope in the cable system by evaluation object, so just this supporting rope of hypothesis has unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage); If should be the displacement component of a direction of a bearing by evaluation object; Just suppose that this bearing is this sense of displacement generation unit displacement (for example 10mm, 20mm, 30mm etc. are unit displacement).For narrating conveniently, the damage and the support displacement of the supporting rope that the present invention will suppose are referred to as unit damage.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, calculate, by the health status numerical value of evaluation object---be called nominal health status vector d
^{i} _{c}(i=1; 2,3 ...); All with respect to this time when beginning circulation, with the health status of Cable Structure as being healthy fully speech, therefore must foundation after the formula that provides of the literary composition nominal health status numerical value that will calculate be converted into true health status numerical value).With a round-robin occur in calculating each time unit damage be different from by evaluation object other time occur in calculating unit damage by evaluation object; And supposition has the unit damage value by evaluation object of unit damage can be different from other by the unit damage value of evaluation object each time, with " nominal unit damage vector D
^{i} _{u}" (shown in (3)) write down that all are by the unit damage of the supposition of evaluation object 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 j of hypothesis had unit damage by evaluation object, 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 be exactly under this condition (with have unit damage by the mark that is numbered of evaluation object) " the numerical value change vector of monitored amount " (when j is had unit damage by evaluation object, use δ C
^{1} _{j}The numerical value change vector of representing monitored amount, δ C
^{1} _{j}Definition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by 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 because that supposition has a unit damage when calculating had unit damage (D for example by evaluation object (for example j by evaluation object)
_{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 N " the numerical value change vector of monitored amount " just to be arranged by evaluation object; The numerical value change vector of each monitored amount has M (general; The individual element of M>=N) is formed " the unit damage monitored numerical quantity transformation matrices Δ 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 suc as formula shown in (8).
${C}_{u}^{i}={\left[\begin{array}{cccccccccc}{C}_{u1}^{i}& {C}_{u2}^{i}& \·& \·& \·& {C}_{\mathrm{uj}}^{i}& \·& \·& \·& {C}_{\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) j of supposition is individual by the unit damage numerical value of evaluation object 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}{cccccccccc}{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; The i time circulation of the expression of M>=N) be because j is individual when by evaluation object unit damage being arranged, 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}_{0}^{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) j of expression had unit damage by evaluation object, D in the formula
^{i} _{Uj}Be vectorial D
^{i} _{u}In j element.Vector δ C
^{i} _{j}Definition suc as formula shown in (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 the change amount of calculating a gained k monitored amount when j has unit damage by evaluation object unit damage D with respect to supposition
^{i} _{Uj}Rate of change, it defines suc as formula shown in (7).
${\mathrm{\δC}}_{j}^{i}={\left[\begin{array}{cccccccccc}{\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 was narrated in front in the formula (7).
${\mathrm{\ΔC}}^{i}=\left[\begin{array}{cccccccccc}{\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, by evaluation object unit damage D is arranged because j is individual
^{i} _{Uj}Cause, the relative value of all monitored amounts changes.Matrix Δ C
^{i}Coding rule and the front vector d of row (subscript j)
^{i} _{o}The coding rule of subscript j of element identical.
The 3rd step: identification is by the current health status of evaluation object.Detailed process is following.
I (i=1,2,3 ...) in the inferior circulation, " 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 numerical quantity transformation matrices Δ 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); By the vectorial d of the current name damage of evaluation object
^{i} _{c}Definition see formula (12).
${C}^{i}={\left[\begin{array}{cccccccccc}{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}{cccccccccc}{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 ...., be that Cable Structure j is individual by the current nominal impairment value of evaluation object in the i time circulation N), vectorial d
^{i} _{c}The coding rule and the matrix Δ C of subscript j of element
^{i}The coding rule of row identical.
When by evaluation object actual damage or support displacement when not too big; 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 (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}" come directly to find the solution to obtain the vectorial d of current name damage
^{i} _{c}Damage vectorial d and obtain current name
^{i} _{c}Acceptable separating (promptly have reasonable error, but can be more accurately from cable system, confirm damaged cable position and degree of injury thereof, confirm the support displacement amount) 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}{cccccccccc}{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), unit damage monitored numerical quantity transformation matrices Δ 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 current name damage
^{i} _{c}Acceptable separating, current actual damage vector d
^{i}The element of (formula (16) is seen in definition) can calculate according to formula (17), has just obtained by the current actual damage vector of evaluation object d
^{i}Thereby, can be by d
^{i}Confirm damaged cable position and degree of injury, confirm the support displacement amount, realized that just damage identification and support displacement discern.
${d}^{i}={\left[\begin{array}{cccccccccc}{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}(i=1,2,3, J=1,2,3 ...., N) j is by the actual damage value of evaluation object in the i time circulation of expression, and formula (17) is seen in its definition, if this was the rope (or pull bar) in the cable system, d so by evaluation object
^{i} _{j}Represent its current damage, d
^{i} _{j}Being to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%, if should be displacement component, a d so of a bearing by evaluation object
^{i} _{j}Represent its present bit shift value.Vector d
^{i}Coding rule and the formula (1) of element 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} _{Cj}Be vectorial d
^{i} _{c}J element.
The 4th step: judge whether to finish this (the i time) circulation, if, then accomplish 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 following.
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 Structure; If mark vector F
^{i}Element be not 0 entirely, then accomplish subsequent step after, get into circulation next time.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, order is d by the health status situation of evaluation object
^{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, through 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}{cccccccccc}{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) representes that j can only be got 0 and 1 two amount by the damage characteristic of evaluation object, 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 current name damage
^{i} _{c}J element (seeing formula (12)), they represent that all j is by the relevant information of evaluation object.
${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 current name damage
^{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 monitored amount monitoring system, signal picker and computing machine etc.Require to monitor in real time or quasi real time each monitored amount.
Software should be used following function by tool: software section should be accomplished the process that first of the present invention sets, and promptly accomplishes functions such as needed among the present invention, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.
The inventive method specifically comprises:
A. for for the purpose of narration is convenient, the present invention is unified to claim that the supporting rope of being assessed is that establishing the quantity of the supporting rope of being assessed and the quantity sum of support displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with the support displacement component; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention representes this numbering with variable j, j=1, and 2,3 ..., N;
B. confirm the measured point with monitored volume coordinate of appointment, give all specified point numberings; Confirmed each measurement point with monitored volume coordinate component, give all measured volume coordinate components numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored spatial data of structure " is made up of above-mentioned all measured volume coordinate components; For simplicity, in the present invention " the monitored spatial data of structure " is called " monitored amount "; The quantity sum of all measured volume coordinate components must not be less than N;
C. utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector d by the Non-Destructive Testing data of evaluation object etc.
^{i} _{o}If during not by the Non-Destructive Testing data of evaluation object, vectorial d
^{i} _{o}Each element numerical value get 0; Vector d
^{i} _{o}The coding rule of element with identical by the coding rule of evaluation object; The present invention representes cycle index with i, i=1, and 2,3 ...; Here be circulation for the first time, i gets 1, the initial health vector d that promptly sets up here
^{i} _{o}Can be embodied as d
^{1} _{o}
D. setting up initial health vector d
^{1} _{o}The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount
^{i} _{o}Here be circulation for the first time, i gets 1, the initial value vector C of the monitored amount of promptly setting up here
^{i} _{o}Can be embodied as C
^{1} _{o}Obtain monitored amount initial value vector C in actual measurement
^{1} _{o}The 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;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure support coordinate data are set up Cable Structure
^{i}Here be circulation for the first time, i gets 1, the Mechanics Calculation benchmark model A of the Cable Structure of promptly setting up here
^{i}Can be embodied as A
^{1}
F. at Mechanics Calculation benchmark model A
^{i}The basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " nominal unit damage vector D
^{i} _{u}";
G. 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 element each vector, that numbering is identical of occurring before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
H. define the vectorial d of current name damage
^{i} _{c}With current actual damage vector d
^{i}, the element number of two damage vectors equals by the quantity of evaluation object, the vectorial d of current name damage
^{i} _{c}Element numerical value represent corresponding to the current nominal degree of injury or the support displacement of evaluation object, current actual damage vector d
^{i}Element numerical value represent corresponding by the current actual damage degree or the support displacement of evaluation object; The element number of the element of two damage vectors equals by the quantity of evaluation object; The element of two damage vectors and be one-to-one relationship between the evaluation object, the coding rules of the element of two damage vectors are with identical by the coding rule of evaluation object;
I. 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 numerical quantity transformation matrices Δ 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
J. 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;
Current actual damage vector d
^{i}Element numerical value represent corresponding to the actual damage degree of evaluation object or actual support displacement, according to current actual damage vector d
^{i}Just can define the impaired and degree of injury of which rope, just can confirm actual support displacement; If a certain element of current actual damage vector is corresponding to being a rope in the cable system; And its numerical value is 0, representes that the pairing rope of this element is intact, does not damage; If its numerical value is 100%; Represent that then the pairing rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this rope; If a certain element of current actual damage vector is corresponding to displacement component, a d so of a bearing
^{i} _{j}Represent its present bit shift value;
K. 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}{cccccccccc}{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 current name damage
^{i} _{c}J element, they all represent j by the relevant information of evaluation object, j=1 in the formula 4,2,3 ..., N;
If mark vector F l.
^{i}Element be 0 entirely, then get back to g step and continue this circulation; If mark vector F
^{i}Element be not 0 entirely, then get into next step, i.e. m step;
M. the initial damage vector d that calculates next time, promptly circulates required the i+1 time according to formula 5
^{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 current name damage
^{i} _{c}J element, F
^{i} _{j}Be mark vector F
^{i}J element, j=1 in the formula 5,2,3 ..., N; Vector d
^{I+1} _{o}The coding rule of element with identical by the coding rule of evaluation object;
N. at Mechanics Calculation benchmark model A
^{i}The basis on, order is d by the health status of evaluation object
^{I+1} _{o}The back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time
^{I+1}
O. pass through Mechanics Calculation benchmark model A
^{I+1}Calculate corresponding to model A
^{I+1}Structure all monitored strains point, with the monitored strain numerical value that should change direction, 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}
P. get back to the f step, beginning is circulation next time.
In step f, at Mechanics Calculation benchmark model A
^{i}The basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C
^{i}" and " nominal unit damage vector D
^{i} _{u}" concrete grammar be:
F1. at the Mechanics Calculation benchmark model A of Cable Structure
^{i}The basis on carry out the several times Mechanics Calculation, equal N on the calculation times numerical value; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or displacement, to be increased unit damage or unit displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases unit displacement again at this sense of displacement; Increase again in calculating each time unit damage or unit displacement be different from by evaluation object other time increase again in calculating unit damage or unit displacement by evaluation object, with " nominal unit damage vector D
^{i} _{u}" unit damage or the unit displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector; Obtain a monitored quantitative changeization vector, have N evaluation object that N monitored quantitative changeization vector just arranged;
F3. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the monitored amount unit change of the Cable Structure matrix Δ C that the N row are arranged successively
^{i}
Beneficial effect: method disclosed by the invention can identify the health status (position and the degree of injury that comprise all damaged cables) of Cable Structure support displacement and cable system simultaneously, and system and method disclosed by the invention is very useful to the safety of Cable Structure.
Embodiment
To the health monitoring of Cable Structure, the invention discloses a kind of system and method that can monitor in the Cable Structure health status of each root rope and each support displacement component in the cable system rationally and effectively simultaneously.The following explanation of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
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 practical implementation, the following step is a kind of in the various steps that can take.
The first step: for the purpose of narrating conveniently, the present invention is unified to claim that supporting rope and the support displacement component assessed are that establishing the quantity of the supporting rope of being assessed and the quantity sum of support displacement component is N, is N by the quantity of evaluation object promptly by evaluation object; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention representes this numbering with variable j, j=1, and 2,3 ..., N.
Confirm that (promptly all characterize the specified point of structure space coordinates for the measured point of appointment; Be provided with K specified point); Each specified point can be exactly a near point the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) of each root rope; This specified point can also be a near point the structural bearings, perhaps directly is exactly the structural bearings fulcrum, gives all specified point numberings; Confirm the measured volume coordinate component (established each measurement point L measured volume coordinate component arranged) of each measurement point, gave the measured volume coordinate component numbering of all appointments.Above-mentioned numbering will be used to generate the vector sum matrix equally in subsequent step." the whole monitored spatial data of structure " described by L volume coordinate component K specified point, that cross each specified point on top definite structure, and the variation of structure space coordinate is exactly variations all specified points, all space specified coordinate components.(individual volume coordinate component measurement value of M=K * L) or calculated value characterize the volume coordinate information of structure to each total M.K and M must not be less than N.For simplicity, in the present invention " the monitored spatial data of structure " is called " monitored amount ".
Second step: utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector d by the Non-Destructive Testing data of evaluation object etc.
^{1} _{o}If during not by the Non-Destructive Testing data of evaluation object, vectorial d
^{1} _{o}Each element numerical value get 0; Vector d
^{1} _{o}The coding rule of element with identical by the coding rule of evaluation object.
The 3rd step: at initial health vector d
^{1} _{o}The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount
^{1} _{o}
The 4th step: the initial value vector C that obtains monitored amount in actual measurement
^{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 confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
According to 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), the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure support coordinate data are set up Cable Structure
_{o}, based on Mechanics Calculation benchmark model A
_{o}The computational data that calculates structure must be very near its measured data, and error generally must not be greater than 5%.
A
_{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 representes cycle index; Alphabetical i is except the place of representing number of steps significantly in the application form of the present invention, and alphabetical i only representes cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A
^{1}, A among the present invention
^{1}Just equal A
_{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 spatial coordinate measuring system, signal conditioner etc.), signal (data) collector, the computing machine and the panalarm of communicating by letter.Each monitored amount 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 when being changed by the health status of evaluation object, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 6th step: establishment and the health monitoring systems software of pass line structure on supervisory control comuter.All move this software at circulation time each time, this software is all the time in operation in other words.This software will be accomplished functions such as the needed monitoring of each task of the present invention, record, control, storage, calculating, notice, warning (being all work that can accomplish with computing machine in this practical implementation method); And can regularly or by the personnel operation health monitoring systems generate Cable Structure health condition form; Can also be according to the condition of setting (for example damage reach a certain value), notice or prompting monitor staff notify specific technician to accomplish necessary evaluation work automatically.
The 7th step: the step begins circulation running thus, for narration conveniently is designated as the i time circulation, and i=1 wherein, 2,3,4,5 ....
The 8th step: the Mechanics Calculation benchmark model in Cable Structure is designated as A
^{i}The basis on carry out the several times Mechanics Calculation, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C
^{i}With nominal unit damage vector D
^{i} _{u}Concrete grammar is:
A. when the i time circulation beginning, at the Mechanics Calculation benchmark model A of Cable Structure
^{i}The basis on carry out the several times Mechanics Calculation, equal N on the calculation times numerical value; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one to be increased on the basis of original damage or displacement by evaluation object again unit damage or unit displacement are arranged; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases unit displacement again at this sense of displacement; Increase again in calculating each time unit damage or unit displacement be different from by evaluation object other time increase again in calculating unit damage or unit displacement by evaluation object, with " nominal unit damage vector D
^{i} _{u}" unit damage or the unit displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed same monitored amount or same by the relevant information of evaluation object object.
B. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector, obtains a monitored quantitative changeization vector δ C
^{i} _{j}There is N N monitored quantitative changeization vector δ C just to be arranged by evaluation object
^{i} _{j}(j=1,2,3 ..., N).
C. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the monitored amount unit change of the Cable Structure matrix Δ C that the N row are arranged successively
^{i}" unit damage monitored quantitative change matrix Δ C
^{i}" the vectorial d of current name damage of coding rule and back definition of row
^{i} _{c}With current actual damage vector d
^{i}The element coding rule identical.
In this step, reach when giving each vectorial element numbering thereafter; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.
The 9th step: set up 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 one by evaluation object on the basis of original damage or displacement, increase again unit damage or unit displacement are arranged in, calculate each time and form a health status vector d
^{i} _{t}, health status vector d
^{i} _{t}Element number equal by the quantity of evaluation object vectorial d
^{i} _{t}All elements in have only the numerical value of an element to get to calculate each time in hypothesis increase unit damage value or the unit displacement value of increase of the rope of unit damage, d
^{i} _{t}The numerical value of other element get 0, what that was not that numbering and the supposition of 0 element increase unit damage or unit displacement is identical by the corresponding relation of evaluation object, with the element of the same numberings 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, formula (13) d
^{i} _{c}Use d
^{i} _{t}Bring into, obtain a linear relationship error vector e
^{i}, calculate a linear relationship error vector e each time
^{i}There is N N calculating just to be arranged, N linear relationship error vector e just arranged by evaluation object
^{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, supply health monitoring systems software to use.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 with the mode of data file.
The tenth 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
^{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 concrete implementation procedure of other algorithm can realize according to the requirement of its specific algorithm in a similar fashion.
According to the goal programming method, formula (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (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 minimum real number γ of an absolute value
^{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{ads}(\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 confirm more exactly damaged cable the separating of position and degree of injury thereof, can confirm separating of all support displacements more exactly.d
^{i}Each element corresponding to one by the health status of evaluation object; If should be the rope (or pull bar) in the cable system by evaluation object; Its current damage of the numeric representation of this element so; The numerical value of this element is to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%; If being somebody's turn to do by evaluation object be a displacement component of a bearing, so its present bit shift value of the numeric representation of this element.
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, 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 tenth step and continue this circulation; If mark vector F
^{i}Element be not 0 entirely, then get into 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 Cable Structure Mechanics Calculation benchmark model A
^{i}The basis on, order is the last vectorial d that calculates of step by the health status of evaluation object
^{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: through 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}, promptly the initial value of monitored amount is vectorial.
The 16 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 17 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 18 step: got back to for the 7th step, beginning is circulation next time.