CN102323097A - Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement - Google Patents

Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement Download PDF

Info

Publication number
CN102323097A
CN102323097A CN201110144482A CN201110144482A CN102323097A CN 102323097 A CN102323097 A CN 102323097A CN 201110144482 A CN201110144482 A CN 201110144482A CN 201110144482 A CN201110144482 A CN 201110144482A CN 102323097 A CN102323097 A CN 102323097A
Authority
CN
China
Prior art keywords
vector
damage
cable
evaluation
rope
Prior art date
Application number
CN201110144482A
Other languages
Chinese (zh)
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 CN201110144482A priority Critical patent/CN102323097A/en
Publication of CN102323097A publication Critical patent/CN102323097A/en

Links

Abstract

The invention discloses a cable force monitoring based progressive method for recognizing a damaged cable and supporting seat generalized displacement. Based on cable force monitoring, i.e., monitoring of the cable force of all support cables and manually-increased cables and aiming at the defect that linear relations between a current numeric vector of a monitored quantity and an initial numeric vector of the monitored quality, a unit damage monitored quantity variance matrix and a current nominal damage vector are approximate, the invention provides a method for segmentally approximating a nonlinear relation by using a linear relation. A large interval is partitioned into continuous small intervals, the linear relation in each small interval is accurate enough, and the supporting seat generalized displacement and the damaged cable can be recognized quickly with appropriate algorithms such as a multi-target optimization algorithm and the like in each small interval.

Description

