CN101832872B - Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring - Google Patents

Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring Download PDF

Info

Publication number
CN101832872B
CN101832872B CN2010101375622A CN201010137562A CN101832872B CN 101832872 B CN101832872 B CN 101832872B CN 2010101375622 A CN2010101375622 A CN 2010101375622A CN 201010137562 A CN201010137562 A CN 201010137562A CN 101832872 B CN101832872 B CN 101832872B
Authority
CN
China
Prior art keywords
vector
evaluation
damage
cable
time
Prior art date
Application number
CN2010101375622A
Other languages
Chinese (zh)
Other versions
CN101832872A (en
Inventor
韩玉林
王芳
Original Assignee
东南大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 东南大学 filed Critical 东南大学
Priority to CN2010101375622A priority Critical patent/CN101832872B/en
Publication of CN101832872A publication Critical patent/CN101832872A/en
Application granted granted Critical
Publication of CN101832872B publication Critical patent/CN101832872B/en

Links

Abstract

The invention discloses a progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring, and provides a method for performing piecewise approximation on nonlinear relationship by using linear relationship to divide a large range into continuous small ranges, considering the fact that the linear relationship between the current value vector and the initial value vector of a monitored volume and linear relationship between the unit damage variation matrix of the monitored volume and the current notional health state vector are approximate, wherein linear relationship in each small range is accurate enough, and each small range can use proper algorithms such as multi-objective optimization algorithm and the like to rapidly identify support displacement, damaged cables and loose cables.

Description

Based on the identification damaged cable of space coordinate monitoring and the progressive method of support displacement
Technical field
The present invention is based on structural health monitoring technology, based on space coordinate monitoring, adopt progressive method to discern the damaged cable in the cable system of support displacement and Cable Structure, belong to the engineering structure health monitoring field.
Background technology
Support 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 support displacement and Cable Structure is a kind of method that has potentiality.When displacement appears in bearing or the health status of cable system when changing, or two kinds of situation when taking place simultaneously; Can cause the variation of the measurable parameter of structure; For example can cause the variation of Suo Li; Can influence the distortion or the strain of Cable Structure; Can influence the shape or the volume coordinate of Cable Structure; Can cause the variation (the for example arbitrarily variation of the angle coordinate of the straight line of any this point of mistake in any section of body structure surface, the perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily) of angle coordinate of any imaginary line of the every bit of Cable Structure, all these change the health status information that has all comprised cable system; In fact the variation of these measurable parameters comprised cable system health status information, comprised support displacement information, that is to say that the measurable parameter that can utilize structure discerns support displacement and damaged cable.
In order reliable monitoring and judgement to be arranged to the health status and the support 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 support 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 support displacement assessment and cable system.
Summary of the invention
Technical matters: the invention discloses a kind of based on the identification damaged cable of space coordinate monitoring and the progressive method of support displacement, this method be based on space coordinate monitoring, adopt health monitor method progressive method, that can discern support displacement and damaged cable rationally and effectively.
Technical scheme: structures such as cable-stayed bridge, suspension bridge, truss-frame structure have a common ground; Be exactly that they have many parts that bear tensile load; Like suspension cable, main push-towing rope, hoist cable, pull bar or the like; The common ground of this class formation is to be support unit with rope, cable or the rod member that only bears tensile load, and the present invention is " Cable Structure " with such structure representation for simplicity.In the military service process of Cable Structure; The supporting system of Cable Structure (refers to all ropeway carrying-ropes, reaches all rod members that only bear tensile load that play supporting role; For simplicity, this patent is called " cable system " with whole support unit unifications of this class formation, but in fact cable system not only refers to support rope; Also comprise the rod member that only bears tensile load) can be impaired; The bearing of Cable Structure also displacement possibly occur simultaneously, and these safety that change to Cable Structure are a kind of threats, and establishing the quantity of rope and the quantity sum of support displacement component is N.For narrate convenient for the purpose of, the unified rope of being assessed of claiming of the present invention is " by evaluation object " with support displacement, gives by the evaluation object serial number, the present invention is with representing this numbering with variable j, j=1,2,3 ..., N, so we can say that N is individual by evaluation object.
The present invention is made up of the two large divisions.Be respectively: one, set up by the method for required knowledge base of evaluation object health monitoring systems and parameter, based on knowledge base (containing parameter) and the volume coordinate of surveying Cable Structure by evaluation object health status appraisal procedure; Two, the software and hardware part of health monitoring systems.
First of the present invention: foundation is used for by the method for the knowledge base of evaluation object health monitoring and parameter.Can be successively circularly as follows, laddering carrying out:
The first step: when beginning circulation each time, in the time of at first need setting up or set up this circulation beginning by evaluation object initial health vector d o i(i=1,2,3 ...), set up the initial Mechanics Calculation benchmark model A of Cable Structure o(for example finite element benchmark model, A in the present invention 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 alphabetical in the present invention i only representes cycle index, i.e. the i time circulation.
The Cable Structure that needs during the i time circulation beginning " initial health vector d o i" (shown in (1)), use d o iCable Structure is (with Mechanics Calculation benchmark model A when representing the i time circulation beginning iThe initial health of Cable Structure expression).
d o i = d o 1 i d o 2 i · · · d oj i · · · d oN i T - - - ( 1 )
D in the formula (1) i Oj(i=1,2,3, J=1,2,3 ...., when N) the i time circulation of expression begins, Mechanics Calculation benchmark model A iIn the j of cable system by the current health status of evaluation object, if should be the rope (or pull bar) in the cable system, d so by evaluation object iRepresent its current damage, d iBeing to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be displacement component, a d so of a bearing by evaluation object iRepresent its present bit shift value.T representes the transposition (back together) of vector in the formula (1).
Set up the initial health vector during circulation beginning for the first time and (be designated as d according to formula (1) 1 o) time, the data and the support displacement that utilize the Non-Destructive Testing data etc. of rope can express the health status of rope are measured and are set up by evaluation object initial health vector d 1 oIf when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vectorial d 1 oIn get 0 with each element numerical value of Suo Xiangguan.
The i time (i=2,3,4,5,6 ...) when beginning circulation need by evaluation object initial health vector d 1 o, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
The Mechanics Calculation benchmark model that need set up during the i time circulation beginning or the Mechanics Calculation benchmark model of having set up are designated as A i
Measured data according to the Cable Structure in the Cable Structure completion (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope) and design drawing, as-constructed drawing, utilize mechanics method (for example finite element method) to set up A oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, structure is surveyed; The measured data that obtains Cable Structure (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope); According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method (for example finite element method) to set up A 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 1Corresponding by the health status of evaluation object by d 1 oDescribe.
The i time (i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins i, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Existing Mechanics Calculation benchmark model A 1With by evaluation object initial health vector d 1 oAfter, model A 1In each by the health status of evaluation object by vectorial d 1 oExpress.At A 1The basis on, all are changed to 0 by the health status numerical value of evaluation object, mechanical model A 1Be updated to one all be that 0 mechanical model (is designated as A by the health status of evaluation object 0), mechanical model A 0Be actually the corresponding mechanical model of Cable Structure of excellent no support displacement.Might as well claim model A 0For the not damaged of Cable Structure does not have support displacement model A 0
" the whole monitored spatial data of structure " described by the volume coordinate specified point of K on the structure, that reach L assigned direction of each specified point, and the variation of structure space coordinate data is exactly the variation of all volume coordinate components of K specified point.(individual volume coordinate measured value of M=K * L) or calculated value characterize the structure space coordinate information to each total M.K and M must not be less than by the quantity N of evaluation object.
For simplicity, in the present invention " the monitored spatial data of structure " abbreviated as " monitored amount ".When mentioning " so-and-so matrix of monitored amount or so-and-so vector " in the back, also can be read as " volume coordinate so-and-so matrix or so-and-so vector ".
Among the present invention with monitored amount initial value vector C i o" (i=1,2,3 ...) initial value (referring to formula (2)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation of expression begins, C i oFull name be " the initial value vector of the i time monitored amount of circulation ".
C o i = C o 1 i C o 2 i · · · C ok i · · · C oM i T - - - ( 2 )
C in the formula (2) i Ok(i=1,2,3, K=1,2,3 ...., M; M>=N; ) k monitored amount when being the i time circulation beginning, in the Cable Structure.Vector C i oBe that M monitored amount by front definition forms according to certain series arrangement, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
During circulation beginning for the first time, " the initial value vector C of the 1st the monitored amount that circulates 1 o" (seeing formula (2)) be made up of measured data, because according to model A 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 ...) when beginning circulation " the vectorial C of the initial value of the i time monitored amount of circulation that needs i o", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) calculate before the loop ends and obtain, concrete grammar is in the back literary composition narration.
Second step: circulation each time needs to set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector ", and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C i, " nominal unit damage vector " that the i time circulation set up is designated as D i u, i=1,2,3 ...
The Cable Structure " unit damage monitored numerical quantity transformation matrices " that circulation is for the first time set up is designated as Δ C 1Set up Δ C 1Process following:
Mechanics Calculation benchmark model A in Cable Structure 1The basis on carry out several times and calculate, equal N on the calculation times numerical value.Calculating hypothesis each time has only one by evaluation object unit damage to be arranged; Concrete, if should be a supporting rope in the cable system by evaluation object, so just this supporting rope of hypothesis has unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage); If should be the displacement component of a direction of a bearing by evaluation object; Just suppose that this bearing is this sense of displacement generation unit displacement (for example 10mm, 20mm, 30mm etc. are unit displacement).For narrating conveniently, the damage and the support displacement of the supporting rope that the present invention will suppose are referred to as unit damage.Calculate for convenient; When setting unit damage in the circulation each time can all be structural health conditions during this time circulation beginning as being healthy fully, and set on this basis unit damage (in subsequent step, calculate, by the health status numerical value of evaluation object---be called nominal health status vector d i c(i=1; 2,3 ...); All with respect to this time when beginning circulation, with the health status of Cable Structure as being healthy fully speech, therefore must foundation after the formula that provides of the literary composition nominal health status numerical value that will calculate be converted into true health status numerical value).With a round-robin occur in calculating each time unit damage be different from by evaluation object other time occur in calculating unit damage by evaluation object; And supposition has the unit damage value by evaluation object of unit damage can be different from other by the unit damage value of evaluation object each time, with " nominal unit damage vector D i u" (shown in (3)) write down that all are by the unit damage of the supposition of evaluation object in each time circulation, circulation time is designated as D for the first time 1 uCalculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, the current calculated value of the M of appointment monitored amount in front; The current calculated value that calculates gained M monitored amount is each time formed one " the current numerical value vector of the calculating of monitored amount ", and (when j of hypothesis had unit damage by evaluation object, available formula (4) was represented the current numerical value vector of the calculating C of M monitored amount of all appointments 1 Tj); The current numerical value vector of the calculating of the monitored amount that calculates each time deducts the initial value vector C of monitored amount 1 o, the gained vector be exactly under this condition (with have unit damage by the mark that is numbered of evaluation object) " the numerical value change vector of monitored amount " (when j is had unit damage by evaluation object, use δ C 1 jThe numerical value change vector of representing monitored amount, δ C 1 jDefinition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vectorial D for formula (4) 1 uJ element D UjGained), the numerical value change of monitored amount vector δ C 1 jEach element representation because that supposition has a unit damage when calculating had unit damage (D for example by evaluation object (for example j by evaluation object) Uj), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage D of supposition UjRate of change; There is N N " the numerical value change vector of monitored amount " just to be arranged by evaluation object; The numerical value change vector of each monitored amount has M (general; The individual element of M>=N) is formed " the unit damage monitored numerical quantity transformation matrices Δ C that M * N element arranged successively by this N " the numerical value change vector of monitored amount " 1" (the capable N row of M), each vectorial δ C 1 j(j=1,2,3 ...., N) be matrix Δ C 1One row, Δ C 1Definition suc as formula shown in (8).
C u i = C u 1 i C u 2 i · · · C uj i · · · C uN i T - - - ( 3 )
Nominal unit damage vector D in the formula (3) i uElement D i Uj(i=1,2,3, J=1,2,3 ...., N) j of supposition is individual by the unit damage numerical value of evaluation object in the i time circulation of expression, vectorial D i uIn the numerical value of each element can be the same or different.
C tj i = C tk 1 i C tk 2 i · · · C tjk i · · · C tjM i T - - - ( 4 )
Elements C in the formula (4) i Tjk(i=1,2,3 ...; J=1,2,3 ...., N; K=1,2,3 ...., M; The i time circulation of the expression of M>=N) be because j is individual when by evaluation object unit damage being arranged, according to the current numerical value of calculating of the monitored amount of pairing k the appointment of coding rule.
δ C j i = C tj i - C 0 i D uj i - - - ( 5 )
The subscript i of each amount in the formula (5) (i=1,2,3 ...) the i time circulation of expression, subscript j (j=1,2,3 ...., N) j of expression had unit damage by evaluation object, D in the formula i UjBe vectorial D i uIn j element.Vector δ C i jDefinition suc as formula shown in (6), δ C i jK (k=1,2,3 ...., M; The individual element δ C of M>=N) i JkRepresent to set up matrix Δ C in the i time circulation iThe time, suppose the change amount of calculating a gained k monitored amount when j has unit damage by evaluation object unit damage D with respect to supposition i UjRate of change, it defines suc as formula shown in (7).
δC j i = δC j 1 i δC j 2 i · · · δC jk i · · · δC jM i T - - - ( 6 )
δC jk i = C tjk i - C ok i D uj i - - - ( 7 )
The definition of each amount was narrated in front in the formula (7).
ΔC i = δC 1 i δ C 2 i · · · δC j i · · · δ C N i - - - ( 8 )
Vectorial δ C in the formula (8) i j(i=1,2,3 ....,, j=1,2,3 ...., N) in the i time circulation of expression, by evaluation object unit damage D is arranged because j is individual i UjCause, the relative value of all monitored amounts changes.Matrix Δ C iCoding rule and the front vector d of row (subscript j) i oThe coding rule of subscript j of element identical.
The 3rd step: identification is by the current health status of evaluation object.Detailed process is following.
I (i=1,2,3 ...) in the inferior circulation, " current (calculating or actual measurement) numerical value vector C of monitored amount i" " the initial value of monitored amount vector C together i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between linear approximate relationship, shown in (9) or formula (10).
C i = C o i + Δ C i · d c i - - - ( 9 )
C i - C o i = Δ C i · d c i - - - ( 10 )
Current (calculating or actual measurement) numerical value vector C of monitored amount in formula (9) and the formula (10) iDefinition be similar to the initial value vector C of monitored amount i oDefinition, see formula (11); By the vectorial d of the current name damage of evaluation object i cDefinition see formula (12).
C i = C 1 i C 2 i · · · C k i · · · C M i T - - - ( 11 )
Elements C in the formula (11) i k(i=1,2,3 ....; K=1,2,3 ...., M; M>=N) be the i time circulation time Cable Structure, according to the current numerical value of the monitored amount of the pairing k of being numbered of coding rule.
d c i = d c 1 i d c 2 i · · · d cj i · · · d cN i T - - - ( 12 )
D in the formula (12) i Cj(i=1,2,3 ....; J=1,2,3 ...., be that Cable Structure j is individual by the current nominal impairment value of evaluation object in the i time circulation N), vectorial d i cThe coding rule and the matrix Δ C of subscript j of element iThe coding rule of row identical.
When by evaluation object actual damage or support displacement when not too big; Because the Cable Structure material still is in the linear elasticity stage; The distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (9) or formula (10) is less with the error of actual conditions, and error can be used error vector e i(formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).
e i = abs ( Δ C i · d c i - C i + C o i ) - - - ( 13 )
Abs () is the function that takes absolute value in the formula (13), and each element of the vector of trying to achieve in the bracket is taken absolute value.
Because there are certain error in formula (9) or the represented linear relationship of formula (10), therefore can not be simply according to formula (9) or formula (10) and " current (actual measurement) numerical value vector C of monitored amount i" come directly to find the solution to obtain the vectorial d of current name damage i cDamage vectorial d and obtain current name i cAcceptable separating (promptly have reasonable error, but can be more accurately from cable system, confirm damaged cable position and degree of injury thereof, confirm the support displacement amount) become a rational solution, available formula (14) is expressed this method.
abs ( ΔC i · d c i - C i + C o i ) ≤ g i - - - ( 14 )
Abs () is the function that takes absolute value in the formula (14), vectorial g iDescription departs from the reasonable deviation of ideal linearity relation (formula (9) or formula (10)), is defined by formula (15).
g i = g 1 i g 2 i · · · g k i · · · g M i T - - - ( 15 )
G in the formula (15) i k(i=1,2,3 ....; K=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (9) or the formula (10) in the i time circulation has been described.Vector g iCan be according to the error vector e of formula (13) definition iTentative calculation is selected.
Initial value vector C in monitored amount i o(survey or calculate), unit damage monitored numerical quantity transformation matrices Δ C iThe current numerical value vector C of (calculating) and monitored amount iWhen (actual measurement obtains) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain the vectorial d of current name damage i cAcceptable separating, current actual damage vector d iThe element of (formula (16) is seen in definition) can calculate according to formula (17), has just obtained by the current actual damage vector of evaluation object d iThereby, can be by d iConfirm damaged cable position and degree of injury, confirm the support displacement amount, realized that just damage identification and support displacement discern.
d i = d 1 i d 2 i · · · d j i · · · d N i T - - - ( 16 )
D in the formula (16) i j(i=1,2,3, J=1,2,3 ...., N) j is by the actual damage value of evaluation object in the i time circulation of expression, and formula (17) is seen in its definition, if this was the rope (or pull bar) in the cable system, d so by evaluation object i jRepresent its current damage, d i jBeing to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%, if should be displacement component, a d so of a bearing by evaluation object i jRepresent its present bit shift value.Vector d iCoding rule and the formula (1) of element in vectorial d i oThe coding rule of element identical.
d j i = 1 - ( 1 - d oj i ) ( 1 - d cj i ) - - - ( 17 )
D in the formula (17) i Oj(i=1,2,3,4, J=1,2,3 ...., N) be vectorial d i oJ element, d i CjBe vectorial d i cJ element.
The 4th step: judge whether to finish this (the i time) circulation, if, then accomplish the preceding tailing in work of this loop ends, for next time (promptly the i+1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following.
In this (the i time) circulation, try to achieve the vectorial d of current name damage i cAfter, at first, set up mark vector F according to formula (18) i, formula (19) has provided mark vector F iThe definition of j element; If mark vector F iElement be 0 entirely, then in this circulation, continue health monitoring and calculating to Cable Structure; If mark vector F iElement be not 0 entirely, then accomplish subsequent step after, get into circulation next time.So-called subsequent step is: at first, according to formula (20) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation I+1 oEach element d I+1 OjThe 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, order is d by the health status situation of evaluation object I+1 oThe back upgrade and to obtain next time (the i+1 time, i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation I+1At last, through to Mechanics Calculation benchmark model A I+1The initial value that calculates monitored amount, by its form next time (promptly the i+1 time, i=1,2,3,4 ...) required " the initial value vector C of monitored amount of circulation I+1 o" (i=1,2,3,4 ...).
F i = F 1 i F 2 i · · · F j i · · · F N i T - - - ( 18 )
Mark vector F in the formula (18) iSubscript i represent the i time the circulation, its element F i j(j=1,2,3 ..., subscript j N) representes that j can only be got 0 and 1 two amount by the damage characteristic of evaluation object, concrete value rule is seen formula (19).
F j i = 0 , if d cj i < D uj i 1 , if d cj i &GreaterEqual; D uj i - - - ( 19 )
Element F in the formula (19) i jBe mark vector F iJ element, D i UjBe nominal unit damage vector D i uJ element (seeing formula (3)), d i CjBe the vectorial d of current name damage i cJ element (seeing formula (12)), they represent that all j is by the relevant information of evaluation object.
d oj i + 1 = 1 - ( 1 - d oj i ) ( 1 - D uj i F j i ) - - - ( 20 )
D in the formula (20) i UjBe nominal unit damage vector D i uJ element (seeing formula (3)), d i CjBe the vectorial d of current name damage i cJ element (seeing formula (12)).
Second portion of the present invention: the software and hardware part of health monitoring systems.
Hardware components comprises monitored amount monitoring system, signal picker and computing machine etc.Require to monitor in real time or quasi real time each monitored amount.
Software should be used following function by tool: software section should be accomplished the process that first of the present invention sets, and promptly accomplishes functions such as needed among the present invention, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.
The inventive method specifically comprises:
A. for for the purpose of narration is convenient, the present invention is unified to claim that the supporting rope of being assessed is that establishing the quantity of the supporting rope of being assessed and the quantity sum of support displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with the support displacement component; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention representes this numbering with variable j, j=1, and 2,3 ..., N;
B. confirm the measured point with monitored volume coordinate of appointment, give all specified point numberings; Confirmed each measurement point with monitored volume coordinate component, give all measured volume coordinate components numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored spatial data of structure " is made up of above-mentioned all measured volume coordinate components; For simplicity, in the present invention " the monitored spatial data of structure " is called " monitored amount "; The quantity sum of all measured volume coordinate components must not be less than N;
C. utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector d by the Non-Destructive Testing data of evaluation object etc. i oIf during not by the Non-Destructive Testing data of evaluation object, vectorial d i oEach element numerical value get 0; Vector d i oThe coding rule of element with identical by the coding rule of evaluation object; The present invention representes cycle index with i, i=1, and 2,3 ...; Here be circulation for the first time, i gets 1, the initial health vector d that promptly sets up here i oCan be embodied as d 1 o
D. setting up initial health vector d 1 oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount i oHere be circulation for the first time, i gets 1, the initial value vector C of the monitored amount of promptly setting up here i oCan be embodied as C 1 oObtain monitored amount initial value vector C in actual measurement 1 oThe time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure support coordinate data of structure;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure support coordinate data are set up Cable Structure iHere be circulation for the first time, i gets 1, the Mechanics Calculation benchmark model A of the Cable Structure of promptly 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, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C i" and " nominal unit damage vector D i u";
G. actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount i"; When numbering to the element of the institute's directed quantity that occurred before this step and this step; Should use same coding rule; Can guarantee element each vector, that numbering is identical of occurring before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
H. define the vectorial d of current name damage i cWith current actual damage vector d i, the element number of two damage vectors equals by the quantity of evaluation object, the vectorial d of current name damage i cElement numerical value represent corresponding to the current nominal degree of injury or the support displacement of evaluation object, current actual damage vector d iElement numerical value represent corresponding by the current actual damage degree or the support displacement of evaluation object; The element number of the element of two damage vectors equals by the quantity of evaluation object; The element of two damage vectors and be one-to-one relationship between the evaluation object, the coding rules of the element of two damage vectors are with identical by the coding rule of evaluation object;
I. according to " the current numerical value vector C of monitored amount i" " the initial value of monitored amount vector C together i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + &Delta; C i &CenterDot; d c i Formula 1
J. the current actual damage vector d that utilizes formula 2 to express iWith initial damage vector d i oWith the vectorial d of current name damage i cElement between relation, calculate current actual damage vector d iAll elements;
d j i = 1 - ( 1 - d Oj i ) ( 1 - d Cj i ) Formula 2
J=1 in the formula 2,2,3 ..., N;
Current actual damage vector d iElement numerical value represent corresponding to the actual damage degree of evaluation object or actual support displacement, according to current actual damage vector d iJust can define the impaired and degree of injury of which rope, just can confirm actual support displacement; If a certain element of current actual damage vector is corresponding to being a rope in the cable system; And its numerical value is 0, representes that the pairing rope of this element is intact, does not damage; If its numerical value is 100%; Represent that then the pairing rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this rope; If a certain element of current actual damage vector is corresponding to displacement component, a d so of a bearing i jRepresent its present bit shift value;
K. try to achieve the vectorial d of current name damage i cAfter, set up mark vector F according to formula 3 i, formula 4 has provided mark vector F iThe definition of j element;
F i = F 1 i F 2 i &CenterDot; &CenterDot; &CenterDot; F j i &CenterDot; &CenterDot; &CenterDot; F N i T Formula 3
F j i = 0 , If d Cj i < D Uj i 1 , If d Cj i &GreaterEqual; D Uj i Formula 4
Element F in the formula 4 i jBe mark vector F iJ element, D i UjBe nominal unit damage vector D i uJ element, d i CjBe the vectorial d of current name damage i cJ element, they all represent j by the relevant information of evaluation object, j=1 in the formula 4,2,3 ..., N;
If mark vector F l. iElement be 0 entirely, then get back to g step and continue this circulation; If mark vector F iElement be not 0 entirely, then get into next step, i.e. m step;
M. the initial damage vector d that calculates next time, promptly circulates required the i+1 time according to formula 5 I+1 oEach element d I+1 Oj
d Oj i + 1 = 1 - ( 1 - d Oj i ) ( 1 - D Uj i F j i ) Formula 5
D in the formula 5 i UjBe nominal unit damage vector D i uJ element, d i CjBe the vectorial d of current name damage i cJ element, F i jBe mark vector F iJ element, j=1 in the formula 5,2,3 ..., N; Vector d I+1 oThe coding rule of element with identical by the coding rule of evaluation object;
N. at Mechanics Calculation benchmark model A iThe basis on, order is d by the health status of evaluation object I+1 oThe back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time I+1
O. pass through Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1Structure all monitored strains point, with the monitored strain numerical value that should change direction, these numerical value are formed next time, the vectorial C of initial value of the required monitored amount that promptly circulates for the i+1 time I+1 o
P. get back to the f step, beginning is circulation next time.
In step f, at Mechanics Calculation benchmark model A iThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C i" and " nominal unit damage vector D i u" concrete grammar be:
F1. at the Mechanics Calculation benchmark model A of Cable Structure iThe basis on carry out the several times Mechanics Calculation, equal N on the calculation times numerical value; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or displacement, to be increased unit damage or unit displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases unit displacement again at this sense of displacement; Increase again in calculating each time unit damage or unit displacement be different from by evaluation object other time increase again in calculating unit damage or unit displacement by evaluation object, with " nominal unit damage vector D i u" unit damage or the unit displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector; Obtain a monitored quantitative changeization vector, have N evaluation object that N monitored quantitative changeization vector just arranged;
F3. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the monitored amount unit change of the Cable Structure matrix Δ C that the N row are arranged successively i
Beneficial effect: method disclosed by the invention can identify the health status (position and the degree of injury that comprise all damaged cables) of Cable Structure support displacement and cable system simultaneously, and system and method disclosed by the invention is very useful to the safety of Cable Structure.
Embodiment
To the health monitoring of Cable Structure, the invention discloses a kind of system and method that can monitor in the Cable Structure health status of each root rope and each support displacement component in the cable system rationally and effectively simultaneously.The following explanation of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
The present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status of the cable system of Cable Structure.During practical implementation, the following step is a kind of in the various steps that can take.
The first step: for the purpose of narrating conveniently, the present invention is unified to claim that supporting rope and the support displacement component assessed are that establishing the quantity of the supporting rope of being assessed and the quantity sum of support displacement component is N, is N by the quantity of evaluation object promptly by evaluation object; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention representes this numbering with variable j, j=1, and 2,3 ..., N.
Confirm that (promptly all characterize the specified point of structure space coordinates for the measured point of appointment; Be provided with K specified point); Each specified point can be exactly a near point the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) of each root rope; This specified point can also be a near point the structural bearings, perhaps directly is exactly the structural bearings fulcrum, gives all specified point numberings; Confirm the measured volume coordinate component (established each measurement point L measured volume coordinate component arranged) of each measurement point, gave the measured volume coordinate component numbering of all appointments.Above-mentioned numbering will be used to generate the vector sum matrix equally in subsequent step." the whole monitored spatial data of structure " described by L volume coordinate component K specified point, that cross each specified point on top definite structure, and the variation of structure space coordinate is exactly variations all specified points, all space specified coordinate components.(individual volume coordinate component measurement value of M=K * L) or calculated value characterize the volume coordinate information of structure to each total M.K and M must not be less than N.For simplicity, in the present invention " the monitored spatial data of structure " is called " monitored amount ".
Second step: utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector d by the Non-Destructive Testing data of evaluation object etc. 1 oIf during not by the Non-Destructive Testing data of evaluation object, vectorial d 1 oEach element numerical value get 0; Vector d 1 oThe coding rule of element with identical by the coding rule of evaluation object.
The 3rd step: at initial health vector d 1 oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount 1 o
The 4th step: the initial value vector C that obtains monitored amount in actual measurement 1 oThe time, can adopt ripe measuring method to carry out cable force measurement, strain measurement, measurement of angle and volume coordinate and measure.Simultaneously; Calculate the initial Suo Li and the Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) of all ropes of Cable Structure after directly measuring or measuring; The original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is just can confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
According to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure (data such as initial Suo Li, structural modal data that comprise structure original geometric form data, strain data, all ropes, to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, strain data, rope force data, bridge), the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure support coordinate data are set up Cable Structure o, based on Mechanics Calculation benchmark model A 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 representes cycle index; Alphabetical i is except the place of representing number of steps significantly in the application form of the present invention, and alphabetical i only representes cycle index, i.e. the i time circulation; The Mechanics Calculation benchmark model of the Cable Structure of setting up when therefore circulation begins for the first time is designated as A 1, A among the present invention 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 spatial coordinate measuring system, signal conditioner etc.), signal (data) collector, the computing machine and the panalarm of communicating by letter.Each monitored amount all must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector with the signal that monitors; Signal is delivered to computing machine through signal picker; Computing machine then is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; When monitoring when being changed by the health status of evaluation object, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 6th step: establishment and the health monitoring systems software of pass line structure on supervisory control comuter.All move this software at circulation time each time, this software is all the time in operation in other words.This software will be accomplished functions such as the needed monitoring of each task of the present invention, record, control, storage, calculating, notice, warning (being all work that can accomplish with computing machine in this practical implementation method); And can regularly or by the personnel operation health monitoring systems generate Cable Structure health condition form; Can also be according to the condition of setting (for example damage reach a certain value), notice or prompting monitor staff notify specific technician to accomplish necessary evaluation work automatically.
The 7th step: the step begins circulation running thus, for narration conveniently is designated as the i time circulation, and i=1 wherein, 2,3,4,5 ....
The 8th step: the Mechanics Calculation benchmark model in Cable Structure is designated as A iThe basis on carry out the several times Mechanics Calculation, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C iWith nominal unit damage vector D i uConcrete 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 N on the calculation times numerical value; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one to be increased on the basis of original damage or displacement by evaluation object again unit damage or unit displacement are arranged; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases unit displacement again at this sense of displacement; Increase again in calculating each time unit damage or unit displacement be different from by evaluation object other time increase again in calculating unit damage or unit displacement by evaluation object, with " nominal unit damage vector D i u" unit damage or the unit displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed same monitored amount or same by the relevant information of evaluation object object.
B. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector, obtains a monitored quantitative changeization vector δ C i jThere is N N monitored quantitative changeization vector δ C just to be arranged by evaluation object i j(j=1,2,3 ..., N).
C. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the monitored amount unit change of the Cable Structure matrix Δ C that the N row are arranged successively i" unit damage monitored quantitative change matrix Δ C i" the vectorial d of current name damage of coding rule and back definition of row i cWith current actual damage vector d iThe element coding rule identical.
In this step, reach when giving each vectorial element numbering thereafter; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.
The 9th step: set up linear relationship error vector e iWith vectorial g iUtilize data (" the initial value vector C of monitored amount of front i o", " unit damage monitored quantitative change matrix Δ C i"); when the 8th step calculated each time; promptly in calculating each time, have only in the hypothesis cable system one by evaluation object on the basis of original damage or displacement, increase again unit damage or unit displacement are arranged in, calculate each time and form a health status vector d i t, health status vector d i tElement number equal by the quantity of evaluation object vectorial d i tAll elements in have only the numerical value of an element to get to calculate each time in hypothesis increase unit damage value or the unit displacement value of increase of the rope of unit damage, d i tThe numerical value of other element get 0, what that was not that numbering and the supposition of 0 element increase unit damage or unit displacement is identical by the corresponding relation of evaluation object, with the element of the same numberings of other vectors with the corresponding relation of this rope; With C i Tj, C i o, Δ C i, d i tBring formula (13) into, formula (13) d i cUse d i tBring into, obtain a linear relationship error vector e i, calculate a linear relationship error vector e each time iThere is N N calculating just to be arranged, N linear relationship error vector e just arranged by evaluation object i, with this N linear relationship error vector e iObtain a vector after the addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e iVector g iEqual final error vector e iWith vectorial g iBe kept on the hard disc of computer of operation health monitoring systems software, supply health monitoring systems software to use.Will " the initial value vector C of monitored amount i o" and " unit damage monitored quantitative change matrix Δ C i" etc. parameter be kept on the hard disc of computer of operation health monitoring systems software with the mode of data file.
The tenth step: actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount i".
The 11 step: according to " the current numerical value vector C of monitored amount i" " the initial value of monitored amount vector C together i o", " unit damage monitored quantitative change matrix Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship (formula (9)) that exists, calculate the vectorial d of the current name damage of cable system according to multi-objective optimization algorithm i cNoninferior solution.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted Sum Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution the vectorial d of current name damage with the goal programming method only i cProcess, the concrete implementation procedure of other algorithm can realize according to the requirement of its specific algorithm in a similar fashion.
According to the goal programming method, formula (9) can transform the multi-objective optimization question shown in an accepted way of doing sth (21) and the formula (22), γ in the formula (21) iBe a real number, R is a real number field, and area of space Ω has limited vectorial d i cSpan (the present embodiment requirements vector d of each element i cEach element be not less than 0, be not more than 1).The meaning of formula (21) is to seek the minimum real number γ of an absolute value i, make formula (22) be met.G (d in the formula (22) i c) by formula (23) definition, weighing vector W in the formula (22) iWith γ iProduct representation formula (22) in G (d i c) and vectorial g iBetween the deviation that allows, g iDefinition referring to formula (15), its value will the 8th the step calculate.Vector W during actual computation iCan with vectorial g iIdentical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of the vectorial d of current name damage i c
minimize?γ i
(21)
γ i∈R, d c i &Element; &Omega;
G ( d c i ) - W i &gamma; i &le; g i - - - ( 22 )
G ( d c i ) = ads ( &Delta; C i &CenterDot; d c i - C i + C o i ) - - - ( 23 )
Try to achieve the vectorial d of current name damage i cAfter, can be according to the vectorial d of the current actual damage that formula (17) obtain iEach element, current actual damage vector d iHave reasonable error exactly but can confirm more exactly damaged cable the separating of position and degree of injury thereof, can confirm separating of all support displacements more exactly.d iEach element corresponding to one by the health status of evaluation object; If should be the rope (or pull bar) in the cable system by evaluation object; Its current damage of the numeric representation of this element so; The numerical value of this element is to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%; If being somebody's turn to do by evaluation object be a displacement component of a bearing, so its present bit shift value of the numeric representation of this element.
The 12 step:, promptly try to achieve the vectorial d of current name damage in the i time circulation in this circulation i cAfter, set up mark vector F according to formula (18), formula (19) iIf mark vector F iElement be 0 entirely, then got back to for the tenth step and continue this circulation; If mark vector F iElement be not 0 entirely, then get into next step, i.e. the 13 step.
The 13 step: according to the initial damage vector d that formula (20) calculates next time, i.e. the i+1 time circulation is required I+1 oEach element d I+1 Oj
The 14 step: at Cable Structure Mechanics Calculation benchmark model A iThe basis on, order is the last vectorial d that calculates of step by the health status of evaluation object I+1 oAfter, obtain new Mechanics Calculation benchmark model, next time promptly the required Mechanics Calculation benchmark model A of (the i+1 time) circulation I+1
The 15 step: through to Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, required vectorial C promptly circulates for the i+1 time I+1 o, promptly the initial value of monitored amount is vectorial.
The 16 step: the computing machine in the health monitoring systems regularly generates cable system health condition form automatically or by the personnel operation health monitoring systems.
The 17 step: under specified requirements, the automatic operation communication panalarm of the computing machine in the health monitoring systems to monitor staff, owner and (or) personnel of appointment report to the police.
The 18 step: got back to for the 7th step, beginning is circulation next time.

