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

Abstract

The invention provides a method for progressively recognizing the damaged cable of a load based on the streamlined strain monitoring process of the generalized displacement. According to the method, based on the strain monitoring process, whether a mechanical calculation benchmark model of a cable structure needs to be updated or not is determined through monitoring the generalized displacement of a bearer, the variation degree of a load and the damaged degree of a damaged cable. After that, the mechanical calculation benchmark model of the cable structure, with the generalized displacement of the bearer, the variation degree of the load and the damaged degree of the damaged cable being taken into account, is obtained. On the basis of the above model, according to the approximately linear relationships between the current numeric vector of a monitored quantity and the current initial numeric vector of the monitored quantity, the numerical value variation matrix of the monitored quantity with unit damage and the current nominal damage vector of an unknown current nominal damage vector, the health status of a core to-be-evaluated object can be recognized.

Description

Simplify generalized displacement strain monitoring load damaged cable progressive-type recognition method

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.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, this method identifies damaged cable based on strain monitoring (monitored strain is called " monitored amount " by this method), belongs 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.Respectively: one, determine evaluation object and monitored amount (also can be called monitored target); Two, the structural health conditions appraisal procedure of the method for knowledge base needed for cable structure health monitoring system and parameter, knowledge based storehouse (containing parameter) and the monitored amount of actual measurement is set up; Three, the software and hardware part of health monitoring systems.

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: determine evaluation object and 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, 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.

If total M in cable system ₁root support cable, i.e. total M ₁individual core evaluation object.

" the whole monitored strain data of structure " can by the specified point of K in structure and the strain of L assigned direction of each specified point describe, the change of structural strain data is exactly the change of all strains of K specified point.Each total individual strain measurement value of M (M=K × L) or calculated value carry out characterisation of structures strain information.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M monitored amount, and 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.

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

The first step: set up initial mechanical Calculation Basis model A _o, the variable quantity of " load likely changed " is setting up initial mechanical Calculation Basis model A _otime be all 0, the variable quantity of " load likely changed " that is identifies below is relative to setting up initial mechanical Calculation Basis model A _otime the structure variable quantity of corresponding load that bears.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, in formula, 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, 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.Vector C _obe to be arranged according to a definite sequence by the monitored amount of M to form, putting in order to this there is no particular/special requirement, only requires all associated vector also array data in this order below.

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.

Second step: circulation starts.When circulation starts each time, first need the evaluation object current initial damage vector d set up or set up when this circulation starts ⁱ _o(i=1,2,3 ...), set up the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _o(such as finite element benchmark model, A in circulation each time ⁱ _oconstantly update), A ⁱ _obearing generalized coordinate " current initial Cable Structure bearing generalized coordinate vector U ⁱ _o" express.Letter i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index in the method, i.e. i-th circulation.

The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d ⁱ _o(such as formula (3) Suo Shi), uses d ⁱ _owhen representing that this circulation starts, Cable Structure is (with current initial mechanical Calculation Basis model A ⁱ _orepresent) the health status of evaluation object.

$d_o^i={\left[\begin{array}{cccccccccc}{d}_{o1}^i& d_{o2}^i& \·& \·& \·& d_{ok}^i& \·& \·& \·& d_{oN}^i\end{array}\right]}^T---\left(3\right)$

D in formula (3) ⁱ _ok(i=1,2,3, K=1,2,3 ...., N) represent i-th time circulation start time, current initial mechanical Calculation Basis model A ⁱ _oin the original state of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d ⁱ _okrepresent its initial damage, d ⁱ _okrepresent 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 one " load that may change ", so d ⁱ _okrepresent that it is relative to setting up initial mechanical Calculation Basis model A _otime the structure variable quantity of the respective components of corresponding load that bears.

Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ⁱ _o, the current initial mechanical Calculation Basis model A of Cable Structure is namely set up for the first time at initial time ⁱ _otime, U ⁱ _ojust equal U _o.

Set up and upgrade d ⁱ _omethod as follows:

When first time, circulation started, (foundation formula (3) is designated as d to set up the current initial damage vector of evaluation object ¹ _o) time, d ¹ _ojust equal d _o.I-th (i=2,3,4,5,6 ...) the secondary evaluation object current initial damage vector d needed when starting that circulates ⁱ _o, be front once (namely the i-th-1 time, i=2,3,4,5,6 ...) circulation terminate before calculate obtain, concrete grammar is described below.

I-th (i=1,2,3,4,5,6 ...) secondary circulation needs the Mechanics Calculation benchmark model set up or the Mechanics Calculation benchmark model of Cable Structure set up to be designated as current initial mechanical Calculation Basis model A when starting ⁱ _o.Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ⁱ _o.Vector U ⁱ _odefinition mode and vectorial U _odefinition mode identical, must set up or set up current initial Cable Structure bearing generalized coordinate vector U when circulation starts each time ⁱ _o.

Set up, upgrade A ⁱ _oand U ⁱ _omethod as follows:

The Mechanics Calculation benchmark model of the Cable Structure set up when first time, circulation started is designated as A ¹ _o, A ¹ _oequal A _o, U ¹ _oequal U _o.A in circulation each time ⁱ _oand U ⁱ _obe constantly update, concrete grammar is described below; At the end of circulation each time, upgrade A ⁱ _oand U ⁱ _othe Mechanics Calculation benchmark model of the Cable Structure required when starting that next time circulated, concrete grammar is described below.

This method " monitored amount current initial value vector C ⁱ _o" (i=1,2,3 ...) initial value (see formula (4)) of all monitored amounts of specifying when representing that i-th time (i=1,2,3,4,5,6 ...) circulation starts, C ⁱ _oalso can be called " the monitored amount that circulates for i-th time current initial value vector ".

$C_o^i={\left[\begin{array}{cccccccccc}{C}_{o1}^i& C_{o2}^i& \·& \·& \·& C_{oj}^i& \·& \·& \·& C_{oM}^i\end{array}\right]}^T---\left(4\right)$

C in formula (2) ⁱ _oj(i=1,2,3, J=1,2,3 ...., M) be jth monitored amount when circulating beginning for i-th time, in Cable Structure.Vector C ⁱ _obe to be arranged according to a definite sequence by the monitored amount of previously defined M to form, putting in order to this there is no particular/special requirement, only requires all associated vector also array data in this order below.

At Modling model A ⁱ _owhile set up " monitored amount current initial value vector C ⁱ _o", monitored amount current initial value vector C ⁱ _orepresent and correspond to A ⁱ _othe concrete numerical value of all monitored amount, C ⁱ _oelement and C _oelement one_to_one corresponding, represent that all monitored amounts are in A in Cable Structure respectively ⁱ _oand A _oconcrete numerical value during two states.

Set up and upgrade C ⁱ _oconcrete grammar as follows:

