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

Abstract

The invention provides a method for progressively recognizing the damaged cable of a load based on the streamlined hybrid monitoring process of the generalized displacement. According to the method, based on the hybrid monitoring process, whether a mechanical calculation benchmark model of a cable structure needs to be updated or not is determined through monitoring the generalized displacement of a bearer, the variation degree of a load and the damaged degree of a damaged cable. After that, the mechanical calculation benchmark model of the cable structure, with the generalized displacement of the bearer, the variation degree of the load and the damaged degree of the damaged cable being taken into account, is obtained. On the basis of the above model, according to the approximately linear relationships between the current numeric vector of a 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 current nominal damage vector, the health status of a core to-be-evaluated object can be recognized.

Description

Simplify generalized displacement hybrid monitoring loading problem rope progressive-type 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.The impaired and lax pair Cable Structure of support cable is a significant threat safely, and damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method.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.During Cable Structure military service, generalized displacement may occur for Cable Structure bearing, the load that Cable Structure is born is it can also happen that change, the health status of Cable Structure may also change simultaneously, in this complex condition, based on hybrid monitoring, (this method judges the health status of Cable Structure to this method by the hybrid monitoring of the change of the measurable parameter to the foregoing different types of Cable Structure of this section, all monitored Cable Structure characteristic parameters are referred to as " monitored amount " by this method, it is made up of due to being now monitored amount the different types of measurable parameter mixing of Cable Structure, this method is referred to as hybrid monitoring) recognize problem cable, belong to engineering structure health monitoring field.

Background technology

Reject the influence of load change, Cable Structure generalized displacement of support to Cable Structure health status recognition result, so that the change of the health status of identification structure exactly, the problem of being in the urgent need to address at present, this method discloses a kind of effective, the cheap method for solving this problem.

The content of the invention

Technical problem：This method discloses a kind of method, under conditions of cost is lower, when bearing has generalized displacement, in the load change that structure is born, the influence of generalized displacement of support and load change to Cable Structure health status recognition result can be rejected, so as to identify the health status of support cable exactly.

During Cable Structure military service, Suo Changdu under support cable free state (it is 0 that now rope tensility, which is also referred to as Suo Li) (is referred to as drift, this method refers exclusively to the drift of that section of rope between the supporting end points of support cable two) it can change, the first purpose of this method seeks to identify the support cable that drift is changed, and identify the knots modification of their drift, this knots modification provides direct basis for the cable force adjustment of the rope, for convenience, the support cable that drift changes is referred to as slack line by this method.

Technical scheme：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 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.

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 support cable and load is " evaluation object ", if the evaluated quantity of support cable and the quantity sum of load is that N, the i.e. quantity of " evaluation object " are N；This method refers exclusively to the evaluated support cable in " evaluation object " with title " core evaluation object ", and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；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；The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, that is, a Q core evaluation object is had, the monitored rope force data of the Cable Structure M in Cable Structure₁The M of individual specified support cable₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change；M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure；M₁It is an integer for being not more than Q not less than 0；Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point；M is had every time₂Individual strain measurement value or calculated value strain to characterize Cable Structure, M₂For K₂And L₂Product；M₂It is no less than 0 integer；Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified；M is had every time₃Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M₃For K₃、L₃And H₃Product；M₃It is an integer not less than 0；The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point；M is had every time₄Individual coordinates measurements or calculated value characterize Cable Structure shape, M₄For K₄And L₄Product；M₄It is an integer not less than 0；The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M₁、M₂、M₃And M₄Sum, it is M to define parameter K, K₁、K₂、K₃And K₄Sum, M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object；For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method；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；While surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, direct measurement calculates the Initial cable force for obtaining all support cables, composition Initial cable force vector F_o；According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cable_oElement coding rule it is identical；

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, 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 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, 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 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 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 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 constitute initial Cable Structure bearing generalized coordinates vector U_o；

