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

Abstract

The invention provides a method for 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. After that, the mechanical calculation benchmark model of the cable structure, with the generalized displacement of the bearer being taken into account, is obtained. On the basis of the above model, a numerical value variation matrix of a monitored quantity with unit damage is obtained through calculation. Finally, according to the approximately linear relationships between the current numeric vector of the monitored quantity and the current initial numeric vector of the monitored quantity, the numerical value variation matrix of the monitored quantity with unit damage and the current nominal damage vector of an unknown to-be-evaluated object, the non-inferior solution of the current nominal damage vector of the to-be-evaluated object is figured out. Therefore, the health status of a core to-be-evaluated object can be recognized.

Description

Simplify the recognition methods of generalized displacement hybrid monitoring loading problem rope

Technical field

Cable-stayed bridge, suspension bridge, the structures such as truss-frame structure have a common ground, be exactly that they have many parts bearing tensile load, as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or the rod member only bearing tensile load are support unit, for simplicity, such structure representation is " Cable Structure " by this method, and by all ropeway carrying-ropes of Cable Structure, carrying cable, and all rod members (being also called two power rod members) only bearing axial tension or axial compression load, be collectively referred to as simplicity " cable system ", ropeway carrying-rope is censured with " support cable " this noun in this method, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word below, two power rod members are just referred to truss-frame structure reality.Impaired and the lax pair Cable Structure of support cable is safely a significant threat, and damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method.In structure military service process, the correct identification of the health status of support cable or cable system is related to the safety of whole Cable Structure.In Cable Structure military service process, may generalized displacement be there is in Cable Structure bearing, the load that Cable Structure is born also may change, the health status of Cable Structure also may change simultaneously, at this complex condition, based on hybrid monitoring, (this method to judge the health status of Cable Structure to this method by the hybrid monitoring of the change of the measurable parameter to the aforementioned dissimilar Cable Structure of this section, all monitored Cable Structure characteristic parameters are referred to as " monitored amount " by this method, because now monitored amount is made up of the dissimilar measurable parameter mixing of Cable Structure, this method claims this to be hybrid monitoring) carry out identification problem rope, belong to engineering structure health monitoring field.

Background technology

Reject load change, Cable Structure generalized displacement of support to the impact of Cable Structure health status recognition result, thus identify the change of the health status of structure exactly, be problem in the urgent need to address at present, this method discloses a kind of effective, the cheap method addressed this problem.

Summary of the invention

Technical matters: this method discloses a kind of method, under the condition that cost is lower, when bearing has generalized displacement, when the load change that structure is born, generalized displacement of support and load change can be rejected on the impact of Cable Structure health status recognition result, thus identify the health status of support cable exactly.

Technical scheme: this method is made up of three parts.Method, knowledge based storehouse (containing parameter) and actual measurement the structural health conditions appraisal procedure of monitored amount, the software and hardware part of health monitoring systems of setting up knowledge base needed for structural healthy monitoring system and parameter respectively.

In the method, the coordinate of bearing about the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing volume coordinate ", alternatively becoming is the volume coordinate of bearing about X, Y, Z axis, bearing is called the volume coordinate component of bearing about this axle about the concrete numerical value of the volume coordinate of some axles, and a volume coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the volume coordinate of some axles; The angular coordinate of bearing about X, Y, Z axis is censured with " bearing angular coordinate ", bearing is called the angular coordinate component of bearing about this axle about the concrete numerical value of the angular coordinate of some axles, and an angular coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the angular coordinate of some axles; Censure bearing angular coordinate and bearing volume coordinate entirety with " bearing generalized coordinate ", a generalized coordinate component also with bearing in this method expresses the volume coordinate of bearing about an axle or the concrete numerical value of angular coordinate; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, and alternatively the change of bearing volume coordinate is called support wire displacement, and a translational component also with bearing in this method expresses the concrete numerical value of bearing about the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, and an angular displacement component also with bearing in this method expresses the concrete numerical value of bearing about the angular displacement of some axles; Generalized displacement of support censures support wire displacement and angular displacement of support is all, and a generalized displacement component also with bearing in this method expresses the displacement of the lines of bearing about some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is support wire displacement or the translational displacement component at gravity direction.

The Part I of this method: the method setting up knowledge base needed for structural healthy monitoring system and parameter.Specific as follows:

1. the external force that object, structure are born can be described as load, and load comprises face load and volume load.Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds.Volume load be continuous distribution in the load of interior of articles each point, as deadweight and the inertial force of object.

Centre-point load is divided into concentrated force and concentrated couple two kinds, in a coordinate system, such as in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component are called a load, the now change of load is embodied as the change of a concentrated force component or a concentrated couple component.