When first time, circulation started, C ¹ _o(i=1, C ⁱ _obe embodied as C ¹ _o) equal C _o; I-th (i=2,3,4,5,6 ...) secondary i-th circulation " the monitored amount current initial value vector C needed when starting that circulates ⁱ _o", be front once (namely the i-th-1 time, i=2,3,4,5,6 ...) circulation terminate before calculate obtain, concrete grammar is described below.I-th (i=1,2,3,4,5,6 ...) in circulation, " monitored amount current initial value vector C ⁱ _o" be constantly update, concrete grammar is described below.Due to according to model A ⁱ _othe initial value calculating gained monitored amount, reliably close to corresponding measured value, in describing below, will represent this calculated value composition of vector and measured value composition of vector with prosign.

U ⁱ _oand d ⁱ _oa ⁱ _ocharacterisitic parameter, C ⁱ _oa ⁱ _omechanics Calculation result composition.

3rd step: in Cable Structure military service process, in circulation each time, in other words in the i-th (i=1,2,3,4,5,6 ...) in secondary circulation, at known A ⁱ _o, U ⁱ _o, C ⁱ _oand d ⁱ _oafter, continuous Actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C ⁱ.C ⁱelement and C _oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.Obtaining vector T ⁱwhile, 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 ⁱ.

In acquisition vector T ⁱafter, upgrade A ⁱ _o, U ⁱ _o, C ⁱ _oand d ⁱ _o:

4th step: circulation time must first be set up " unit damage monitored numerical quantity transformation matrices " and " evaluation object unit change vector " each time, and " unit damage monitored numerical quantity transformation matrices " that i-th circulation is set up is designated as Δ C ⁱ(i=1,2,3 ...)." evaluation object unit change vector " that i-th circulation is set up is designated as D ⁱ _u.Δ C in circulation each time ⁱand D ⁱ _uneed according to circumstances to constantly update, namely at the current initial mechanical Calculation Basis model A of renewal ⁱ _o, current initial Cable Structure bearing generalized coordinate vector U ⁱ _owith monitored amount current initial value vector C ⁱ _oafter, upgrade unit damage monitored numerical quantity transformation matrices Δ C ⁱwith evaluation object unit change vector D ⁱ _u.

First unit damage monitored numerical quantity transformation matrices Δ C is set up in the steps below when circulation starts each time ⁱwith evaluation object unit change vector D ⁱ _u; If have updated A in the third step ⁱ _o, (namely upgrading) unit damage monitored numerical quantity transformation matrices Δ C so must be re-established in this step ⁱwith evaluation object unit change vector D ⁱ _u; If do not upgrade A in the third step ⁱ _o, unit damage monitored numerical quantity transformation matrices Δ C so need not be re-established in this step ⁱwith evaluation object unit change vector D ⁱ _u; Set up and re-establish (namely upgrading) Δ C ⁱand D ⁱ _udetailed process identical, arrange as follows:

At the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _obasis on carry out several times calculating, calculation times numerically equals the quantity of all ropes.Calculate hypothesis each time to only have an evaluation object (original load change amount can be 0 in original damage or load change amount, also can not be 0) 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 increases unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again, 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).For convenience of calculating, can be all structural health conditions when this circulation is started when in circulation, setting increases unit damage or load unit change each time as being completely healthy, and set on this basis unit damage or load unit change (in subsequent step, calculate, the damage numerical value of evaluation object or load change amount---be called nominal fatigue d ⁱ _c(i=1,2,3 ...), all relative to when this circulation is started, by the health status of evaluation object as being completely healthy speech, therefore must the formula that hereinafter provides of foundation the nominal fatigue calculated is converted into true damage).With occurring in the calculating each time once circulated 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, and the unit damage value or the load unit that suppose there is the evaluation object of unit damage or load unit change each time change unit damage value or the load unit change numerical value that numerical value can be different from other evaluation objects, with " evaluation object unit change vector D ⁱ _u" (such as formula (5) Suo Shi) record the unit damage of supposition or the load unit change of all evaluation objects in each circulation, first time circulation time be designated as D ¹ _ucalculate the current calculated value all utilizing M that mechanics method (such as finite element method) calculates Cable Structure, that specified above monitored amount each time, the current calculated value calculating a gained M monitored amount each time forms one " monitored amount calculation current vector ", and (when supposing that a kth evaluation object has unit damage, available formula (6) represents the monitored amount calculation current vector C of all M specified a monitored amount ⁱ _tk); The monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C ⁱ _o, gained vector is exactly that " the numerical value change vector of monitored amount " of (to have the position of the rope of unit damage or numbering etc. for mark) (when a kth evaluation object has unit damage, uses δ C under this condition ⁱ _krepresent the numerical value change vector of monitored amount, δ C ⁱ _kdefinition see formula (7), formula (8) and formula (9), formula (7) deducts after formula (4) again divided by vectorial D for formula (6) ⁱ _ua kth element D ⁱ _ukgained), the numerical value change vector δ C of monitored amount ⁱ _keach element representation change (such as D due to the unit damage of that evaluation object (such as a kth evaluation object) or load unit that suppose there is unit damage or load unit change when calculating ⁱ _uk), and the numerical value knots modification of monitored amount corresponding to this element caused is relative to the unit damage of supposition or load unit change numerical value D ⁱ _ukrate of change; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ", the numerical value change vector of each monitored amount has M element, forms by this N number of " numerical value change vector of monitored amount " " the unit damage monitored numerical quantity transformation matrices Δ C having M × N number of element successively ⁱ" (the capable N row of M), each vectorial δ C ⁱ _k(k=1,2,3 ...., N) be matrix Δ C ⁱrow, Δ C ⁱdefinition such as formula shown in (10).

$D_u^i={\left[\begin{array}{cccccccccc}{D}_{u1}^i& D_{u2}^i& \·& \·& \·& D_{uk}^i& \·& \·& \·& D_{uN}^i\end{array}\right]}^T---\left(5\right)$

Evaluation object unit change vector D in formula (5) ⁱ _uelement D ⁱ _uk(i=1,2,3, K=1,2,3 ...., N) represent that the unit damage of a kth evaluation object of supposition in i-th circulation or load unit change numerical value, vectorial D ⁱ _uin the numerical value of each element can be the same or different.

$C_{tk}^i={\left[\begin{array}{cccccccccc}{C}_{tk1}^i& C_{tk2}^i& \·& \·& \·& C_{tkj}^i& \·& \·& \·& C_{tkM}^i\end{array}\right]}^T---\left(6\right)$

Elements C in formula (6) ⁱ _tkj(i=1,2,3, K=1,2,3 ...., N; J=1,2,3 ...., M) represent i-th circulation due to a kth evaluation object have unit damage or a load unit change time, according to the calculating current value of the individual monitored amount of specifying of the jth corresponding to coding rule.