D. the physical and mechanical properties parameter of various materials, initial Cable Structure bearing generalized coordinates vector U 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, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure_oAll Cable Structure data obtained with preceding step, 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_oCable Structure bearing generalized coordinates data be exactly initial Cable Structure bearing generalized coordinates vector U_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；U_oAnd 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_o、d_oAnd U_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；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure bearing generalized coordinates data " with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent, vectorial Uⁱ _oDefinition mode and vector U_oDefinition mode it is identical, Uⁱ _oElement and U_oElement correspond；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；Uⁱ _oAnd 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, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_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 q is walked by f；During structure military service, constantly actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas constitute current Cable Structure actual measurement bearing angular coordinate vector Uⁱ, vectorial UⁱDefinition mode and vector U_oDefinition mode it is identical, UⁱElement and U_oElement correspond；Vectorial U is obtained in actual measurementⁱWhile, 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；Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, actual measurement obtains the rope force datas of all Q roots support cables in Cable Structure, all these rope force datas composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, Actual measurement obtains the space coordinate of two supporting end points of all Q roots support cables, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vector_oElement coding rule it is identical；

G. bearing generalized coordinates vector U is surveyed according to current Cable Structureⁱ, current initial mechanical calculating benchmark model A is updated according to step g1 to g3ⁱ _o, the monitored current initial value vector C of amountⁱ _oWith current initial Cable Structure bearing generalized coordinates vector Uⁱ _o, and the current initial damage vector d of evaluation objectⁱ _oKeep constant；

G1. U is comparedⁱWith Uⁱ _oIf, UⁱEqual to Uⁱ _o, then Aⁱ _o、Cⁱ _oAnd Uⁱ _oKeep constant, otherwise need to follow these steps to Aⁱ _o、Cⁱ _oAnd Uⁱ _oIt is updated；

G2. U is calculatedⁱWith U_oDifference, UⁱWith U_oDifference be exactly generalized displacement of support of the Cable Structure bearing on initial position, represent generalized displacement of support with generalized displacement of support vector V, V is equal to UⁱSubtract U_o；

G3. to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, to A_oIn Cable Structure bearing apply the current initial mechanical calculating benchmark model A that is updated after generalized displacement of support constraintⁱ _o, update Aⁱ _oWhile, Uⁱ _oAll elements numerical value also uses UⁱAll elements numerical value correspondence is replaced, that is, have updated Uⁱ _o, thus obtained properly corresponding to Aⁱ _oUⁱ _o, now dⁱ _oKeep constant；As renewal Aⁱ _oAfterwards, Aⁱ _oRope the health status current initial damage vector d of evaluation objectⁱ _oRepresent, Aⁱ _oBearing generalized coordinates with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent；Update Cⁱ _oMethod be：As renewal Aⁱ _oAfterwards, A is obtained by Mechanics Calculationⁱ _oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Cⁱ _o；

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 load unit change on the basis of original damage 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 a load, it is assumed that the load is further added by load unit change, D is usedⁱ _ukThis increased unit damage or load unit change is recorded, wherein k represents the numbering for increasing unit damage or the evaluation object of 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 load unit change are further added by calculating each time is different from the evaluation object that unit damage 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 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 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 same monitored amount different evaluation objects increase unit damage or load unit change when different unit change amplitudes；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 correspondence evaluation object nominal fatigue degree or nominal load variable quantity, 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；

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ⁱ _kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, dⁱ _kThe numerical response degree of relaxation or the damage of the support cable；Damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line, the currently practical injury vector d of evaluation objectⁱIn correspond to slack line element numerical expression be with slack line relax degree mechanic equivalent currently practical equivalent damage degree；If the evaluation object is a load, then dⁱ _kRepresent the actual change amount of the load；By the currently practical injury vector d of evaluation objectⁱIn the Q element related to Q root support cables take out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical；The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., Q；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line；Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity；Identify after slack line, utilize the currently practical injury vector d of support cable^ciThese slack lines, the degree that relaxed with it mechanic equivalent the currently practical equivalent damage degree of expression, utilizes current cable force vector FⁱEnd points horizontal range vector is supported with current support cable two, using the initial drift that the support cable obtained is walked in b is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial, Initial cable force vector F_oUtilize the physical and mechanical properties parameter that various materials used in the Cable Structure of acquisition are walked in b, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is：First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical；2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope；When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa；Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined；So it is achieved that relaxation identification and the non-destructive tests of support cable；Institute's demand power is by current cable force vector F during calculatingⁱCorresponding element is provided；Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far the problem of this method realizes rejecting load change and the influence of generalized displacement of support, Cable Structure rope is recognized, generalized displacement of support and the change influence of support cable health status, load change amount identification are rejected while realizing；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

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；

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, first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then makes the health status of rope for 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. remove once, i.e. i+1 time circulates required current initial Cable Structure bearing generalized coordinates vector Uⁱ⁺¹ _oThe current initial Cable Structure bearing generalized coordinates vector U circulated equal to ithⁱ _o；

