CN101793622B - Method for distinguishing slack supporting cable based on cable force monitoring during support settlement - Google Patents

Abstract

The invention relates to a method for distinguishing a slack supporting cable based on cable force monitoring during support settlement; the coordinates of a support of a monitoring structure are used for deciding whether a mechanical calculation benchmark module of the structure needs to be updated or not; only when the coordinates of the support of the structure changes, the mechanical calculation benchmark module of the structure is updated, thus obtaining the new mechanical calculation benchmark module of the structure which records the support settlement of the structure; and based on the module, a unit damage monitored quantity transformation matrix is obtained by calculation. According to the approximately linear relationship of the current vector of the strain with the initial vector of the strain, the virtual unit damage strain transformation matrix and the current virtual damage vector, the virtual damaged cable is calculated and distinguished. After a nondestructive testing method is used for distinguishing the true damaged cable, the remained virtual damaged cable is the slack cable, i.e., the cable force needs to be regulated; and the length of the cable needing to be regulated can be determined according to the relationship between the slack degree and the virtual damage degree.

Description

When support settlement is arranged based on the method for the lax support cable of the identification of cable force monitoring

Technical field

When support settlement is arranged, rope supporting structure (particularly large-scale Cable Structure is discerned in the monitoring that the present invention is based on Suo Li equivalent, for example large-scale cable-stayed bridge, suspension bridge) cable system (referring to all support cables) in need adjust the support cable of Suo Li, and provide the long adjustment amount of concrete rope, belong to the engineering structure security fields.

Background technology

Cable system is Cable Structure (particularly large-scale Cable Structure normally; for example large-scale cable-stayed bridge, suspension bridge) key components; owing to reason such as lax; new construction is completed, and the Suo Li of support cable can change usually after a period of time; the lax variation that also can cause the supporting cable force of its support cable behind the structure long service; these change the variation that all will cause structural internal force; safety to structure causes harmful effect; will cause the inefficacy of structure when serious, therefore accurately and timely discern the support cable that needs to adjust Suo Li and be very important.

The health status of support cable system changes and (for example takes place lax; damage etc.) after; except meeting causes the variation of Suo Li; also can cause the variation of other measurable parameter of structure; for example can cause the variation of supporting cable force; in fact the variation of Suo Li has comprised the health status information of cable system; that is to say the health status that to utilize the rope force data to judge structure; can (the present invention be called monitored Suo Li " monitored amount " based on cable force monitoring; the back is mentioned " monitored amount " and just is meant monitored Suo Li) discern damaged cable; monitored amount is except the influence that is subjected to the cable system health status; also can be subjected to the influence of Cable Structure support settlement (usually can take place), also not have a kind of disclosed at present; effectively health monitoring systems and method have solved this problem.Therefore can monitor based on monitored amount and discern the rope that needs to adjust Suo Li, like this when support settlement is arranged, a method of can rational and effective setting up the relation of (specifically the characteristic parameter according to rope characterizes the rope that needs to adjust Suo Li) between the characteristic parameter of monitored amount with all ropes just must be arranged, and the need of setting up based on this method are adjusted the recognition result of the support cable of Suo Li just can be more credible.

Summary of the invention

Technical matters: the objective of the invention is when the Cable Structure bearing has sedimentation, at the identification problem of support cable in the cable system in the Cable Structure, that need to adjust Suo Li, disclose a kind of monitoring based on Suo Li equivalent, can discern rationally and effectively need to adjust Suo Li support settlement is arranged the time based on the method for the lax support cable of the identification of cable force monitoring.

Technical scheme: the reason according to the Suo Li of support cable changes can change the three kinds of situations that be divided into the Suo Li of support cable: the one, and support cable has been subjected to damage, and for example localized cracks and corrosion or the like have appearred in support cable; The 2nd, support cable and not damaged, but variation has also taken place in Suo Li, the one of the main reasons that this variation occurs is that variation has taken place the Suo Changdu (be called drift, the present invention specially refers to the drift of that section rope between support cable two supporting end points) under the support cable free state (this moment, Suo Zhangli claimed that also Suo Li is 0); The 3rd, support cable and not damaged, but the Cable Structure bearing has had displacement (wherein the component at gravity direction just is called as sedimentation), also can cause the variation of structural internal force, also will cause the variation of Suo Li certainly.One of fundamental purpose of the present invention is exactly when support displacement is arranged, and identify drift the support cable that changes has taken place, and identify the change amount of their drift, and this change amount provides direct foundation for the Suo Li adjustment of this rope.The reason that the support cable drift changes is not single, and for convenience, the present invention is referred to as slack line with the support cable that drift changes.

The present invention is made up of the two large divisions.Be respectively: one, set up required knowledge base of the health monitoring systems be used for discerning support cable cable system, that need to adjust Suo Li and parameter method, based on knowledge base (containing parameter), based on actual measurement Cable Structure support coordinate, adjust the method for the support cable of Suo Li based on need monitoring, the identification Cable Structure of monitored amount equivalent.Two, the software and hardware part of health monitoring systems.

First of the present invention: set up required knowledge base of the health monitoring systems be used for discerning support cable cable system, that need to adjust Suo Li and parameter method, based on knowledge base (containing parameter), based on actual measurement Cable Structure support coordinate, adjust the method for the support cable of Suo Li based on need monitoring, the identification Cable Structure of monitored amount equivalent.Can carry out as follows, to obtain the health status assessment of cable system more accurately.

The first step: at first set up the initial virtual lesion vector of cable system d _o(, claim here " virtual lesion ", back with) because in fact support cable may be lax and not damage is the expression difference, set up the initial mechanical calculating benchmark model A of Cable Structure _o(for example finite element benchmark model, A in the present invention _oBe constant).

If total N root rope in the cable system, " initial virtual lesion vector is designated as d to cable system _o" (as the formula (1)), use d _oThe expression Cable Structure is (with the initial mechanical calculating benchmark model A of Cable Structure _oThe health status of cable system expression).

d _o=[d _O1d _O2D _OjD _ON] ^T(1) d in the formula (1) _Oj(j=1,2,3 ...., N) expression A _oIn the initial virtual lesion value of j root rope of cable system, d _OjBeing that to represent that j root rope not damaged did not have at 0 o'clock lax, 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%.T represents the transposition (back together) of vector in the formula (1).