${\mathrm{\δC}}_k^i=\frac{C_{tk}^i-C_o^i}{D_{uk}^i}---\left(7\right)$

In formula (7) each amount subscript i (i=1,2,3 ...) represent i-th circulation, subscript k (k=1,2,3 ...., N) represent the unit damage that a kth evaluation object increases or load unit change, D in formula ⁱ _ukvectorial D ⁱ _uin a kth element.Vector δ C ⁱ _kdefinition such as formula shown in (7) and formula (8), δ C ⁱ _kjth (j=1,2,3 ...., M) individual element δ C ⁱ _kj(defining such as formula (9) Suo Shi) represents in i-th circulation, sets up matrix Δ C ⁱtime, assuming that a kth evaluation object has the knots modification calculating a gained jth monitored amount when unit damage or load unit change relative to the unit damage supposed or load unit change D ⁱ _ukrate of change.

${\mathrm{\δC}}_k^i={\left[\begin{array}{cccccccccc}{\mathrm{\δC}}_{k1}^i& {\mathrm{\δC}}_{k2}^i& \·& \·& \·& {\mathrm{\δC}}_{kj}^i& \·& \·& \·& {\mathrm{\δC}}_{kM}^i\end{array}\right]}^T---\left(8\right)$

${\mathrm{\δC}}_{kj}^i=\frac{C_{tkj}^i-C_{oj}^i}{D_{uk}^i}---\left(9\right)$

${\mathrm{\ΔC}}^i=\[\begin{array}{cccccccccc}{\mathrm{\δC}}_1^i& {\mathrm{\δC}}_2^i& \·& \·& \·& {\mathrm{\δC}}_k^i& \·& \·& \·& {\mathrm{\δC}}_N^i\end{array}\]---\left(10\right)$

Vectorial δ C in formula (10) ⁱ _k(i=1,2,3 ...., k=1,2,3 ...., N) represent in i-th circulation, because a kth evaluation object increases unit damage or load unit change D ⁱ _ukcause, the change of the relative value of all monitored amounts.Matrix Δ C ⁱthe coding rule of row (subscript k) and vectorial d above ⁱ _othe coding rule of subscript k of element identical.

5th step: the current health state identifying Cable Structure.Detailed process is as follows.

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

$C^i=C_o^i+{\mathrm{\ΔC}}^i\·d_c^i---\left(11\right)$

$C^i-C_o^i={\mathrm{\ΔC}}^i\·d_c^i---\left(12\right)$

Monitored amount current value vector C in formula (11) and formula (12) ⁱdefinition be similar to monitored amount current initial value vector C ⁱ _odefinition, see formula (13); Evaluation object current nominal fatigue vector d ⁱ _cdefinition see formula (14).

$C^i={\left[\begin{array}{cccccccccc}{C}_1^i& C_2^i& \·& \·& \·& C_j^i& \·& \·& \·& C_M^i\end{array}\right]}^T---\left(13\right)$

Elements C in formula (13) ⁱ _j(i=1,2,3 ....; J=1,2,3 ...., M) be i-th circulation time Cable Structure, the current value that be numbered the monitored amount of j of foundation corresponding to coding rule.

$d_c^i={\left[\begin{array}{cccccccccc}{d}_{c1}^i& d_{c2}^i& \·& \·& \·& d_{ck}^i& \·& \·& \·& d_{cN}^i\end{array}\right]}^T---\left(14\right)$

D in formula (14) ⁱ _ck(i=1,2,3 ....; K=1,2,3 ...., N) be the current nominal fatigue of a kth evaluation object in i-th circulation or current nominal load changing value, vectorial d ⁱ _cthe coding rule of subscript k of element and matrix Δ C ⁱthe coding rule of row identical.

When support cable actual damage is not too large, because Cable Structure material is still in the linear elasticity stage, the distortion of Cable Structure is also less, formula (11) or a kind of like this linear relationship represented by formula (12) less with the error of actual conditions, error can use error vector e ⁱ(formula (15)) define, the error of expression (11) or the shown linear relationship of formula (12).

$e^i=abs({\mathrm{\ΔC}}^i\·d_c^i-C^i+C_o^i)---\left(15\right)$

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

There is certain error in the linear relationship represented by formula (11) or formula (12), therefore obtains evaluation object current nominal fatigue vector d ⁱ _cacceptable solution (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more exactly from cable system, also can determine load change numerical value more exactly) become a rational solution, available formula (16) expresses this method.

$abs({\mathrm{\ΔC}}^i\·d_c^i-C^i+C_o^i)\≤g^i---\left(16\right)$

In formula (16), abs () is the function that takes absolute value, vectorial g ⁱdescription departs from the legitimate skew of ideal linearity relation (formula (11) or formula (12)), is defined by formula (17).

$g^i={\left[\begin{array}{cccccccccc}{g}_1^i& g_2^i& \·& \·& \·& g_j^i& \·& \·& \·& g_M^i\end{array}\right]}^T---\left(17\right)$

In formula (17) describe the maximum allowable offset departing from formula (11) or the ideal linearity relation shown in formula (12) in i-th circulation.Vector g ⁱthe error vector e that can define according to formula (15) ⁱtentative calculation is selected.

At monitored amount current initial value vector C ⁱ _o, unit damage monitored numerical quantity transformation matrices Δ C ⁱwith monitored amount current value vector C ⁱtime known, suitable algorithm (such as multi-objective optimization algorithm) can be utilized to solve formula (16), obtain evaluation object current nominal fatigue vector d ⁱ _cacceptable solution, evaluation object current actual damage vector d ⁱthe element of (formula (18) is shown in definition) can calculate according to formula (19), thus can by d ⁱdetermine the health status of evaluation object.

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

D in formula (18) ⁱ _k(i=1,2,3, K=1,2,3 ...., N) represent the current actual health status of a kth evaluation object in i-th circulation, formula (19) is shown in its definition, if this evaluation object is a support cable (or pull bar) in cable system, so d ⁱ _krepresent its current actual damage, d ⁱ _krepresent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; If this evaluation object is load, so a d ⁱ _krepresent the current actual change numerical value of the load of its correspondence, vectorial d ⁱthe coding rule of element and formula (1) in vectorial d _othe coding rule of element identical.So far this method achieves two kinds of functions that existing method can not possess, be respectively, one, when the load change that structural bearings generation generalized displacement and structure are born, generalized displacement of support and load change can be rejected on the impact of Cable Structure health status recognition result, thus identify the structure health monitoring method of damaged cable exactly; Two, this method is while identifying damaged cable, can also identify the change of load simultaneously, and namely this method can reject generalized displacement of support and the impact with the change of support cable health status, realizes the correct identification of load change degree.

D in formula (19) ⁱ _ok(i=1,2,3,4, K=1,2,3 ...., N) be evaluation object current initial damage vector d ⁱ _oa kth element, d ⁱ _ckevaluation object current nominal fatigue vector d ⁱ _ca kth element.

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.