Q. 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 sum of the quantity and JZW " being possible to the load changed " of the support cable of Cable Structure is 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.

Monitored multiclass parameter can include：Suo Li, strain, angle and space coordinate, are described below respectively：

If having Q root support cables in cable system, that is, have a Q core evaluation object, the monitored rope force data of the Cable Structure M in Cable Structure₁The M of individual specified rope₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified ropes change.M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M₁It is an integer not less than 0.

The monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point.M is had every time₂(M₂=K₂×L₂) individual strain measurement value or calculated value characterize Cable Structure strain.M₂It is an integer not less than 0.

The monitored angle-data of Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified.M is had every time₃(M₃=K₃×L₃×H₃) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M₃It is an integer not less than 0.

The monitored shape data of Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point.M is had every time₄(M₄=K₄×L₄) individual coordinates measurements or calculated value characterize Cable Structure shape.M₄It is an integer not less than 0.

The monitored amount of summary, whole Cable Structure has M (M=M₁+M₂+M₃+M₄) individual monitored amount, define parameter K (K=M₁+K₂+K₃+K₄), the quantity that M cannot be less than core evaluation object plus 4, 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.

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；The physical parameter (such as density) and mechanical property parameters (such as modulus of elasticity, Poisson's ratio) of various materials used in Cable Structure are obtained using conventional method (consult reference materials or survey) simultaneously.

While surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, support cable initial elastic modulus vector E is constituted using the modulus of elasticity data of support cable therein_o；While surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, direct measurement calculates the Initial cable force for obtaining all support cables, composition Initial cable force vector F_o；According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li be 0 when length, in free state when cross-sectional area and the unit length in free state weight, successively composition support cable initial drift vector l_o, initial free cross-sectional area vector A_oWith the weight vector ω of initial free unit length_o, support cable initial elastic modulus vector E_o, support cable initial drift vector l_o, initial free cross-sectional area vector A_oWith the weight vector ω of initial free unit length_oElement coding rule and Initial cable force vector F_oElement coding rule it is identical.

Monitored amount initial value vector C is obtained in Actual measurement_oWhile, the Actual measurement data of Cable Structure are obtained using conventional method Actual measurement.The Actual measurement data of Cable Structure, which include Non-destructive Testing Data of support cable etc., can express the measured datas such as the initial geometric data of data, Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinates data, 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 structure space measurement of coordinates data of health status of rope.Initial Cable Structure bearing generalized coordinates data constitute initial Cable Structure bearing generalized coordinates vector U_o.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.The data and Cable Structure load measurement data that Non-destructive Testing Data using support cable etc. can express the health status of support cable 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 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.Utilize 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 and initial Cable Structure bearing generalized coordinates vector U_o, 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, model A_oMiddle bearing generalized coordinates vector U_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.Corresponding to A_oCable Structure bearing generalized coordinates data with initial Cable Structure bearing generalized coordinates vector U_oRepresent；U_oAnd 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；Current initial mechanical calculating benchmark model A that ith circulation needs when starting, corresponding to Cable Structureⁱ _oThe current initial Cable Structure bearing generalized coordinates vector U of Cable Structure bearing generalized coordinates data compositionⁱ _o, the current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timeⁱ _oWhen, Uⁱ _oIt is equal to U_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；Uⁱ _oAnd 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, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_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 subsystem containing angular surveying, cable force measurement subsystem, strain measurement subsystem, space coordinate measurement subsystem, signal conditioner), Cable Structure bearing generalized coordinates monitoring system (containing total powerstation, angle measuring sensor, signal conditioner etc.), support cable cable force monitoring system, space coordinate monitoring system, signal (data) collector, computer and the communication alert equipment of the supporting end points of support cable.Each monitored amount, the Suo Li of each support cable, the space coordinate of the supporting end points of each support cable, the bearing generalized coordinates of each Cable Structure 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 continuous 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.

