CN105067323A - Method for progressively recognizing damaged cable generalized displacement of load based on streamlined angle monitoring process - Google Patents

Abstract

The invention provides a method for progressively recognizing the damaged cable generalized displacement of a load based on the streamlined angle monitoring process. According to the method, based on the angle monitoring process, a mechanical calculation benchmark model of a cable structure is established. After that, a numerical value variation matrix of a monitored quantity with unit damage is obtained through calculation. Finally, according to the approximately linear relationships between the current numeric vector of the monitored quantity and the current initial numeric vector of the monitored quantity, the numerical value variation matrix of the monitored quantity with unit damage and the current nominal damage vector of an unknown to-be-evaluated object, the non-inferior solution of the current nominal damage vector of the to-be-evaluated object is figured out. In this way, the health status of a core to-be-evaluated object can be recognized.

Description

Simplify angle monitor load damaged cable generalized displacement progressive recognition method

Technical field

Cable-stayed bridge, suspension bridge, the structures such as truss structure have a common ground, it is exactly that they have many parts for bearing tensile load, such as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of the class formation is with rope, cable is subjected only to the rod member of tensile load for supporting member, for convenience, such structure representation is " Cable Structure " by this method, and by all carrying ropes of Cable Structure, carry cable, and all rod members (also known as two power rod members) for being subjected only to axial tension or axial compression load, it is collectively referred to as convenience " cable system ", in this method carrying rope is censured with " support cable " this noun, carrying cable and the rod member for being subjected only to axial tension or axial compression load, sometimes referred to simply as " rope ", so when " rope " this word is used below, two power rod members are just actually referred to truss structure.During structure military service, the correct identification to support cable or the health status of cable system is related to the safety of whole Cable Structure.When support cable is impaired, during Cable Structure military service, generalized displacement may occur for Cable Structure bearing, the load that Cable Structure is born may also be while change, even if the health status of actually Cable Structure does not change, the load that Cable Structure is born may also individually change, in this complex condition, this method recognizes generalized displacement of support and damaged cable, category engineering structure health monitoring field based on angle monitor (monitored angle is referred to as " monitored amount " by this method).

Background technology

The influence to Cable Structure health status recognition result of load change, Cable Structure generalized displacement of support is rejected, so that the change of the health status of identification structure exactly, the problem of be in the urgent need to address at present；Reject load change, Cable Structure health status and change influence to Cable Structure generalized displacement of support recognition result, so as to recognize Cable Structure generalized displacement of support exactly, the problem of be also in the urgent need to address at present；This method discloses a kind of effective, cheap method for solving the two problems.

The content of the invention

Technical problem：This method discloses a kind of method, under conditions of cost is lower, realizes two kinds of functions, is respectively, and one, the influence to Cable Structure health status recognition result of generalized displacement of support and load change is rejected, so as to identify the health status of support cable exactly；2nd, this method can also reject the influence of load change and the change of Cable Structure health status to Cable Structure generalized displacement of support recognition result, so as to recognize Cable Structure generalized displacement of support exactly.

Technical scheme：This method is made up of two parts.It is respectively：First, evaluation object and monitored amount (being referred to as monitored target) are determined；2nd, knowledge base and the method for parameter, knowledge based storehouse (containing parameter) and the monitored structural health conditions appraisal procedure measured of actual measurement set up needed for cable structure health monitoring system；3rd, the software and hardware part of health monitoring systems.

In the method, coordinate of the bearing on the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing space coordinate ", space coordinate of the bearing on X, Y, Z axis can also be said to be, bearing is referred to as in space coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the space coordinate of some axle with a space coordinate subscale of bearing on the concrete numerical value of the space coordinate of some axle；Angular coordinate of the bearing on X, Y, Z axis is censured with " bearing angular coordinate ", bearing is referred to as in angular coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the angular coordinate of some axle with an angular coordinate subscale of bearing on the concrete numerical value of the angular coordinate of some axle；Censure that bearing angular coordinate and bearing space coordinate are all with " bearing generalized coordinates ", in this method also with a generalized coordinates subscale of bearing up to bearing on the space coordinate of axle or the concrete numerical value of angular coordinate；Bearing is referred to as bearing displacement of the lines on the change of the coordinate of X, Y, Z axis, it may also be said to which the change of bearing space coordinate is referred to as the concrete numerical value in bearing displacement of the lines, this method also with the translational component expression bearing of bearing on the displacement of the lines of some axle；Bearing is referred to as the concrete numerical value in angular displacement of support, this method also with the angular displacement component expression bearing of bearing on the angular displacement of some axle on the change of the angular coordinate of X, Y, Z axis；Generalized displacement of support censures bearing displacement of the lines and angular displacement of support is all, and also bearing is reached on the displacement of the lines of some axle or the concrete numerical value of angular displacement with a generalized displacement subscale of bearing in this method；Bearing displacement of the lines is alternatively referred to as translational displacement, and support settlement is the component of bearing displacement of the lines or translational displacement in gravity direction.

The Part I of this method：Determine evaluation object and monitored amount.

The external force that object, structure are born can be described as load, and load includes face load and volume load.Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load.Volume load is the continuously distributed load in interior of articles each point, the deadweight of such as object and inertia force.

Concentrfated load is divided into two kinds of concentrated force and concentrated couple, in a coordinate system, for example in Descartes's rectangular coordinate system, one concentrated force can resolve into three components, likewise, a concentrated couple can also resolve into three components, if load is actually concentrfated load, force component or a concentrated couple component is concentrated to be referred to as a load by one in the method, the now change of load is embodied as a change for concentrating force component or a concentrated couple component.