Set up the initial virtual lesion vector of cable system and (be designated as d according to formula (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 the initial virtual lesion vector of cable system d _oIf when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vectorial d _oEach element numerical value get 0.

Set up the initial mechanical calculating benchmark model A of Cable Structure _o(for example finite element benchmark model) and current Mechanics Calculation benchmark model A ^t _oThe method of (for example finite element benchmark model).A in the present invention _oBe constant.A ^t _oBring in constant renewal in.Set up A _oAnd A ^t _oMethod as follows:

Design drawing, as-constructed drawing and the measured data of the Cable Structure in being completed according to Cable Structure (comprises that the Non-Destructive Testing data etc. of rope can express measured datas such as the data of the health status of rope, Cable Structure shape data, structure angle-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), utilize mechanics method (for example finite element method) to set up A _oIf there is not the measured data of the structure in the Cable Structure completion, so just before setting up health monitoring systems, structure is surveyed, the measured data that obtains Cable Structure (comprises the Cable Structure shape data, the rope force data, the draw-bar pull data, Cable Structure support coordinate data, measured datas such as Cable Structure modal data, to cable-stayed bridge, suspension bridge and the bridge type data of Yan Shiqiao, the rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope), design drawing according to these data and Cable Structure, as-constructed drawing utilizes mechanics method (for example finite element method) to set up A _oNo matter which kind of method to obtain A with _o, based on A _oThe Cable Structure computational data that calculates (to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge) must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this _oStrain computational data, Suo Li computational data, Cable Structure shape computational data and displacement computational data, Cable Structure angle-data etc. under the analog case of calculating gained, the measured data when truly taking place near institute's analog case reliably.Corresponding to A _oCable Structure support coordinate data form initial Cable Structure support coordinate vector U _oA _oAnd U _oBe constant.

If total N root support cable in the cable system, structure rope force data is just described by the Suo Li of N root support cable.For simplicity, in the present invention " the monitored rope force data of structure " abbreviated as " monitored amount ".When mentioning " monitored amount so-and-so matrix or so-and-so vector " in the back, also can be read as " Suo Li so-and-so matrix or so-and-so vector ".

Among the present invention with monitored amount initial vector C _oThe vector (seeing formula (1)) that the initial value of all monitored amounts of expression Cable Structure is formed.Requirement is obtaining A _oThe time obtain C _oBecause of subject to the foregoing, the monitored amount of calculating gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial monitored amount reliably, in the narration of back, will represent this calculated value and measured value with prosign.

C _o=[C _O1C _O2C _OjC _OM] ^T(1) C in the formula (1) _Oj(j=1,2,3 ...., M; M 〉=N) is the original bulk of j monitored amount in the Cable Structure, this component according to coding rule corresponding to specific j monitored amount.T represents the transposition (back together) of vector.

The vector of forming by the currency of all monitored amounts in the Cable Structure with the current numerical value vector of monitored amount C among the present invention (formula (2) is seen in definition).

C=[C ₁C ₂C _jC _M] ^T(2) C in the formula (2) _j(j=1,2,3 ...., M; M 〉=N) is the currency of j monitored amount in the Cable Structure, this component C _jAccording to coding rule and C _OjCorresponding to same " monitored amount ".

Second step: the current Mechanics Calculation benchmark model A that sets up Cable Structure ^t _o(finite element benchmark model for example, A in the health monitoring systems operational process ^t _oBring in constant renewal in) and current cable structure actual measurement support coordinate vector U ^tIn Cable Structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data (all data composition current cable structure actual measurement support coordinate vector U ^t, vectorial U ^tDefinition mode and vectorial U _oIdentical).Cable Structure support coordinate current data when for simplicity, being upgraded current Mechanics Calculation benchmark model the last time is designated as current cable structural bearings coordinate vector U ^t _oSet up and renewal A ^t _oMethod be: in the moment that health monitoring systems is started working for the first time, the current Mechanics Calculation benchmark model A of Cable Structure ^t _oJust equal A _oIn Cable Structure military service process, constantly actual measurement obtains Cable Structure support coordinate data and obtains current cable structure actual measurement support coordinate vector U ^tIf, U ^tEqual U ^t _o, then do not need A ^t _oUpgrade; If U ^tBe not equal to U ^t _o, then need A ^t _oUpgrade, at this moment U ^tWith U _oDifference be exactly the Cable Structure bearing about initial position (corresponding to A _o) support displacement (representing support displacement with support displacement vector V, is exactly support settlement in the displacement of gravity direction).Upgrade A ^t _oMethod be: to A _oIn the Cable Structure bearing apply the 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 _oIn 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 ^t _o, upgrade A ^t _oAfter, U ^t _oAll elements numerical value U ^tAll elements numerical value replaces, and has promptly upgraded U ^t _o, so just obtained correctly corresponding to A ^t _oU ^t _o

The 3rd step: set up " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and " nominal virtual unit damage vector D _u", Δ C and D _uBring in constant renewal in, promptly upgrading current Mechanics Calculation benchmark model A ^t _oThe time, upgrade virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _u

Set up and renewal virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _uProcess as follows:

Current Mechanics Calculation benchmark model A in Cable Structure ^t _oThe basis on carry out several times and calculate, equal the quantity of all ropes on the calculation times numerical value.Calculating each time in the hypothesis cable system has only a rope to increase virtual unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is virtual unit damage) again on the basis of original virtual lesion (original virtual lesion can be 0, can not be 0 also).Being convenient and calculating, can all be the bar structure health status when setting virtual unit damage as being healthy fully, and set on this basis virtual unit damage (in subsequent step, damage numerical value that calculate, rope---be called nominal virtual lesion d _c, all be with respect to the health status of rope as being healthy fully speech, the nominal virtual lesion that the formula that therefore must foundation hereinafter provides will calculate be converted into actual and virtual and damage).The rope that occurs virtual lesion in calculating each time is different from the rope of the virtual lesion that occurs in other time calculating, and supposition each time has the virtual unit damage value of the rope of virtual lesion can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector D _u" (as the formula (3)) write down the virtual unit damage of the supposition of all ropes, is designated as D _uCalculate 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 _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 _o, the gained vector is exactly that " the numerical value change vector of monitored amount " of (is mark with the position of rope that virtual unit damage is arranged or numbering etc.) (when j root rope has virtual unit damage, used δ C under this condition _jThe numerical value change vector of representing monitored amount, δ C _jDefinition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vectorial D for formula (4) _uJ element D _UjGained), the numerical value change of monitored amount vector δ C _jEach element representation since when calculating supposition the virtual unit damage (D for example of the Na Gensuo (for example j root rope) of virtual 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 virtual unit damage D of supposition _UjRate 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), form " unit damage monitored quantitative change matrix Δ C " (capable N row of M) that M * N element arranged, each vectorial δ C successively by this N " the numerical value change vector of monitored amount " _j(j=1,2,3 ...., be the row of matrix Δ C N), the definition of Δ C is as the formula (8).

D _u=[D _U1D _U2D _UjD _UN] ^T(3) nominal virtual unit damage vector D in the formula (3) _uElement D _Uj(j=1,2,3 ...., N) the virtual unit damage numerical value of the j root rope of expression supposition, vectorial D _uIn the numerical value of each element can be the same or different.

C _Tj=[C _Tk1C _Tk2C _TjkC _TjM] ^T(4) Elements C in the formula (4) _Tjk(j=1,2,3 ...., N; K=1,2,3 ...., M; When the expression of M 〉=N) j root rope has virtual unit damage, according to the current numerical value of calculating of the monitored amount of pairing k the appointment of coding rule.

$\mathrm{\δ}{C}_j=\frac{C_{\mathrm{tj}}-C_o}{D_{\mathrm{uj}}}---\left(5\right)$

Subscript j in the formula (5) (j=1,2,3 ...., N) expression j root rope has virtual unit damage, D in the formula _UjBe vectorial D _uIn j element.Vector δ C _jDefinition as the formula (6), δ C _jK (k=1,2,3 ...., M; The individual element δ C of M 〉=N) _JkWhen matrix Δ C is set up in expression, suppose that the change amount of calculating a gained k monitored amount when j root rope has virtual unit damage is with respect to the virtual unit damage D that supposes _UjRate of change, it defines as the formula (7).

δC _j＝[δC _j1?δC _j2…δC _jk…δC _jM] ^T (6)

$\mathrm{\δ}{C}_{\mathrm{jk}}=\frac{C_{\mathrm{tjk}}-C_{\mathrm{ok}}}{D_{\mathrm{uj}}}---\left(7\right)$

The definition of each amount has been previously described in the formula (7).

Δ C=[δ C ₁δ C ₂δ C _jδ C _N] vectorial δ C in (8) formula (8) _j(j=1,2,3 ...., N) expression is because j root rope has virtual unit damage D _UjCause, the relative value of all monitored amounts changes.The coding rule of the row (subscript j) of matrix Δ C and front vector d _oThe coding rule of subscript j of element identical.

In Cable Structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, in case monitor U ^tBe not equal to U ^t _o, then needed to get back to second step to A ^t _oUpgrade, to A ^t _oEntering this step after upgrading again upgrades Δ C.In fact Δ C brings in constant renewal in, and is promptly upgrading current Mechanics Calculation benchmark model A ^t _oAfterwards, upgrade the monitored numerical quantity transformation matrices of virtual unit damage Δ C.

The 4th step: the current health status of identification cable system.Detailed process is as follows.

Cable system " current (calculating or actual measurement) numerical value vector C of monitored amount " is " the initial value vector C of monitored amount together _o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and " current nominal virtual lesion vector d _c" between linear approximate relationship, shown in (9) or formula (10).

C＝C _o+ΔC·d _c (9)

C-C _o=Δ Cd _c(10) definition of current (calculate or actual measurement) numerical value vector C of monitored amount is similar to the initial value vector C of monitored amount in formula (9) and the formula (10) _oDefinition, see formula (11); Cable system " current nominal virtual lesion vector d _c" definition see formula (12).

C=[C ₁C ₂C _kC _M] ^T(11) Elements C in the formula (11) _k(k=1,2,3 ...., M; M 〉=N) be Cable Structure, according to the current numerical value of the monitored amount of the pairing k of being numbered of coding rule.

d _c=[d _C1d _C2D _CjD _CN] ^T(12) d in the formula (12) _Cj(j=1,2,3 ...., N) be the current nominal virtual lesion value of cable system j root rope, vectorial d _cThe coding rule of subscript j of element identical with the coding rule of the row of matrix Δ C.

When the rope actual damage is not too big, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure is also less, the represented a kind of like this linear relationship of formula (9) or formula (10) is less with the error of actual conditions, error can be used error vector e (formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).

E=abs (Δ Cd _c-C+C _o) abs () is the function that takes absolute value in (13) formula (13), 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 directly find the solution and obtain " the vectorial d of current nominal virtual lesion according to formula (9) or formula (10) and " the vectorial C of current (actual measurement) numerical value of monitored amount " _c".If done like this, the vectorial d that obtains _cIn element in addition bigger negative value can appear, just negative damage, this obviously is irrational.Therefore obtain vectorial d _cAcceptable separating (promptly have reasonable error, but can determine the position and the virtual lesion degree thereof of virtual damaged cable more exactly) become a rational solution, available formula (14) is expressed this method.

Abs (Δ Cd _c-C+C _oAbs () is the function that takes absolute value in)≤g (14) formula (14), and vectorial g describes the reasonable deviation that departs from ideal linearity relation (formula (9) or formula (10)), is defined by formula (15).

G=[g ₁g ₂G _kG _M] ^T(15) g in the formula (15) _k(k=1,2,3 ...., M) maximum allowable offset of the ideal linearity relation that departs from shown in formula (9) or the formula (10) has been described.Vector g can be selected according to the error vector e tentative calculation of formula (13) definition.

At " the initial value vector C of monitored amount _o" (survey or calculate), " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " (calculating) and " the current numerical value vector C of monitored amount " be when (actual measurement obtains) is known; can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain " current nominal virtual lesion vector d _c" acceptable separating; the element of " current actual virtual lesion vector d " (formula (16) is seen in definition) can calculate according to formula (17) then; just obtained " current actual virtual lesion vector d "; thus can determine the position and the virtual lesion degree of virtual damaged cable by d; then according to position and the relax level of below the method for narration being determined slack line, just determined to need to adjust rope and the long adjustment amount of rope thereof of Suo Li.

D=[d ₁d ₂D _jD _N] ^T(16) d in the formula (16) _j(j=1,2,3 ...., N) the actual virtual lesion value of expression j root rope, formula (17), d are seen in its definition _jBeing that to represent that j root rope not damaged did not have at 0 o'clock lax, 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 in the coding rule of the element of vectorial d and the formula (1) _oThe coding rule of element identical.

d _j=1-(1-d _Oj) (1-d _Cj) the middle d of (17) formula (17) _Oj(j=1,2,3 ...., N) be vectorial d _oJ element, d _CjBe vectorial d _cJ element.

After narration has obtained the actual virtual lesion vector of Suo Dangqian d below, how to determine the position and the relax level of slack line.

If total N root support cable in the cable system, structure rope force data is described by the Suo Li of N root support cable.Available " initial rope force vector F _o" the initial Suo Li (formula (18) is seen in definition) of all support cables in the expression Cable Structure.Because the initial Suo Li that calculates gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial Suo Li reliably, in the narration of back, will represent this calculated value and measured value with prosign.

F _o=[F _O1F _O2F _OjF _ON] ^T(18) F in the formula (18) _o(j=1,2,3 ...., N) being the initial Suo Li of j root support cable in the Cable Structure, this element is according to the Suo Li of coding rule corresponding to the appointment support cable.Vector F _oIt is constant.Setting up the initial mechanical calculating benchmark model A of Cable Structure _oThe time used vectorial F _o

The current cable power (formula (19) is seen in definition) of all support cables in the Cable Structure that obtains with " current cable force vector F " expression actual measurement among the present invention.

F=[F ₁F ₂F _jF _N] ^T(19) F in the formula (19) _j(j=1,2,3 ...., N) be the current cable power of j root support cable in the Cable Structure.

Among the present invention, under support cable original state (not damaged, do not have lax), and support cable is when being in free state (free state refers to that Suo Li is 0, back with), and the length of support cable is called initial drift, with " initial drift vector l _o" the initial drift (formula (20) is seen in definition) of all support cables in the expression Cable Structure.

l _o=[l _O1l _O2L _OjL _ON] ^T(20) l in the formula (20) _Oj(j=1,2,3 ...., N) be the initial drift of j root support cable in the Cable Structure.Vector l _oBe constant, after when beginning, determining, just no longer change.

Among the present invention, with the current drift (formula (21) is seen in definition) of all support cables in " current drift vector l " expression Cable Structure.

L=[l ₁l ₂L _jL _N] ^T(21) l in the formula (21) _j(j=1,2,3 ...., N) be the current drift of j root support cable in the Cable Structure.

Among the present invention, represent the change amount (formula (22) and formula (23) are seen in definition) of the drift of all support cables in the Cable Structure with " drift changes vectorial Δ l " (or claiming support cable current relax level vector).

Δ l=[Δ l ₁Δ l ₂Δ l _jΔ l _N] ^T(22) Δ l in the formula (22) _j(j=1,2,3 ...., N) being the change amount of the drift of j root support cable in the current cable structure, formula (23), Δ l are seen in its definition _jBe not that 0 rope is a slack line, Δ l _jNumerical value be the slack of rope, and the current relax level of expression cable system j root support cable also is the long adjustment amount of rope of this rope when adjusting Suo Li.

Δl _j＝l _j-l _oj (23)

By slack line is carried out the relax level identification that slack line is carried out in the mechanics equivalence with damaged cable, the mechanical condition of equivalence is in the present invention:

The mechanics parameters of initial drift, geometrical property parameter and material when one, the nothing of the rope of two equivalences relaxes with not damaged is identical;

Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.

When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing slack line with the damaged cable of equivalence, Cable Structure any variation can not take place, vice versa.

Among the present invention, with j support cable (its current relax level Δ l _jDefinition) carries out the current actual virtual lesion degree d of virtual impaired support cable of equivalence _jExpression (d _jDefinition see formula (16) and formula (17)).The current relax level Δ l of j lax support cable _j(Δ l _jDefinition see formula (22)) with the current actual virtual lesion degree d of damaged cable of equivalence _jBetween relation determine by aforementioned two mechanics equivalent conditions.Δ l _jSame d _jBetween physical relationship can adopt accomplished in many ways, for example can directly determine (referring to formula (24)) according to aforementioned equivalent condition, also can adopt based on the Ernst equivalent elastic modulus to replace the E in the formula (24) to revise back definite (referring to formula (25)), also can adopt and determine based on other methods such as trial and error procedure of finite element method.

$\mathrm{\Δ}{l}_j=\frac{d_j}{1-d_j}\frac{F_j}{\mathrm{EA}+F_j}{l}_{\mathrm{oj}}---\left(24\right)$

$\mathrm{\Δ}{l}_i=\frac{d_i}{1-d_i}\frac{F_i}{\left[\frac{E}{1+\frac{{\left({\mathrm{\ω}}_i{l}_{\mathrm{ix}}\right)}^2\mathrm{AE}}{12{\left(F_i\right)}^3}}\right]A+F_i}{l}_{\mathrm{oi}}---\left(25\right)$

E is the elastic modulus of this support cable in formula (24) and the formula (25), and A is the cross-sectional area of this support cable, F _jBe the current cable power of this support cable, d _jBe the current actual virtual lesion degree of this support cable, ω is the weight of the unit length of this support cable, l _JxIt is the horizontal range of two supporting end points of this support cable.Item in the formula (25) in [] is the Ernst equivalent elastic modulus of this support cable, can just can determine the current relax level vector of support cable Δ l by formula (24) or formula (25).Formula (25) is the correction to formula (24).

Second portion of the present invention: the software and hardware part of health monitoring systems.Hardware components comprises monitoring system (monitoring the horizontal range of monitored amount, monitoring Cable Structure support coordinate, monitoring Suo Li, monitoring support cable two supporting end points), signal picker and computing machine etc.Requirement is monitored each monitored amount in real time or quasi real time, monitors the Suo Li of each support cable, is monitored the horizontal range that each support cable two supports end points.Software should the following function of tool: 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;

B. structure rope force data is just described by the Suo Li of N root support cable; For simplicity, in the present invention " the monitored rope force data of structure " abbreviated as " monitored amount "; When mentioning " monitored amount so-and-so matrix or so-and-so vector " in the back, also can be read as " Suo Li so-and-so matrix or so-and-so vector ";

C. the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up initial virtual lesion vector d _oIf when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, vectorial d ¹ _oEach element numerical value get 0.

D. setting up initial virtual lesion vector d _oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount _o

E. setting up initial virtual lesion vector d _oInitial value vector C with monitored amount _oThe time, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F _oSimultaneously, obtain the initial drift of all support cables, form initial drift vector l according to structural design data, completion data _oSimultaneously, obtain the initial geometric data of Cable Structure according to structural design data, completion data or actual measurement; Simultaneously, survey or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information;

F. set up the initial mechanical calculating benchmark model A of Cable Structure _o, set up initial Cable Structure support coordinate vector U _o, set up the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe measured data of the Cable Structure in Cable Structure completion, this measured data comprises measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, all ropes, 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 _oIf there is not the measured data of the structure in the Cable Structure completion, so just before setting up health monitoring systems, this Cable Structure is surveyed, obtain the measured data of Cable Structure equally, according to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure equally _oNo matter which kind of method to obtain A with _o, based on A _oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A _oCable Structure support coordinate data form initial Cable Structure support coordinate vector U _oA _oAnd U _oBe constant; For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A ^t _o", A in structure military service process ^t _oCan bring in constant renewal in as required, during beginning, A ^t _oEqual A _oEqually for sake of convenience, name " Cable Structure actual measurement support coordinate vector U ^t", in structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, and all Cable Structure support coordinate current datas are formed " current cable structure actual measurement support coordinate vector U ^t", vectorial U ^tElement and vectorial U _oThe coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, the last time is upgraded A ^t _oThe time Cable Structure support coordinate current data be designated as current cable structural bearings coordinate vector U ^t _oDuring beginning, A ^t _oEqual A _o, U ^t _oEqual U _oA _oThe health status of corresponding rope is by d _oDescribe;

When g. health monitoring systems is started working, make A ^t _oEqual A _oConstantly actual measurement obtains Cable Structure support coordinate current data in structure military service process, and all Cable Structure support coordinate current datas are formed current cable structure actual measurement support coordinate vector U ^t, according to current cable structure actual measurement support coordinate vector U ^t, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary ^t _oWith current cable structural bearings coordinate vector U ^t _o

H. at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _u

I. actual measurement obtains the current cable power of all support cables of Cable Structure, forms current cable force vector F; Simultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount ".Actual measurement calculates the volume coordinate of two supporting end points of all support cables, and the volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.When numbering to the element of the institute's directed quantity that occurred before this step and this step, should use same coding rule, each vector that can guarantee before this step and this step like this and occur afterwards, number identical element, represent same monitored amount, corresponding to vectorial defined relevant information under this element;

J. define current nominal virtual lesion vector d to be asked _cWith current actual virtual lesion vector d.Damage vectorial d _o, d _cEqualing the quantity of rope with the element number of d, 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 virtual lesion degree or the health status of corresponding rope;

K. according to " the current numerical value vector C of monitored amount " " the vectorial C of the initial value of monitored amount together _o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and " current nominal virtual lesion vector d _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 _cOther outer amount is known, finds the solution formula 1 and just can calculate current nominal virtual lesion vector d _c

C=C _o+ Δ Cd _cFormula 1

1. utilize the element d of the current actual virtual lesion vector d of formula 2 expression _jWith initial virtual lesion vector d _oElement d _OjWith current nominal virtual lesion vector d _cElement d _CjBetween relation, calculate all elements of current actual virtual lesion vector d.

d _j=1-(1-d _Oj) (1-d _Cj) formula 2

J=1 in the formula 2,2,3 ..., N.

Because the element numerical value of current actual virtual lesion vector d is represented the current actual virtual lesion degree of corresponding rope, be actual relax level or actual damage degree, numerical value is not that the support cable of 0 element correspondence is exactly problematic support cable among the current actual virtual lesion vector d, problematic support cable may be slack line, also may be damaged cable, its numerical response the degree of lax or damage;

M. identify damaged cable in the problematic support cable that identifies from the 1st step, remaining is exactly slack line.

N. utilize the current actual virtual lesion vector d that obtains in the 1st step to obtain the current actual virtual lesion degree of slack line, the current cable force vector F that utilization obtained in the i step, utilization utilizes the vectorial l of the initial drift that obtains in the e step in the volume coordinates of two supporting end points of all support cables of i step acquisition _oUtilization is in elastic modulus, density, the initial cross sectional area data of all ropes of e step acquisition, by with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual virtual lesion degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa.Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has just determined the long adjustment amount of rope of the support cable that those need adjust Suo Li.The lax identification and the damage identification of support cable have so just been realized.Institute's demand power is provided by current cable force vector F corresponding element during calculating.

In step g, according to current cable structure actual measurement support coordinate vector U ^t, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary ^t _oWith current cable structural bearings coordinate vector U ^t _oConcrete grammar be:

G1. actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tEqual U ^t _o, then do not need A ^t _oUpgrade;

G2. actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tBe not equal to U ^t _o, then need A ^t _oUpgrade, update method is: calculate U earlier ^tWith U _oPoor, U ^tWith U _oDifference be exactly that the current cable structural bearings is about setting up A _oThe 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 ^t _oMethod be: to A _oIn the Cable Structure bearing apply the constraint of current support displacement, the numerical value of current support displacement constraint is just taken from the numerical value of corresponding element among the current support displacement vector V, to A _oIn 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 ^t _o, upgrade A ^t _oThe time, U ^t _oAll elements numerical value is also used U ^tAll elements numerical value replaces, and has promptly upgraded U ^t _o, so just obtained correctly corresponding to A ^t _oU ^t _o

In step h, at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _uConcrete grammar be:

When h1. health monitoring systems was started working for the first time, directly h2 obtained Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of step h4 set by step _uAfter, if in the step g to A ^t _oUpgrade, directly h2 obtains Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of step h4 set by step _uIf, in step g not to A ^t _oUpgrade, then directly change step I herein over to and carry out follow-up work;

H2. at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, there is N root rope that N calculating is just arranged, calculating each time in the hypothesis cable system has only a rope to increase virtual unit damage again on the basis of original virtual lesion, the rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating, and supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector D _u" write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is 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 _Tj" the current evaluation vector C of monitored amount that expression is corresponding _Tj".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 evaluation vector C of monitored amount _Tj" deduct " initial value of monitored amount vector C _o" obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; There is N root rope that N " the numerical value change vector of monitored amount " just arranged;

H4. form " the monitored numerical quantity transformation matrices of the virtual unit damage Δ C " that the N row are arranged successively by this N " the numerical value change vector of monitored amount "; Each row of " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " are corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector d _cIdentical with the element coding rule of current actual virtual lesion vector d.

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 or lax condition of more rope monitoring and evaluation very exactly go out the health status (position and relax level or the degree of injury that comprise all slack lines and damaged cable) of cable system.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 present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status (relax level and the extent of damage that comprise rope) of the cable system of Cable Structure.During 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.

Cable Structure rope force data is just described by the Suo Li of N root support cable.For simplicity, when in the present invention " the monitored rope force data of structure " being abbreviated as " monitored amount " and mentioning " monitored amount so-and-so matrix or so-and-so vector " in the back, also can be read as " Suo Li so-and-so matrix or so-and-so vector ".

Second step: the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up initial virtual lesion vector d _oIf 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 not damaged, no relaxed state vectorial d _oEach element numerical value get 0.

The 3rd step: setting up initial virtual lesion vector d _oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms " the initial value vector C of monitored amount _o"; Simultaneously, directly measure the initial Suo Li of all support cables that calculate Cable Structure, form " initial rope force vector F _o"; Simultaneously, obtain the initial drift of all ropes, form " the initial drift vector of support cable l according to structural design data, completion data _o"; Simultaneously, survey or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information.

The 4th step: setting up initial virtual lesion vector d ¹ _oThe time, can adopt ripe measuring method to carry out cable force measurement, strain measurement, measurement of angle and volume coordinate and measure.Calculate Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) after directly measuring or measuring, the original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is to 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, set up the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe measured data of the Cable Structure in Cable Structure completion, this measured data comprises measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, all ropes, 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 _oIf there is not the measured data of the structure in the Cable Structure completion, so just before setting up health monitoring systems, this Cable Structure is surveyed, obtain the measured data of Cable Structure equally, according to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure equally _oNo matter which kind of method to obtain A with _o, based on A _oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A _oCable Structure support coordinate data form initial Cable Structure support coordinate vector U _oA _oAnd U _oBe constant; For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A ^t _o", A in structure military service process ^t _oCan bring in constant renewal in as required, during beginning, A ^t _oEqual A _oEqually for sake of convenience, name " Cable Structure actual measurement support coordinate vector U ^t", in structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, and all Cable Structure support coordinate current datas are formed " current cable structure actual measurement support coordinate vector U ^t", vectorial U ^tElement and vectorial U _oThe coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, the last time is upgraded A ^t _oThe time Cable Structure support coordinate current data be designated as current cable structural bearings coordinate vector U ^t _oDuring beginning, A ^t _oEqual A _o, U ^t _oEqual U _oA _oThe health status of corresponding rope is by d _oDescribe;

