CN102706646A - Approximant type identification method of damaged cable based on space coordinate monitoring during temperature variation - Google Patents

Abstract

The invention discloses an approximant type identification method of a damaged cable based on space coordinate monitoring during temperature variation. The method is based on space coordinate monitoring and comprises the following steps of: judging whether needing to update a mechanical calculation benchmark model of a cable structure by monitoring cable structure temperature, environment temperature and damage degree; updating the mechanical calculation benchmark model of the cable structure only when temperature varies and (or) the damage degree is high, thus obtaining the mechanical calculation benchmark model of the cable structure with the temperature being taken into account; obtaining a unit damage monitored quantity numerical variation matrix by calculation based on the mechanical calculation benchmark model; rapidly calculating the non-inferior solution of a current nominal damage vector by algorithms such as a multiple-object optimization algorithm, according to the approximately linear relationships between a current numeric vector of the monitored quantity and a current initial numeric vector of the monitored quantity, the unit damage monitored quantity numerical variation matrix and the current nominal damage vector to be solved; and correctly determining the position and the damage degree of the damaged cable during temperature variation according to the non-inferior solution.

Description

During temperature variation based on the approximant recognition methods of the damaged cable of space coordinate monitoring

Technical field

Structures such as cable-stayed bridge, suspension bridge, truss-frame structure have a common ground; Be exactly that they have many parts that bear tensile load; Like suspension cable, main push-towing rope, hoist cable, pull bar or the like; The common ground of this class formation is to be support unit with rope, cable or the rod member that only bears tensile load, and this method is " Cable Structure " with such structure representation for simplicity.Variation along with environment temperature; The temperature of Cable Structure also can change; When the Cable Structure temperature changed, the supporting system that this method is discerned Cable Structure based on space coordinate monitoring (referred to all ropeway carrying-ropes, reaches all rod members that only bear tensile load that play supporting role, for simplicity; This patent is called " cable system " with whole support unit unifications of this class formation; But in fact cable system not only refers to support rope, also comprises the rod member that only bears tensile load, equally for ease; Censure all ropeway carrying-ropes and the rod member that only bears tensile load that all play supporting role with " supporting rope " this noun in this method; Equally for ease, censure all ropeway carrying-ropes and the rod member that only bears tensile load that all play supporting role with " supporting rope " this noun in this method) in damaged cable (truss-frame structure just is meant the impaired rod member that only bears tensile load), belong to the engineering structure health monitoring field.

Background technology

Cable system is the key components of Cable Structure normally; Its inefficacy usually brings the inefficacy of total, and the damaged cable of discerning based on structural health monitoring technology in the cable system of Cable Structure (also referring to only bear the rod member of tensile load as previously mentioned) is a kind of method that has potentiality.After the health status of cable system changes; Can cause the variation of the measurable parameter of structure; For example the volume coordinate of Cable Structure can change; In fact the variation of volume coordinate has comprised the health status information of cable system; That is to say the health status that to utilize the structure space coordinate data to judge structure, can discern damaged cable based on space coordinate monitoring (this method is called monitored volume coordinate " monitored amount ", and the back is mentioned " monitored amount " and just is meant monitored volume coordinate).Monitored amount is except the influence that receives the cable system health status; Also can receive the influence of Cable Structure temperature variation (usually can take place); Under the condition that the Cable Structure temperature changes; If can realize damaged cable identification based on monitoring to monitored amount, the safety of Cable Structure is had significant values, also there are not a kind of disclosed, effective health monitoring systems and method to solve this problem at present.

When Cable Structure has temperature variation; In order reliable monitoring and judgement to be arranged to the health status of the cable system of Cable Structure; Must there be one can rational and effective set up the method for each monitored quantitative changeization, can provides the health evaluating of more believable cable system based on the health monitoring systems of this method foundation with the relation between the health status of all ropes in the cable system.

Summary of the invention

Technical matters: the purpose of this method is when Cable Structure has temperature variation, to the health monitoring problem of cable system in the Cable Structure, discloses a kind of based on health monitor method space coordinate monitoring, that can monitor cable system in the Cable Structure rationally and effectively.

Technical scheme: this method is made up of three parts.Be respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; When two, setting up temperature variation based on the method for required knowledge base of the approximant recognition system of the damaged cable of space coordinate monitoring and parameter, during based on the temperature variation of the volume coordinate of knowledge base (containing parameter) and actual measurement Cable Structure and actual measurement Cable Structure temperature based on the approximant recognition methods of the damaged cable of space coordinate monitoring; Three, during temperature variation based on the software and hardware part of the approximant recognition system of damaged cable of space coordinate monitoring.

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

At first confirm " the temperature survey calculating method of the Cable Structure of this method ".Because the temperature of Cable Structure possibly change; For example the temperature of the different parts of Cable Structure change along with the variation of intensity of sunshine, along with the variation of environment temperature changes; The surface of Cable Structure and temperature inside possibly be time dependent sometimes; The surface of Cable Structure possibly be different with temperature inside, and the surface of Cable Structure and temperature inside difference are time dependent, and the Mechanics Calculation of the Cable Structure when this just makes the account temperature condition is quite complicated with monitoring; For cost is measured in simplification problem, minimizing calculated amount and reduction; In order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", and is specific as follows especially:

The first step; Inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; Statistics obtains in T cloudy day 0 the highest temperature and the lowest temperature between back 30 minutes of the moment of sunrise next day at each cloudy day; Sunrise be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunrise constantly, the sunrise that can inquire about data or calculate each required day through conventional meteorology constantly, each cloudy day 0 up to sunrise next day constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature; T cloudy day arranged; The maximum temperature difference of daily temperature that T cloudy day just arranged, the maximal value of getting in the maximum temperature difference of daily temperature at T cloudy day is with reference to temperature difference per day, is designated as Δ T with reference to temperature difference per day _rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation _h, for Δ T is got in convenient narration _hUnit be ℃/m.On the surface of Cable Structure, get " R Cable Structure surface point "; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data ".When on the surface of Cable Structure, getting " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " is narrated with the condition that must satisfy that distributes in the back.From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least 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 of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is 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 and dark color, the sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day; The non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate.Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly.

Second step; Monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T _EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T _PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; Have R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise, maximal value wherein is designated as Δ T _SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes.Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " 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 "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _Tmax

In the 3rd step, measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes 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 centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front _EmaxBe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _PmaxBe not more than Δ T after deducting 2 degrees centigrade _Emax, and Cable Structure surface maximum temperature difference Δ T _SmaxBe not more than Δ T _PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains 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 per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains 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 the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the 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 that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T _hThe numerical value that obtains is for Δ T is got in convenient narration _hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have 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 in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure.

The second portion of this method: when setting up temperature variation based on the method for required knowledge base of the approximant recognition system of the damaged cable of space coordinate monitoring and parameter, during based on the temperature variation of the volume coordinate of knowledge base (containing parameter) and actual measurement Cable Structure and actual measurement Cable Structure temperature based on the approximant recognition methods of the damaged cable of space coordinate monitoring.Can carry out successively as follows, to obtain the health status assessment of cable system more accurately.

The first step: establish total N root supporting rope, at first confirm the coding rule of supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule.

" the whole monitored spatial data of structure " described by the volume coordinate specified point of K on the structure, that reach L assigned direction of each specified point, and the variation of structure space coordinate data is exactly the variation of all volume coordinate components of K specified point.(individual volume coordinate measured value of M=K * L) or calculated value characterize the structure space coordinate information to each total M.K and M must not be less than the quantity N of supporting rope.

For simplicity, in the method " the monitored spatial data of structure " abbreviated as " monitored amount ".When mentioning " so-and-so matrix of monitored amount or so-and-so vector " in the back, also can be read as " volume coordinate so-and-so matrix or so-and-so vector ".

Set up initial Mechanics Calculation benchmark model A _oThe time; In Cable Structure completion; Perhaps before setting up health monitoring (damaged cable identification) system; Calculate " Cable Structure steady temperature data " (can use conventional thermometry to measure, for example use thermal resistance to measure) according to " the temperature survey calculating method of the Cable Structure of this method " measurement, this moment " Cable Structure steady temperature data " are used vector T _oExpression is called initial Cable Structure steady temperature data vector T _oObtain T in actual measurement _oThe time, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.Calculate initial Cable Structure steady temperature data vector T in actual measurement _oThe time, use conventional method (consult reference materials or survey) to obtain the temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Calculate initial Cable Structure steady temperature data vector T in actual measurement _oThe time, use the conventional method actual measurement to calculate the actual measurement computational data of Cable Structure.The actual measurement computational data of Cable Structure comprises that Non-Destructive Testing data of supporting rope etc. can express measured datas such as the initial geometric data of data, Cable Structure of the health status of rope, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, structural strain data, structure angle measurement data, structure space measurement of coordinates data.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 on the structure, and purpose is to confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on the bridge two ends, so-called bridge type data that Here it is.The Non-Destructive Testing data of utilization supporting rope etc. can be expressed the data of the health status of rope and set up cable system initial damage vector d _oD is used in (shown in (1)) _oThe expression Cable Structure is (with initial Mechanics Calculation benchmark model A _oThe initial health of cable system expression).If when not having Non-Destructive Testing data and other of rope can express the data of health status of supporting rope, perhaps can think when the structure original state is the not damaged state vectorial d _oEach element numerical value get 0.Utilize the measured data of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, the Non-Destructive Testing data of supporting rope, temperature variant physics and the mechanical property parameters and the initial Cable Structure steady temperature data vector T of the employed various materials of Cable Structure _oUtilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial Mechanics Calculation benchmark model A _o

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

D in the formula (1) _Oj(j=1,2,3 ...., N) the initial Mechanics Calculation benchmark model A of expression _oIn the initial damage value of j root rope of cable system, d _OjBeing to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, and T representes the transposition (back with) of vector.

Obtain T in actual measurement _oThe time, use conventional method directly to measure the initial value of all monitored amounts of the Cable Structure that calculates, form monitored amount initial value vector C _o(seeing formula (2)).Requirement is obtaining A _oThe time obtain C _o, monitored amount initial value vector C _oExpression is corresponding to A _oThe concrete numerical value of " monitored amount ".Because of subject to the foregoing, the monitored amount of calculating gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial monitored amount reliably, in the narration of back, will represent this calculated value and measured value with prosign.

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

C in the formula (2) _Ok(k=1,2,3 ...., M) be k monitored amount in the Cable Structure.Vector C _oBe to be formed according to certain series arrangement by the monitored amount of M, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.