Claims (2)

1. one kind based on the identification damaged cable of space coordinate monitoring and the progressive method of support displacement, it is characterized in that said method comprises:
A. claim that supporting rope and the support displacement component assessed are that establishing the quantity of the supporting rope of being assessed and the quantity sum of support displacement component is N, is N by the quantity of evaluation object promptly by evaluation object; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; This method is represented this numbering with variable j, j=1, and 2,3 ..., N;
B. confirm the measured point with monitored volume coordinate of appointment, give all specified point numberings; Confirmed each measurement point with monitored volume coordinate component, give all measured volume coordinate components numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored spatial data of structure " is made up of above-mentioned all measured volume coordinate components; For simplicity, in the method " the whole monitored spatial data of structure " is called " monitored amount "; The quantity sum of all measured volume coordinate components must not be less than N;
C. utilize by the Non-Destructive Testing data of evaluation object and can express by the foundation of the data of the health status of evaluation object by evaluation object initial health vector d i oIf during not by the Non-Destructive Testing data of evaluation object, vectorial d i oEach element numerical value get 0; Vector d i oThe coding rule of element with identical by the coding rule of evaluation object; This method is represented cycle index with i, i=1, and 2,3 ...; Here be circulation for the first time, i gets 1, the initial health vector d that promptly sets up here i oBe embodied as d 1 o
D. setting up initial health vector d 1 oThe time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector C of monitored amount i oHere be circulation for the first time, i gets 1, the initial value vector C of the monitored amount of promptly setting up here i oBe embodied as C 1 oObtain monitored amount initial value vector C in actual measurement 1 oThe time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure support coordinate data of structure;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure support coordinate data are set up Cable Structure iHere be circulation for the first time, i gets 1, the Mechanics Calculation benchmark model A of the Cable Structure of promptly setting up here iBe embodied as A 1
F. at Mechanics Calculation benchmark model A iThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C i" and " nominal unit damage vector D i u";
G. actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector C of monitored amount i"; When numbering to the element of the institute's directed quantity that occurred before this step and this step; Should use same coding rule; Can guarantee element each vector, that numbering is identical of occurring before this step and this step like this, represent same monitored amount, corresponding to vectorial defined relevant information under this element;
H. define the vectorial d of current name damage i cWith current actual damage vector d i, the element number of two damage vectors equals by the quantity of evaluation object, the vectorial d of current name damage i cElement numerical value represent corresponding to the current nominal degree of injury or the support displacement of evaluation object, current actual damage vector d iElement numerical value represent corresponding by the current actual damage degree or the support displacement of evaluation object; The element number of the element of two damage vectors equals by the quantity of evaluation object; The element of two damage vectors and be one-to-one relationship between the evaluation object, the coding rules of the element of two damage vectors are with identical by the coding rule of evaluation object;
I. according to " the current numerical value vector C of monitored amount i" " the initial value of monitored amount vector C together i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + &Delta; C i &CenterDot; d c i Formula 1
J. the current actual damage vector d that utilizes formula 2 to express iWith initial health vector d i oWith the vectorial d of current name damage i cElement between relation, calculate current actual damage vector d iAll elements;
d j i = 1 - ( 1 - d Oj i ) ( 1 - d Cj i ) Formula 2
J=1 in the formula 2,2,3 ..., N;
Current actual damage vector d iElement numerical value represent corresponding to the actual damage degree of evaluation object or actual support displacement, according to current actual damage vector d iJust can define the impaired and degree of injury of which rope, just can confirm actual support displacement; If a certain element of current actual damage vector is corresponding to a rope in the cable system; And its numerical value is 0, representes that the pairing rope of this element is intact, does not damage; If its numerical value is 100%; Represent that then the pairing rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this rope; If a certain element of current actual damage vector is corresponding to displacement component, a d so of a bearing i jRepresent its present bit shift value;
K. try to achieve the vectorial d of current name damage i cAfter, set up mark vector F according to formula 3 i, formula 4 has provided mark vector F iThe definition of j element;
F i = F 1 i F 2 i &CenterDot; &CenterDot; &CenterDot; F j i &CenterDot; &CenterDot; &CenterDot; F N i T Formula 3
F j i = 0 , If d Cj i < D Uj i 1 , If d Cj i &GreaterEqual; D Uj i Formula 4
Element F in the formula 4 i jBe mark vector F iJ element, D i UjBe nominal unit damage vector D i uJ element, d i CjBe the vectorial d of current name damage i cJ element, they all represent j by the relevant information of evaluation object, j=1 in the formula 4,2,3 ..., N;
1. if mark vector F iElement be 0 entirely, then get back to g step and continue this circulation; If mark vector F iElement be not 0 entirely, then get into next step, i.e. m step;
M. the initial damage vector d that calculates next time, promptly circulates required the i+1 time according to formula 5 I+1 oEach element d I+1 Oj
d Oj i + 1 = 1 - ( 1 - d Oj i ) ( 1 - D Uj i F j i ) Formula 5
D in the formula 5 i UjBe nominal unit damage vector D i uJ element, F i jBe mark vector F iJ element, j=1 in the formula 5,2,3 ..., N; Vector d I+1 oThe coding rule of element with identical by the coding rule of evaluation object;
N. at Mechanics Calculation benchmark model A iThe basis on, order is d by the health status of evaluation object I+1 oThe back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time I+1
O. pass through Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1Structure all monitored strains point, with the monitored strain numerical value that should change direction, these numerical value are formed next time, the vectorial C of initial value of the required monitored amount that promptly circulates for the i+1 time I+1 o
P. get back to the f step, beginning is circulation next time.
2. according to claim 1 based on the identification damaged cable of space coordinate monitoring and the progressive method of support displacement, it is characterized in that in step f, at Mechanics Calculation benchmark model A iThe basis on carry out the several times Mechanics Calculation, through calculate obtaining " unit damage monitored numerical quantity transformation matrices Δ C i" and " nominal unit damage vector D i u" concrete grammar be:
F1. at the Mechanics Calculation benchmark model A of Cable Structure iThe basis on carry out the several times Mechanics Calculation, equal N on the calculation times numerical value; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or displacement, to be increased unit damage or unit displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases unit displacement again at this sense of displacement; Increase again in calculating each time unit damage or unit displacement be different from by evaluation object other time increase again in calculating unit damage or unit displacement by evaluation object, with " nominal unit damage vector D i u" write down the unit damage or the unit displacement that increase again of all supposition; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector; Obtain a monitored quantitative changeization vector, have N evaluation object that N monitored quantitative changeization vector just arranged;
F3. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged successively i
CN2010101375622A 2010-03-31 2010-03-31 Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring CN101832872B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101375622A CN101832872B (en) 2010-03-31 2010-03-31 Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101375622A CN101832872B (en) 2010-03-31 2010-03-31 Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring

Publications (2)

Publication Number Publication Date
CN101832872A CN101832872A (en) 2010-09-15
CN101832872B true CN101832872B (en) 2012-05-23

Family

ID=42717019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101375622A CN101832872B (en) 2010-03-31 2010-03-31 Progressive method for identifying damaged cables and support displacement based on spatial coordinate monitoring

Country Status (1)

Country Link
CN (1) CN101832872B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102297768A (en) * 2011-05-13 2011-12-28 东南大学 Progressive method for identifying damaged cable and bearer angular displacement monitored by space coordinate
CN102288424A (en) * 2011-05-13 2011-12-21 东南大学 Progressive method for recognizing damaged cable and angular displacement of support based on hybrid monitoring
CN102252861A (en) * 2011-05-31 2011-11-23 东南大学 Progressive method for identifying support generalized displacement of damaged cable based on space coordinate monitoring
CN102323077A (en) * 2011-05-31 2012-01-18 东南大学 Method for identifying damaged cable, slack cable and support generalized displacement on the basis of spatial coordinate monitoring
CN102323084A (en) * 2011-05-31 2012-01-18 东南大学 Progressive method for identifying damaged cable and support generalized displacement on the basis of mixed monitoring
CN102706579B (en) * 2012-05-29 2015-05-06 东南大学 Method for identifying damaged cable and support translation based on space coordinate monitoring during temperature change
CN102706596B (en) * 2012-05-29 2015-06-10 东南大学 Problem cable and support translation progressive identification method used on basis of strain monitoring during temperature variation
CN102706592B (en) * 2012-05-29 2015-06-10 东南大学 Problem cable and support translation progressive identification method on basis of space coordinate monitoring during temperature variation
CN102721551B (en) * 2012-05-29 2015-06-10 东南大学 Problematic cable and support translation identification method based on space coordinate monitoring during temperature change
CN102706575B (en) * 2012-05-29 2015-04-15 东南大学 Damaged cable and supporting seat translation progressive-type identification method based on space coordinate monitoring at moment of temperature variation
CN102706678B (en) * 2012-05-30 2015-10-07 东南大学 The approximant recognition methods of angular displacement of support temperature variation space coordinate monitoring damaged cable

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101334275A (en) * 2008-07-29 2008-12-31 东南大学 Cable system health monitoring method based on strain monitoring cable structure
CN101387600A (en) * 2008-10-17 2009-03-18 东南大学 Cable system health monitoring method in cable structure based on mixed monitor
CN101476990A (en) * 2009-01-20 2009-07-08 东南大学 Progressive health monitoring method for cable system in cable structure based on angle monitoring
CN101482434A (en) * 2009-01-20 2009-07-15 东南大学 Progressive health monitoring method for cable system in cable structure based on cable force monitoring
CN101504333A (en) * 2009-01-09 2009-08-12 东南大学 Progressive health monitoring method for cable structure median cord system based on strain monitoring

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101334275A (en) * 2008-07-29 2008-12-31 东南大学 Cable system health monitoring method based on strain monitoring cable structure
CN101387600A (en) * 2008-10-17 2009-03-18 东南大学 Cable system health monitoring method in cable structure based on mixed monitor
CN101504333A (en) * 2009-01-09 2009-08-12 东南大学 Progressive health monitoring method for cable structure median cord system based on strain monitoring
CN101476990A (en) * 2009-01-20 2009-07-08 东南大学 Progressive health monitoring method for cable system in cable structure based on angle monitoring
CN101482434A (en) * 2009-01-20 2009-07-15 东南大学 Progressive health monitoring method for cable system in cable structure based on cable force monitoring