The 5th step: the hardware components of pass line structural healthy monitoring system.Hardware components comprises at least: horizontal range monitoring system (for example measuring with total powerstation), signal (data) collector, the computing machine and the panalarm of communicating by letter of monitored amount monitoring system (for example containing acceleration transducer, signal conditioner etc.), cable force monitoring system (for example containing acceleration transducer, signal conditioner etc.), each support cable two supporting end points.The horizontal range of the Suo Li of each monitored amount, each support cable and each root support cable two supporting end points all must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector with the signal that monitors; Signal is delivered to computing machine through signal picker; Computing machine then is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; Have lax or during damage when monitoring rope, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.

The 6th step: establishment and the cable system health monitoring systems software of installation and operation Cable Structure on supervisory control comuter.This software will be finished functions such as monitoring that the present invention's method of the lax support cable of the identification of cable force 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: when health monitoring systems is started working, make A ^t _oEqual A _oConstantly actual measurement (for example measuring with total powerstation) obtains Cable Structure support coordinate current data in structure military service process, and all Cable Structure support coordinate current datas are formed current cable structure actual measurement support coordinate vector U ^t, according to current cable structure actual measurement support coordinate vector U ^t, upgrade the current Mechanics Calculation benchmark model of Cable Structure A where necessary ^t _oWith current cable structural bearings coordinate vector U ^t _oConcrete grammar is: actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tEqual U ^t _o, then do not need A ^t _oUpgrade; Actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tBe not equal to U ^t _o, then need A ^t _oUpgrade.Upgrade A ^t _oMethod be: calculate earlier U ^tWith U _oPoor, U ^tWith U _oDifference be exactly that the current cable structural bearings is about setting up A _oThe 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 ^t _oMethod be: to A _oIn the Cable Structure bearing apply the constraint of current support displacement, the numerical value of current support displacement constraint is just taken from the numerical value of corresponding element among the current support displacement vector V, to A _oIn 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 ^t _o, upgrade A ^t _oThe time, U ^t _oAll elements numerical value is also used U ^tAll elements numerical value replaces, and has promptly upgraded U ^t _o, so just obtained correctly corresponding to A ^t _oU ^t _o