No matter which kind of method to obtain initial Mechanics Calculation benchmark model A with _o, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T _o), based on A _oThe Cable Structure computational data that calculates must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this _oSuo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of calculating gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when truly taking place near institute's analog case reliably.Model A _oThe health status of middle supporting rope is with cable system initial damage vector d _oExpression, Cable Structure Cable Structure steady temperature data are with initial Cable Structure steady temperature data vector T _oExpression.Because based on A _oThe evaluation that calculates all monitored amounts is very near the initial value (actual measurement obtains) of all monitored amounts, so also can be used in A _oThe basis on, carry out Mechanics Calculation obtains, A _oThe evaluation of each monitored amount form monitored amount initial value vector C _oWe can say T _oAnd d _oBe A _oParameter, C _oBy A _oMechanics Calculation result form.

Second step: circulation beginning.During circulation beginning each time, the current initial damage vector of the cable system d in the time of at first need setting up or set up this circulation beginning ⁱ _o(i=1,2,3 ...), set up the current initial Mechanics Calculation benchmark model A of Cable Structure ⁱ _o(finite element benchmark model for example, A in circulation each time ⁱ _oBring in constant renewal in), A ⁱ _oTemperature Distribution with " current initial Cable Structure steady temperature data vector T ⁱ _o" express.Letter i is except the place of representing number of steps significantly, and alphabetical in the method i only representes cycle index, i.e. the i time circulation.A _oAnd A ⁱ _oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure.

The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d ⁱ _oD is used in (shown in (3)) ⁱ _oCable Structure is (with current initial Mechanics Calculation benchmark model A when representing this time circulation beginning ⁱ _oThe health status of cable system expression).

$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 the formula (3) ⁱ _Oj(i=1,2,3, J=1,2,3 ...., when N) the i time circulation of expression begins, current initial Mechanics Calculation benchmark model A ⁱ _oIn the initial damage value of j root rope of cable system, d ⁱ _OjBeing to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%.

Set up and renewal d ⁱ _oMethod following:

During circulation beginning for the first time, set up the current initial damage vector of cable system and (be designated as d according to formula (3) ¹ _o) time, d ¹ _oJust equal d _oI (i=2,3,4,5,6 ...) the current initial damage vector of the cable system of needs d when inferior circulation begins ⁱ _o, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.

I (i=1,2,3,4,5,6 ...) the Mechanics Calculation benchmark model of the Mechanics Calculation benchmark model that need set up when inferior circulation begins or the Cable Structure of having set up is designated as current initial Mechanics Calculation benchmark model A ⁱ _oCorresponding to A ⁱ _o" Cable Structure steady temperature data " use vector T ⁱ _oExpression is called current initial Cable Structure steady temperature data vector T ⁱ _oVector T ⁱ _oDefinition mode and vector T _oDefinition mode identical, when beginning circulation must be set up or set up and be called current initial Cable Structure steady temperature data vector T each time ⁱ _o

Set up, upgrade A ⁱ _oAnd T ⁱ _oMethod following:

The Mechanics Calculation benchmark model of the Cable Structure of setting up during circulation beginning for the first time is designated as A ¹ _o, A ¹ _oEqual A _o, T ¹ _oEqual T _oA in circulation each time ⁱ _oAnd T ⁱ _oBring in constant renewal in, concrete grammar is the literary composition narration in the back; When loop ends each time, upgrade A ⁱ _oAnd T ⁱ _oThe Mechanics Calculation benchmark model of required Cable Structure when next time being circulated beginning, concrete grammar is the literary composition narration in the back.

This method is with " the current initial value vector of monitored amount C ⁱ _o" (i=1,2,3 ...) initial value (referring to formula (4)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation of expression begins, C ⁱ _oAlso can be called " the i time current initial value of the monitored amount of circulation 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 the formula (2) ⁱ _Ok(i=1,2,3, K=1,2,3 ...., k monitored amount when M) being the i time circulation beginning, in the Cable Structure.Vector C ⁱ _oBe that M monitored amount by front definition forms according to certain series arrangement, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.

Setting up model A ⁱ _oThe time set up " the current initial value vector of monitored amount C ⁱ _o", the current initial value vector of monitored amount C ⁱ _oExpression is corresponding to A ⁱ _oThe concrete numerical value of all monitored amounts, C ⁱ _oElement and C _oElement corresponding one by one, represent that respectively all monitored amounts are in A in Cable Structure ⁱ _oAnd A _oConcrete numerical value during two states.

Set up and renewal C ⁱ _oConcrete grammar following:

During circulation beginning for the first time, C ¹ _o(i=1, C ⁱ _oBe embodied as C ¹ _o) equal C _oI (i=2,3,4,5,6 ...) " the current initial value vector of the monitored amount C of the i time circulation of needs when inferior circulation begins ⁱ _o", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) calculate before the loop ends and obtain, concrete grammar is in the back literary composition narration.The i time (i=1,2,3,4,5,6 ...) in the circulation, " the current initial value vector of monitored amount C ⁱ _o" bring in constant renewal in, concrete grammar is the literary composition narration in the back.Because according to model A ⁱ _oThe initial value of calculating the monitored amount of gained approaches corresponding measured value reliably, in the narration of back, will represent this calculated value composition of vector and measured value composition of vector with prosign.

We can say T ⁱ _oAnd d ⁱ _oBe A ⁱ _oCharacterisitic parameter, C ⁱ _oBe A ⁱ _oAt T ⁱ _oAnd d ⁱ _oMechanics Calculation result under the condition forms.

The 3rd step: in Cable Structure military service process, in circulation each time, in other words in i (i=1,2,3,4,5,6 ...) in the inferior circulation, at known A ⁱ _o, T ⁱ _o, C ⁱ _oAnd d ⁱ _oAfter; Constantly survey the current data of calculating acquisition " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed " current cable structure steady temperature data vector T ⁱ", vector T ⁱDefinition mode and vector T _oDefinition mode identical; In the actual measurement vector T ⁱThe time, actual measurement obtains the currency of all monitored amounts in the Cable Structure, and all these numerical value are formed the current numerical value vector of monitored amount C ⁱC ⁱElement and C _oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.

Obtaining vector T ⁱAfter, upgrade A according to following concrete grammar ⁱ _o, T ⁱ _o, C ⁱ _oAnd d ⁱ _o:

Compare T ⁱAnd T ⁱ _oIf, T ⁱEqual T ⁱ _o, then need be to A ⁱ _oUpgrade; If T ⁱBe not equal to T ⁱ _o, then need be to A ⁱ _oUpgrade, at this moment T ⁱWith T _oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ⁱWith T _oThe difference represent that with steady temperature change vector S S equals T ⁱDeduct T _o, S representes the variation of Cable Structure steady temperature data.Upgrade A ⁱ _oMethod be: at A _oThe basis on make that the health status of rope is the current initial damage vector of cable system d ⁱ _o, more further to A _oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A _oIn Cable Structure to apply what obtain after the temperature variation be exactly the current initial Mechanics Calculation benchmark model A that upgrades ⁱ _o, upgrade A ⁱ _oThe time, T ⁱ _oAll elements numerical value is also used T ⁱAll elements numerical value replaces, and has promptly upgraded current initial Cable Structure steady temperature data vector T ⁱ _o, so just obtained correctly corresponding to A ⁱ _oT ⁱ _o, this moment d ⁱ _oRemain unchanged, so far just realized A ⁱ _oRenewal.When upgrading A ⁱ _oAfter, this moment A ⁱ _oThe health status of rope with the current initial damage of cable system vector d ⁱ _oExpression, A ⁱ _oThe Cable Structure steady temperature with current cable structure steady temperature data vector T ⁱExpression obtains A through Mechanics Calculation ⁱ _oIn concrete numerical value all monitored amounts, current, with these concrete numerical value replacements C ⁱ _oMiddle corresponding element has so just been realized the current initial value vector of monitored amount C ⁱ _oRenewal.

The 4th step: circulation time must be set up " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector " earlier each time, and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C ⁱ(i=1,2,3 ...)." nominal unit damage vector " that the i time circulation set up is designated as D ⁱ _uΔ C in circulation each time ⁱAnd D ⁱ _uNeed according to circumstances to bring in constant renewal in, promptly upgrading current initial Mechanics Calculation benchmark model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oAfter, upgrade unit damage monitored numerical quantity transformation matrices Δ C ⁱWith nominal unit damage vector D ⁱ _u

Earlier set up unit damage monitored numerical quantity transformation matrices Δ C during circulation beginning each time by following step ⁱWith nominal unit damage vector D ⁱ _uIf in the 3rd step, upgraded A ⁱ _o, in this step, must rebulid (promptly upgrading) unit damage monitored numerical quantity transformation matrices Δ C so ⁱWith nominal unit damage vector D ⁱ _uIf in the 3rd step, do not upgrade A ⁱ _o, in this step, needn't rebulid unit damage monitored numerical quantity transformation matrices Δ C so ⁱWith nominal unit damage vector D ⁱ _uSet up and rebulid (promptly upgrading) Δ C ⁱAnd D ⁱ _uDetailed process identical, row as follows:

Current initial Mechanics Calculation benchmark model A in Cable Structure ⁱ _oThe basis on carry out several times and calculate, equal the quantity of all ropes on the calculation times numerical value.Calculating each time in the hypothesis cable system has only a rope on the basis of original damage (original damage can be 0, can not be 0 also), to increase unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again.Calculate for convenient; When setting unit damage in the circulation each time can all be structural health conditions during this time circulation beginning as being healthy fully, and set on this basis unit damage (in subsequent step, damage numerical value that calculate, rope---be called name damage d ⁱ _c(i=1,2,3 ...), all with respect to this time when beginning circulation, with the health status of rope as being healthy fully speech, therefore must according to after the formula that provides of the literary composition name that will calculate damage be converted into true damage).The rope of appearance damage was different from the rope that appearance damages in other time calculating during a round-robin calculated each time together, and supposed that each time the unit damage value of the rope that damage is arranged can be different from the unit damage value of other ropes, with " the vectorial D of nominal unit damage ⁱ _u" (shown in (5)) write down the unit damage of the supposition of all ropes in each time circulation, circulation time is designated as D for the first time ¹ _uCalculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, the current calculated value of the M of appointment monitored amount in front; The current calculated value that calculates gained M monitored amount is each time formed one " monitored amount is calculated current numerical value vector ", and (when hypothesis j root rope had unit damage, available formula (6) represented that the monitored amount of M monitored amount of all appointments calculates current numerical value vector C ⁱ _Tj); The monitored amount that calculates is each time calculated current numerical value vector and is deducted the current initial value vector of monitored amount C ⁱ _o, the gained vector is exactly that " the numerical value change vector of monitored amount " of (is mark with the position of rope that unit damage is arranged or numbering etc.) (when j root rope has unit damage, used δ C under this condition ⁱ _jThe numerical value change vector of representing monitored amount, δ C ⁱ _jDefinition see formula (7), formula (8) and formula (9), formula (7) deducts after the formula (4) again divided by vectorial D for formula (6) ⁱ _uJ element D ⁱ _UjGained), the numerical value change of monitored amount vector δ C ⁱ _jEach element representation since when calculating supposition the unit damage (D for example of the Na Gensuo (for example j root rope) of unit damage is arranged ⁱ _Uj), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage D of supposition ⁱ _UjRate of change; There is N root rope that N " the numerical value change vector of monitored amount " just arranged; The numerical value change vector of each monitored amount has M element, forms " the unit damage monitored numerical quantity transformation matrices Δ C that M * N element arranged successively by this N " the numerical value change vector of monitored amount " ⁱ" (the capable N row of M), each vectorial δ C ⁱ _j(j=1,2,3 ...., N) be matrix Δ C ⁱOne row, Δ C ⁱDefinition suc as formula shown in (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 the formula (5) ⁱ _uElement D ⁱ _Uj(i=1,2,3, J=1,2,3 ...., N) the unit damage numerical value of the j root rope of supposition in the i time circulation of expression, 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 the formula (6) ⁱ _Tjk(i=1,2,3, J=1,2,3 ...., N; K=1,2,3 ...., M) the i time circulation of expression be owing to j root rope when unit damage is arranged, according to the current numerical value of calculating of the monitored amount of pairing k the appointment of coding rule.

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

The subscript i of each amount in the formula (7) (i=1,2,3 ...) the i time circulation of expression, subscript j (j=1,2,3 ...., N) expression j root rope has unit damage, D in the formula ⁱ _UjBe vectorial D ⁱ _uIn j element.Vector δ C ⁱ _jDefinition suc as formula shown in (8), δ C ⁱ _jK (k=1,2,3 ...., M) individual element δ C ⁱ _JkRepresent to set up matrix Δ C in the i time circulation ⁱThe time, suppose that the change amount of calculating a gained k monitored amount when j root rope has unit damage is with respect to the unit damage D that supposes ⁱ _UjRate of change, it defines suc 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)$

The definition of each amount was narrated in front in the formula (9).

${\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 the formula (10) ⁱ _j(i=1,2,3 ....,, j=1,2,3 ...., N) in the i time circulation of expression, because j root rope has unit damage D ⁱ _UjCause, the relative value of all monitored amounts changes.Matrix Δ C ⁱCoding rule and the front vector d of row (subscript j) ⁱ _oThe coding rule of subscript j of element identical.

The 5th step: the current health status of identification cable system.Detailed process is following.

I (i=1,2,3 ...) in the inferior circulation, be utilized in " the current numerical value vector of the monitored amount C that the actual measurement of the 3rd step obtains ⁱ" " the current initial value of monitored amount vector C together ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " the vectorial d of current name damage ⁱ _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)$

The current numerical value vector of monitored amount C in formula (11) and the formula (12) ⁱDefinition be similar to the current initial value of monitored amount vector C ⁱ _oDefinition, see formula (13); The vectorial d of the current name damage of cable system ⁱ _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 the formula (13) ⁱ _k(i=1,2,3 ....; K=1,2,3 ...., M) be the i time circulation time Cable Structure, according to the current numerical value of the monitored amount of the pairing k of being numbered of 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 the formula (14) ⁱ _Cj(i=1,2,3 ....; J=1,2,3 ...., N) be the current nominal impairment value of cable system j root rope in the i time circulation, vectorial d ⁱ _cThe coding rule and the matrix Δ C of subscript j of element ⁱThe coding rule of row identical.

When the rope actual damage was not too big, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure was also less, and the represented a kind of like this linear relationship of formula (11) or formula (12) is less with the error of actual conditions, and error can be used error vector e ⁱ(formula (15)) definition, the error of linear relationship shown in expression (11) or the formula (12).

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

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

Because there are certain error in formula (11) or the represented linear relationship of formula (12), therefore can not be simply according to formula (11) or formula (12) and " the current numerical value vector of monitored amount C ⁱ" come directly to find the solution to obtain the vectorial d of Suo Dangqian name damage ⁱ _cIf done like this, the damage vector d that obtains ⁱ _cIn element in addition bigger negative value can appear, just negative damage, this obviously is irrational.Therefore obtaining rope damages vectorial d ⁱ _cAcceptable separating (promptly have reasonable error, but can be more accurately from cable system, confirm the position and the degree of injury thereof of damaged cable) become a rational solution, available formula (16) is expressed this method.

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

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

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

