Summary of the invention
technical matters:the objective of the invention is when the Cable Structure bearing has angular displacement, for in cable system in Cable Structure, identification problem that need the support cable of adjustment Suo Li, structure health monitoring method a kind of hybrid monitoring based on to the multiclass parameter, that can identify rationally and effectively the support cable that needs adjustment Suo Li is disclosed.
The reason changed according to the Suo Li of support cable can change the three kinds of situations that be divided into by the Suo Li of support cable: the one, and support cable has been subject to damage, and for example localized cracks and corrosion etc. have appearred in support cable; The 2nd, support cable not damaged, but variation has also occurred in Suo Li, the one of the main reasons that this variation occurs is that variation has occurred the Suo Changdu (be called drift, the present invention specially refers to the drift of that section rope between support cable two supporting end points) under support cable free state (now Suo Zhangli also claims that Suo Li is 0); The 3rd, support cable not damaged, but the Cable Structure bearing has had rotation, also can cause the variation of structural internal force, certainly also will cause the variation of Suo Li.One of fundamental purpose of the present invention is exactly when angular displacement of support, identify drift the support cable changed has occurred, and identify the change amount of their drift, and the Suo Li that this change amount is this rope adjusts provides direct basis.The reason that the support cable drift changes is not single, and for convenient, the support cable that the present invention changes drift is referred to as slack line.
technical scheme:the present invention is comprised of the two large divisions.Respectively: one, set up the required knowledge base of health monitoring systems for identifying support cable cable system, that need to adjust Suo Li and parameter method, adjust the method for the support cable of Suo Li based on knowledge base (containing parameter), need angular displacement of support, monitoring based on monitored amount equivalent based on the actual measurement Cable Structure, the identification Cable Structure; Two, the software and hardware part of health monitoring systems.
First of the present invention: set up the required knowledge base of health monitoring systems for identifying support cable cable system, that need to adjust Suo Li and parameter method, adjust the method for the support cable of Suo Li based on knowledge base (containing parameter), need angular displacement of support, monitoring based on monitored amount equivalent based on the actual measurement Cable Structure, the identification Cable Structure.Can be successively circularly as follows, laddering carrying out, to obtain the health status assessment of cable system more accurately.
The first step: when circulation starts each time, the initial virtual lesion vector of cable system while at first needing to set up or set up this circulation beginning
d o i (
i=1,2,3 ...Because in fact support cable may be lax and not damage is to mean difference, claim here " virtual lesion ", rear with), set up the initial mechanical calculating benchmark model A of Cable Structure
o(for example finite element benchmark model, A in the present invention
oconstant), set up the current Mechanics Calculation benchmark model A of Cable Structure
ti o(finite element benchmark model for example, A in circulation each time
ti oconstantly update), 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 that means significantly number of steps, and alphabetical i only means cycle index in the present invention, i.e. the i time circulation.
If total in cable system
nthe root rope, " initial virtual lesion vector is designated as the cable system that the i time circulation needs while starting
d o i " (as the formula (1)), use
d o i cable Structure while meaning this circulation beginning is (with Mechanics Calculation benchmark model A
ithe health status of cable system expression).
In formula (1)
d i oj (
i=1,2,3,
; j=1,2,3 ....,
n) while meaning that the i time circulation starts, Mechanics Calculation benchmark model A
iin cable system
jthe initial virtual lesion value of root rope,
d i oj be to mean at 0 o'clock
jroot rope not damaged is without lax, means while being 100% that this rope thoroughly loses load-bearing capacity, means the in the time of between 0 and 100%
jthe root rope is lost the load-bearing capacity of corresponding proportion.In formula (1)
tmean vectorial transposition (rear same).
Setting up the initial virtual lesion vector of cable system when circulation starts for the first time (is designated as according to formula (1)
d 1 o ) time, Non-destructive Testing Data that utilizes rope etc. can be expressed the data of the health status of rope and be set up the initial virtual lesion vector of cable system
d 1 o .If while there is no the data of the Non-destructive Testing Data of rope and other health status that can express rope, or can think that the structure original state is not damaged during without relaxed state, vector
d 1 o each element numerical value get 0.
The i time (
i=2,3,4,5,6 ...) the initial virtual lesion vector of cable system of needs while circulating beginning
d i o , be front once (the i-1 time,
i=2,3,4,5,6 ...) the front calculating acquisition of circulation end, concrete grammar is described below.
The Mechanics Calculation benchmark model that need to set up when the i time circulation starts or the Mechanics Calculation benchmark model of having set up are designated as A
i.
Design drawing, as-constructed drawing and the measured data of the Cable Structure in being completed according to Cable Structure (comprises that the Non-destructive Testing Data etc. of rope can express the measured datas such as the data of the health status of rope, Cable Structure shape data, structure angle degrees of 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 the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge), utilize mechanics method (for example finite element method) to set up A
o, if there is no the measured data of the structure in 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, the 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, utilize mechanics method (for example finite element method) to set up A
o.No matter which kind of method to obtain A by
o, based on A
othe Cable Structure computational data calculated (to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge) must approach its measured data very much, and error generally must not be greater than 5%.But utility A like this
ocalculate strain computational data, Suo Li computational data, Cable Structure shape computational data and displacement computational data, Cable Structure angle-data etc. under the analog case of gained, the measured data when approaching reliably institute's analog case and truly occurring.Corresponding to A
ocable Structure bearing angular data form initial Cable Structure bearing angular coordinate vector
u o.A
owith
u obe constant, only set up when circulation beginning for the first time.
The Mechanics Calculation benchmark model of the Cable Structure of setting up when circulation starts for the first time is designated as A
1, A
1just equal A
o.A
1the health status of corresponding rope by
d 1 o describe.
The i time (
i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs while circulating beginning
i, be front once (the i-1 time,
i=2,3,4,5,6 ...) the front calculating acquisition of circulation end, concrete grammar is described below.
Existing Mechanics Calculation benchmark model A
1with the initial virtual lesion vector of cable system
d 1 o after, model A
1in the virtual lesion of each rope by vector
d 1 o express.At A
1basis on, the virtual lesion value of all ropes is changed to 0, mechanical model A
1the virtual lesion that is updated to all ropes is all that 0 mechanical model (is designated as A
0), mechanical model A
0be actually the excellent mechanical model corresponding without lax Cable Structure.Might as well claim model A
0for the not damaged of Cable Structure without relaxation model A
0.
Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate are described below respectively:
If total in cable system
nthe root rope, the monitored rope force data of structure is by structure
m 1 individual appointment rope
m 1 individual rope force data is described, and the variation of structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.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 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 all tested strains 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 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 the variation of angle coordinate components appointment straight lines all specified points, all, 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 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 rope
n.Due to
mindividual monitored amount is dissimilar, so the present invention is called " method of identifying the support cable that needs adjustment Suo Li based on hybrid monitoring ".
For simplicity, in the present invention by " monitored all parameters of structure " referred to as " monitored amount ".
" the initial value vector of monitored amount for the present invention
c i o " (
i=1,2,3 ...) mean the i time (
i=1,2,3,4,5,6 ...) initial value (referring to formula (2)) of the monitored amount of all appointments while circulating beginning,
c i o full name be the initial value vector of monitored amount " the i time circulation ".
In formula (2)
c i ok (
i=1,2,3,
k=1,2,3,
., M; M>=N;) be that the i time circulation is while starting, in Cable Structure the
kindividual monitored amount.Vector
c i o by previously defined
mindividual monitored amount is arranged and is formed according to a definite sequence, this is put in order and there is no specific (special) requirements, only require all associated vector in back also in this order array data get final product.
When circulation starts for the first time, " the initial value vector of the 1st the monitored amount that circulates
c 1 o " (seeing formula (2)) be comprised of measured data, due to according to model A
1calculate the initial value of the monitored amount of gained reliably close to corresponding measured value, in the narration of back, will mean 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 while circulating beginning
c i o ", be front once (the i-1 time,
i=2,3,4,5,6 ...) the front calculating acquisition of circulation end, concrete grammar is described below.
Second step: in Cable Structure military service process, in circulation each time, constantly actual measurement obtains Cable Structure bearing angular coordinate current data (all data composition current cable structure actual measurement bearing angular coordinate vectors
u ti , vector
u ti definition mode with the vector
u oidentical).For simplicity, for the i time circulation, the Cable Structure bearing angular coordinate current data when last time is upgraded to current Mechanics Calculation benchmark model is designated as current cable structural bearings angular coordinate vector
u ti o.Set up and renewal A
ti omethod be: each time the circulation the zero hour, the current Mechanics Calculation benchmark model A of Cable Structure
ti ojust equal A
i(
i=1,2,3,4,5,6 ...).In Cable Structure military service process, constantly actual measurement obtains Cable Structure bearing angular data and obtains current cable structure actual measurement bearing angular coordinate vector
u ti if,
u ti equal
u ti o, do not need A
ti oupgraded; If
u ti be not equal to
u ti o, need A
ti oupgraded, now
u ti with
u odifference be exactly the Cable Structure bearing about initial position (corresponding to A
o) angular displacement of support (use the angular displacement of support vector
vmean angular displacement of support).Upgrade A
ti omethod be: at A
obasis on to make the health status of rope be cable system initial damage vector
d i o , more further to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support, the numerical value of current angular displacement of support constraint is just taken from current angular displacement of support vector
vthe numerical value of middle corresponding element, to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
ti o, upgrade A
ti oafter,
u ti oall elements numerical value is used
u ti all elements numerical value replaces, and has upgraded
u ti o, so just obtained correctly corresponding to A
ti o's
u ti o.
The 3rd step: circulation each time needs to set up " the monitored numerical quantity transformation matrices of virtual unit damage " and " nominal virtual unit damage vector ",
," the monitored numerical quantity transformation matrices of virtual unit damage " that the i time circulation set up is designated as
Δ C i (
i=1,2,3 ...)." nominal virtual unit damage vector " that the i time circulation set up is designated as
d i u .In circulation each time
Δ C i with
d i u constantly update, upgrading current Mechanics Calculation benchmark model A
ti othe time, upgrade the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i with nominal virtual unit damage vector
d i u .
Set up and the monitored numerical quantity transformation matrices of renewal virtual unit damage
Δ C i with nominal virtual unit damage vector
d i u process as follows:
Current Mechanics Calculation benchmark model A in Cable Structure
ti obasis on carry out several times calculating, equal the quantity of all ropes on calculation times numerical value.Calculating each time in the hypothesis cable system only has a rope to increase virtual unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is virtual unit damage) on the basis of original virtual lesion (original virtual lesion can be 0, can not be also 0) again.For convenience of calculating, while setting virtual unit damage in the circulation each time, can be all structural health conditions when this circulation started as being fully healthy, and set on this basis virtual unit damage (in subsequent step, damage numerical value that calculate, rope---be called nominal virtual lesion
d i c (
i=1,2,3 ...), all when this circulation started, by the health status of rope as being fully healthy speech, the formula therefore must foundation hereinafter provided be converted into actual and virtual by the nominal virtual lesion calculated and damage).Occur in calculating each time with once circulation that the rope of virtual lesion is different from the rope of the virtual lesion occurred in other calculating, and supposition each time has the virtual unit damage value of the rope of virtual lesion can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
d i u " (as the formula (3)) record the virtual unit damage of supposition of all ropes in each 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, calculate gained each time
mthe current calculated value of individual monitored amount forms one " the calculating current value vector of monitored amount " (when hypothesis the
jwhen the root rope has unit damage, available formula (4) means all appointments
mthe calculating current value vector of individual monitored amount
c 1 tj ); The calculating current value vector of the monitored amount calculated 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 under this condition, (take and the position of the rope of virtual unit damage or numbering etc. be arranged as mark) is (when
jwhen the root rope has virtual unit damage, use
δ C 1 j the numerical value change vector that means monitored amount,
δ C 1 j definition see formula (5), formula (6) and formula (7), formula (5) deducts after 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 the supposition owing to calculating of each element representation the Na Gensuo (for example the of virtual unit damage is arranged
jthe root rope) virtual unit damage (for example
d uj ), and the numerical value change amount of the corresponding monitored amount of this element caused is with respect to the virtual unit damage of supposition
d uj rate of change; Have
nthe root rope 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 " forms and has successively
m * N" the unit damage monitored quantitative change matrix of individual element
Δ C 1 "
(MoK
nrow
), each vector
δ C 1 j (
j=1,2,3 ....,
n) be matrix
Δ C 1 row,
Δ C 1 definition as the formula (8).
(3)
Nominal virtual unit damage vector in formula (3)
d i u element
d i uj (
i=1,2,3,
j=1,2,3 ....,
n) mean
iin inferior circulation, suppose
jthe virtual unit damage numerical value of root rope, vector
d i u in the numerical value of each element can be the same or different.
Element in formula (4)
c i tjk (
i=1,2,3,
; j=1,2,3 ....,
n; k=1,2,3 ....,
m; M>=N) mean
iinferior circulation is due to
jwhen the root rope has virtual unit damage, according to coding rule corresponding
kthe calculating current value of the monitored amount of individual appointment.
The subscript of each amount in formula (5)
i(
i=1,2,3 ...) mean
iinferior circulation, subscript
j(
j=1,2,3 ....,
n) mean
jthe root rope has virtual unit damage, in 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 mean the
iin inferior circulation, set up matrix
Δ C i the time, suppose
jwhen having virtual unit damage, the root rope calculates gained the
kthe change amount of individual monitored amount is with respect to the virtual unit damage of supposition
d i uj rate of change, it defines suc as formula shown in (7).
In formula (7), the definition of each amount has been previously described.
(8)
Vector in formula (8)
δ C i j (
i=1,2,3 ....,,
j=1,2,3 ....,
n) mean
iin inferior circulation, due to
jthe root rope has virtual 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.
In Cable Structure military service process, in circulation each time, constantly actual measurement obtains Cable Structure bearing angular coordinate current data, once monitor
u ti be not equal to
u ti o, need to get back to second step to A
ti oupgraded, to A
ti oenter again this step after being upgraded right
Δ C i upgraded.In fact in circulation each time
Δ C i constantly update, upgrading current Mechanics Calculation benchmark model A
ti oafterwards, upgrade the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i .
The 4th step: the current health status of identification cable system.Detailed process is as follows.
The
i(
i=1,2,3 ...) in inferior circulation, cable system " current (calculating or actual measurement) numerical value vector of monitored amount
c i " " the initial value vector of monitored amount together
c i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion 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 formula (10)
c i definition be similar to the initial value vector of monitored amount
c i o definition, see formula (11); Cable system " current nominal virtual lesion vector
d i c " definition see formula (12).
Element in formula (11)
c i k (
i=1,2,3 ....;
k=1,2,3 ....,
m; M>=N) be
iinferior circulation time Cable Structure, according to coding rule is corresponding, be numbered
kthe current value of monitored amount.
In formula (12)
d i cj (
i=1,2,3 ....;
j=1,2,3 ....,
n) be
iin inferior circulation, cable system
jthe current nominal virtual lesion value of root rope, vector
d i c the subscript of element
jcoding rule and matrix
Δ C i the coding rule of row identical.
When the rope actual damage is not too large, 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 formula (10).
In formula (13)
abs ()be the function that takes absolute value, each element of the vector of trying to achieve in 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 " come direct solution to obtain " current nominal virtual lesion vector
d i c ".If done like this, the vector obtained
d i c in element even there will be larger negative value, namely negative damage, this is obviously irrational.Therefore obtain vector
d i c acceptable solution (with reasonable error, but can determine more exactly position and the virtual lesion degree thereof of virtual damaged cable) become a rational solution, available formula (14) is expressed this method.
(14)
In formula (14)
abs ()the function that takes absolute value, vector
g i description departs from ideal linearity relation (formula (9) or formula (10))
Legitimate skew, by formula (15), defined.
In 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 formula (10) in inferior circulation.Vector
g i can be according to the error vector of formula (13) definition
e i tentative calculation is selected.
At " the initial value vector of monitored amount
c i o " (survey or calculate), " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " (calculating) and " the current value vector of monitored amount
c i " (actual measurement obtains) when known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to solve formula (14), obtains " current nominal virtual lesion vector
d i c " acceptable solution, " current actual virtual lesion vector then
d i " element of (formula (16) is shown in definition) can calculate according to formula (17), namely obtained " " current actual virtual lesion vector
d i thereby, can by
d i determine position and the virtual lesion degree of virtual damaged cable, then according to below the method for narration being determined to position and the relax level of slack line, namely determined the rope and the long adjustment amount of rope thereof that need to adjust Suo Li.
In formula (16)
d i j (
i=1,2,3,
; j=1,2,3 ....,
n) mean
iin inferior circulation
jthe actual virtual lesion value of root rope, formula (17) is shown in its definition,
d i j be to mean at 0 o'clock
jroot rope not damaged is without lax, means while being 100% that this rope thoroughly loses load-bearing capacity, means the in the time of between 0 and 100%
jthe root rope is lost the load-bearing capacity of corresponding proportion, vector
d i the coding rule of element and formula (1) in vector
d i o the coding rule of element identical.
In 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.
Below narration has obtained the actual virtual lesion vector of Suo Dangqian
d i after, the position of how to confirm slack line and relax level.
If total in cable system
nthe root support cable, structure rope force data by
nthe Suo Li of root support cable describes.Available " initial rope force vector
f o " mean the initial Suo Li (formula (18) is shown in definition) of all support cables in Cable Structure.Because the calculating benchmark model based on Cable Structure calculates the initial Suo Li of gained reliably close to the measured data of initial Suo Li, in the narration of back, will mean this calculated value and measured value with prosign.
In formula (18)
f oi (
j=1,2,3,
., N) be in Cable Structure
jthe initial Suo Li of root support cable, this element is the Suo Li corresponding to the appointment support cable according to coding rule.Vector
f o constant, irrelevant with cycle index, after determining when circulation starts for the first time, just no longer change.Setting up the Mechanics Calculation benchmark model A of Cable Structure
1the time used vector
f o .
Use " current cable force vector in the present invention
f i " mean
ithe current cable power (formula (19) is shown in definition) of all support cables in the Cable Structure that inferior circulation time actual measurement obtains.
In formula (19)
f i j (
i=1,2,3,4,
; j=1,2,3,
., N) be
iin inferior circulation time Cable Structure
jthe current cable power of root support cable.
In the present invention, under support cable original state (not damaged, without lax), and support cable is when free state (free state refers to that Suo Li is 0, rear with), and the length of support cable is called initial drift, with " initial drift vector
l o " mean the initial drift (formula (20) is shown in definition) of all support cables in Cable Structure.
In formula (20)
l oj (
j=1,2,3,
., N) be in Cable Structure
jthe initial drift of root support cable.Vector
l o constant, irrelevant with cycle index, after determining when circulation starts for the first time, just no longer change.
In the present invention, with " current drift vector
l i " mean
ithe current drift of all support cables in inferior circulation time Cable Structure (formula (21) is shown in definition).
In formula (21)
l i j (
i=1,2,3,4,
; j=1,2,3,
., N) be
iin inferior circulation time Cable Structure
jthe current drift of root support cable.
In the present invention, with " drift changes vectorial Δ
l i " (or claim support cable current relax level vector) mean the
ithe change amount of the drift of all support cables in inferior circulation time Cable Structure (formula (22) and formula (23) are shown in definition).
(22)
Δ in formula (22)
l i j (
i=1,2,3,4,
; j=1,2,3,
., N) be current (
iinferior circulation time) in Cable Structure
jthe change amount of the drift of root support cable, formula (23), Δ are shown in its definition
l i j be not that 0 rope is slack line, Δ
l i j the numerical value slack that is rope, and mean cable system the
jthe current relax level of root support cable is also the long adjustment amount of rope of this rope while adjusting Suo Li.
By slack line is carried out to the relax level identification that slack line is carried out in the mechanics equivalence with damaged cable, the mechanical condition of equivalence is in the present invention:
The mechanics parameters of initial drift, geometrical property parameter and material when one, the nothing of the rope of two equivalences relaxes with not damaged is identical;
Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.
While meeting above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if, after with equivalent damaged cable, replacing slack line, Cable Structure any variation can not occur, vice versa.
In the present invention, the
iinferior circulation time, with the
jindividual support cable (its current relax level Δ
l i j the current actual virtual lesion degree of definition) carrying out the virtual impaired support cable of equivalence is used
d i j mean (
d i j definition see formula (16) and formula (17)).Lax the
jthe current relax level Δ of individual support cable
l i j (Δ
l i j definition see formula (22)) with the current actual virtual lesion degree of equivalent damaged cable
d i j between relation by aforementioned two mechanics equivalent conditions, determined.Δ
l i j with
d i j between physical relationship can adopt accomplished in many ways, for example can be directly according to aforementioned equivalent condition, determine (referring to formula (24)), also can adopt based on the Ernst equivalent elastic modulus and replace in formula (24)
edetermine after revising (referring to formula (25)), also can adopt other methods such as trial and error procedure based on finite element method to determine.
In formula (24) and formula (25)
ethe elastic modulus of this support cable,
athe cross-sectional area of this support cable,
f i j the current cable power of this support cable,
d i j the current actual virtual lesion degree of this support cable,
ω i the weight of the unit length of this support cable,
l i jx it is the horizontal range of two supporting end points of this support cable.Item in formula (25) in [] is the Ernst equivalent elastic modulus of this support cable, by formula (24) or formula (25), can just can determine the current relax level vector of support cable Δ
l i .Formula (25) is the correction to formula (24).
The 5th step: judge whether to finish this (
iinferior) circulation, if so, complete the tailing in work before this circulation finishes, for next time (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) circulation in try to achieve " current nominal virtual lesion vector
d i c " after, at first, according to formula (26), set up " mark vector
b i ", formula (27) has provided " mark vector
b i "
jthe definition of individual element; If " mark vector
b i " element be 0 entirely, in this circulation, continue health monitoring and the calculating to cable system; If " mark vector
b i " element be not 0 entirely, after completing subsequent step, enter next time circulation.So-called subsequent step is: at first, according to formula (28), calculate next time (
i+ 1 time,
i=1,2,3,4 ...) required " initial virtual lesion vector circulates
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
0basis on, make the health status of rope be
d i+
1 o rear renewal obtain next time (the i+1 time,
i=1,2,3,4 ...) circulate Mechanics Calculation benchmark model A required
i+1; Finally, by Mechanics Calculation benchmark model A
i+1the initial value that calculates monitored amount, consisting of next time (the i+1 time,
i=1,2,3,4 ...) required " the initial value vector of monitored amount circulates
c i+1 o " (
i=1,2,3,4 ...).
Mark vector in formula (26)
b i subscript
imean the
iinferior circulation, its element
b i j (
j=1,2,3 ..., subscript N)
jmean the
jthe damage characteristic of root rope, can only get 0 and 1 two amount, and concrete value rule is shown in formula (27).
Element in formula (27)
b i j it is " mark vector
b i "
jindividual element,
d i uj it is " nominal virtual unit damage vector
d i u "
jindividual element (seeing formula (3)),
d i cj it is " current nominal virtual lesion vector
d i c "
jindividual element (seeing formula (12)), they all mean
jthe relevant information of root rope.
In formula (28)
d i uj it is " nominal virtual unit damage vector
d i u "
jindividual element (seeing formula (3)),
d i cj it is " current nominal virtual lesion 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 monitoring system (monitoring the horizontal range of monitored amount, monitoring Cable Structure angular displacement of support, monitoring Suo Li, monitoring support cable two supporting end points), signal picker and computing machine etc.Requirement is monitored in real time or quasi real time each monitored amount, monitors the Suo Li of each support cable, is monitored the horizontal range that each support cable two supports end points.Software should the following function of tool: software section should complete the process that first of the present invention sets, complete needed in the present invention, can be by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
The inventive method specifically comprises:
A. establish total N root rope, at first determine the coding rule of rope, by this rule, by rope numberings all in Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step;
While b. determining hybrid monitoring appointment by the support cable of monitored Suo Li, establish in cable system total
nthe root rope, the monitored rope force data of structure is by structure
m 1 individual appointment rope
m 1 individual rope force data is described, and the variation of structure Suo Li is exactly the variation of the Suo Li of all appointment ropes; 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; While determining hybrid monitoring appointment by the measured point of monitored strain, the monitored strain data of structure can be by 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 all tested strains 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 to be not less than 0 integer; While determining hybrid monitoring appointment by the measured point of monitored angle, the monitored angle-data of structure is by 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 the variation of angle coordinate components appointment straight lines all specified points, all, 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; While determining hybrid monitoring appointment by monitored shape data, the monitored shape data of structure is by 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 rope
n; Due to
mindividual monitored amount is dissimilar, so the present invention is called " health monitor method of the cable system based on hybrid monitoring during angular displacement of support "; For simplicity, in the present invention by this step listed " monitored all parameters of structure during hybrid monitoring " referred to as " monitored amount ";
C. Non-destructive Testing Data that utilizes rope etc. can be expressed the data of the health status of rope and be set up initial virtual lesion vector
d i o , wherein i means cycle index, back i and subscript i mean cycle index, and i=1,2,3, Circulation time for the first time
d i o be designated as
d 1 o .If while there is no the data of the Non-destructive Testing Data of rope and other health status that can express rope, or can think that the structure original state is for without lax, not damaged state the time, vector
d 1 o each element numerical value get 0.
D. setting up initial virtual lesion 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 ;
E. setting up initial virtual lesion vector
d 1 o initial value vector with monitored amount
c 1 o the time, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector
f o ; Simultaneously, obtain the initial drift of all support cables according to structural design data, completion data, form initial drift vector
l o ; Simultaneously, obtain the initial geometric data of Cable Structure according to structural design data, completion data or actual measurement; Simultaneously, survey or obtain according to structural design, completion information elastic modulus, density, the initial cross sectional area of all ropes;
F. set up the initial mechanical calculating benchmark model A of Cable Structure
o, set up initial Cable Structure bearing angular coordinate vector
u o, the Mechanics Calculation benchmark model A of the Cable Structure needed while setting up circulation beginning for the first time
1; The measured data of the Cable Structure in being completed according to Cable Structure, this measured data comprises the measured datas such as the elastic modulus, density, initial cross sectional area of Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure bearing angular data, Cable Structure modal data, all ropes, and the Non-destructive Testing Data of rope etc. can be expressed the data of the health status of rope, according to design drawing and as-constructed drawing, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure
o; If there is no the measured data of the structure in Cable Structure completion, so just before setting up health monitoring systems, this Cable Structure is surveyed, obtain equally the measured data of Cable Structure, according to design drawing, the as-constructed drawing of these data and Cable Structure, utilize equally mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure
o; No matter which kind of method to obtain A by
o, based on A
othe Cable Structure computational data calculated must approach its measured data very much, and difference therebetween must not be greater than 5%; Corresponding to A
ocable Structure bearing angular data form initial Cable Structure bearing angular coordinate vector
u o; A
owith
u obe constant, only set up when circulation beginning for the first time; The Mechanics Calculation benchmark model of the Cable Structure of setting up when the i time circulation starts is designated as A
i, wherein i means cycle index; In application form of the present invention, alphabetical i is except the place that means significantly number of steps, and alphabetical i only means cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation starts for the first time is designated as A
1, A in the present invention
1just equal A
o; For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A
ti o", A in circulation each time
ti ocan constantly update as required, when circulation starts each time, A
ti oequal A
i; Equally for sake of convenience, name " Cable Structure actual measurement bearing angular coordinate vector
u ti ", in circulation each time, constantly actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas form current cable structure actual measurement bearing angular coordinate vector
u ti , vector
u ti element with the vector
u othe angular coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, for the i time circulation, the last time is upgraded to A
ti othe time Cable Structure bearing angular coordinate current data be designated as current cable structural bearings angular coordinate vector
u ti o; When circulation starts for the first time, A
t1 oequal A
1,
u t1 oequal
u o; A
1the health status of corresponding rope by
d 1 o describe; Mechanics Calculation benchmark model A
ithe health status of corresponding rope by
d i o describe;
When g. circulation starts each time, make A
ti oequal A
i; Actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas form current cable structure actual measurement bearing angular coordinate vector
u ti , according to current cable structure actual measurement bearing angular coordinate vector
u ti , upgrade where necessary the current Mechanics Calculation benchmark model of Cable Structure A
ti owith current cable structural bearings angular coordinate vector
u ti o;
H. at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on carry out the several times Mechanics Calculation, obtain the monitored numerical quantity transformation matrices of Cable Structure virtual unit damage by calculating
Δ C i with nominal virtual unit damage vector
d i u ;
I. actual measurement obtains the current cable power of all support cables of Cable Structure, forms the current cable force vector
f i ; Simultaneously, actual measurement obtains the current measured value of the monitored amount of all appointments of Cable Structure, forms " the current value vector of monitored amount
c i ".Actual measurement obtains the volume coordinate of two supporting end points of all support cables, and the volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.While numbering to the element of the institute's directed quantity occurred before this step and this step, should use same coding rule, each vector that can guarantee like this before this step and this step and occur afterwards, number identical element, mean same monitored amount, corresponding to vectorial defined relevant information under this element;
J. define current nominal virtual lesion vector to be asked
d i c with current actual virtual lesion vector
d i .The damage vector
d i o ,
d i c with
d i element number equal the quantity of rope, damaging between vectorial element and rope is one-to-one relationship, damages virtual lesion degree or health status that vectorial element numerical value represents corresponding rope;
K. according to " the current value vector of monitored amount
c i " " the initial value vector of monitored amount together
c i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, in formula 1, removes
d i c other outer amount is known, solves formula 1 and just can calculate current nominal virtual lesion vector
d i c ;
L. the current actual virtual lesion vector that utilizes formula 2 to express
d i element
d i j with initial virtual lesion vector
d i o element
d i oj with current nominal virtual lesion vector
d i c element
d i cj between relation, calculate current actual virtual lesion vector
d i all elements.
In formula 2
j=1,2,3 ..., N.
Due to current actual virtual lesion vector
d i element numerical value represent the current actual virtual lesion degree of corresponding rope, i.e. actual relax level or actual damage degree, current actual virtual lesion vector
d i middle numerical value is not that 0 support cable corresponding to element is exactly problematic support cable, and problematic support cable may be slack line, also may be damaged cable, its numerical response the degree of lax or damage;
M. identify damaged cable in the problematic support cable identified from the l step, remaining is exactly slack line.
N. utilize the current actual virtual lesion vector obtained in the l step
d i obtain the current actual virtual lesion degree of slack line, utilize the current cable force vector obtained in the i step
f i , utilize two volume coordinates that support end points in all support cables of i step acquisition, utilize the initial drift vector obtained in the e step
l o elastic modulus, density, the initial cross sectional area data of all ropes that utilization obtains in the e step, by by slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual virtual lesion degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.While meeting above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if, after with equivalent slack line, replacing damaged cable, Cable Structure any variation can not occur, vice versa.Try to achieve according to aforementioned mechanics equivalent condition the relax level that those are judged as slack line, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li.Lax identification and the damage identification of support cable have so just been realized.During calculating, institute's demand power is by the current cable force vector
f i corresponding element provides.
O. try to achieve current nominal virtual lesion vector
d i c after, set up mark vector according to formula 3
b i , formula 4 has provided mark vector
b i jthe definition of individual element;
Element in formula 3, formula 4
b i j it is mark vector
b i jindividual element,
d i uj it is nominal virtual unit damage vector
d i u jindividual element,
d i cj it is current nominal virtual lesion vector
d i c jindividual element, they all mean
jthe relevant information of root rope.In formula 4
j=1,2,3 ..., N.
If mark vector p.
b i element be 0 entirely, get back to g step and continue this circulation; If mark vector
b i element be not 0 entirely, enter next step, i.e. q step.
Q. according to formula 5 calculate next time,
i+ 1 required initial virtual lesion vector of circulation
d i+
1 o each element
d i+
1 oj ;
In formula 5
d i uj ithe nominal virtual unit damage vector of inferior circulation
d i u jindividual element,
d i cj ithe current nominal virtual lesion vector of inferior circulation
d i c jindividual element,
b i j iinferior circulation mark vector
b i jindividual element.In formula 5
j=1,2,3 ..., N.
R. at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on, make the health status of rope be
d i+
1 o rear renewal obtains next time, required Mechanics Calculation benchmark model A circulates for the i+1 time
i+1;
S. pass through Mechanics Calculation benchmark model A
i+1calculate corresponding to model A
i+1the numerical value of all monitored amounts of structure, these numerical value form next time, the initial value vector of required monitored amount circulates for the i+1 time
c i+1 o ;
T. set up the required current Mechanics Calculation benchmark model of the Cable Structure A that next time, circulates for the i+1 time
ti+1 o, get A
ti+1 oequal A
i+1;
U. set up the required current cable structural bearings angular coordinate vector that next time, circulates for the i+1 time
u ti+
1 o, get
u ti+
1 oequal
u ti o;
V. get back to step g, start circulation next time.
In step g, according to current cable structure actual measurement bearing angular coordinate vector
u ti , upgrade where necessary the current Mechanics Calculation benchmark model of Cable Structure A
ti owith current cable structural bearings angular coordinate vector
u ti oconcrete grammar be:
G1. actual measurement obtains current cable structure actual measurement bearing angular coordinate vector
u ti after, relatively
u ti with
u ti oif,
u ti equal
u ti o, do not need A
ti oupgraded;
G2. actual measurement obtains current cable structure actual measurement bearing angular coordinate vector
u ti after, relatively
u ti with
u ti oif,
u ti be not equal to
u ti o, need A
ti oupgraded, update method is: first calculate
u ti with
u opoor,
u ti with
u odifference be exactly that the current cable structural bearings is about setting up A
othe time the current angular displacement of support of Cable Structure bearing, with current angular displacement of support vector
vmean angular displacement of support, current angular displacement of support vector
vin element and the angular displacement of support component between be one-to-one relationship, current angular displacement of support vector
vin the numerical value of an element corresponding to the rotation of an assigned direction of an appointment bearing; Upgrade A
ti omethod be: at A
obasis on to make the health status of rope be cable system initial damage vector
d i o , more further to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support, the numerical value of current angular displacement of support constraint is just taken from current angular displacement of support vector
vthe numerical value of middle corresponding element, to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
ti o, upgrade A
ti othe time,
u ti oall elements numerical value is also used
u ti all elements numerical value replaces, and has upgraded
u ti o, so just obtained correctly corresponding to A
ti o's
u ti o.
In step h, at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on carry out the several times Mechanics Calculation, obtain the monitored quantitative change matrix of Cable Structure virtual unit damage by calculating
Δ C i with nominal virtual unit damage vector
d i u concrete grammar be:
H1. when the i time circulation starts, directly press step h2 to the monitored quantitative change matrix of the listed method acquisition Cable Structure virtual unit damage of step h4
Δ C i with nominal virtual unit damage vector
d i u ; The moment started in non-the i time circulation, in step g to A
ti oafter being upgraded, directly press step h2 to the monitored quantitative change matrix of the listed method acquisition Cable Structure virtual unit damage of step h4
Δ C i with nominal virtual unit damage vector
d i u ; The moment started in non-the i time circulation, if in step g not to A
ti oupgraded, directly proceed to herein step I and carry out follow-up work;
H2. at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on calculation times numerical value, have
nthe root rope just has
ninferior calculating, calculating each time in the hypothesis cable system only has a rope to increase virtual unit damage on the basis of original virtual lesion again, the rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other calculating, and supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
d i u " record the unit damage of the supposition of all ropes, calculate each time the current value of all monitored amounts, the current value of all monitored amounts that calculate each time forms one " the calculating current value vector of monitored amount ".When hypothesis the
jwhen the root rope has unit damage, available
c i tj mean corresponding " the current evaluation vector of monitored amount
c i tj ".While giving in this step each vectorial element numbering, should use same coding rule with other vector in the present invention, can guarantee like this any one element in each vector in this step, with in other vector, number identical element, expressed the relevant information of same monitored amount or same target;
H3. that calculating each time " current evaluation vector of monitored amount
c i tj " deduct " the initial value vector of monitored amount
c i o " obtain a vector, then obtain one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition during each element of this vector is calculated divided by this; Have
nthe root rope just has
nindividual " the numerical value change vector of monitored amount ";
H4. by this
nindividual " the numerical value change vector of monitored amount " forms and has successively
n" the monitored numerical quantity transformation matrices of virtual unit damage of row
Δ C i "; " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " each row corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector
d i c with current actual virtual lesion vector
d i the element coding rule identical.
beneficial effect:system and method disclosed by the invention in the situation that the Cable Structure bearing angular displacement occurs, having under the synchronous impaired or lax condition of more rope unusual monitor evaluate the health status (position and relax level or the degree of injury that comprise all slack lines and damaged cable) of cable system.This is due to " the current value vector of monitored amount
c i " " the initial value vector of monitored amount together
c i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between linear relationship be similar to; be actually nonlinear relation; when particularly the or extent of damage more at damaged cable is larger; the nonlinear characteristic of the relation between above-mentioned amount is more obvious; for overcoming this obstacle, the invention discloses a kind of health monitor method that approaches this nonlinear relationship in minizone by linear relationship.In fact the present invention has used the method by linear relationship piecewise approximation nonlinear relationship, large interval is divided into to minizone one by one, in each minizone internal linear relation, be enough accurately, the health status of the cable system obtained according to its judgement is also reliably, and therefore system and method disclosed by the invention is very useful to effective health monitoring of cable system.
Embodiment
When angular displacement of support, for the health monitoring of the cable system of Cable Structure, the invention discloses a kind of system and method for health status of each root rope of the cable system that can monitor rationally and effectively Cable Structure.The following describes of embodiments of the invention is in fact only exemplary, and purpose never is to limit application of the present invention or use.
In the situation that angular displacement appears in the Cable Structure bearing, the present invention adopts a kind of algorithm, and this algorithm is for monitoring the health status (relax level and the extent of damage that comprise rope) of cable system of Cable Structure.During concrete enforcement, the following step is a kind of in the various steps that can take.
The first step: determine type, position and the quantity of monitored amount, and numbering.Detailed process is:
If total N root rope, the coding rule of at first definite rope, by this rule, by rope numberings all in Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step.
While determining hybrid monitoring appointment by the support cable of monitored Suo Li, establish in cable system total
nthe root rope, the monitored rope force data of structure is by structure
m 1 individual appointment rope
m 1 individual rope force data is described, and the variation of structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.Each total
m 1 individual cable force measurement value or calculated value characterize the rope force information of structure.
m 1 be one and be not less than 0 integer.When reality is selected the rope of monitored Suo Li, the rope that can select those Suo Li to be easy to measure is monitored rope.
While determining hybrid monitoring appointment by the measured point of monitored strain, the monitored strain data of structure can be by 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 all tested strains 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.Near for example, a point each fixed endpoint that can be exactly each root rope by the measured point of monitored strain (being the stiff end of drag-line on bridge of cable-stayed bridge), this point should not be generally stress concentration point, to avoid 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.
While determining hybrid monitoring appointment by the measured point of monitored angle, the monitored angle-data of structure is by 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 the variation of angle coordinate components appointment straight lines all specified points, all, 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 fixed endpoint (being for example the stiff end of drag-line on bridge floor of cable-stayed bridge) or a near point it of each root rope, and the point of monitored angle-data generally should all not be chosen as " fixed endpoint of the rope of the monitored Suo Li of appointment or point in its vicinity in hybrid monitoring " and " point of the monitored strain of appointment or point in its vicinity in 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, be exactly in fact measurement of dip angle here.
While determining hybrid monitoring appointment by monitored shape data, the monitored shape data of structure is by 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 (being for example the stiff end of drag-line on bridge of cable-stayed bridge) of each root rope; Here " fixed endpoint of the rope of the monitored Suo Li of appointment or point in its vicinity in hybrid monitoring ", " point of the monitored strain of appointment or point in its vicinity in hybrid monitoring " and " point of the monitored angle-data of appointment or point in its vicinity in hybrid monitoring " should all not selected in selected point being monitored.
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 rope
n.Due to
mindividual monitored amount is dissimilar, so the present invention is called " method of identifying the support cable that needs adjustment Suo Li based on hybrid monitoring ".For simplicity, in the present invention by this step listed " monitored all parameters of structure during hybrid monitoring " referred to as " monitored amount ".
Second step: Non-destructive Testing Data that utilizes rope etc. can be expressed the data of the health status of rope and be set up initial virtual lesion vector
d 1 o .If while there is no the data of the Non-destructive Testing Data of rope and other health status that can express rope, or can think the structure original state be not damaged, during without relaxed state, vector
d 1 o each element numerical value get 0.
The 3rd step: setting up initial virtual lesion 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 "; Simultaneously, directly measure the initial Suo Li of all support cables that calculate Cable Structure, form " initial rope force vector
f o "; Simultaneously, obtain the initial drift of all ropes according to structural design data, completion data, form " the initial drift vector of support cable
l o "; Simultaneously, survey or obtain according to structural design, completion information elastic modulus, density, the initial cross sectional area of all ropes.
The 4th step: setting up initial virtual lesion vector
d 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.Calculate Cable Structure original geometric form data (being exactly its initial bridge type data for cable-stayed bridge) after directly measuring or measuring, the original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on structure, and purpose is to determine according to these coordinate datas the geometric properties of Cable Structure.For cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
Set up the initial mechanical calculating benchmark model A of Cable Structure
o, set up initial Cable Structure bearing angular coordinate vector
u o, the Mechanics Calculation benchmark model A of the Cable Structure needed while setting up circulation beginning for the first time
1; The measured data of the Cable Structure in being completed according to Cable Structure, this measured data comprises the 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, and the Non-destructive Testing Data of rope etc. can be expressed the data of the health status of rope, according to design drawing and as-constructed drawing, utilize mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure
o; If there is no the measured data of the structure in Cable Structure completion, so just before setting up health monitoring systems, this Cable Structure is surveyed, obtain equally the measured data of Cable Structure, according to design drawing, the as-constructed drawing of these data and Cable Structure, utilize equally mechanics method to set up the initial mechanical calculating benchmark model A of Cable Structure
o; No matter which kind of method to obtain A by
o, based on A
othe Cable Structure computational data calculated must approach its measured data very much, and difference therebetween generally must not be greater than 5%; Corresponding to A
ocable Structure bearing angular data form initial Cable Structure bearing angular coordinate vector
u o; A
owith
u obe constant, only set up when circulation beginning for the first time; The Mechanics Calculation benchmark model of the Cable Structure of setting up when the i time circulation starts is designated as A
i, wherein i means cycle index; In application form of the present invention, alphabetical i is except the place that means significantly number of steps, and alphabetical i only means cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation starts for the first time is designated as A
1, A in the present invention
1just equal A
o; For sake of convenience, name " the current Mechanics Calculation benchmark model of Cable Structure A
ti o", A in circulation each time
ti ocan constantly update as required, when circulation starts each time, A
ti oequal A
i; Equally for sake of convenience, name " Cable Structure actual measurement bearing angular coordinate vector
u ti ", in circulation each time, constantly actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas form current cable structure actual measurement bearing angular coordinate vector
u ti , vector
u ti element with the vector
u othe angular coordinate of the equidirectional of the element representation same abutment of same position; For sake of convenience, for the i time circulation, the last time is upgraded to A
ti othe time Cable Structure bearing angular coordinate current data be designated as current cable structural bearings angular coordinate vector
u ti o; When circulation starts for the first time, A
t1 oequal A
1,
u t1 oequal
u o.A
1the health status of corresponding rope by
d 1 o describe; Mechanics Calculation benchmark model A
ithe health status of corresponding rope by
d i o describe.
The 5th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system is (for example, containing the measurement of angle subsystem, the cable force measurement subsystem, the strain measurement subsystem, volume coordinate is measured subsystem, signal conditioner etc.), Cable Structure bearing angular coordinate monitoring system is (containing angle measuring sensor, signal conditioner etc.), the cable force monitoring system is (for example, containing acceleration transducer, signal conditioner etc.), the horizontal range monitoring system of each support cable two supporting end points, signal (data) collector, computing machine and the panalarm of communicating by letter.The horizontal range of the Suo Li of each monitored amount, each support cable and each root support cable two supporting end points must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector by the signal monitored; Signal is delivered to computing machine through signal picker; Computing machine is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; When monitoring rope lax or damage arranged, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 6th step: establishment the cable system health monitoring systems software of pass line structure on supervisory control comuter.All move this software at circulation time each time, this software is all the time in operation in other words.This software will complete the functions (all work that can complete with computing machine in this specific implementation method) such as monitoring that the present invention's " progressive method of the identification slack line based on hybrid monitoring during angular displacement of support " required by task wants, record, control, storage, calculating, notice, warning, and can regularly or by personnel's operational health monitoring system generate cable system health condition form, can also for example, according to the condition of setting (damage reach a certain value), notice or prompting monitor staff notify specific technician to complete necessary evaluation work automatically.
The 7th step: step starts circulation running thus, is designated as for sake of convenience the i time circulation, i=1 wherein, and 2,3,4,5 ...Actual measurement (containing angle measuring sensor, signal conditioner etc.) obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas form current cable structure actual measurement bearing angular coordinate vector
u ti , according to current cable structure actual measurement bearing angular coordinate vector
u ti , upgrade where necessary the current Mechanics Calculation benchmark model of Cable Structure A
ti owith current cable structural bearings angular coordinate vector
u ti o.Concrete grammar is:
When circulation starts each time, make A
ti oequal A
i;
Actual measurement obtains current cable structure actual measurement bearing angular coordinate vector
u ti after, relatively
u ti with
u ti oif,
u ti equal
u ti o, do not need A
ti oupgraded;
Actual measurement obtains current cable structure actual measurement bearing angular coordinate vector
u ti after, relatively
u ti with
u ti oif,
u ti be not equal to
u ti o, need A
ti oupgraded, update method is: first calculate
u ti with
u opoor,
u ti with
u odifference be exactly that the current cable structural bearings is about setting up A
othe time the current angular displacement of support of Cable Structure bearing, with current angular displacement of support vector
vmean angular displacement of support, current angular displacement of support vector
vin element and the angular displacement of support component between be one-to-one relationship, current angular displacement of support vector
vin the numerical value of an element corresponding to the rotation of an assigned direction of an appointment bearing; Upgrade A
ti omethod be: at A
obasis on to make the health status of rope be cable system initial damage vector
d i o , more further to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support, the numerical value of current angular displacement of support constraint is just taken from current angular displacement of support vector
vthe numerical value of middle corresponding element, to A
oin the Cable Structure bearing apply the constraint of current angular displacement of support after, that finally obtain is exactly the current Mechanics Calculation benchmark model A of renewal
ti o, upgrade A
ti othe time,
u ti oall elements numerical value is also used
u ti all elements numerical value replaces, and has upgraded
u ti o, so just obtained correctly corresponding to A
ti o's
u ti o.
The 8th step: at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on carry out the several times Mechanics Calculation, obtain the monitored numerical quantity transformation matrices of Cable Structure virtual unit damage by calculating
Δ C i with nominal virtual unit damage vector
d i u .Concrete grammar is:
A. when the i time circulation starts, directly press step b to the monitored quantitative change matrix of the listed method acquisition Cable Structure virtual unit damage of steps d
Δ C i with nominal virtual unit damage vector
d i u ; The moment started in non-the i time circulation, in the 7th step to A
ti oafter being upgraded, directly press step b to the monitored quantitative change matrix of the listed method acquisition Cable Structure virtual unit damage of steps d
Δ C i with nominal virtual unit damage vector
d i u ; The moment started in non-the i time circulation, if in the 7th step not to A
ti oupgraded, directly proceed to herein the 9th step and carry out follow-up work;
B. at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on calculation times numerical value, have
nthe root rope just has
ninferior calculating, calculating each time in the hypothesis cable system only has a rope to increase virtual unit damage on the basis of original virtual lesion again, the rope that occurs virtual unit damage in calculating each time is different from the rope that occurs virtual unit damage in other calculating, and supposition each time has the virtual unit damage value of the rope of virtual unit damage can be different from the virtual unit damage value of other ropes, uses " nominal virtual unit damage vector
d i u " record the unit damage of the supposition of all ropes, calculate each time the current value of all monitored amounts, the current value of all monitored amounts that calculate each time forms one " the calculating current value vector of monitored amount ".When hypothesis the
jwhen the root rope has unit damage, available
c i tj mean corresponding " the current evaluation vector of monitored amount
c i tj ".While giving in this step each vectorial element numbering, should use same coding rule with other vector in the present invention, can guarantee like this any one element in each vector in this step, with in other vector, number identical element, expressed the relevant information of same monitored amount or same target;
C. that calculating each time " current evaluation vector of monitored amount
c i tj " deduct " the initial value vector of monitored amount
c i o " obtain a vector, then obtain one " the numerical value change vector of monitored amount " after the virtual unit damage value of supposition during each element of this vector is calculated divided by this; Have
nthe root rope just has
nindividual " the numerical value change vector of monitored amount ";
D. by this
nindividual " the numerical value change vector of monitored amount " forms and has successively
n" the monitored numerical quantity transformation matrices of virtual unit damage of row
Δ C i "; " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " each row corresponding to one " the numerical value change vector of monitored amount "; The coding rule of the row of " the monitored quantitative change matrix of virtual unit damage " and current nominal virtual lesion vector
d i c with current actual virtual lesion vector
d i the element coding rule identical.
While reaching in this step thereafter to each vectorial element numbering, should use same coding rule with other vector in the present invention, can guarantee like this any one element in each vector in this step, with in other vector, number identical element, 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 ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i "), when the 8th step is calculated each time, when in calculating each time, in the hypothesis cable system, only having on the basis of rope at original virtual lesion and increase virtual unit damage again, calculate each time and form " virtual lesion a vector
d i t ", the virtual lesion vector
d i t element number equal the quantity of rope, virtual lesion vector
d i t all elements in only have the numerical value of an element to get to calculate each time in hypothesis increase the virtual unit damage value of the rope of virtual unit damage,
d i t the numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the rope of virtual unit damage, is 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 bringing formula (13) into (should be noted, in formula (13)
c i with
c i tj bring into,
d i c with
d i t bring into), obtain a linear relationship error vector
e i , calculate each time a linear relationship error vector
e i ; Have
nthe root rope just has
ninferior calculating, just have
nindividual linear relationship error vector
e i , by this
nindividual linear relationship error vector
e i obtain a vector after addition, by 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 .By vector
g i be kept on the hard disc of computer of operation health monitoring systems software, for the health monitoring systems software application.
Will " initial rope force vector
f o ", " the initial value vector of monitored amount
c i o ", " nominal virtual unit damage vector
d i u ", " initial drift vector
l o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and the parameters such as unit weight of the elastic modulus of all ropes, initial cross sectional area, rope in the mode of data file, be kept on the hard disc of computer of operation health monitoring systems software.
The tenth step: actual measurement obtains the current cable power of all support cables of Cable Structure, forms the current cable force vector
f i ; Simultaneously, actual measurement obtains the current measured value of the monitored amount of all appointments of Cable Structure, forms " the current value vector of monitored amount
c i ".Actual measurement obtains the volume coordinate of two supporting end points of all support cables, and the volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges.
The 11 step: according to " current (calculating or actual measurement) numerical value vector of monitored amount
c i " " the initial value vector of monitored amount together
c i o ", " the monitored numerical quantity transformation matrices of virtual unit damage
Δ C i " and " current nominal virtual lesion vector
d i c " between the linear approximate relationship (seeing formula (9)) that exists, calculate the current nominal virtual lesion vector of cable system 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: the multiple-objection optimization based on genetic algorithm, the multiple-objection optimization based on artificial neural network, the multi-objective optimization algorithm based on population, the multiple-objection optimization based on ant group algorithm, leash law (Constrain Method), weighted method (Weighted Sum Method), Objective Programming (Goal Attainment Method) etc.Because various multi-objective optimization algorithms are all conventional algorithms, can realize easily, this implementation step only be take Objective Programming and is solved current nominal virtual lesion vector as example provides
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 Objective Programming, formula (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (29) and formula (30), in formula (29)
γ i 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 (29) is to find the real number of an absolute value minimum
γ i , make formula (30) be met.In formula (30)
g(d i c )by formula (31) definition, weighing vector in formula (30)
w i with
γ i product representation formula (30) in
g(d i c )with vector
g i between the deviation that allows,
g i definition referring to formula (15), its value will calculate in the 8th step.Vector during actual computation
w i can with vector
g i identical.The concrete programming of Objective Programming realizes having had universal program directly to adopt.Just can be in the hope of current name damage vector according to Objective Programming
d i c .
Try to achieve current nominal virtual lesion vector
d i c after
,the current actual virtual lesion vector that can obtain according to formula (17)
d i each element, current actual virtual lesion vector
d i exactly with reasonable error but can be more exactly from all ropes, determine the position of problematic rope (be virtual damaged cable, may be impaired may be also lax) and the solution of virtual lesion degree thereof.If the current actual virtual lesion vector solved
d i the numerical value of a certain element be 0, mean that the corresponding rope of this element is intact, not damage or lax; If its numerical value is 100%, mean that the corresponding rope of this element has completely lost load-bearing capacity; If its numerical value between 0 and 100%, means this rope and has lost the load-bearing capacity of corresponding proportion.
The 12 step: due to current actual virtual lesion vector
d i element numerical value represent the virtual lesion degree of corresponding rope, so it is impaired or relaxed and possible degree of injury or relax level just to define which Suo Keneng according to current actual virtual lesion vector, but damage has occurred actually or has occurred to relax in these ropes, need be differentiated.The method of differentiating is varied; can be by removing the protective seam of support cable; support cable is carried out to visual discriminating; perhaps by optical imaging apparatus, carry out visual discriminating; also can to whether support cable is impaired, be differentiated by lossless detection method, UT (Ultrasonic Testing) is exactly a kind of now widely used lossless detection method.After differentiating, those do not find that damage and virtual lesion degree are not that 0 support cable is exactly that lax rope has occurred, and need exactly to adjust the rope of Suo Li, can be completed by software the solving of relax level (being the long adjustment amount of rope) of slack line.
The 13 step: in this circulation,
itry to achieve current nominal virtual lesion vector in inferior circulation
d i c after, at first, according to formula (26), formula (27), set up mark vector
b i .
The 14 step: if mark vector
b i element be 0 entirely, get back to the 7th step and continue this circulation; If mark vector
b i element be not 0 entirely, enter next step, i.e. the 15 step.
The 15 step: according to formula (28) calculate next time,
i+ 1 required initial virtual lesion vector of circulation
d i+
1 o each element
d i+
1 oj .
The 16 step: at the current Mechanics Calculation benchmark model of Cable Structure A
ti obasis on, the health status that makes rope is the vector that previous step calculates
d i+
1 o after, obtain new Mechanics Calculation benchmark model, the Mechanics Calculation benchmark model A that (the i+1 time) circulates required next time
i+1.
The 17 step: by Mechanics Calculation benchmark model A
i+1calculate corresponding to model A
i+1the numerical value of all monitored amounts of structure, these numerical value form next time, required vector circulates for the i+1 time
c i+1 o , i.e. the initial value vector of monitored amount
.
The 18 step: set up next time, the required current Mechanics Calculation benchmark model of Cable Structure A circulates for the i+1 time
ti+1 o, get A
ti+1 oequal A
i+1.
The 19 step: set up next time, required current cable structural bearings angular coordinate vector circulates for the i+1 time
u ti+
1 o, get
u ti+
1 oequal
u ti o.
The 20 step: get back to the 7th step, start circulation next time.