The 8th step: at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _uConcrete grammar is as follows:

When a. health monitoring systems was started working for the first time, directly b obtained Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of steps d set by step _uAfter, if the 7th the step in to A ^t _oUpgrade, directly b obtains Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of steps d set by step _uIf, the 7th the step in not to A ^t _oUpgrade, then directly changing for the 9th step herein over to carries out follow-up work;

B. at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, there is N root rope that N calculating is just arranged, calculating each time in the hypothesis cable system has only a rope to increase virtual unit damage again on the basis of original virtual lesion, the rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating, and supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector D _u" write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is 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 _Tj" the current evaluation vector C of monitored amount that expression is corresponding _Tj".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 evaluation vector C of monitored amount _Tj" deduct " initial value of monitored amount vector C _o" obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; There is N root rope that N " the numerical value change vector of monitored amount " just arranged;

D. form " the monitored numerical quantity transformation matrices of the virtual unit damage Δ C " that the N row are arranged successively by this N " the numerical value change vector of monitored amount "; Each row of " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " are corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector d _cIdentical with the element coding rule of current actual virtual lesion vector d.

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 and vectorial g.Utilize data (" the initial value vector C of monitored amount of front _o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C "); when the 8th step calculated each time; promptly in calculating each time, have only in the hypothesis cable system increase virtual unit damage again on the basis of rope at original virtual lesion in, calculate each time and form " virtual lesion vector a d _t", virtual lesion vector d _tElement number equal the quantity of rope, virtual lesion vector d _tAll elements in have only the numerical value of an element to get to calculate each time in hypothesis increase the virtual unit damage value of the rope of virtual unit damage, d _tThe numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the rope of virtual unit damage, be identical with the element of the same numbering of other vectors with the corresponding relation of this rope; With C _Tj, C _o, Δ C, d _tBring formula (13) into and (note C C in the formula (13) _TjBring d into _cUse d _tBring into), obtain a linear relationship error vector e, calculate a linear relationship error vector e each time; Have N root rope that N calculating is just arranged, N linear relationship error vector e just arranged, will obtain a vector after this N the linear relationship error vector e addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e.Vector g equals final error vector e.Vectorial g is kept on the hard disc of computer of operation health monitoring systems software, uses for health monitoring systems software.

