CN102706607B - Based on the approximant recognition methods of the damaged cable of angle monitor during temperature variation - Google Patents

Abstract

During temperature variation based on the approximant recognition methods of the damaged cable of angle monitor based on angle monitor, determined whether the Mechanics Calculation benchmark model that needs to upgrade Cable Structure by monitoring Cable Structure temperature, environment temperature and damaged cable degree of injury, only has the Mechanics Calculation benchmark model just upgrading Cable Structure when temperature changes and (or) damaged cable degree of injury is larger, thus obtain the new Mechanics Calculation benchmark model counting the Cable Structure of temperature, the basis of this model calculates and obtains unit damage monitored numerical quantity transformation matrices.According to the current value vector of monitored amount with monitored amount current initial value vector, the linear approximate relationship that exists between unit damage monitored numerical quantity transformation matrices and current nominal fatigue vector to be asked, multi-objective optimization algorithm scheduling algorithm can be utilized to calculate the noninferior solution of current nominal fatigue vector fast, accordingly can when there being temperature variation, determine position and the degree of injury thereof of damaged cable more exactly.

Description

Based on the approximant recognition methods of the damaged cable of angle monitor during temperature variation

Technical field

The structures such as cable-stayed bridge, suspension bridge, 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 only bears the rod member of tensile load for support unit, and such structure representation is " Cable Structure " by this method for simplicity.Along with the change of environment temperature, the temperature of Cable Structure also can change, when Cable Structure temperature changes, based on angle monitor, this method identifies that the supporting system of Cable Structure (refers to all ropeway carrying-ropes, and all rod members only bearing tensile load play supporting role, for simplicity, whole support units of this class formation are collectively referred to as " cable system " by this patent, but in fact cable system not only refers to support cable, also the rod member only bearing tensile load is comprised, also for ease of conveniently, all ropeway carrying-ropes and all rod members only bearing tensile load play supporting role is censured with " support cable " this noun in this method, also for ease of conveniently, censure all ropeway carrying-ropes and all rod members only bearing tensile load play supporting role with " support cable " this noun in this method) in damaged cable (the impaired rod member only bearing tensile load is just referred to truss-frame structure), belong to engineering structure health monitoring field.

Background technology

The key components of cable system normally Cable Structure, its inefficacy usually brings the inefficacy of total, and the damaged cable (also referring to only bear the rod member of tensile load as previously mentioned) that structure based health monitoring technique identifies in the cable system of Cable Structure is a kind of method of great potential.After the health status of cable system changes, the change of the measurable parameter of structure can be caused, such as caused the change of the angle coordinate of any imaginary line of the every bit of the Cable Structure (change of the angle coordinate of the straight line of any this point of mistake in the section of such as body structure surface any point, or the change of the angle coordinate of the normal of body structure surface any point), in fact the change of structural point contains the health status information of cable system, that is structural point data can be utilized to judge the health status of structure, can (monitored angle-data be called " monitored amount " by this method based on angle monitor, mention " monitored amount " below and just refer to monitored angle-data) identify damaged cable.Monitored amount is except being subject to the impact of cable system health status; also can be subject to the impact of Cable Structure temperature variation (usually can occur); under the condition that Cable Structure temperature changes; if damaged cable identification can be realized based on to the monitoring of monitored amount; to the safety of Cable Structure, there is important value, also do not have a kind of disclosed, effective health monitoring systems and method to solve this problem at present.

When Cable Structure has temperature variation, in order to can monitoring reliably be had to the health status of the cable system of Cable Structure and judge, must have one can rationally effectively set up each monitored amount change with ropes all in cable system health status between the method for relation, the health monitoring systems set up based on the method can provide the health evaluating of more believable cable system.

Summary of the invention

Technical matters: the object of this method is when Cable Structure has temperature variation, for the health monitoring problem of cable system in Cable Structure, disclose a kind of based on angle monitor, health monitor method that cable system in Cable Structure can be monitored rationally and effectively.

Technical scheme: this method is made up of three parts.Respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; Based on the approximant recognition methods of damaged cable of angle monitor during temperature variation based on the angle of the method for the knowledge base needed for the approximant recognition system of damaged cable of angle monitor and parameter, knowledge based storehouse (containing parameter) and actual measurement Cable Structure and actual measurement Cable Structure temperature when two, setting up temperature variation; Three, during temperature variation based on the software and hardware part of the approximant recognition system of the damaged cable of angle monitor.

The Part I of this method: " the temperature survey calculating method of the Cable Structure of this method ".

First determine " the temperature survey calculating method of the Cable Structure of this method ".Temperature due to Cable Structure may be change, the temperature of the different parts of such as Cable Structure changes along with the change of intensity of sunshine, along with the change of environment temperature changes, the surface of Cable Structure may be time dependent with inner temperature sometimes, the surface of Cable Structure may be different from inner temperature, the surface of Cable Structure is time dependent with inner temperature difference, this just makes the Mechanics Calculation of Cable Structure when considering temperature conditions and monitoring quite complicated, for simplifying problem, reduce calculated amount and reduce and measure cost, especially in order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", specific as follows:

The first step, inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day _r.Be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of describing _hunit be DEG C/m.The surface of Cable Structure is got " R Cable Structure surface point ", the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " describes later with the condition that must meet that distributes.From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, be called reference plate, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate.Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, the moment of survey record data is called the physical record data moment.

Second step, Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T _emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T _pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T _smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _tmax.

3rd step, survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, reference plate maximum temperature difference Δ T _pmaxwith Cable Structure surface maximum temperature difference Δ T _smaxall be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxΔ T is not more than after deducting 2 degrees Celsius _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _tmaxbe not more than 1 degree Celsius, this method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtained, gets Δ T for convenience of describing _hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure.

The Part II of this method: based on the approximant recognition methods of damaged cable of angle monitor during temperature variation based on the angle of the method for the knowledge base needed for the approximant recognition system of damaged cable of angle monitor and parameter, knowledge based storehouse (containing parameter) and actual measurement Cable Structure and actual measurement Cable Structure temperature when setting up temperature variation.Can carry out successively as follows, to obtain the health state evaluation of cable system more accurately.

The first step: establish total N root support cable, first determine the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step.

" the whole monitored angle-data of structure " by the specified point of K in structure, cross L of each specified point appointment straight line, H angle coordinate component of each appointment straight line describe, the change of structural point is exactly the change of all angle coordinate components of specifying of all specified points, all appointment straight line.Each total M (M=K × L × H) individual angle coordinate component measurement value or calculated value carry out the angle information of characterisation of structures.K and M must not be less than the quantity N of support cable.

For simplicity, in the method by " the monitored angle-data of structure " referred to as " monitored amount ".When mentioning " so-and-so matrix of monitored amount or so-and-so vector " later, also can be read as " monitored angle so-and-so matrix or so-and-so vector ".

Set up initial mechanical Calculation Basis model A _otime, when Cable Structure is completed, or before setting up health monitoring (damaged cable identification) system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T _orepresent, be called initial Cable Structure steady temperature data vector T _o.T is obtained in actual measurement _owhile, use the direct survey calculation of conventional method to obtain the initial number of all monitored amount of Cable Structure.Initial Cable Structure steady temperature data vector T is obtained at Actual measurement _owhile, use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained at Actual measurement _owhile, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, structural strain data, structural point measurement data, structure space measurement of coordinates data that the Non-destructive Testing Data of support cable etc. can express the health status of rope.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.The data utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope set up cable system initial damage vector d _o(as the formula (1)), d is used _orepresent that Cable Structure is (with initial mechanical Calculation Basis model A _orepresent) the initial health of cable system.If there is no the Non-destructive Testing Data of rope and other are when can express the data of the health status of support cable, or when can think that structure original state is not damaged state, vectorial d _oeach element numerical value get 0.The temperature variant 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 steady temperature data vector T _o, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A _o.

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