Progressive method based on the generalized displacement of cable force monitoring identification damaged cable bearing
Technical field
Structures such as cable-stayed bridge, suspension bridge, truss-frame structure have a common ground; Be exactly that they have many parts that bear tensile load; Like suspension cable, main push-towing rope, hoist cable, pull bar or the like; The common ground of this class formation is to be support unit with rope, cable or the rod member that only bears tensile load, and the present invention is " Cable Structure " with such structure representation for simplicity.In the military service process of Cable Structure; The supporting system of Cable Structure (refers to all ropeway carrying-ropes, reaches all rod members that only bear tensile load that play supporting role; For simplicity; This patent is called " cable system " with whole support unit unifications of this class formation; But in fact cable system not only refers to support rope, also comprises the rod member that only bears tensile load) can be impaired, the bearing of Cable Structure generalized displacement also possibly occur (for example the bearing generalized displacement refers to that bearing is along the angular displacement around X, Y, Z axle of the displacement of the lines of X, Y, Z axle and bearing simultaneously; Corresponding to the bearing generalized displacement; The bearing generalized coordinate refers to that bearing is about the coordinate of X, Y, Z axle and the bearing angular coordinate about X, Y, Z axle); These safety that change to Cable Structure are a kind of threats; The present invention is based on structural health monitoring technology, based on cable force monitoring, adopt progressive method to discern the damaged cable in the cable system of bearing generalized displacement and Cable Structure, belong to the engineering structure health monitoring field.
Background technology
The bearing generalized displacement is a significant threat to Cable Structure safety; Same; Cable system is the key components of Cable Structure normally; Its inefficacy usually brings the inefficacy of total, and the damaged cable of discerning based on structural health monitoring technology in the cable system of bearing generalized displacement and Cable Structure is a kind of method that has potentiality.When generalized displacement appears in bearing or the health status of cable system when changing, or two kinds of situation when taking place simultaneously; Can cause the variation of the measurable parameter of structure; For example can cause the variation of Suo Li; Can influence the distortion or the strain of Cable Structure; Can influence the shape or the volume coordinate of Cable Structure; Can cause the variation (the for example arbitrarily variation of the angle coordinate of the straight line of any this point of mistake in any section of body structure surface, the perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily) of angle coordinate of any imaginary line of the every bit of Cable Structure, all these change the health status information that has all comprised cable system; In fact the variation of these measurable parameters comprised cable system health status information, comprised bearing generalized displacement information, that is to say that the measurable parameter that can utilize structure discerns bearing generalized displacement and damaged cable.
In order reliable monitoring and judgement to be arranged to the health status and the bearing generalized displacement of the cable system of Cable Structure; The method of the variation of a measurable parameter of can rational and effective setting up Cable Structure with the relation between the health status of all ropes in bearing generalized displacement and the cable system must be arranged, and the health monitoring systems of setting up based on this method can provide the health evaluating of more believable bearing generalized displacement assessment and cable system.
Summary of the invention
Technical matters:The invention discloses a kind of based on cable force monitoring, adopt health monitor method progressive method, that can discern bearing generalized displacement and damaged cable rationally and effectively.
Technical scheme:If the quantity sum of the quantity of rope and bearing generalized displacement component does NFor the purpose of narrating conveniently, the present invention is unified to claim that rope and the bearing generalized displacement assessed are " by evaluation object ", gives by the evaluation object serial number, and the present invention is with using variable jRepresent this numbering, j=1,2,3 ..., N, therefore we can say NIndividual by evaluation object.
If it is total in the cable system M 1 Root supporting rope, structure rope force data comprises this M 1 The Suo Li of root supporting rope, obviously M 1 Less than by the quantity of evaluation object NOnly pass through M 1 Individual supporting rope M 1 Individual rope force data is found the solution unknown NIndividual state by evaluation object is impossible, and the present invention in monitoring all M 1 On the basis of root supporting cable force, increase to be no less than ( N-M 1 ) individual other monitored amounts.
Increase be no less than ( N-M 1 ) other individual monitored amounts remain Suo Li, narration as follows:
Structurally artificially increase M 2 (M 2 Be not less than N-M 1 )The root rope increases newly M 2 The rigidity of root rope is compared with the rigidity of any supporting rope of Cable Structure, can be little a lot, for example little 10 times, increase newly M 2 The Suo Li of root rope should be less, and for example its xsect normal stress should be less than its fatigue limit, and these requirements can guarantee to increase newly M 2 Fatigue damage can not take place in the root rope, increases newly M 2 The fully anchoring of the two ends of root rope guarantees can not occur relaxing, and increases newly M 2 The root rope should obtain sufficient anti-corrosion protection, and assurance increases newly M 2 Damage and lax can not take place in the root rope, will monitor in the monitoring structural health conditions process that this increases newly M 2 The Suo Li of root rope.
Comprehensive above-mentioned monitored amount, total is total M (M=M 1 + M 2 )The root rope MIndividual monitored amount, MShould be greater than by the quantity of evaluation object NBecause MIndividual monitored amount all is Suo Li, so the present invention is called " based on the progressive method of cable force monitoring identification damaged cable bearing generalized displacement ".
For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " monitored amount ".Give MIndividual monitored amount serial number, this numbering will be used to generate the vector sum matrix in subsequent step.The present invention is with using variable kRepresent this numbering, k=1,2,3 ..., M
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 rope force data 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 is alphabetical in the present invention iOnly represent cycle index, promptly iInferior circulation.
The iThe Cable Structure that needs when inferior circulation begins " initial health vector d o i " (shown in (1)), use d o i Expression the iCable Structure when inferior circulation begins is (with Mechanics Calculation benchmark model A iThe initial health of Cable Structure expression).
(1)
In the formula (1) d i Oj ( i=1,2,3, ; j=1,2,3 ...., N) expression circulates for the i time when beginning, Mechanics Calculation benchmark model A iIn cable system jIndividual by the current health status of evaluation object, if should be the rope (or pull bar) in the cable system, so by evaluation object d i Represent its current damage, d i Being to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be a generalized displacement component of a bearing by evaluation object, so d i Represent its current generalized displacement numerical value.In the formula (1) TThe transposition of expression vector (back together).
Setting up the initial health vector during circulation beginning for the first time (is designated as according to formula (1) d 1 o ) time, the data and the bearing generalized displacement that utilize the Non-Destructive Testing data etc. of rope can express the health status of rope are measured and are set up by evaluation object initial health vector d 1 o If when not having the data of Non-Destructive Testing data and other health status that can express rope of rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vector d 1 o In get 0 with each element numerical value of Suo Xiangguan.
The i time ( i=2,3,4,5,6 ...) needs was vectorial by the evaluation object initial health when circulation began d i o , be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Set up the Mechanics Calculation benchmark model A of Cable Structure oThe method of (for example finite element benchmark model) is following:
At first on Cable Structure, increase M 2 (M 2 Be not less than N-M 1 )The root rope increases newly M 2 The rigidity of the isostructural any supporting rope of the rigidity of root rope is compared, and is can be little a lot, for example little 10 times, will monitor in the monitoring structural health conditions process that this increases newly M 2 The Suo Li of root rope.Actual measurement obtains that this increases newly before structural healthy monitoring system is started working M 2 The Suo Li of root rope.Measure simultaneously and to increase newly M 2 The geometric parameter of root rope and mechanics parameter measure and increase newly M 2 Two coordinates that end points is installed on Cable Structure of root rope.Deserving to be called the information of stating is what increase newly M 2 All information of root rope.
Increase newly M 2 Behind all ten-fours of root rope, set up A again oSet up A oThe time, increase newly according to known M 2 All information of root rope; Measured data according to the Cable Structure in the Cable Structure completion (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope) and design drawing, as-constructed drawing, utilize mechanics method (for example finite element method) to set up A oIf there is not the measured data of the structure in the Cable Structure completion; So just before setting up health monitoring systems, structure is surveyed; The measured data that obtains Cable Structure (comprises measured datas such as Cable Structure shape data, rope force data, draw-bar pull data, Cable Structure bearing generalized coordinate data, Cable Structure modal data; To cable-stayed bridge, suspension bridge and bridge type data of Yan Shiqiao, rope force data, the modal data of bridge, the Non-Destructive Testing data of rope etc. can be expressed the data of the health status of rope); According to design drawing, the as-constructed drawing of these data and Cable Structure, utilize mechanics method (for example finite element method) to set up A oNo matter which kind of method to obtain A with o, based on A oThe Cable Structure computational data that calculates (to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, rope force data, bridge) must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this oStrain computational data, Suo Li computational data, Cable Structure shape computational data and displacement computational data, Cable Structure angle-data etc. under the analog case of calculating gained, the measured data when truly taking place near institute's analog case reliably.A oBe constant, only when circulation beginning for the first time, set up.
The Mechanics Calculation benchmark model 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
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 o Describe.
The i time ( i=2,3,4,5,6 ...) the Mechanics Calculation benchmark model A of needs when circulation begins i, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Existing Mechanics Calculation benchmark model A 1With vectorial by the evaluation object initial health d 1 o After, model A 1In each by the health status of evaluation object by vector d 1 o Express.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 bearing generalized displacement.Might as well claim model A 0For the not damaged of Cable Structure does not have bearing generalized displacement model A 0
Vectorial among the present invention with monitored amount initial value C i o " ( i=1,2,3 ...) expression the i time ( i=1,2,3,4,5,6 ...) initial value (referring to formula (2)) of the monitored amount of all appointments when circulation begins, C i o Full name be " the initial value vector of the i time monitored amount of circulation ".
(2)
In the formula (2) C i Ok ( i=1,2,3, k=1,2,3, ., M; M>=N;) be in the i time when beginning circulation, the Cable Structure the kIndividual monitored amount.Vector C i o Define by the front MIndividual monitored amount forms according to certain series arrangement, and this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
During circulation beginning for the first time, " the initial value vector of the 1st the monitored amount that circulates C 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 need " initial value of the i time monitored amount of circulation is vectorial C i o ", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
Second step: circulation each time needs to set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector ", and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C i , " nominal unit damage vector " that the i time circulation set up is designated as D i u , i=1,2,3 ...
The Cable Structure " unit damage monitored numerical quantity transformation matrices " that circulation is for the first time set up is designated as Δ C 1 Set up Δ C 1 Process following:
Mechanics Calculation benchmark model A in Cable Structure 1The basis on carry out several times and calculate, equal on the calculation times numerical value NCalculating hypothesis each time has only one by evaluation object unit damage to be arranged; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis has unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage); If should be by evaluation object be the generalized displacement component of a direction of a bearing, just suppose that this bearing is this generalized displacement direction generation unit generalized displacement (for example get 1 millimeter, 2 millimeters, 3 millimeters etc. and be the unit line displacement, get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. and be the unit angular displacement).For narrating conveniently, the damage and the bearing generalized displacement of the supporting rope that the present invention will suppose are referred to as unit damage.Calculate for convenient; When setting unit damage in the circulation each time can all be structural health conditions during this time circulation beginning as being healthy fully, and set on this basis unit damage (in subsequent step, calculate, by the health status numerical value of evaluation object---be called nominal health status vector d 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 for the first time D 1 u Calculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, appointment in front MThe current calculated value of individual monitored amount calculates gained each time MThe current calculated value of individual monitored amount is formed one " the current numerical value vector of the calculating of monitored amount " (when hypothesis the jIndividual when by evaluation object unit damage being arranged, available formula (4) is represented all appointments MThe current numerical value vector of the calculating of individual monitored amount C 1 Tj ); The current numerical value vector of the calculating of the monitored amount that calculates each time deducts the initial value vector of monitored amount C 1 o , the gained vector be exactly under this condition (with have unit damage by the mark that is numbered of evaluation object) " the numerical value change vector of monitored amount " (when the jIndividual when unit damage being arranged by evaluation object, use δ C 1 j The numerical value change vector of representing monitored amount, δ C 1 j Definition see formula (5), formula (6) and formula (7), formula (5) deducts after the formula (2) again divided by vector for formula (4) D 1 u jIndividual element D Uj Gained), the numerical value change of monitored amount vector δ C 1 j Each element representation since that supposition has a unit damage when calculating by evaluation object (for example the jIndividual by evaluation object) there is unit damage (for example D Uj ), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage of supposition D Uj Rate of change; Have NIndividual just had by evaluation object NIndividual " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has M(it is general, M>=N) individual element, by this NIndividual " the numerical value change vector of monitored amount " formed successively to be had M * N" the unit damage monitored numerical quantity transformation matrices of individual element Δ C 1 " (MOK NRow ), each vector δ C 1 j ( j=1,2,3 ...., N) be matrix Δ C 1 One row, Δ C 1 Definition suc as formula shown in (8).
(3)
Nominal unit damage vector in the formula (3) D i u Element D i Uj ( i=1,2,3, j=1,2,3 ...., N) expression the iSuppose in the inferior circulation jIndividual by the unit damage numerical value of evaluation object, vector D i u In the numerical value of each element can be the same or different.
(4)
Element in the formula (4) C i Tjk ( i=1,2,3, j=1,2,3 ...., N; k=1,2,3 ...., M; M>=N) expression the iInferior circulation is because the jIndividual when unit damage being arranged by evaluation object, according to coding rule pairing kThe current numerical value of the calculating of the monitored amount of individual appointment.
(5)
The subscript of each amount in the formula (5) i( i=1,2,3 ...) expression the iInferior circulation, subscript j( j=1,2,3 ...., N) expression the jIndividual had unit damage by evaluation object, in the formula D i Uj It is vector D i u In jIndividual element.Vector δ C i j Definition suc as formula shown in (6), δ C i j k( k=1,2,3 ...., M; M>=N) individual element δ C i Jk Expression the iIn the inferior circulation, set up matrix Δ C i The time, suppose jThe individual gained that calculates when unit damage being arranged by evaluation object kThe change amount of individual monitored amount is with respect to the unit damage of supposition D i Uj Rate of change, it defines suc as formula shown in (7).
(6)
(7)
The definition of each amount was narrated in front in the formula (7).
(8)
Vector in the formula (8) δ C i j ( i=1,2,3 ....,, j=1,2,3 ...., N) expression the iIn the inferior circulation, because the jIndividual had unit damage by evaluation object D i Uj Cause, the relative value of all monitored amounts changes.Matrix Δ C i Row (subscript j) coding rule and front vector d i o The subscript of element jCoding rule identical.
The 3rd step: identification is by the current health status of evaluation object.Detailed process is following.
The i( i=1,2,3 ...) in the inferior circulation, " current (calculating or actual measurement) numerical value vector of monitored amount C i " " initial value of monitored amount is vectorial together C i o ", " unit damage monitored numerical quantity transformation matrices Δ C i " and " current name damage vector d i c " between linear approximate relationship, shown in (9) or formula (10).
(9)
(10)
Current (calculating or actual measurement) numerical value vector of monitored amount in formula (9) and the formula (10) C i Definition be similar to the initial value vector of monitored amount C i o Definition, see formula (11); By the current name damage of evaluation object vector d i c Definition see formula (12).
(11)
Element in the formula (11) C i k ( i=1,2,3 ....; k=1,2,3 ...., M; M>=N) be iInferior circulation time Cable Structure, be numbered according to coding rule is pairing kThe current numerical value of monitored amount.
(12)
In the formula (12) d i Cj ( i=1,2,3 ....; j=1,2,3 ...., N) be iCable Structure in the inferior circulation jIndividual by the current nominal impairment value of evaluation object, vector d i c The subscript of element jCoding rule and matrix Δ C i The coding rule of row identical.
When by evaluation object actual damage or bearing generalized displacement when not too big; Because the Cable Structure material still is in the linear elasticity stage; The distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (9) or formula (10) is less with the error of actual conditions, and error can be used error vector e i (formula (13)) definition, the error of linear relationship shown in expression (9) or the formula (10).
(13)
In the formula (13) Abs ()Be the function that takes absolute value, each element of the vector of trying to achieve in the bracket is taken absolute value.
Because there are certain error in formula (9) or the represented linear relationship of formula (10), therefore can not be simply according to formula (9) or formula (10) and " current (actual measurement) numerical value vector of monitored amount C i " come directly to find the solution to obtain current name damage vector d i c And it is vectorial to obtain current name damage d i c Acceptable separating (promptly have reasonable error, but can be more accurately from cable system, confirm damaged cable position and degree of injury thereof, confirm bearing generalized displacement amount) become a rational solution, available formula (14) is expressed this method.
(14)
In the formula (14) Abs ()Be the function that takes absolute value, vector g i Description departs from ideal linearity relation (formula (9) or formula (10))
Reasonable deviation, define by formula (15).
(15)
In the formula (15) g i k ( i=1,2,3 ....; k=1,2,3 ...., M) described iThe maximum allowable offset that departs from the ideal linearity relation shown in formula (9) or the formula (10) in the inferior circulation.Vector g i Can be according to the error vector of formula (13) definition e i Tentative calculation is selected.
Initial value vector in monitored amount C i o (survey or calculate), unit damage monitored numerical quantity transformation matrices Δ C i The current numerical value vector of (calculating) and monitored amount C i When (actual measurement obtains) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (14), obtain current name damage vector d i c Acceptable separating, current actual damage vector d i The element of (formula (16) is seen in definition) can calculate according to formula (17), has just obtained by the current actual damage vector of evaluation object d i Thereby, can by d i Confirm damaged cable position and degree of injury, confirm bearing generalized displacement amount, realized that just damage identification and bearing generalized displacement discern.
(16)
In the formula (16) d i j ( i=1,2,3, ; j=1,2,3 ...., N) expression the iIn the inferior circulation jIndividual by the actual damage value of evaluation object, formula (17) is seen in its definition, if should be the rope (or pull bar) in the cable system by evaluation object, so d i j Represent its current damage, d i j Being to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%, if should be a generalized displacement component of a bearing by evaluation object, so d i j Represent its current generalized displacement numerical value.Vector d i Coding rule and the formula (1) of element in vector d i o The coding rule of element identical.
(17)
In the formula (17) d i Oj ( i=1,2,3,4, ; j=1,2,3 ...., N) be vector d i o jIndividual element, d i Cj It is vector d i c jIndividual element.
The 4th step: judge whether to finish this (the iInferior) circulation, if then accomplish the tailing in work before this loop ends, for next time (promptly the i+ 1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following.
This ( iInferior) try to achieve current name damage vector in the circulation d i c After, at first, set up mark vector according to formula (18) F i , formula (19) has provided mark vector F i jThe definition of individual element; If mark vector F i Element be 0 entirely, then in this circulation, continue health monitoring and calculating to Cable Structure; If mark vector F i Element be not 0 entirely, then accomplish subsequent step after, get into circulation next time.So-called subsequent step is: at first, calculate next time (promptly according to formula (20) i+ 1 time, i=1,2,3,4 ...) the required initial damage vector of circulation d i+ 1 o Each element d i+ 1 Oj The second, at Mechanics Calculation benchmark model A i( i=1,2,3,4 ...) or the not damaged model A of Cable Structure 0The basis on, the order by the health status situation of evaluation object done d i+ 1 o The back upgrade and to obtain next time (the i+1 time, i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation I+1At last, through to Mechanics Calculation benchmark model A I+1The initial value that calculates monitored amount, by its form next time (promptly the i+1 time, i=1,2,3,4 ...) required " the initial value vector of monitored amount of circulation C I+1 o " ( i=1,2,3,4 ...).
(18)
Mark vector in the formula (18) F i Subscript iExpression the iInferior circulation, its element F i j ( j=1,2,3 ..., subscript N) jExpression the jIndividual by the damage characteristic of evaluation object, can only get 0 and 1 two amount, concrete value rule is seen formula (19).
(19)
Element in the formula (19) F i j It is mark vector F i jIndividual element, D i Uj It is nominal unit damage vector D i u jIndividual element (seeing formula (3)), d i Cj It is current name damage vector d i c jIndividual element (seeing formula (12)), they all represent jIndividual by the relevant information of evaluation object.
(20)
In the formula (20) D i Uj It is nominal unit damage vector D i u jIndividual element (seeing formula (3)), d i Cj It is current name damage vector d i c jIndividual element (seeing formula (12)).
 
Second portion of the present invention: the software and hardware part of health monitoring systems.
Hardware components comprises monitored amount monitoring system, signal picker and computing machine etc.Require to monitor in real time or quasi real time each monitored amount.
Software should be used following function by tool: software section should be accomplished the process that first of the present invention sets, and promptly accomplishes functions such as needed among the present invention, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.
 
The inventive method specifically comprises:
A. for for the purpose of narration is convenient, the present invention is unified to claim that the supporting rope of being assessed is by evaluation object with bearing generalized displacement component, establishes the quantity of the supporting rope of being assessed and the quantity sum of bearing generalized displacement component and does N, promptly done by the quantity of evaluation object NConfirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention uses variable jRepresent this numbering, j=1,2,3 ..., N
B. establish in the cable system total M 1 Root supporting rope, structure rope force data comprises this M 1 The Suo Li of root supporting rope, obviously M 1 Less than by the quantity of evaluation object NOnly pass through M 1 Individual supporting rope M 1 Individual rope force data is found the solution unknown NIndividual state by evaluation object is impossible, and the present invention in monitoring all M 1 On the basis of root supporting cable force, structurally artificially increase M 2 The root rope will monitor in the monitoring structural health conditions process that this increases newly M 2 The Suo Li of root rope; Comprehensive above-mentioned monitored amount, total is total MThe root rope MIndividual Suo Li is monitored, promptly has MIndividual monitored amount, M whereinFor M 1 With M 2 Sum; MShould be greater than by the quantity of evaluation object NIncrease newly M 2 The rigidity of root rope is compared with the rigidity of any supporting rope of Cable Structure, should be little many; Increase newly M 2 The Suo Li of root rope should than the Suo Lixiao of any supporting rope of Cable Structure many, even can guarantee what this increased newly like this M 2 Damage or lax has appearred in the root rope, and is very little to the influence of the stress of other members of Cable Structure, strain, distortion; Increase newly M 2 Normal stress should be less than its fatigue limit on the xsect of root rope, and these requirements can guarantee to increase newly M 2 Fatigue damage can not take place in the root rope; Increase newly M 2 The fully anchoring of the two ends of root rope guarantees can not occur relaxing; Increase newly M 2 The root rope should obtain sufficient anti-corrosion protection, and assurance increases newly M 2 Damage and lax can not take place in the root rope; For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " monitored amount "; Give MIndividual monitored amount serial number, the present invention is with using variable jRepresent this numbering, j=1,2,3 ..., M, this numbering will be used to generate the vector sum matrix in subsequent step;
C. utilize to be expressed by the data of the health status of evaluation object and set up by evaluation object initial health vector by the Non-Destructive Testing data of evaluation object etc. d i o If during not by the Non-Destructive Testing data of evaluation object, vector d i o Each element numerical value get 0; Vector d i o The coding rule of element with identical by the coding rule of evaluation object; The present invention representes cycle index with i, i=1, and 2,3, Here be circulation for the first time, i gets 1, the initial health vector of promptly setting up here d i o Can be embodied as d 1 o
D. setting up the initial health vector d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount C i o Here be circulation for the first time, i gets 1, the initial value vector of the monitored amount of promptly setting up here C i o Can be embodied as C 1 o Obtain monitored amount initial value vector in actual measurement C 1 o The time, actual measurement obtains the initial rope force data of all ropes of Cable Structure, the initial geometric data and the initial Cable Structure bearing generalized coordinate data of structure; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure bearing generalized coordinate data are set up Cable Structure 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 of monitored amount C i "; When numbering to the element of the institute's directed quantity that occurred before this step and this step; Should use same coding rule; Can guarantee 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. in the monitoring structural health conditions process, to what increase newly M 2 Root Suo Jinhang Non-Destructive Testing therefrom identifies the rope that damage occurs or relax;
I. according to monitored amount coding rule, from the initial value vector of monitored amount C i o The middle element of removing appearance damage that identifies among the step h or the rope correspondence that relaxes; According to monitored amount coding rule, from unit damage monitored numerical quantity transformation matrices Δ C i The middle row of removing appearance damage that identifies among the step h or the rope correspondence that relaxes; According to monitored amount coding rule, from the current numerical value vector of monitored amount C i The middle element of removing appearance damage that identifies among the step h or the rope correspondence that relaxes;
J. define current name damage vector d i c With current actual damage vector d i , the element number of two damage vectors equals by the quantity of evaluation object, current name damage vector d i c Element numerical value represent corresponding to the current nominal degree of injury of evaluation object or bearing generalized displacement, current actual damage vector d i Element numerical value represent corresponding by the current actual damage degree or the bearing generalized displacement of evaluation object; The element number of the element of two damage vectors equals by the quantity of evaluation object; The element of two damage vectors and be one-to-one relationship between the evaluation object, the coding rules of the element of two damage vectors are with identical by the coding rule of evaluation object;
K. according to " the current numerical value vector of monitored amount C i " " initial value of monitored amount is vectorial together C i o ", " unit damage monitored numerical quantity transformation matrices Δ C i " and " current name damage vector d i c " between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes in the formula 1 d i c Other outer amount is known, finds the solution formula 1 and just can calculate current name damage vector d i c
formula 1
L. the current actual damage vector that utilizes formula 2 to express d i With the initial damage vector d i o With current name damage vector d i c Element between relation, calculate current actual damage vector d i All elements;
formula 2
In the formula 2 j=1,2,3 ..., N;
Current actual damage vector d i Element numerical value represent corresponding to the actual damage degree of evaluation object or actual bearing generalized displacement, according to current actual damage vector d i Just can define the impaired and degree of injury of which rope, just can confirm actual bearing generalized displacement; If a certain element of current actual damage vector is corresponding to being a rope in the cable system; And its numerical value is 0, representes that the pairing rope of this element is intact, does not damage; If its numerical value is 100%; Represent that then the pairing rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this rope; If a certain element of current actual damage vector is corresponding to a generalized displacement component of a bearing, so d i j Represent its current generalized displacement numerical value;
M. try to achieve current name damage vector d i c After, set up mark vector according to formula 3 F i , formula 4 has provided mark vector F i jThe definition of individual element;
formula 3
formula 4
Element in the formula 4 F i j It is mark vector F i jIndividual element, D i Uj It is nominal unit damage vector D i u jIndividual element, d i Cj It is current name damage vector d i c jIndividual element, they all represent jIndividual by the relevant information of evaluation object, in the formula 4 j=1,2,3 ..., N;
If mark vector n. F i Element be 0 entirely, then get back to g step and continue this circulation; If mark vector F i Element be not 0 entirely, then get into next step, i.e. o step;
O. according to formula 5 calculate next time, promptly iThe initial damage vector that+1 circulation is required d i+ 1 o Each element d i+ 1 Oj
formula 5
In the formula 5 D i Uj It is nominal unit damage vector D i u jIndividual element, d i Cj It is current name damage vector d i c jIndividual element, F i j It is mark vector F i jIndividual element is in the formula 5 j=1,2,3 ..., N; Vector d i+ 1 o The coding rule of element with identical by the coding rule of evaluation object;
P. at Mechanics Calculation benchmark model A iThe basis on, the order by the health status of evaluation object done d i+ 1 o The back renewal obtains next time, required Mechanics Calculation benchmark model A promptly circulates for the i+1 time I+1
Q. pass through Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1Structure all monitored strains point, with the monitored strain numerical value that should change direction, these numerical value are formed next time, the initial value of the required monitored amount that promptly circulates for the i+1 time is vectorial C I+1 o
R. 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 on the calculation times numerical value NAccording to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or generalized displacement, to be increased unit damage or unit generalized displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the generalized displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases the unit generalized displacement again in this generalized displacement direction; Increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D i u " unit damage or the unit generalized displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit generalized displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector, obtains a monitored quantitative changeization vector, has NIndividual just had by evaluation object NIndividual monitored quantitative changeization vector;
F3. by this NIndividual monitored quantitative change vector according to NIndividual by the coding rule of evaluation object, forming successively has NThe monitored amount unit change of the Cable Structure matrix of row Δ C 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 generalized displacement of Cable Structure bearing and cable system simultaneously, and system and method disclosed by the invention is very useful to the safety of Cable Structure.
Embodiment
To the health monitoring of Cable Structure, the invention discloses a kind of system and method that can monitor in the Cable Structure health status of each root rope and each bearing generalized displacement component in the cable system rationally and effectively simultaneously.The following explanation of embodiments of the invention in fact only is exemplary, and purpose never is to limit application of the present invention or use.
The present invention adopts a kind of algorithm, and this algorithm is used for monitoring the health status of the cable system of Cable Structure.During practical implementation, the following step is a kind of in the various steps that can take.
The first step: for the purpose of narrating conveniently, the present invention is unified to claim that supporting rope and the bearing generalized displacement component assessed are by evaluation object, establishes the quantity of the supporting rope of being assessed and the quantity sum of bearing generalized displacement component and does N, promptly done by the quantity of evaluation object NConfirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The present invention uses variable jRepresent this numbering, j=1,2,3 ..., N
If it is total in the cable system M 1 Root supporting rope, structure rope force data comprises this M 1 The Suo Li of root supporting rope, obviously M 1 Less than by the quantity of evaluation object NOnly pass through M 1 Individual supporting rope M 1 Individual rope force data is found the solution unknown NIndividual state by evaluation object is impossible, and the present invention in monitoring all M 1 On the basis of root supporting cable force, increase to be no less than ( N-M 1 ) individual other monitored amounts.
Increase be no less than ( N-M 1 ) other individual monitored amounts remain Suo Li, narration as follows:
Structurally artificially increase M 2 (M 2 Be not less than N-M 1 )The root rope increases newly M 2 The rigidity of root rope is compared with the extensional rigidity of any supporting rope of Cable Structure, can be little a lot, for example little 10 times; Increase newly M 2 The Suo Li of root rope should than the Suo Lixiao of any supporting rope of Cable Structure many, even can guarantee what this increased newly like this M 2 Damage or lax has appearred in the root rope, and is very little to the influence of the stress of other members of Cable Structure, strain, distortion; Increase newly M 2 Normal stress should for example have only 1/2nd of fatigue limit less than its fatigue limit on the xsect of root rope, and these requirements can guarantee to increase newly M 2 Fatigue damage can not take place in the root rope; Increase newly M 2 The fully anchoring of the two ends of root rope guarantees can not occur relaxing; Increase newly M 2 The root rope should obtain sufficient anti-corrosion protection, and assurance increases newly M 2 Damage and lax can not take place in the root rope; Can also adopt the modes that increase rope to guarantee the reliability of health monitoring more, for example make M 2 Be not less than N-M 1 2 times, the rope force data of in the monitoring structural health conditions process, only selecting intact rope wherein (is called actual operable monitored amount, writes down its quantity and do K, KMust not less than N) and the corresponding monitored amount unit change of Cable Structure matrix Δ CCarry out the health status assessment, because M 2 Be not less than N-M 1 2 times, can guarantee actual operable; To monitor in the monitoring structural health conditions process that this increases newly M 2 The Suo Li of root rope.Increase newly M 2 The position that the root rope should be installed structurally, personnel are easy to arrive is convenient to personnel it is carried out Non-Destructive Testing.
Comprehensive above-mentioned monitored amount, total is total M (M=M 1 + M 2 )The root rope MIndividual monitored amount, MMust not be less than by the quantity of evaluation object NBecause MIndividual monitored amount all is Suo Li, so the present invention is called " based on the identification damaged cable of cable force monitoring and the health monitor method of bearing generalized displacement ".Give MIndividual monitored amount serial number, the present invention is with using variable jRepresent this numbering, j=1,2,3 ..., M, this numbering will be used to generate the vector sum matrix in subsequent step.
For simplicity, in the present invention " all monitored parameters of structure " are abbreviated as " 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 by the Non-Destructive Testing data of evaluation object etc. d 1 o If during not by the Non-Destructive Testing data of evaluation object, vector d 1 o Each element numerical value get 0; Vector d 1 o The coding rule of element with identical by the coding rule of evaluation object.
The 3rd step: at the initial health vector d 1 o The time, directly measurement calculates the initial value of all monitored amounts of Cable Structure, forms the initial value vector of monitored amount C 1 o
The 4th step: the initial value vector that obtains monitored amount in actual measurement C 1 o The time, can adopt ripe measuring method to carry out cable force measurement, strain measurement, measurement of angle and volume coordinate and measure.Simultaneously; Calculate the initial Suo Li and the Cable Structure original geometric form data (is exactly its initial bridge type data for cable-stayed bridge) of all ropes of Cable Structure after directly measuring or measuring; The original geometric form data of Cable Structure can be the spatial datas that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is just can confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the original geometric form data can be the spatial datas that the spatial data of the end points of all ropes adds some points on the bridge two ends, so-called bridge type data that Here it is.
According to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure (data such as initial Suo Li, structural modal data that comprise structure original geometric form data, strain data, all ropes, to cable-stayed bridge, suspension bridge and the modal data of the bridge type data of Yan Shiqiao, strain data, rope force data, bridge), the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure bearing generalized coordinate data are set up Cable Structure o, based on Mechanics Calculation benchmark model A oThe computational data that calculates structure must be very near its measured data, and error generally must not be greater than 5%.
A oBe constant, only when circulation beginning for the first time, set up; The Mechanics Calculation benchmark model of the Cable Structure of setting up during the i time circulation beginning is designated as A 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 cable force measurement 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 i With nominal unit damage vector D i u Concrete grammar is:
A. when the i time circulation beginning, at the Mechanics Calculation benchmark model A of Cable Structure iThe basis on carry out the several times Mechanics Calculation, equal on the calculation times numerical value NAccording to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one to be increased on the basis of original damage or generalized displacement by evaluation object again unit damage or unit generalized displacement are arranged; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again; If should be the generalized displacement component of a direction of a bearing, just suppose that this bearing increases the unit generalized displacement again in this generalized displacement direction and (if should be the translational component of the x direction of a bearing by evaluation object for example, just supposes that this bearing has the unit line displacement in the x direction by evaluation object; If should be the angular displacement component around the x axle of a bearing by evaluation object; Just suppose that this bearing has the unit angular displacement around the x axle), increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D i u " unit damage or the unit generalized displacement that increase again of all supposition of record record, wherein iExpression the iThe current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure is calculated in inferior circulation each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount and calculated current numerical value vector; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed same monitored amount or same by the relevant information of evaluation object object.
B. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit generalized displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector, obtains a monitored quantitative changeization vector δ C i j Have NIndividual just had by evaluation object NIndividual monitored quantitative changeization vector δ C i j ( J=1,2,3 ..., N).
C. by this NIndividual monitored quantitative change vector according to NIndividual by the coding rule of evaluation object, forming successively has NThe monitored amount unit change of the Cable Structure matrix of row Δ C i " unit damage monitored quantitative change matrix Δ C i " the current name damage vector of coding rule and back definition of row d i c With current actual damage vector d i The element coding rule identical.
In this step, reach when giving each vectorial element numbering thereafter; Should use same coding rule with other vector among the present invention; Can guarantee any element in each vector in this step like this; With element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target.
The 9th step: set up the linear relationship error vector e i And vector g i Utilize data (" the initial value vector of monitored amount of front C i o ", " unit damage monitored quantitative change matrix Δ C i "); when the 8th step calculated each time; promptly in calculating each time, have only in the hypothesis cable system one by evaluation object on the basis of original damage or generalized displacement, increase again unit damage or unit generalized displacement are arranged in, calculate each time and form a health status vector d i t , the health status vector d i t Element number equal by the quantity of evaluation object, vector d i t All elements in have only the numerical value of an element to get to calculate each time in hypothesis increase unit damage value or the unit generalized displacement value of increase of the rope of unit damage, d i t The numerical value of other element get 0, what that was not that numbering and the supposition of 0 element increase unit damage or unit generalized displacement is identical by the corresponding relation of evaluation object, with the element of the same numberings of other vectors with the corresponding relation of this rope; Will C i Tj , C i o , Δ C i , d i t Bring formula (13) into, formula (13) d i c Use d i t Bring into, obtain a linear relationship error vector e i , calculate a linear relationship error vector each time e i Have NIndividual just had by evaluation object NInferior calculating just has NIndividual linear relationship error vector e i , with this NIndividual linear relationship error vector e i Obtain a vector after the addition, with this vector each element divided by NAfter the new vector that obtains be exactly final linear relationship error vector e i Vector g i Equal final error vector e i With vector g i Be kept on the hard disc of computer of operation health monitoring systems software, supply health monitoring systems software to use.Will " the initial value vector of monitored amount C i o " and " unit damage monitored quantitative change matrix Δ C i " etc. parameter be kept on the hard disc of computer of operation health monitoring systems software with the mode of data file.
The tenth step: actual measurement obtain Cable Structure all specify the current measured value of monitored amount, form " the current numerical value vector of monitored amount C i ".
The 11 step: to what increase newly M 2 Root Suo Jinhang Non-Destructive Testing, for example UT (Ultrasonic Testing), visual examination, infrared imaging inspection therefrom identify the rope that damage occurs or relax.
The 12 step: according to monitored amount coding rule, from the initial value vector of monitored amount C i o The middle element of removing appearance damage that identifies in the 11 step or the rope correspondence that relaxes; According to monitored amount coding rule, from unit damage monitored quantitative change matrix Δ C i The middle row of removing appearance damage that identifies in the 11 step or the rope correspondence that relaxes; Current numerical value vector from monitored amount C i The middle element of removing appearance damage that identifies in the 11 step or the rope correspondence that relaxes; According to monitored amount coding rule, from vector gThe middle element of removing appearance damage that identifies in the 11 step or the rope correspondence that relaxes.
The 13 step: according to " the current numerical value vector of monitored amount C i " " initial value of monitored amount is vectorial together C i o ", " unit damage monitored quantitative change matrix Δ C i " and " current name damage vector d i c " between the linear approximate relationship (formula (9)) that exists, calculate the current name damage of cable system vector according to multi-objective optimization algorithm d i c Noninferior solution.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted Sum Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution current name damage vector with the goal programming method only d i c Process, the 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) γ i Be a real number, RBe real number field, area of space Ω has limited vector d i c Span (the present embodiment requirements vector of each element d i c Each element be not less than 0, be not more than 1).The meaning of formula (21) is to seek the minimum real number of an absolute value γ i , make formula (22) be met.In the formula (22) G (d i c )By formula (23) definition, weighing vector in the formula (22) W i With γ i Product representation formula (22) in G (d i c )With vector g i Between the deviation that allows, g i Definition referring to formula (15), its value will the 8th the step calculate.Vector during actual computation W i Can with vector g i Identical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Just can be according to the goal programming method in the hope of current name damage vector d i c
(21)
(22)
(23)
Try to achieve current name damage vector d i c After ,Can be vectorial according to the current actual damage that formula (17) obtain d i Each element, current actual damage vector d i Have reasonable error exactly but can confirm more exactly damaged cable the separating of position and degree of injury thereof, can confirm separating of all bearing generalized displacements more exactly. d i Each element corresponding to one by the health status of evaluation object; If should be the rope (or pull bar) in the cable system by evaluation object; Its current damage of the numeric representation of this element so; The numerical value of this element is to represent this rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of this rope forfeiture corresponding proportion in the time of between 0 and 100%; If being somebody's turn to do by evaluation object be a generalized displacement component of a bearing, so its current generalized displacement numerical value of the numeric representation of this element.
The 14 step: in this circulation, promptly the iTry to achieve current name damage vector in the inferior circulation d i c After, set up mark vector according to formula (18), formula (19) F i If mark vector F i Element be 0 entirely, then got back to for the tenth step and continue this circulation; If mark vector F i Element be not 0 entirely, then get into next step, i.e. the 15 step.
The 15 step: according to formula (20) calculate next time, promptly the iThe initial damage vector that+1 circulation is required d i+ 1 o Each element d i+ 1 Oj
The 16 step: at Cable Structure Mechanics Calculation benchmark model A iThe basis on, order is the last vector that calculates of step by the health status of evaluation object d i+ 1 o After, obtain new Mechanics Calculation benchmark model, next time promptly the required Mechanics Calculation benchmark model A of (the i+1 time) circulation I+1
The 17 step: through to Mechanics Calculation benchmark model A I+1Calculate corresponding to model A I+1The numerical value of all monitored amounts of structure, these numerical value are formed next time, required vector promptly circulates for the i+1 time C I+1 o , promptly the initial value of monitored amount is vectorial
The 18 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 19 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 20 step: got back to for the 7th step, beginning is circulation next time.

Claims (2)

1. progressive method based on the generalized displacement of cable force monitoring identification damaged cable bearing is characterized in that said method comprises:
A. for for the purpose of narration is convenient, unifiedly claim that the supporting rope of being assessed is that establishing the quantity of the supporting rope of being assessed and the quantity sum of bearing generalized displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with bearing generalized displacement component; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; Represent this numbering with variable j, j=1,2,3 ..., N;
B. establish total M in the cable system 1Root supporting rope, structure rope force data comprises this M 1The Suo Li of root supporting rope, obviously M 1Less than by the quantity N of evaluation object; Only pass through M 1The M of individual supporting rope 1It is impossible by the state of evaluation object that individual rope force data is found the solution unknown N, at the whole M of monitoring 1On the basis of root supporting cable force, structurally artificially increase M 2The root rope will be monitored the M that this increases newly in the monitoring structural health conditions process 2The Suo Li of root rope; Comprehensive above-mentioned monitored amount, M Suo Li of the total M root rope of total is monitored, and M monitored amount promptly arranged, and wherein M is M 1With M 2Sum; M should be greater than by the quantity N of evaluation object; The M that increases newly 2The rigidity of root rope is compared with the rigidity of any supporting rope of Cable Structure, should be little many; The M that increases newly 2The Suo Li of root rope should than the Suo Lixiao of any supporting rope of Cable Structure many, even can guarantee the M that this increases newly like this 2Damage or lax has appearred in the root rope, and is very little to the influence of the stress of other members of Cable Structure, strain, distortion; The M that increases newly 2Normal stress should be less than its fatigue limit on the xsect of root rope, and these requirements can guarantee the M that increases newly 2Fatigue damage can not take place in the root rope; The M that increases newly 2The fully anchoring of the two ends of root rope guarantees can not occur relaxing; The M that increases newly 2The root rope should obtain sufficient anti-corrosion protection, guarantees the M that increases newly 2Damage and lax can not take place in the root rope; For simplicity, " all monitored parameters of structure " are abbreviated as " monitored amount "; Give M monitored amount serial number, the present invention is with representing this numbering with variable j, j=1, and 2,3 ..., M, this numbering will be used to generate the vector sum matrix in subsequent step;
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; Represent cycle index with i, i=1,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 bearing generalized coordinate data of structure; The bearing generalized coordinate comprises two kinds of line amount and angle amounts;
E. according to the measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, the Non-Destructive Testing data of rope and the Mechanics Calculation benchmark model A that initial Cable Structure bearing generalized coordinate data are set up Cable Structure 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. in the monitoring structural health conditions process, to the M that increases newly 2Root Suo Jinhang Non-Destructive Testing therefrom identifies the rope that damage occurs or relax;
I. according to monitored amount coding rule, from the initial value vector C of monitored amount i oThe middle element of removing appearance damage that identifies among the step h or the rope correspondence that relaxes; According to monitored amount coding rule, from unit damage monitored numerical quantity transformation matrices Δ C iThe middle row of removing appearance damage that identifies among the step h or the rope correspondence that relaxes; According to monitored amount coding rule, from the current numerical value vector C of monitored amount iThe middle element of removing appearance damage that identifies among the step h or the rope correspondence that relaxes;
J. 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 of evaluation object or bearing generalized displacement, current actual damage vector d iElement numerical value represent corresponding by the current actual damage degree or the bearing generalized displacement of evaluation object; The element number of the element of two damage vectors equals by the quantity of evaluation object; The element of two damage vectors and be one-to-one relationship between the evaluation object, the coding rules of the element of two damage vectors are with identical by the coding rule of evaluation object;
K. 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 + Δ C i · d c i Formula 1
L. 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 bearing generalized displacement, according to current actual damage vector d iJust can define the impaired and degree of injury of which rope, just can confirm actual bearing generalized displacement; If a certain element of current actual damage vector is corresponding to being a rope in the cable system; And its numerical value is 0, representes that the pairing rope of this element is intact, does not damage; If its numerical value is 100%; Represent that then the pairing rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this rope; If a certain element of current actual damage vector is corresponding to a generalized displacement component of a bearing, so Represent its current generalized displacement numerical value;
M. 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 · · · F j i · · · 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 in the formula 4 Be 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 n. 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. o step;
O. 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;
P. 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
Q. 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
R. get back to the f step, beginning is circulation next time.
2. the progressive method based on the generalized displacement of cable force monitoring identification damaged cable bearing according to claim 1 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 generalized displacement, to be increased unit damage or unit generalized displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage again, if should be the generalized displacement component of a direction of a bearing by evaluation object, just supposes that this bearing increases the unit generalized displacement again in this generalized displacement direction; Increase again in calculating each time unit damage or unit generalized displacement be different from by evaluation object other time increase again in calculating unit damage or unit generalized displacement by evaluation object, with " nominal unit damage vector D i u" unit damage or the unit generalized displacement that increase again of all supposition of record record; wherein i representes the i time circulation; calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, and the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount, and to calculate current numerical value vectorial;
F2. the monitored amount that calculates is each time calculated and is calculated unit damage or the unit generalized displacement numerical value of being supposed divided by this time again after current numerical value vector deducts monitored amount initial value vector; Obtain a monitored quantitative changeization vector, have N evaluation object that N monitored quantitative changeization vector just arranged;
F3. individual by the coding rule of evaluation object by this N monitored quantitative change vector according to N, form the monitored amount unit change of the Cable Structure matrix Δ C that the N row are arranged successively i
CN201110144482A 2011-05-31 2011-05-31 Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement CN102323097A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110144482A CN102323097A (en) 2011-05-31 2011-05-31 Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110144482A CN102323097A (en) 2011-05-31 2011-05-31 Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement

Publications (1)

Publication Number Publication Date
CN102323097A true CN102323097A (en) 2012-01-18

Family

ID=45450886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110144482A CN102323097A (en) 2011-05-31 2011-05-31 Cable force monitoring based progressive method for recognizing damaged cable and supporting seat generalized displacement

Country Status (1)

Country Link
CN (1) CN102323097A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706670A (en) * 2012-05-30 2012-10-03 东南大学 Damaged cable and support generalized displacement identification method based on cable force monitoring of temperature change
CN102706651A (en) * 2012-05-29 2012-10-03 东南大学 Damaged cable and support translation progressive identification method on basis of cable force monitoring during temperature variation

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998057166A1 (en) * 1997-06-11 1998-12-17 Pure Technologies Ltd. Method and apparatus for monitoring of tensioned cables
CN101782946A (en) * 2010-03-17 2010-07-21 东南大学 Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement
CN101788403A (en) * 2010-03-17 2010-07-28 东南大学 Progressive method for identifying loose support cable based on strain monitoring during support settlement
CN101806670A (en) * 2010-04-06 2010-08-18 苏州大学 Method and device for detecting missing parts of forging press

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998057166A1 (en) * 1997-06-11 1998-12-17 Pure Technologies Ltd. Method and apparatus for monitoring of tensioned cables
CN101782946A (en) * 2010-03-17 2010-07-21 东南大学 Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement
CN101788403A (en) * 2010-03-17 2010-07-28 东南大学 Progressive method for identifying loose support cable based on strain monitoring during support settlement
CN101806670A (en) * 2010-04-06 2010-08-18 苏州大学 Method and device for detecting missing parts of forging press

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102706651A (en) * 2012-05-29 2012-10-03 东南大学 Damaged cable and support translation progressive identification method on basis of cable force monitoring during temperature variation
CN102706651B (en) * 2012-05-29 2015-09-30 东南大学 The damaged cable of cable force monitoring and support translation progressive identification method during temperature variation
CN102706670A (en) * 2012-05-30 2012-10-03 东南大学 Damaged cable and support generalized displacement identification method based on cable force monitoring of temperature change
CN102706670B (en) * 2012-05-30 2015-10-07 东南大学 The damaged cable of temperature variation cable force monitoring and generalized displacement of support recognition methods

Similar Documents

Publication Publication Date Title
CN100582719C (en) Health monitoring method for cable system in cable structure
CN101387600B (en) Cable system health monitoring method in cable structure based on mixed monitor
CN100595514C (en) Cable system health monitoring method based on strain monitoring cable structure
CN101673316B (en) Method for identifying supporting cables needing to adjust cable force based on hybrid monitoring
CN101587026B (en) Method of identifying support cable in need of regulating cable tension based on cable tension monitoring
CN101476989B (en) Progressive health monitoring method for cable system in cable structure based on mixed monitoring
CN102706659A (en) Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change
CN101476990B (en) Progressive health monitoring method for cable system in cable structure based on angle monitoring
CN101477340B (en) Progressive health monitoring method for cable system in cable structure based on cable force monitoring
CN105956216B (en) Correction method for finite element model greatly across steel bridge based on uniform temperature response monitor value
CN101334338B (en) Cable structure cable system healthy monitoring method based on space coordinates monitoring
CN103604642A (en) Problematic cable/concentrated load progressive recognition method on basis of generalized displacement cable power monitoring
CN103616222A (en) Hybrid-monitoring linear displacement progressive identification method for defective cable and concentrated loads
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
CN103604649A (en) Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring
CN103604623A (en) Progressive type damaged rope and concentrated load identification method based on angular displacement strain monitoring
CN103604628A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement angle monitoring
CN103868731A (en) Generalized displacement strain monitoring damaged cable load identification method
CN101504333B (en) Progressive health monitoring method for cable structure median cord system based on strain monitoring
CN102735459A (en) Progressive generalized displacement identification method for problem cables and support based on temperature change tension force monitoring
CN101788403B (en) Progressive method for identifying loose support cable based on strain monitoring during support settlement
CN103604644A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement strain monitoring
CN101782946B (en) Progressive type method for identifying loose supporting ropes based on space coordinate monitoring during support settlement
CN103616231A (en) Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20120118

C12 Rejection of a patent application after its publication