Will " initial rope force vector F _o", " initial value of monitored amount vector C _o", " nominal virtual unit damage vector D _u", " initial drift vector l _o", the parameters such as unit weight of the elastic modulus of " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and all ropes, initial cross sectional area, rope are kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.

The tenth step: actual measurement obtains the current cable power of all support cables of Cable Structure, forms current cable force vector F; Simultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount ".Actual measurement calculates the volume coordinate of two supporting end points of all support cables, and the volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.

The 11 step: foundation " current (calculating or actual measurement) numerical value vector C of monitored amount " is " the initial value vector C of monitored amount together _o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and " current nominal virtual lesion vector d _c" between the linear approximate relationship (seeing formula (9)) that exists, calculate the current nominal virtual lesion vector d of cable system according to multi-objective optimization algorithm _cNoninferior 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 (Constran Method), weighted method (Weghted Sum Method), goal programming method (Goal Attanment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution current nominal virtual lesion vector d with the goal programming method only _cProcess, 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 (29) and the formula (30), and γ is a real number in the formula (29), and R is a real number field, and area of space Ω has limited vectorial d _cSpan (the present embodiment requirements vector d of each element _cEach element be not less than 0, be not more than 1).The meaning of formula (29) is to seek the real number γ of an absolute value minimum, makes formula (30) be met.G (d in the formula (30) _c) by formula (31) definition, G (d in the product representation formula (30) of weighing vector W and γ in the formula (30) _c) and vectorial g between the deviation that allows, the definition of g is referring to formula (15), its value will calculate in the 8th step.Vector W can be identical with vectorial g during actual computation.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of the vectorial d of current name damage _c