Also Published As

Publication number Publication date
CN101832872A (en) 2010-09-15

Similar Documents

Publication Publication Date Title
CN101673316B (en) Method for identifying supporting cables needing to adjust cable force based on hybrid monitoring
CN101387600B (en) Cable system health monitoring method in cable structure based on mixed monitor
CN101587026B (en) Method of identifying support cable in need of regulating cable tension based on cable tension monitoring
CN102706659A (en) Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change
CN100595514C (en) Cable system health monitoring method based on strain monitoring cable structure
CN101476989B (en) Progressive health monitoring method for cable system in cable structure based on mixed monitoring
CN101476990B (en) Progressive health monitoring method for cable system in cable structure based on angle monitoring
CN103604549A (en) Linear displacement progressive type identification method for strain monitoring of damaged cable and intensive load
CN103616222A (en) Hybrid-monitoring linear displacement progressive identification method for defective cable and concentrated loads
CN103604649A (en) Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring
CN102735459A (en) Progressive generalized displacement identification method for problem cables and support based on temperature change tension force monitoring
CN101477340B (en) Progressive health monitoring method for cable system in cable structure based on cable force monitoring
CN103604628A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement angle monitoring
CN103604644A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement strain monitoring
CN103604634A (en) Linear displacement progressive type identification method for space coordinate monitoring of damaged cable and intensive load
CN103616249A (en) Strain-monitoring damaged cable centralized load generalized displacement recognition method
CN101782945B (en) Method for identifying loose supporting ropes based on space coordinate monitoring during support settlement
CN103604623A (en) Progressive type damaged rope and concentrated load identification method based on angular displacement strain monitoring
CN101788403B (en) Progressive method for identifying loose support cable based on strain monitoring during support settlement
CN101782946B (en) Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement
CN103604631A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement cable power monitoring
CN103604636A (en) Linear displacement progressive type identification method for cable force monitoring of damaged cable and intensive load
CN101334338B (en) Cable structure cable system healthy monitoring method based on space coordinates monitoring
CN103616231A (en) Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads
CN103604638A (en) Problematic cable/concentrated load/linear displacement recognition method on basis of strain monitoring

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
TR01 Transfer of patent right

Effective date of registration: 20131018

Address after: 226600 No. 8 Yingbin Road, software park, Haian County, Jiangsu Province

Patentee after: Haian Changzhou University Technology Transfer Center Co., Ltd.

Patentee after: Southeast University

Address before: 210009 No. 2, Southeast University Road, Jiangning Development Zone, Nanjing, Jiangsu

Patentee before: Southeast University

ASS Succession or assignment of patent right

Owner name: SOWTHEAST UNIV.

Effective date: 20131018

Owner name: HAIAN CHANGDA TECHNOLOGY TRANSFER CENTER CO., LTD.

Free format text: FORMER OWNER: SOWTHEAST UNIV.

Effective date: 20131018

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 210009 NANJING, JIANGSU PROVINCE TO: 226600 NANTONG, JIANGSU PROVINCE