Distributed load is divided into line distributed load and EDS maps load, the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, the size distribution intensity of distributed load is expressed, distribution intensity distribution characteristics is (such as uniform, sine function equal distribution feature) and amplitude is expressed, and (such as two distributed loads are all uniform, but its amplitude is different, can well-distributed pressure be example so that the concept of amplitude to be described: same structure bears two different well-distributed pressures, two distributed loads are all uniformly distributed loads, but the amplitude of a distributed load is 10MPa, the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, in fact 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, a 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 all not identical, so in the method the component of these several distributed loads is regarded as the independently distributed load of same quantity, now load just represents the component of a distributed load, also component identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized a distributed load or be called a load.

Volume load is that continuous distribution is in the load of interior of articles each point, as deadweight and the inertial force of object, the description of volume load at least comprises the zone of action of volume load and the size of volume load, the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics is (such as uniform, linear function equal distribution feature) and amplitude is expressed, and (such as two individual stow lotuses are all uniform, but its amplitude is different, can conduct oneself with dignity for example is to illustrate the concept of amplitude: the material of two parts of same structure is different, therefore density is different, so although this volume load suffered by two parts is all uniform, but the amplitude of the volume load suffered by a part may be 10kN/m ³, the amplitude of the volume load suffered by another part is 50kN/m ³).If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of 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 stow lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles about coordinate system, that is, in Descartes's rectangular coordinate system, volume load can resolve into three components), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is all not identical, so in the method the component of this several body stow lotus is regarded as the independently load of same quantity, also the volume sharing part of the load identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized an individual stow lotus or be called a load.

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

When load is embodied as distributed load, in the method, " load unit change " in fact 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 " specifically refers to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " specifically refers to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " specifically refers to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " specifically refers to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that changes of amplitude of distribution intensity ", 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 " in fact 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 changes of amplitude of distribution intensity ", briefly, " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.

First the quantity of the load that may change that Cable Structure is born is confirmed.According to the feature of the load that Cable Structure is born, confirm wherein " load likely changed ", or all load is considered as " load likely changed ", if total JZW the load that may change, i.e. total JZW secondary evaluation object.

If the quantity sum of the quantity of the support cable of Cable Structure and JZW " load likely changed " is N, i.e. total N number of evaluation object.To evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.

Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, be described below respectively:

If total Q root support cable in cable system, i.e. total Q core evaluation object, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment rope ₁individual rope force data describes, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment ropes.Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M ₁it is an integer being not less than 0.

The monitored strain data of Cable Structure can by K in Cable Structure ₂the L of individual specified point and each specified point ₂the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K ₂the change of all tested strain of individual specified point.Each total M ₂(M ₂=K ₂× L ₂) individual strain measurement value or calculated value characterize Cable Structure strain.M ₂it is an integer being not less than 0.

The monitored angle-data of Cable Structure is by K in Cable Structure ₃individual specified point, cross the L of each specified point ₃the H of individual appointment straight line, each appointment straight line ₃individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying.Each total M ₃(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 being not less than 0.

The monitored shape data of Cable Structure is by K in Cable Structure ₄the L of individual specified point and each specified point ₄the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K ₄the change of all coordinate components of individual specified point.Each total M ₄(M ₄=K ₄× L ₄) individual coordinates measurements or calculated value characterize Cable Structure shape.M ₄it is an integer being not less than 0.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M (M=M ₁+ M ₂+ M ₃+ M ₄) individual monitored amount, definition parameter K (K=M ₁+ K ₂+ K ₃+ K ₄), M is greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object.

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

2. set up the initial mechanical Calculation Basis model A of Cable Structure _o(such as finite element benchmark model) and current initial mechanical Calculation Basis model A ^t _othe method of (such as finite element benchmark model), sets up and A _ocorresponding monitored amount initial value vector C _omethod, set up and A ^t _ocorresponding monitored amount current initial value vector C ^t _omethod.A in the method _o, C _o, A ^t _oand C ^t _oconstantly update.Set up and upgrade A _o, C _o, A ^t _oand C ^t _omethod as follows.Monitored amount initial value vector C _othe coding rule of coding rule and M monitored amount identical.

Set up initial mechanical Calculation Basis model A _otime, the variable quantity data of the data and " load likely changed " that utilize the Non-destructive Testing Data etc. of support cable can express the health status of support cable set up evaluation object initial damage vector d _o(such as formula (1) Suo Shi), uses d _orepresent that Cable Structure is (with initial mechanical Calculation Basis model A _orepresent) the initial health of evaluation object.If there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d _oin each element numerical value relevant to support cable get 0.Vector d _oin each element numerical value relevant to the variable quantity of load get 0.The physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure to use, utilizes mechanics method (such as finite element method) to set up initial mechanical Calculation Basis model A _o.Corresponding to A _ocable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U _o.

d _o＝[d _o1d _o2· · · d _ok· · · d _oN] ^T(1)

D in formula (1) _ok(k=1,2,3 ...., N) represent initial mechanical Calculation Basis model A _oin the original state of a kth evaluation object, subscript T represents the transposition (afterwards with) of vector.

The initial value of all monitored amount of the Cable Structure using the direct survey calculation of conventional method to obtain before utilization, forms monitored amount initial value vector C _o(see formula (2)).Require at acquisition A _owhile obtain C _o, monitored amount initial value vector C _orepresent and correspond to A _othe concrete numerical value of " monitored amount ".Because of subject to the foregoing, the Calculation Basis model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in describing below, will represent this calculated value and measured value with prosign.

C _o＝[C _o1C _o2· · · C _oj· · · C _oM] ^T(2)

C in formula (2) _oj(j=1,2,3 ...., M) be the original bulk of jth monitored amount in Cable Structure, this component corresponds to a specific jth monitored amount according to coding rule.

No matter which kind of method to obtain initial mechanical Calculation Basis model A by _o, based on A _othe Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A _othe Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, 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 evaluation object initial damage vector d of middle evaluation object _orepresent, model A _othe vectorial U of middle bearing generalized coordinate _orepresent.Due to based on A _othe initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A _obasis on, carry out Mechanics Calculation obtains, A _othe evaluation of each monitored amount form monitored amount initial value vector C _o.U _o, d _oa _oparameter, alternatively C _oby A _omechanics Calculation result composition.

Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ^t _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ^t _o, Cable Structure is in A ^t _oduring state, this method monitored amount current initial value vector C ^t _orepresent the concrete numerical value of all monitored amounts.

Set up and upgrade current initial mechanical Calculation Basis model A ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _o, concrete grammar sees below.

In Cable Structure, the currency of all monitored amounts forms monitored amount current value vector C (formula (3) is shown in definition).

C＝[C ₁C ₂· · · C _j· · · C _M] ^T(3)

C in formula (3) _j(j=1,2,3 ...., M) be the currency of jth monitored amount in Cable Structure, this component C _jaccording to coding rule and C _ojcorresponding to same " monitored amount ".In Cable Structure military service process, constantly actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C.

3. set up and upgrade the method for Cable Structure unit damage monitored numerical quantity transformation matrices Δ C.

Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal ^t _owith monitored amount current initial value vector C ^t _owhile, upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C.Concrete grammar is as follows:

At the current initial mechanical Calculation Basis model A of Cable Structure ^t _obasis on carry out several times calculating, calculation times numerically equals the quantity of all evaluation objects.Calculating hypothesis each time only has an evaluation object (to use vectorial d at initial damage _ocorresponding element represent) basis on increase again unit damage or load unit change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable has unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is a load, just suppose that this load is at vectorial d _othe basis that this load represented has a variable quantity increases again load unit change (if this load is distributed load, and this distributed load is line distributed load, load unit change can get 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. for unit change; If this load is distributed load, and this distributed load is EDS maps load, and load unit change can get 1MPa, 2MPa, 3MPa or 1kNm/m ², 2kNm/m ², 3kNm/m ²deng be unit change; If this load is centre-point load, and this centre-point load is couple, and load unit change can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centre-point load, and this centre-point load is concentrated force, and load unit change can get 1kN, 2kN, 3kN etc. for unit change; If this load is volume load, load unit change can get 1kN/m ³, 2kN/m ³, 3kN/m ³deng be unit change), use D _ukrecord this unit damage or load unit change, wherein k represents the numbering of the evaluation object that unit damage or load unit change occur.With " evaluation object unit change vector D _u" (such as formula (4) Suo Shi) records all unit damage or load unit changes.Occur in calculating each time that the evaluation object of unit damage or load unit change is different from during other time calculates the evaluation object occurring unit damage or load unit change, calculate the current calculated value all utilizing mechanics method (such as finite element method) to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, and (when supposing that a kth evaluation object has unit damage or load unit to change, available formula (5) represents monitored amount calculation current vector C _t ^k); Calculate monitored amount calculation current vector C each time _t ^kdeduct monitored amount current initial value vector C ^t _oafter calculate the unit damage supposed or load unit change numerical value D divided by this time again _uk, gained vector is exactly that the monitored amount unit change vector of (being numbered mark with what have an evaluation object of unit damage or load unit change) (when a kth evaluation object has unit damage or load unit changes, uses δ C under this condition _krepresent monitored amount unit change vector, formula (6) is shown in definition), the Unit alteration amount of the monitored amount corresponding to this element that each element representation of monitored amount unit change vector causes owing to suppose there is the unit damage of that evaluation object of unit damage or load unit change or load unit change when calculating; N number of evaluation object is had just to have N number of monitored amount unit change vector, owing to there being M monitored amount, so each monitored amount unit change vector has M element, the definition of monitored amount unit change matrix Δ C, the Δ C that have M × N number of element is made up of successively such as formula shown in (6) this N number of monitored amount unit change vector.

D _u＝[D _u1D _u2· · · D _uk· · · D _uN] ^T(4)

Evaluation object unit change vector D in formula (4) _uelement D _uk(k=1,2,3 ...., N) represent that the unit damage of a kth evaluation object or load unit change numerical value.

$C_t^k={\left[\begin{array}{cccccccccc}{C}_{t1}^k& C_{t2}^k& \·& \·& \·& C_{tj}^k& \·& \·& \·& C_{tM}^k\end{array}\right]}^T---\left(5\right)$

Elements C in formula (5) _tj ^k(k=1,2,3 ...., N; J=1,2,3 ...., M) represent due to a kth evaluation object have unit damage or a load unit change time, according to the current calculated amount of the individual monitored amount of the jth corresponding to coding rule.

${\mathrm{\δC}}_k=\frac{C_t^k-C_o^t}{D_{uk}}---\left(6\right)$

$\mathrm{\Δ}C=\left[\begin{array}{cccccc}{\mathrm{\ΔC}}_{1,1}& {\mathrm{\ΔC}}_{1,2}& \·& {\mathrm{\ΔC}}_{1,k}& \·& {\mathrm{\ΔC}}_{1,N}\\ {\mathrm{\ΔC}}_{2,1}& {\mathrm{\ΔC}}_{2,2}& \·& {\mathrm{\ΔC}}_{2,k}& \·& {\mathrm{\ΔC}}_{2,N}\\ \·& \·& \·& \·& \·& \·\\ {\mathrm{\ΔC}}_{j,1}& {\mathrm{\ΔC}}_{j,2}& \·& {\mathrm{\ΔC}}_{j,k}& \·& {\mathrm{\ΔC}}_{j,N}\\ \·& \·& \·& \·& \·& \·\\ {\mathrm{\ΔC}}_{M,1}& {\mathrm{\ΔC}}_{M,2}& \·& {\mathrm{\ΔC}}_{M,k}& \·& {\mathrm{\ΔC}}_{M,N}\end{array}\right]---\left(7\right)$

Δ C in formula (7) _j,k(k=1,2,3 ...., N; J=1,2,3 ...., M) represent only due to a kth evaluation object have unit damage or load unit change to cause, according to the unit change (algebraic value) of the calculating current value of the monitored amount of the jth corresponding to coding rule, the vectorial δ C of monitored amount unit change _kbe actually the row in matrix Δ C.

4. monitored amount current value vector C (calculating or actual measurement) is with monitored amount current initial value vector C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D _uand the linear approximate relationship between evaluation object current nominal fatigue vector d, shown in (8) or formula (9).The definition of evaluation object current nominal fatigue vector d is see formula (10).

$C=C_o^t+\mathrm{\Δ}C\·d---\left(8\right)$

$C-C_o^t=\mathrm{\Δ}C\·d---\left(9\right)$

d＝[d ₁d ₂· · · d _k· · · d _N] ^T(10)

D in formula (10) _k(k=1,2,3 ...., N) be the current health state of a kth evaluation object in Cable Structure, if this evaluation object is a rope (or pull bar) in cable system, so d _krepresent its current damage, d _krepresent not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, represent the load-bearing capacity losing corresponding proportion time between 0 and 100%, if this evaluation object is load, so a d _krepresent its variable quantity.

The error of the linear relationship error vector e expression (8) that available formula (11) defines or the shown linear relationship of formula (9).

$e=abs(\mathrm{\Δ}C\·d-C+C_o^t)---\left(11\right)$

In formula (11), abs () is the function that takes absolute value, and takes absolute value to each element of the vector of trying to achieve in bracket.

The Part II of this method: the structural health conditions appraisal procedure of knowledge based storehouse (containing parameter) and the monitored amount of actual measurement.

There is certain error in the linear relationship represented by formula (8) or formula (9), therefore simply can not carry out direct solution according to formula (8) or formula (9) and actual measurement monitored amount current value vector C and obtain evaluation object current nominal fatigue vector d.If this has been doned, the element in the evaluation object obtained current nominal fatigue vector d even there will be larger negative value, namely negative damage, and this is obviously irrational.Therefore the acceptable solution of evaluation object current nominal fatigue vector d is obtained (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more accurately from cable system) become a rational solution, available formula (12) expresses this method.

$abs(\mathrm{\Δ}C\·d-C+C_o^t)\≤g---\left(12\right)$

In formula (12), abs () is the function that takes absolute value, and vectorial g describes the legitimate skew departing from ideal linearity relation (formula (8) or formula (9)), is defined by formula (13).

g＝[g ₁g ₂· · · g _j· · · g _M] ^T(13)

G in formula (13) _j(j=1,2,3 ...., M) describe the maximum allowable offset departing from formula (8) or the ideal linearity relation shown in formula (9).The error vector e tentative calculation that vector g can define according to formula (11) is selected.

At monitored amount current initial value vector C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, survey monitored amount current value vector C known time, suitable algorithm (such as multi-objective optimization algorithm) can be utilized to solve formula (12), obtain the acceptable solution of evaluation object current nominal fatigue vector d.

Definition evaluation object current actual damage vector d ^a(see formula (14)), can by d ^adetermine the health status of evaluation object.

$d^a={\left[\begin{array}{cccccccccc}{d}_1^a& d_2^a& \·& \·& \·& d_k^a& \·& \·& \·& d_N^a\end{array}\right]}^T---\left(14\right)$

D in formula (14) ^a _k(k=1,2,3,., N) represent eliminate Cable Structure generalized displacement of support and load change on after the impact of health status recognition result, the current actual health status of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d ^a _krepresent its current actual health status, its definition is shown in formula (15), d ^a _krepresent not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; If this evaluation object is a load, its definition is shown in formula (15), so d ^a _krepresent that it is relative to setting up initial mechanical Calculation Basis model A _otime the structure variable quantity of corresponding load that bears.Vector d ^athe coding rule of element and formula (1) in vectorial d _othe coding rule of element identical.

D in formula (15) _ok(k=1,2,3 ...., N) be vectorial d _oa kth element, d _kit is a kth element of vectorial d.

Describe below and obtain evaluation object current actual damage vector d ^aafter, how to identify damaged cable and slack line.

Total Q root support cable in Cable Structure, Cable Structure rope force data is described by the Suo Li of Q root support cable.Available " Initial cable force vector F _o" represent the Initial cable force (formula (16) is shown in definition) of all support cables in Cable Structure.

F _o＝[F _o1F _o2· · · F _oh· · · F _oQ] ^T(16)

F in formula (16) _o(h=1,2,3 ...., Q) be the Initial cable force of h root support cable in Cable Structure, this element corresponds to the Suo Li specifying support cable according to coding rule.Vector F _oit is constant.Setting up the initial mechanical Calculation Basis model A of Cable Structure _osynchronization, use conventional method direct survey calculation to obtain the rope force data of all support cables, all these rope force datas composition Initial cable force vector F _o.Setting up the initial mechanical Calculation Basis model A of Cable Structure _otime in fact employ vectorial F _o.

By current for evaluation object actual damage vector d ^ain relevant to support cable Q element take out, composition support cable current actual damage vector d ^ca, support cable current actual damage vector d ^cathe coding rule of element and Initial cable force vector F _othe coding rule of element identical.Support cable current actual damage vector d ^cah element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., Q; Current actual damage vector d ^camiddle numerical value be not 0 element correspond to the support cable of unsoundness problem, carry out Non-Destructive Testing to the support cable of these unsoundness problems, after Non-Destructive Testing finds out that this support cable is not damaged, so this element numerical value (uses d ^ca _hrepresent) represent this support cable and d ^ca _hrelaxing of impairment value mechanic equivalent, just determine slack line thus, the computing method of concrete slack are described below.

Represent the current cable power (formula (17) is shown in definition) of all support cables in the Cable Structure of surveying and obtaining with " current cable force vector F " in this method.

F＝[F ₁F ₂· · · F _h· · · F _Q] ^T(17)

F in formula (17) _h(h=1,2,3 ...., Q) be the current cable power of h root support cable in Cable Structure.Obtain the synchronization of monitored amount current value vector C in actual measurement, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F.The element of vector F and vectorial F _othe coding rule of element identical.

In this method, under support cable original state, (namely Cable Structure is in initial mechanical Calculation Basis model A _oduring the state expressed), when support cable is in free state (free state refers to that Suo Li is 0, rear same), the length of support cable is called initial drift, with " initial drift vector l _o" represent the initial drift (formula (18) is shown in definition) of all support cables in Cable Structure.

l _o＝[l _o1l _o2· · · l _oh· · · l _oQ] ^T(18)

L in formula (18) _oh(h=1,2,3 ...., Q) be the initial drift of h root support cable in Cable Structure.Vector l _obeing constant, after determining when starting, just no longer changing.

Similar, under support cable original state, (namely Cable Structure is in initial mechanical Calculation Basis model A _oduring the state expressed), and when support cable is in free state, the cross-sectional area of support cable is called initial free cross-sectional area, with " initial free cross-sectional area vector A _o" representing the initial free cross-sectional area (formula (19) is shown in definition) of all support cables in Cable Structure, the weight of the unit length of support cable is called the weight of initial free unit length, with " the weight vector ω of initial free unit length _o" represent the weight (formula (20) is shown in definition) of the initial free unit length of all support cables in Cable Structure.

A _o＝[A _o1A _o2· · · A _oh· · · A _oQ] ^T(19)

A in formula (19) _oh(h=1,2,3 ...., Q) be the initial free cross-sectional area of h root support cable in Cable Structure.Vector A _obeing constant, after determining when starting, just no longer changing.

ω _o＝[ω _o1ω _o2· · · ω _oh· · · ω _oQ] ^T(20)

ω in formula (20) _oh(h=1,2,3 ...., Q) be the weight of the free unit length of initial freedom of h root support cable in Cable Structure.Vector ω _obeing constant, after determining when starting, just no longer changing.

Vector d ^caelement, the element of vectorial F, vectorial l _oelement, vectorial A _oelement, vectorial ω _oelement and vectorial F _othe coding rule of element identical, the different information of the same support cable of element representation of the identical numbering of these vectors.

In this method, the synchronization of monitored amount current value vector C is obtained in actual measurement, in Cable Structure, with " current drift vector l ", the current drift of all support cables represents that (formula (21) is shown in definition, now support cable may be intact, also may be impaired, also may be lax).

l＝[l ₁l ₂· · · l _h· · · l _Q] ^T(21)

L in formula (21) _h(h=1,2,3 ...., Q) be the current drift of h root support cable in Cable Structure.

In this method, represent the knots modification (formula (22) and formula (23) are shown in definition) of the drift of all support cables in Cable Structure with " drift changes vectorial Δ l " (or claiming support cable current slack degree vector).

Δl＝[Δl ₁Δl ₂· · · Δl _h· · · Δl _Q] ^T(22)

Δ l in formula (22) _h(h=1,2,3 ...., Q) be the knots modification of the drift of h root support cable in current cable structure, its definition is shown in formula (23), Δ l _hbe not 0 rope be slack line, Δ l _hnumerical value be the slack of rope, and representing the current slack degree of cable system h root support cable, is also the long adjustment amount of rope of this rope during adjustment Suo Li.

Δl _h＝l _h-l _oh(23)

In the method by slack line is carried out with damaged cable the relax level identification that mechanic equivalent carries out slack line, mechanic equivalent condition is:

One, the nothing of the rope of two equivalences is lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter and material;

Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical.

When meeting above-mentioned two mechanic equivalent conditions, such two support cables mechanics function is in the structure exactly identical, if after namely replacing slack line with the damaged cable of equivalence, any change can not occur Cable Structure, and vice versa.

Obtain support cable current actual damage vector d ^caafter, d ^cah element d ^ca _h(h=1,2,3 ...., Q) represent the actual damage value of h root support cable, although by d ^ca _hbe called the actual damage value of h root rope or the actual damage degree of h root rope, but also may be lax because h root support cable may be impaired, so d ^cah element d ^ca _hthe actual damage value of the h root support cable represented is actually the actual equivalent damage value of h root support cable, when h root support cable is actually impaired, and d ^ca _hthe actual damage value of the h root support cable just represented, when h root support cable is actually lax, d ^a _hthe h root support cable just represented with the actual damage value of lax equivalence, for sake of convenience, claim d in the method ^a _hrepresent h root support cable not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, time between 0 and 100%, represent that h root support cable loses the load-bearing capacity of corresponding proportion, by support cable current actual damage vector d ^cajust can identify the support cable that health status goes wrong afterwards, but in the support cable that goes wrong of these health status, some is impaired, some relaxes, if h support cable is actually relax (its current slack degree Δ l _hdefinition), the current slack degree Δ l of h so lax support cable _h(Δ l _hdefinition see formula (23)) with the current actual damage degree d of damaged cable of equivalence ^ca _hbetween relation determined by aforementioned two mechanic equivalent conditions.Δ l _hsame d ^ca _hbetween physical relationship can adopt accomplished in many ways, such as can directly determine (see formula (24)) according to aforementioned equivalent condition, also can adopt after replacing the E in formula (24) to revise based on Ernst equivalent elastic modulus and determine (see formula (25)), other method such as trial and error procedure based on finite element method also can be adopted to determine.

${\mathrm{\Δl}}_h=\frac{d_h^{ca}}{1-d_h^{ca}}\frac{F_h}{E_{oh}{A}_{oh}+F_h}{l}_{oh}---\left(24\right)$

${\mathrm{\Δl}}_h=\frac{d_h^{ca}}{1-d_h^{ca}}\frac{F_h}{\[\frac{E_{oh}}{1+\frac{{\left({\mathrm{\ω}}_{oh}{l}_{xh}\right)}^2{A}_{oh}{E}_{oh}}{12{\left(F_h\right)}^3}}\]A_{oh}+F_h}{l}_{oh}---\left(25\right)$

Formula (24) and the middle E of formula (25) _ohvectorial E _oh element, be setting up the initial mechanical Calculation Basis model A of Cable Structure _osynchronization, the elastic modulus of h root support cable, A _ohsetting up the initial mechanical Calculation Basis model A of Cable Structure _osynchronization, the cross-sectional area of h root support cable, ω _ohsetting up the initial mechanical Calculation Basis model A of Cable Structure _osynchronization, the weight of the unit length of h root support cable, F _hthe synchronization obtaining monitored amount current value vector C in actual measurement, the current cable power of h root support cable, d ^ca _hthe current actual damage degree of h root support cable, l _xhthe horizontal range of the synchronization obtaining monitored amount current value vector C in actual measurement, two supporting end points of h root support cable, l _xhthat current support cable two supports end points horizontal range vector l _xan element, current support cable two support end points horizontal range vector l _xthe coding rule of element and initial drift vector l _othe coding rule of element identical, E _hcan obtain according to the characteristic material data looking into or survey h root support cable, A _hand ω _hcan according to A _oh, ω _oh, F _hand T _oobtained by Typical physical and Mechanics Calculation.Item in formula (25) in [] is the Ernst equivalent elastic modulus of this support cable, just can determine support cable current slack degree vector Δ l by formula (24) or formula (25).Formula (25) is the correction to formula (24).

So far this method achieves the accurate identification of the health status of core evaluation object with a kind of effective, cheap method.May exact value be departed to the recognition result of the health status of secondary evaluation object more, only require the correct health status identifying core evaluation object in the method.

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

Hardware components comprises monitoring system (comprising the space coordinate monitoring system of supporting end points of monitored amount monitoring system, Cable Structure bearing generalized coordinate monitoring system, support cable cable force monitoring system, support cable), signal picker and computing machine etc.Require that the volume coordinate of each bearing generalized coordinate of each monitored amount, Cable Structure, the Suo Li of each root support cable, the supporting end points of each root support cable must monitored system real-time monitor.

Software should to complete in this method required, can by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.

This method specifically comprises:

Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated support cable and load to be " evaluation object ", if the quantity sum of the quantity of evaluated support cable and load is N, namely the quantity of " evaluation object " is N; This method title " core evaluation object " specially refers to the evaluated support cable in " evaluation object ", and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; Specify when determining hybrid monitoring by the support cable of monitored Suo Li, if total Q root support cable in cable system, i.e. total Q core evaluation object, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment support cable ₁individual rope force data describes, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment support cables; Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M ₁it is an integer being not less than that 0 is not more than Q; Specify when determining hybrid monitoring by the measured point of monitored strain, the monitored strain data of Cable Structure can by K in Cable Structure ₂the L of individual specified point and each specified point ₂the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K ₂the change of all tested strain of individual specified point; Each total M ₂individual strain measurement value or calculated value characterize Cable Structure strain, M ₂for K ₂and L ₂long-pending; M ₂be be not less than 0 integer; Specify when determining hybrid monitoring by the measured point of monitored angle, the monitored angle-data of Cable Structure is by K in Cable Structure ₃individual specified point, cross the L of each specified point ₃the H of individual appointment straight line, each appointment straight line ₃individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying; Each total M ₃individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M ₃for K ₃, L ₃and H ₃long-pending; M ₃it is an integer being not less than 0; Specify when determining hybrid monitoring by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure ₄the L of individual specified point and each specified point ₄the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K ₄the change of all coordinate components of individual specified point; Each total M ₄individual coordinates measurements or calculated value characterize Cable Structure shape, M ₄for K ₄and L ₄long-pending; M ₄it is an integer being not less than 0; The monitored amount of comprehensive above-mentioned hybrid monitoring, whole Cable Structure has M monitored amount, and M is M ₁, M ₂, M ₃and M ₄sum, definition parameter K, K is M ₁, K ₂, K ₃and K ₄sum, M must be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; For simplicity, in the method by the monitored amount of the M listed by this step referred to as " monitored amount "; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centre-point load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated 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 comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact 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; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the one-component of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;

B. survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses;

C. actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the measured data of initial Cable Structure, the measured data of initial Cable Structure comprises Cable Structure centre-point 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 coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Non-destructive Testing Data comprising support cable, the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C _o, monitored amount initial value vector C _othe coding rule of coding rule and M monitored amount identical, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d _o, vectorial d _orepresent with initial mechanical Calculation Basis model A _othe initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d _oelement number equal N, d _oelement and evaluation object be one-to-one relationship, vectorial d _othe coding rule of element identical with the coding rule of evaluation object, if d _oevaluation object corresponding to some elements be support cable, so a d in cable system _othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses load-bearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the load-bearing capacity of corresponding proportion, if d _oevaluation object corresponding to some elements be some load, get d in this method _othis element numerical value be 0, the initial value representing the change of this load is 0, if there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d _oin each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U _o, actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F _o, according to comprising Cable Structure design data, the data of completion data obtain the length of all support cables when free state and Suo Li are 0, in free state time cross-sectional area and in free state time the weight of unit length, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum and Initial cable force vector F _othe coding rule of element identical,

Physical and mechanical properties parameter, the initial Cable Structure bearing generalized coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, support cable initial health data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure _owith all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A of Cable Structure _o, based on A _othe Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A _oevaluation object health status with evaluation object initial damage vector d _orepresent; Corresponding to A _othe initial value monitored amount initial value vector C of all monitored amount _orepresent; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _owith monitored amount current initial value vector C ^t _o; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _owith monitored amount current initial value vector C ^t _otime, the current initial mechanical Calculation Basis model A of Cable Structure ^t _ojust equal the initial mechanical Calculation Basis model A of Cable Structure _o, monitored amount current initial value vector C ^t _ojust equal monitored amount initial value vector C _o; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ^t _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ^t _o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _otime, U ^t _ojust equal U _o; A ^t _othe initial health of evaluation object and A _othe health status of evaluation object identical, also use evaluation object initial damage vector d _orepresent, A in cyclic process below ^t _othe initial health of evaluation object use evaluation object initial damage vector d all the time _orepresent; U _oand d _oa _oparameter, by A _othe initial value of all monitored amount that obtains of Mechanics Calculation result and C _othe initial value of all monitored amount represented is identical, therefore alternatively C _oby A _omechanics Calculation result composition; U ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition;

E. from entering the circulation being walked to m step by e here; In structure military service process, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, all Cable Structure bearing generalized coordinate current data composition current cable structure actual measurement bearing generalized coordinate vector U ^t, vectorial U ^tdefinition mode and vectorial U _odefinition mode identical;

F. according to current cable structure actual measurement bearing generalized coordinate vector U ^t, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3 ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _o;

F1. U is compared ^twith U ^t _oif, U ^tequal U ^t _o, then A ^t _o, U ^t _oand C ^t _oremain unchanged, otherwise need to follow these steps to A ^t _o, U ^t _oand C ^t _oupgrade;

F2. U is calculated ^twith U _odifference, U ^twith U _odifference be exactly the generalized displacement of support of Cable Structure bearing about initial position, with generalized displacement of support vector V represent generalized displacement of support, V equals U ^tdeduct U _o, be one-to-one relationship between the element in generalized displacement of support vector V and generalized displacement of support component, in generalized displacement of support vector V, the numerical value of an element corresponds to the generalized displacement of an assigned direction of an appointment bearing;

F3. to A _oin Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint just takes from the numerical value of corresponding element in generalized displacement of support vector V, to A _omiddle Cable Structure bearing obtains the current initial mechanical Calculation Basis model A upgraded after applying generalized displacement of support constraint ^t _o, upgrade A ^t _owhile, U ^t _oall elements numerical value also uses U ^tall elements numerical value correspondence replaces, and namely have updated U ^t _o, so just obtain and correctly correspond to A ^t _ou ^t _o; Upgrade C ^t _omethod be: when renewal A ^t _oafter, obtain A by Mechanics Calculation ^t _oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C ^t _o; A ^t _othe initial health of support cable use evaluation object initial damage vector d all the time _orepresent;

G. at current initial mechanical Calculation Basis model A ^t _obasis on carry out several times Mechanics Calculation according to step g 1 to g4, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D _u;

G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _oafterwards, Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D must then be upgraded _u;

G2. at the current initial mechanical Calculation Basis model A of Cable Structure ^t _obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object on the basis of original damage or load, increase unit damage or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just supposes that this support cable is at vectorial d _othe basis that this support cable represented has a damage increases unit damage again, if this evaluation object is a load, just supposes that this load is at vectorial d _othe basis that this load represented has a variable quantity increases load unit change again, use D _ukrecord unit damage or the load unit change of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or load unit change, D _ukevaluation object unit change vector D _uan element, evaluation object unit change vector D _uthe coding rule of element and vectorial d _othe coding rule of element identical; The evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector _oelement number rule identical;

G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C ^t _oobtain a vector, then each element of this vector is calculated the unit damage or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;

G4. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d _othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;

H. current cable structure actual measurement bearing generalized coordinate vector U is obtained in actual measurement ^twhile, actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C ^t _owith monitored amount initial value vector C _odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time; Obtain the synchronization of monitored amount current value vector C in actual measurement, actual measurement obtains the rope force data of all Q root support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F _othe coding rule of element identical; The synchronization of monitored amount current value vector C is obtained in actual measurement, Actual measurement obtains the volume coordinate of two supporting end points of all Q root support cables, the difference of the volume 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 form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector _othe coding rule of element identical;

I. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be one-to-one relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object; The coding rule of the element of vector d and vectorial d _othe coding rule of element identical;

J. according to monitored amount current value vector C with monitored amount current initial value vector C ^t _o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;

$C=C_o^t+\mathrm{\Δ}C\·d$ Formula 1

K. evaluation object current actual damage vector d is defined ^a, evaluation object current actual damage vector d ^aelement number equal the quantity of evaluation object, evaluation object current actual damage vector d ^aelement and evaluation object between be one-to-one relationship, evaluation object current actual damage vector d ^aelement numerical value represent actual damage degree or the real load variable quantity of corresponding evaluation object; Vector d ^athe coding rule of element and vectorial d _othe coding rule of element identical;

L. the evaluation object utilizing formula 2 to express current actual damage vector d ^aa kth element d ^a _kwith evaluation object initial damage vector d _oa kth element d _okwith a kth element d of evaluation object current nominal fatigue vector d _kbetween relation, calculate evaluation object current actual damage vector d ^aall elements;

K=1 in formula 2,2,3 ...., N, d ^a _krepresent the current actual health status of a kth evaluation object, d ^a _krepresent when being 0 that a kth evaluation object is without health problem, d ^a _knumerical value represents when not being 0 that a kth evaluation object is the evaluation object of unsoundness problem, if this evaluation object is support cable, so a d in cable system ^a _krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d ^a _kthe degree of the lax or damage of this support cable of numerical response; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, evaluation object current actual damage vector d ^ain correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent; If this evaluation object is load, so a d ^a _krepresent the actual change amount of this load; After identifying slack line, utilize evaluation object current actual damage vector d ^athese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize the current cable force vector F that obtains in f step and current support cable two to support end points horizontal range vector, utilize the initial drift vector of the support cable obtained in b step, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the vectorial F of Initial cable force _othe physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in b step to use, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanic equivalent condition is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting above-mentioned two mechanic equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the non-destructive tests of support cable; Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so according to evaluation object current actual damage vector d ^athe health status of core evaluation object can be determined;

M. get back to e step, start the circulation next time being walked to m step by e.

Beneficial effect: structural healthy monitoring system is first by using sensor to carry out long-term on-line monitoring to structural response, after obtaining Monitoring Data, (or off-line) analysis is online carried out to it and obtain structural health conditions data, due to the complicacy of structure, structural healthy monitoring system needs to use a large amount of sensor equipment to carry out monitoring structural health conditions, therefore its cost is usually quite high, and 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 ingredient of the correct identification of structural health conditions, or even they are whole, and the impact of correct identification on the correct identification of the health status of Cable Structure of the change (such as by the change of the quality and quantity of the automobile of cable-stayed bridge) of secondary evaluation object (load that such as structure is born) is very little, or even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is normally suitable, the quantity of secondary evaluation object is also usually greater than the quantity of core evaluation object, and the quantity of such evaluation object is usually many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, conventional method requires that the quantity of monitored amount (using sensor device measuring to obtain) must be more than or equal to the quantity of evaluation object, when the number ratio of the secondary evaluation object changed is larger (in fact often so), the quantity of the sensor equipment required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, unacceptablely even high.Inventor studies discovery, in the secondary evaluation object (normal load that such as structure is born, the normal load of structure refers to that the load that structure is being born is no more than the structure allowable load limited according to structural design book or structure completion book) change less time (be exactly that structure only bears normal load for load, whether the load that structure is born is normal load, can be observed by methods such as naked eyes and determine, if find that the load that structure is born is not normal load, so artificially remove, after removing improper load, structure just only bears normal load), the amplitude of variation (this instructions is called " secondary response ") of the structural response caused by them much smaller than core evaluation object change (such as support cable is impaired) caused by the amplitude of variation (this instructions is called " core response ") of structural response, secondary response and core respond total change (this instructions is called " global response ") that sum is structural response, obvious core response dominate in global response, based on this, find to choose when determining monitored amount quantity to be a bit larger tham core evaluation object quantity even if inventor studies, but much smaller than the numerical value (this method is exactly do like this) of evaluation object quantity, even if that is adopt the relatively few a lot of sensor equipment 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 this method cost of structural healthy monitoring system of advising is more much lower than the cost of the structural healthy monitoring system required by conventional method apparently, that is this method can realize to the health status of the core evaluation object of Cable Structure with the much lower condition of cost assessment, can this benefit be used structural health monitoring technology is very important.

Embodiment

The following describes of embodiment of this method is in fact only exemplary, and object is never the application or the use that limit this method.

The first step: the quantity first confirming the load that may change that Cable Structure is born.According to the feature of the load that Cable Structure is born, confirm wherein " load likely changed ", or all load is considered as " load likely changed ", if total JZW the load that may change, i.e. total JZW secondary evaluation object.

If the quantity sum of the quantity of the support cable of Cable Structure and JZW " load likely changed " is N, i.e. total N number of evaluation object.To evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.

Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, be described below respectively:

If total Q root support cable in cable system, i.e. total Q core evaluation object, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment rope ₁individual rope force data describes, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment ropes.Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M ₁it is an integer being not less than 0.

The monitored strain data of Cable Structure can by K in Cable Structure ₂the L of individual specified point and each specified point ₂the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K ₂the change of all tested strain of individual specified point.Each total M ₂(M ₂=K ₂× L ₂) individual strain measurement value or calculated value characterize Cable Structure strain.M ₂it is an integer being not less than 0.

The monitored angle-data of Cable Structure is by K in Cable Structure ₃individual specified point, cross the L of each specified point ₃the H of individual appointment straight line, each appointment straight line ₃individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying.Each total M ₃(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 being not less than 0.

The monitored shape data of Cable Structure is by K in Cable Structure ₄the L of individual specified point and each specified point ₄the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K ₄the change of all coordinate components of individual specified point.Each total M ₄(M ₄=K ₄× L ₄) individual coordinates measurements or calculated value characterize Cable Structure shape.M ₄it is an integer being not less than 0.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M (M=M ₁+ M ₂+ M ₃+ M ₄) individual monitored amount, definition parameter K (K=M ₁+ K ₂+ K ₃+ K ₄), the quantity that M must not be less than core evaluation object adds the quantity that 4, M is less than evaluation object.

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

Second step: set up initial mechanical Calculation Basis model A _o.

When Cable Structure is completed, or before setting up health monitoring systems, conventional method (consult reference materials or survey) is used to obtain physical parameter (such as density) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses, use the direct survey calculation of conventional method to obtain the initial value of all monitored amount of Cable Structure simultaneously, form monitored amount initial value vector C _o.

Monitored amount initial value vector C is obtained at Actual measurement _osynchronization, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinate data, Cable Structure centre-point 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, structure space measurement of coordinates data that the Non-destructive Testing Data of support cable etc. can express the health status of rope.Initial Cable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U _o.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is the geometric properties according to these coordinate data determination Cable Structure.For cable-stayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, so-called bridge type data that Here it is.Utilize the Non-destructive Testing Data etc. of support cable can express the data of the health status of support cable and Cable Structure load measurement data set up evaluation object initial damage vector d _o(such as formula (1) Suo Shi), uses d _orepresent that Cable Structure is (with initial mechanical Calculation Basis model A _orepresent) the initial health of evaluation object.If there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d _oin each element numerical value relevant to support cable get 0; If d _oevaluation object corresponding to some elements be some load, get d in this method _othis element numerical value be 0, the initial value representing the change of this load is 0.The physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure to use and initial Cable Structure bearing generalized coordinate vector U _o, utilize mechanics method (such as finite element method) to set up initial mechanical Calculation Basis model A _o.

Actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, utilize the elastic modulus data composition support cable initial elastic modulus vector E of support cable wherein _o; Obtaining monitored amount initial value vector C _osynchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F _o; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time cross-sectional area and in free state time the weight of unit length, form the initial drift vector l of support cable successively _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 _othe coding rule of element and Initial cable force vector F _othe coding rule of element identical.

3rd step: first time sets up current initial mechanical Calculation Basis model A ^t _o, monitored amount current initial value vector C ^t _o" current initial Cable Structure bearing generalized coordinate vector U ^t _o", concrete grammar is: at initial time, and namely first time sets up current initial mechanical Calculation Basis model A ^t _owith monitored amount current initial value vector C ^t _otime, A ^t _ojust equal A _o, C ^t _ojust equal C _o, corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ^t _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ^t _o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _otime, U ^t _ojust equal U _o.A ^t _othe health status of evaluation object and A _oevaluation object health status (evaluation object initial damage vector d _orepresent) identical, A in cyclic process ^t _othe health status of evaluation object use evaluation object initial damage vector d all the time _orepresent.U ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition.

4th step: in Cable Structure military service process, continuous actual measurement obtains the current measured value of all monitored amount of Cable Structure, composition " monitored amount current value vector C ", actual measurement simultaneously obtains Cable Structure bearing generalized coordinate current data, all data composition current cable structure actual measurement bearing generalized coordinate vector U ^t, vectorial U ^tdefinition mode and vectorial U _odefinition mode identical.

Obtaining the synchronization of monitored amount current value vector C, actual measurement obtains the rope force data of all Q root support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F _othe coding rule of element identical; Obtaining the synchronization of monitored amount current value vector C, Actual measurement obtains the volume coordinate of two supporting end points of all Q root support cables, the difference of the volume 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 root support cables form current support cable two and support end points horizontal range vector l _x, current support cable two supports end points horizontal range vector l _xthe coding rule of element and Initial cable force vector F _othe coding rule of element identical.

5th step: according to current cable structure actual measurement bearing generalized coordinate vector U ^t, upgrade current initial mechanical Calculation Basis model A where necessary ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _o.Current cable structure actual measurement bearing generalized coordinate vector U is obtained in the 4th step actual measurement ^tafter, compare U ^tand U ^t _oif, U ^tequal U ^t _o, then do not need A ^t _o, U ^t _oand C ^t _oupgrade, otherwise need A ^t _o, U ^t _oand C ^t _oupgrade, update method is undertaken by the method specified in technical scheme.

6th step: at current initial mechanical Calculation Basis model A ^t _obasis on carry out several times Mechanics Calculation, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D _u.Concrete grammar described in technical scheme and claims.

7th step: set up linear relationship error vector e and vectorial g.Utilize (the monitored amount current initial value vector C of data above ^t _o, unit damage monitored numerical quantity transformation matrices Δ C), while the 6th step calculates each time, namely only have the increase unit damage of an evaluation object or load unit change D calculating each time in hypothesis evaluation object _ukthe evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current value all utilizing mechanics method (such as adopting finite element method) to calculate all monitored amounts in Cable Structure each time, while calculating the monitored amount calculation current vector C of composition one each time, calculate composition injury vector d each time, originally walk out of existing injury vector d only to use in this step, in all elements of injury vector d, only have the numerical value of an element to get D _uk, the numerical value of other element gets 0, the coding rule of the element of injury vector d and vectorial d _othe coding rule of element identical; By C, C ^t _o, Δ C, D _u, d brings formula (12) into, obtain a linear relationship error vector e, calculate a linear relationship error vector e each time; N number of evaluation object is had just to have N calculating, just there is N number of linear relationship error vector e, obtaining a vector after being added by this N number of linear relationship error vector e, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N.Vector g equals final error vector e.

8th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system is (such as containing measurement of angle subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate measures subsystem, signal conditioner etc.), Cable Structure bearing generalized coordinate monitoring system is (containing total powerstation, angle measuring sensor, signal conditioner etc.), support cable cable force monitoring system, the space coordinate monitoring system of the supporting end points of support cable, signal (data) collector, computing machine and communication alert equipment.The volume coordinate of the Suo Li of each monitored amount, each root support cable, the supporting end points of each root support cable, the bearing generalized coordinate of each Cable Structure must arrive by monitored system monitoring, monitoring system by the Signal transmissions that monitors to signal (data) collector; Signal is delivered to computing machine through signal picker; The health monitoring software of the evaluation object running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring evaluation object health status and changing, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.

9th step: by current for monitored amount initial value vector C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D _uparameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.

Tenth step: establishment installation and operation this method system software on computers, the function (i.e. all work that can complete with computing machine in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method required by task and wants

11 step: according to monitored amount current value vector C with monitored amount current initial value vector C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D _uand the linear approximate relationship (formula (8)) existed between evaluation object current nominal fatigue vector d (being made up of all Suo Dangqian nominal fatigue amounts), calculate the noninferior solution of evaluation object current nominal fatigue vector d according to multi-objective optimization algorithm, namely can determine the position of damaged cable and the solution of nominal fatigue degree thereof more exactly with reasonable error from all ropes.

The Objective Programming in multi-objective optimization algorithm (Goal Attainment Method) can be adopted to solve current injury vector d, according to Objective Programming, formula (8) can transform the multi-objective optimization question shown in an accepted way of doing sth (26) and formula (27), in formula (26), γ is a real number, R is real number field, and area of space Ω limits the span of each element of vectorial d.Formula (26) be meant to the minimum real number γ of searching one, formula (27) is met.In formula (27), G (d) is defined by formula (28), the middle deviation allowed between G (d) and vectorial g of the product representation formula (27) of weighing vector W and γ in formula (27), the definition of g is see formula (13), and its value calculates in the 5th step.During actual computation, vector W can be identical with vectorial g.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of evaluation object current nominal fatigue vector d.

minimize γ

(26)

γ∈R,d∈Ω

G(d)-Wγ≤g (27)

G(d)＝abs(ΔC·d-C+C _o) (28)

The element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be one-to-one relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load intensity of variation of corresponding evaluation object; The coding rule of the element of vector d and vectorial d _othe coding rule of element identical.

12 step: definition evaluation object current actual damage vector d ^a, evaluation object current actual damage vector d ^aelement number equal the quantity of evaluation object, evaluation object current actual damage vector d ^aelement and evaluation object between be one-to-one relationship, evaluation object current actual damage vector d ^aelement numerical value represent actual damage degree or the real load intensity of variation of corresponding evaluation object; Vector d ^athe coding rule of element and vectorial d _othe coding rule of element identical.The evaluation object utilizing formula (15) to express current actual damage vector d ^aa kth element d ^a _kwith evaluation object initial damage vector d _oa kth element d _okwith a kth element d of evaluation object current nominal fatigue vector d _kbetween relation, calculate evaluation object current actual damage vector d ^aall elements.

D ^a _krepresent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system ^a _krepresent its current actual damage, d ^a _krepresent when being 0 that the support cable of its correspondence is without health problem, d ^a _knumerical value represents when not being 0 that the support cable of its correspondence is the support cable of unsoundness problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, the degree of the lax or damage of its numerical response.

D ^a _krepresent the current actual health status of a kth evaluation object, if this evaluation object is a load, its definition is shown in formula (15), so d ^a _krepresent that it is relative to setting up initial mechanical Calculation Basis model A _otime the structure variable quantity of corresponding load that bears.

By current for evaluation object actual damage vector d ^ain relevant to support cable Q element take out, composition support cable current actual damage vector d ^ca, support cable current actual damage vector d ^cathe coding rule of element and Initial cable force vector F _othe coding rule of element identical.Support cable current actual damage vector d ^cah element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., Q; Support cable current actual damage vector d ^camiddle numerical value be not 0 element correspond to the support cable of unsoundness problem, adopt lossless detection method to distinguish damaged cable and slack line from the support cable of these unsoundness problems.Damaged cable is at support cable current actual damage vector d ^cathe numerical value of the element of middle correspondence just represents its degree of injury, the numerical value of corresponding element represents when being 100% that this support cable thoroughly loses load-bearing capacity, represent time between 0 and 100% that this support cable loses the load-bearing capacity of corresponding proportion, so far just have identified damaged cable and degree of injury thereof.

Can in the hope of the relax level of these ropes (i.e. the long adjustment amount of rope) according to formula (24) or formula (25).So just achieve the lax identification of support cable.So far just damaged cable and slack line is all identified.

So far this method achieves the accurate identification of the health status of core evaluation object with a kind of effective, cheap method.May exact value be departed to the recognition result of the health status of secondary evaluation object more, only require the correct health status identifying core evaluation object in the method.

13 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.

14 step: under specified requirements, the computing machine automatic operation communication alert equipment in health monitoring systems is reported to the police to monitor staff, owner and (or) the personnel that specify.

15 step: get back to the 4th step, starts by the circulation of the 4th step to the 15 step.

Claims (1)

1. simplify the recognition methods of generalized displacement hybrid monitoring loading problem rope, it is characterized in that described method comprises:

Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated support cable and load to be " evaluation object ", if the quantity sum of the quantity of evaluated support cable and load is N, namely the quantity of " evaluation object " is N; This method title " core evaluation object " specially refers to the evaluated support cable in " evaluation object ", and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; Specify when determining hybrid monitoring by the support cable of monitored Suo Li, if total Q root support cable in cable system, i.e. total Q core evaluation object, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment support cable ₁individual rope force data describes, and the change of Cable Structure Suo Li is exactly the change of the Suo Li of all appointment support cables; Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M ₁it is an integer being not less than that 0 is not more than Q; Specify when determining hybrid monitoring by the measured point of monitored strain, the monitored strain data of Cable Structure can by K in Cable Structure ₂the L of individual specified point and each specified point ₂the strain of individual assigned direction describes, and the change of Cable Structure strain data is exactly K ₂the change of all tested strain of individual specified point; Each total M ₂individual strain measurement value or calculated value characterize Cable Structure strain, M ₂for K ₂and L ₂long-pending; M ₂be be not less than 0 integer; Specify when determining hybrid monitoring by the measured point of monitored angle, the monitored angle-data of Cable Structure is by K in Cable Structure ₃individual specified point, cross the L of each specified point ₃the H of individual appointment straight line, each appointment straight line ₃individual angle coordinate component describes, and the change of Cable Structure angle is exactly the change of all specified points, all appointment straight line, all angle coordinate component of specifying; Each total M ₃individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M ₃for K ₃, L ₃and H ₃long-pending; M ₃it is an integer being not less than 0; Specify when determining hybrid monitoring by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure ₄the L of individual specified point and each specified point ₄the volume coordinate of individual assigned direction describes, and the change of Cable Structure shape data is exactly K ₄the change of all coordinate components of individual specified point; Each total M ₄individual coordinates measurements or calculated value characterize Cable Structure shape, M ₄for K ₄and L ₄long-pending; M ₄it is an integer being not less than 0; The monitored amount of comprehensive above-mentioned hybrid monitoring, whole Cable Structure has M monitored amount, and M is M ₁, M ₂, M ₃and M ₄sum, definition parameter K, K is M ₁, K ₂, K ₃and K ₄sum, M must be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; For simplicity, in the method by the monitored amount of the M listed by this step referred to as " monitored amount "; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centre-point load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated 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 comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact 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; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the one-component of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;

B. survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses;

C. actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the measured data of initial Cable Structure, the measured data of initial Cable Structure comprises Cable Structure centre-point 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 coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Non-destructive Testing Data comprising support cable, the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C _o, monitored amount initial value vector C _othe coding rule of coding rule and M monitored amount identical, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d _o, vectorial d _orepresent with initial mechanical Calculation Basis model A _othe initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d _oelement number equal N, d _oelement and evaluation object be one-to-one relationship, vectorial d _othe coding rule of element identical with the coding rule of evaluation object, if d _oevaluation object corresponding to some elements be support cable, so a d in cable system _othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses load-bearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the load-bearing capacity of corresponding proportion, if d _oevaluation object corresponding to some elements be some load, get d in this method _othis element numerical value be 0, the initial value representing the change of this load is 0, if there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d _oin each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U _o, actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F _o, according to comprising Cable Structure design data, the data of completion data obtain the length of all support cables when free state and Suo Li are 0, in free state time cross-sectional area and in free state time the weight of unit length, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum and Initial cable force vector F _othe coding rule of element identical,

Physical and mechanical properties parameter, the initial Cable Structure bearing generalized coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, support cable initial health data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure _owith all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A of Cable Structure _o, based on A _othe Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A _oevaluation object health status with evaluation object initial damage vector d _orepresent; Corresponding to A _othe initial value monitored amount initial value vector C of all monitored amount _orepresent; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _owith monitored amount current initial value vector C ^t _o; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _owith monitored amount current initial value vector C ^t _otime, the current initial mechanical Calculation Basis model A of Cable Structure ^t _ojust equal the initial mechanical Calculation Basis model A of Cable Structure _o, monitored amount current initial value vector C ^t _ojust equal monitored amount initial value vector C _o; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ^t _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ^t _o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ^t _otime, U ^t _ojust equal U _o; A ^t _othe initial health of evaluation object and A _othe health status of evaluation object identical, also use evaluation object initial damage vector d _orepresent, A in cyclic process below ^t _othe initial health of evaluation object use evaluation object initial damage vector d all the time _orepresent; U _oand d _oa _oparameter, by A _othe initial value of all monitored amount that obtains of Mechanics Calculation result and C _othe initial value of all monitored amount represented is identical, therefore alternatively C _oby A _omechanics Calculation result composition; U ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition;

E. from entering the circulation being walked to m step by e here; In structure military service process, constantly actual measurement obtains Cable Structure bearing generalized coordinate current data, all Cable Structure bearing generalized coordinate current data composition current cable structure actual measurement bearing generalized coordinate vector U ^t, vectorial U ^tdefinition mode and vectorial U _odefinition mode identical;

F. according to current cable structure actual measurement bearing generalized coordinate vector U ^t, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3 ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _o;

F1. U is compared ^twith U ^t _oif, U ^tequal U ^t _o, then A ^t _o, U ^t _oand C ^t _oremain unchanged, otherwise need to follow these steps to A ^t _o, U ^t _oand C ^t _oupgrade;

F2. U is calculated ^twith U _odifference, U ^twith U _odifference be exactly the generalized displacement of support of Cable Structure bearing about initial position, with generalized displacement of support vector V represent generalized displacement of support, V equals U ^tdeduct U _o, be one-to-one relationship between the element in generalized displacement of support vector V and generalized displacement of support component, in generalized displacement of support vector V, the numerical value of an element corresponds to the generalized displacement of an assigned direction of an appointment bearing;

F3. to A _oin Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint just takes from the numerical value of corresponding element in generalized displacement of support vector V, to A _omiddle Cable Structure bearing obtains the current initial mechanical Calculation Basis model A upgraded after applying generalized displacement of support constraint ^t _o, upgrade A ^t _owhile, U ^t _oall elements numerical value also uses U ^tall elements numerical value correspondence replaces, and namely have updated U ^t _o, so just obtain and correctly correspond to A ^t _ou ^t _o; Upgrade C ^t _omethod be: when renewal A ^t _oafter, obtain A by Mechanics Calculation ^t _oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C ^t _o; A ^t _othe initial health of support cable use evaluation object initial damage vector d all the time _orepresent;

G. at current initial mechanical Calculation Basis model A ^t _obasis on carry out several times Mechanics Calculation according to step g 1 to g4, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D _u;

G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal ^t _o, current initial Cable Structure bearing generalized coordinate vector U ^t _owith monitored amount current initial value vector C ^t _oafterwards, Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D must then be upgraded _u;

G2. at the current initial mechanical Calculation Basis model A of Cable Structure ^t _obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object on the basis of original damage or load, increase unit damage or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just supposes that this support cable is at vectorial d _othe basis that this support cable represented has a damage increases unit damage again, if this evaluation object is a load, just supposes that this load is at vectorial d _othe basis that this load represented has a variable quantity increases load unit change again, use D _ukrecord unit damage or the load unit change of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or load unit change, D _ukevaluation object unit change vector D _uan element, evaluation object unit change vector D _uthe coding rule of element and vectorial d _othe coding rule of element identical; The evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector _oelement number rule identical;

G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C ^t _oobtain a vector, then each element of this vector is calculated the unit damage or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;

G4. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d _othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;

H. current cable structure actual measurement bearing generalized coordinate vector U is obtained in actual measurement ^twhile, actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C ^t _owith monitored amount initial value vector C _odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time; Obtain the synchronization of monitored amount current value vector C in actual measurement, actual measurement obtains the rope force data of all Q root support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F _othe coding rule of element identical; The synchronization of monitored amount current value vector C is obtained in actual measurement, Actual measurement obtains the volume coordinate of two supporting end points of all Q root support cables, the difference of the volume 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 form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector _othe coding rule of element identical;

I. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be one-to-one relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object; The coding rule of the element of vector d and vectorial d _othe coding rule of element identical;

J. according to monitored amount current value vector C with monitored amount current initial value vector C ^t _o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;

$C=C_o^t+\mathrm{\Δ}C\·d$ Formula 1

K. evaluation object current actual damage vector d is defined ^a, evaluation object current actual damage vector d ^aelement number equal the quantity of evaluation object, evaluation object current actual damage vector d ^aelement and evaluation object between be one-to-one relationship, evaluation object current actual damage vector d ^aelement numerical value represent actual damage degree or the real load variable quantity of corresponding evaluation object; Vector d ^athe coding rule of element and vectorial d _othe coding rule of element identical;

L. the evaluation object utilizing formula 2 to express current actual damage vector d ^aa kth element d ^a _kwith evaluation object initial damage vector d _oa kth element d _okwith a kth element d of evaluation object current nominal fatigue vector d _kbetween relation, calculate evaluation object current actual damage vector d ^aall elements;

K=1 in formula 2,2,3 ...., N, d ^a _krepresent the current actual health status of a kth evaluation object, d ^a _krepresent when being 0 that a kth evaluation object is without health problem, d ^a _knumerical value represents when not being 0 that a kth evaluation object is the evaluation object of unsoundness problem, if this evaluation object is support cable, so a d in cable system ^a _krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d ^a _kthe degree of the lax or damage of this support cable of numerical response; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, evaluation object current actual damage vector d ^ain correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent; If this evaluation object is load, so a d ^a _krepresent the actual change amount of this load; After identifying slack line, utilize evaluation object current actual damage vector d ^athese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize the current cable force vector F that obtains in f step and current support cable two to support end points horizontal range vector, utilize the initial drift vector of the support cable obtained in b step, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the vectorial F of Initial cable force _othe physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in b step to use, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanic equivalent condition is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting above-mentioned two mechanic equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the non-destructive tests of support cable; Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so according to evaluation object current actual damage vector d ^athe health status of core evaluation object can be determined;

M. get back to e step, start the circulation next time being walked to m step by e.