CN102323084A - Progressive method for identifying damaged cable and support generalized displacement on the basis of mixed monitoring - Google Patents

Abstract

The invention provides a progressive method for identifying a damaged cable and a support generalized displacement on the basis of mixed monitoring. Taking into consideration the fact that linear relationships between the current value vector of the monitored amount and an initial value vector of the monitored amount, a unit damage monitored amount variation matrix and the current nominal damage vector are approximate with each other, the invention provides a method for approaching nonlinear relation in a sectional way by using the linear relationships in order to overcome the defects. A big interval is divided into continuous small intervals; the linear relation is accurate enough within each small interval; and the support generalized displacement and the damaged cable can be rapidly identified within each small interval by using a multi-target optimization algorithm and other appropriate algorithms.

Description

Progressive method based on the generalized displacement of hybrid monitoring identification damaged cable bearing

Technical field

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 comprises the rod member that only bears tensile load) can be impaired, the bearing of Cable Structure generalized displacement also possibly occur (for example the bearing generalized displacement refers to that bearing is along the angular displacement around X, Y, Z axle of the displacement of the lines of X, Y, Z axle and bearing simultaneously; Corresponding to the bearing generalized displacement; The bearing generalized coordinate refers to that bearing is about the coordinate of X, Y, Z axle and the bearing angular coordinate about X, Y, Z axle); These safety that change to Cable Structure are a kind of threats; The present invention is based on structural health monitoring technology, based on hybrid monitoring, adopt progressive method to discern the damaged cable in the cable system of bearing generalized displacement and Cable Structure, belong to the engineering structure health monitoring field.

Background technology

The bearing generalized displacement is a significant threat to Cable Structure safety; Same; Cable system is the key components of Cable Structure normally; Its inefficacy usually brings the inefficacy of total, and the damaged cable of discerning based on structural health monitoring technology in the cable system of bearing generalized displacement and Cable Structure is a kind of method that has potentiality.When generalized displacement appears in bearing or the health status of cable system when changing, or two kinds of situation when taking place simultaneously; Can cause the variation of the measurable parameter of structure; For example can cause the variation of Suo Li; Can influence the distortion or the strain of Cable Structure; Can influence the shape or the volume coordinate of Cable Structure; Can cause the variation (the for example arbitrarily variation of the angle coordinate of the straight line of any this point of mistake in any section of body structure surface, the perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily) of angle coordinate of any imaginary line of the every bit of Cable Structure, all these change the health status information that has all comprised cable system; In fact the variation of these measurable parameters comprised cable system health status information, comprised bearing generalized displacement information; Therefore can be through the hybrid monitoring of the variation of the characteristic parameter of these dissimilar structures being judged the health status of structure, the present invention's architectural feature parameter that all are monitored is referred to as " monitored amount ", because monitored amount was mixed by the dissimilar measurable parameter of structure and formed this moment; The present invention claims that this is a hybrid monitoring, that is to say and can utilize hybrid monitoring to discern bearing generalized displacement and damaged cable.

In order reliable monitoring and judgement to be arranged to the health status and the bearing generalized displacement of the cable system of Cable Structure; The method of the variation of a measurable parameter of can rational and effective setting up Cable Structure with the relation between the health status of all ropes in bearing generalized displacement and the cable system must be arranged, and the health monitoring systems of setting up based on this method can provide the health evaluating of more believable bearing generalized displacement assessment and cable system.

Summary of the invention

Technical matters:The invention discloses a kind of based on to the hybrid monitoring of multiclass parameter, adopt health monitor method progressive method, that can discern bearing generalized displacement and damaged cable rationally and effectively.

Technical scheme:If the quantity sum of the quantity of rope and bearing generalized displacement component does NFor the purpose of narrating conveniently, the present invention is unified to claim that rope and the bearing generalized displacement assessed are " by evaluation object ", gives by the evaluation object serial number, and the present invention is with using variable jRepresent this numbering, j=1,2,3 ..., N, therefore we can say NIndividual 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 to the hybrid monitoring data of 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=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 _oBe constant), set up the Mechanics Calculation benchmark model A of Cable Structure ⁱ(finite element benchmark model for example, i=1,2,3 ...).Letter i is except the place of representing number of steps significantly, and is alphabetical in the present invention iOnly represent cycle index, promptly iInferior circulation.

The iThe Cable Structure that needs when inferior circulation begins " initial health vector d _o ⁱ " (shown in (1)), use d _o ⁱ Expression the iCable Structure when inferior circulation begins is (with Mechanics Calculation benchmark model A ⁱThe initial health of Cable Structure expression).

(1)

In the formula (1) d ⁱ _Oj ( i=1,2,3, ; j=1,2,3 ...., N) expression circulates for the i time when beginning, Mechanics Calculation benchmark model A ⁱIn cable system jIndividual by the current health status of evaluation object, if should be the rope (or pull bar) in the cable system, so by evaluation object d _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 a generalized displacement component of a bearing by evaluation object, so d _i Represent its current generalized displacement numerical value.In the formula (1) TThe transposition of expression vector (back together).

Setting up the initial health vector during circulation beginning for the first time (is designated as according to formula (1) d ¹ _o ) time, the data and the bearing generalized 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 ¹ _o If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vector d ¹ _o In get 0 with each element numerical value of Suo Xiangguan.