G in the formula (17) ⁱ _k(i=1,2,3 ....; K=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (11) or the formula (12) in the i time circulation has been described.Vector g ⁱCan be according to the error vector e of formula (15) definition ⁱTentative calculation is selected.

At the current initial value vector of monitored amount C ⁱ _o, unit damage monitored numerical quantity transformation matrices Δ C ⁱWith the current numerical value vector of monitored amount C ⁱWhen known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (16), obtain the vectorial d of the current name damage of cable system ⁱ _cAcceptable separating, the current actual damage of cable system vector d ⁱThe element of (formula (18) is seen in definition) can calculate according to formula (19), thereby can be by d ⁱConfirm the position and the degree of injury of damaged cable, just realized the health monitoring of cable system, realized 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 the formula (18) ⁱ _j(i=1,2,3, J=1,2,3 ...., N) the actual damage value of j root rope in the i time circulation of expression, formula (19), d are seen in its definition ⁱ _jBeing to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, vectorial d ⁱCoding rule and the formula (1) of element 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 the formula (19) ⁱ _Oj(i=1,2,3,4, J=1,2,3 ...., N) be the current initial damage vector of cable system d ⁱ _oJ element, d ⁱ _CjBe the vectorial d of the current name damage of cable system ⁱ _cJ element.

The 6th step: judge whether to finish this (the i time) circulation, if, then accomplish the preceding tailing in work of this loop ends, for next time (promptly the i+1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following:

In this (the i time) circulation, try to achieve the vectorial d of current name damage ⁱ _cAfter, at first, set up mark vector B according to formula (20) ⁱ, formula (21) has provided mark vector B ⁱThe definition of j element; If mark vector B ⁱElement be 0 entirely, then got back to for the 3rd step and proceed health monitoring and calculating cable system; If mark vector B ⁱElement be not 0 entirely, then accomplish subsequent step after, get into circulation next time.

So-called subsequent step is: at first, according to formula (22) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation ^I+1 _oEach element d ^I+1 _OjThe second, at Mechanics Calculation benchmark model A _oThe basis on, make A _oIn the health status of rope be d ^I+1 _oRather than be d _oAfter, more further to A _oIn Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation that applies just taken from steady temperature change vector S, and steady temperature change vector S equals T ⁱDeduct T _o), so just obtained next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Mechanics Calculation benchmark mould A of circulation ^I+1 _o, next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation ^I+1 _oEqual T ⁱ _o, to A ^I+1 _oCarrying out Mechanics Calculation obtains corresponding to A ^I+1 _oConcrete numerical value all monitored amounts, current, these concrete numerical value form next time (promptly the i+1 time, i=1,2,3,4 ...) the current initial value vector C of the required monitored amount of circulation ^I+1 _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 the formula (20) ⁱSubscript i represent the i time the circulation, its element B ⁱ _j(j=1,2,3 ..., subscript j N) representes the damage characteristic of j root rope, can only get 0 and 1 two amount, concrete value rule is seen 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 the formula (21) ⁱ _jBe mark vector B ⁱJ element, D ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element (seeing formula (3)), d ⁱ _CjBe the vectorial d of the current name damage of cable system ⁱ _cJ element (seeing formula (14)), they all represent the relevant information of j 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 the formula (22) ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element (seeing formula (5)), d ⁱ _OjBe the current initial damage vector of cable system d ⁱ _oJ element (seeing formula (3)).

The third part of this method: during temperature variation based on the software and hardware part of the approximant recognition system of damaged cable of space coordinate monitoring.

Hardware components comprises monitoring system (comprising monitored amount monitoring system, temperature monitoring system), signal picker and computing machine etc.Require the temperature required measured data of monitoring acquisition in real time, require simultaneously each monitored amount of monitoring in real time.

Software section should be accomplished the process that this method sets, and promptly accomplishes functions such as needed in this method, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.

This method specifically comprises:

A. establish total N root supporting rope, at first confirm the coding rule of supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; Confirm the monitored point with monitored volume coordinate of appointment, give all monitored some numberings; Confirmed each monitored point with monitored volume coordinate component, give all monitored volume coordinate components numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored spatial data of Cable Structure " is made up of above-mentioned all monitored volume coordinate components; For simplicity, in the method " the monitored spatial data of Cable Structure " is called " monitored amount "; The quantity of monitored point must not be less than the quantity of supporting rope; The quantity sum of all monitored volume coordinate components must not be less than the quantity of supporting rope; Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly;

B. this method definition " the temperature survey calculating method of the Cable Structure of this method " set by step b1 to b3 carry out;

B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure; Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; One be called the cloudy day all day with what can not see the sun daytime in the method; Statistics obtain in T cloudy day each cloudy day 0 up to the sunrise next day highest temperature and the lowest temperature between back 30 minutes constantly, sunrise is meant the sunrise moment on the meteorology that base area revolutions and revolution rule confirm constantly, does not represent necessarily can see the sun same day; The sunrise that can inquire about data or calculate each required day through conventional meteorology constantly; Each cloudy day 0 up to next day sunrise constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, T cloudy day arranged, the maximum temperature difference of T cloudy daily temperature is just arranged; Get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day _rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation _h, for Δ T is got in convenient narration _hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data "; From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least 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 of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness "; Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is called reference plate; Reference plate can not contact with ground; Reference plate overhead distance is not less than 1.5 meters, and the one side of this reference plate is called sunny slope on the sunny side; The sunny slope of reference plate is coarse and dark color; The sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day, the non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate;

B2: monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains Cable Structure place environment; Be called the environment maximum temperature difference, be designated as Δ T _EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day of temperature that minimum temperature obtains the sunny slope of reference plate; Be called the reference plate maximum temperature difference, be designated as Δ T _PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; There is R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise; Maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T _SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes; Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " 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 "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _Tmax

B3: measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes 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 centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front _EmaxBe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _PmaxBe not more than Δ T after deducting 2 degrees centigrade _Emax, and Cable Structure surface maximum temperature difference Δ T _SmaxBe not more than Δ T _PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains 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 per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains 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 the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the 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 that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T _hThe numerical value that obtains is for Δ T is got in convenient narration _hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have 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 in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure;

C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " and calculate the Cable Structure steady temperature data under the original state; Cable Structure steady temperature data under the 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 physics and the mechanical property parameters of the employed various materials of Cable Structure; Obtain T in actual measurement _oThe time, just obtaining initial Cable Structure steady temperature data vector T _oThe synchronization in the moment; Directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure comprises the Non-Destructive Testing data of the health status of expressing the supporting rope, the initial value of all monitored amounts, the initial rope force data of all supporting ropes, 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 is formed monitored amount initial value vector C _oUtilize the Non-Destructive Testing data of the health status that can express the supporting rope to set up cable system initial damage vector d _o, cable system initial damage vector d _oElement number equal N, d _oElement be one-to-one relationship with the supporting rope, cable system initial damage vector d _oElement numerical value represent the degree of injury of corresponding supporting rope, if cable system initial damage vector d _oThe numerical value of a certain element be 0, represent that the pairing supporting rope of this element is intact, not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%; Represent that then this supporting rope lost the load-bearing capacity of corresponding proportion; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, perhaps think when the Cable Structure original state is the not damaged state vectorial d _oEach element numerical value get 0;

D. according to the temperature variant physics of the Non-Destructive Testing data of the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, supporting rope, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure steady temperature data vector T _oWith all Cable Structure data that preceding step obtains, set up the initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " _o, based on A _oThe Cable Structure computational data that calculates must be very near 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 _oSupporting rope health status with cable system initial damage vector d _oExpression; Corresponding to A _oThe initial value of all monitored amounts with monitored amount initial value vector C _oExpression; Set up the current initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " for the first time ⁱ _o, the current initial value of monitored amount vector C ⁱ _o" current initial Cable Structure steady temperature data vector T ⁱ _o"; Set up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oThe time, the current initial Mechanics Calculation benchmark model A of Cable Structure ⁱ _oJust equal the initial Mechanics Calculation benchmark model A of Cable Structure _o, the current initial value vector of monitored amount C ⁱ _oJust equal monitored amount initial value vector C _oA ⁱ _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 Mechanics Calculation benchmark model A of Cable Structure the first time ⁱ _oThe time, T ⁱ _oJust equal T _oA ⁱ _oThe initial health and the A of supporting rope _oThe health status of supporting rope identical, also use cable system initial damage vector d _oExpression, A in the cyclic process of back ⁱ _oThe initial health of supporting rope use cable system initial damage vector d all the time _oExpression; T _oAnd d _oBe A _oParameter, by A _oThe initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result _oThe initial value of all monitored amounts of expression is identical, therefore also we can say C _oBy A _oMechanics Calculation result form, work as T ⁱ _oAnd d _oBe A ⁱ _oParameter the time, C ^t _oBy A ⁱ _oMechanics Calculation result form; A in the method _o, C _o, d _oAnd T _oBe constant;

E. in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning ⁱ _o, A _oAnd A ⁱ _oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _oExpression, vector T ⁱ _oDefinition mode and vector T _oDefinition mode identical, T ⁱ _oElement and T _oElement corresponding one by one; The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d ⁱ _o, d ⁱ _oCable Structure A when representing this time circulation beginning ⁱ _oThe health status of cable system, d ⁱ _oDefinition mode and d _oDefinition mode identical, d ⁱ _oElement and d _oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C ⁱ _oExpression, vectorial C ⁱ _oDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _oElement and C _oElement corresponding one by one, the current initial value vector of monitored amount C ⁱ _oExpression is corresponding to A ⁱ _oThe concrete numerical value of all monitored amounts; T ⁱ _oAnd d ⁱ _oBe A ⁱ _oCharacterisitic parameter; C ⁱ _oBy A ⁱ _oMechanics Calculation result form; During circulation beginning for the first time, A ⁱ _oBe designated as A ¹ _o, set up A ¹ _oMethod for making A ¹ _oEqual A _oDuring circulation beginning for the first time, T ⁱ _oBe designated as T ¹ _o, set up T ¹ _oMethod for making T ¹ _oEqual T _oDuring circulation beginning for the first time, d ⁱ _oBe designated as d ¹ _o, set up d ¹ _oMethod for making d ¹ _oEqual d _oDuring circulation beginning for the first time, C ⁱ _oBe designated as C ¹ _o, set up C ¹ _oMethod for making C ¹ _oEqual C _o

F. go on foot the q circulation in step from getting into here by f; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T ⁱ, vector T ⁱDefinition mode and vector T _oDefinition mode identical, T ⁱElement and T _oElement corresponding one by one; Obtain vector T in actual measurement ⁱThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T ⁱThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value are formed the current numerical value vector of monitored amount C ⁱ, vectorial C ⁱDefinition mode and vectorial C _oDefinition mode identical, C ⁱElement and C _oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly;

G. according to current cable structure steady temperature data vector T ⁱ, upgrade current initial Mechanics Calculation benchmark model A according to step g 1 to g2 ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _o, and the current initial damage vector of cable system d ⁱ _oRemain unchanged;

G1. compare T ⁱAnd T ⁱ _oIf, T ⁱEqual T ⁱ _o, A then ⁱ _o, T ⁱ _o, C ⁱ _oAnd d ⁱ _oRemain unchanged;

G2. compare T ⁱAnd T ⁱ _oIf, T ⁱBe not equal to T ⁱ _o, then need be to current initial Mechanics Calculation benchmark model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oUpgrade, update method is: at first calculate T ⁱWith T _oPoor, T ⁱWith T _oThe difference represent that with steady temperature change vector S S equals T ⁱDeduct T _oAt initial Mechanics Calculation benchmark model A _oThe basis on the health status of order supporting rope be the current initial damage vector of cable system d ⁱ _o, more further to A _oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A _oIn Cable Structure to apply what obtain after the temperature variation be exactly the current initial Mechanics Calculation benchmark model A that upgrades ⁱ _oUpgrading A ⁱ _oThe time, T ⁱ _oAll elements numerical value is used T ⁱAll elements numerical value replaces, and has promptly upgraded T ⁱ _o, so just obtained correctly corresponding to A ⁱ _oT ⁱ _o, this moment d ⁱ _oRemain unchanged; Upgrade A ⁱ _oAfter, obtain A through Mechanics Calculation ⁱ _oIn concrete numerical value all monitored amounts, current, replace vectorial C with these concrete numerical value ⁱ _oMiddle corresponding element has so just been realized vectorial C ⁱ _oRenewal;

H. at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity transformation matrices Δ C through calculating ⁱWith nominal unit damage vector D ⁱ _u

H1. the i time when beginning circulation, directly h2 obtains Δ C to the listed method of step h4 set by step ⁱAnd D ⁱ _uAt other constantly, when in step g to A ⁱ _oAfter upgrading, h2 regains Δ C to the listed method of step h4 set by step ⁱAnd D ⁱ _uIf, in step g not to A ⁱ _oUpgrade, then directly change step I herein over to and carry out follow-up work;

H2. at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on carry out the several times Mechanics Calculation; Equal the quantity of all supporting ropes on the calculation times numerical value; There is N root supporting rope that N calculating is just arranged; Calculating each time in the hypothesis cable system has only a supporting rope on the basis of original damage, to increase unit damage again; The supporting rope of appearance damage is different from the supporting rope that occurs damage in other time calculating in calculating each time, and supposes that each time the unit damage value of the supporting rope that damage is arranged can be different from the unit damage value of other supporting ropes, uses " nominal unit damage vector D ⁱ _u" write down the unit damage of the supposition of all ropes, vectorial D ⁱ _uElement coding rule and vectorial d _oThe coding rule of element identical, calculate the current numerical value of all monitored amounts in the Cable Structure each time, the current numerical value of the monitored amount of all that calculate is each time formed one " monitored amount is calculated current numerical value vector "; When hypothesis j root supporting rope has unit damage, available C ⁱ _Tj" monitored amount is calculated current numerical value vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C ⁱ _TjDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _TjElement and C _oElement corresponding one by one;

H3. the vectorial C that calculates each time ⁱ _TjDeduct vectorial C ⁱ _oObtain a vector, obtain " numerical value change vector δ a C of monitored amount during each element that again should vector all calculates divided by this after the unit damage value of supposition ⁱ _j"; There is N root supporting rope that N " the numerical value change vector of monitored amount " just arranged;

H4. form " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount " ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row corresponding 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 _oThe element coding rule identical;

I. define the vectorial d of current name damage ⁱ _cWith current actual damage vector d ⁱ, d ⁱ _cAnd d ⁱElement number equal to support the quantity of rope, d ⁱ _cAnd d ⁱElement and supporting be one-to-one relationship between the rope, d ⁱ _cAnd d ⁱElement numerical value represent the degree of injury or the health status of corresponding supporting rope, d ⁱ _cAnd d ⁱWith cable system initial damage vector d _oThe element coding rule identical, d ⁱ _cElement, d ⁱElement and d _oElement be one-to-one relationship;

J. according to the current numerical value vector of monitored amount C ⁱWith " the current initial value vector of monitored amount C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " the vectorial d of current name damage ⁱ _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 ⁱ _cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage ⁱ _c

$C^i=C_o^i+\mathrm{\&Delta;}{C}^i\&CenterDot;d_c^i$ Formula 1

K. the current actual damage vector d that utilizes formula 2 to express ⁱJ element d ⁱ _jWith the current initial damage vector of cable system d ⁱ _oJ element d ⁱ _OjWith the vectorial d of current name damage ⁱ _cJ 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 the formula 2,2,3 ..., N; Current actual damage vector d ⁱElement numerical value represent the degree of injury of corresponding supporting rope, so according to current actual damage vector d ⁱJust can define the impaired and degree of injury of which supporting rope, promptly realize the health monitoring of cable system in the Cable Structure, realize damaged cable identification; If the numerical value of a certain element of current actual damage vector is 0, represent that the pairing supporting rope of this element is intact, do not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity; If its numerical value between 0 and 100%, is then represented this supporting rope and has been lost the load-bearing capacity of corresponding proportion;

L. try to achieve the vectorial d of current name damage ⁱ _cAfter, set up mark vector B according to formula 3 ⁱ, formula 4 has provided mark vector B ⁱThe definition of j element;

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

formula 4

Element B in the formula 4 ⁱ _jBe mark vector B ⁱJ element, D ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element, d ⁱ _CjBe the vectorial d of the current name damage of cable system ⁱ _cJ element, they all represent the relevant information of j root supporting rope, j=1 in the formula 4,2,3 ..., N;

If mark vector B m. ⁱElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B ⁱElement be not 0 entirely, then get into next step, be step n;

N. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required cable system of circulation d ^I+1 _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 the formula 5 ^I+1 _OjBe the current initial damage vector of the required cable system d that next time, promptly circulates for the i+1 time ^I+1 _oJ element, d ⁱ _OjBe this, i.e. the i time current initial damage vector of round-robin cable system d ⁱ _oJ element, D ⁱ _UjBe the i time round-robin name unit damage vector D ⁱ _uJ element, B ⁱ _jBe the i time round-robin mark vector B ⁱJ element, j=1 in the formula 5,2,3 ..., N;

O. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation ^I+1 _oEqual the current initial Cable Structure steady temperature data vector T of round-robin the i time ⁱ _o

P. at initial Mechanics Calculation benchmark model A _oThe basis on, the health status that makes rope is d ^I+1 _oAfter, to initial Mechanics Calculation benchmark model A _oIn Cable Structure apply temperature variation, the steady temperature change vector S that the numerical value of the temperature variation that applies is just got this, be i.e. used in the i time circulation is to A _oIn Cable Structure apply obtain after the temperature variation be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation ^I+1Obtain A ^I+1After, obtain A through Mechanics Calculation ^I+1In concrete numerical value all monitored amounts, current, these concrete numerical value are formed next time, the vectorial C of the current initial value of required monitored amount that promptly circulates for the i+1 time ^I+1 _o

Q. get back to step f, beginning is circulation next time.

Beneficial effect: when receiving the influencing of factors such as sunshine and environment temperature when the temperature field of Cable Structure; The temperature field of Cable Structure is constantly to change; The change of temperature field of Cable Structure must influence the monitored amount of Cable Structure; Have only and partly rejected by the influence in temperature field monitored amount could to carry out rational cable structure health monitoring based on monitored amount; And the temperature field measurement of Cable Structure is very complicated with calculating; This method discloses and has comprised a kind of simple, economic, feasible, cable structure health monitoring method of Cable Structure temperature field computing method efficiently that is suitable for cable structure health monitoring, adopts this method synchronous when impaired at many ropes of Cable Structure, and the temperature of Cable Structure is when changing along with the time; Monitoring and evaluation identifies the health status (position and the degree of injury that comprise all damaged cables) of cable system very exactly, and the disclosed system and method for this method is very useful to effective health monitoring of cable system.

Embodiment

When temperature variation was arranged, to the health monitoring of the cable system of Cable Structure, this method disclosed a kind of system and method for the health status of each root rope in the cable system of can monitoring rationally and effectively in the identification Cable Structure.The following explanation of the embodiment of this method in fact only is exemplary, and purpose never is to limit the application or the use of this method.

This method adopts a kind of algorithm, and this algorithm is used for monitoring the health status of the cable system of discerning Cable Structure.During practical implementation, the following step is a kind of in the various steps that can take.

The first step: confirm " the temperature survey calculating method of the Cable Structure of this method " that these method concrete steps are following:

The a step: inquiry or actual measurement (can use conventional thermometry to measure; For example use thermal resistance to measure) obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure (for example finite element model).Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; Statistics obtains in T cloudy day 0 the highest temperature and the lowest temperature between back 30 minutes of the moment of sunrise next day at each cloudy day; Sunrise be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunrise constantly, the sunrise that can inquire about data or calculate each required day through conventional meteorology constantly, each cloudy day 0 up to sunrise next day constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature; T cloudy day arranged; The maximum temperature difference of daily temperature that T cloudy day just arranged, the maximal value of getting in the maximum temperature difference of daily temperature at T cloudy day is with reference to temperature difference per day, is designated as Δ T with reference to temperature difference per day _rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation _h, for Δ T is got in convenient narration _hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual observation record; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data ".From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level,, can choose height above sea level 0m, 50m, 100m and height above sea level 200m so if for example the sea level elevation of Cable Structure is between 0m to 200m; Intersect with imaginary surface level and Cable Structure surface at each sea level elevation place that chooses; Obtain intersection, surface level and Cable Structure are crossing to obtain friendship face, and intersection is the outer edge line of friendship face; Intersection place on surface level and Cable Structure surface chooses 6 points; From the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", and it is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure.In the measurement Cable Structure of choosing along in 6 directions of the Temperature Distribution of wall thickness; At first according to the sunny slope of definite Cable Structure such as the physical dimension of the meteorological data throughout the year in Cable Structure position zone and Cable Structure, volume coordinate, Cable Structure surrounding environment 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 each sea level elevation place that chooses; Aforementioned intersection respectively has one section in sunny slope and in the shade; This of intersection respectively has a mid point for two sections, crosses these two mid points and gets the outer normal of Cable Structure, and this method is called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure with these two outer normals; This method is called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure with these two outer normal directions; The outer normal of obvious sunny slope and in the shade outer normal all intersect with aforementioned intersection, and two intersection points are also just arranged, and these two intersection points are divided into two line segments with intersection; On two line segments, get 2 points respectively; Totally 4 points, each line segment is divided into 3 sections of equal in length in two line segments of the intersection of naming a person for a particular job of getting, and gets the outer normal on Cable Structure surface at these 4 some places; Just chosen the outer normal on 6 Cable Structure surfaces 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, one at these two intersection points are on the Cable Structure outside surface; Another is on inside surface, if Cable Structure is solid, these two intersection points are all on the Cable Structure outside surface; Connect these two intersection points and obtain a straight-line segment, on straight-line segment, choose three points again, these three these straight-line segments of naming a person for a particular job are divided into four sections; Three points measuring that Cable Structure chooses at this and two end points of straight-line segment, the temperature of 5 points altogether; Concrete can hole on Cable Structure earlier, and temperature sensor is embedded in this 5 some places, special; Can not on the supporting rope, hole; The supporting rope is only measured the temperature that supports the rope surface point, and in any case, the temperature that records all is called this place " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " " Cable Structure is along temperature profile data of thickness " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measurement acquisition are called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method.If chosen H different altitude above sea level; At each sea level elevation place, chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness, measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness; Wherein H and E are not less than 3; B is not less than 2, and is special, and E equals 1 for the supporting rope; That " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE; The back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement, claim to survey the temperature data that obtains and will be " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure; Obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual observation record; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground getable this day fullest sunshine (as long as the same day sunrise arranged, this position just should by solar radiation to), at the flat board (the square flat board that for example the wide 3mm of 30cm is thick) of carbon steel material of this position of sound production (for example No. 45 carbon steels); Be called reference plate; Reference plate can not contact with ground, and reference plate overhead distance is not less than 1.5 meters, and reference plate can place the top of the wooden thermometer screen that meets meteorology temperature measurement requirement; The one side of this reference plate on the sunny side; (for example, in the time of on the Northern Hemisphere, sunny slope faces up towards south to be called sunny slope; Full daytime is all by sunshine; Sunny slope should have the suitable gradient to make snow not accumulate or behind snow, clear up sunny slope), the sunny slope of reference plate is coarse and (helping accepting solar irradiation) dark color, the sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day; The non-sunny slope of reference plate is covered with insulation material (the for example thick lime carbonate insulation material of 5mm), real-time monitoring record is obtained the temperature of the sunny slope of reference plate.

The b step; Monitoring in real time (can use conventional thermometry to measure; For example use thermal resistance to measure; For example every at a distance from temperature data of 10 minutes survey records) write down R the Cable Structure surface temperature measured data that obtains above-mentioned R Cable Structure surface point, monitoring in real time simultaneously (can be used conventional thermometry to measure, for example use the thermal resistance measurement; For example every at a distance from temperature data of 10 minutes survey records) obtain Cable Structure that the front defines temperature profile data along thickness; Monitoring in real time simultaneously (can use conventional thermometry to measure, for example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature, be for example every at a distance from temperature data of 10 minutes survey records) record obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; (can use conventional thermometry to measure through real-time monitoring; For example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature; For example every at a distance from temperature data of 10 minutes survey records) record obtain being carved at sunrise the same day sunrise next day constantly the Cable Structure between back 30 minutes belong to the temperature measured data sequence of environment; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T _Emax(the temperature measured data sequence that for example earlier Cable Structure is belonged to environment carries out curve fitting through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment; Then through ask curve to the derivative of time or through ask with numerical method on the curve each corresponding to the point of survey record data time to the change of time rate) temperature that obtains Cable Structure place environment is about the change of time rate, this rate of change is also along with the time changes; (can use conventional thermometry to measure through real-time monitoring; For example use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark; For example every at a distance from temperature data of 10 minutes survey records) obtain being carved at sunrise the same day measured data sequence of temperature of the sunny slope of the reference plate between back 30 minutes of the moment of sunrise next day; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T _Pmax(can use conventional thermometry to measure through real-time monitoring; For example use thermal resistance to measure the Cable Structure surface point; For example every at a distance from temperature data of 10 minutes survey records) record obtains being carved at sunrise the same day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes of the moment of sunrise next day; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; Have R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise, maximal value wherein is designated as Δ T _Smax(for example earlier each Cable Structure surface temperature measured data sequence is carried out curve fitting through conventional mathematical computations by each Cable Structure surface temperature measured data sequence; Then through ask curve to the derivative of time or through ask with numerical method on the curve each corresponding to the point of survey record data time to the change of time rate) temperature that obtains each Cable Structure surface point is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes.Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " 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 "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _Tmax

In the c step, measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, Δ T _PmaxWith Δ T _SmaxAll be not more than 5 degrees centigrade; Second b condition that must satisfy be the same day be carved at sunrise sunrise next day constantly between back 30 minutes during this period of time in, measure the Δ T that calculates in front _EmaxBe not more than with reference to temperature difference per day Δ T _r, and measure the Δ T that calculates in front _PmaxDeduct 2 degrees centigrade and be not more than Δ T _Emax, and measure the Δ T that calculates in front _SmaxBe not more than Δ T _PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains 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 per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains 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 the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _TmaxBe not more than 1 degree centigrade.This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly the moment of the acquisition Cable Structure steady temperature data of this method constantly; Then; According to the Cable Structure heat transfer characteristic; Utilize to obtain R Cable Structure surface temperature measured data and " HBE Cable Structure is along the thickness temperature measured data " in the moment of Cable Structure steady temperature data; Utilize the thermal conduction study computation model (for example finite element model) of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer (for example finite element method) in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; The accounting temperature that also comprises HBE " measuring the point of Cable Structure " that Cable Structure is selected in front along the temperature profile data of thickness; The accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data "; When R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data in the Cable Structure in moment of acquisition Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method; This moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T _hThe numerical value that obtains is for Δ T is got in convenient narration _hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have 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 in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure.

Second step: set up initial Mechanics Calculation benchmark model A _o

If total N root supporting rope, the coding rule of at first definite rope, with rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule.Confirm the measured point (promptly all characterize the specified point of Cable Structure volume coordinate, are provided with K specified point) of appointment, give all specified point numberings; Confirm the measured volume coordinate component (established each measurement point L measured volume coordinate component arranged) of each measurement point, gave the measured volume coordinate component numbering of all appointments.Above-mentioned numbering will be used to generate the vector sum matrix equally in subsequent step." the whole monitored spatial data of Cable Structure " described by L volume coordinate component K specified point, that cross each specified point on top definite Cable Structure, and the variation of Cable Structure volume coordinate is exactly variations all specified points, all space specified coordinate components.(individual volume coordinate component measurement value of M=K * L) or calculated value characterize the volume coordinate information of Cable Structure to each total M.K and M must not be less than the quantity N of supporting rope.For simplicity, in the method " the monitored spatial data of Cable Structure " abbreviated as " monitored amount ".

In Cable Structure completion; Perhaps before setting up health monitoring (damaged cable identification) system; Calculating " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " measurement (can use conventional thermometry to measure; For example use thermal resistance to measure), this moment " Cable Structure steady temperature data " are used vector T _oExpression is called initial Cable Structure steady temperature data vector T _oObtain T in actual measurement _oThe time, just at the synchronization in the moment that obtains initial Cable Structure steady temperature data vector, use conventional method directly to measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial value vector C _o

Can be specifically in this method according to the synchronization of following method in the moment that obtains so-and-so Cable Structure steady temperature data vector such as (for example initial or current); Use so-and-so method measurement to calculate the data of the monitored amount of so-and-so measured amount (for example all monitored amounts of Cable Structure): in survey record temperature (comprising that Cable Structure belongs to the temperature and the Cable Structure surface temperature of the sunny slope of the temperature of environment, reference plate); For example every at a distance from temperature of 10 minutes survey records, equally also every so simultaneously at a distance from 10 minutes the monitored amount of so-and-so measured amount of survey record (for example all monitored amounts of Cable Structure) data.In case confirmed to obtain the moment of Cable Structure steady temperature data; Data with the monitored amount of so-and-so measured amount (for example all monitored amounts of Cable Structure) of the moment synchronization that obtains Cable Structure steady temperature data just are called the synchronization in the moment that obtains Cable Structure steady temperature data so, use so-and-so method to measure the data of the monitored amount of so-and-so measured amount that computing method obtain.

Use conventional method (consult reference materials or survey) to obtain the temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Calculate initial Cable Structure steady temperature data vector T in actual measurement _oThe time, just obtaining initial Cable Structure steady temperature data vector T _oThe synchronization in the moment, use the conventional method actual measurement to calculate the actual measurement computational data of Cable Structure.The actual measurement computational data of Cable Structure comprises that Non-Destructive Testing data of supporting rope etc. can express measured datas such as the initial geometric data of data, Cable Structure of the health status of rope, 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.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 on the Cable Structure, and purpose is to confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on the bridge two ends, so-called bridge type data that Here it is.The Non-Destructive Testing data of utilization supporting rope etc. can be expressed the data of the health status of rope and set up cable system initial damage vector d _oIf when not having Non-Destructive Testing data and other of rope can express the data of health status of supporting rope, perhaps can think when the Cable Structure original state is the not damaged state vectorial d _oEach element numerical value get 0.Utilize the measured data of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, the Non-Destructive Testing data of supporting rope, temperature variant physics and the mechanical property parameters and the initial Cable Structure steady temperature data vector T of the employed various materials of Cable Structure _o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial Mechanics Calculation benchmark model A _o

No matter which kind of method to obtain initial Mechanics Calculation benchmark model A with _o, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T _o), based on A _oThe Cable Structure computational data that calculates must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this _oSuo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of calculating gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when truly taking place near institute's analog case reliably.Model A _oThe health status of middle supporting rope is with cable system initial damage vector d _oExpression, Cable Structure Cable Structure steady temperature data are with initial Cable Structure steady temperature data vector T _oExpression.Because based on A _oThe evaluation that calculates all monitored amounts is very near the initial value (actual measurement obtains) of all monitored amounts, so also can be used in A _oThe basis on, carry out Mechanics Calculation obtains, A _oThe evaluation of each monitored amount form monitored amount initial value vector C _oCorresponding to A _o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T _o"; Corresponding to A _oSupporting rope health status with cable system initial damage vector d _oExpression; Corresponding to A _oThe initial value of all monitored amounts with monitored amount initial value vector C _oExpression.T _oAnd d _oBe A _oParameter, C _oBy A _oMechanics Calculation result form.