minimize γ

(29)

γ∈R，d _c∈Ω

G(d _c)-Wγ≤g (30)

G(d _c)＝abs(ΔC·d _c-C+C _o) (31)

Try to achieve current nominal virtual lesion vector d _cAfter, can be according to vectorial each element of d of the current actual virtual lesion that formula (17) obtain, current actual virtual lesion vector d has reasonable error exactly but can determine the position of problematic rope (be virtual damaged cable, may be impaired also may be lax) more exactly and separating of virtual lesion degree from all ropes.If the numerical value of a certain element of the current actual virtual lesion that solves vector d is 0, represent that the pairing rope of this element is intact, not damage or lax; If its numerical value is 100%, represent that then the pairing rope of this element has completely lost load-bearing capacity; If its numerical value between 0 and 100%, is then represented this rope and has been lost the load-bearing capacity of corresponding proportion.

The 12 step: because the element numerical value of current actual virtual lesion vector d is represented the virtual lesion degree of corresponding rope, so it is impaired or relaxed and possible degree of injury or relax level just to define which Suo Keneng according to current actual virtual lesion vector, but damage has taken place actually or has taken place to relax in these ropes, need differentiate.The method of differentiating is varied; can be by removing the protective seam of support cable; support cable is carried out visual discriminating; perhaps carry out visual discriminating by optical imaging apparatus; also can be by lossless detection method to support cable impaired discriminating the whether, UT (Ultrasonic Testing) is exactly a kind of present widely used lossless detection method.Those do not find to damage and the virtual lesion degree is not that 0 support cable is exactly that lax rope has taken place to differentiate the back, need adjust the rope of Suo Li exactly, can be in the hope of the relax level (being the long adjustment amount of rope) of these ropes according to formula (24) or formula (25).So just realized comprising the health monitoring of the cable system of damage identification and the lax Cable Structure of discerning.