The i time ( i=2,3,4,5,6 ...) needs was vectorial by the evaluation object initial health when circulation began d ⁱ _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 ⁱ

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

The Mechanics Calculation benchmark model of the Cable Structure of setting up during circulation beginning for the first time is designated as A ¹, A ¹Just equal A _oA ¹Corresponding by the health status of evaluation object by d ¹ _o Describe.

The i time ( i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins ⁱ, 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 ¹With vectorial by the evaluation object initial health d ¹ _o After, model A ¹In each by the health status of evaluation object by vector d ¹ _o Express.At A ¹The basis on, all are changed to 0 by the health status numerical value of evaluation object, mechanical model A ¹Be updated to one all be that 0 mechanical model (is designated as A by the health status of evaluation object ⁰), mechanical model A ⁰Be actually the corresponding mechanical model of Cable Structure of excellent no bearing generalized displacement.Might as well claim model A ⁰For the not damaged of Cable Structure does not have bearing generalized displacement model A ⁰

Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, and narration is as follows respectively:

If it is total in the cable system QThe root rope, the monitored rope force data of structure is by on the structure M ₁ Individual appointment rope M ₁ Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope.Each total M ₁ Individual cable force measurement value or calculated value characterize the rope force information of structure. M ₁ Be one and be not less than 0 integer.

The monitored strain data of structure can be by on the structure K ₂ Individual specified point, and each specified point L ₂ The strain of individual assigned direction is described, and the variation of structural strain data is exactly K ₂ The variation of the tested strain of all of individual specified point.Each total M ₂ (M ₂ =K ₂ * L ₂ )Individual strain measurement value or calculated value characterize structural strain. M ₂ Be one and be not less than 0 integer.

The monitored angle-data of structure is by on the structure K ₃ Individual specified point, cross each specified point L ₃ Individual appointment straight line, each specifies straight line H ₃ Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M ₃ (M ₃ =K ₃ * L ₃ * H ₃ )Individual angle coordinate component measurement value or calculated value characterize the angle information of structure. M ₃ Be one and be not less than 0 integer.

The monitored shape data of structure is by on the structure K ₄ Individual specified point, and each specified point L ₄ The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K ₄ The variation of all coordinate components of individual specified point.Each total M ₄ (M ₄ =K ₄ * L ₄ )Individual measurement of coordinates value or calculated value characterize planform. M ₄ Be one and be not less than 0 integer.

Comprehensive above-mentioned monitored amount, total is total M (M=M ₁ + M ₂ + M ₃ + M ₄ )Individual monitored amount, the definition parameter K( K=M ₁ + K ₂ + K ₃ + K ₄ ), KWith MMust not be less than by the quantity of evaluation object NBecause MIndividual monitored amount is dissimilar, so the present invention is called " based on the progressive method of hybrid monitoring identification damaged cable bearing generalized displacement ".

For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " monitored amount ".

Vectorial among the present invention with monitored amount initial value C ⁱ _o " ( i=1,2,3 ...) expression the i time ( i=1,2,3,4,5,6 ...) initial value (referring to formula (2)) of the monitored amount of all appointments when circulation begins, C ⁱ _o Full name be " the initial value vector of the i time monitored amount of circulation ".

(2)

In the formula (2) C ⁱ _Ok ( i=1,2,3, k=1,2,3, ., M; M>=N;) be in the i time when beginning circulation, the Cable Structure the kIndividual monitored amount.Vector C ⁱ _o Define by the front MIndividual monitored amount forms according to certain series arrangement, and this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.

During circulation beginning for the first time, " the initial value vector of the 1st the monitored amount that circulates C ¹ _o " (seeing formula (2)) be made up of measured data, because according to model A ¹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 need " initial value of the i time monitored amount of circulation is vectorial C ⁱ _o ", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.

Second step: 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 ⁱ , " nominal unit damage vector " that the i time circulation set up is designated as D ⁱ _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 ¹ Set up Δ C ¹ Process following:

Mechanics Calculation benchmark model A in Cable Structure ¹The basis on carry out several times and calculate, equal on the calculation times numerical value NCalculating 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 by evaluation object be the generalized displacement component of a direction of a bearing, just suppose that this bearing is this generalized displacement direction generation unit generalized displacement (for example get 1 millimeter, 2 millimeters, 3 millimeters etc. and be the unit line displacement, get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. and be the unit angular displacement).For narrating conveniently, the damage and the bearing generalized 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 ⁱ _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 ⁱ _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 for the first time D ¹ _u Calculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, appointment in front MThe current calculated value of individual monitored amount calculates gained each time MThe current calculated value of individual monitored amount is formed one " the current numerical value vector of the calculating of monitored amount " (when hypothesis the jIndividual when by evaluation object unit damage being arranged, available formula (4) is represented all appointments MThe current numerical value vector of the calculating of individual monitored amount C ¹ _Tj ); The current numerical value vector of the calculating of the monitored amount that calculates each time deducts the initial value vector of monitored amount C ¹ _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 the jIndividual when unit damage being arranged by evaluation object, use δ C ¹ _j The numerical value change vector of representing monitored amount, δ C ¹ _j Definition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vector for formula (4) D ¹ _u jIndividual element D _Uj Gained), the numerical value change of monitored amount vector δ C ¹ _j Each element representation since that supposition has a unit damage when calculating by evaluation object (for example the jIndividual by evaluation object) there is unit damage (for example D _Uj ), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage of supposition D _Uj Rate of change; Have NIndividual just had by evaluation object NIndividual " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has M(it is general, M>=N) individual element, by this NIndividual " the numerical value change vector of monitored amount " formed successively to be had M * N" the unit damage monitored numerical quantity transformation matrices of individual element Δ C ¹ " (MOK NRow ), each vector δ C ¹ _j ( j=1,2,3 ...., N) be matrix Δ C ¹ One row, Δ C ¹ Definition suc as formula shown in (8).