The 3rd step: in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning ⁱ _o, A _oAnd A ⁱ _oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _oExpression, vector T ⁱ _oDefinition mode and vector T _oDefinition mode identical, T ⁱ _oElement and T _oElement corresponding one by one; The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d ⁱ _o, d ⁱ _oCable Structure A when representing this time circulation beginning ⁱ _oThe health status of cable system, d ⁱ _oDefinition mode and d _oDefinition mode identical, d ⁱ _oElement and d _oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C ⁱ _oExpression, vectorial C ⁱ _oDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _oElement and C _oElement corresponding one by one, the current initial value vector of monitored amount C ⁱ _oExpression is corresponding to A ⁱ _oThe concrete numerical value of all monitored amounts; T ⁱ _oAnd d ⁱ _oBe A ⁱ _oCharacterisitic parameter; C ⁱ _oBy A ⁱ _oMechanics Calculation result form; During circulation beginning for the first time, A ⁱ _oBe designated as A ¹ _o, set up A ¹ _oMethod for making A ¹ _oEqual A _oDuring circulation beginning for the first time, T ⁱ _oBe designated as T ¹ _o, set up T ¹ _oMethod for making T ¹ _oEqual T _oDuring circulation beginning for the first time, d ⁱ _oBe designated as d ¹ _o, set up d ¹ _oMethod for making d ¹ _oEqual d _oDuring circulation beginning for the first time, C ⁱ _oBe designated as C ¹ _o, set up C ¹ _oMethod for making C ¹ _oEqual C _o