Claims (2)

1. One kind when support settlement is arranged based on the method for the lax support cable of the identification of cable force monitoring, it is characterized in that described method 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;

   B. structure rope force data is just described by the Suo Li of N root support cable; For simplicity, " the monitored rope force data of structure " abbreviated as " monitored amount "; When mentioning " monitored amount so-and-so matrix or so-and-so vector " in the back, also can be read as " Suo Li so-and-so matrix or so-and-so vector ";

   C. the data of utilizing the Non-Destructive Testing data etc. of rope can express the health status of rope are set up initial virtual lesion vector d _oIf when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, vectorial d _oEach element numerical value get 0;

   D. setting up initial virtual lesion vector d _oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount _o

   E. setting up initial virtual lesion vector d _oInitial value vector C with monitored amount _oThe time, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F _oSimultaneously, obtain the initial drift of all support cables, form initial drift vector l according to structural design data, completion data _oSimultaneously, obtain the initial geometric data of Cable Structure according to structural design data, completion data or actual measurement; Simultaneously, survey or obtain elastic modulus, density, the initial cross sectional area of all ropes according to structural design, completion information;

   F. set up the initial mechanical calculating benchmark model A of Cable Structure _o, set up initial Cable Structure support coordinate vector U _o, set up the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe measured data of the Cable Structure in Cable Structure completion, this measured data comprises measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, all ropes, and the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope, according to design drawing and as-constructed drawing, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure _oIf there is not the measured data of the structure in the Cable Structure completion, so just before setting up health monitoring systems, this Cable Structure is surveyed, obtain the measured data of Cable Structure equally, according to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure equally _oNo matter which kind of method to obtain A with _o, based on A _oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A _oCable Structure support coordinate data form initial Cable Structure support coordinate vector U _oA _oAnd U _oBe constant; For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A ^t _o", A in structure military service process ^t _oCan bring in constant renewal in as required, during beginning, A ^t _oEqual A _oEqually for sake of convenience, name " Cable Structure actual measurement support coordinate vector U ^t", in structure military service process, constantly actual measurement obtains Cable Structure support coordinate current data, and all Cable Structure support coordinate current datas are formed " current cable structure actual measurement support coordinate vector U ^t", vectorial U ^tElement and vectorial U _oThe coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, the last time is upgraded A ^t _oThe time Cable Structure support coordinate current data be designated as current cable structural bearings coordinate vector U ^t _oDuring beginning, A ^t _oEqual A _o, U ^t _oEqual U _oA _oThe health status of corresponding rope is by d _oDescribe;

   When g. health monitoring systems is started working, make A ^t _oEqual A _oConstantly actual measurement obtains Cable Structure support coordinate current data in structure military service process, and all Cable Structure support coordinate current datas are formed current cable structure actual measurement support coordinate vector U ^t, according to current cable structure actual measurement support coordinate vector U ^t, upgrade the current Mechanics Calculation benchmark model of Cable Structure A according to step g 1 to g2 ^t _oWith current cable structural bearings coordinate vector U ^t _o

   G1. actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tEqual U ^t _o, then do not need A ^t _oUpgrade;

   G2. actual measurement obtains current cable structure actual measurement support coordinate vector U ^tAfter, compare U ^tAnd U ^t _oIf, U ^tBe not equal to U ^t _o, then need A ^t _oUpgrade, update method is: calculate U earlier ^tWith U _oPoor, U ^tWith U _oDifference be exactly that the current cable structural bearings is about setting up A _oThe 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 ^t _oMethod be: to A _oIn the Cable Structure bearing apply the constraint of current support displacement, the numerical value of current support displacement constraint is just taken from the numerical value of corresponding element among the current support displacement vector V, to A _oIn 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 ^t _o, upgrade A ^t _oThe time, U ^t _oAll elements numerical value is also used U ^tAll elements numerical value replaces, and has promptly upgraded U ^t _o, so just obtained correctly corresponding to A ^t _oU ^t _o

   H. at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _u

   I. actual measurement obtains the current cable power of all support cables of Cable Structure, forms current cable force vector F; Simultaneously, actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount "; Actual measurement calculates the volume coordinate of two supporting end points of all support cables, and the volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges; When numbering to the element of the institute's directed quantity that occurred before this step and this step, should use same coding rule, each vector that can guarantee before this step and this step like this and occur afterwards, number identical element, represent same monitored amount, corresponding to vectorial defined relevant information under this element;

   J. define current nominal virtual lesion vector d to be asked _cWith current actual virtual lesion vector d; Damage vectorial d _o, d _cEqualing the quantity of rope with the element number of d, 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 virtual lesion degree or the health status of corresponding rope;

   K. according to " the current numerical value vector C of monitored amount " " the vectorial C of the initial value of monitored amount together _o", " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " and " current nominal virtual lesion vector d _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 _cOther outer amount is known, finds the solution formula 1 and just can calculate current nominal virtual lesion vector d _c

   C=C _o+ Δ Cd _cFormula 1

   L. utilize the element d of the current actual virtual lesion vector d of formula 2 expression _iWith initial virtual lesion vector d _oElement d _OjWith current nominal virtual lesion vector d _cElement d _CjBetween relation, calculate all elements of current actual virtual lesion vector d;

   d _j=1-(1-d _Oj) (1-d _Cj) formula 2

   J=1 in the formula 2,2,3 ..., N;

   Because the element numerical value of current actual virtual lesion vector d is represented the current actual virtual lesion degree of corresponding rope, be actual relax level or actual damage degree, numerical value is not that the support cable of 0 element correspondence is exactly problematic support cable among the current actual virtual lesion vector d, problematic support cable may be slack line, also may be damaged cable, its numerical response the degree of lax or damage;

   M. identify damaged cable in the problematic support cable that identifies from the 1st step, remaining is exactly slack line;

   N. utilize the current actual virtual lesion vector d that obtains in the 1st step to obtain the current actual virtual lesion degree of slack line, the current cable force vector F that utilization obtained in the i step, utilization utilizes the vectorial l of the initial drift that obtains in the e step in the volume coordinates of two supporting end points of all support cables of i step acquisition _oUtilization is in elastic modulus, density, the initial cross sectional area data of all ropes of e step acquisition, by with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual virtual lesion degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, just determine the long adjustment amount of rope of the support cable that those need adjust Suo Li, so just realized the lax identification and the damage identification of support cable; Institute's demand power is provided by current cable force vector F corresponding element during calculating.