(3)

Nominal unit damage vector in the formula (3) D ⁱ _u Element D ⁱ _Uj ( i=1,2,3, j=1,2,3 ...., N) expression the iSuppose in the inferior circulation jIndividual by the unit damage numerical value of evaluation object, vector D ⁱ _u In the numerical value of each element can be the same or different.

(4)

Element in the formula (4) C ⁱ _Tjk ( i=1,2,3, j=1,2,3 ...., N; k=1,2,3 ...., M; M>=N) expression the iInferior circulation is because the jIndividual when unit damage being arranged by evaluation object, according to coding rule pairing kThe current numerical value of the calculating of the monitored amount of individual appointment.

(5)

The subscript of each amount in the formula (5) i( i=1,2,3 ...) expression the iInferior circulation, subscript j( j=1,2,3 ...., N) expression the jIndividual had unit damage by evaluation object, in the formula D ⁱ _Uj It is vector D ⁱ _u In jIndividual element.Vector δ C ⁱ _j Definition suc as formula shown in (6), δ C ⁱ _j k( k=1,2,3 ...., M; M>=N) individual element δ C ⁱ _Jk Expression the iIn the inferior circulation, set up matrix Δ C ⁱ The time, suppose jThe individual gained that calculates when unit damage being arranged by evaluation object kThe change amount of individual monitored amount is with respect to the unit damage of supposition D ⁱ _Uj Rate of change, it defines suc as formula shown in (7).

(6)

(7)

The definition of each amount was narrated in front in the formula (7).

(8)

Vector in the formula (8) δ C ⁱ _j ( i=1,2,3 ....,, j=1,2,3 ...., N) expression the iIn the inferior circulation, because the jIndividual had unit damage by evaluation object D ⁱ _Uj Cause, the relative value of all monitored amounts changes.Matrix Δ C ⁱ Row (subscript j) coding rule and front vector d ⁱ _o The subscript of element jCoding rule identical.

The 3rd step: identification is by the current health status of evaluation object.Detailed process is following.

The i( i=1,2,3 ...) in the inferior circulation, " current (calculating or actual measurement) numerical value vector of monitored amount C ⁱ " " initial value of monitored amount is vectorial together C ⁱ _o ", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ " and " current name damage vector d ⁱ _c " between linear approximate relationship, shown in (9) or formula (10).

(9)

(10)

Current (calculating or actual measurement) numerical value vector of monitored amount in formula (9) and the formula (10) C ⁱ Definition be similar to the initial value vector of monitored amount C ⁱ _o Definition, see formula (11); By the current name damage of evaluation object vector d ⁱ _c Definition see formula (12).

(11)

Element in the formula (11) C ⁱ _k ( i=1,2,3 ....; k=1,2,3 ...., M; M>=N) be iInferior circulation time Cable Structure, be numbered according to coding rule is pairing kThe current numerical value of monitored amount.

(12)

In the formula (12) d ⁱ _Cj ( i=1,2,3 ....; j=1,2,3 ...., N) be iCable Structure in the inferior circulation jIndividual by the current nominal impairment value of evaluation object, vector d ⁱ _c The subscript of element jCoding rule and matrix Δ C ⁱ The coding rule of row identical.

When by evaluation object actual damage or bearing generalized 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 ⁱ (formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).

(13)

In the formula (13) Abs ()Be the function that takes absolute value, each element of the vector of trying to achieve in the bracket is taken absolute value.

Because there are certain error in formula (9) or the represented linear relationship of formula (10), therefore can not be simply according to formula (9) or formula (10) and " current (actual measurement) numerical value vector of monitored amount C ⁱ " come directly to find the solution to obtain current name damage vector d ⁱ _c And it is vectorial to obtain current name damage d ⁱ _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 bearing generalized displacement amount) become a rational solution, available formula (14) is expressed this method.

(14)

In the formula (14) Abs ()Be the function that takes absolute value, vector g ⁱ Description departs from ideal linearity relation (formula (9) or formula (10))

Reasonable deviation, define by formula (15).

(15)

In the formula (15) g ⁱ _k ( i=1,2,3 ....; k=1,2,3 ...., M) described iThe maximum allowable offset that departs from the ideal linearity relation shown in formula (9) or the formula (10) in the inferior circulation.Vector g ⁱ Can be according to the error vector of formula (13) definition e ⁱ Tentative calculation is selected.