Monitored amount current value vector C is obtained in actual measurementⁱWhile, actual measurement obtains Cable Structure bearing generalized coordinates current data, all data composition current cable structure actual measurement bearing generalized coordinates vector Uⁱ。

Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, actual measurement obtains the rope force datas of all Q roots support cables in Cable Structure, all these rope force datas composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, Actual measurement obtains the space coordinate of two supporting end points of all Q roots support cables, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, and two supporting end points horizontal range data of all Q roots support cables constitute the current supporting end points horizontal range vector of support cables two lⁱ _x, the current supporting end points horizontal range vector of support cable two lⁱ _xElement coding rule and Initial cable force vector F_oElement coding rule it is identical.

7th step：Obtaining current Cable Structure actual measurement bearing generalized coordinates vector UⁱAfterwards, U is comparedⁱAnd Uⁱ _oIf, UⁱEqual to Uⁱ _o, then need not be to Aⁱ _o、Cⁱ _oAnd Uⁱ _oIt is updated, otherwise needs to current initial mechanical calculating benchmark model Aⁱ _o, current initial Cable Structure bearing generalized coordinates vector Uⁱ _oWith monitored amount current initial value vector Cⁱ _oIt is updated, and the current initial damage vector d of evaluation objectⁱ _oKeep constant, update method follows these steps a to step c progress：

A. U is calculatedⁱWith U_oDifference, UⁱWith U_oDifference be exactly generalized displacement of support of the Cable Structure bearing on initial position, represent generalized displacement of support with generalized displacement of support vector V, V is equal to UⁱSubtract U_o, it is the generalized displacement that the numerical value of an element in one-to-one relationship, generalized displacement of support vector V corresponds to an assigned direction of a specified bearing between the element and generalized displacement of support component in generalized displacement of support vector V.

B. to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, to A_oIn Cable Structure bearing apply the current initial mechanical calculating benchmark model A that is updated after generalized displacement of support constraintⁱ _o, update Aⁱ _oWhile, Uⁱ _oAll elements numerical value also uses UⁱAll elements numerical value correspondence is replaced, that is, have updated Uⁱ _o, thus obtained properly corresponding to Aⁱ _oUⁱ _o, now dⁱ _oKeep constant；As renewal Aⁱ _oAfterwards, Aⁱ _oRope the health status current initial damage vector d of evaluation objectⁱ _oRepresent, Aⁱ _oBearing generalized coordinates with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent.

C. C is updatedⁱ _oMethod be：As renewal Aⁱ _oAfterwards, A is obtained by Mechanics Calculationⁱ _oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Cⁱ _o；

8th 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 directly as listed by step b to step dⁱAnd Dⁱ _u；At other moment, when in the 7th step to Aⁱ _oAfter being updated, it is necessary to which the method as listed by step b to step d regains Δ CⁱAnd Dⁱ _uIf, not to A in the 7th stepⁱ _oIt is updated, then is directly transferred to the 9th step here 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 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 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 load unit change are recorded, wherein k represents the numbering for occurring unit damage or the evaluation object of 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 load unit change in calculating each time is different from the evaluation object for occurring unit damage 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 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 or load unit that assume in each element of the vector divided by this calculating 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 same monitored amount different evaluation objects increase unit damage or load unit change when different unit change amplitudes；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.

9th 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 load unit change is calculated each time; when assuming that kth (k=1,2,3; ...; N) individual evaluation object increase unit damage or during 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 or load unit change evaluation object unit damage 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 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 (1) 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 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.

Tenth 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 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.

11st 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 (2), current nominal fatigue vector d calculated according to multi-objective optimization algorithmⁱ _cNoninferior solution, that is, position and its solution of nominal fatigue degree of damaged cable can be determined with reasonable error but relatively accurately from all ropes.