6th step: judge whether to terminate this (i-th time) circulation, if so, then completes this tailing in work before terminating that circulates, for next time (namely the i-th+1 time, i=1,2,3,4 ...) circulating prepares Mechanics Calculation benchmark model and necessary vector.Detailed process is as follows:

Current nominal fatigue vector d is tried to achieve in this (i-th time) circulation ⁱ _cafter, first, set up mark vector B according to formula (20) ⁱ, formula (21) gives mark vector B ⁱthe definition of a kth element; If mark vector B ⁱelement be 0 entirely, then get back to the 3rd step and proceed health monitoring to Cable Structure and calculating; If mark vector B ⁱelement be not 0, then after completing subsequent step entirely, enter and circulate next time.

So-called subsequent step is: first, according to formula (22) calculate next time (namely the i-th+1 time, i=1,2,3,4 ...) needed for circulation initial damage vector d ⁱ⁺¹ _oeach element d ⁱ⁺¹ _ok; The second, at Mechanics Calculation benchmark model A _obasis on, make A _oin the health status of evaluation object be d ⁱ⁺¹ _oinstead of be d _oafter, more further to A _oin Cable Structure apply generalized displacement of support constraint (as previously mentioned, the numerical value of the generalized displacement of support of applying constraint just take from generalized displacement of support vector V, and generalized displacement of support vector V equals U ⁱdeduct U _o), so just obtain next time (namely the i-th+1 time, i=1,2,3,4 ...) current initial mechanical Calculation Basis mould A needed for circulation ⁱ⁺¹ _o, next time (namely the i-th+1 time, i=1,2,3,4 ...) current initial Cable Structure bearing generalized coordinate vector U needed for circulation ⁱ⁺¹ _oequal U ⁱ _o, to A ⁱ⁺¹ _ocarry out Mechanics Calculation to obtain corresponding to A ⁱ⁺¹ _oall monitored amount, current concrete numerical value, these concrete numerical value compositions next time (namely the i-th+1 time, i=1,2,3,4 ...) the current initial value vector C of monitored amount needed for circulation ⁱ⁺¹ _o.

$B^i={\left[\begin{array}{cccccccccc}{B}_1^i& B_2^i& \·& \·& \·& B_k^i& \·& \·& \·& B_N^i\end{array}\right]}^T---\left(20\right)$

Mark vector B in formula (20) ⁱsubscript i represent i-th circulation, its element B ⁱ _k(k=1,2,3 ..., N) subscript k represent the health status feature of a kth evaluation object, can only get 0 and 1 two amount, concrete value rule is shown in formula (21).

Element B in formula (21) ⁱ _kmark vector B ⁱa kth element, D ⁱ _ukevaluation object unit change vector D ⁱ _ua kth element (see formula (5)), d ⁱ _ckevaluation object current nominal fatigue vector d ⁱ _ca kth element (see formula (14)), they all represent the relevant information of a kth evaluation object.

D in formula (22) ⁱ _ukevaluation object unit change vector D ⁱ _ua kth element (see formula (5)), d ⁱ _okevaluation object current initial damage vector d ⁱ _oa kth element (see formula (3)).

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

Hardware components comprises monitoring system (comprising monitored amount monitoring system, Cable Structure bearing generalized coordinate monitoring system), signal picker and computing machine etc.Require that Real-Time Monitoring obtains the measured data of required Cable Structure bearing generalized coordinate, require each monitored amount of Real-Time Monitoring simultaneously.

Software section should complete the process set by this method, namely 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; If total M in cable system ₁root support cable, the quantity of obvious core evaluation object is exactly M ₁; Determine the point being monitored of specifying, namely point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified points numbering; Determine monitored should the changing direction of point being monitored, and give all monitored strain numberings of specifying, " monitored strain numbering " will be used for generating vector sum matrix in subsequent step, and " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; The external force that object, structure are born can be described as load, and load 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, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular 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,

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 _ocable Structure bearing generalized coordinate data be exactly initial Cable Structure bearing generalized coordinate vector U _o; 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; 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, A in the method _o, C _o, d _oand U _oconstant;

E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting ⁱ _o; When i-th circulation starts, corresponding to A ⁱ _o" Cable Structure bearing generalized coordinate data " with current initial Cable Structure bearing generalized coordinate vector U ⁱ _orepresent, vectorial U ⁱ _odefinition mode and vectorial U _odefinition mode identical, U ⁱ _oelement and U _oelement one_to_one corresponding; The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d ⁱ _o, d ⁱ _ocable Structure A when representing that this circulation starts ⁱ _othe health status of evaluation object, d ⁱ _odefinition mode and d _odefinition mode identical, d ⁱ _oelement and d _oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C ⁱ _orepresent, vectorial C ⁱ _odefinition mode and vectorial C _odefinition mode identical, C ⁱ _oelement and C _oelement one_to_one corresponding, monitored amount current initial value vector C ⁱ _orepresent and correspond to A ⁱ _othe concrete numerical value of all monitored amount; U ⁱ _oand d ⁱ _oa ⁱ _ocharacterisitic parameter, C ⁱ _oby A ⁱ _omechanics Calculation result composition; When first time, circulation started, A ⁱ _obe designated as A ¹ _o, set up A ¹ _omethod for making A ¹ _oequal A _o; When first time, circulation started, U ⁱ _obe designated as U ¹ _o, set up U ¹ _omethod for making U ¹ _oequal U _o; When first time, circulation started, d ⁱ _obe designated as d ¹ _o, set up d ¹ _omethod for making d ¹ _oequal d _o; When first time, circulation started, C ⁱ _obe designated as C ¹ _o, set up C ¹ _omethod for making C ¹ _oequal C _o;

F. from entering the circulation being walked to q step by f 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 ⁱ, vectorial U ⁱdefinition mode and vectorial U _odefinition mode identical, U ⁱelement and U _oelement one_to_one corresponding; Vectorial U is obtained in actual measurement ⁱwhile, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C ⁱ, vectorial C ⁱdefinition mode and vectorial C _odefinition mode identical, C ⁱelement and C _oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time;

G. according to current cable structure actual measurement bearing generalized coordinate vector U ⁱ, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g3 ⁱ _o, monitored amount current initial value vector C ⁱ _owith current initial Cable Structure bearing generalized coordinate vector U ⁱ _o, and evaluation object current initial damage vector d ⁱ _oremain unchanged;

G1. U is compared ⁱwith U ⁱ _oif, U ⁱequal U ⁱ _o, then A ⁱ _o, C ⁱ _oand U ⁱ _oremain unchanged, otherwise need to follow these steps to A ⁱ _o, C ⁱ _oand U ⁱ _oupgrade;