Initial value vector in monitored amount C ⁱ _o (survey or calculate), unit damage monitored numerical quantity transformation matrices Δ C ⁱ The current numerical value vector of (calculating) and monitored amount C ⁱ When (actual measurement obtains) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain current name damage vector d ⁱ _c Acceptable separating, current actual damage vector d ⁱ 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 ⁱ Thereby, can by d ⁱ Confirm damaged cable position and degree of injury, confirm bearing generalized displacement amount, realized that just damage identification and bearing generalized displacement discern.

(16)

In the formula (16) d ⁱ _j ( i=1,2,3, ; j=1,2,3 ...., N) expression the iIn the inferior circulation jIndividual by the actual damage value of evaluation object, formula (17) is seen in its definition, if should be the rope (or pull bar) in the cable system by evaluation object, so d ⁱ _j Represent its current damage, d ⁱ _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 a generalized displacement component of a bearing by evaluation object, so d ⁱ _j Represent its current generalized displacement numerical value.Vector d ⁱ Coding rule and the formula (1) of element in vector d ⁱ _o The coding rule of element identical.

(17)

In the formula (17) d ⁱ _Oj ( i=1,2,3,4, ; j=1,2,3 ...., N) be vector d ⁱ _o jIndividual element, d ⁱ _Cj It is vector d ⁱ _c jIndividual element.