The 4th step: the hardware components that the damaged cable recognition system is installed.Hardware components comprises at least: monitored amount monitoring system (for example containing spatial coordinate measuring 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, the computing machine and the panalarm of communicating by letter.Each monitored amount, each temperature all must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector with the signal that monitors; Signal is delivered to computing machine through signal picker; Computing machine then is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; When monitoring rope when damage is arranged, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.

The 5th step: establishment and on computers during the installation and operation temperature variation based on the approximant recognition methods system software of the damaged cable of space coordinate monitoring, this software will be accomplished functions such as monitoring that this method " during temperature variation based on the approximant recognition methods of the damaged cable of space coordinate monitoring " required by task wants, record, control, storage, calculating, notice, warning (in this practical implementation method all can with the work of computing machine completion).

The 6th step: the step begins the circulation running thus; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T ⁱ, vector T ⁱDefinition mode and vector T _oDefinition mode identical, T ⁱElement and T _oElement corresponding one by one; In the actual measurement vector T ⁱThe time, just obtaining current cable structure steady temperature data vector T ⁱThe synchronization in the moment, actual measurement obtains the currency of all monitored amounts in the Cable Structure, all these numerical value are formed the current numerical value vector of monitored amount C ⁱ, vectorial C ⁱDefinition mode and vectorial C _oDefinition mode identical, C ⁱElement and C _oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.

The 7th step: according to current cable structure steady temperature data vector T ⁱ, upgrade current initial Mechanics Calculation benchmark model A according to step a to b ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _o, and the current initial damage vector of cable system d ⁱ _oRemain unchanged;

A. be utilized in the T that obtains in the 6th step ⁱ, compare T ⁱAnd T ⁱ _oIf, T ⁱEqual T ⁱ _o, A then ⁱ _o, T ⁱ _o, C ⁱ _oAnd d ⁱ _oRemain unchanged;

B. be utilized in the T that obtains in the 6th step ⁱ, compare T ⁱAnd T ⁱ _oIf, T ⁱEqual T ⁱ _o, then need be to current initial Mechanics Calculation benchmark model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oUpgrade, update method is: at first calculate T ⁱWith T _oPoor, T ⁱWith T _oThe difference represent that with steady temperature change vector S S equals T ⁱDeduct T _oAt initial Mechanics Calculation benchmark model A _oThe basis on the health status of order supporting rope be the current initial damage vector of cable system d ⁱ _o, more further to A _oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A _oIn Cable Structure to apply what obtain after the temperature variation be exactly the current initial Mechanics Calculation benchmark model A that upgrades ⁱ _oUpgrading A ⁱ _oThe time, T ⁱ _oAll elements numerical value is used T ⁱAll elements numerical value replaces, and has promptly upgraded T ⁱ _o, so just obtained correctly corresponding to A ⁱ _oT ⁱ _o, this moment d ⁱ _oRemain unchanged; Upgrade A ⁱ _oAfter, obtain A through Mechanics Calculation ⁱ _oIn concrete numerical value all monitored amounts, current, replace vectorial C with these concrete numerical value ⁱ _oMiddle corresponding element has so just been realized vectorial C ⁱ _oRenewal, and the current initial damage of cable system vector d ⁱ _oRemain unchanged.

The 8th step: at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on, carry out the several times Mechanics Calculation according to step a to steps d, set up unit damage monitored numerical quantity transformation matrices Δ C through calculating ⁱWith nominal unit damage vector D ⁱ _u

A. the i time when beginning circulation, directly b obtains Δ C to the listed method of steps d set by step ⁱAnd D ⁱ _uAt other constantly, when in the 7th step to A ⁱ _oAfter upgrading, b regains Δ C to the listed method of steps d set by step ⁱAnd D ⁱ _uIf, the 7th the step in not to A ⁱ _oUpgrade, then directly changing for the 9th step herein over to carries out follow-up work.

B. at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on carry out the several times Mechanics Calculation; Equal the quantity of all supporting ropes on the calculation times numerical value; There is N root supporting rope that N calculating is just arranged; Calculating each time in the hypothesis cable system has only a supporting rope on the basis of original damage, to increase unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again; The supporting rope of appearance damage is different from the supporting rope that occurs damage in other time calculating in calculating each time, and supposes that each time the unit damage value of the supporting rope that damage is arranged can be different from the unit damage value of other supporting ropes, uses " nominal unit damage vector D ⁱ _u" write down the unit damage of the supposition of all ropes, vectorial D ⁱ _uElement coding rule and vectorial d _oThe coding rule of element identical, calculate the current numerical value of all monitored amounts in the Cable Structure each time, the current numerical value of the monitored amount of all that calculate is each time formed one " monitored amount is calculated current numerical value vector "; When the hypothesis j (j=1,2,3 ..., when N) root supporting rope has unit damage, available C ⁱ _Tj" monitored amount is calculated current numerical value vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C ⁱ _TjDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _TjElement and C _oElement corresponding one by one.

C. the vectorial C that calculates each time ⁱ _TjDeduct vectorial C ⁱ _oObtain a vector, obtain " numerical value change vector δ a C of monitored amount during each element that again should vector all calculates divided by this after the unit damage value of supposition ⁱ _j"; There is N root supporting rope that N " the numerical value change vector of monitored amount " just arranged.

D. form " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount " ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row corresponding 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 _oThe element coding rule identical.

The 9th step: set up linear relationship error vector e ⁱWith vectorial g ⁱUtilize data (" the current initial value vector of the monitored amount C of front ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ"), when the 8th step calculated each time, promptly in calculating each time, have only in the hypothesis cable system increase unit damage again on the basis of rope in original damage in; when hypothesis j (j=1,2,3; ...; when N) root supporting rope has unit damage, calculate each time and form a damage vector, use d ⁱ _TiExpression should damage vector, and corresponding monitored amount is calculated current numerical value vector and is C ⁱ _Tj(referring to the 8th step), damage vectorial d ⁱ _TjElement number equal the quantity of rope, vectorial d ⁱ _TjAll elements in have only the numerical value of an element to get to calculate each time in hypothesis increase the unit damage value of the rope of unit damage, d ⁱ _TjThe numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the rope of unit damage, be identical with the element of the same numberings of other vectors with the corresponding relation of this rope; d ⁱ _TjWith cable system initial damage vector d _oThe element coding rule identical, d ⁱ _TjElement and d _oElement be one-to-one relationship.With 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 ⁱ _je ⁱ _jSubscript j represent j (j=1,2,3 ..., N) root supporting rope has unit damage.There is N root rope that N calculating is just arranged, N linear relationship error vector e just arranged ⁱ _j, with this N linear relationship error vector e ⁱ _jObtain a vector after the addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e ⁱVector g ⁱEqual final error vector e ⁱWith vectorial g ⁱBe kept on the hard disc of computer of operation health monitoring systems software, supply health monitoring systems software to use.

$e_j^i=\mathrm{abs}(\mathrm{\&Delta;}{C}^i\&CenterDot;d_{\mathrm{tj}}^i-C_{\mathrm{tj}}^i+C_o^i)---\left(23\right)$

The tenth step: define the vectorial d of current name damage ⁱ _cWith current actual damage vector d ⁱ, d ⁱ _cAnd d ⁱElement number equal to support the quantity of rope, d ⁱ _cAnd d ⁱElement and supporting be one-to-one relationship between the rope, d ⁱ _cAnd d ⁱElement numerical value represent the degree of injury or the health status of corresponding supporting rope, d ⁱ _cAnd d ⁱWith cable system initial damage vector d _oThe element coding rule identical, d ⁱ _cElement, d ⁱElement and d _oElement be one-to-one relationship.

The 11 step: according to the current numerical value vector of monitored amount C ⁱWith " the current initial value vector of monitored amount C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " the vectorial d of current name damage ⁱ _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula (11), calculates the vectorial d of current name damage according to multi-objective optimization algorithm ⁱ _cNoninferior solution, just have reasonable error but can be more exactly from all ropes, confirm the position of damaged cable and separating of nominal degree of injury thereof.

The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted SUm Method), goal programming method (GoalAttainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution the vectorial d of current name damage with the goal programming method only ⁱ _cProcess, the concrete implementation procedure of other algorithm can realize according to the requirement of its specific algorithm in a similar fashion.

According to the goal programming method, formula (11) can transform the multi-objective optimization question shown in an accepted way of doing sth (24) and the formula (25), and γ is a real number in the formula (24), and R is a real number field, and area of space Ω has limited vectorial d ⁱ _cSpan (the present embodiment requirements vector d of each element ⁱ _cEach element be not less than 0, be not more than 1).The meaning of formula (24) is to seek the real number γ of a minimum, makes formula (25) be met.G (d in the formula (25) ⁱ _c) by formula (25) definition, G (d in the product representation formula (25) of weighing vector W and γ in the formula (25) ⁱ _c) and vectorial g ⁱBetween the deviation that allows, g ⁱDefinition referring to formula (17), its value the 9th the step calculate.During actual computation the vector W can with vectorial g ⁱIdentical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Use the goal programming method just can damage vectorial d in the hope of current name ⁱ _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)$