G2. U is calculated ⁱwith U _odifference, U ⁱwith 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 ⁱdeduct U _o;

G3. 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 _oin Cable Structure bearing apply generalized displacement of support constraint after obtain upgrade current initial mechanical Calculation Basis model A ⁱ _o, upgrade A ⁱ _owhile, U ⁱ _oall elements numerical value also uses U ⁱall elements numerical value correspondence replaces, and namely have updated U ⁱ _o, so just obtain and correctly correspond to A ⁱ _ou ⁱ _o, now d ⁱ _oremain unchanged; As renewal A ⁱ _oafter, A ⁱ _othe health status evaluation object of rope current initial damage vector d ⁱ _orepresent, A ⁱ _othe current initial Cable Structure bearing generalized coordinate vector U of bearing generalized coordinate ⁱ _orepresent; Upgrade C ⁱ _omethod be: when renewal A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C ⁱ _o;

H. at current initial mechanical Calculation Basis model A ⁱ _obasis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C ⁱwith evaluation object unit change vector D ⁱ _u;

H1., when i-th circulation starts, directly Δ C is obtained by method listed by step h2 to step h4 ⁱand D ⁱ _u; In other moment, when in step g to A ⁱ _oafter upgrading, Δ C must be regained by method listed by step h2 to step h4 ⁱand D ⁱ _uif, not to A in step g ⁱ _oupgrade, then directly proceed to step I herein and carry out follow-up work;

H2. at current initial mechanical Calculation Basis model A ⁱ _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 suppose that this support cable increases unit damage again, if this evaluation object is a load, just suppose that this load 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 again different from during other time calculates the evaluation object increasing unit damage or load unit change again, 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; When supposing that a kth evaluation object increases unit damage or load unit change again, use C ⁱ _tkrepresent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target; C ⁱ _tkdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tkelement and C _oelement one_to_one corresponding;

H3. the vectorial C calculated each time ⁱ _tkdeduct vectorial C ⁱ _oobtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated divided by this unit damage or load unit change numerical value supposed ⁱ _k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ";

H4. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively ⁱ"; Unit damage monitored numerical quantity transformation matrices Δ C ⁱeach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ⁱevery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of row and vectorial d _othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of coding rule and M monitored amount of row identical;

I. current nominal fatigue vector d is defined ⁱ _cwith current actual damage vector d ⁱ, d ⁱ _cand d ⁱelement number equal the quantity of evaluation object, d ⁱ _cand d ⁱelement and evaluation object between be one-to-one relationship, d ⁱ _celement numerical value represent nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object, d ⁱ _cand d ⁱwith evaluation object initial damage vector d _oelement number rule identical, d ⁱ _celement, d ⁱelement and d _oelement be one-to-one relationship;

J. according to monitored amount current value vector C ⁱwith " monitored amount current initial value vector C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " current nominal fatigue vector d ⁱ _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1 ⁱ _cother outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d ⁱ _c;

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

K. the current actual damage vector d utilizing formula 2 to express ⁱa kth element d ⁱ _kwith evaluation object current initial damage vector d ⁱ _oa kth element d ⁱ _okwith current nominal fatigue vector d ⁱ _ca kth element d ⁱ _ckbetween relation, calculate current actual damage vector d ⁱall elements;

K=1 in formula 2,2,3 ..., N; d ⁱ _krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system ⁱ _krepresent its current actual damage, d ⁱ _krepresent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; So according to evaluation object current actual damage vector d ^athe health status of core evaluation object can be determined;

L. current nominal fatigue vector d is tried to achieve ⁱ _cafter, set up mark vector B according to formula 3 ⁱ, formula 4 gives mark vector B ⁱthe definition of a kth element;

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

Element B in formula 4 ⁱ _kmark vector B ⁱa kth element, D ⁱ _ukevaluation object unit change vector D ⁱ _ua kth element, d ⁱ _ckevaluation object current nominal fatigue vector d ⁱ _ca kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;

If m. mark vector B ⁱelement be 0 entirely, then get back to step f continue this circulation; If mark vector B ⁱelement be not 0 entirely, then enter next step, i.e. step n;

N. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the i-th+1 time circulation ⁱ⁺¹ _oeach element;

D in formula 5 ⁱ⁺¹ _okthe evaluation object current initial damage vector d next time, namely needed for the i-th+1 time circulation ⁱ⁺¹ _oa kth element, d ⁱ _okthis, i.e. the evaluation object of i-th circulation current initial damage vector d ⁱ _oa kth element, D ⁱ _ukthe evaluation object unit change vector D of i-th circulation ⁱ _ua kth element, B ⁱ _kthe mark vector B of i-th circulation ⁱa kth element, k=1 in formula 5,2,3 ..., N;

O. at initial mechanical Calculation Basis model A _obasis on, first 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, then makes the health status of rope be d ⁱ⁺¹ _oafter obtain be exactly next time, namely the i-th+1 time circulation needed for Mechanics Calculation benchmark model A ⁱ⁺¹; Obtain A ⁱ⁺¹after, obtain A by Mechanics Calculation ⁱ⁺¹in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C ⁱ⁺¹ _o;

P. take off once, current initial Cable Structure bearing generalized coordinate vector U namely needed for the i-th+1 time circulation ⁱ⁺¹ _oequal the current initial Cable Structure bearing generalized coordinate vector U of i-th circulation ⁱ _o;

Q. get back to step f, start to circulate next time.

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.

If total M in cable system ₁root support cable, i.e. total M ₁individual core evaluation object.

" the whole monitored strain data of structure " can by the specified point of K in structure and the strain of L assigned direction of each specified point describe, the change of structural strain data is exactly the change of all strains of K specified point.Each total individual strain measurement value of M (M=K × L) or calculated value carry out characterisation of structures strain information.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M monitored amount, and the quantity that M must not be less than core evaluation object adds the quantity N 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, use the direct survey calculation of conventional method to obtain the initial value of all monitored amount of Cable Structure, form monitored amount initial value vector C _o; Use conventional method (consult reference materials or survey) 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 simultaneously.