The 4th step: judge whether to finish this (the iInferior) circulation, if then accomplish the tailing in work before this loop ends, for next time (promptly the i+ 1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following.

This ( iInferior) try to achieve current name damage vector in the circulation d ⁱ _c After, at first, set up mark vector according to formula (18) F ⁱ , formula (19) has provided mark vector F ⁱ jThe definition of individual element; If mark vector F ⁱ Element be 0 entirely, then in this circulation, continue health monitoring and calculating to Cable Structure; If mark vector F ⁱ Element be not 0 entirely, then accomplish subsequent step after, get into circulation next time.So-called subsequent step is: at first, calculate next time (promptly according to formula (20) i+ 1 time, i=1,2,3,4 ...) the required initial damage vector of circulation d ^{i+ 1} _o Each element d ^{i+ 1} _Oj The second, at Mechanics Calculation benchmark model A ⁱ( i=1,2,3,4 ...) or the not damaged model A of Cable Structure ⁰The basis on, the order by the health status situation of evaluation object done d ^{i+ 1} _o The back upgrade and to obtain next time (the i+1 time, i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation ^I+1At last, through to Mechanics Calculation benchmark model A ^I+1The initial value that calculates monitored amount, by its form next time (promptly the i+1 time, i=1,2,3,4 ...) required " the initial value vector of monitored amount of circulation C ^I+1 _o " ( i=1,2,3,4 ...).

(18)

Mark vector in the formula (18) F ⁱ Subscript iExpression the iInferior circulation, its element F ⁱ _j ( j=1,2,3 ..., subscript N) jExpression the jIndividual by the damage characteristic of evaluation object, can only get 0 and 1 two amount, concrete value rule is seen formula (19).

(19)

Element in the formula (19) F ⁱ _j It is mark vector F ⁱ jIndividual element, D ⁱ _Uj It is nominal unit damage vector D ⁱ _u jIndividual element (seeing formula (3)), d ⁱ _Cj It is current name damage vector d ⁱ _c jIndividual element (seeing formula (12)), they all represent jIndividual by the relevant information of evaluation object.

(20)

In the formula (20) D ⁱ _Uj It is nominal unit damage vector D ⁱ _u jIndividual element (seeing formula (3)), d ⁱ _Cj It is current name damage vector d ⁱ _c jIndividual 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 by evaluation object with bearing generalized displacement component, establishes the quantity of the supporting rope of being assessed and the quantity sum of bearing generalized displacement component and does N, promptly done by the quantity of evaluation object NConfirm 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 uses variable jRepresent this numbering, j=1,2,3 ..., N

The supporting rope with monitored Suo Li of appointment when b. confirming hybrid monitoring is established in the cable system total QThe root rope, the monitored rope force data of structure is by on the structure M ₁ Individual appointment rope M ₁ Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope; Each total M ₁ Individual cable force measurement value or calculated value characterize the rope force information of structure; M ₁ Be one and be not less than 0 integer; The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of structure can be by on the structure K ₂ Individual specified point, and each specified point L ₂ The strain of individual assigned direction is described, and the variation of structural strain data is exactly K ₂ The variation of the tested strain of all of individual specified point; Each total M ₂ Individual strain measurement value or calculated value characterize structural strain, M ₂ For K ₂ With L ₂ Long-pending; M ₂ Be to be not less than 0 integer; The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of structure is by on the structure K ₃ Individual specified point, cross each specified point L ₃ Individual appointment straight line, each specifies straight line H ₃ Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total M ₃ Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M ₃ For K ₃ , L ₃ With H ₃ Long-pending; M ₃ Be one and be not less than 0 integer; When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by on the structure K ₄ Individual specified point, and each specified point L ₄ The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K ₄ The variation of all coordinate components of individual specified point; Each total M ₄ Individual measurement of coordinates value or calculated value characterize planform, M ₄ For K ₄ With L ₄ Long-pending; M ₄ Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total is total MIndividual monitored amount , MFor M ₁ , M ₂ , M ₃ With M ₄ Sum, the definition parameter K, KFor M ₁ , K ₂ , K ₃ With K ₄ Sum, KWith MMust not be less than by the quantity of evaluation object NBecause MIndividual monitored amount is dissimilar, so the present invention is called " based on the progressive method of hybrid monitoring identification damaged cable bearing generalized displacement "; For simplicity, in the present invention that this step is listed " all monitored parameters of structure during hybrid monitoring " abbreviate " monitored amount " as;

C. utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector by the Non-Destructive Testing data of evaluation object etc. d ⁱ _o If during not by the Non-Destructive Testing data of evaluation object, vector d ⁱ _o Each element numerical value get 0; Vector d ⁱ _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 of promptly setting up here d ⁱ _o Can be embodied as d ¹ _o

D. setting up the initial health vector d ¹ _o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount C ⁱ _o Here be circulation for the first time, i gets 1, the initial value vector of the monitored amount of promptly setting up here C ⁱ _o Can be embodied as C ¹ _o Obtain monitored amount initial value vector in actual measurement C ¹ _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 bearing generalized coordinate data of structure; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;

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 bearing generalized coordinate data are set up Cable Structure ⁱ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 ⁱCan be embodied as A ¹

F. at Mechanics Calculation benchmark model A ⁱThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C ⁱ " and " nominal unit damage vector D ⁱ _u ";

G. actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount C ⁱ "; 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 current name damage vector d ⁱ _c With current actual damage vector d ⁱ , the element number of two damage vectors equals by the quantity of evaluation object, current name damage vector d ⁱ _c Element numerical value represent corresponding to the current nominal degree of injury of evaluation object or bearing generalized displacement, current actual damage vector d ⁱ Element numerical value represent corresponding by the current actual damage degree or the bearing generalized 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 of monitored amount C ⁱ " " initial value of monitored amount is vectorial together C ⁱ _o ", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ " and " current name damage vector d ⁱ _c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes in the formula 1 d ⁱ _c Other outer amount is known, finds the solution formula 1 and just can calculate current name damage vector d ⁱ _c

formula 1

J. the current actual damage vector that utilizes formula 2 to express d ⁱ With the initial damage vector d ⁱ _o With current name damage vector d ⁱ _c Element between relation, calculate current actual damage vector d ⁱ All elements;

formula 2

In the formula 2 j=1,2,3 ..., N;

Current actual damage vector d ⁱ Element numerical value represent corresponding to the actual damage degree of evaluation object or actual bearing generalized displacement, according to current actual damage vector d ⁱ Just can define the impaired and degree of injury of which rope, just can confirm actual bearing generalized 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 a generalized displacement component of a bearing, so d ⁱ _j Represent its current generalized displacement numerical value;

K. try to achieve current name damage vector d ⁱ _c After, set up mark vector according to formula 3 F ⁱ , formula 4 has provided mark vector F ⁱ jThe definition of individual element;

formula 3

formula 4

Element in the formula 4 F ⁱ _j It is mark vector F ⁱ jIndividual element, D ⁱ _Uj It is nominal unit damage vector D ⁱ _u jIndividual element, d ⁱ _Cj It is current name damage vector d ⁱ _c jIndividual element, they all represent jIndividual by the relevant information of evaluation object, in the formula 4 j=1,2,3 ..., N;

If mark vector l. F ⁱ Element be 0 entirely, then get back to g step and continue this circulation; If mark vector F ⁱ Element be not 0 entirely, then get into next step, i.e. m step;

M. according to formula 5 calculate next time, promptly iThe initial damage vector that+1 circulation is required d ^{i+ 1} _o Each element d ^{i+ 1} _Oj

formula 5

In the formula 5 D ⁱ _Uj It is nominal unit damage vector D ⁱ _u jIndividual element, d ⁱ _Cj It is current name damage vector d ⁱ _c jIndividual element, F ⁱ _j It is mark vector F ⁱ jIndividual element is in the formula 5 j=1,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 ⁱThe basis on, the order by the health status of evaluation object done d ^{i+ 1} _o The back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time ^I+1

O. pass through Mechanics Calculation benchmark model A ^I+1Calculate corresponding to model A ^I+1Structure all monitored strains point, with the monitored strain numerical value that should change direction, these numerical value are formed next time, the initial value of the required monitored amount that promptly circulates for the i+1 time is vectorial C ^I+1 _o

P. get back to the f step, beginning is circulation next time.

 

In step f, at Mechanics Calculation benchmark model A ⁱThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C ⁱ " and " nominal unit damage vector D ⁱ _u " concrete grammar be:

F1. at the Mechanics Calculation benchmark model A of Cable Structure ⁱThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value NAccording to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or generalized displacement, to be increased unit damage or unit generalized 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 generalized displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases the unit generalized displacement again in this generalized displacement direction; Increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D ⁱ _u " unit damage or the unit generalized 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 generalized 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, has NIndividual just had by evaluation object NIndividual monitored quantitative changeization vector;

F3. by this NIndividual monitored quantitative change vector according to NIndividual by the coding rule of evaluation object, forming successively has NThe monitored amount unit change of the Cable Structure matrix of row Δ C ⁱ

Beneficial effect:Method disclosed by the invention can identify the health status (position and the degree of injury that comprise all damaged cables) of generalized displacement of Cable Structure bearing 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 bearing generalized 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 bearing generalized displacement component assessed are by evaluation object, establishes the quantity of the supporting rope of being assessed and the quantity sum of bearing generalized displacement component and does N, promptly done by the quantity of evaluation object NConfirm 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 uses variable jRepresent this numbering, j=1,2,3 ..., N

The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of structure can be by on the structure K ₂ Individual specified point, and each specified point L ₂ The strain of individual assigned direction is described, and the variation of structural strain data is exactly K ₂ The variation of the tested strain of all of individual specified point.Each total M ₂ Individual strain measurement value or calculated value characterize structural strain, M ₂ For K ₂ With L ₂ Long-pending. M ₂ Be one and be not less than 0 integer.A near point each fixed endpoint that can be exactly each root rope (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge) with the measured point of monitored strain; This specified point can also be a near point the structural bearings; This point generally should not be a stress concentration point; Avoiding occurring excessive strain measurement value, near the fixed endpoint of the rope of the monitored Suo Li of appointment or it when these points generally should all be not hybrid monitoring yet.

The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of structure is by on the structure K ₃ Individual specified point, cross each specified point L ₃ Individual appointment straight line, each specifies straight line H ₃ Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M ₃ Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M ₃ For K ₃ , L ₃ With H ₃ Long-pending. M ₃ Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) or a near point it of each root rope; This specified point can also be a near point the structural bearings, and the point of monitored angle-data generally should all not be chosen as " fixed endpoint of the rope of the monitored Suo Li of appointment or near the point it in the hybrid monitoring " and " point of the monitored strain of appointment or near the point it in the hybrid monitoring "; Can only measure an angle coordinate of specifying straight line at each specified point, for example measure body structure surface normal or the tangent line of the specified point angle coordinate with respect to the acceleration of gravity direction, in fact be exactly measurement of dip angle here.

When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by on the structure K ₄ Individual specified point, and each specified point L ₄ The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K ₄ The variation of all coordinate components of individual specified point.Each total M ₄ Individual measurement of coordinates value or calculated value characterize planform, M ₄ For K ₄ With L ₄ Long-pending. M ₄ Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (for example being the stiff end of drag-line on bridge of cable-stayed bridge) of each root rope, and this specified point can also be a near point the structural bearings, perhaps directly is exactly the structural bearings fulcrum; Here selected monitored point should all not selected " fixed endpoint of the rope of the monitored Suo Li of appointment or near the point it in the hybrid monitoring ", " point of the monitored strain of appointment or near the point it in the hybrid monitoring " and " point of the monitored angle-data of appointment or near the point it in the hybrid monitoring " for use.

Comprehensive above-mentioned monitored amount, total is total with regard to hybrid monitoring MIndividual monitored amount , MFor M ₁ , M ₂ , M ₃ With M ₄ Sum, the definition parameter K, KFor M ₁ , K ₂ , K ₃ With K ₄ Sum, KWith MMust not be less than by the quantity of evaluation object NBecause MIndividual monitored amount is dissimilar, so the present invention is called " based on the progressive method of hybrid monitoring identification damaged cable bearing generalized displacement ".For simplicity, in the present invention that this step is listed " all monitored parameters of structure during hybrid monitoring " abbreviate " monitored amount " as.

Second step: utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector by the Non-Destructive Testing data of evaluation object etc. d ¹ _o If during not by the Non-Destructive Testing data of evaluation object, vector d ¹ _o Each element numerical value get 0; Vector d ¹ _o The coding rule of element with identical by the coding rule of evaluation object.

The 3rd step: at the initial health vector d ¹ _o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount C ¹ _o

The 4th step: the initial value vector that obtains monitored amount in actual measurement C ¹ _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 bearing generalized coordinate data are set up Cable Structure _o, based on Mechanics Calculation benchmark model A _oThe computational data that calculates structure must be very near its measured data, and error generally must not be greater than 5%.

A _oBe constant, only when circulation beginning for the first time, set up; The Mechanics Calculation benchmark model of the Cable Structure of setting up during the i time circulation beginning is designated as A ⁱ, 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 ¹, A among the present invention ¹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 measurement of angle subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate measurement subsystem, signal conditioner etc.), signal (data) collector, the computing machine and the panalarm of communicating by letter.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 ⁱThe basis on carry out the several times Mechanics Calculation, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C ⁱ With nominal unit damage vector D ⁱ _u Concrete grammar is:

A. when the i time circulation beginning, at the Mechanics Calculation benchmark model A of Cable Structure ⁱThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value NAccording 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 generalized displacement by evaluation object again unit damage or unit generalized 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 generalized displacement component of a direction of a bearing, just suppose that this bearing increases the unit generalized displacement again in this generalized displacement direction and (if should be the translational component of the x direction of a bearing by evaluation object for example, just supposes that this bearing has the unit line displacement in the x direction by evaluation object; If should be the angular displacement component around the x axle of a bearing by evaluation object; Just suppose that this bearing has the unit angular displacement around the x axle), increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D ⁱ _u " unit damage or the unit generalized displacement that increase again of all supposition of record record, wherein iExpression the iThe current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure is calculated in inferior circulation each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount and calculated current numerical value vector; 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 generalized 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 ⁱ _j Have NIndividual just had by evaluation object NIndividual monitored quantitative changeization vector δ C ⁱ _j ( J=1,2,3 ..., N).

C. by this NIndividual monitored quantitative change vector according to NIndividual by the coding rule of evaluation object, forming successively has NThe monitored amount unit change of the Cable Structure matrix of row Δ C ⁱ " unit damage monitored quantitative change matrix Δ C ⁱ " the current name damage vector of coding rule and back definition of row d ⁱ _c With current actual damage vector d ⁱ The element coding rule identical.