The 12 step: according to the current actual damage vector of cable system d ⁱThe definition (seeing formula (19)) of definition (seeing formula (18)) and its element calculate current actual damage vector d ⁱEach element, thereby can be by d ⁱConfirm the position and the degree of injury of damaged cable, just realized the health monitoring of cable system, realized damaged cable identification.d ⁱ _j(i=1,2,3, J=1,2,3 ...., N) the actual damage value of j root rope in the i time circulation of expression, formula (19), d are seen in its definition ⁱ _jBeing to represent j root rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%.

The 13 step: the computing machine in the health monitoring systems regularly generates cable system health condition form automatically or by the personnel operation health monitoring systems.

The 14 step: under specified requirements, the automatic operation communication panalarm of the computing machine in the health monitoring systems to monitor staff, owner and (or) personnel of appointment report to the police.

The 15 step: set up mark vector B according to formula (20) ⁱ, formula (21) has provided mark vector B ⁱThe definition of j element; If mark vector B ⁱElement be 0 entirely, then got back to for the 6th step and proceed health monitoring and calculating cable system; If mark vector B ⁱElement be not 0 entirely, then accomplish subsequent step after, get into circulation next time.

The 16 the step: according to formula (22) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation ^I+1 _oEach element d ^I+1 _Oj(j=1,2,3 ..., N); The second, at Mechanics Calculation benchmark model A _oThe basis on, make A _oIn the health status of rope be d ^I+1 _oRather than be d _oAfter, more further to A _oIn Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation that applies just taken from steady temperature change vector S, and steady temperature change vector S equals T ⁱDeduct T _o), so just obtained next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Mechanics Calculation benchmark mould A of circulation ^I+1 _o, next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation ^I+1 _oEqual T ⁱ _o, to A ^I+1 _oCarrying out Mechanics Calculation obtains corresponding to A ^I+1 _oConcrete numerical value all monitored amounts, current, these concrete numerical value form next time (promptly the i+1 time, i=1,2,3,4 ...) the current initial value vector C of the required monitored amount of circulation ^I+1 _o

The 17 the step: got back to for the 6th step, the beginning by the 6th go on foot the 17 the step circulation.

Claims (1)