D in formula (1) _oj(j=1,2,3 ...., N) represent initial mechanical Calculation Basis model A _oin the initial damage value of jth root rope of cable system, d _ojrepresent jth root rope not damaged when being 0, represent that this rope thoroughly loses load-bearing capacity when being 100%, represent that the load-bearing capacity of corresponding proportion lost by jth root rope time between 0 and 100%, T represents the transposition (same afterwards) of vector.

T is obtained in actual measurement _owhile, 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 _ok···C _oM] ^T(2)

C in formula (2) _ok(k=1,2,3 ...., M) be a kth monitored amount in Cable Structure.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, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T _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 cable system initial damage vector d of middle support cable _orepresent, the initial Cable Structure steady temperature data vector T of Cable Structure Cable Structure steady temperature data _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.T can be said _oand d _oa _oparameter, C _oby A _omechanics Calculation result composition.

Second step: circulation starts.When circulation starts each time, first need the cable system 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 ⁱ _otemperature Distribution " current initial Cable Structure steady temperature data vector T ⁱ _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.A _oand A ⁱ _ocount temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure.

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

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

D in formula (3) ⁱ _oj(i=1,2,3, J=1,2,3 ...., N) represent i-th time circulation start time, current initial mechanical Calculation Basis model A ⁱ _oin the initial damage value of jth root rope of cable system, d ⁱ _ojrepresent jth root rope not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, time between 0 and 100%, represent that the load-bearing capacity of corresponding proportion lost by jth root rope.

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 cable system ¹ _o) time, d ¹ _ojust equal d _o.I-th (i=2,3,4,5,6 ...) the secondary cable system 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 A ⁱ _o" Cable Structure steady temperature data " use vector T ⁱ _orepresent, be called current initial Cable Structure steady temperature data vector T ⁱ _o.Vector T ⁱ _odefinition mode and vector T _odefinition mode identical, must set up or set up and be called current initial Cable Structure steady temperature data vector T when circulation starts each time ⁱ _o.

Set up, upgrade A ⁱ _oand T ⁱ _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, T ¹ _oequal T _o.A in circulation each time ⁱ _oand T ⁱ _obe constantly update, concrete grammar is described below; At the end of circulation each time, upgrade A ⁱ _oand T ⁱ _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_{\mathrm{ok}}^i& \·& \·& \·& C_{\mathrm{oM}}^i\end{array}\right]}^T---\left(4\right)$

C in formula (2) ⁱ _ok(i=1,2,3, K=1,2,3 ...., M) be kth 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.

T can be said ⁱ _oand d ⁱ _oa ⁱ _ocharacterisitic parameter, C ⁱ _oa ⁱ _oat T ⁱ _oand d ⁱ _omechanics Calculation result composition under condition.

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, T ⁱ _o, C ⁱ _oand d ⁱ _oafter, the current data of " Cable Structure steady temperature data " is obtained, current data composition " the current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement ⁱ", vector T ⁱdefinition mode and vector T _odefinition mode identical; In actual measurement vector T ⁱ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 ⁱ.C ⁱelement and C _oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.

In acquisition vector T ⁱafter, upgrade A according to following concrete grammar ⁱ _o, T ⁱ _o, C ⁱ _oand d ⁱ _o:

Relatively T ⁱand T ⁱ _oif, T ⁱequal T ⁱ _o, then do not need A ⁱ _oupgrade; If T ⁱbe not equal to T ⁱ _o, then need A ⁱ _oupgrade, now T ⁱwith T _odifference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ⁱwith T _odifference represent with steady temperature change vector S, S equals T ⁱdeduct T _o, S represents the change of Cable Structure steady temperature data.Upgrade A ⁱ _omethod be: at A _obasis on make the health status of rope be cable system current initial damage vector d ⁱ _o, more further to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A _oin Cable Structure applies to obtain after temperature variation is exactly the current initial mechanical Calculation Basis model A of renewal ⁱ _o, upgrade A ⁱ _owhile, T ⁱ _oall elements numerical value also uses T ⁱall elements numerical value replaces, and namely have updated current initial Cable Structure steady temperature data vector T ⁱ _o, so just obtain and correctly correspond to A ⁱ _ot ⁱ _o, now d ⁱ _oremain unchanged, so far just achieve A ⁱ _orenewal.As renewal A ⁱ _oafter, now A ⁱ _othe health status cable system of rope current initial damage vector d ⁱ _orepresent, A ⁱ _ocable Structure steady temperature current cable structure steady temperature data vector T ⁱrepresent, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, replace C with these concrete numerical value ⁱ _othe element of middle correspondence, so just achieves monitored amount current initial value vector C ⁱ _orenewal.

4th step: circulation time must first be set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage 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 ...)." nominal unit damage 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 steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _oafter, upgrade unit damage monitored numerical quantity transformation matrices Δ C ⁱwith nominal unit damage 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 nominal unit damage 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 nominal unit damage 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 nominal unit damage 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 each time in hypothesis cable system and only have a rope on the basis of original damage (original damage can be 0, also can not be 0), increase unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again.For convenience of calculating, can be all structural health conditions when this circulation is started each time when setting unit damage in circulation as being completely healthy, and set on this basis unit damage (in subsequent step, calculate, the damage numerical value of rope---be called nominal fatigue d ⁱ _c(i=1,2,3 ...), all relative to when this circulation is started, by the health status of rope 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 rope damaged is different from during other time calculates the rope occurring damaging, and the unit damage value that suppose there is the rope of damage each time can be different from the unit damage value of other ropes, with " the vectorial D of nominal unit damage ⁱ _u" (as the formula (5)) record the unit damage of supposition of all ropes 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 gauge calculates current value vector ", and (when supposing that jth root rope has unit damage, available formula (6) represents that the monitored gauge of all M specified a monitored amount calculates current value vector C ⁱ _tj); The monitored gauge calculated each time is calculated current value vector and is deducted 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 jth root rope has unit damage, uses δ C under this condition ⁱ _jrepresent the numerical value change vector of monitored amount, δ C ⁱ _jdefinition see formula (7), formula (8) and formula (9), formula (7) deducts after formula (4) again divided by vectorial D for formula (6) ⁱ _ua jth element D ⁱ _ujgained), the numerical value change vector δ C of monitored amount ⁱ _jeach element representation owing to suppose there is unit damage (the such as D of the Na Gensuo (such as jth root rope) of unit damage when calculating ⁱ _uj), and the numerical value knots modification of monitored amount corresponding to this element caused is relative to the unit damage D of supposition ⁱ _ujrate of change; N root rope 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 ⁱ _j(j=1,2,3 ...., N) be matrix Δ C ⁱrow, Δ C ⁱdefinition as the formula (10).

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

Nominal unit damage vector D in formula (5) ⁱ _uelement D ⁱ _uj(i=1,2,3, J=1,2,3 ...., N) represent the unit damage numerical value of jth root rope of supposition in i-th circulation, vectorial D ⁱ _uin the numerical value of each element can be the same or different.

$C_{\mathrm{tj}}^i={\left[\begin{array}{cccccccccc}{C}_{\mathrm{tk}1}^i& C_{\mathrm{tk}2}^i& \·& \·& \·& C_{\mathrm{tjk}}^i& \·& \·& \·& C_{\mathrm{tjM}}^i\end{array}\right]}^T---\left(6\right)$

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