Monitored amount initial value vector C is obtained at Actual measurement _owhile, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Non-destructive Testing Data etc. that the Actual measurement data of Cable Structure comprise support cable can express the data of the health status of rope, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinate data (comprise bearing about Descartes rectangular coordinate system X, Y, the volume coordinate of Z axis and angular coordinate and initial Cable Structure bearing spatial data and initial Cable Structure bearing angular 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, the measured datas such as structure space measurement of coordinates data.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, use 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.

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 support cable _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.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.Corresponding to A _ocable Structure bearing generalized coordinate data initial Cable Structure bearing generalized coordinate vector U _orepresent; U _oand d _oa _oparameter, C _oby A _omechanics Calculation result composition.

3rd step: in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting ⁱ _o; That i-th circulation needs when starting, corresponding to the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _ocable Structure bearing generalized coordinate data composition current initial Cable Structure bearing generalized coordinate vector U ⁱ _o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ⁱ _otime, U ⁱ _ojust equal U _o.The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d ⁱ _o, d ⁱ _ocable Structure A when representing that this circulation starts ⁱ _othe health status of evaluation object, d ⁱ _odefinition mode and d _odefinition mode identical, d ⁱ _oelement and d _oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C ⁱ _orepresent, vectorial C ⁱ _odefinition mode and vectorial C _odefinition mode identical, C ⁱ _oelement and C _oelement one_to_one corresponding, monitored amount current initial value vector C ⁱ _orepresent and correspond to A ⁱ _othe concrete numerical value of all monitored amount; U ⁱ _oand d ⁱ _oa ⁱ _ocharacterisitic parameter; C ⁱ _oby A ⁱ _omechanics Calculation result composition; When first time, circulation started, A ⁱ _obe designated as A ¹ _o, set up A ¹ _omethod for making A ¹ _oequal A _o; When first time, circulation started, U ⁱ _obe designated as U ¹ _o, set up U ¹ _omethod for making U ¹ _oequal U _o; When first time, circulation started, d ⁱ _obe designated as d ¹ _o, set up d ¹ _omethod for making d ¹ _oequal d _o; When first time, circulation started, C ⁱ _obe designated as C ¹ _o, set up C ¹ _omethod for making C ¹ _oequal C _o.

4th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (such as containing strain measurement system, signal conditioner etc.), Cable Structure bearing angular coordinate monitoring system (containing total powerstation, angle measuring sensor, signal conditioner etc.), signal (data) collector, computing machine and communication alert equipment.The bearing angular coordinate of each monitored amount, 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.

5th step: establishment the system software of installation and operation this method 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.

6th step: step starts circulation running thus, and in structure military service process, continuous Actual measurement obtains the currency of all monitored amounts in Cable Structure, all these numerical value form monitored amount current value vector C ⁱ, vectorial C ⁱdefinition mode and vectorial C _odefinition mode identical, C ⁱelement and C _oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.

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

7th step: obtaining current cable structure actual measurement bearing generalized coordinate vector U ⁱafter, compare U ⁱand U ⁱ _oif, U ⁱequal U ⁱ _o, then do not need A ^t _o, U ^t _oand C ^t _oupgrade, otherwise need current initial mechanical Calculation Basis model A ⁱ _o, current initial Cable Structure bearing generalized coordinate vector U ⁱ _owith monitored amount current initial value vector C ⁱ _oupgrade, and evaluation object current initial damage vector d ⁱ _oremain unchanged, update method follows these steps to a to step c and carries out:

A. U is calculated ⁱwith U _odifference, U ⁱwith 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 ⁱdeduct 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.

B. 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 _oin Cable Structure bearing apply generalized displacement of support constraint after obtain upgrade current initial mechanical Calculation Basis model A ⁱ _o, upgrade A ⁱ _owhile, U ⁱ _oall elements numerical value also uses U ⁱall elements numerical value correspondence replaces, and namely have updated U ⁱ _o, so just obtain and correctly correspond to A ⁱ _ou ⁱ _o, now d ⁱ _oremain unchanged; As renewal A ⁱ _oafter, A ⁱ _othe health status evaluation object of rope current initial damage vector d ⁱ _orepresent, A ⁱ _othe current initial Cable Structure bearing generalized coordinate vector U of bearing generalized coordinate ⁱ _orepresent.

C. C is upgraded ⁱ _omethod be: when renewal A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C ⁱ _o;

8th step: at current initial mechanical Calculation Basis model A ⁱ _obasis on, carry out several times Mechanics Calculation according to step a to steps d, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C ⁱwith evaluation object unit change vector D ⁱ _u.

A., when i-th circulation starts, directly Δ C is obtained by method listed by step b to steps d ⁱand D ⁱ _u; In other moment, when in the 7th step to A ⁱ _oafter upgrading, Δ C must be regained by method listed by step b to steps d ⁱand D ⁱ _uif, not to A in the 7th step ⁱ _oupgrade, then directly proceed to the 9th step herein and carry out follow-up work.

B. at current initial mechanical Calculation Basis model A ⁱ _obasis on carry out several times Mechanics Calculation, vectorial d ⁱ _orepresent A ⁱ _othe health status of evaluation object, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; Calculate hypothesis each time and only have an evaluation object at vectorial d ⁱ _othe basis of the health status of the evaluation object represented there is 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 is at vectorial d ⁱ _othe basis that this support cable represented has a damage there is unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again, 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, 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; 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 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; When supposing that a kth evaluation object has unit damage or load unit change, available C ⁱ _tkrepresent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target; C ⁱ _tkdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tkelement and C _oelement one_to_one corresponding.

C. the vectorial C calculated each time ⁱ _tkdeduct vectorial C ⁱ _oobtain a vector, then by each element of this vector divided by the unit damage supposed in this calculating or load unit change numerical value D ⁱ _ukafter obtain " the numerical value change of a monitored amount vector δ C ⁱ _k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ".

D. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively ⁱ"; Unit damage monitored numerical quantity transformation matrices Δ C ⁱeach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ⁱevery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of row and vectorial d _othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of coding rule and M monitored amount of row identical.

9th step: set up linear relationship error vector e ⁱwith vectorial g ⁱ.Utilize (" the monitored amount current initial value vector C of data above ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ"); while the 8th step calculates each time; namely while calculating each time in hypothesis evaluation object and only having the increase unit damage of an evaluation object or load unit change; when hypothesis kth (k=1,2,3; ..., N), when individual evaluation object increases unit damage or load unit change, calculate composition injury vector each time, use d ⁱ _tkrepresent this injury vector, corresponding monitored amount calculation current vector is C ⁱ _tk(see the 8th step), injury vector d ⁱ _tkelement number equal the quantity of evaluation object, vectorial d ⁱ _tkall elements in only have the numerical value of an element to get to calculate each time in hypothesis increase the unit damage of evaluation object or the load unit changing value of unit damage or load unit change, d ⁱ _tkthe numerical value of other element get 0, that be not numbering and the supposition of the element of 0 increase the evaluation object that unit damage or load unit change corresponding relation, be identical with the element of the same numbering of other vectors with the corresponding relation of this evaluation object; d ⁱ _tkwith evaluation object initial damage vector d _oelement number rule identical, d ⁱ _tkelement and d _oelement be one-to-one relationship.By C ⁱ _tk, C ⁱ _o, Δ C ⁱ, d ⁱ _tkbring formula (23) into, obtain a linear relationship error vector e ⁱ _k, calculate a linear relationship error vector e each time ⁱ _k; e ⁱ _ksubscript k represent kth (k=1,2,3 ..., N) and individual evaluation object increases unit damage or load unit changes.There is N number of evaluation object just to have N calculating, just have N number of linear relationship error vector e ⁱ _k, by this N number of linear relationship error vector e ⁱ _kobtaining a vector after addition, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N ⁱ.Vector g ⁱequal final error vector e ⁱ.By vectorial g ⁱbe kept on the hard disc of computer of operation health monitoring systems software, for health monitoring systems software application.

$e_k^i=abs({\mathrm{\ΔC}}^i\·d_{tk}^i-C_{tk}^i+C_o^i)---\left(23\right)$

Tenth step: define current nominal fatigue vector d ⁱ _cwith current actual damage vector d ⁱ, d ⁱ _cand d ⁱelement number equal the quantity of evaluation object, d ⁱ _cand d ⁱelement and evaluation object between be one-to-one relationship, d ⁱ _cand d ⁱelement numerical value represent degree of injury or the load change degree of corresponding evaluation object, d ⁱ _cand d ⁱwith evaluation object initial damage vector d _oelement number rule identical, d ⁱ _celement, d ⁱelement and d _oelement be one-to-one relationship.

11 step: according to monitored amount current value vector C ⁱwith " monitored amount current initial value vector C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " current nominal fatigue vector d ⁱ _c" between exist linear approximate relationship, this linear approximate relationship can be expressed as formula (11), according to multi-objective optimization algorithm calculate current nominal fatigue vector d ⁱ _cnoninferior solution, 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 of multi-objective optimization algorithm (Goal Attainment Method) can be adopted to solve current nominal fatigue vector d ⁱ _c, according to Objective Programming, formula (11) can transform the multi-objective optimization question shown in an accepted way of doing sth (24) and formula (25), and in formula (24), γ is a real number, and R is real number field, and area of space Ω limits vectorial d ⁱ _cspan (the present embodiment requirements vector d of each element ⁱ _ceach element be not less than 0, be not more than 1).Formula (24) be meant to the minimum real number γ of searching one, formula (25) is met.G (d in formula (25) ⁱ _c) defined by formula (25), the middle G (d of the product representation formula (25) of weighing vector W and γ in formula (25) ⁱ _c) and vectorial g ⁱbetween allow deviation, g ⁱdefinition see formula (17), its value calculates in the 9th step.During actual computation vector W can with vectorial g ⁱidentical.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of current nominal fatigue vector d ⁱ _c.

minimize γ

(24)

$\mathrm{\γ}\∈R,d_c^i\∈\mathrm{\Ω}$

$G\left(d_c^i\right)-W\mathrm{\γ}\≤g^i---\left(25\right)$

$G\left(d_c^i\right)=abs({\mathrm{\ΔC}}^i\·d_c^i-C^i+C_o^i)---\left(26\right)$

12 step: according to cable system current actual damage vector d ⁱdefinition (see formula (18)) and the definition (see formula (19)) of its element calculate current actual damage vector d ⁱeach element, thus can by d ⁱdetermine the health status of evaluation object.Current actual damage vector d ⁱa kth element d ⁱ _krepresent the current actual health status of a kth evaluation object in i-th circulation.

D ⁱ _krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system ⁱ _krepresent its current actual damage, d ⁱ _krepresent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion.

D ⁱ _krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is load, so a d ⁱ _krepresent its current real load change numerical value, so according to evaluation object current actual damage vector d ⁱimpaired and the degree of injury of which support cable can be defined, define which load and there occurs change and numerical value thereof.

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: set up mark vector B according to formula (20) ⁱ, formula (21) gives mark vector B ⁱthe definition of a kth element; If mark vector B ⁱelement be 0 entirely, then get back to the 6th step and proceed health monitoring to cable system and calculating; If mark vector B ⁱelement be not 0, then after completing subsequent step entirely, enter and circulate next time.

16 step: according to formula (22) calculate next time (namely the i-th+1 time, i=1,2,3,4 ...) needed for circulation initial damage vector d ⁱ⁺¹ _oeach element d ⁱ⁺¹ _ok(k=1,2,3 ..., N); The second, at initial mechanical Calculation Basis model A _obasis on, 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, then makes the health status of rope be d ⁱ⁺¹ _oafter obtain be exactly next time, namely the i-th+1 time (i=1,2,3,4 ...) Mechanics Calculation benchmark model A needed for circulation ⁱ⁺¹; Next time (namely the i-th+1 time, i=1,2,3,4 ...) current initial Cable Structure bearing generalized coordinate vector U needed for circulation ⁱ⁺¹ _oequal U ⁱ _o.Obtain A ⁱ⁺¹, d ⁱ⁺¹ _oand U ⁱ⁺¹ _oafter, obtain A by Mechanics Calculation ⁱ⁺¹in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C ⁱ⁺¹ _o.

17 step: get back to the 6th step, starts by the circulation of the 6th step to the 17 step.

Claims (1)

1. simplify generalized displacement strain monitoring load damaged cable progressive-type recognition method, 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; If total M in cable system ₁root support cable, the quantity of obvious core evaluation object is exactly M ₁; Determine the point being monitored of specifying, namely point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified points numbering; Determine monitored should the changing direction of point being monitored, and give all monitored strain numberings of specifying, " monitored strain numbering " will be used for generating vector sum matrix in subsequent step, and " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; The external force that object, structure are born can be described as load, and load 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, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular 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,

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 _ocable Structure bearing generalized coordinate data be exactly initial Cable Structure bearing generalized coordinate vector U _o; 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; 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, A in the method _o, C _o, d _oand U _oconstant;

E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting ⁱ _o; When i-th circulation starts, corresponding to A ⁱ _o" Cable Structure bearing generalized coordinate data " with current initial Cable Structure bearing generalized coordinate vector U ⁱ _orepresent, vectorial U ⁱ _odefinition mode and vectorial U _odefinition mode identical, U ⁱ _oelement and U _oelement one_to_one corresponding; The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d ⁱ _o, d ⁱ _ocable Structure A when representing that this circulation starts ⁱ _othe health status of evaluation object, d ⁱ _odefinition mode and d _odefinition mode identical, d ⁱ _oelement and d _oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C ⁱ _orepresent, vectorial C ⁱ _odefinition mode and vectorial C _odefinition mode identical, C ⁱ _oelement and C _oelement one_to_one corresponding, monitored amount current initial value vector C ⁱ _orepresent and correspond to A ⁱ _othe concrete numerical value of all monitored amount; U ⁱ _oand d ⁱ _oa ⁱ _ocharacterisitic parameter, C ⁱ _oby A ⁱ _omechanics Calculation result composition; When first time, circulation started, A ⁱ _obe designated as A ¹ _o, set up A ¹ _omethod for making A ¹ _oequal A _o; When first time, circulation started, U ⁱ _obe designated as U ¹ _o, set up U ¹ _omethod for making U ¹ _oequal U _o; When first time, circulation started, d ⁱ _obe designated as d ¹ _o, set up d ¹ _omethod for making d ¹ _oequal d _o; When first time, circulation started, C ⁱ _obe designated as C ¹ _o, set up C ¹ _omethod for making C ¹ _oequal C _o;

F. from entering the circulation being walked to q step by f 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 ⁱ, vectorial U ⁱdefinition mode and vectorial U _odefinition mode identical, U ⁱelement and U _oelement one_to_one corresponding; Vectorial U is obtained in actual measurement ⁱwhile, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C ⁱ, vectorial C ⁱdefinition mode and vectorial C _odefinition mode identical, C ⁱelement and C _oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time;

G. according to current cable structure actual measurement bearing generalized coordinate vector U ⁱ, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g3 ⁱ _o, monitored amount current initial value vector C ⁱ _owith current initial Cable Structure bearing generalized coordinate vector U ⁱ _o, and evaluation object current initial damage vector d ⁱ _oremain unchanged;

G1. U is compared ⁱwith U ⁱ _oif, U ⁱequal U ⁱ _o, then A ⁱ _o, C ⁱ _oand U ⁱ _oremain unchanged, otherwise need to follow these steps to A ⁱ _o, C ⁱ _oand U ⁱ _oupgrade;

G2. U is calculated ⁱwith U _odifference, U ⁱwith 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 ⁱdeduct U _o;

G3. 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 _oin Cable Structure bearing apply generalized displacement of support constraint after obtain upgrade current initial mechanical Calculation Basis model A ⁱ _o, upgrade A ⁱ _owhile, U ⁱ _oall elements numerical value also uses U ⁱall elements numerical value correspondence replaces, and namely have updated U ⁱ _o, so just obtain and correctly correspond to A ⁱ _ou ⁱ _o, now d ⁱ _oremain unchanged; As renewal A ⁱ _oafter, A ⁱ _othe health status evaluation object of rope current initial damage vector d ⁱ _orepresent, A ⁱ _othe current initial Cable Structure bearing generalized coordinate vector U of bearing generalized coordinate ⁱ _orepresent; Upgrade C ⁱ _omethod be: when renewal A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C ⁱ _o;

H. at current initial mechanical Calculation Basis model A ⁱ _obasis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C ⁱwith evaluation object unit change vector D ⁱ _u;

H1., when i-th circulation starts, directly Δ C is obtained by method listed by step h2 to step h4 ⁱand D ⁱ _u; In other moment, when in step g to A ⁱ _oafter upgrading, Δ C must be regained by method listed by step h2 to step h4 ⁱand D ⁱ _uif, not to A in step g ⁱ _oupgrade, then directly proceed to step I herein and carry out follow-up work;

H2. at current initial mechanical Calculation Basis model A ⁱ _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 suppose that this support cable increases unit damage again, if this evaluation object is a load, just suppose that this load 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 again different from during other time calculates the evaluation object increasing unit damage or load unit change again, 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; When supposing that a kth evaluation object increases unit damage or load unit change again, use C ⁱ _tkrepresent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target; C ⁱ _tkdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tkelement and C _oelement one_to_one corresponding;

H3. the vectorial C calculated each time ⁱ _tkdeduct vectorial C ⁱ _oobtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated divided by this unit damage or load unit change numerical value supposed ⁱ _k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ";

H4. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively ⁱ"; Unit damage monitored numerical quantity transformation matrices Δ C ⁱeach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ⁱevery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of row and vectorial d _othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ⁱthe coding rule of coding rule and M monitored amount of row identical;

I. current nominal fatigue vector d is defined ⁱ _cwith current actual damage vector d ⁱ, d ⁱ _cand d ⁱelement number equal the quantity of evaluation object, d ⁱ _cand d ⁱelement and evaluation object between be one-to-one relationship, d ⁱ _celement numerical value represent nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object, d ⁱ _cand d ⁱwith evaluation object initial damage vector d _oelement number rule identical, d ⁱ _celement, d ⁱelement and d _oelement be one-to-one relationship;

J. according to monitored amount current value vector C ⁱwith " monitored amount current initial value vector C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " current nominal fatigue vector d ⁱ _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1 ⁱ _cother outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d ⁱ _c;

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

K. the current actual damage vector d utilizing formula 2 to express ⁱa kth element d ⁱ _kwith evaluation object current initial damage vector d ⁱ _oa kth element d ⁱ _okwith current nominal fatigue vector d ⁱ _ca kth element d ⁱ _ckbetween relation, calculate current actual damage vector d ⁱall elements;

formula 2

K=1 in formula 2,2,3 ..., N; d ⁱ _krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system ⁱ _krepresent its current actual damage, d ⁱ _krepresent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; So according to evaluation object current actual damage vector d ^athe health status of core evaluation object can be determined;

L. current nominal fatigue vector d is tried to achieve ⁱ _cafter, set up mark vector B according to formula 3 ⁱ, formula 4 gives mark vector B ⁱthe definition of a kth element;

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

formula 4

Element B in formula 4 ⁱ _kmark vector B ⁱa kth element, D ⁱ _ukevaluation object unit change vector D ⁱ _ua kth element, d ⁱ _ckevaluation object current nominal fatigue vector d ⁱ _ca kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;

If m. mark vector B ⁱelement be 0 entirely, then get back to step f continue this circulation; If mark vector B ⁱelement be not 0 entirely, then enter next step, i.e. step n;

N. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the i-th+1 time circulation ⁱ⁺¹ _oeach element;

formula 5

D in formula 5 ⁱ⁺¹ _okthe evaluation object current initial damage vector d next time, namely needed for the i-th+1 time circulation ⁱ⁺¹ _oa kth element, d ⁱ _okthis, i.e. the evaluation object of i-th circulation current initial damage vector d ⁱ _oa kth element, D ⁱ _ukthe evaluation object unit change vector D of i-th circulation ⁱ _ua kth element, B ⁱ _kthe mark vector B of i-th circulation ⁱa kth element, k=1 in formula 5,2,3 ..., N;

O. at initial mechanical Calculation Basis model A _obasis on, first 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, then makes the health status of rope be d ⁱ⁺¹ _oafter obtain be exactly next time, namely the i-th+1 time circulation needed for Mechanics Calculation benchmark model A ⁱ⁺¹; Obtain A ⁱ⁺¹after, obtain A by Mechanics Calculation ⁱ⁺¹in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C ⁱ⁺¹ _o;

P. take off once, current initial Cable Structure bearing generalized coordinate vector U namely needed for the i-th+1 time circulation ⁱ⁺¹ _oequal the current initial Cable Structure bearing generalized coordinate vector U of i-th circulation ⁱ _o;

Q. get back to step f, start to circulate next time.