2. 2. according to claim 1 when support settlement is arranged based on the method for the lax support cable of the identification of cable force monitoring, it is characterized in that in step h, at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure virtual unit damage monitored numerical quantity transformation matrices Δ C and nominal virtual unit damage vector D _uConcrete grammar be:

   When h1. health monitoring systems was started working for the first time, directly h2 obtained Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of step h4 set by step _uAfter, if in the step g to A ^t _oUpgrade, directly h2 obtains Cable Structure virtual unit damage monitored quantitative change matrix Δ C and nominal virtual unit damage vector D to the listed method of step h4 set by step _uIf, in step g not to A ^t _oUpgrade, then directly change step I herein over to and carry out follow-up work;

   H2. at the current Mechanics Calculation benchmark model of Cable Structure A ^t _oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, there is N root rope that N calculating is just arranged, calculating each time in the hypothesis cable system has only a rope to increase virtual unit damage again on the basis of original virtual lesion, the rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other time calculating, and supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector D _u" write down the unit damage of the supposition of all ropes, calculate the current numerical value of all monitored amounts each time, the current numerical value of the monitored amount of all that calculate is 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 _Tj" the current evaluation vector C of monitored amount that expression is corresponding _Tj"; When giving each vectorial element numbering in this step, should use same coding rule with other vector, can guarantee any one element in each vector in this step like this,, express the relevant information of same monitored amount or same target with element in other vector, that numbering is identical;

   H3. that calculates each time " the current evaluation vector C of monitored amount _Tj" deduct " initial value of monitored amount vector C _o" obtain a vector, during all calculating divided by this, each element that again should vector obtains one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition; There is N root rope that N " the numerical value change vector of monitored amount " just arranged;

   H4. form " the monitored numerical quantity transformation matrices of the virtual unit damage Δ C " that the N row are arranged successively by this N " the numerical value change vector of monitored amount "; Each row of " the monitored numerical quantity transformation matrices of virtual unit damage Δ C " are corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector d _cIdentical with the element coding rule of current actual virtual lesion vector d.