$\mathrm{\δ}{C}_j^i=\frac{C_{\mathrm{tj}}^i-C_o^i}{D_{\mathrm{uj}}^i}---\left(7\right)$

In formula (7) each amount subscript i (i=1,2,3 ...) and represent i-th circulation, subscript j (j=1,2,3 ...., N) represent that jth root rope has unit damage, D in formula ⁱ _ujvectorial D ⁱ _uin a jth element.Vector δ C ⁱ _jdefinition such as formula shown in (8), δ C ⁱ _jkth (k=1,2,3 ...., M) individual element δ C ⁱ _jkrepresent in i-th circulation, set up matrix Δ C ⁱtime, assuming that jth root rope has the unit damage D of knots modification relative to supposition calculating a gained kth monitored amount during unit damage ⁱ _ujrate of change, its definition such as formula shown in (9).

$\mathrm{\δ}{C}_j^i={\left[\begin{array}{cccccccccc}\mathrm{\δ}{C}_{j1}^i& \mathrm{\δ}{C}_{j2}^i& \·& \·& \·& \mathrm{\δ}{C}_{\mathrm{jk}}^i& \·& \·& \·& {\mathrm{\δC}}_{\mathrm{jM}}^i\end{array}\right]}^T---\left(8\right)$

$\mathrm{\δ}{C}_{\mathrm{jk}}^i=\frac{C_{\mathrm{tjk}}^i-C_{\mathrm{ok}}^i}{D_{\mathrm{uj}}^i}---\left(9\right)$

In formula (9), the definition of each amount has been previously described.

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

Vectorial δ C in formula (10) ⁱ _j(i=1,2,3 ...., j=1,2,3 ...., N) represent in i-th circulation, because jth root rope has unit damage D ⁱ _ujcause, the change of the relative value of all monitored amounts.Matrix Δ C ⁱthe coding rule of row (subscript j) and vectorial d above ⁱ _othe coding rule of subscript j of element identical.

5th step: the current health state identifying cable system.Detailed process is as follows.

I-th (i=1,2,3 ...) in secondary circulation, utilize " the monitored amount current value vector C obtained in the 3rd 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); Cable system current nominal fatigue vector d ⁱ _cdefinition see formula (14).

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

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

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

D in formula (14) ⁱ _cj(i=1,2,3 ....; J=1,2,3 ...., N) be the current nominal fatigue value of cable system jth root rope in i-th circulation, vectorial d ⁱ _cthe coding rule of subscript j of element and matrix Δ C ⁱthe coding rule of row identical.

When rope 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=\mathrm{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 can not simply according to formula (11) or formula (12) and " monitored amount current value vector C ⁱ" come direct solution obtain Suo Dangqian nominal fatigue vector d ⁱ _c.If this has been doned, the injury vector d obtained ⁱ _cin element even there will be larger negative value, namely negative damage, this is obviously irrational.Therefore rope injury vector d is obtained ⁱ _cacceptable solution (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more accurately from cable system) become a rational solution, available formula (16) expresses this method.

$\mathrm{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& \text{\·}& \·& \·& g_k^i& \·& \·& \·& g_M^i\end{array}\right]}^T---\left(17\right)$

G in formula (17) ⁱ _k(i=1,2,3 ....; K=1,2,3 ...., M) 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 cable system current nominal fatigue vector d ⁱ _cacceptable solution, cable system 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 position and the degree of injury of damaged cable, namely achieve the health monitoring of cable system, achieve damaged cable identification.

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

D in formula (18) ⁱ _j(i=1,2,3, J=1,2,3 ...., N) represent the actual damage value of jth root rope in i-th circulation, its definition is shown in formula (19), d ⁱ _jrepresent jth root rope not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, represent time between 0 and 100% that the load-bearing capacity of corresponding proportion lost by jth root rope, vectorial d ⁱthe coding rule of element and formula (1) in vectorial d _othe coding rule of element identical.

$d_j^i=1-(1-d_{\mathrm{oj}}^i)(1-d_{\mathrm{cj}}^i)---\left(19\right)$

D in formula (19) ⁱ _oj(i=1,2,3,4, J=1,2,3 ...., N) be cable system current initial damage vector d ⁱ _oa jth element, d ⁱ _cjcable system current nominal fatigue vector d ⁱ _ca jth element.

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 jth element; If mark vector B ⁱelement be 0 entirely, then get back to the 3rd 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.

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 ⁱ⁺¹ _oj; The second, at Mechanics Calculation benchmark model A _obasis on, make A _oin the health status of rope be d ⁱ⁺¹ _oinstead of be d _oafter, more further to A _oin Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, and steady temperature change vector S equals T ⁱdeduct T _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 steady temperature data vector T needed for circulation ⁱ⁺¹ _oequal T ⁱ _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_j^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 ⁱ _j(j=1,2,3 ..., N) subscript j represent the damage characteristic of jth root rope, can only get 0 and 1 two amount, concrete value rule is shown in formula (21).