In this step, reach when giving each vectorial element numbering thereafter; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.

The 9th step: set up the linear relationship error vector e ⁱ And vector g ⁱ Utilize data (" the initial value vector of monitored amount of front C ⁱ _o ", " unit damage monitored quantitative change matrix Δ C ⁱ "); 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 generalized displacement, increase again unit damage or unit generalized displacement are arranged in, calculate each time and form a health status vector d ⁱ _t , the health status vector d ⁱ _t Element number equal by the quantity of evaluation object, vector d ⁱ _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 generalized displacement value of increase of the rope of unit damage, d ⁱ _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 generalized 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; Will C ⁱ _Tj , C ⁱ _o , Δ C ⁱ , d ⁱ _t Bring formula (13) into, formula (13) d ⁱ _c Use d ⁱ _t Bring into, obtain a linear relationship error vector e ⁱ , calculate a linear relationship error vector each time e ⁱ Have NIndividual just had by evaluation object NInferior calculating just has NIndividual linear relationship error vector e ⁱ , with this NIndividual linear relationship error vector e ⁱ Obtain a vector after the addition, with this vector each element divided by NAfter the new vector that obtains be exactly final linear relationship error vector e ⁱ Vector g ⁱ Equal final error vector e ⁱ With vector g ⁱ 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 of monitored amount C ⁱ _o " and " unit damage monitored quantitative change matrix Δ C ⁱ " 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 of monitored amount C ⁱ ".

The 11 step: according to " the current numerical value vector of monitored amount C ⁱ " " initial value of monitored amount is vectorial together C ⁱ _o ", " unit damage monitored quantitative change matrix Δ C ⁱ " and " current name damage vector d ⁱ _c " between the linear approximate relationship (formula (9)) that exists, calculate the current name damage of cable system vector according to multi-objective optimization algorithm d ⁱ _c Noninferior solution.