Distributed load is divided into line distributed load and EDS maps load, the size of the zone of action and distributed load of the description of distributed load at least including distributed load, the size of distributed load is expressed with distribution intensity, being distributed intensity, to express, (such as two distributed loads are all uniform with distribution characteristics (such as uniform, SIN function equal distribution feature) and amplitude, but its amplitude is different, and the concept of amplitude can be illustrated by taking well-distributed pressure as an example：Same structure bears two different well-distributed pressures, and two distributed loads are all uniform loads, but the amplitude of a distributed load is 10MPa, and the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, one distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is not all identical, the component of this several distributed load is so regarded as same amount of independent distributed load in the method, now a load just represents the component of a distributed load, the amplitude changing ratio identical component of wherein distribution intensity can also be synthesized into a distributed load or is load.

Volume load is the continuously distributed load in interior of articles each point, deadweight and inertia force such as object, the size of the zone of action and volume load of the description of volume load at least including volume load, the size of volume load is expressed with distribution intensity, with distribution characteristics (such as uniform, linear function equal distribution feature) and amplitude, to express, (it is all uniform that lotuses are storaged in such as two individuals to distribution intensity, but its amplitude is different, and the concept of amplitude is illustrated exemplified by can conducting oneself with dignity：The material of two parts of same structure is different, therefore density is different, so while the volume load suffered by the two parts is all uniform, but the amplitude of the volume load suffered by a part is probably 10kN/m³, the amplitude of the volume load suffered by another part is 50kN/m³).If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes.In a coordinate system, one individual stowage lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles on coordinate system, that is, volume load can resolve into three components in Descartes's rectangular coordinate system), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is not all identical, the component that lotus so is storaged in this several body in the method regards same amount of independent load as, the amplitude changing ratio identical volume sharing part of the load of wherein distribution intensity can also be synthesized to an individual stowage lotus or it is load.

When load is embodied as concentrfated load, in the method, " load unit change " actually refers to " unit change of concentrfated load ", similar, " load change " is referred specifically to " change of the size of concentrfated load ", " load change amount " is referred specifically to " variable quantity of the size of concentrfated load ", " load change degree " is referred specifically to " intensity of variation of the size of concentrfated load ", " the actual change amount of load " refers to " the actual change amount of the size of concentrfated load ", " load changed " refers to " concentrfated load that size changes ", briefly, now " so-and-so load so-and-so change " refers to " size of so-and-so concentrfated load so-and-so change ".

When load is embodied as distributed load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, similar, " load change " is referred specifically to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " is referred specifically to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " is referred specifically to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " is referred specifically to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that the amplitude of distribution intensity changes ", briefly, now " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so distributed load so-and-so change ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.

When load is embodied as volume load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of volume load ", similar, " load change " refers to " change of the amplitude of the distribution intensity of volume load ", " load change amount " refers to " variable quantity of the amplitude of the distribution intensity of volume load ", " load change degree " refers to " intensity of variation of the amplitude of the distribution intensity of volume load ", " the actual change amount of load " refers to " the actual change amount of the amplitude of the distribution intensity of volume load ", " load changed " refers to " the volume load that the amplitude of distribution intensity changes ", briefly, " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.

First confirm that the quantity for the load that the possibility that Cable Structure is born changes.The characteristics of load born according to Cable Structure, confirmation wherein " is possible to the load changed ", or all load is considered as " being possible to the load changed ", if the shared JZW load that may be changed, that is, have JZW secondary evaluation objects.

If the quantity of the generalized displacement of support component of Cable Structure, the quantity of the support cable of Cable Structure and JZW " being possible to the load changed " quantity sum are N, that is, have N number of evaluation object.Evaluation object serial number is given, the numbering will be used to generate vector sum matrix in subsequent step.

If evaluated support cable and the quantity sum of generalized displacement of support component are that P, the i.e. quantity of core evaluation object are P, if the quantity of evaluated generalized displacement of support component is Z, if the quantity of evaluated support cable is M₁。

" the monitored angle-data of the whole of structure " H angle coordinate components of L specified straight lines, each specified straight lines of K specified point, excessively each specified point in structure are described, the change for changing exactly all specified points, all specified straight lines all angle coordinate components specified of structural point.M (M=K × L × H) individual angle coordinate component measurement value or calculated value is had every time to characterize the angle information of structure.

The monitored amount of summary, whole Cable Structure has M monitored amounts, and M should be more than the quantity of core evaluation object, and M is less than the quantity of evaluation object.

For convenience, " monitored all parameters of Cable Structure " are referred to as in the method " monitored amount ".To M monitored amount serial numbers, the numbering will be used to generate vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 ..., M with variable j.

The Part II of this method：Knowledge base and the method for parameter, knowledge based storehouse (containing parameter) and the monitored structural health conditions appraisal procedure measured of actual measurement set up needed for cable structure health monitoring system.It can carry out successively as follows, to obtain the health state evaluation of more accurately evaluation object.

The first step：Set up initial mechanical calculating benchmark model A_oWhen Cable Structure is completed, or before health monitoring systems are set up, the initial value for all monitored amounts for obtaining Cable Structure is calculated using conventional method direct measurement, the physical parameter of various materials used in Cable Structure and mechanical property parameters (such as modulus of elasticity, Poisson's ratio) are obtained using conventional method (consult reference materials or survey) simultaneously, while obtaining the Actual measurement data of Cable Structure using conventional method Actual measurement.The Actual measurement data of Cable Structure are the data of the health status that can express rope including the Non-destructive Testing Data of support cable first, the measured data of the Actual measurement data of Cable Structure still including the initial generalized displacement measurement data of Cable Structure bearing, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinates data, Cable Structure modal data, structural strain data, structural point measurement data, structure space measurement of coordinates data, load data.The initial geometric data of Cable Structure can be spatial data of the spatial data plus a series of point in structure of the end points of all ropes, it is therefore intended that the geometric properties of Cable Structure are determined according to these coordinate datas.The initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model A_oWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs.For cable-stayed bridge, if initial geometric data, which can be the spatial data of the end points of all ropes, adds the spatial data done on bridge two ends, here it is so-called bridge type data.The variable quantity of " being possible to the load changed " is setting up initial mechanical calculating benchmark model A_oWhen all 0, that is to say, that the variable quantity of " being possible to the load changed " for identifying below is relative to setting up initial mechanical calculating benchmark model A_oWhen the variable quantity of corresponding load that is born of structure.Non-destructive Testing Data using support cable etc. can express the initial generalized displacement measurement data of the data of the health status of support cable, Cable Structure bearing and the variable quantity data of " being possible to the load changed " set up evaluation object initial damage vector d_o(as shown in formula (1)), uses d_oRepresent Cable Structure (with initial mechanical calculating benchmark model A_oRepresent) evaluation object initial health.If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0.Vectorial d_oIn each element numerical value related to the variable quantity of load take 0.If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial d_oIn each element numerical value related to Cable Structure generalized displacement of support take 0.Using the physical and mechanical properties parameter of various materials used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure, initial mechanical calculating benchmark model A is set up using mechanics method (such as FInite Element)_o。

d_o=[d_o1d_o2···d_ok···d_oN]^T(1)

D in formula (1)_ok(k=1,2,3 ..., N) represent initial mechanical calculating benchmark model A_oIn k-th of evaluation object original state, if the evaluation object is a rope (or pull bar) in cable system, then d_okRepresent its initial damage, d_okFor 0 when represent not damaged, represent that the rope thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity；If the evaluation object is a generalized displacement component of a bearing, then d_okRepresent its initial displacement numerical value；If the evaluation object is one " load that may be changed ", then d_okRepresent its initial value, d_okFor 0, that is to say, that the variable quantity of " being possible to the load changed " for identifying below is relative to setting up initial mechanical calculating benchmark model A_oWhen the variable quantity of corresponding load that is born of structure.Subscript T represents the transposition (rear same) of vector in formula.

The initial value of all monitored amounts of obtained Cable Structure, the monitored amount initial value vector C of composition are calculated using conventional method direct measurement_o(see formula (2)).It is required that obtaining A_oWhile obtain C_o, it is monitored amount initial value vector C_oRepresent to correspond to A_o" monitored amount " concrete numerical value.Because subject to the foregoing, the monitored amount obtained by the calculating benchmark model calculating based on Cable Structure in narration below, will represent the calculated value and measured value reliably close to the measured data of initial monitored amount with same symbol.

C_o=[C_o1C_o2···C_oj···C_oM]^T(2)

C in formula (2)_oj(j=1,2,3 ... ..., M) are the primary quantities of j-th of monitored amount in Cable Structure, and the component corresponds to specific j-th monitored amount according to coding rule.Vectorial C_oIt is to be formed by M monitored amounts according to certain order arrangement, this, which is put in order, has no particular/special requirement, only requires that all associated vectors also arrange data in this order below.

No matter which kind of method to obtain initial mechanical calculating benchmark model A with_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, error typically cannot be greater than 5%.So can utility A_oSuo Li under analog case obtained by calculating calculates that data, strain calculation data, Cable Structure shapometer count evidence and displacement meter counts evidence, Cable Structure angle-data, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A_oThe health status of middle evaluation object evaluation object initial damage vector d_oRepresent.Due to based on A_oThe initial value (actual measurement is obtained) for obtaining the evaluations of all monitored amounts closely all monitored amounts is calculated, so A can also be used in_oOn the basis of, carry out Mechanics Calculation obtain, A_oEach monitored amount the monitored amount initial value vector C of evaluation composition_o。d_oIt is A_oParameter, it may also be said to C_oBy A_oMechanics Calculation result composition.

Second step：Circulation starts.When circulation starts each time, it is necessary first to set up or set up the current initial damage vector d of evaluation object when this circulation startsⁱ _o(i=1,2,3 ...), the current initial mechanical calculating benchmark model A for setting up Cable Structureⁱ _o(such as finite element benchmark model, the A in circulating each timeⁱ _oIt is to constantly update).Alphabetical i is in addition to the place for clearly indicating that number of steps, and letter i only represents that cycle-index, i.e. ith are circulated in the method.

The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o(as shown in formula (3)), uses dⁱ _oRepresent that Cable Structure is (with current initial mechanical calculating benchmark model A during this time circulation beginningⁱ _oRepresent) evaluation object health status.

In formula (3)When representing that ith circulation starts, current initial mechanical calculating benchmark model Aⁱ _oIn k-th of evaluation object original state, if the evaluation object is a rope (or pull bar) in cable system, then dⁱ _okRepresent its initial damage, dⁱ _okFor 0 when represent not damaged, represent that the rope thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity, if the evaluation object is one " load that may be changed ", then dⁱ _okRepresent it relative to setting up initial mechanical calculating benchmark model A_oWhen the variable quantity of corresponding load that is born of structure, if the evaluation object is a generalized displacement component of a bearing, then dⁱ _okRepresent its initial displacement numerical value.

Set up and update dⁱ _oMethod it is as follows：

When circulation starts for the first time, setting up the current initial damage vector of evaluation object, (foundation formula (3) is designated as d¹ _o) when, d¹ _oIt is equal to d_o.The current initial damage vector d of evaluation object that i-th (i=2,3,4,5,6 ...) secondary circulations need when startingⁱ _o, it is to circulate to terminate preceding calculating acquisition in preceding once (i.e. the i-th -1 time, i=2,3,4,5,6 ...), specific method is described below.

I-th (i=1,2,3,4,5,6 ...) secondary circulations need the Mechanics Calculation benchmark model for the Mechanics Calculation benchmark model or Cable Structure set up set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o。

Set up, update Aⁱ _oMethod it is as follows：

The Mechanics Calculation benchmark model for the Cable Structure set up when circulation starts for the first time is designated as A¹ _o, A¹ _oEqual to A_o.The A in circulating each timeⁱ _oIt is to constantly update, specific method is described below；At the end of circulating each time, A is updatedⁱ _oThe Mechanics Calculation benchmark model of Cable Structure required during beginning is circulated next time, and specific method is described below.

This method " monitored current initial value vector C of amountⁱ _o" (i=1,2,3 ...) represent ith (i=1,2,3,4,5,6 ...) circulation start when all monitored amounts specified initial value (referring to formula (4)), Cⁱ _oBe referred to as " the current initial value vector of the monitored amount of ith circulation ".

In formula (2)J-th of monitored amount when being ith circulation beginning, in Cable Structure.Vectorial Cⁱ _oIt is to be formed by previously defined M monitored amounts according to certain order arrangement, this, which is put in order, has no particular/special requirement, only requires that all associated vectors also arrange data in this order below.

Setting up model Aⁱ _oWhile set up " the monitored current initial value vector C of amountⁱ _o", it is monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value, Cⁱ _oElement and C_oElement correspond, represent respectively all monitored amounts Cable Structure be in Aⁱ _oAnd A_oConcrete numerical value during two states.

Set up and update Cⁱ _oSpecific method it is as follows：

When circulation starts for the first time, C¹ _o(i=1, Cⁱ _oIt is embodied as C¹ _o) it is equal to C_o；Ith circulation " the monitored current initial value vector C of amount that i-th (i=2,3,4,5,6 ...) secondary circulations need when startingⁱ _o", it is to circulate to terminate preceding calculating acquisition in preceding once (i.e. the i-th -1 time, i=2,3,4,5,6 ...), specific method is described below.In ith (i=1,2,3,4,5,6 ...) circulation, " the monitored current initial value vector C of amountⁱ _o" it is to constantly update, specific method is described below.Due to according to model Aⁱ _oThe initial value of the monitored amount of gained is calculated reliably close to corresponding measured value, in narration below, the calculated value composition of vector and measured value composition of vector will be represented with same symbol.

dⁱ _oIt is Aⁱ _oCharacterisitic parameter, Cⁱ _oIt is Aⁱ _oIn dⁱ _oUnder the conditions of Mechanics Calculation result composition.

3rd step：During Cable Structure military service, in circulating each time, in other words in i-th (i=1,2,3,4,5,6 ...) secondary circulation, in known Aⁱ _o、Cⁱ _oAnd dⁱ _oAfterwards, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositionⁱ。CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time.

4th step：" unit damage monitored numerical quantity transformation matrices " and " evaluation object unit change vector " must be first set up when circulating each time, " the unit damage monitored numerical quantity transformation matrices " that ith circulation is set up are designated as Δ Cⁱ(i=1,2,3 ...)." the evaluation object unit change vector " that ith circulation is set up is designated as Dⁱ _u.The Δ C in circulating each timeⁱAnd Dⁱ _uNeed according to circumstances to constantly update, that is, updating current initial mechanical calculating benchmark model Aⁱ _oWith monitored amount current initial value vector Cⁱ _oAfterwards, unit damage monitored numerical quantity transformation matrices Δ C is updatedⁱWith evaluation object unit change vector Dⁱ _u。

Unit damage monitored numerical quantity transformation matrices Δ C is first set up in the steps below when circulation starts each timeⁱWith evaluation object unit change vector Dⁱ _u, detailed process arranges as follows：

In the current initial mechanical calculating benchmark model A of Cable Structureⁱ _oOn the basis of calculated several times, calculation times are numerically equal to the quantity of all ropes.Calculate each time and assume that (original load change amount can be 0 to only one of which evaluation object in original damage or generalized displacement or load change amount, can not also be to be further added by unit generalized displacement, unit damage or load unit change on the basis of 0), specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage (such as taking 5%, 10%, 20% or 30% equivalent damage to be damaged for unit), if the evaluation object is a load, it is assumed that the load is in vectorial dⁱ _oLoad unit is further added by the basis of the existing variable quantity of the load represented to change (if the load is distributed load, and the distributed load is line distributed load, load unit change can take 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. to change for unit；If the load is distributed load, and the distributed load is EDS maps load, and load unit change can take 1MPa, 2MPa, 3MPa or 1kNm/m²、2kNm/m²、3kNm/m²Change Deng for unit；If the load is concentrfated load, and the concentrfated load is couple, and load unit change can take 1kNm, 2kNm, 3kNm etc. to change for unit；If the load is concentrfated load, and the concentrfated load is concentrated force, and load unit change can take 1kN, 2kN, 3kN etc. to change for unit；If the load is volume load, load unit change can take 1kN/m³、2kN/m³、3kN/m³Deng for unit change), if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing in the direction of displacement in vectorial dⁱ _oUnit generalized displacement occurs again (if such as evaluation object is the translational component in the x directions of a bearing on the basis of the existing generalized displacement of the bearing represented, it is assumed that the bearing has unit displacement of the lines in x directions, for example take 1mm, if the evaluation object is the angular displacement component around x-axis of a bearing, it is assumed that the bearing has unit angular displacement around x-axis, for example, take ten a ten thousandth radians).For convenience of calculating, can all be that this time is circulated structural health conditions when starting as being complete health when setting increase unit generalized displacement, unit damage or load unit change in circulating each time, and on this basis the generalized displacement of setting increase unit, unit damage or load unit change (in subsequent step, calculate, the generalized displacement numerical value of evaluation object, damage numerical value or load change amount --- be referred to as nominal fatigue dⁱ _c(i=1,2,3 ...), be relative to by this time circulate it is when starting, by the health status of evaluation object as being complete health, it is therefore necessary to the nominal fatigue calculated is converted into true damage according to formula given hereinlater).With there is unit generalized displacement in the calculating each time of one cycle, unit damage or the evaluation object of load unit change are different from unit generalized displacement occur in other calculating, unit damage or the evaluation object of load unit change, and unit generalized displacement is suppose there is each time, the unit generalized displacement of unit damage or the evaluation object of load unit change, unit damage value or load unit change numerical value can be differently configured from the unit generalized displacement of other evaluation objects, unit damage value or load unit change numerical value, with " evaluation object unit change vector Dⁱ _u" (as shown in formula (5)) record the unit generalized displacement of hypothesis of all evaluation objects in each circulation, unit damage or load unit change, and D is designated as when circulating for the first time¹ _uCalculate each time and all calculate Cable Structure, current calculated values in the monitored amount of the M above specified using mechanics method (such as FInite Element), the current calculated value of M monitored amounts obtained by calculating each time constitutes one " being monitored amount calculation current vector " (when assuming that k-th of evaluation object has unit damage, can represent all the M specified the monitored monitored amount calculation current vector C measured with formula (6)ⁱ _tk)；Obtained monitored amount calculation current vector is calculated each time subtracts the monitored current initial value vector C of amountⁱ _o, " the numerical value change vector of monitored amount " (to have the position of the rope of unit damage or numbering etc. to be mark) (when k-th of evaluation object has unit damage, uses δ C under the conditions of gained vector is exactly thisⁱ _kRepresent the numerical value change vector of monitored amount, δ Cⁱ _kDefinition see formula (7), formula (8) and formula (9), formula (7) is that formula (6) is subtracted after formula (4) again divided by vector Dⁱ _uK-th of element Dⁱ _ukGained), it is monitored the numerical value change vector δ C of amountⁱ _kEach element representation due to suppose there is unit generalized displacement, the unit generalized displacement of that evaluation object (such as k-th evaluation object) of unit damage or load unit change, unit damage or load unit change (such as D when calculatingⁱ _uk), caused by monitored amount corresponding to the element numerical value knots modification relative to the unit generalized displacement of hypothesis, unit damage or load unit change numerical value Dⁱ _ukRate of change；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has M element, is made up of " the unit damage monitored numerical quantity transformation matrices Δ C of M × N number of element successively this N number of " numerical value change vector of monitored amount "ⁱ" (M rows N row), each vectorial δ Cⁱ _k(k=1,2,3 ..., N) it is matrix Δ CⁱA row, Δ CⁱDefinition such as formula (10) shown in.

Evaluation object unit change vector D in formula (5)ⁱ _uElement Dⁱ _uk(i=1,2,3 ...；K=1,2,3 ... ..., N) represent that unit generalized displacement, unit damage or the load unit of k-th of the evaluation object assumed in ith circulation change numerical value, vectorial Dⁱ _uIn the numerical value of each element can be the same or different.

Elements C in formula (6)ⁱ _tkj(i=1,2,3 ...；K=1,2,3 ..., N；J=1,2,3 ... ..., M) ith circulation is represented due to when k-th of evaluation object has unit generalized displacement, unit damage or load unit change, the calculating current value for the monitored amount specified according to j-th corresponding to coding rule.

The subscript i (i=1,2,3 ...) respectively measured in formula (7) represents ith circulation, subscript k (k=1,2,3, ..., N) represent D in k-th of evaluation object increase unit generalized displacement, unit damage or load unit change, formulaⁱ _ukIt is vectorial Dⁱ _uIn k-th of element.Vectorial δ Cⁱ _kDefinition such as formula (7) and formula (8) shown in, δ Cⁱ _kJth (j=1,2,3 ..., M) individual element δ Cⁱ _kj(definition is as shown in formula (9)) is represented in ith circulation, sets up matrix Δ CⁱWhen, it is assumed that k-th of evaluation object calculates j-th of monitored amount of gained knots modification when having unit generalized displacement, unit damage or load unit change changes D relative to the unit generalized displacement of hypothesis, unit damage or load unitⁱ _ukRate of change.

Vector δ C in formula (10)ⁱ _k(i=1,2,3 ... ..., k=1,2,3 ... ..., N) represent in ith circulation, due to k-th of evaluation object increase unit generalized displacement, unit damage or load unit change Dⁱ _ukCaused by, the change of the relative values of all monitored amounts.Matrix Δ CⁱRow (subscript k) coding rule and above vector dⁱ _oElement subscript k coding rule it is identical.

5th step：Recognize the current health state of Cable Structure.Detailed process is as follows.

In i-th (i=1,2,3 ...) secondary circulation, " the monitored amount current value vector C obtained in second step actual measurement is utilizedⁱ" " monitored to measure current initial value vector C togetherⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between linear approximate relationship, as shown in formula (11) or formula (12).

Amount current value vector C is monitored in formula (11) and formula (12)ⁱBe defined similarly as the monitored current initial value vector C of amountⁱ _oDefinition, see formula (13)；The current nominal fatigue vector d of evaluation objectⁱ _cDefinition see formula (14).

Elements C in formula (13)ⁱ _j(i=1,2,3 ...；J=1,2,3 ..., M) it is Cable Structure, the current value for the monitored amount for being j according to the numbering corresponding to coding rule when ith is circulated.

D in formula (14)ⁱ _ck(i=1,2,3 ...；K=1,2,3 ... ..., N) be ith circulation in k-th of evaluation object current nominal generalized displacement numerical value, current nominal fatigue or current nominal load changing value, vectorial dⁱ _cElement subscript k coding rule and matrix Δ CⁱRow coding rule it is identical.

When support cable actual damage is less big, because Cable Structure material remains at the linear elasticity stage, the deformation of Cable Structure is also smaller, and such a linear relationship represented by formula (11) or formula (12) is smaller with the error of actual conditions, and error can use error vector eⁱ(formula (15)) are defined, the error of linear relationship shown in expression (11) or formula (12).

Abs () is the function that takes absolute value in formula (15), and each the vectorial element tried to achieve in bracket is taken absolute value.

Because the linear relationship represented by formula (11) or formula (12) has certain error, therefore obtain the current nominal fatigue vector d of evaluation objectⁱ _cAcceptable solution (i.e. with reasonable error, but the position of damaged cable can be determined relatively accurately from cable system and its degree of injury, load change numerical value can also be relatively accurately determined) turn into a rational solution, this method can be expressed with formula (16).

Abs () is the function that takes absolute value in formula (16), vectorial gⁱDescription deviates the legitimate skew of ideal linearity relation (formula (11) or formula (12)), is defined by formula (17).

G in formula (17)ⁱ _j(i=1,2,3 ...；J=1,2,3 ..., M) describe ith circulation in deviate formula (11) or formula (12) shown in ideal linearity relation maximum allowable offset.Vectorial gⁱThe error vector e that can be defined according to formula (15)ⁱTentative calculation is selected.

In the monitored current initial value vector C of amountⁱ _o, unit damage monitored numerical quantity transformation matrices Δ CⁱWith monitored amount current value vector CⁱWhen known, it is possible to use suitable algorithm (such as multi-objective optimization algorithm) solves formula (16), obtains the current nominal fatigue vector d of evaluation objectⁱ _cAcceptable solution, the currently practical injury vector d of cable systemⁱThe element of (formula (18) is shown in definition) can be calculated according to formula (19) and obtained, so as to by dⁱIt is determined that determining the health status of evaluation object.

D in formula (18)ⁱ _k(i=1,2,3 ...；K=1,2,3 ... ..., N) the currently practical health status of k-th of evaluation object in ith circulation is represented, its definition is shown in formula (19).

D in formula (19)ⁱ _ok(i=1,2,3,4 ...；K=1,2,3 ..., N) it is the current initial damage vector d of evaluation objectⁱ _oK-th of element, dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element.dⁱ _k(i=1,2,3 ...；K=1,2,3 ... ..., N) represent ith circulation in k-th of evaluation object currently practical health status, if the evaluation object is a support cable (or pull bar) in cable system, then dⁱ _kRepresent its currently practical damage, dⁱ _kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity；If the evaluation object is a load, then dⁱ _kRepresent the currently practical change numerical value of its corresponding load；If the evaluation object is a generalized displacement component of a bearing, then dⁱ _kRepresent its currently practical generalized displacement numerical value；Vectorial dⁱElement coding rule and formula (1) in vector d_oElement coding rule it is identical.

So far this method realizes accurately identifying for the health status of core evaluation object in a kind of effective, cheap method.It may deviate the recognition result of the health status of secondary evaluation object that exact value is more, the health status of correct identification core evaluation object is required nothing more than in the method.

6th step：Judge whether to terminate this (ith) circulation, be that (i.e. i+1 time, i=1,2,3,4 ...) circulation prepares Mechanics Calculation benchmark model and necessary vector next time if it is, completing the round-off work before this circulation terminates.Detailed process is as follows：

Current nominal fatigue vector d is tried to achieve in this (ith) circulationⁱ _cAfterwards, first, mark vector B is set up according to formula (20)ⁱ, formula (21) gives mark vector BⁱK-th of element definition；If mark vector BⁱElement be all 0, then return to the 3rd step and proceed this circulation, continue health monitoring and calculating to Cable Structure；If mark vector BⁱElement be not all 0, then complete after subsequent step, into circulating next time.

So-called subsequent step is：First, the initial damage vector d obtained next time needed for (i.e. i+1 time, i=1,2,3,4 ...) circulation is calculated according to formula (22)ⁱ⁺¹ _oEach element dⁱ⁺¹ _ok；Second, in Mechanics Calculation benchmark model A_oOn the basis of, make A_oIn evaluation object health status be dⁱ⁺¹ _oRather than be d_oAfterwards, the current initial mechanical calculating benchmark mould A next time needed for (i.e. i+1 time, i=1,2,3,4 ...) circulation has thus been obtainedⁱ⁺¹ _o, to Aⁱ⁺¹ _oMechanics Calculation is carried out to obtain corresponding to Aⁱ⁺¹ _oAll monitored amounts, current concrete numerical values, these concrete numerical values constitute the current initial value vector C of required monitored amount of (i.e. i+1 time, i=1,2,3,4 ...) circulation next timeⁱ⁺¹ _o。

Mark vector B in formula (20)ⁱSubscript i represent ith circulate, its element Bⁱ _kThe subscript k of (k=1,2,3 ..., N) represents the health status feature of k-th of evaluation object, can only take 0 and 1 two amount, and specific value rule is shown in formula (21).

Element B in formula (21)ⁱ _kIt is mark vector BⁱK-th of element, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element (see formula (5)), dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element (see formula (14)), they all represent the relevant information of k-th of evaluation object.

D in formula (22)ⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element (see formula (5)), dⁱ _okIt is the current initial damage vector d of evaluation objectⁱ _oK-th of element (see formula (3)).

The Part III of this method：The software and hardware part of health monitoring systems.

Hardware components include monitoring system (including monitored amount monitoring system), signal picker and computer etc..It is required that monitoring each monitored amount in real time.

Software section should can complete the process set by this method, that is, complete it is required in this method, can with computer implemented monitoring, record, control, storage, calculate, notify, the function such as alarm.

This method is specifically included：

A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable；The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation；This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, i.e. evaluation object is N；The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step；This method represents this numbering, k=1,2,3 ..., N with variable k；This method refers exclusively to evaluated support cable and generalized displacement of support component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of generalized displacement of support component are P, i.e. the quantity of core evaluation object is P, and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；If having M in cable system₁Root support cable；It is determined that the measured point specified, gives all specified points numbering；It is determined that the measured straight line of excessively each measurement point, to all measured straight line numberings specified；The measured angle coordinate component of each measured straight line is determined, is numbered to all measured angle coordinate components；Above-mentioned numbering will be used to generate vector sum matrix in subsequent step；" the monitored angle-data of the whole of Cable Structure " is made up of above-mentioned all measured angle coordinate components；For convenience, " the monitored angle-data of Cable Structure " is referred to as in the method " monitored amount "；The quantity sum of all monitored amounts is designated as M, and M should be more than the quantity of core evaluation object, and M is less than the quantity of evaluation object；The external force that object, structure are born can be described as load, and load includes face load and volume load；Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load；Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force；Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component；Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant；Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load；Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes；Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method；

B. survey or consult reference materials and obtain the physical and mechanical properties parameter of various materials used in Cable Structure；

C. while surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, direct measurement calculates the measured data for obtaining initial Cable Structure, and the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, and the data of the health status that can express support cable now are referred to as support cable initial health data；The monitored amount initial value vector C of initial value composition of all monitored amounts_o, it is monitored amount initial value vector C_oCoding rule and M monitored amounts coding rules it is identical；Evaluation object initial damage vector d is set up using support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement data_o, vectorial d_oRepresent with initial mechanical calculating benchmark model A_oThe initial health of the evaluation object of the Cable Structure of expression；Evaluation object initial damage vector d_oElement number be equal to N, d_oElement and evaluation object be one-to-one relationship, vectorial d_oElement coding rule it is identical with the coding rule of evaluation object；If d_oThe corresponding evaluation object of some element be a support cable in cable system, then d_oThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion；If d_oThe corresponding evaluation object of some element be some bearing some generalized displacement component, then d_oThe element numerical value represent this bearing the generalized displacement component initial value；If d_oThe corresponding evaluation object of some element be to take d in some load, this method_oThe element numerical value be 0, the initial value for representing the change of this load is 0；If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial d_oIn each element numerical value related to Cable Structure generalized displacement of support take 0；If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0；Initial Cable Structure bearing generalized coordinates data refer to the bearing generalized coordinates data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model A_oWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs；

D. all Cable Structure data that the physical and mechanical properties parameter and preceding step of various materials according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure are obtained, set up the initial mechanical calculating benchmark model A of Cable Structure_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot be greater than 5%；Corresponding to A_oEvaluation object health status evaluation object initial damage vector d_oRepresent；Corresponding to A_oAll monitored amounts initial value with monitored amount initial value vector C_oRepresent；d_oIt is A_oParameter, by A_oThe obtained initial value of all monitored amounts of Mechanics Calculation result and C_oThe initial value of all monitored amounts represented is identical, therefore alternatively C_oBy A_oMechanics Calculation result composition, A in the method_o、C_oAnd d_oIt is constant；

E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o；The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；dⁱ _oIt is Aⁱ _oCharacterisitic parameter, Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o；

F. enter from here and the circulation walked to pth is walked by f；

G. during structure military service, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time；

H. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u；

H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4ⁱAnd Dⁱ _u；At other moment, when in step g to Aⁱ _oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CⁱAnd Dⁱ _uIf, not to A in step gⁱ _oIt is updated, then is directly transferred to step i here and carries out follow-up work；

H2. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating；According to the coding rule of evaluation object, calculated successively；Calculate each time and assume that only one of which evaluation object is further added by unit damage or unit generalized displacement or load unit change on the basis of original damage or generalized displacement or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing is further added by unit generalized displacement in the direction of displacement, if the evaluation object is a load, it is assumed that the load is further added by load unit change, use Dⁱ _ukThis increased unit damage or unit generalized displacement or load unit change are recorded, wherein k represents to increase the numbering of the evaluation object of unit damage or unit generalized displacement or load unit change, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uAn element, evaluation object unit change vector Dⁱ _uElement coding rule and vector d_oElement coding rule it is identical；The evaluation object that unit damage or unit generalized displacement or load unit change are further added by calculating each time is different from the evaluation object that unit damage or unit generalized displacement or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector；When assuming that k-th of evaluation object is further added by unit damage or unit generalized displacement or load unit change, C is usedⁱ _tkRepresent corresponding " monitored amount calculation current vector "；When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed；Cⁱ _tkDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _tkElement and C_oElement correspond；

H3. obtained vectorial C is calculated each timeⁱ _tkSubtract vectorial Cⁱ _oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or unit generalized displacement or load unit change numerical valueⁱ _k”；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount "；

H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "ⁱ”；Unit damage monitored numerical quantity transformation matrices Δ CⁱEach row correspond to a monitored amount unit change vector；Unit damage monitored numerical quantity transformation matrices Δ CⁱEvery a line correspond to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit generalized displacement or load unit change；Unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and vector d_oElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and M monitored amounts coding rules it is identical；

I. current nominal fatigue vector d is definedⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cElement numerical value represent the nominal fatigue degree or nominal generalized displacement or nominal load variable quantity of correspondence evaluation object, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship；

J. according to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1ⁱ _cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector dⁱ _c；

Formula 1

K. the currently practical injury vector d expressed using formula 2ⁱK-th of element dⁱ _kWith the current initial damage vector d of evaluation objectⁱ _oK-th of element dⁱ _okWith current nominal fatigue vector dⁱ _cK-th of element dⁱ _ckBetween relation, calculating obtain currently practical injury vector dⁱAll elements；

Formula 2

K=1,2,3 in formula 2 ..., N；dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kRepresent its currently practical damage, dⁱ _kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity；If the evaluation object is a load, then dⁱ _kRepresent the actual change amount of the load；If the evaluation object is a generalized displacement component of a bearing, then dⁱ _kRepresent its currently practical generalized displacement numerical value；So far this method, which is realized, rejects generalized displacement of support, load change influence, Cable Structure damaged cable identification, realize simultaneously and reject generalized displacement of support and the change influence of support cable health status, load change amount identification, also achieve and reject load change and the change influence of support cable health status, generalized displacement of support identification；So far this method realizes accurately identifying for the health status of core evaluation object in a kind of effective, cheap method；Recognition result deviation exact value to the health status of secondary evaluation object is more, therefore not accepts and believe, and the health status of correct identification core evaluation object is required nothing more than in the method；

L. current nominal fatigue vector d is tried to achieveⁱ _cAfterwards, mark vector B is set up according to formula 3ⁱ, formula 4 gives mark vector BⁱK-th of element definition；

Formula 3

Formula 4

Element B in formula 4ⁱ _kIt is mark vector BⁱK-th of element, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element, dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4；

If m. mark vector BⁱElement be all 0, then return to step f continue this circulation；If mark vector BⁱElement be not all 0, then enter next step, i.e. step n；

N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5ⁱ⁺¹ _oEach element；

Formula 5

D in formula 5ⁱ⁺¹ _okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationⁱ⁺¹ _oK-th of element, dⁱ _okThis, i.e. the current initial damage vector d of evaluation object of ith circulationⁱ _oK-th of element, Dⁱ _ukIt is the evaluation object unit change vector D of ith circulationⁱ _uK-th of element, Bⁱ _kIt is the mark vector B of ith circulationⁱK-th of element, k=1,2,3 ... ..., N in formula 5；

O. in initial mechanical calculating benchmark model A_oOn the basis of, the health status for making rope is dⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationⁱ⁺¹；Obtain Aⁱ⁺¹Afterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o；

P. step f is returned to, starts to circulate next time.

Beneficial effect：Structural healthy monitoring system is monitored on-line for a long time by using sensor to structural response first, (or offline) analysis online is carried out after acquisition Monitoring Data to it and obtains structural health conditions data, due to the complexity of structure, structural healthy monitoring system needs to use the equipment such as substantial amounts of sensor to carry out monitoring structural health conditions, therefore its cost is generally at a relatively high, therefore cost problem is a subject matter of limit structural health monitoring technique application.On the other hand, the correct identification of the health status of core evaluation object (such as suspension cable) is the indispensable part of the correct identification of structural health conditions, even its whole, and influence of the correct identification of the change of secondary evaluation object (load that such as structure is born) (such as by the change of the quality and quantity of the automobile of cable-stayed bridge) to the correct identification of the health status of Cable Structure be it is very little, it is even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is typically suitable, the quantity of secondary evaluation object is also frequently more than the quantity of core evaluation object, and the quantity of such evaluation object is often many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, the quantity of the monitored amount of conventional method requirement (being obtained using device measurings such as sensors) have to be larger than the quantity equal to evaluation object, when the secondary evaluation object changed quantity than it is larger when (practically always such), the quantity of the equipment such as the sensor required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, or even unacceptablely high.Inventor's research is found,In secondary evaluation object (such as normal load that structure is born,The normal load of structure refers to that the load that structure bearing is no more than the structure allowable load limited according to structure design book or structure completion book) change it is smaller when (be exactly that structure is only subjected only to normal load for load,Whether the load that structure is born is normal load,It can be determined by the observation of the methods such as naked eyes,If it find that the load that structure is born is not normal load,It is so artificial to remove,Remove after improper load,Structure is just solely subjected to normal load),Amplitude of variation of the amplitude of variation (this specification is called " secondary response ") of structural response caused by them much smaller than the structural response caused by the change (such as support cable is damaged) of core evaluation object (this specification is called " core response "),Secondary response responds total change that sum is structural response with core (this specification is called " global response "),Obvious core response occupies leading position in global response,Based on this,Even if inventor's research hair is chosen when now determining that monitored amount quantity is slightly larger than core evaluation object quantity,But much smaller than the numerical value of evaluation object quantity (this method is exactly so to do),Even if that is using equipment such as the relatively few many sensors of quantity,Still the state of health data of core evaluation object can accurately be obtained,Meet the core demand of structural health conditions monitoring,Therefore cost of the cost of the structural healthy monitoring system proposed by this method apparently than the structural healthy monitoring system required by conventional method is much lower,That is this method can realize the assessment of the health status to the core evaluation object of Cable Structure with the much lower condition of cost,This benefit is that can structural health monitoring technology be used is very important.

Embodiment

Illustrate what is be substantially merely exemplary below the embodiment of this method, and purpose is never to limit the application of this method or used.

The first step：First confirm that the quantity for the load that the possibility that Cable Structure is born changes.The characteristics of load born according to Cable Structure, confirmation wherein " is possible to the load changed ", or all load is considered as " being possible to the load changed ", if the shared JZW load that may be changed, that is, have JZW secondary evaluation objects.

If the quantity of the generalized displacement of support component of Cable Structure, the quantity of the support cable of Cable Structure and JZW " being possible to the load changed " quantity sum are N, that is, have N number of evaluation object.Evaluation object serial number is given, the numbering will be used to generate vector sum matrix in subsequent step.

If evaluated support cable and the quantity sum of generalized displacement of support component are that P, the i.e. quantity of core evaluation object are P, if the quantity of evaluated generalized displacement of support component is Z, if the quantity of evaluated support cable is M₁。

" the monitored angle-data of the whole of structure " H angle coordinate components of L specified straight lines, each specified straight lines of K specified point, excessively each specified point in structure are described, the change for changing exactly all specified points, all specified straight lines all angle coordinate components specified of structural point.M (M=K × L × H) individual angle coordinate component measurement value or calculated value is had every time to characterize the angle information of structure.

The monitored amount of summary, whole Cable Structure has M monitored amounts, and the quantity that M cannot be less than core evaluation object plus 4, and M is less than the quantity N of evaluation object.

For convenience, " monitored all parameters of Cable Structure " are referred to as in the method " monitored amount ".To M monitored amount serial numbers, the numbering will be used to generate vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 ..., M with variable j.

Second step：Set up initial mechanical calculating benchmark model A_o。

When Cable Structure is completed, or before health monitoring systems are set up, the initial value for all monitored amounts for obtaining Cable Structure, the monitored amount initial value vector C of composition are calculated using conventional method direct measurement_o。

Obtaining monitored amount initial value vector C_oWhile, obtain the physical parameter of various materials used in Cable Structure and mechanical property parameters (such as modulus of elasticity, Poisson's ratio) using conventional method (consult reference materials or survey).

Obtaining monitored amount initial value vector C_oWhile, the Actual measurement data of Cable Structure are obtained using conventional method Actual measurement.The Actual measurement data of Cable Structure include the data that Non-destructive Testing Data of support cable etc. can express the health status of rope, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure bearing generalized coordinates data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structural point measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing generalized coordinates data refer to the bearing generalized coordinates data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model A_oWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs.The initial geometric data of Cable Structure can be spatial data of the spatial data plus a series of point in structure of the end points of all ropes, it is therefore intended that the geometric properties of Cable Structure are determined according to these coordinate datas.For cable-stayed bridge, if initial geometric data, which can be the spatial data of the end points of all ropes, adds the spatial data done on bridge two ends, here it is so-called bridge type data.Non-destructive Testing Data using support cable etc. can express the initial generalized displacement measurement data of the data of the health status of support cable, Cable Structure bearing and Cable Structure load measurement data set up evaluation object initial damage vector d_o, use d_oRepresent Cable Structure (with initial mechanical calculating benchmark model A_oRepresent) evaluation object initial health.If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0；If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial d_oIn each element numerical value related to Cable Structure generalized displacement of support take 0；If d_oThe corresponding evaluation object of some element be to take d in some load, this method_oThe element numerical value be 0, the initial value for representing the change of this load is 0.Using the physical and mechanical properties parameter of various materials used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure, initial mechanical calculating benchmark model A is set up using mechanics method (such as FInite Element)_o。

No matter which kind of method to obtain initial mechanical calculating benchmark model A with_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, error typically cannot be greater than 5%.So can utility A_oSuo Li under analog case obtained by calculating calculates that data, strain calculation data, Cable Structure shapometer count evidence and displacement meter counts evidence, Cable Structure angle-data, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A_oThe health status of middle support cable evaluation object initial damage vector d_oRepresent.Due to based on A_oThe initial value (actual measurement is obtained) for obtaining the evaluations of all monitored amounts closely all monitored amounts is calculated, so A can also be used in_oOn the basis of, carry out Mechanics Calculation obtain, A_oEach monitored amount the monitored amount initial value vector C of evaluation composition_o.Corresponding to A_oEvaluation object health status evaluation object initial damage vector d_oRepresent；Corresponding to A_oAll monitored amounts initial value with monitored amount initial value vector C_oRepresent.d_oIt is A_oParameter, C_oBy A_oMechanics Calculation result composition.

3rd step：In the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o；The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；dⁱ _oIt is Aⁱ _oCharacterisitic parameter；Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o。

4th step：The hardware components of cable structure health monitoring system are installed.Hardware components at least include：Monitored amount monitoring system (such as containing angle measurement system, signal conditioner), signal (data) collector, computer and communication alert equipment.Each monitored amount must be monitored system monitoring and arrive, and the signal monitored is transferred to signal (data) collector by monitoring system；Signal is delivered to computer through signal picker；Computer is then responsible for the health monitoring software of the evaluation object of operation Cable Structure, including the signal that the transmission of tracer signal collector comes；When monitoring that evaluation object health status is changed, computer control communication warning device is alarmed to monitoring personnel, owner and (or) the personnel specified.

5th step：Establishment and the on computers system software of installation and operation this method, the software will complete the functions (all work that can be completed with computer i.e. in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, alarm that this method required by task is wanted.

6th step：Thus step starts the cycle over running, during structure military service, and actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time.

7th step：In current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out according to step a to step d, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u。

A. when ith circulates beginning, method obtains Δ C as listed by step b to step dⁱAnd Dⁱ _u；At other moment, directly it is transferred to the 8th step and carries out follow-up work.

B. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of carry out Mechanics Calculation several times, vectorial dⁱ _oRepresent Aⁱ _oEvaluation object health status, calculation times are numerically equal to the quantity N of all evaluation objects, have N number of assessment object just to have n times calculating；Calculate each time and assume only one of which evaluation object in vectorial dⁱ _oUnit damage or unit generalized displacement or load unit change occur on the basis of the health status of the evaluation object of expression, if specifically, the evaluation object is a support cable in cable system, then it is assumed that the support cable is in vectorial dⁱ _oThe support cable represented has unit damage (for example taking 5%, 10%, 20% or 30% equivalent damage to be damaged for unit) again on the basis of having damaged, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing in the direction of displacement in vectorial dⁱ _oUnit generalized displacement occurs again (if such as evaluation object is the translational component in the x directions of a bearing on the basis of the existing generalized displacement of the bearing represented, it is assumed that the bearing has unit displacement of the lines in x directions, for example take 1mm, if the evaluation object is the angular displacement component around x-axis of a bearing, it is assumed that the bearing has unit angular displacement around x-axis, for example take ten a ten thousandth radians), if the evaluation object is a load, it is assumed that the load is in vectorial dⁱ _oLoad unit is further added by the basis of the existing variable quantity of the load represented to change (if the load is distributed load, and the distributed load is line distributed load, load unit change can take 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. to change for unit；If the load is distributed load, and the distributed load is EDS maps load, and load unit change can take 1MPa, 2MPa, 3MPa or 1kNm/m²、2kNm/m²、3kNm/m²Change Deng for unit；If the load is concentrfated load, and the concentrfated load is couple, and load unit change can take 1kNm, 2kNm, 3kNm etc. to change for unit；If the load is concentrfated load, and the concentrfated load is concentrated force, and load unit change can take 1kN, 2kN, 3kN etc. to change for unit；If the load is volume load, load unit change can take 1kN/m³、2kN/m³、 3kN/m³Change Deng for unit), use Dⁱ _ukThis unit damage or unit generalized displacement or load unit change are recorded, wherein k represents the numbering for occurring the evaluation object of unit damage or unit generalized displacement or load unit change, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uAn element, evaluation object unit change vector Dⁱ _uElement coding rule and vector d_oElement coding rule it is identical；The evaluation object for occurring unit damage or unit generalized displacement or load unit change in calculating each time is different from the evaluation object for occurring unit damage or unit generalized displacement or load unit change in other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector；When assuming that k-th of evaluation object has unit damage or unit generalized displacement or load unit change, C can be usedⁱ _tkRepresent corresponding " monitored amount calculation current vector "；When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed；Cⁱ _tkDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _tkElement and C_oElement correspond.

C. obtained vectorial C is calculated each timeⁱ _tkSubtract vectorial Cⁱ _oA vector is obtained, then the unit damage assumed in each element of the vector divided by this calculating or unit generalized displacement or load unit are changed into numerical value Dⁱ _uk" numerical value change vector δ C for monitored amount is obtained afterwardsⁱ _k”；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount ".

D. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "ⁱ”；Unit damage monitored numerical quantity transformation matrices Δ CⁱEach row correspond to a monitored amount unit change vector；Unit damage monitored numerical quantity transformation matrices Δ CⁱEvery a line correspond to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit generalized displacement or load unit change；Unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and vector d_oElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and M monitored amounts coding rules it is identical.

8th step：Set up linear relationship error vector eⁱWith vectorial gⁱ.Utilize data (" the monitored current initial value vector C of amount aboveⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ"); while the 8th step is calculated each time; i.e. while the increase unit damage for assuming only one of which evaluation object in evaluation object or unit generalized displacement or load unit change is calculated each time; when assuming that kth (k=1,2,3; ...; N) individual evaluation object increase unit damage or when unit generalized displacement or load unit change, one injury vector of composition is calculated each time, d is usedⁱ _tkThe injury vector is represented, corresponding monitored amount calculation current vector is Cⁱ _tk(referring to the 8th step), injury vector dⁱ _tkElement number be equal to evaluation object quantity, vectorial dⁱ _tkAll elements in only one of which element numerical value take each time calculate in assume increase unit damage unit generalized displacement or load unit change evaluation object unit damage or unit generalized displacement or load unit changing value, dⁱ _tkThe numerical value of other elements take 0, that for 0 element evaluation object of the numbering with assuming increase unit damage or unit generalized displacement or load unit change corresponding relation, with the elements of the same numberings of other vectors with the corresponding relation of the evaluation object be identical；dⁱ _tkWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _tkElement and d_oElement be one-to-one relationship.By Cⁱ _tk、Cⁱ _o、ΔCⁱ、dⁱ _tkBring formula (23) into, obtain a linear relationship error vector eⁱ _k, calculate obtain a linear relationship error vector e each timeⁱ _k；eⁱ _kSubscript k represent the individual evaluation object increase unit damage of kth (k=1,2,3 ..., N) or unit generalized displacement or load unit change.There is N number of evaluation object just to have n times calculating, just there is N number of linear relationship error vector eⁱ _k, by this N number of linear relationship error vector eⁱ _kA vector is obtained after addition, is exactly final linear relationship error vector e by the new vector obtained after each element divided by N of this vectorⁱ.Vectorial gⁱEqual to final error vector eⁱ.By vectorial gⁱOn the hard disc of computer for being stored in operation health monitoring systems software, used for health monitoring systems software.

9th step：Define current nominal fatigue vector dⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cAnd dⁱElement numerical value represent correspondence evaluation object degree of injury or unit generalized displacement or load change degree, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship.

Tenth step：According to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as formula (11), current nominal fatigue vector d calculated according to multi-objective optimization algorithmⁱ _cNoninferior solution, that is, with reasonable error but can relatively accurately reflect evaluation object health status change solution.

To solve current nominal fatigue vector d using the Objective Programming (GoalAttainmentMethod) of multi-objective optimization algorithmⁱ _c, according to Objective Programming, it is a real number that formula (11), which can convert γ in an accepted way of doing sth (24) and the multi-objective optimization question shown in formula (25), formula (24), and R is real number field, and area of space Ω limits vectorial dⁱ _cThe span of each element (the present embodiment requires vector dⁱ _cEach element be not less than 0, no more than 1).Formula (24) means one minimum real number γ of searching so that formula (25) is met.G (d in formula (25)ⁱ _c) defined by formula (25), G (d in weighing vector W and γ product representation formula (25) in formula (25)ⁱ _c) and vector gⁱBetween the deviation that allows, gⁱDefinition referring to formula (17), its value the 9th step calculate obtain.Vector W can be with vectorial g during actual calculatingⁱIt is identical.The specific programming realization of Objective Programming has had general program directly to use.Can be in the hope of current nominal fatigue vector d using Objective Programmingⁱ _c。

minimizeγ

11st step：According to the currently practical injury vector d of cable systemⁱDefinition (see formula (18)) and its element definition (see formula (19)) calculate obtain currently practical injury vector dⁱEach element, so as to by dⁱDetermine the health status of evaluation object.Currently practical injury vector dⁱK-th of element dⁱ _kRepresent the currently practical health status of k-th of evaluation object in ith circulation.

dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kRepresent its currently practical damage, dⁱ _kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity.

dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a generalized displacement component of a bearing, then dⁱ _kRepresent its currently practical generalized displacement numerical value.

dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a load, then dⁱ _kIts currently practical load change numerical value is represented, so according to the currently practical injury vector d of evaluation objectⁱIt can determine which support cable is damaged and its degree of injury, it is determined that there is which load to be changed and its numerical value.

So far this method realizes accurately identifying for the health status of core evaluation object in a kind of effective, cheap method.It may deviate the recognition result of the health status of secondary evaluation object that exact value is more, the health status of correct identification core evaluation object is required nothing more than in the method.

12nd step：Computer in health monitoring systems is periodically automatic or generates cable system health condition form by human users' health monitoring systems.

13rd step：Under specified requirements, the computer in health monitoring systems is automatically brought into operation communication alert equipment and alarmed to monitoring personnel, owner and (or) the personnel specified.

14th step：Mark vector B is set up according to formula (20)ⁱ, formula (21) gives mark vector BⁱK-th of element definition；If mark vector BⁱElement be all 0, then return to the 6th step and proceed health monitoring and calculating to cable system；If mark vector BⁱElement be not all 0, then complete after subsequent step, into circulating next time.

15th step：The initial damage vector d obtained next time needed for (i.e. i+1 time, i=1,2,3,4 ...) circulation is calculated according to formula (22)ⁱ⁺¹ _oEach element dⁱ⁺¹ _ok(k=1,2,3 ..., N)；Second, in initial mechanical calculating benchmark model A_oOn the basis of, the health status for making rope is dⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time (i=1,2,3,4 ...) circulationⁱ⁺¹.Obtain Aⁱ⁺¹、dⁱ⁺¹ _oAfterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o。

16th step：The 6th step is returned to, starts the circulation by the 6th step to the 16th step.

Claims (1)

1. simplify angle monitor load damaged cable generalized displacement progressive recognition method, it is characterised in that methods described includes：

A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable；The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation；This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, i.e. evaluation object is N；The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step；This method represents this numbering, k=1,2,3 ..., N with variable k；This method refers exclusively to evaluated support cable and generalized displacement of support component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of generalized displacement of support component are P, i.e. the quantity of core evaluation object is P, and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；If having M in cable system₁Root support cable；It is determined that the measured point specified, gives all specified points numbering；It is determined that the measured straight line of excessively each measurement point, to all measured straight line numberings specified；The measured angle coordinate component of each measured straight line is determined, is numbered to all measured angle coordinate components；Above-mentioned numbering will be used to generate vector sum matrix in subsequent step；" the monitored angle-data of the whole of Cable Structure " is made up of above-mentioned all measured angle coordinate components；For convenience, " the monitored angle-data of Cable Structure " is referred to as in the method " monitored amount "；The quantity sum of all monitored amounts is designated as M, and M should be more than the quantity of core evaluation object, and M is less than the quantity of evaluation object；The external force that object, structure are born can be described as load, and load includes face load and volume load；Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load；Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force；Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component；Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant；Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load；Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes；Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method；

B. survey or consult reference materials and obtain the physical and mechanical properties parameter of various materials used in Cable Structure；

C. while surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, direct measurement calculates the measured data for obtaining initial Cable Structure, and the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, and the data of the health status that can express support cable now are referred to as support cable initial health data；The monitored amount initial value vector C of initial value composition of all monitored amounts_o, it is monitored amount initial value vector C_oCoding rule and M monitored amounts coding rules it is identical；Evaluation object initial damage vector d is set up using support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement data_o, vectorial d_oRepresent with initial mechanical calculating benchmark model A_oThe initial health of the evaluation object of the Cable Structure of expression；Evaluation object initial damage vector d_oElement number be equal to N, d_oElement and evaluation object be one-to-one relationship, vectorial d_oElement coding rule it is identical with the coding rule of evaluation object；If d_oThe corresponding evaluation object of some element be a support cable in cable system, then d_oThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion；If d_oThe corresponding evaluation object of some element be some bearing some generalized displacement component, then d_oThe element numerical value represent this bearing the generalized displacement component initial value；If d_oThe corresponding evaluation object of some element be to take d in some load, this method_oThe element numerical value be 0, the initial value for representing the change of this load is 0；If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial d_oIn each element numerical value related to Cable Structure generalized displacement of support take 0；If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0；Initial Cable Structure bearing generalized coordinates data refer to the bearing generalized coordinates data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model A_oWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs；

D. all Cable Structure data that the physical and mechanical properties parameter and preceding step of various materials according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure are obtained, set up the initial mechanical calculating benchmark model A of Cable Structure_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot be greater than 5%；Corresponding to A_oEvaluation object health status evaluation object initial damage vector d_oRepresent；Corresponding to A_oAll monitored amounts initial value with monitored amount initial value vector C_oRepresent；d_oIt is A_oParameter, by A_oThe obtained initial value of all monitored amounts of Mechanics Calculation result and C_oThe initial value of all monitored amounts represented is identical, therefore alternatively C_oBy A_oMechanics Calculation result composition, A in the method_o、C_oAnd d_oIt is constant；

E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o；The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；dⁱ _oIt is Aⁱ _oCharacterisitic parameter, Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o；

F. enter from here and the circulation walked to pth is walked by f；

G. during structure military service, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time；

H. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u；

H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4ⁱAnd Dⁱ _u；At other moment, when in step g to Aⁱ _oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CⁱAnd Dⁱ _uIf, not to A in step gⁱ _oIt is updated, then is directly transferred to step i here and carries out follow-up work；

H2. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating；According to the coding rule of evaluation object, calculated successively；Calculate each time and assume that only one of which evaluation object is further added by unit damage or unit generalized displacement or load unit change on the basis of original damage or generalized displacement or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing is further added by unit generalized displacement in the direction of displacement, if the evaluation object is a load, it is assumed that the load is further added by load unit change, use Dⁱ _ukThis increased unit damage or unit generalized displacement or load unit change are recorded, wherein k represents to increase the numbering of the evaluation object of unit damage or unit generalized displacement or load unit change, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uAn element, evaluation object unit change vector Dⁱ _uElement coding rule and vector d_oElement coding rule it is identical；The evaluation object that unit damage or unit generalized displacement or load unit change are further added by calculating each time is different from the evaluation object that unit damage or unit generalized displacement or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector；When assuming that k-th of evaluation object is further added by unit damage or unit generalized displacement or load unit change, C is usedⁱ _tkRepresent corresponding " monitored amount calculation current vector "；When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed；Cⁱ _tkDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _tkElement and C_oElement correspond；

H3. obtained vectorial C is calculated each timeⁱ _tkSubtract vectorial Cⁱ _oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or unit generalized displacement or load unit change numerical valueⁱ _k”；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount "；

H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "ⁱ”；Unit damage monitored numerical quantity transformation matrices Δ CⁱEach row correspond to a monitored amount unit change vector；Unit damage monitored numerical quantity transformation matrices Δ CⁱEvery a line correspond to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit generalized displacement or load unit change；Unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and vector d_oElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and M monitored amounts coding rules it is identical；

I. current nominal fatigue vector d is definedⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cElement numerical value represent the nominal fatigue degree or nominal generalized displacement or nominal load variable quantity of correspondence evaluation object, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship；

J. according to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1ⁱ _cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector dⁱ _c；

$C^i=C_o^i+\mathrm{\Δ}C\·d_c^i$ Formula 1

K. the currently practical injury vector d expressed using formula 2ⁱK-th of element dⁱ _kWith the current initial damage vector d of evaluation objectⁱ _oK-th of element dⁱ _okWith current nominal fatigue vector dⁱ _cK-th of element dⁱ _ckBetween relation, calculating obtain currently practical injury vector dⁱAll elements；

Formula 2

K=1,2,3 in formula 2 ..., N；dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kRepresent its currently practical damage, dⁱ _kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity；If the evaluation object is a load, then dⁱ _kRepresent the actual change amount of the load；If the evaluation object is a generalized displacement component of a bearing, then dⁱ _kRepresent its currently practical generalized displacement numerical value；So far this method, which is realized, rejects generalized displacement of support, load change influence, Cable Structure damaged cable identification, realize simultaneously and reject generalized displacement of support and the change influence of support cable health status, load change amount identification, also achieve and reject load change and the change influence of support cable health status, generalized displacement of support identification；So far this method realizes accurately identifying for the health status of core evaluation object in a kind of effective, cheap method；Recognition result deviation exact value to the health status of secondary evaluation object is more, therefore not accepts and believe, and the health status of correct identification core evaluation object is required nothing more than in the method；

L. current nominal fatigue vector d is tried to achieveⁱ _cAfterwards, mark vector B is set up according to formula 3ⁱ, formula 4 gives mark vector BⁱK-th of element definition；

$B^i={\left[\begin{array}{cccccccccc}{B}_1^i& B_2^i& \·& \·& \·& B_k^i& \·& \·& \·& B_N^i\end{array}\right]}^T$ Formula 3

Formula 4

Element B in formula 4ⁱ _kIt is mark vector BⁱK-th of element, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element, dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4；

If m. mark vector BⁱElement be all 0, then return to step f continue this circulation；If mark vector BⁱElement be not all 0, then enter next step, i.e. step n；

N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5ⁱ⁺¹ _oEach element；

Formula 5

D in formula 5ⁱ⁺¹ _okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationⁱ⁺¹ _oK-th of element, dⁱ _okThis, i.e. the current initial damage vector d of evaluation object of ith circulationⁱ _oK-th of element, Dⁱ _ukIt is the evaluation object unit change vector D of ith circulationⁱ _uK-th of element, Bⁱ _kIt is the mark vector B of ith circulationⁱK-th of element, k=1,2,3 ... ..., N in formula 5；

O. in initial mechanical calculating benchmark model A_oOn the basis of, the health status for making rope is dⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationⁱ⁺¹；Obtain Aⁱ⁺¹Afterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o；

P. step f is returned to, starts to circulate next time.