The Objective Programming (GoalAttainmentMethod) in multi-objective optimization algorithm etc. can be used to solve current nominal fatigue vector dⁱ _c, according to Objective Programming, it is a real number that formula (2), which can convert γ in an accepted way of doing sth (3) and the multi-objective optimization question shown in formula (4), formula (3), and R is real number field, and area of space Ω limits vectorial dⁱ _cEach element span.Formula (3) means one minimum real number γ of searching so that formula (4) is met.G (d in formula (4)ⁱ _c) defined by formula (5), G (d in weighing vector W and γ product representation formula (4) in formula (4)ⁱ _c) and vector gⁱBetween the deviation that allows.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。

12nd step：According to the currently practical injury vector d of cable systemⁱDefinition (see formula (6)) and its element definition (see formula (7)) 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 in formula (6)ⁱ _k(i=1,2,3 ...；K=1,2,3 ... ..., N) represent ith circulation in k-th of evaluation object currently practical health status, formula (7) is shown in its definition, 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, vectorial dⁱElement coding rule and formula (1) in vector d_oElement coding rule it is identical.

D in formula (7)ⁱ _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ⁱ _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 its corresponding support cable without health problem, dⁱ _kNumerical value represents that its corresponding support cable is the support cable of unsoundness problem when not being 0, and the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, the degree of the relaxation of its numerical response or damage.

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ⁱ _kRepresent its currently practical load change numerical value.

By the currently practical injury vector d of evaluation objectⁱIn the Q element related to support cable take out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical.The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., Q；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line.The method for distinguishing that reflects is varied; can be by the protective layer of the support cable for removing unsoundness problem; support cable is visually differentiated; or visually differentiated by optical imaging apparatus; it can also differentiate by the way that whether lossless detection method is impaired to support cable, ultrasonic examination is exactly a kind of now widely used lossless detection method.Those do not find that damage, unsoundness problem support cable is exactly that there occurs loose rope after discriminating, exactly need to adjust Suo Li rope, are exactly slack line, these need to adjust Suo Li rope in the currently practical injury vector d of support cable^ciIn corresponding element numerical value (such as one of element can use d^ci _hRepresent) degree of injury with the relaxation degree mechanic equivalents of these support cables is represented, slack line is thus determined that, the computational methods of specific slack illustrate below.Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity；According to the currently practical injury vector d of support cable^ci, identified from the support cable of unsoundness problem after slack line, remaining is exactly damaged cable, these damaged cables are in the currently practical injury vector d of support cable^ciIn the numerical value of corresponding element mean that its degree of injury, the numerical value of corresponding element represents that the support cable thoroughly loses bearing capacity when being 100%, represent that the support cable loses the bearing capacity of corresponding proportion when between 0 and 100%, so far just have identified damaged cable and its degree of injury.

Utilize the currently practical injury vector d of support cable^ciThe currently practical equivalent damage degree with its degree mechanic equivalent that relaxes of the slack line of expression, monitored amount current value vector C is obtained using in the 6th step, in actual measurementⁱSynchronization obtain current cable force vector FⁱWith current support cable two supporting end points horizontal range vector lⁱ _x, (can also set up the initial mechanical calculating benchmark model A of Cable Structure using in second step, while surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure_oSynchronization) obtain support cable initial elastic modulus vector E_o, support cable initial drift vector l_o, initial free cross-sectional area vector A_oWith the weight vector ω of initial free unit length_oUtilize the physical and mechanical properties parameter of various materials used in the Cable Structure obtained in second step, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, specifically can according to formula (8) or formula (9) can in the hope of these ropes relaxation degree (i.e. the long adjustment amount of rope).So it is achieved that the relaxation identification of support cable.So far damaged cable and slack line all just are identified.