The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted Sum Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution current name damage vector with the goal programming method only d ⁱ _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) γ ⁱ Be a real number, RBe real number field, area of space Ω has limited vector d ⁱ _c Span (the present embodiment requirements vector of each element d ⁱ _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 γ ⁱ , make formula (22) be met.In the formula (22) G (d ⁱ _c )By formula (23) definition, weighing vector in the formula (22) W ⁱ With γ ⁱ Product representation formula (22) in G (d ⁱ _c )With vector g ⁱ Between the deviation that allows, g ⁱ Definition referring to formula (15), its value will the 8th the step calculate.Vector during actual computation W ⁱ Can with vector g ⁱ Identical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of current name damage vector d ⁱ _c

(21)

(22)

(23)

Try to achieve current name damage vector d ⁱ _c After ,Can be vectorial according to the current actual damage that formula (17) obtain d ⁱ Each element, current actual damage vector d ⁱ 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 bearing generalized displacements more exactly. d ⁱ 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 generalized displacement component of a bearing, so its current generalized displacement numerical value of the numeric representation of this element.

The 12 step: in this circulation, promptly the iTry to achieve current name damage vector in the inferior circulation d ⁱ _c After, set up mark vector according to formula (18), formula (19) F ⁱ If mark vector F ⁱ Element be 0 entirely, then got back to for the tenth step and continue this circulation; If mark vector F ⁱ Element be not 0 entirely, then get into next step, i.e. the 13 step.

The 13 step: according to formula (20) calculate next time, promptly the iThe initial damage vector that+1 circulation is required d ^{i+ 1} _o Each element d ^{i+ 1} _Oj

The 14 step: at Cable Structure Mechanics Calculation benchmark model A ⁱThe basis on, order is the last vector that calculates of step by the health status of evaluation object d ^{i+ 1} _o After, obtain new Mechanics Calculation benchmark model, next time promptly the required Mechanics Calculation benchmark model A of (the i+1 time) circulation ^I+1

The 15 step: through to Mechanics Calculation benchmark model A ^I+1Calculate corresponding to model A ^I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, required vector promptly circulates for the i+1 time C ^I+1 _o , promptly the initial value of monitored amount is vectorial

The 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.

Claims (2)

1. progressive method based on the generalized displacement of hybrid monitoring identification damaged cable bearing is characterized in that said method comprises:

A. for for the purpose of narration is convenient, unifiedly 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 bearing generalized displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with bearing generalized 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; Represent this numbering with variable j, j=1,2,3 ..., N;

The supporting rope with monitored Suo Li of appointment when b. confirming hybrid monitoring is established total Q root rope in the cable system, and the monitored rope force data of structure is by M on the structure ₁The M of individual appointment rope ₁Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope; Each total M ₁Individual cable force measurement value or calculated value characterize the rope force information of structure; M ₁Be one and be not less than 0 integer; The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of structure can be by K on the structure ₂L individual specified point, that reach each specified point ₂The strain of individual assigned direction is described, and the variation of structural strain data is exactly K ₂The variation of the tested strain of all of individual specified point; Each total M ₂Individual strain measurement value or calculated value characterize structural strain, M ₂Be K ₂And L ₂Long-pending; M ₂Be to be not less than 0 integer; The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of structure is by K on the structure ₃L individual specified point, that cross each specified point ₃H individual appointment straight line, each appointment straight line ₃Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total M ₃Individual angle coordinate component measurement value or calculated value characterize the angle information of structure, M ₃Be K ₃, L ₃And H ₃Long-pending; M ₃Be one and be not less than 0 integer; When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by K on the structure ₄L individual specified point, that reach each specified point ₄The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K ₄The variation of all coordinate components of individual specified point; Each total M ₄Individual measurement of coordinates value or calculated value characterize planform, M ₄Be K ₄And L ₄Long-pending; M ₄Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of total, M is M ₁, M ₂, M ₃And M ₄Sum, definition parameter K, K is M ₁, K ₂, K ₃And K ₄Sum, K and M must not be less than by the quantity N of evaluation object; Because M monitored amount is dissimilar, so be called " based on the progressive method of hybrid monitoring identification damaged cable bearing generalized displacement "; For simplicity, this step is listed " all monitored parameters of structure during hybrid monitoring " abbreviate " monitored amount " as;

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. ⁱ _oIf during not by the Non-Destructive Testing data of evaluation object, vectorial d ⁱ _oEach element numerical value get 0; Vector d ⁱ _oThe 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 ⁱ _oCan be embodied as d ¹ _o

D. setting up initial health 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 ⁱ _oHere be circulation for the first time, i gets 1, the initial value vector C of the monitored amount of promptly setting up here ⁱ _oCan be embodied as C ¹ _oObtain monitored amount initial value vector C in actual measurement ¹ _oThe time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure bearing generalized coordinate data of structure; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;

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 bearing generalized coordinate data are set up Cable Structure ⁱ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 ⁱCan be embodied as A ¹

F. at Mechanics Calculation benchmark model A ⁱThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " nominal unit damage vector D ⁱ _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 ⁱ"; 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 ⁱ _cWith current actual damage vector d ⁱ, the element number of two damage vectors equals by the quantity of evaluation object, the vectorial d of current name damage ⁱ _cElement numerical value represent corresponding to the current nominal degree of injury of evaluation object or bearing generalized displacement, current actual damage vector d ⁱElement numerical value represent corresponding by the current actual damage degree or the bearing generalized 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 ⁱ" " the initial value of monitored amount vector C together ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " the vectorial d of current name damage ⁱ _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 the vectorial d of current name damage ⁱ _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 ⁱWith initial damage vector d ⁱ _oWith the vectorial d of current name damage ⁱ _cElement between relation, calculate current actual damage vector d ⁱ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 ⁱElement numerical value represent corresponding to the actual damage degree of evaluation object or actual bearing generalized displacement, according to current actual damage vector d ⁱJust can define the impaired and degree of injury of which rope, just can confirm actual bearing generalized 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 generalized displacement component, a d so of a bearing ⁱ _jRepresent its current generalized displacement numerical value;

K. try to achieve the vectorial d of current name damage ⁱ _cAfter, set up mark vector F according to formula 3 ⁱ, formula 4 has provided mark vector F ⁱ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\&GreaterEqual;D_{\mathrm{Uj}}^i\end{array}\right.$ Formula 4

Element F in the formula 4 ⁱ _jBe mark vector F ⁱJ element, D ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element, d ⁱ _CjBe the vectorial d of current name damage ⁱ _cJ 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. ⁱElement be 0 entirely, then get back to g step and continue this circulation; If mark vector F ⁱ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 _oEach 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 ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element, d ⁱ _CjBe the vectorial d of current name damage ⁱ _cJ element, F ⁱ _jBe mark vector F ⁱJ element, j=1 in the formula 5,2,3 ..., N; Vector d ^I+1 _oThe coding rule of element with identical by the coding rule of evaluation object;

N. at Mechanics Calculation benchmark model A ⁱThe basis on, order is d by the health status of evaluation object ^I+1 _oThe 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+1Calculate corresponding to model A ^I+1Structure 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.

2. the progressive method based on the generalized displacement of hybrid monitoring identification damaged cable bearing according to claim 1 is characterized in that in step f, at Mechanics Calculation benchmark model A ⁱThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " nominal unit damage vector D ⁱ _u" concrete grammar be:

F1. at the Mechanics Calculation benchmark model A of Cable Structure ⁱ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 generalized displacement, to be increased unit damage or unit generalized 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 generalized displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases the unit generalized displacement again in this generalized displacement direction; Increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D ⁱ _u" unit damage or the unit generalized 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 generalized 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 ⁱ