CN102288424A - Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring - Google Patents
Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring Download PDFInfo
- Publication number
- CN102288424A CN102288424A CN2011101226143A CN201110122614A CN102288424A CN 102288424 A CN102288424 A CN 102288424A CN 2011101226143 A CN2011101226143 A CN 2011101226143A CN 201110122614 A CN201110122614 A CN 201110122614A CN 102288424 A CN102288424 A CN 102288424A
- Authority
- CN
- China
- Prior art keywords
- vector
- individual
- monitored
- evaluation object
- damage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention provides a progressive method for recognizing a damaged cable and the angular displacement of a support based on hybrid monitoring. In the method, the defect that the linear relationships between the current value vector of a monitored value and the initial value vector of the monitored value, the change matrix of a unit damage monitored value and a current nominal damage vector are approximate to each other is taken into account. In order to overcome the defect, the invention provides a method of approximating a nonlinear relationship by using linear relationship segments, wherein a large interval is divided into continuous individual small intervals; the linear relationship within each small interval is accurate enough; and the angular displacement of the support and the damaged cable can be recognized quickly through an appropriate algorithm such as a multi-target optimization algorithm and the like within each small interval.
Description
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, as 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, and 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 cable, also comprise the rod member that only bears tensile load) can be impaired, angular displacement also may appear in the bearing of while Cable Structure, and (for example bearing is around coordinate axis X, Y, the rotation of Z, in fact be exactly that bearing is around coordinate axis X, Y, the angular displacement of Z), 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 damaged cable in the cable system of bearing angular displacement and Cable Structure, belong to the engineering structure health monitoring field.
Background technology
The bearing angular 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 angular displacement and Cable Structure is a kind of method that has potentiality.When angular displacement appears in bearing, or the health status of cable system is when changing, when perhaps two kinds of situations take 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 variation (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 of angle coordinate of any imaginary line of the every bit of Cable Structure, the perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily), all these change the health status information that has all comprised cable system, in fact the variation of these measurable parameters has comprised the health status information of cable system, comprised the bearing angular displacement information, therefore can be by 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 this moment, monitored amount was made up of the dissimilar measurable parameter mixing of structure, the present invention claims that this is a hybrid monitoring, that is to say and can utilize hybrid monitoring to discern bearing angular displacement and damaged cable.
In order reliable monitoring and judgement to be arranged to the health status and the bearing angular 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 angular 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 angular 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 angular displacement and damaged cable rationally and effectively.
Technical scheme:If the quantity sum of the quantity of rope and bearing angular displacement component is
NFor sake of convenience, the present invention unitedly calls evaluated rope and bearing angular displacement to be " evaluation object ", gives 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 evaluation object.
The present invention is made up of the two large divisions.Be respectively: one, set up the method for required knowledge base of evaluation object health monitoring systems and parameter, based on knowledge base (containing parameter) with to the evaluation object health status appraisal procedure of the hybrid monitoring data of Cable Structure; Two, the software and hardware part of health monitoring systems.
First of the present invention: foundation is used for the knowledge base of evaluation object health monitoring and the method for parameter.Can be successively circularly as follows, laddering carrying out:
The first step: during circulation beginning each time, the evaluation object initial health vector when at first needing to set up or set up this circulation beginning
d o i (
i=1,2,3 ...), set up the initial mechanical calculating benchmark model A of Cable Structure
o(for example finite element benchmark model, A in the present invention
oBe constant), set up the Mechanics Calculation benchmark model A of Cable Structure
i(finite element benchmark model for example,
i=1,2,3 ...).Letter i is except the place of representing number of steps significantly, and is alphabetical in the present invention
iOnly represent cycle index, promptly
iInferior circulation.
The
iCable Structure " the initial health vector that needs when inferior circulation begins
d o i " (as the formula (1)), use
d o i Expression the
iCable Structure when inferior circulation begins is (with Mechanics Calculation benchmark model A
iThe initial health of Cable Structure expression).
In the formula (1)
d i Oj (
i=1,2,3,
; j=1,2,3 ....,
N) expression circulates for the i time when beginning, Mechanics Calculation benchmark model A
iIn cable system
jThe current health status of individual evaluation object is if this evaluation object is the rope (or pull bar) in the cable system, so
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, represents to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if this evaluation object is an angular displacement component of a bearing, so
d i Represent its current angular 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 1 o ) time, Non-Destructive Testing data of utilizing rope etc. can be expressed the data and the bearing angle displacement measurement of the health status of rope and set up evaluation object initial health vector
d 1 o If when not having the data of the Non-Destructive Testing data of rope and other health status that can express rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vector
d 1 o In get 0 with each element numerical value of Suo Xiangguan.
The i time (
i=2,3,4,5,6 ...) the evaluation object initial health vector of needs when circulation begins
d i o , be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is described below.
The Mechanics Calculation benchmark model that need set up during the i time circulation beginning or the Mechanics Calculation benchmark model of having set up are designated as A
i
Measured data according to the Cable Structure in the Cable Structure completion (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure bearing angular data, Cable Structure modal data, to cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope) and design drawing, as-constructed drawing, utilize mechanics method (for example finite element method) to set up A
oIf there is not the measured data of the structure in the Cable Structure completion, so just before setting up health monitoring systems, structure is surveyed, the measured data that obtains Cable Structure (comprises the Cable Structure shape data, the rope force data, the draw-bar pull data, Cable Structure support coordinate data, Cable Structure bearing angular 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.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
1, A
1Just equal A
oA
1The health status of corresponding evaluation object by
d 1 o Describe.
The i time (
i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins
i, be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is described below.
Existing Mechanics Calculation benchmark model A
1With evaluation object initial health vector
d 1 o After, model A
1In the health status of each evaluation object by vector
d 1 o Express.At A
1The basis on, the health status numerical value of all evaluation objects is changed to 0, mechanical model A
1The health status that is updated to all evaluation objects all is that 0 mechanical model (is designated as A
0), mechanical model A
0Be actually the mechanical model of the Cable Structure correspondence of excellent no bearing angular displacement.Might as well claim model A
0For the not damaged of Cable Structure does not have bearing angular displacement model A
0
Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate are described below 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 1 Individual appointment rope
M 1 Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope.Each total
M 1 Individual cable force measurement value or calculated value characterize the rope force information of structure.
M 1 Be one and be not less than 0 integer.
The monitored strain data of structure can be by on the structure
K 2 Individual specified point, and each specified point
L 2 The strain of individual assigned direction is described, and the variation of structural strain data is exactly
K 2 The variation of the tested strain of all of individual specified point.Each total
M 2 (M 2 =K 2 * L 2 )Individual strain measurement value or calculated value characterize structural strain.
M 2 Be one and be not less than 0 integer.
The monitored angle-data of structure is by on the structure
K 3 Individual specified point, cross each specified point
L 3 Individual appointment straight line, each specifies straight line
H 3 Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total
M 3 (M 3 =K 3 * L 3 * H 3 )Individual angle coordinate component measurement value or calculated value characterize the angle information of structure.
M 3 Be one and be not less than 0 integer.
The monitored shape data of structure is by on the structure
K 4 Individual specified point, and each specified point
L 4 The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly
K 4 The variation of all coordinate components of individual specified point.Each total
M 4 (M 4 =K 4 * L 4 )Individual measurement of coordinates value or calculated value characterize planform.
M 4 Be one and be not less than 0 integer.
Comprehensive above-mentioned monitored amount, total is total
M(M=M 1 + M 2 + M 3 + M 4 )Individual monitored amount, the definition parameter
K(
K=M 1 + K 2 + K 3 + K 4 ),
KWith
MMust not be less than the quantity of evaluation object
NBecause
MIndividual monitored amount is dissimilar, so the present invention is called " based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular displacement ".
For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " monitored amount ".
Among the present invention with monitored amount initial value vector
C i 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 i o Full name be the initial value vector of monitored amount " the i time circulation ".
In the formula (2)
C i 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 i o Be by previously defined
MIndividual monitored amount forms according to certain series arrangement, and this is put in order there is no specific (special) requirements, only require all associated vector of back also in this order array data get final product.
During circulation beginning for the first time, " the initial value vector of the 1st the monitored amount that circulates
C 1 o " (seeing formula (2)) be made up of measured data, because according to model A
1The 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 ...) " the initial value vector of the i time monitored amount of circulation of needs when circulation begins
C i o ", be preceding once (promptly the i-1 time,
i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is described below.
Second step: circulation each time needs to set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector ", and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as
Δ C i , " nominal unit damage vector " that the i time circulation set up is designated as
D i u ,
i=1,2,3 ...
The Cable Structure " unit damage monitored numerical quantity transformation matrices " that circulation is for the first time set up is designated as
Δ C 1 Set up
Δ C 1 Process as follows:
Mechanics Calculation benchmark model A in Cable Structure
1The basis on carry out several times and calculate, equal on the calculation times numerical value
NCalculating hypothesis each time has only an evaluation object that unit damage is arranged, concrete, if this evaluation object is a support cable in the cable system, so just suppose that this support cable has unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the angular displacement component of a direction of a bearing, just suppose that this bearing is this angular displacement direction generation unit angular displacement (for example get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. and be the unit angular displacement).For sake of convenience, the present invention is referred to as unit damage with the damage and the bearing angular displacement of the support cable of supposition.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, health status numerical value that calculate, evaluation object---be called nominal health status vector
d i c (
i=1,2,3 ...), all with respect to this time when beginning circulation, with the health status of Cable Structure as being healthy fully speech, the nominal health status numerical value that the formula that therefore must foundation hereinafter provides will calculate be converted into true health status numerical value).The evaluation object that occurs unit damage during a round-robin calculates each time together is different from the evaluation object that occurs unit damage in other time calculating, and supposition each time has the unit damage value of the evaluation object of unit damage can be different from the unit damage value of other evaluation objects, uses " nominal unit damage vector
D i u " (as the formula (3)) write down the unit damage of the supposition of all evaluation objects in each time circulation, circulation time is designated as for the first time
D 1 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
jWhen individual evaluation object had unit damage, available formula (4) was represented all appointments
MThe current numerical value vector of the calculating of individual monitored amount
C 1 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 1 o , the gained vector is exactly that " the numerical value change vector of monitored amount " of (with the mark that is numbered of evaluation object that unit damage is arranged) is (when under this condition
jWhen individual evaluation object has unit damage, use
δ C 1 j The numerical value change vector of representing monitored amount,
δ C 1 j Definition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vector for formula (4)
D 1 u jIndividual element
D Uj Gained), the numerical value change vector of monitored amount
δ C 1 j Each element representation since when calculating supposition that evaluation object (for example the of unit damage is arranged
jIndividual 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 evaluation object just has
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 1 "
(MOK
NRow
), each vector
δ C 1 j (
j=1,2,3 ....,
N) be matrix
Δ C 1 One row,
Δ C 1 Definition as the formula (8).
Nominal unit damage vector in the formula (3)
D i u Element
D i Uj (
i=1,2,3,
j=1,2,3 ....,
N) expression the
iSuppose in the inferior circulation
jThe unit damage numerical value of individual evaluation object, vector
D i u In the numerical value of each element can be the same or different.
Element in the formula (4)
C i Tjk (
i=1,2,3,
j=1,2,3 ....,
N; k=1,2,3 ....,
M; M 〉=N) expression the
iInferior circulation is because the
jWhen individual evaluation object has unit damage, according to coding rule pairing
kThe current numerical value of the calculating of the monitored amount of individual appointment.
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 evaluation object has unit damage, in the formula
D i Uj It is vector
D i u In
jIndividual element.Vector
δ C i j Definition as the formula (6),
δ C i j k(
k=1,2,3 ....,
M; M 〉=N) individual element
δ C i Jk Expression the
iIn the inferior circulation, set up matrix
Δ C i The time, suppose
jWhen having unit damage, individual evaluation object calculates gained the
kThe change amount of individual monitored amount is with respect to the unit damage of supposition
D i Uj Rate of change, it defines as the formula (7).
The definition of each amount has been previously described in the formula (7).
Vector in the formula (8)
δ C i j (
i=1,2,3 ....,,
j=1,2,3 ....,
N) expression the
iIn the inferior circulation, because the
jIndividual evaluation object has unit damage
D i Uj Cause, the relative value of all monitored amounts changes.Matrix
Δ C i Row (subscript
j) coding rule and front vector
d i o The subscript of element
jCoding rule identical.
The 3rd step: the current health status of identification evaluation object.Detailed process is as follows.
The
i(
i=1,2,3 ...) in the inferior circulation, " current (calculating or actual measurement) numerical value vector of monitored amount
C i " " the initial value vector of monitored amount together
C i o ", " unit damage monitored numerical quantity transformation matrices
Δ C i " and " current name damage vector
d i c " between linear approximate relationship, shown in (9) or formula (10).
Current (calculating or actual measurement) numerical value vector of monitored amount in formula (9) and the formula (10)
C i Definition be similar to the initial value vector of monitored amount
C i o Definition, see formula (11); The current name damage of evaluation object vector
d i c Definition see formula (12).
Element in the formula (11)
C i 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.
In the formula (12)
d i Cj (
i=1,2,3 ....;
j=1,2,3 ....,
N) be
iCable Structure in the inferior circulation
jThe current nominal impairment value of individual evaluation object, vector
d i c The subscript of element
jCoding rule and matrix
Δ C i The coding rule of row identical.
When evaluation object actual damage or bearing angular displacement are not too big, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (9) or formula (10) is less with the error of actual conditions, and error can be used error vector
e i (formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).
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 i " directly find the solution and obtain current name damage vector
d i c And it is vectorial to obtain current name damage
d i c Acceptable separating (promptly have reasonable error, but can be more accurately from cable system, determine damaged cable position and degree of injury thereof, determine the bearing angular displacement) 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 i Description departs from ideal linearity relation (formula (9) or formula (10))
Reasonable deviation, define by formula (15).
In the formula (15)
g i 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 i Can be according to the error vector of formula (13) definition
e i Tentative calculation is selected.
Initial value vector in monitored amount
C i o (survey or calculate), unit damage monitored numerical quantity transformation matrices
Δ C i The current numerical value vector of (calculating) and monitored amount
C i 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 i c Acceptable separating, current actual damage vector
d i The element of (formula (16) is seen in definition) can calculate according to formula (17), has just obtained the current actual damage vector of evaluation object
d i Thereby, can by
d i Determine damaged cable position and degree of injury, determine the bearing angular displacement, realized that just damage identification and bearing angular displacement discern.
In the formula (16)
d i j (
i=1,2,3,
; j=1,2,3 ....,
N) expression the
iIn the inferior circulation
jThe actual damage value of individual evaluation object, its definition are seen formula (17), if this evaluation object is the rope (or pull bar) in the cable system, so
d i j Represent its current damage,
d i j Being to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, represents the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%, if this evaluation object is an angular displacement component of a bearing, so
d i j Represent its current angular displacement numerical value.Vector
d i The coding rule of element and formula (1) in vector
d i o The coding rule of element identical.
In the formula (17)
d i Oj (
i=1,2,3,4,
; j=1,2,3 ....,
N) be vector
d i o jIndividual element,
d i Cj It is vector
d i c jIndividual element.
The 4th step: judge whether to finish this (the
iInferior) circulation, if then finish 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 as follows.
This (
iInferior) try to achieve current name damage vector in the circulation
d i c After, at first, set up mark vector according to formula (18)
F i , formula (19) has provided mark vector
F i jThe definition of individual element; If mark vector
F i Element be 0 entirely, then in this circulation, continue health monitoring and calculating to Cable Structure; If mark vector
F i Element be not 0 entirely, then finish subsequent step after, enter next time circulation.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(
i=1,2,3,4 ...) or the not damaged model A of Cable Structure
0The basis on, make the health status situation of evaluation object be
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, by 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 ...).
Mark vector in the formula (18)
F i Subscript
iExpression the
iInferior circulation, its element
F i j (
j=1,2,3 ..., subscript N)
jExpression the
jThe damage characteristic of individual evaluation object can only be got 0 and 1 two amount, and concrete value rule is seen formula (19).
Element in the formula (19)
F i j It is mark vector
F i jIndividual element,
D i Uj It is nominal unit damage vector
D i u jIndividual element (seeing formula (3)),
d i Cj It is current name damage vector
d i c jIndividual element (seeing formula (12)), they all represent
jThe relevant information of individual evaluation object.
(20)
In the formula (20)
D i Uj It is nominal unit damage vector
D i u jIndividual element (seeing formula (3)),
d i Cj It is current name damage vector
d i 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 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. for sake of convenience, it is evaluation object that the present invention unitedly calls evaluated support cable and bearing angular displacement component, establishes the quantity of evaluated support cable and the quantity sum of bearing angular displacement component to be
N, promptly the quantity of evaluation object is
NDetermine the coding rule of evaluation object, with evaluation object numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention's variable
jRepresent this numbering,
j=1,2,3 ...,
N
The support cable with monitored Suo Li of appointment when b. determining 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 1 Individual appointment rope
M 1 Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope; Each total
M 1 Individual cable force measurement value or calculated value characterize the rope force information of structure;
M 1 Be one and be not less than 0 integer; The measured point with monitored strain of appointment when determining hybrid monitoring, the monitored strain data of structure can be by on the structure
K 2 Individual specified point, and each specified point
L 2 The strain of individual assigned direction is described, and the variation of structural strain data is exactly
K 2 The variation of the tested strain of all of individual specified point; Each total
M 2 Individual strain measurement value or calculated value characterize structural strain,
M 2 For
K 2 With
L 2 Long-pending;
M 2 Be to be not less than 0 integer; The measured point with monitored angle of appointment when determining hybrid monitoring, the monitored angle-data of structure is by on the structure
K 3 Individual specified point, cross each specified point
L 3 Individual appointment straight line, each specifies straight line
H 3 Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total
M 3 Individual angle coordinate component measurement value or calculated value characterize the angle information of structure,
M 3 For
K 3 , L 3 With
H 3 Long-pending;
M 3 Be one and be not less than 0 integer; When determining hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by on the structure
K 4 Individual specified point, and each specified point
L 4 The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly
K 4 The variation of all coordinate components of individual specified point; Each total
M 4 Individual measurement of coordinates value or calculated value characterize planform,
M 4 For
K 4 With
L 4 Long-pending;
M 4 Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total is total
MIndividual monitored amount
, MFor
M 1 , M 2 , M 3 With
M 4 Sum, the definition parameter
K,
KFor
M 1 , K 2 , K 3 With
K 4 Sum,
KWith
MMust not be less than the quantity of evaluation object
NBecause
MIndividual monitored amount is dissimilar, so the present invention is called " based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular 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. Non-Destructive Testing data of utilizing evaluation object etc. can be expressed the data of the health status of evaluation object and set up evaluation object initial health vector
d i o If when not having the Non-Destructive Testing data of evaluation object, vector
d i o Each element numerical value get 0; Vector
d i o The coding rule of element and the coding rule of evaluation object identical; The present invention represents cycle index with i, i=1, and 2,3, Here be circulation for the first time, i gets 1, i.e. the initial health vector of setting up here
d i o Can be embodied as
d 1 o
D. setting up the initial health vector
d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount
C i o Here be circulation for the first time, i gets 1, i.e. the initial value vector of the monitored amount of setting up here
C i o Can be embodied as
C 1 o Obtain monitored amount initial value vector in actual measurement
C 1 o The time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure bearing angular data of structure;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure bearing angular data is set up Cable Structure
iHere be circulation for the first time, i gets 1, i.e. the Mechanics Calculation benchmark model A of the Cable Structure of setting up here
iCan be embodied as A
1
F. at Mechanics Calculation benchmark model A
iThe basis on carry out the several times Mechanics Calculation, by calculate obtaining " unit damage monitored numerical quantity transformation matrices
Δ C i " and " nominal unit damage vector
D i u ";
G. actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount
C i "; When numbering to the element of the institute's directed quantity that occurred before this step and this step, should use same coding rule, can guarantee the element each vector, that numbering is identical that occurs before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
H. define current name damage vector
d i c With current actual damage vector
d i , the element number of two damage vectors equals the quantity of evaluation object, current name damage vector
d i c Element numerical value represent the current nominal degree of injury or the bearing angular displacement of corresponding evaluation object, current actual damage vector
d i Element numerical value represent the current actual damage degree or the bearing angular displacement of corresponding evaluation object, the element number of the element of two damage vectors equals the quantity of evaluation object, between the element of two damage vectors and the evaluation object is one-to-one relationship, and the coding rule of the element of two damage vectors is identical with the coding rule of evaluation object;
I. according to " the current numerical value vector of monitored amount
C i " " the initial value vector of monitored amount together
C i o ", " unit damage monitored numerical quantity transformation matrices
Δ C i " and " current name damage vector
d i c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes in the formula 1
d i c Other outer amount is known, finds the solution formula 1 and just can calculate current name damage vector
d i c
J. the current actual damage vector that utilizes formula 2 to express
d i With the initial damage vector
d i o With current name damage vector
d i c Element between relation, calculate current actual damage vector
d i All elements;
In the formula 2
j=1,2,3 ..., N;
Current actual damage vector
d i Element numerical value represent the actual damage degree or the actual bearing angular displacement of corresponding evaluation object, according to current actual damage vector
d i Just can define the impaired and degree of injury of which rope, just can determine actual bearing angular 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, represent that the pairing rope of this element is intact, do not damage, if its numerical value is 100%, represent that then the pairing rope of this element has completely lost load-bearing capacity, 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 an angular displacement component of a bearing, so
d i j Represent its current angular displacement numerical value;
K. try to achieve current name damage vector
d i c After, set up mark vector according to formula 3
F i , formula 4 has provided mark vector
F i jThe definition of individual element;
Element in the formula 4
F i j It is mark vector
F i jIndividual element,
D i Uj It is nominal unit damage vector
D i u jIndividual element,
d i Cj It is current name damage vector
d i c jIndividual element, they all represent
jThe relevant information of individual evaluation object is in the formula 4
j=1,2,3 ..., N;
If mark vector l.
F i Element be 0 entirely, then get back to g step and continue this circulation; If mark vector
F i Element be not 0 entirely, then enter next step, i.e. m step;
M. according to formula 5 calculate next time, promptly
i+ 1 required initial damage vector of circulation
d i+
1 o Each element
d i+
1 Oj
In the formula 5
D i Uj It is nominal unit damage vector
D i u jIndividual element,
d i Cj It is current name damage vector
d i c jIndividual element,
F i j It is mark vector
F i jIndividual element is in the formula 5
j=1,2,3 ..., N; Vector
d i+
1 o The coding rule of element and the coding rule of evaluation object identical;
N. at Mechanics Calculation benchmark model A
iThe basis on, make the health status of evaluation object be
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 vector of required monitored amount promptly circulates for the i+1 time
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
iThe basis on carry out the several times Mechanics Calculation, by calculate obtaining " unit damage monitored numerical quantity transformation matrices
Δ C i " and " nominal unit damage vector
D i u " concrete grammar be:
F1. at the Mechanics Calculation benchmark model A of Cable Structure
iThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value
NCoding rule according to evaluation object calculates successively; Calculating hypothesis each time has only an evaluation object to increase unit damage or unit angular displacement again on the basis of original damage or angular displacement, concrete, if this evaluation object is a support cable in the cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is the angular displacement component of a direction of a bearing, just suppose that this bearing increases the unit angular displacement again in this angular displacement direction, the evaluation object that increases unit damage or unit angular displacement in calculating each time again is different from the evaluation object that increases unit damage or unit angular displacement in other time calculating again, with " nominal unit damage vector
D i u " unit damage or the unit angular displacement that increase again of all supposition of record record; wherein i represents 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 monitored amount of current calculated value composition of the monitored amount of all that calculate is calculated current numerical value vector each time;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit angular 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 change vector, has
NIndividual evaluation object just has
NIndividual monitored quantitative change vector;
F3. by this
NIndividual monitored quantitative change vector according to
NThe coding rule of individual evaluation object, forming successively has
NThe monitored amount unit change of the Cable Structure matrix of row
Δ C i
Beneficial effect:Method disclosed by the invention can identify the health status (position and the degree of injury that comprise all damaged cables) of angular 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
At 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 angular displacement component in the cable system rationally and effectively simultaneously.The following describes of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
The present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status of the cable system of Cable Structure.During concrete enforcement, the following step is a kind of in the various steps that can take.
The first step: for sake of convenience, it is evaluation object that the present invention unitedly calls evaluated support cable and bearing angular displacement component, establishes the quantity of evaluated support cable and the quantity sum of bearing angular displacement component to be
N, promptly the quantity of evaluation object is
NDetermine the coding rule of evaluation object, with evaluation object numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention's variable
jRepresent this numbering,
j=1,2,3 ...,
N
The measured point with monitored strain of appointment when determining hybrid monitoring, the monitored strain data of structure can be by on the structure
K 2 Individual specified point, and each specified point
L 2 The strain of individual assigned direction is described, and the variation of structural strain data is exactly
K 2 The variation of the tested strain of all of individual specified point.Each total
M 2 Individual strain measurement value or calculated value characterize structural strain,
M 2 For
K 2 With
L 2 Long-pending.
M 2 Be one and be not less than 0 integer.A near point each fixed endpoint that can be exactly each root rope (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge) with the measured point of monitored strain, this 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, the fixed endpoint of the rope of the monitored Suo Li of appointment or in its vicinity when these points generally should all be not hybrid monitoring yet.
The measured point with monitored angle of appointment when determining hybrid monitoring, the monitored angle-data of structure is by on the structure
K 3 Individual specified point, cross each specified point
L 3 Individual appointment straight line, each specifies straight line
H 3 Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total
M 3 Individual angle coordinate component measurement value or calculated value characterize the angle information of structure,
M 3 For
K 3 , L 3 With
H 3 Long-pending.
M 3 Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) or a near point it of each root rope, 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 point in its vicinity in the hybrid monitoring " and " point of the monitored strain of appointment or point in its vicinity in the hybrid monitoring "; Can only measure an angle coordinate of specifying straight line at each specified point, for example measure the body structure surface normal of specified point or the tangent line angle coordinate with respect to the acceleration of gravity direction, in fact be exactly measurement of dip angle here.
When determining hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by on the structure
K 4 Individual specified point, and each specified point
L 4 The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly
K 4 The variation of all coordinate components of individual specified point.Each total
M 4 Individual measurement of coordinates value or calculated value characterize planform,
M 4 For
K 4 With
L 4 Long-pending.
M 4 Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (for example being the stiff end of drag-line on bridge of cable-stayed bridge) of each root rope, and this specified point can also be a near point the structural bearings, perhaps directly is exactly the structural bearings fulcrum; Here Xuan Ding monitored point should all not selected " fixed endpoint of the rope of the monitored Suo Li of appointment or point in its vicinity in the hybrid monitoring ", " point of the monitored strain of appointment or point in its vicinity in the hybrid monitoring " and " point of the monitored angle-data of appointment or point in its vicinity in the hybrid monitoring " for use.
Comprehensive above-mentioned monitored amount, total is total with regard to hybrid monitoring
MIndividual monitored amount
, MFor
M 1 , M 2 , M 3 With
M 4 Sum, the definition parameter
K,
KFor
M 1 , K 2 , K 3 With
K 4 Sum,
KWith
MMust not be less than the quantity of evaluation object
NBecause
MIndividual monitored amount is dissimilar, so the present invention is called " based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular 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: the data of utilizing the Non-Destructive Testing data etc. of evaluation object can express the health status of evaluation object are set up evaluation object initial health vector
d 1 o If when not having the Non-Destructive Testing data of evaluation object, vector
d 1 o Each element numerical value get 0; Vector
d 1 o The coding rule of element and the coding rule of evaluation object identical.
The 3rd step: at the initial health vector
d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount
C 1 o
The 4th step: the initial value vector that obtains monitored amount in actual measurement
C 1 o The time, can adopt ripe measuring method to carry out cable force measurement, strain measurement, measurement of angle and volume coordinate and measure.Simultaneously, calculate the initial Suo Li and the Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) of all ropes of Cable Structure after directly measuring or measuring, the original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is just can determine according to these coordinate datas the geometric properties of Cable Structure.For cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
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 angular data is 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
i, wherein i represents cycle index; Alphabetical i is except the place of representing number of steps significantly in the application form of the present invention, and alphabetical i only represents cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A
1, A among the present invention
1Just 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 the health status that monitors evaluation object changes, 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 finished functions such as the needed monitoring of each task of the present invention, record, control, storage, calculating, notice, warning (being 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 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 finish necessary evaluation work automatically.
The 7th step: the step begins circulation running thus, is designated as the i time circulation for sake of convenience, i=1 wherein, and 2,3,4,5 ...
The 8th step: the Mechanics Calculation benchmark model in Cable Structure is designated as A
iThe basis on carry out the several times Mechanics Calculation, by calculate obtaining Cable Structure unit damage monitored quantitative change matrix
Δ C i With nominal unit damage vector
D i u Concrete grammar is:
A. when the i time circulation beginning, at the Mechanics Calculation benchmark model A of Cable Structure
iThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value
NCoding rule according to evaluation object calculates successively; Calculating hypothesis each time has only an evaluation object to increase on the basis of original damage or angular displacement again unit damage or unit angular displacement is arranged, concrete, if this evaluation object is a support cable in the cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is the angular displacement component of a direction of a bearing, just suppose that this bearing increases the unit angular displacement again in this angular displacement direction, the evaluation object that increases unit damage or unit angular displacement in calculating each time again is different from the evaluation object that increases unit damage or unit angular displacement in other time calculating again, with " nominal unit damage vector
D i u " unit damage or the unit angular 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 monitored amount of current calculated value composition of the monitored amount of all that calculate is calculated current numerical value vector each time; 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 evaluation object object.
B. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit angular 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 change vector
δ C i j Have
NIndividual evaluation object just has
NIndividual monitored quantitative change vector
δ C i j (
J=1,2,3 ..., N).
C. by this
NIndividual monitored quantitative change vector according to
NThe coding rule of individual evaluation object, forming successively has
NThe monitored amount unit change of the Cable Structure matrix of row
Δ C i " unit damage monitored quantitative change matrix
Δ C i " the coding rule of row and the current name damage vector that defines later
d i c With current actual damage vector
d i The element coding rule identical.
Reach in this step when giving each vectorial element numbering thereafter, should use same coding rule with other vector among the present invention, can guarantee any one element in each vector in this step like this, with element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.
The 9th step: set up the linear relationship error vector
e i And vector
g i Utilize data (" the initial value vector of monitored amount of front
C i o ", " unit damage monitored quantitative change matrix
Δ C i "); when the 8th step calculated each time; when supposing in calculating each time that promptly having only an evaluation object to increase in the cable system has unit damage or unit angular displacement on the basis of original damage or angular displacement again, calculate each time and form a health status vector
d i t , the health status vector
d i t Element number equal the quantity of evaluation object, vector
d i t All elements in have only the numerical value of an element to get to calculate each time in hypothesis increase the unit damage value of rope of unit damage or the unit value of angular displacement of increase,
d i t The numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the evaluation object of unit damage or unit angular displacement, be identical with the element of the same numbering of other vectors with the corresponding relation of this rope; Will
C i Tj , C i o ,
Δ C i ,
d i t Bring formula (13) into, formula (13)
d i c With
d i t Bring into, obtain a linear relationship error vector
e i , calculate a linear relationship error vector each time
e i Have
NIndividual evaluation object just has
NInferior calculating just has
NIndividual linear relationship error vector
e i , with this
NIndividual linear relationship error vector
e i Obtain a vector after the addition, with each element of this vector divided by
NAfter the new vector that obtains be exactly final linear relationship error vector
e i Vector
g i Equal final error vector
e i With vector
g i Be kept on the hard disc of computer of operation health monitoring systems software, use for health monitoring systems software.With " the initial value vector of monitored amount
C i o " and " unit damage monitored quantitative change matrix
Δ C i " etc. parameter be kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
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 i ".
The 11 step: according to " the current numerical value vector of monitored amount
C i " " the initial value vector of monitored amount together
C i o ", " unit damage monitored quantitative change matrix
Δ C i " and " current name damage vector
d i 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 i c Noninferior solution.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted Sum Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution current name damage vector with the goal programming method only
d i c Process, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.
According to the goal programming method, formula (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (21) and the formula (22), in the formula (21)
γ i Be a real number,
RBe real number field, area of space Ω has limited vector
d i c Span (the present embodiment requirements vector of each element
d i c Each element be not less than 0, be not more than 1).The meaning of formula (21) is to seek the real number of an absolute value minimum
γ i , make formula (22) be met.In the formula (22)
G (d i c )By formula (23) definition, weighing vector in the formula (22)
W i With
γ i Product representation formula (22) in
G (d i c )With vector
g i Between the deviation that allows,
g i Definition referring to formula (15), its value will the 8th the step calculate.Vector during actual computation
W i Can with vector
g i Identical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of current name damage vector
d i c
(21)
(22)
(23)
Try to achieve current name damage vector
d i c After
,The current actual damage vector that can obtain according to formula (17)
d i Each element, current actual damage vector
d i Have reasonable error exactly but can determine the position of damaged cable more exactly and the separating of degree of injury, can determine separating of all bearing angular displacements more exactly.
d i Each element corresponding to the health status of an evaluation object, if this evaluation object is the rope (or pull bar) in the cable system, 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, be to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, the load-bearing capacity of representing this rope forfeiture corresponding proportion in the time of between 0 and 100%, if this evaluation object is an angular displacement component of a bearing, its current angular displacement numerical value of the numeric representation of this element so.
The 12 step: in this circulation, promptly the
iTry to achieve current name damage vector in the inferior circulation
d i c After, set up mark vector according to formula (18), formula (19)
F i If mark vector
F i Element be 0 entirely, then got back to for the tenth step and continue this circulation; If mark vector
F i Element be not 0 entirely, then enter next step, i.e. the 13 step.
The 13 step: according to formula (20) calculate next time, promptly the
i+ 1 required initial damage vector of circulation
d i+
1 o Each element
d i+
1 Oj
The 14 step: at Cable Structure Mechanics Calculation benchmark model A
iThe basis on, the health status that makes evaluation object is the vector that previous step calculates
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: by 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 , i.e. the initial value vector of monitored amount
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. one kind based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular displacement, it is characterized in that described method comprises:
A. for sake of convenience, unitedly calling evaluated support cable and bearing angular displacement component is evaluation object, establishes the quantity of evaluated support cable and the quantity sum of bearing angular displacement component to be
N, promptly the quantity of evaluation object is
NDetermine the coding rule of evaluation object, with evaluation object numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; Use variable
jRepresent this numbering,
j=1,2,3 ...,
N
The support cable with monitored Suo Li of appointment when b. determining 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 1 Individual appointment rope
M 1 Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope; Each total
M 1 Individual cable force measurement value or calculated value characterize the rope force information of structure;
M 1 Be one and be not less than 0 integer; The measured point with monitored strain of appointment when determining hybrid monitoring, the monitored strain data of structure can be by on the structure
K 2 Individual specified point, and each specified point
L 2 The strain of individual assigned direction is described, and the variation of structural strain data is exactly
K 2 The variation of the tested strain of all of individual specified point; Each total
M 2 Individual strain measurement value or calculated value characterize structural strain,
M 2 For
K 2 With
L 2 Long-pending;
M 2 Be to be not less than 0 integer; The measured point with monitored angle of appointment when determining hybrid monitoring, the monitored angle-data of structure is by on the structure
K 3 Individual specified point, cross each specified point
L 3 Individual appointment straight line, each specifies straight line
H 3 Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total
M 3 Individual angle coordinate component measurement value or calculated value characterize the angle information of structure,
M 3 For
K 3 , L 3 With
H 3 Long-pending;
M 3 Be one and be not less than 0 integer; When determining hybrid monitoring appointment with monitored shape data, the monitored shape data of structure is by on the structure
K 4 Individual specified point, and each specified point
L 4 The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly
K 4 The variation of all coordinate components of individual specified point; Each total
M 4 Individual measurement of coordinates value or calculated value characterize planform,
M 4 For
K 4 With
L 4 Long-pending;
M 4 Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total is total
MIndividual monitored amount
, MFor
M 1 , M 2 , M 3 With
M 4 Sum, the definition parameter
K,
KFor
M 1 , K 2 , K 3 With
K 4 Sum,
KWith
MMust not be less than the quantity of evaluation object
NBecause
MIndividual monitored amount is dissimilar, so the present invention is called " based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular 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. Non-Destructive Testing data of utilizing evaluation object etc. can be expressed the data of the health status of evaluation object and set up evaluation object initial health vector
d i o If when not having the Non-Destructive Testing data of evaluation object, vector
d i o Each element numerical value get 0; Vector
d i o The coding rule of element and the coding rule of evaluation object identical; The present invention represents cycle index with i, i=1, and 2,3, Here be circulation for the first time, i gets 1, i.e. the initial health vector of setting up here
d i o Can be embodied as
d 1 o
D. setting up the initial health vector
d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount
C i o Here be circulation for the first time, i gets 1, i.e. the initial value vector of the monitored amount of setting up here
C i o Can be embodied as
C 1 o Obtain monitored amount initial value vector in actual measurement
C 1 o The time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure bearing angular data of structure;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure bearing angular data is set up Cable Structure
iHere be circulation for the first time, i gets 1, i.e. the Mechanics Calculation benchmark model A of the Cable Structure of setting up here
iCan be embodied as A
1
F. at Mechanics Calculation benchmark model A
iThe basis on carry out the several times Mechanics Calculation, by calculate obtaining " unit damage monitored numerical quantity transformation matrices
Δ C i " and " nominal unit damage vector
D i u ";
G. actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount
C i "; When numbering to the element of the institute's directed quantity that occurred before this step and this step, should use same coding rule, can guarantee the element each vector, that numbering is identical that occurs before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
H. define current name damage vector
d i c With current actual damage vector
d i , the element number of two damage vectors equals the quantity of evaluation object, current name damage vector
d i c Element numerical value represent the current nominal degree of injury or the bearing angular displacement of corresponding evaluation object, current actual damage vector
d i Element numerical value represent the current actual damage degree or the bearing angular displacement of corresponding evaluation object, the element number of the element of two damage vectors equals the quantity of evaluation object, between the element of two damage vectors and the evaluation object is one-to-one relationship, and the coding rule of the element of two damage vectors is identical with the coding rule of evaluation object;
I. according to " the current numerical value vector of monitored amount
C i " " the initial value vector of monitored amount together
C i o ", " unit damage monitored numerical quantity transformation matrices
Δ C i " and " current name damage vector
d i c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes in the formula 1
d i c Other outer amount is known, finds the solution formula 1 and just can calculate current name damage vector
d i c
J. the current actual damage vector that utilizes formula 2 to express
d i With the initial damage vector
d i o With current name damage vector
d i c Element between relation, calculate current actual damage vector
d i All elements;
In the formula 2
j=1,2,3 ..., N;
Current actual damage vector
d i Element numerical value represent the actual damage degree or the actual bearing angular displacement of corresponding evaluation object, according to current actual damage vector
d i Just can define the impaired and degree of injury of which rope, just can determine actual bearing angular 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, represent that the pairing rope of this element is intact, do not damage, if its numerical value is 100%, represent that then the pairing rope of this element has completely lost load-bearing capacity, 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 an angular displacement component of a bearing, so
d i j Represent its current angular displacement numerical value;
K. try to achieve current name damage vector
d i c After, set up mark vector according to formula 3
F i , formula 4 has provided mark vector
F i jThe definition of individual element;
Element in the formula 4
F i j It is mark vector
F i jIndividual element,
D i Uj It is nominal unit damage vector
D i u jIndividual element,
d i Cj It is current name damage vector
d i c jIndividual element, they all represent
jThe relevant information of individual evaluation object is in the formula 4
j=1,2,3 ..., N;
If mark vector l.
F i Element be 0 entirely, then get back to g step and continue this circulation; If mark vector
F i Element be not 0 entirely, then enter next step, i.e. m step;
M. according to formula 5 calculate next time, promptly
i+ 1 required initial damage vector of circulation
d i+
1 o Each element
d i+
1 Oj
In the formula 5
D i Uj It is nominal unit damage vector
D i u jIndividual element,
d i Cj It is current name damage vector
d i c jIndividual element,
F i j It is mark vector
F i jIndividual element is in the formula 5
j=1,2,3 ..., N; Vector
d i+
1 o The coding rule of element and the coding rule of evaluation object identical;
N. at Mechanics Calculation benchmark model A
iThe basis on, make the health status of evaluation object be
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 vector of required monitored amount promptly circulates for the i+1 time
C I+1 o
P. get back to the f step, beginning is circulation next time.
2. according to claim 1ly it is characterized in that in step f, at Mechanics Calculation benchmark model A based on the identification damaged cable of hybrid monitoring and the progressive method of bearing angular displacement
iThe basis on carry out the several times Mechanics Calculation, by calculate obtaining " unit damage monitored numerical quantity transformation matrices
Δ C i " and " nominal unit damage vector
D i u " concrete grammar be:
F1. at the Mechanics Calculation benchmark model A of Cable Structure
iThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value
NCoding rule according to evaluation object calculates successively; Calculating hypothesis each time has only an evaluation object to increase unit damage or unit angular displacement again on the basis of original damage or angular displacement, concrete, if this evaluation object is a support cable in the cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is the angular displacement component of a direction of a bearing, just suppose that this bearing increases the unit angular displacement again in this angular displacement direction, the evaluation object that increases unit damage or unit angular displacement in calculating each time again is different from the evaluation object that increases unit damage or unit angular displacement in other time calculating again, with " nominal unit damage vector
D i u " unit damage or the unit angular displacement that increase again of all supposition of record record; wherein i represents 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 monitored amount of current calculated value composition of the monitored amount of all that calculate is calculated current numerical value vector each time;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit angular 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 change vector, has
NIndividual evaluation object just has
NIndividual monitored quantitative change vector;
F3. by this
NIndividual monitored quantitative change vector according to
NThe coding rule of individual evaluation object, forming successively has
NThe monitored amount unit change of the Cable Structure matrix of row
Δ C i
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011101226143A CN102288424A (en) | 2011-05-13 | 2011-05-13 | Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011101226143A CN102288424A (en) | 2011-05-13 | 2011-05-13 | Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102288424A true CN102288424A (en) | 2011-12-21 |
Family
ID=45335024
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2011101226143A Pending CN102288424A (en) | 2011-05-13 | 2011-05-13 | Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102288424A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102706630A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Damaged cable and support angular displacement progressive identification method on basis of hybrid monitoring during temperature variation |
| CN102706637A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Damaged cable and support translation progressive identification method on basis of hybrid monitoring during temperature variation |
| CN102706658A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Defective cable and support angular displacement progressive identification method based on mixed monitoring of temperature change |
| CN102706623A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Slack cable progressive identification method on basis of hybrid monitoring during generalized displacement of support and temperature variation |
| CN102735481A (en) * | 2012-05-30 | 2012-10-17 | 东南大学 | Approximation type injured cable identification method based on hybrid monitoring on support angular displacement and temperature change |
| CN103852313A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Angular displacement mixed-monitoring load identification method for problematic cable |
| CN105004552A (en) * | 2015-07-23 | 2015-10-28 | 东南大学 | Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined linear displacement |
| CN105067358A (en) * | 2015-07-23 | 2015-11-18 | 东南大学 | Method for progressively recognizing load of damaged cable based on hybrid monitoring process of streamlined linear displacement |
| CN105115773A (en) * | 2015-07-23 | 2015-12-02 | 东南大学 | Identification method for problematic cable and load through simplified cable force monitoring |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006337144A (en) * | 2005-06-01 | 2006-12-14 | Kawasaki Heavy Ind Ltd | Fatigue life diagnostic method and diagnostic support device of bridge |
| JP2007297777A (en) * | 2006-04-27 | 2007-11-15 | Nippon Steel Engineering Co Ltd | Cable for suspension structure and measurement system |
| CN101476990A (en) * | 2009-01-20 | 2009-07-08 | 东南大学 | Progressive health monitoring method for cable system in cable structure based on angle monitoring |
| CN101819098A (en) * | 2010-03-31 | 2010-09-01 | 东南大学 | Mixed monitoring based progressive method for identifying damaged cables and support displacement |
| CN101832872A (en) * | 2010-03-31 | 2010-09-15 | 东南大学 | Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring |
-
2011
- 2011-05-13 CN CN2011101226143A patent/CN102288424A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006337144A (en) * | 2005-06-01 | 2006-12-14 | Kawasaki Heavy Ind Ltd | Fatigue life diagnostic method and diagnostic support device of bridge |
| JP2007297777A (en) * | 2006-04-27 | 2007-11-15 | Nippon Steel Engineering Co Ltd | Cable for suspension structure and measurement system |
| CN101476990A (en) * | 2009-01-20 | 2009-07-08 | 东南大学 | Progressive health monitoring method for cable system in cable structure based on angle monitoring |
| CN101819098A (en) * | 2010-03-31 | 2010-09-01 | 东南大学 | Mixed monitoring based progressive method for identifying damaged cables and support displacement |
| CN101832872A (en) * | 2010-03-31 | 2010-09-15 | 东南大学 | Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102706637B (en) * | 2012-05-29 | 2015-10-07 | 东南大学 | The damaged cable of hybrid monitoring and support translation progressive identification method during temperature variation |
| CN102706637A (en) * | 2012-05-29 | 2012-10-03 | 东南大学 | Damaged cable and support translation progressive identification method on basis of hybrid monitoring during temperature variation |
| CN102706623A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Slack cable progressive identification method on basis of hybrid monitoring during generalized displacement of support and temperature variation |
| CN102706630A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Damaged cable and support angular displacement progressive identification method on basis of hybrid monitoring during temperature variation |
| CN102735481A (en) * | 2012-05-30 | 2012-10-17 | 东南大学 | Approximation type injured cable identification method based on hybrid monitoring on support angular displacement and temperature change |
| CN102706658B (en) * | 2012-05-30 | 2015-07-08 | 东南大学 | Defective cable and support angular displacement progressive identification method based on mixed monitoring of temperature change |
| CN102706630B (en) * | 2012-05-30 | 2015-07-29 | 东南大学 | The damaged cable of temperature variation hybrid monitoring and support angular displacement progressive identification method |
| CN102735481B (en) * | 2012-05-30 | 2015-09-30 | 东南大学 | The damaged cable of hybrid monitoring approximant recognition methods during angular displacement of support temperature variation |
| CN102706658A (en) * | 2012-05-30 | 2012-10-03 | 东南大学 | Defective cable and support angular displacement progressive identification method based on mixed monitoring of temperature change |
| CN103852313A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Angular displacement mixed-monitoring load identification method for problematic cable |
| CN105004552A (en) * | 2015-07-23 | 2015-10-28 | 东南大学 | Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined linear displacement |
| CN105067358A (en) * | 2015-07-23 | 2015-11-18 | 东南大学 | Method for progressively recognizing load of damaged cable based on hybrid monitoring process of streamlined linear displacement |
| CN105115773A (en) * | 2015-07-23 | 2015-12-02 | 东南大学 | Identification method for problematic cable and load through simplified cable force monitoring |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101832874B (en) | Progressive cable structure health monitoring method based on angle monitoring | |
| CN102323084A (en) | Progressive method for identifying damaged cable and support generalized displacement on the basis of mixed monitoring | |
| CN101832872B (en) | Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring | |
| CN102288424A (en) | Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring | |
| CN102221478A (en) | Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support | |
| CN102297766A (en) | Progressive method for identifying damaged rope and angular displacement of support based on strain monitoring | |
| CN101806667B (en) | Progressive method for identifying damaged cable and support displacement based on cable force monitoring | |
| CN102252861A (en) | Progressive method for identifying support generalized displacement of damaged cable based on space coordinate monitoring | |
| CN101819098B (en) | Mixed monitoring based progressive method for identifying damaged cables and support displacement | |
| CN102221479A (en) | Progressive method for identifying damaged cable, loose cable and generalized displacement of supporting seat based on space coordinate monitoring | |
| CN102323075B (en) | Progressive method for identifying damaged cable, slack cable and support generalized displacement on basis of angle monitoring | |
| CN101793627B (en) | Health monitoring method of increment cable structure based on hybrid monitoring | |
| CN101832875B (en) | Progressive cable structure health monitoring method based on cable force monitoring | |
| CN102323085A (en) | Angle monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement | |
| CN102288423B (en) | Progressive method for recognizing damaged cable and angular displacement of support based on angle monitoring | |
| CN102323086A (en) | Strain monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement | |
| CN101819097B (en) | Progressive method for identifying damaged cable and support displacement based on strain monitoring | |
| CN102323074A (en) | The method of going forward one by one based on the generalized displacement of hybrid monitoring identification damaged cable slack line bearing | |
| CN102221480A (en) | Method for identifying damaged cable, loose cable and generalized displacement of supporting seat based on hybrid monitoring | |
| CN102297780A (en) | Progressive method for identifying damaged cables, loose cables and supporting-base displacement based on hybrid monitoring | |
| CN101832876A (en) | Progressive cable structure health monitoring method based on spatial coordinate monitoring | |
| CN101793629A (en) | Method for monitoring health of progressive cable structure based on strain monitor | |
| CN102288438B (en) | Method for recognizing damaged cable, slack cable and angular displacement of support based on hybrid monitoring | |
| CN102288440A (en) | Progressive method for recognizing damaged cable, slack cable and angular displacement of support based on space coordinate monitoring | |
| CN102297781A (en) | Progressive method for identifying damaged cable, slacking cable and bearer angular displacement based on strain monitoring |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20111221 |