E in formula (8) and formula (9)_ohIt is to set up the initial mechanical calculating benchmark model A of Cable Structure_oSynchronization, the modulus of elasticity of h root support cables, A_ohIt is to set up the initial mechanical calculating benchmark model A of Cable Structure_oSynchronization, the cross-sectional area of h root support cables, ω_ohIt is to set up the initial mechanical calculating benchmark model A of Cable Structure_oSynchronization, the weight of the unit length of h root support cables, l_ohIt is to set up the initial mechanical calculating benchmark model A of Cable Structure_oSynchronization, the drift of h root support cables, E in Cable Structure_ohIt is vectorial E_oH-th of element, be to set up the initial mechanical calculating benchmark model A of Cable Structure_oSynchronization, the modulus of elasticity of h root support cables, F in Cable Structureⁱ _hIt is to obtain monitored amount current value vector C in actual measurementⁱSynchronization, the current cable power of h root support cables, d^ci _hIt is the currently practical degree of injury of h root support cables, lⁱ _xhIt is to obtain monitored amount current value vector C in actual measurementⁱSynchronization, the horizontal range of two of h root support cables supporting end points, lⁱ _xhIt is the current supporting end points horizontal range vector of support cable two lⁱ _xAn element, current support cable two supporting end points horizontal range vector lⁱ _xElement coding rule and initial drift vector l_oElement coding rule it is identical.Item in formula (9) in [] is the Ernst equivalent elastic modulus of the support cable, and the vectorial Δ l of support cable current slack degree can be assured that by formula (8) or formula (9)ⁱ.Formula (9) is the amendment to formula (8).

So according to the currently practical injury vector d of evaluation objectⁱProblem cable and its health degree are can determine, 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.

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

14th 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.

15th step：Set up mark vector Bⁱ, formula (11) 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.

Mark vector B in formula (10)ⁱ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 (11)ⁱ _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 the relevant information of k-th of evaluation object.

16th 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 (12)ⁱ⁺¹ _oEach element dⁱ⁺¹ _ok(k=1,2,3 ..., N)；Second, in initial mechanical calculating benchmark model A_oOn the basis of, to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then makes the health status of rope for 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ⁱ⁺¹；Current initial Cable Structure bearing generalized coordinates vector U needed for (i.e. i+1 time, i=1,2,3,4 ...) is circulated next timeⁱ⁺¹ _oEqual to Uⁱ _o.Obtain Aⁱ⁺¹、dⁱ⁺¹ _oAnd Uⁱ⁺¹ _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。

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

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

Claims (1)

1. simplify generalized displacement hybrid monitoring loading problem rope progressive-type 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 support cable and load is " evaluation object ", if the evaluated quantity of support cable and the quantity sum of load is that N, the i.e. quantity of " evaluation object " are N；This method refers exclusively to the evaluated support cable in " evaluation object " with title " core evaluation object ", and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；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；The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, that is, a Q core evaluation object is had, the monitored rope force data of the Cable Structure M in Cable Structure₁The M of individual specified support cable₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change；M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure；M₁It is an integer for being not more than Q not less than 0；Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point；M is had every time₂Individual strain measurement value or calculated value strain to characterize Cable Structure, M₂For K₂And L₂Product；M₂It is no less than 0 integer；Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified；M is had every time₃Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M₃For K₃、L₃And H₃Product；M₃It is an integer not less than 0；The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point；M is had every time₄Individual coordinates measurements or calculated value characterize Cable Structure shape, M₄For K₄And L₄Product；M₄It is an integer not less than 0；The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M₁、M₂、M₃And M₄Sum, it is M to define parameter K, K₁、K₂、K₃And K₄Sum, M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object；For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method；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；While surveying or consulting reference materials the physical and mechanical properties parameter for obtaining various materials used in Cable Structure, direct measurement calculates the Initial cable force for obtaining all support cables, composition Initial cable force vector F_o；According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cable_oElement coding rule it is identical；

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, 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 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, 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 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 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 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 constitute initial Cable Structure bearing generalized coordinates vector U_o；

D. the physical and mechanical properties parameter of various materials, initial Cable Structure bearing generalized coordinates vector U 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, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure_oAll Cable Structure data obtained with preceding step, 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_oCable Structure bearing generalized coordinates data be exactly initial Cable Structure bearing generalized coordinates vector U_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；U_oAnd 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_o、d_oAnd U_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；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure bearing generalized coordinates data " with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent, vectorial Uⁱ _oDefinition mode and vector U_oDefinition mode it is identical, Uⁱ _oElement and U_oElement correspond；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；Uⁱ _oAnd 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, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_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 q is walked by f；During structure military service, constantly actual measurement obtains Cable Structure bearing angular coordinate current data, and all Cable Structure bearing angular coordinate current datas constitute current Cable Structure actual measurement bearing angular coordinate vector Uⁱ, vectorial UⁱDefinition mode and vector U_oDefinition mode it is identical, UⁱElement and U_oElement correspond；Vectorial U is obtained in actual measurementⁱWhile, 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；Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, actual measurement obtains the rope force datas of all Q roots support cables in Cable Structure, all these rope force datas composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；Monitored amount current value vector C is obtained in actual measurementⁱSynchronization, Actual measurement obtains the space coordinate of two supporting end points of all Q roots support cables, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vector_oElement coding rule it is identical；