1. During a temperature variation based on the approximant recognition methods of the damaged cable of space coordinate monitoring, it is characterized in that said method comprises:

   A. establish total N root supporting rope, at first confirm the coding rule of supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; Confirm the monitored point with monitored volume coordinate of appointment, give all monitored some numberings; Confirmed each monitored point with monitored volume coordinate component, give all monitored volume coordinate components numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored spatial data of Cable Structure " is made up of above-mentioned all monitored volume coordinate components; For simplicity, in the method " the monitored spatial data of Cable Structure " is called " monitored amount "; The quantity of monitored point must not be less than the quantity of supporting rope; The quantity sum of all monitored volume coordinate components must not be less than the quantity of supporting rope; Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly;

   B. this method definition " the temperature survey calculating method of the Cable Structure of this method " set by step b1 to b3 carry out;

   B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure; Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; One be called the cloudy day all day with what can not see the sun daytime in the method; Statistics obtain in T cloudy day each cloudy day 0 up to the sunrise next day highest temperature and the lowest temperature between back 30 minutes constantly, sunrise is meant the sunrise moment on the meteorology that base area revolutions and revolution rule confirm constantly, does not represent necessarily can see the sun same day; The sunrise that can inquire about data or calculate each required day through conventional meteorology constantly; Each cloudy day 0 up to next day sunrise constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, T cloudy day arranged, the maximum temperature difference of T cloudy daily temperature is just arranged; Get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day _rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation _h, for Δ T is got in convenient narration _hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data "; From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least 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 of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness "; Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is called reference plate; Reference plate can not contact with ground; Reference plate overhead distance is not less than 1.5 meters, and the one side of this reference plate is called sunny slope on the sunny side; The sunny slope of reference plate is coarse and dark color; The sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day, the non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate;

   B2: monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains Cable Structure place environment; Be called the environment maximum temperature difference, be designated as Δ T _EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day of temperature that minimum temperature obtains the sunny slope of reference plate; Be called the reference plate maximum temperature difference, be designated as Δ T _PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; There is R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise; Maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T _SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes; Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " 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 "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _Tmax

   B3: measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes 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 centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front _EmaxBe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _PmaxBe not more than Δ T after deducting 2 degrees centigrade _Emax, and Cable Structure surface maximum temperature difference Δ T _SmaxBe not more than Δ T _PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains 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 per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains 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 the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T _TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the 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 that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T _hThe numerical value that obtains is for Δ T is got in convenient narration _hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have 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 in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure;

   C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " and calculate the Cable Structure steady temperature data under the original state; Cable Structure steady temperature data under the 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 physics and the mechanical property parameters of the employed various materials of Cable Structure; Obtain T in actual measurement _oThe time, just obtaining initial Cable Structure steady temperature data vector T _oThe synchronization in the moment; Directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure comprises the Non-Destructive Testing data of the health status of expressing the supporting rope, the initial value of all monitored amounts, the initial rope force data of all supporting ropes, 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 is formed monitored amount initial value vector C _oUtilize the Non-Destructive Testing data of the health status that can express the supporting rope to set up cable system initial damage vector d _o, cable system initial damage vector d _oElement number equal N, d _oElement be one-to-one relationship with the supporting rope, cable system initial damage vector d _oElement numerical value represent the degree of injury of corresponding supporting rope, if cable system initial damage vector d _oThe numerical value of a certain element be 0, represent that the pairing supporting rope of this element is intact, not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%; Represent that then this supporting rope lost the load-bearing capacity of corresponding proportion; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, perhaps think when the Cable Structure original state is the not damaged state vectorial d _oEach element numerical value get 0;

   D. according to the temperature variant physics of the Non-Destructive Testing data of the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, supporting rope, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure steady temperature data vector T _oWith all Cable Structure data that preceding step obtains, set up the initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " _o, based on A _oThe Cable Structure computational data that calculates must be very near 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 _oSupporting rope health status with cable system initial damage vector d _oExpression; Corresponding to A _oThe initial value of all monitored amounts with monitored amount initial value vector C _oExpression; Set up the current initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " for the first time ⁱ _o, the current initial value of monitored amount vector C ⁱ _o" current initial Cable Structure steady temperature data vector T ⁱ _o"; Set up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oThe time, the current initial Mechanics Calculation benchmark model A of Cable Structure ⁱ _oJust equal the initial Mechanics Calculation benchmark model A of Cable Structure _o, the current initial value vector of monitored amount C ⁱ _oJust equal monitored amount initial value vector C _oA ⁱ _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 Mechanics Calculation benchmark model A of Cable Structure the first time ⁱ _oThe time, T ⁱ _oJust equal T _oA ⁱ _oThe initial health and the A of supporting rope _oThe health status of supporting rope identical, also use cable system initial damage vector d _oExpression, A in the cyclic process of back ⁱ _oThe initial health of supporting rope use cable system initial damage vector d all the time _oExpression; T _oAnd d _oBe A _oParameter, by A _oThe initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result _oThe initial value of all monitored amounts of expression is identical, therefore also we can say C _oBy A _oMechanics Calculation result form, work as T ⁱ _oAnd d _oBe A ⁱ _oParameter the time, C ^t _oBy A ⁱ _oMechanics Calculation result form; A in the method _o, C _o, d _oAnd T _oBe constant;

   E. in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning ⁱ _o, A _oAnd A ⁱ _oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A ⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T ⁱ _oExpression, vector T ⁱ _oDefinition mode and vector T _oDefinition mode identical, T ⁱ _oElement and T _oElement corresponding one by one; The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d ⁱ _o, d ⁱ _oCable Structure A when representing this time circulation beginning ⁱ _oThe health status of cable system, d ⁱ _oDefinition mode and d _oDefinition mode identical, d ⁱ _oElement and d _oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C ⁱ _oExpression, vectorial C ⁱ _oDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _oElement and C _oElement corresponding one by one, the current initial value vector of monitored amount C ⁱ _oExpression is corresponding to A ⁱ _oThe concrete numerical value of all monitored amounts; T ⁱ _oAnd d ⁱ _oBe A ⁱ _oCharacterisitic parameter; C ⁱ _oBy A ⁱ _oMechanics Calculation result form; During circulation beginning for the first time, A ⁱ _oBe designated as A ¹ _o, set up A ¹ _oMethod for making A ¹ _oEqual A _oDuring circulation beginning for the first time, T ⁱ _oBe designated as T ¹ _o, set up T ¹ _oMethod for making T ¹ _oEqual T _oDuring circulation beginning for the first time, d ⁱ _oBe designated as d ¹ _o, set up d ¹ _oMethod for making d ¹ _oEqual d _oDuring circulation beginning for the first time, C ⁱ _oBe designated as C ¹ _o, set up C ¹ _oMethod for making C ¹ _oEqual C _o

   F. go on foot the q circulation in step from getting into here by f; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T ⁱ, vector T ⁱDefinition mode and vector T _oDefinition mode identical, T ⁱElement and T _oElement corresponding one by one; Obtain vector T in actual measurement ⁱThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T ⁱThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value are formed the current numerical value vector of monitored amount C ⁱ, vectorial C ⁱDefinition mode and vectorial C _oDefinition mode identical, C ⁱElement and C _oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly;

   G. according to current cable structure steady temperature data vector T ⁱ, upgrade current initial Mechanics Calculation benchmark model A according to step g 1 to g2 ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _o, and the current initial damage vector of cable system d ⁱ _oRemain unchanged;

   G1. compare T ⁱAnd T ⁱ _oIf, T ⁱEqual T ⁱ _o, A then ⁱ _o, T ⁱ _o, C ⁱ _oAnd d ⁱ _oRemain unchanged;

   G2. compare T ⁱAnd T ⁱ _oIf, T ⁱBe not equal to T ⁱ _o, then need be to current initial Mechanics Calculation benchmark model A ⁱ _o, current initial Cable Structure steady temperature data vector T ⁱ _oWith the current initial value vector of monitored amount C ⁱ _oUpgrade, update method is: at first calculate T ⁱWith T _oPoor, T ⁱWith T _oThe difference represent that with steady temperature change vector S S equals T ⁱDeduct T _oAt initial Mechanics Calculation benchmark model A _oThe basis on the health status of order supporting rope be the current initial damage vector of cable system d ⁱ _o, more further to A _oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A _oIn Cable Structure to apply what obtain after the temperature variation be exactly the current initial Mechanics Calculation benchmark model A that upgrades ⁱ _oUpgrading A ⁱ _oThe time, T ⁱ _oAll elements numerical value is used T ⁱAll elements numerical value replaces, and has promptly upgraded T ⁱ _o, so just obtained correctly corresponding to A ⁱ _oT ⁱ _o, this moment d ⁱ _oRemain unchanged; Upgrade A ⁱ _oAfter, obtain A through Mechanics Calculation ⁱ _oIn concrete numerical value all monitored amounts, current, replace vectorial C with these concrete numerical value ⁱ _oMiddle corresponding element has so just been realized vectorial C ⁱ _oRenewal;

   H. at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity transformation matrices Δ C through calculating ⁱWith nominal unit damage vector D ⁱ _u

   H1. the i time when beginning circulation, directly h2 obtains Δ C to the listed method of step h4 set by step ⁱAnd D ⁱ _uAt other constantly, when in step g to A ⁱ _oAfter upgrading, h2 regains Δ C to the listed method of step h4 set by step ⁱAnd D ⁱ _uIf, in step g not to A ⁱ _oUpgrade, then directly change step I herein over to and carry out follow-up work;

   H2. at current initial Mechanics Calculation benchmark model A ⁱ _oThe basis on carry out the several times Mechanics Calculation; Equal the quantity of all supporting ropes on the calculation times numerical value; There is N root supporting rope that N calculating is just arranged; Calculating each time in the hypothesis cable system has only a supporting rope on the basis of original damage, to increase unit damage again; The supporting rope of appearance damage is different from the supporting rope that occurs damage in other time calculating in calculating each time, and supposes that each time the unit damage value of the supporting rope that damage is arranged can be different from the unit damage value of other supporting ropes, uses " nominal unit damage vector D ⁱ _u" write down the unit damage of the supposition of all ropes, vectorial D ⁱ _uElement coding rule and vectorial d _oThe coding rule of element identical, calculate the current numerical value of all monitored amounts in the Cable Structure each time, the current numerical value of the monitored amount of all that calculate is each time formed one " monitored amount is calculated current numerical value vector "; When hypothesis j root supporting rope has unit damage, available C ⁱ _Tj" monitored amount is calculated current numerical value vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C ⁱ _TjDefinition mode and vectorial C _oDefinition mode identical, C ⁱ _TjElement and C _oElement corresponding one by one;

   H3. the vectorial C that calculates each time ⁱ _TjDeduct vectorial C ⁱ _oObtain a vector, obtain " numerical value change vector δ a C of monitored amount during each element that again should vector all calculates divided by this after the unit damage value of supposition ⁱ _j"; There is N root supporting rope that N " the numerical value change vector of monitored amount " just arranged;

   H4. form " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged successively by this N " the numerical value change vector of monitored amount " ⁱ"; " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" each row corresponding 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 _oThe element coding rule identical;

   I. define the vectorial d of current name damage ⁱ _cWith current actual damage vector d ⁱ, d ⁱ _cAnd d ⁱElement number equal to support the quantity of rope, d ⁱ _cAnd d ⁱElement and supporting be one-to-one relationship between the rope, d ⁱ _cAnd d ⁱElement numerical value represent the degree of injury or the health status of corresponding supporting rope, d ⁱ _cAnd d ⁱWith cable system initial damage vector d _oThe element coding rule identical, d ⁱ _cElement, d ⁱElement and d _oElement be one-to-one relationship;

   J. according to the current numerical value vector of monitored amount C ⁱWith " the current initial value vector of monitored amount C ⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ C ⁱ" and " the vectorial d of current name damage ⁱ _c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 ⁱ _cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage ⁱ _c

   $C^i=C_o^i+\mathrm{\&Delta;}{C}^i\&CenterDot;d_c^i$ Formula 1

   K. the current actual damage vector d that utilizes formula 2 to express ⁱJ element d ⁱ _jWith the current initial damage vector of cable system d ⁱ _oJ element d ⁱ _OjWith the vectorial d of current name damage ⁱ _cJ 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 the formula 2,2,3 ..., N; Current actual damage vector d ⁱElement numerical value represent the degree of injury of corresponding supporting rope, so according to current actual damage vector d ⁱJust can define the impaired and degree of injury of which supporting rope, promptly realize the health monitoring of cable system in the Cable Structure, realize damaged cable identification; If the numerical value of a certain element of current actual damage vector is 0, represent that the pairing supporting rope of this element is intact, do not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity; If its numerical value between 0 and 100%, is then represented this supporting rope and has been lost the load-bearing capacity of corresponding proportion;

   L. try to achieve the vectorial d of current name damage ⁱ _cAfter, set up mark vector B according to formula 3 ⁱ, formula 4 has provided mark vector B ⁱThe definition of j element;

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

   

   formula 4

   Element B in the formula 4 ⁱ _jBe mark vector B ⁱJ element, D ⁱ _UjBe nominal unit damage vector D ⁱ _uJ element, d ⁱ _CjBe the vectorial d of the current name damage of cable system ⁱ _cJ element, they all represent the relevant information of j root supporting rope, j=1 in the formula 4,2,3 ..., N;

   If mark vector B m. ⁱElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B ⁱElement be not 0 entirely, then get into next step, be step n;

   N. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required cable system of circulation d ^I+1 _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 the formula 5 ^I+1 _OjBe the current initial damage vector of the required cable system d that next time, promptly circulates for the i+1 time ^I+1 _oJ element, d ⁱ _OjBe this, i.e. the i time current initial damage vector of round-robin cable system d ⁱ _oJ element, D ⁱ _UjBe the i time round-robin name unit damage vector D ⁱ _uJ element, B ⁱ _jBe the i time round-robin mark vector B ⁱJ element, j=1 in the formula 5,2,3 ..., N;

   O. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation ^I+1 _oEqual the current initial Cable Structure steady temperature data vector T of round-robin the i time ⁱ _o

   P. at initial Mechanics Calculation benchmark model A _oThe basis on, the health status that makes rope is d ^I+1 _oAfter, to initial Mechanics Calculation benchmark model A _oIn Cable Structure apply temperature variation, the steady temperature change vector S that the numerical value of the temperature variation that applies is just got this, be i.e. used in the i time circulation is to A _oIn Cable Structure apply obtain after the temperature variation be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation ^I+1Obtain A ^I+1After, obtain A through Mechanics Calculation ^I+1In concrete numerical value all monitored amounts, current, these concrete numerical value are formed next time, the vectorial C of the current initial value of required monitored amount that promptly circulates for the i+1 time ^I+1 _o

   Q. get back to step f, beginning is circulation next time.