$B_j^i=\left\{\begin{array}{ccc}0,& \mathrm{if}& d_{\mathrm{cj}}^i<D_{\mathrm{uj}}^i\\ 1,& \mathrm{if}& d_{\mathrm{cj}}^i\&GreaterEqual;D_{\mathrm{uj}}^i\end{array}\right.---\left(21\right)$

Element B in formula (21) ⁱ _jmark vector B ⁱa jth element, D ⁱ _ujnominal unit damage vector D ⁱ _ua jth element (see formula (3)), d ⁱ _cjcable system current nominal fatigue vector d ⁱ _ca jth element (see formula (14)), they all represent the relevant information of jth root rope.

$d_{\mathrm{oj}}^{i+1}=1-(1-d_{\mathrm{oj}}^i)(1-D_{\mathrm{uj}}^i{B}_j^i)---\left(22\right)$

D in formula (22) ⁱ _ujnominal unit damage vector D ⁱ _ua jth element (see formula (5)), d ⁱ _ojcable system current initial damage vector d ⁱ _oa jth element (see formula (3)).

The Part III of this method: based on the software and hardware part of the approximant recognition system of the damaged cable of angle monitor during temperature variation.

Hardware components comprises monitoring system (comprising monitored amount monitoring system, temperature monitoring system), signal picker and computing machine etc.Require that Real-Time Monitoring obtains temperature required measured data, 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:

A. establish total N root support cable, first determine the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; Determine point being monitored, number to all points being monitored; Determined the monitored straight line of each point being monitored, to all monitored straight line numberings; Determine the monitored angle coordinate component of each monitored straight line, to all monitored angle coordinate component numberings; Above-mentioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored angle-data of Cable Structure " is made up of above-mentioned all monitored angle coordinate components; For simplicity, in the method by " the monitored angle-data of Cable Structure " referred to as " monitored amount "; The quantity of point being monitored must not be less than the quantity of support cable; The quantity sum of all monitored angle coordinate components must not be less than the quantity of support cable; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, the moment of survey record data is called the physical record data moment;

B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;

B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day _r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of describing _hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate,

B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T _emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T _pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T _smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _tmax,

B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, reference plate maximum temperature difference Δ T _pmaxwith Cable Structure surface maximum temperature difference Δ T _smaxall be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxΔ T is not more than after deducting 2 degrees Celsius _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _tmaxbe not more than 1 degree Celsius, this method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtained, gets Δ T for convenience of describing _hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure,

C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T _o"; Survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses; T is obtained in actual measurement _owhile, namely at the initial Cable Structure steady temperature data vector T of acquisition _othe synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Non-destructive Testing Data of the health status expressing support cable, 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 support coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C _o; The Non-destructive Testing Data that utilization can express the health status of support cable sets up cable system initial damage vector d _o, cable system initial damage vector d _oelement number equal N, d _oelement and support cable be one-to-one relationship, cable system initial damage vector d _oelement numerical value represent the degree of injury of corresponding support cable, if cable system initial damage vector d _othe numerical value of a certain element be 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 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 when thinking that Cable Structure original state is not damaged state, vectorial d _oeach element numerical value get 0;

Temperature variant physical and mechanical properties parameter, the initial Cable Structure steady temperature data vector T 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, the Non-destructive Testing Data of support cable, Cable Structure _owith all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " _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 _o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T _o"; Corresponding to A _osupport cable health status with cable system initial damage vector d _orepresent; Corresponding to A _othe initial value monitored amount initial value vector C of all monitored amount _orepresent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " ⁱ _o, monitored amount current initial value vector C ⁱ _o" current initial Cable Structure steady temperature data vector T ⁱ _o"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ⁱ _owith monitored amount current initial value vector C ⁱ _otime, the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _ojust equal the initial mechanical Calculation Basis model A of Cable Structure _o, monitored amount current initial value vector C ⁱ _ojust equal monitored amount initial value vector C _o; A ⁱ _ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T ⁱ _o", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ⁱ _otime, T ⁱ _ojust equal T _o; A ⁱ _othe initial health of support cable and A _othe health status of support cable identical, also use cable system initial damage vector d _orepresent, A in cyclic process below ⁱ _othe initial health of support cable use cable system initial damage vector d all the time _orepresent; T _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, work as T ⁱ _oand d _oa ⁱ _oparameter time, C ^t _oby A ⁱ _omechanics Calculation result composition; A in the method _o, C _o, d _oand T _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, A _oand A ⁱ _ocount temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When i-th circulation starts, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _orepresent, vector T ⁱ _odefinition mode and vector T _odefinition mode identical, T ⁱ _oelement and T _oelement one_to_one corresponding; The current initial damage vector of cable system 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 cable system, 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; T ⁱ _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, T ⁱ _obe designated as T ¹ _o, set up T ¹ _omethod for making T ¹ _oequal T _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 Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement ⁱ, vector T ⁱdefinition mode and vector T _odefinition mode identical, T ⁱelement and T _oelement one_to_one corresponding; Vector T is obtained in actual measurement ⁱwhile, actual measurement obtains at acquisition current cable structure steady temperature data vector T ⁱmoment synchronization Cable Structure in the currency of all monitored amounts, 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 steady temperature data vector T ⁱ, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g2 ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _o, and cable system current initial damage vector d ⁱ _oremain unchanged;

G1. T is compared ⁱand T ⁱ _oif, T ⁱequal T ⁱ _o, then A ⁱ _o, T ⁱ _o, C ⁱ _oand d ⁱ _oremain unchanged;

G2. T is compared ⁱand T ⁱ _oif, T ⁱbe not equal to T ⁱ _o, then need current initial mechanical Calculation Basis model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _oupgrade, update method is: first calculate T ⁱwith T _odifference, T ⁱwith T _odifference represent with steady temperature change vector S, S equals T ⁱdeduct T _o; At initial mechanical Calculation Basis model A _obasis on make the health status of support cable be cable system current initial damage vector d ⁱ _o, more further to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A _oin Cable Structure applies to obtain after temperature variation is exactly the current initial mechanical Calculation Basis model A of renewal ⁱ _o; At renewal A ⁱ _owhile, T ⁱ _oall elements numerical value T ⁱall elements numerical value replaces, and namely have updated T ⁱ _o, so just obtain and correctly correspond to A ⁱ _ot ⁱ _o, now d ⁱ _oremain unchanged; Upgrade A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, with these concrete numerical value substituting vector C ⁱ _othe element of middle correspondence, so just achieves vectorial C ⁱ _orenewal;

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 nominal unit damage 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 of all support cables, N root support cable is had just to have N calculating, calculating each time in hypothesis cable system only has a support cable to increase unit damage again on the basis of original damage, occur in calculating each time that the support cable damaged is different from during other time calculates the support cable occurring damaging, and the unit damage value that suppose there is the support cable of damage each time can be different from the unit damage value of other support cables, with " nominal unit damage vector D ⁱ _u" record the unit damage of the supposition of all ropes, vectorial D ⁱ _uelement number rule with vectorial d _othe coding rule of element identical, calculate the current value of all monitored amounts in Cable Structure each time, the current value of all monitored amount calculated each time forms one " monitored gauge calculates current value vector "; When supposing that jth root support cable has unit damage, available C ⁱ _tjrepresent corresponding " monitored gauge calculates current value 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 ⁱ _tjdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tjelement and C _oelement one_to_one corresponding;

H3. the vectorial C calculated each time ⁱ _tjdeduct vectorial C ⁱ _oobtain a vector, then after the unit damage value of supposition in being calculated divided by this by each element of this vector, obtain " numerical value change vector δ a C for monitored amount ⁱ _j"; N root support cable is had just to have N number of " the numerical value change vector of monitored amount ";

H4. form by this N number of " numerical value change vector of monitored amount " " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange successively ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row correspond to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d _oelement number rule identical;

I. current nominal fatigue vector d is defined ⁱ _cwith current actual damage vector d ⁱ, d ⁱ _cand d ⁱelement number equal the quantity of support cable, d ⁱ _cand d ⁱelement and support cable between be one-to-one relationship, d ⁱ _cand d ⁱelement numerical value represent degree of injury or the health status of corresponding support cable, d ⁱ _cand d ⁱwith cable system 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{\&Delta;}{C}^i\&CenterDot;d_c^i$ Formula 1

K. the current actual damage vector d utilizing formula 2 to express ⁱa jth element d ⁱ _jwith cable system current initial damage vector d ⁱ _oa jth element d ⁱ _ojwith current nominal fatigue vector d ⁱ _ca jth element d ⁱ _cjbetween relation, calculate current actual damage vector d ⁱall elements;

$d_j^i=1-(1-d_{\mathrm{oj}}^i)(1-d_{\mathrm{cj}}^i)$ Formula 2

J=1 in formula 2,2,3 ..., N; Current actual damage vector d ⁱelement numerical value represent the degree of injury of corresponding support cable, so according to current actual damage vector d ⁱjust can define the impaired and degree of injury of which support cable, namely achieve the health monitoring of cable system in Cable Structure, achieve damaged cable identification; If the numerical value of a certain element of current actual damage vector 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;

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 jth element;

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

formula 4

Element B in formula 4 ⁱ _jmark vector B ⁱa jth element, D ⁱ _ujnominal unit damage vector D ⁱ _ua jth element, d ⁱ _cjcable system current nominal fatigue vector d ⁱ _ca jth element, they all represent the relevant information of jth root support cable, j=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, cable system current initial damage vector d namely needed for the i-th+1 time circulation ⁱ⁺¹ _oeach element;

$d_{\mathrm{oj}}^{i+1}=1-(1-d_{\mathrm{oj}}^i)(1-D_{\mathrm{ju}}^i{B}_j^i)$ Formula 5

D in formula 5 ⁱ⁺¹ _ojthe cable system current initial damage vector d next time, namely needed for the i-th+1 time circulation ⁱ⁺¹ _oa jth element, d ⁱ _ojthis, i.e. the cable system of i-th circulation current initial damage vector d ⁱ _oa jth element, D ⁱ _ujthe nominal unit damage vector D of i-th circulation ⁱ _ua jth element, B ⁱ _jthe mark vector B of i-th circulation ⁱa jth element, j=1 in formula 5,2,3 ..., N;

O. take off once, namely the i-th+1 time circulation needed for current initial Cable Structure steady temperature data vector T ⁱ⁺¹ _oequal the current initial Cable Structure steady temperature data vector T of i-th circulation ⁱ _o;

P. at initial mechanical Calculation Basis model A _obasis on, make the health status of rope be d ⁱ⁺¹ _oafter, to initial mechanical Calculation Basis model A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just get this, the steady temperature change vector S that namely uses in i-th circulation, to A _oin Cable Structure applies to obtain after temperature variation is exactly Mechanics Calculation benchmark model A next time, namely needed for the i-th+1 time circulation ⁱ⁺¹; 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;

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

Beneficial effect: when the temperature field of Cable Structure is subject to affecting of the factor such as sunshine and environment temperature, the temperature field of Cable Structure is constantly change, the change of temperature field of Cable Structure must affect the monitored amount of Cable Structure, only have and monitored amount is rejected could carry out rational cable structure health monitoring based on monitored amount by temperature profile effect part, and the temperature field measurement of Cable Structure and calculating are very complicated, this method discloses to comprise and is a kind ofly suitable for the simple of cable structure health monitoring, economical, feasible, the cable structure health monitoring method of efficient Cable Structure Calculation Method of Temperature Field, adopt this method many ropes of Cable Structure synchronous impaired time, and the temperature of Cable Structure along with time variations time, very monitor assessment can identify the health status (comprising position and the degree of injury of all damaged cables) of cable system, the effective health monitoring of system and method disclosed in this method to cable system is highly profitable.

Embodiment

When there being temperature variation, for the health monitoring of the cable system of Cable Structure, this method discloses a kind of system and method can monitoring the health status identifying each root rope in cable system in Cable Structure rationally and effectively.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.

This method adopts a kind of algorithm, and this algorithm is for monitoring the health status of the cable system identified in Cable Structure.During concrete enforcement, the following step is the one in the various steps that can take.

The first step: determine " the temperature survey calculating method of the Cable Structure of this method ", the method concrete steps are as follows:

A walks: inquiry or actual measurement (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure) obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model (such as finite element model) of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day _r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of describing _hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual observation record, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, if the sea level elevation of such as Cable Structure is between 0m to 200m, so can choose height above sea level 0m, 50m, 100m and height above sea level 200m, crossing with Cable Structure surface with imaginary surface level at the sea level elevation place that each is chosen, obtain intersection, surface level is crossing with Cable Structure obtains cross surface, intersection is the outer edge line of cross surface, 6 points are chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness.In 6 directions of the measurement Cable Structure chosen along the Temperature Distribution of wall thickness, first according to the meteorological data throughout the year in region, Cable Structure position and the physical dimension of Cable Structure, volume coordinate, Cable Structure surrounding environment etc. determines the sunny slope of Cable Structure and in the shade, the sunny slope of Cable Structure and in the shade face are the parts on the surface of Cable Structure, at the sea level elevation place that each is chosen, aforementioned intersection respectively has one section in sunny slope and in the shade, these two sections of intersection respectively have a mid point, cross the outer normal that these two mid points get Cable Structure, these two outer normals are called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure by this method, these two outer normal directions are called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure by this method, the outer normal of obvious sunny slope and the outer normal of in the shade all crossing with aforementioned intersection, also two intersection points are just had, intersection is divided into two line segments by these two intersection points, 2 points are got respectively on two line segments, totally 4 points, each line segment in two of intersection line segments is divided into equal 3 sections of length by taken point, the outer normal on Cable Structure surface is got at these 4 some places, the outer normal on 6 Cable Structure surfaces is just have chosen so altogether at each selected sea level elevation place, the direction of 6 outer normals is exactly " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures the direction of Cable Structure along the Temperature Distribution of wall thickness " line and Cable Structure, if Cable Structure is hollow, these, two intersection points are on Cable Structure outside surface, another on an internal surface, if Cable Structure is solid, these two intersection points are all on Cable Structure outside surface, connect these two intersection points and obtain a straight-line segment, straight-line segment is chosen three points again, this straight-line segment is divided into four sections by these three points, measure two end points of Cable Structure at these three points chosen and straight-line segment, amount to the temperature of 5 points, concrete can first hole in Cable Structure, temperature sensor is embedded in this 5 some places, especially, can not hole in support cable, support cable is only measured to the temperature of support cable surface point, in any case, the temperature recorded all is called this place " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".If have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measuring temperature in Cable Structure location according to meteorology to require to choose a position, meeting obtaining in this position actual observation record the temperature that meteorology measures the Cable Structure place environment of temperature requirement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable (as long as there was sunrise the same day, this position just should by sunlight), at the flat board (square plate that the wide 3mm of such as 30cm is thick) of this position of sound production one piece of carbon steel material (such as No. 45 carbon steels), be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, reference plate can be placed in and meet the top that meteorology temperature measures the wooden thermometer screen required, the one side of this reference plate on the sunny side, be called that sunny slope (such as, time on the Northern Hemisphere, sunny slope faces up towards south, full daytime is all by sunshine, sunny slope should have the suitable gradient to make snow can not accumulate or clear up sunny slope after snow), the sunny slope of reference plate is coarse with dark color (being conducive to accepting solar irradiation), the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the non-sunny slope of reference plate is covered with insulation material (the thick calcium carbonate insulation material of such as 5mm), Real-Time Monitoring record is obtained the temperature of the sunny slope of reference plate.

B walks, Real-Time Monitoring (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) record R the Cable Structure surface temperature measured data obtaining above-mentioned R Cable Structure surface point, Real-Time Monitoring (can be measured by ordinary temperature measuring method simultaneously, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) obtain the temperature profile data of previously defined Cable Structure along thickness, Real-Time Monitoring (can be measured by ordinary temperature measuring method simultaneously, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record the temperature record obtaining the Cable Structure place environment meeting the requirement of meteorology measurement temperature, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes the same day, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T _emax, (such as first the temperature measured data sequence of Cable Structure place environment is carried out curve fitting by Conventional mathematical calculating by the temperature measured data sequence of Cable Structure place environment, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of Cable Structure place environment, this rate of change is also along with time variations, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, such as use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, such as every 10 minutes survey records temperature data) obtain being carved at sunrise the same day measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T _pmax, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, thermal resistance is such as used to measure Cable Structure surface point, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes the same day, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T _smax, (such as first each Cable Structure surface temperature measured data sequence is carried out curve fitting by Conventional mathematical calculating by each Cable Structure surface temperature measured data sequence, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of each Cable Structure surface point, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _tmax.

C walks, and survey calculation obtains Cable Structure steady temperature data; First, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology; The a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, Δ T _pmaxwith Δ T _smaxall be not more than 5 degrees Celsius; The b condition that Section 2 must meet be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the Δ T that survey calculation obtains above _emaxbe not more than with reference to temperature difference per day Δ T _r, and the Δ T that survey calculation obtains above _pmaxdeduct 2 degrees Celsius and be not more than Δ T _emax, and the Δ T that survey calculation obtains above _smaxbe not more than Δ T _pmax; One of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition; Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise; Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _tmaxbe not more than 1 degree Celsius.This method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly moment of the acquisition Cable Structure steady temperature data of this method, then, according to Cable Structure heat transfer characteristic, utilize R Cable Structure surface temperature measured data and " HBE Cable Structure is along the thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model (such as finite element model) of Cable Structure, the Temperature Distribution that (such as finite element method) obtains the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%; Cable Structure surface comprises support cable surface; Second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtained, gets Δ T for convenience of describing _hunit be DEG C/m, be m for convenience of describing the unit getting Δ h; " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two; 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure.

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

If total N root support cable, first determines the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step.Determine the point being monitored (namely the specified point of all sign Cable Structure angle displacements, is provided with K specified point) of specifying, number to all specified points; Determined the monitored straight line (establishing each point being monitored to have L monitored straight line) of each point being monitored, to all monitored straight line numberings; Determine the monitored angle coordinate component (establishing each monitored straight line to have H angle coordinate component) of each monitored straight line, to all monitored angle coordinate component numberings.Above-mentioned numbering will be used for equally generating vector sum matrix in subsequent step." the whole monitored angle-data of Cable Structure " by K specified point in the Cable Structure determined above, cross L of each specified point appointment straight line, H angle coordinate component of each appointment straight line describe, the change of Cable Structure angle is exactly the change of all angle coordinate components of specifying of all specified points, all appointment straight line.Each total M (M=K × L × H) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.K and M must not be less than the quantity N of support cable.For simplicity, in the method " the monitored angle-data of Cable Structure " is called " monitored amount ".Only can measure an angle coordinate of a monitored straight line in each point being monitored, such as, measure the angle coordinate of Cable Structure surface normal relative to acceleration of gravity direction of specified point, be in fact exactly measurement of dip angle here.

When Cable Structure is completed, or before setting up health monitoring (damaged cable identification) system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T _orepresent, be called initial Cable Structure steady temperature data vector T _o.T is obtained in actual measurement _owhile, namely at the synchronization in the moment of the initial Cable Structure steady temperature data vector of acquisition, 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.

Specifically the synchronization in moment of so-and-so Cable Structure steady temperature data vector such as (such as initial or current) can be being obtained according to following method in this method, so-and-so method survey calculation is used to obtain the data of the monitored amount of so-and-so measured amount (all monitored amount of such as Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) while, such as every 10 minutes survey records temperature, so simultaneously equally also every 10 minutes the monitored amount of so-and-so measured amount of survey record (all monitored amount of such as Cable Structure) data.Once determine the moment obtaining Cable Structure steady temperature data, so just be called the synchronization in the moment obtaining Cable Structure steady temperature data, the data of the monitored amount of so-and-so measured amount using so-and-so method survey calculation method to obtain with the data of the monitored amount of so-and-so measured amount (all monitored amount of such as Cable Structure) of the moment synchronization of acquisition Cable Structure steady temperature data.

Conventional method (consult reference materials or survey) is used to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained at Actual measurement _owhile, namely at the initial Cable Structure steady temperature data vector T of acquisition _othe synchronization in moment, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, Cable Structure strain data, Cable Structure angle measurement data, Cable Structure volume coordinate measurement data that the Non-destructive Testing Data of support cable etc. can express the health status of rope.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 Cable 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.The data utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope set up cable system initial damage vector d _o.If there is no the Non-destructive Testing Data of rope and other are when can express the data of the health status of support cable, or when can think that Cable Structure original state is not damaged state, vectorial d _oeach element numerical value get 0.The temperature variant 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 steady temperature data vector T _o, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and 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, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T _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 cable system initial damage vector d of middle support cable _orepresent, the initial Cable Structure steady temperature data vector T of Cable Structure Cable Structure steady temperature data _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 _o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T _o"; Corresponding to A _osupport cable health status with cable system initial damage vector d _orepresent; Corresponding to A _othe initial value monitored amount initial value vector C of all monitored amount _orepresent.T _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, A _oand A ⁱ _ocount temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When i-th circulation starts, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _orepresent, vector T ⁱ _odefinition mode and vector T _odefinition mode identical, T ⁱ _oelement and T _oelement one_to_one corresponding; The current initial damage vector of cable system 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 cable system, 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; T ⁱ _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, T ⁱ _obe designated as T ¹ _o, set up T ¹ _omethod for making T ¹ _oequal T _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 installing damaged cable recognition system.Hardware components at least comprises: monitored amount monitoring system (such as containing angle measuring sensor system, signal conditioner etc.), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, computing machine and communication alert equipment.Each monitored amount, each temperature 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 cable system running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring rope and having damage, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.

5th step: establishment and on computers installation and operation temperature variation time based on the approximant recognition methods system software of damaged cable of angle monitor, the function (in this specific implementation method all work that can complete with computing machine) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method " based on the approximant recognition methods of the damaged cable of angle monitor during temperature variation " required by task and wants.

6th step: step starts circulation running thus, in Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement ⁱ, vector T ⁱdefinition mode and vector T _odefinition mode identical, T ⁱelement and T _oelement one_to_one corresponding; In actual measurement vector T ⁱwhile, namely at acquisition current cable structure steady temperature data vector T ⁱthe synchronization in moment, 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.

7th step: according to current cable structure steady temperature data vector T ⁱ, upgrade current initial mechanical Calculation Basis model A according to step a to b ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _o, and cable system current initial damage vector d ⁱ _oremain unchanged;

A. the T obtained in the 6th step is utilized ⁱ, compare T ⁱand T ⁱ _oif, T ⁱequal T ⁱ _o, then A ⁱ _o, T ⁱ _o, C ⁱ _oand d ⁱ _oremain unchanged;

B. the T obtained in the 6th step is utilized ⁱ, compare T ⁱand T ⁱ _oif, T ⁱbe not equal to T ⁱ _o, then need current initial mechanical Calculation Basis model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _oupgrade, update method is: first calculate T ⁱwith T _odifference, T ⁱwith T _odifference represent with steady temperature change vector S, S equals T ⁱdeduct T _o; At initial mechanical Calculation Basis model A _obasis on make the health status of support cable be cable system current initial damage vector d ⁱ _o, more further to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A _oin Cable Structure applies to obtain after temperature variation is exactly the current initial mechanical Calculation Basis model A of renewal ⁱ _o; At renewal A ⁱ _owhile, T ⁱ _oall elements numerical value T ⁱall elements numerical value replaces, and namely have updated T ⁱ _o, so just obtain and correctly correspond to A ⁱ _ot ⁱ _o, now d ⁱ _oremain unchanged; Upgrade A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, with these concrete numerical value substituting vector C ⁱ _othe element of middle correspondence, so just achieves vectorial C ⁱ _orenewal, and cable system current initial damage vector d ⁱ _oremain unchanged.

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 nominal unit damage 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, calculation times numerically equals the quantity of all support cables, N root support cable is had just to have N calculating, calculate each time and suppose that only having a support cable on the basis of original damage, increase unit damage again in cable system (such as gets 5%, 10%, 20% or 30% equivalent damage is unit damage), occur in calculating each time that the support cable damaged is different from during other time calculates the support cable occurring damaging, and the unit damage value that suppose there is the support cable of damage each time can be different from the unit damage value of other support cables, with " nominal unit damage vector D ⁱ _u" record the unit damage of the supposition of all ropes, vectorial D ⁱ _uelement number rule with vectorial d _othe coding rule of element identical, calculate the current value of all monitored amounts in Cable Structure each time, the current value of all monitored amount calculated each time forms one " monitored gauge calculates current value vector ", when hypothesis jth (j=1,2,3 ..., N) and root support cable is when having a unit damage, available C ⁱ _tjrepresent corresponding " monitored gauge calculates current value 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 ⁱ _tjdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tjelement and C _oelement one_to_one corresponding.

C. the vectorial C calculated each time ⁱ _tjdeduct vectorial C ⁱ _oobtain a vector, then after the unit damage value of supposition in being calculated divided by this by each element of this vector, obtain " numerical value change vector δ a C for monitored amount ⁱ _j"; N root support cable is had just to have N number of " the numerical value change vector of monitored amount ".

D. form by this N number of " numerical value change vector of monitored amount " " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange successively ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row correspond to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d _oelement number rule 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 each time calculate in suppose in cable system, to only have a rope to increase unit damage again on the basis of original damage while; when hypothesis jth (j=1,2,3; ..., N), when root support cable has a unit damage, calculate composition injury vector each time, use d ⁱ _tjrepresent this injury vector, it is C that corresponding monitored gauge calculates current value vector ⁱ _tj(see the 8th step), injury vector d ⁱ _tjelement number equal the quantity of rope, vectorial d ⁱ _tjall elements in only have the numerical value of an element to get to calculate each time in hypothesis increase the unit damage value of rope of unit damage, d ⁱ _tjthe numerical value of other element get 0, that be not numbering and the supposition of the element of 0 increase the rope of unit damage corresponding relation, be identical with the element of the same numbering of other vectors with the corresponding relation of this rope; d ⁱ _tjwith cable system initial damage vector d _oelement number rule identical, d ⁱ _tjelement and d _oelement be one-to-one relationship.By C ⁱ _tj, C ⁱ _o, Δ C ⁱ, d ⁱ _tjbring formula (23) into, obtain a linear relationship error vector e ⁱ _j, calculate a linear relationship error vector e each time ⁱ _j; e ⁱ _jsubscript j represent jth (j=1,2,3 ..., N) and root support cable has unit damage.There is N root rope just to have N calculating, just have N number of linear relationship error vector e ⁱ _j, by this N number of linear relationship error vector e ⁱ _jobtaining 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_j^i=\mathrm{abs}(\mathrm{\&Delta;}{C}^i\&CenterDot;d_{\mathrm{tj}}^i-C_{\mathrm{tj}}^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 support cable, d ⁱ _cand d ⁱelement and support cable between be one-to-one relationship, d ⁱ _cand d ⁱelement numerical value represent degree of injury or the health status of corresponding support cable, d ⁱ _cand d ⁱwith cable system 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 multi-objective optimization algorithm that can adopt has a variety of, such as: the multiple-objection optimization based on genetic algorithm, the multiple-objection optimization based on artificial neural network, the multi-objective optimization algorithm based on population, the multiple-objection optimization based on ant group algorithm, leash law (Constrain Method), weighted method (Weighted SUm Method), Objective Programming (Goal Attainment Method) etc.Because various multi-objective optimization algorithm is all conventional algorithm, can realize easily, this implementation step only provides for Objective Programming and solves current nominal fatigue vector d ⁱ _cprocess, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.

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.

$\begin{array}{cc}\min \mathrm{imize}& \mathrm{\&gamma;}\\ \mathrm{\&gamma;}\&Element;R,& d_c^i\&Element;\mathrm{\&Omega;}\end{array}---\left(24\right)$

$G\left(d_c^i\right)-\mathrm{W\&gamma;}\&le;g^i---\left(25\right)$

$G\left(d_c^i\right)=\mathrm{abs}(\mathrm{\&Delta;}{C}^i\&CenterDot;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 position and the degree of injury of damaged cable, namely achieve the health monitoring of cable system, achieve damaged cable identification.D ⁱ _j(i=1,2,3, J=1,2,3 ...., N) represent the actual damage value of jth root rope in i-th circulation, its definition is shown in formula (19), d ⁱ _jrepresent jth root rope not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, time between 0 and 100%, represent that the load-bearing capacity of corresponding proportion lost by jth root rope.

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 jth 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 ⁱ⁺¹ _oj(j=1,2,3 ..., N); The second, at Mechanics Calculation benchmark model A _obasis on, make A _oin the health status of rope be d ⁱ⁺¹ _oinstead of be d _oafter, more further to A _oin Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, and steady temperature change vector S equals T ⁱdeduct T _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 steady temperature data vector T needed for circulation ⁱ⁺¹ _oequal T ⁱ _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.

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

Claims (1)

1. during temperature variation based on the approximant recognition methods of the damaged cable of angle monitor, it is characterized in that described method comprises:

A. establish total N root support cable, first determine the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; Determine point being monitored, number to all points being monitored; Determined the monitored straight line of each point being monitored, to all monitored straight line numberings; Determine the monitored angle coordinate component of each monitored straight line, to all monitored angle coordinate component numberings; Above-mentioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored angle-data of Cable Structure " is made up of above-mentioned all monitored angle coordinate components; For simplicity, in the method by " the whole monitored angle-data of Cable Structure " referred to as " monitored amount "; The quantity of point being monitored must not be less than the quantity of support cable; The quantity sum of all monitored angle coordinate components must not be less than the quantity of support cable; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, the moment of survey record data is called the physical record data moment;

B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;

B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day _r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of describing _hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, only measure the temperature of the surface point of support cable, measure all temperature be selected a little, the temperature recorded is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate,

B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T _emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T _pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T _smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _tmax,

B3: survey calculation obtains Cable Structure steady temperature data; First, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology; The a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, reference plate maximum temperature difference Δ T _pmaxwith Cable Structure surface maximum temperature difference Δ T _smaxall be not more than 5 degrees Celsius; The b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxΔ T is not more than after deducting 2 degrees Celsius _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax; One of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition; Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise; Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _tmaxbe not more than 1 degree Celsius; This method utilizes above-mentioned six conditions, any one in following two kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, and simultaneously the second moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition; When the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data; If the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data; The amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method; This method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method; then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtained, gets Δ T for convenience of describing _hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure,

C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T _o"; Survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses; T is obtained in actual measurement _owhile, namely at the initial Cable Structure steady temperature data vector T of acquisition _othe synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Non-destructive Testing Data of the health status expressing support cable, 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 support coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C _o; The Non-destructive Testing Data that utilization can express the health status of support cable sets up cable system initial damage vector d _o, cable system initial damage vector d _oelement number equal N, d _oelement and support cable be one-to-one relationship, cable system initial damage vector d _oelement numerical value represent the degree of injury of corresponding support cable, if cable system initial damage vector d _othe numerical value of a certain element be 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 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 when thinking that Cable Structure original state is not damaged state, vectorial d _oeach element numerical value get 0;

Temperature variant physical and mechanical properties parameter, the initial Cable Structure steady temperature data vector T 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, the Non-destructive Testing Data of support cable, Cable Structure _owith all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " _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 _o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T _o"; Corresponding to A _osupport cable health status with cable system initial damage vector d _orepresent; Corresponding to A _othe initial value monitored amount initial value vector C of all monitored amount _orepresent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " ⁱ _o, monitored amount current initial value vector C ⁱ _o" current initial Cable Structure steady temperature data vector T ⁱ _o"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ⁱ _owith monitored amount current initial value vector C ⁱ _otime, the current initial mechanical Calculation Basis model A of Cable Structure ⁱ _ojust equal the initial mechanical Calculation Basis model A of Cable Structure _o, monitored amount current initial value vector C ⁱ _ojust equal monitored amount initial value vector C _o; A ⁱ _ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T ⁱ _o", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time ⁱ _otime, T ⁱ _ojust equal T _o; A ⁱ _othe initial health of support cable and A _othe health status of support cable identical, also use cable system initial damage vector d _orepresent, A in cyclic process below ⁱ _othe initial health of support cable use cable system initial damage vector d all the time _orepresent; T _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, work as T ⁱ _oand d _oa ⁱ _oparameter time, C ^t _oby A ⁱ _omechanics Calculation result composition; A in the method _o, C _o, d _oand T _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, A _oand A ⁱ _ocount temperature parameter, can the Effect on Mechanical Properties of accounting temperature change to Cable Structure; When i-th circulation starts, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _orepresent, vector T ⁱ _odefinition mode and vector T _odefinition mode identical, T ⁱ _oelement and T _oelement one_to_one corresponding; The current initial damage vector of cable system 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 cable system, 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; T ⁱ _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, T ⁱ _obe designated as T ¹ _o, set up T ¹ _omethod for making T ¹ _oequal T _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 Cable Structure military service process, the current data of Cable Structure steady temperature data is obtained, the current data composition current cable structure steady temperature data vector T of all " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement ⁱ, vector T ⁱdefinition mode and vector T _odefinition mode identical, T ⁱelement and T _oelement one_to_one corresponding; Vector T is obtained in actual measurement ⁱwhile, actual measurement obtains at acquisition current cable structure steady temperature data vector T ⁱmoment synchronization Cable Structure in the currency of all monitored amounts, 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 steady temperature data vector T ⁱ, upgrade current initial mechanical Calculation Basis model A according to step g 1 to g2 ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _o, and cable system current initial damage vector d ⁱ _oremain unchanged;

G1. T is compared ⁱand T ⁱ _oif, T ⁱequal T ⁱ _o, then A ⁱ _o, T ⁱ _o, C ⁱ _oand d ⁱ _oremain unchanged;

G2. T is compared ⁱand T ⁱ _oif, T ⁱbe not equal to T ⁱ _o, then need current initial mechanical Calculation Basis model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _owith monitored amount current initial value vector C ⁱ _oupgrade, update method is: first calculate T ⁱwith T _odifference, T ⁱwith T _odifference represent with steady temperature change vector S, S equals T ⁱdeduct T _o; At initial mechanical Calculation Basis model A _obasis on make the health status of support cable be cable system current initial damage vector d ⁱ _o, more further to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A _oin Cable Structure applies to obtain after temperature variation is exactly the current initial mechanical Calculation Basis model A of renewal ⁱ _o; At renewal A ⁱ _owhile, T ⁱ _oall elements numerical value T ⁱall elements numerical value replaces, and namely have updated T ⁱ _o, so just obtain and correctly correspond to A ⁱ _ot ⁱ _o, now d ⁱ _oremain unchanged; Upgrade A ⁱ _oafter, obtain A by Mechanics Calculation ⁱ _oin all monitored amounts, current concrete numerical value, with these concrete numerical value substituting vector C ⁱ _othe element of middle correspondence, so just achieves vectorial C ⁱ _orenewal;

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 nominal unit damage 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 of all support cables, N root support cable is had just to have N calculating, calculating each time in hypothesis cable system only has a support cable to increase unit damage again on the basis of original damage, occur in calculating each time that the support cable damaged is different from during other time calculates the support cable occurring damaging, and the unit damage value that suppose there is the support cable of damage each time can be different from the unit damage value of other support cables, with " nominal unit damage vector D ⁱ _u" record the unit damage of the supposition of all ropes, vectorial D ⁱ _uelement number rule with vectorial d _othe coding rule of element identical, calculate the current value of all monitored amounts in Cable Structure each time, the current value of all monitored amount calculated each time forms one " monitored gauge calculates current value vector "; When supposing that jth root support cable has unit damage, available C ⁱ _tjrepresent corresponding " monitored gauge calculates current value 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 ⁱ _tjdefinition mode and vectorial C _odefinition mode identical, C ⁱ _tjelement and C _oelement one_to_one corresponding;

H3. the vectorial C calculated each time ⁱ _tjdeduct vectorial C ⁱ _oobtain a vector, then after the unit damage value of supposition in being calculated divided by this by each element of this vector, obtain " numerical value change vector δ a C for monitored amount ⁱ _j"; N root support cable is had just to have N number of " the numerical value change vector of monitored amount ";

H4. form by this N number of " numerical value change vector of monitored amount " " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange successively ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row correspond to one the numerical value change of the monitored amount " vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d _oelement number rule identical;

I. current nominal fatigue vector d is defined ⁱ _cwith current actual damage vector d ⁱ, d ⁱ _cand d ⁱelement number equal the quantity of support cable, d ⁱ _cand d ⁱelement and support cable between be one-to-one relationship, d ⁱ _cand d ⁱelement numerical value represent degree of injury or the health status of corresponding support cable, d ⁱ _cand d ⁱwith cable system 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{\&Delta;}{C}^i\&CenterDot;d_c^i$ Formula 1

K. the current actual damage vector d utilizing formula 2 to express ⁱa jth element d ⁱ _jwith cable system current initial damage vector d ⁱ _oa jth element d ⁱ _ojwith current nominal fatigue vector d ⁱ _ca jth element d ⁱ _cjbetween relation, calculate current actual damage vector d ⁱall elements;

$d_j^i=1-(1-d_{\mathrm{ij}}^i)(1-d_{\mathrm{cj}}^i)$ Formula 2

J=1 in formula 2,2,3 ..., N; Current actual damage vector d ⁱelement numerical value represent the degree of injury of corresponding support cable, so according to current actual damage vector d ⁱjust can define the impaired and degree of injury of which support cable, namely achieve the health monitoring of cable system in Cable Structure, achieve damaged cable identification; If the numerical value of a certain element of current actual damage vector 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;

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 jth element;

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

Element B in formula 4 ⁱ _jmark vector B ⁱa jth element, D ⁱ _ujnominal unit damage vector D ⁱ _ua jth element, d ⁱ _cjcable system current nominal fatigue vector d ⁱ _ca jth element, they all represent the relevant information of jth root support cable, j=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, cable system current initial damage vector d namely needed for the i-th+1 time circulation ⁱ⁺¹ _oeach element;

$d_{\mathrm{oj}}^{i+1}=1-(1-d_{\mathrm{oj}}^i)(1-D_{\mathrm{uj}}^i{B}_j^i)$ Formula 5

D in formula 5 ⁱ⁺¹ _ojthe cable system current initial damage vector d next time, namely needed for the i-th+1 time circulation ⁱ⁺¹ _oa jth element, d ⁱ _ojthis, i.e. the cable system of i-th circulation current initial damage vector d ⁱ _oa jth element, D ⁱ _ujthe nominal unit damage vector D of i-th circulation ⁱ _ua jth element, B ⁱ _jthe mark vector B of i-th circulation ⁱa jth element, j=1 in formula 5,2,3 ..., N;

O. take off once, namely the i-th+1 time circulation needed for current initial Cable Structure steady temperature data vector T ⁱ⁺¹ _oequal the current initial Cable Structure steady temperature data vector T of i-th circulation ⁱ _o;

P. at initial mechanical Calculation Basis model A _obasis on, make the health status of rope be d ⁱ⁺¹ _oafter, to initial mechanical Calculation Basis model A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just get this, the steady temperature change vector S that namely uses in i-th circulation, to A _oin Cable Structure applies to obtain after temperature variation is exactly Mechanics Calculation benchmark model A next time, namely needed for the i-th+1 time circulation ⁱ⁺¹; 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;

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