G. bearing generalized coordinates vector U is surveyed according to current Cable Structureⁱ, current initial mechanical calculating benchmark model A is updated according to step g1 to g3ⁱ _o, the monitored current initial value vector C of amountⁱ _oWith current initial Cable Structure bearing generalized coordinates vector Uⁱ _o, and the current initial damage vector d of evaluation objectⁱ _oKeep constant；

G1. U is comparedⁱWith Uⁱ _oIf, UⁱEqual to Uⁱ _o, then Aⁱ _o、Cⁱ _oAnd Uⁱ _oKeep constant, otherwise need to follow these steps to Aⁱ _o、Cⁱ _oAnd Uⁱ _oIt is updated；

G2. U is calculatedⁱWith U_oDifference, UⁱWith U_oDifference be exactly generalized displacement of support of the Cable Structure bearing on initial position, represent generalized displacement of support with generalized displacement of support vector V, V is equal to UⁱSubtract U_o；

G3. to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, to A_oIn Cable Structure bearing apply the current initial mechanical calculating benchmark model A that is updated after generalized displacement of support constraintⁱ _o, update Aⁱ _oWhile, Uⁱ _oAll elements numerical value also uses UⁱAll elements numerical value correspondence is replaced, that is, have updated Uⁱ _o, thus obtained properly corresponding to Aⁱ _oUⁱ _o, now dⁱ _oKeep constant；As renewal Aⁱ _oAfterwards, Aⁱ _oRope the health status current initial damage vector d of evaluation objectⁱ _oRepresent, Aⁱ _oBearing generalized coordinates with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent；Update Cⁱ _oMethod be：As renewal Aⁱ _oAfterwards, A is obtained by Mechanics Calculationⁱ _oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Cⁱ _o；

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 load unit change on the basis of original damage 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 a load, it is assumed that the load is further added by load unit change, D is usedⁱ _ukThis increased unit damage or load unit change is recorded, wherein k represents the numbering for increasing unit damage or the evaluation object of 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 load unit change are further added by calculating each time is different from the evaluation object that unit damage 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 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 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 same monitored amount different evaluation objects increase unit damage or load unit change when different unit change amplitudes；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 correspondence evaluation object nominal fatigue degree or nominal load variable quantity, 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ⁱ _kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, dⁱ _kThe numerical response degree of relaxation or the damage of the support cable；Damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line, the currently practical injury vector d of evaluation objectⁱIn correspond to slack line element numerical expression be with slack line relax degree mechanic equivalent currently practical equivalent damage degree；If the evaluation object is a load, then dⁱ _kRepresent the actual change amount of the load；By the currently practical injury vector d of evaluation objectⁱIn the Q element related to Q root support cables take out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical；The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., Q；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line；Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity；Identify after slack line, utilize the currently practical injury vector d of support cable^ciThese slack lines, the degree that relaxed with it mechanic equivalent the currently practical equivalent damage degree of expression, utilizes current cable force vector FⁱEnd points horizontal range vector is supported with current support cable two, using the initial drift that the support cable obtained is walked in b is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial, Initial cable force vector F_oUtilize the physical and mechanical properties parameter that various materials used in the Cable Structure of acquisition are walked in b, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is：First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical；2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope；When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa；Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined；So it is achieved that relaxation identification and the non-destructive tests of support cable；Institute's demand power is by current cable force vector F during calculatingⁱCorresponding element is provided；Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far the problem of this method realizes rejecting load change and the influence of generalized displacement of support, Cable Structure rope is recognized, generalized displacement of support and the change influence of support cable health status, load change amount identification are rejected while realizing；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, first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then makes the health status of rope for 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. remove once, i.e. i+1 time circulates required current initial Cable Structure bearing generalized coordinates vector Uⁱ⁺¹ _oThe current initial Cable Structure bearing generalized coordinates vector U circulated equal to ithⁱ _o；

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