CN103776644A - Method for identifying damaged cable, concentration load and support linear displacement through mixed monitoring - Google Patents

Abstract

The invention discloses a method for identifying a damaged cable, a concentration load and support linear displacement through mixed monitoring. The method is based on the mixed monitoring and comprises the following steps: judging whether it is required to update a mechanical calculation benchmark model of a cable structure by monitoring cable structure temperature and environment temperature; obtaining the mechanical calculation benchmark model of the cable structure, with the cable structure temperature and the environment temperature being taken into account; obtaining a unit damage monitored quantity value-variation matrix by calculation based on the mechanical calculation benchmark model; calculating the non-inferior solution of a current nominal damage vector of an evaluated object, according to the approximately linear relations between a current numeric vector of a monitored quantity and a current initial numeric vector of the monitored quantity, the unit damage monitored quantity value-variation matrix, and the to-be-solved current nominal damage vector of the evaluated object; and according to the non-inferior solution, identifying the support linear displacement, concentration load variance and the damaged cable during temperature variation.

Description

The damaged cable centre-point load support wire displacement recognition methods of hybrid monitoring

Technical field

Cable-stayed bridge, suspension bridge, the structures such as truss-frame structure have a common ground, be exactly that they have many parts that bear tensile load, as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or the rod member that only bears tensile load are support unit, for simplicity, this method is " Cable Structure " by such structure representation, and by all ropeway carrying-ropes of Cable Structure, carrying cable, and all rod members (being called again two power rod members) that only bear axial tension or axial compression load, unified be called " cable system " for simplicity, in this method, censure ropeway carrying-rope with " support cable " this noun, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word in the back, truss-frame structure reality is just referred to two power rod members.In structure military service process, the correct identification of the health status to support cable or cable system is related to the safety of whole Cable Structure.In the time that environment temperature changes, the temperature of Cable Structure generally also can be along with changing, in the time that Cable Structure temperature changes, may there is displacement of the lines in Cable Structure bearing, the centre-point load that Cable Structure is born also may change, the health status of Cable Structure also may change simultaneously, at this complex condition, based on hybrid monitoring, (this method judges the health status of Cable Structure to this method by the hybrid monitoring of the variation of the measurable parameter to the aforementioned dissimilar Cable Structure of this section, all monitored Cable Structure characteristic parameters are referred to as " monitored amount " by this method, because monitored amount is now mixed and formed by the dissimilar measurable parameter of Cable Structure, this method claims that this is hybrid monitoring) identify support wire displacement, the variable quantity of the centre-point load that damaged cable and Cable Structure are born, belong to engineering structure health monitoring field.

Background technology

Reject load change, the displacement of Cable Structure support wire and structure temperature and change the impact on Cable Structure health status recognition result, thereby the variation of the health status of recognition structure is exactly current problem in the urgent need to address; Rejecting load change, the variation of Cable Structure health status and structure temperature and change the impact on Cable Structure support wire displacement recognition result, thereby identify exactly the displacement of Cable Structure support wire, is also current problem in the urgent need to address; Same, the variation of rejecting structure temperature, the displacement of Cable Structure support wire and structural health conditions change the impact of the recognition result of the variable quantity of the centre-point load that structure is born, significant equally to structural safety.Based on structural health monitoring technology, this method discloses a kind of effective ways that solve these three problems.

When changing appears in the centre-point load of bearing when Cable Structure, or Cable Structure generation support wire displacement, or the temperature of Cable Structure is while changing, for example, or the health status of cable system is while changing (damaging), or when four kinds of situations occur simultaneously, can cause the variation of the measurable parameter of Cable Structure, for example can cause the variation of Suo Li, can affect distortion or the strain of Cable Structure, can affect shape or the volume coordinate of Cable Structure, can cause variation (for example variation of the angle coordinate of the straight line of any this point of mistake in the section of body structure surface any point of the angle coordinate of any imaginary line of the every bit of Cable Structure, or the variation of the angle coordinate of the normal of body structure surface any point), all these change the health status information that has all comprised cable system, also comprised the variable quantity information of support wire displacement and centre-point load, that is to say and can utilize the measurable parameter of Cable Structure to identify support wire displacement, the variable quantity of damaged cable and centre-point load.

In the time that bearing has displacement of the lines, current published technology, in method, some only can be in the time that other all conditions be constant the variation of (load of only only having structure to bear changes) recognition structure bearing load, the variation of some recognition structure health status of only can (only only having structural health conditions to change) in the time that other all conditions is constant, the variation of some only can (only only have structure temperature and structural health conditions to change) in the time that other all conditions is constant recognition structure (environment) temperature and structural health conditions, also do not have at present a kind of disclosed, effective method is recognition structure bearing load simultaneously, the variation of structure (environment) temperature and structural health conditions, when the load of bearing in structure in other words and structure (environment) temperature changes simultaneously, also there is no the variation that effective method can recognition structure health status, and the load that structure is born and structure (environment) temperature usually changes, so when the load of how to bear in structure and structure (environment) temperature variation, reject load change and structure temperature and change the impact on Cable Structure health status recognition result, thereby the variation of the health status of recognition structure exactly, it is current problem in the urgent need to address, this method discloses a kind of method, in the time that bearing has displacement of the lines, when the centre-point load that can bear in Cable Structure and structure (environment) temperature changes, reject support wire displacement, load change and structure temperature change the impact on Cable Structure health status recognition result, monitor to identify damaged cable based on monitored amount, the safety of Cable Structure is had to important value.

Same, in current disclosed method, thereby also do not occur rejecting the correct knowledge method for distinguishing of realizing centre-point load intensity of variation of support wire displacement, structure temperature variation and the impact of support cable health status, and concerning structure, the identification of load change is also very important; This method, in identifying damaged cable, can also identify the variation of centre-point load simultaneously, and this method can be rejected the impact that support wire displacement, structure temperature variation and support cable health status change, and realizes the correct identification of centre-point load intensity of variation; This method can also be rejected load change, Cable Structure health status changes and structure temperature changes the impact on Cable Structure support wire displacement recognition result, thereby identifies exactly the displacement of Cable Structure support wire.

That is to say, this method has realized three kinds of functions that existing method can not possess.

Summary of the invention

Technical matters: this method discloses a kind of method, three kinds of functions that existing method can not possess are realized, be respectively, one, reject support wire displacement, centre-point load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the health status of support cable; Two, this method, in identifying damaged cable, can also identify the variation of centre-point load simultaneously, and this method can be rejected the impact that support wire displacement, structure temperature variation and support cable health status change, and realizes the correct identification of centre-point load intensity of variation; Three, this method can also be rejected load change, Cable Structure health status changes and structure temperature changes the impact on Cable Structure support wire displacement recognition result, thereby identifies exactly the displacement of Cable Structure support wire.

Technical scheme: this method is made up of three parts.Be respectively set up the method for the required knowledge base of structural healthy monitoring system and parameter, based on knowledge base (containing parameter) with survey the structural health conditions appraisal procedure of monitored amount, the software and hardware part of health monitoring systems.

In the method, censure the coordinate of bearing about the X, Y, Z axis of Descartes's rectangular coordinate system with " bearing volume coordinate ", also can be said to is the volume coordinate of bearing about X, Y, Z axis, bearing is called the volume coordinate component of bearing about this axle about the concrete numerical value of the volume coordinate of some axles, and in this method, also a volume coordinate component with bearing is expressed the concrete numerical value of bearing about the volume coordinate of some axles; Censure the angular coordinate of bearing about X, Y, Z axis with " bearing angular coordinate ", bearing is called the angular coordinate component of bearing about this axle about the concrete numerical value of the angular coordinate of some axles, and in this method, also an angular coordinate component with bearing is expressed the concrete numerical value of bearing about the angular coordinate of some axles; All by " bearing generalized coordinate " denotion bearing angular coordinate and bearing volume coordinate, in this method, also a generalized coordinate component with bearing is expressed the concrete numerical value of bearing about volume coordinate or the angular coordinate of an axle; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, also can say that the change of bearing volume coordinate is called support wire displacement, and in this method, also a translational component with bearing is expressed the concrete numerical value of bearing about the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, and in this method, also an angular displacement component with bearing is expressed the concrete numerical value of bearing about the angular displacement of some axles; Generalized displacement of support denotion support wire displacement and angular displacement of support are all, and in this method, also a generalized displacement component with bearing is expressed bearing about the displacement of the lines of some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is support wire displacement or the translational displacement component at gravity direction.

First confirm the quantity of the centre-point load that may change that Cable Structure bears.The feature of the centre-point load of bearing according to Cable Structure, confirm wherein " centre-point load likely changing ", or all centre-point load is considered as " centre-point load likely changing ", establishes total JZW the centre-point load that may change.

Centre-point load is divided into two kinds of concentrated force and concentrated couples, in coordinate system, for example, in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, in the method a concentrated force component or a concentrated couple component is called to a centre-point load.

If the quantity sum of the quantity of the support wire displacement component of the quantity of the support cable of Cable Structure, Cable Structure and JZW " centre-point load likely changing " is N.For sake of convenience, it is " evaluation object " that this method unitedly calls evaluated support cable, support wire displacement and " centre-point load likely changing ", total N evaluation object.Give evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.

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

If total Q root support cable in cable system, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment rope ₁individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M ₁be one and be not less than 0 integer.

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

The monitored angle-data of Cable Structure is by K in Cable Structure ₃l individual specified point, that cross each specified point ₃the H of appointment straight line individual appointment straight line, each ₃individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments.Each total M ₃(M ₃=K ₃× L ₃× H ₃) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M ₃be one and be not less than 0 integer.

The monitored shape data of Cable Structure is by K in Cable Structure ₄l individual specified point and each specified point ₄the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K ₄the variation of all coordinate components of individual specified point.Each total M ₄(M ₄=K ₄× L ₄) individual measurement of coordinates value or calculated value characterize Cable Structure shape.M ₄be one and be not less than 0 integer.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M(M=M ₁+ M ₂+ M ₃+ M ₄) individual monitored amount, definition parameter K(K=M ₁+ K ₂+ K ₃+ K ₄), K and M must not be less than N.

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

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

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

The first step, inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of 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 of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, 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 _r.Between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of narration _hunit be ℃/m.On the surface of Cable Structure, get " R Cable Structure surface point ", the temperature of this R Cable Structure surface point will be obtained by actual measurement below, 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 by Calculation of Heat Transfer, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ".Get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " is narrated in the back with the condition that distribution must be satisfied.From the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at least choose two points 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 ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, measure Cable Structure along each and be no less than three points along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, 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 by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method.Measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual measurement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day, at the flat board of a carbon steel material of this position of sound production, be called reference plate, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained to the temperature of the sunny slope of reference plate.In this method, must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring, the moment of survey record data is called the physical record data moment.

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

The 3rd step, measures and calculates acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 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 Celsius, the b condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxafter deducting 2 degrees Celsius, be not more than Δ T _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax, only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition, Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 5 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 to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day, Section 6 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 Celsius, this method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to 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, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in 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 Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of 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 of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtaining, gets Δ T for convenience of narration _hunit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of 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 between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure.

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

Set up initial mechanical calculating benchmark model A _otime, in Cable Structure completion, or setting up before structural healthy monitoring system, calculating " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " measurement (can measure by conventional thermometry, for example use thermal resistance to measure), " Cable Structure steady temperature data " now use vector T _orepresent, be called initial Cable Structure steady temperature data vector T _o.Obtain T in actual measurement _otime, namely, obtaining the synchronization in moment of Cable Structure steady temperature data, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.Use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Obtain initial Cable Structure steady temperature data vector T at Actual measurement _otime, namely, obtaining the synchronization in moment of Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.First the Actual measurement data of Cable Structure are the data of the health status that can express rope including the Non-destructive Testing Data of support cable, and the Actual measurement data of Cable Structure are the measured data including Cable Structure bearing initial line displacement measurement data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinate data, Cable Structure modal data, structural strain data, structure angle measurement data, structure space measurement of coordinates data and load data still.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is to determine according to these coordinate datas the geometric properties of Cable Structure.Cable Structure bearing initial line displacement measurement data refer to setting up initial mechanical calculating benchmark model A _otime, the displacement of the lines that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point.For 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 bridge two ends, so-called bridge type data that Here it is.The variable quantity of " centre-point load likely changing " is being set up initial mechanical calculating benchmark model A _otime be all 0, the variable quantity that is to say " centre-point load likely changing " that identifies is below with respect to setting up initial mechanical calculating benchmark model A _otime the structure corresponding centre-point load of bearing variable quantity.The variable quantity data of utilizing the Non-destructive Testing Data etc. of support cable can express the data of the health status of support cable, Cable Structure bearing initial line displacement measurement data and " centre-point load likely changing " are set up evaluation object initial damage vector d _o(as the formula (1)), use d _orepresent that Cable Structure is (with initial mechanical calculating benchmark model A _orepresent) the initial health of evaluation object.If while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d _oin the each element numerical value relevant to support cable get 0.Vector d _oin the each element numerical value relevant to the variable quantity of centre-point load get 0.If there is no Cable Structure bearing initial line displacement measurement data or can think that the displacement of Cable Structure bearing initial line is at 0 o'clock, vectorial d _oin the each element numerical value relevant to the displacement of Cable Structure support wire get 0.Utilize temperature variant physical and mechanical properties parameter and the initial Cable Structure steady temperature data vector T of the various materials that Non-destructive Testing Data, Cable Structure bearing initial line displacement measurement data and the Cable Structure of the design drawing, as-constructed drawing of Cable Structure and the measured data of initial Cable Structure, support cable use _o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical calculating benchmark model A _o.

d _o＝[d _o1 d _o2 · · · d _ok · · · d _oN］ ^T (1)

D in formula (1) _ok(k=1,2,3 ...., N) represent initial mechanical calculating benchmark model A _oin the original state of k evaluation object, if this evaluation object is the rope (or pull bar) in cable system, d so _okrepresent its initial damage, d _okbe to represent not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion; If this evaluation object is translational component, so a d of a bearing _okrepresent its initial displacement numerical value; If this evaluation object is one " centre-point load that may change ", d so _okrepresent its initial value, d _okbe 0, the variable quantity that is to say " centre-point load likely changing " that identifies is below with respect to setting up initial mechanical calculating benchmark model A _otime the structure corresponding centre-point load of bearing variable quantity.Subscript T represents vectorial transposition (same afterwards).

Obtain T in actual measurement _otime, namely, obtaining the synchronization in moment of Cable Structure steady temperature data, use conventional method directly to measure the initial value of all monitored amounts of the Cable Structure calculating, form monitored amount initial value vector C _o(seeing formula (2)).Require obtaining A _otime obtain C _o, monitored amount initial value vector C _orepresent corresponding to A _othe concrete numerical value of " monitored amount ".Because of subject to the foregoing, the calculating benchmark model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in narration below, will represent this calculated value and measured value with prosign.

C _o＝[C _o1 C _o2 · · · C _oj · · · C _oM］ ^T (2)

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

No matter which kind of method to obtain initial mechanical calculating benchmark model A by _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 calculating must approach its measured data very much, and error generally must not be greater than 5%.Like this can utility A _ocalculate Suo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when approaching reliably institute's analog case and truly occurring.Model A _oevaluation object initial damage vector d for the health status of middle evaluation object _orepresent initial Cable Structure steady temperature data vector T for Cable Structure steady temperature data _orepresent.Due to based on A _othe evaluation that calculates all monitored amounts approaches the initial value (actual measurement obtains) of all monitored amounts very much, so also can be used in A _obasis on, carry out Mechanics Calculation obtains, A _othe evaluation of each monitored amount form monitored amount initial value vector C _o.T _oand d _oa _oparameter, also can say C _oby A _omechanics Calculation result composition.

Set up and upgrade current initial mechanical calculating benchmark model A ^t _omethod be: (namely set up for the first time A at initial time ^t _otime), A ^t _ojust equal A _o, A ^t _ocorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T ^t _o", at initial time, T ^t _ojust equal T _o, vector T ^t _odefinition mode and vector T _odefinition mode identical.A ^t _oinitial health and the A of evaluation object _othe health status of evaluation object identical, also use evaluation object initial damage vector d _orepresent A in cyclic process below ^t _othe initial health of evaluation object use all the time evaluation object initial damage vector d _orepresent; Cable Structure is in A ^t _owhen state, the current initial value vector of monitored amount C for this method ^t _orepresent the concrete numerical value of all monitored amounts, C ^t _oelement and C _oelement corresponding one by one, represent respectively all monitored amounts in Cable Structure in A ^t _oand A _oconcrete numerical value when two states.At initial time, C ^t _ojust equal C _o, T ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition; In Cable Structure military service process, the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T ^t", vector T ^tdefinition mode and vector T _odefinition mode identical); If T ^tequal T ^t _o, do not need A ^t _oupgrade, otherwise need to be to A ^t _oand T ^t _oupgrade, update method is: the first step is calculated T ^twith T _opoor, T ^twith T _odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ^twith T _opoor represent with steady temperature change vector S, S equals T ^tdeduct T _o, S represents the variation of Cable Structure steady temperature data; Second step is to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A _oin the temperature variation that applies of Cable Structure after obtain the current initial mechanical calculating benchmark model A that upgrades ^t _o, upgrade A ^t _otime, T ^t _oall elements numerical value is also used T ^tcorresponding replacement of all elements numerical value, upgraded T ^t _o, so just obtained correctly corresponding to A ^t _ot ^t _o; Upgrade C ^t _omethod be: when upgrade A ^t _oafter, obtain A by Mechanics Calculation ^t _oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C ^t _o.

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

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

C in formula (3) _j(j=1,2,3 ...., M) be the currency of j monitored amount in Cable Structure, this component C _jaccording to coding rule and C _ojcorresponding to same " monitored amount ".Obtain current cable structure steady temperature data vector T in actual measurement ^tsynchronization, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, forms monitored amount current value vector C.

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

Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A ^t _owith the current initial value vector of monitored amount C ^t _otime, upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C.Concrete grammar is as follows:

At the current initial mechanical calculating benchmark model A of Cable Structure ^t _obasis on carry out several times calculating, on calculation times numerical value, equal the quantity of all evaluation objects.Calculating each time hypothesis only has an evaluation object (to use vectorial d in original damage or displacement of the lines or centre-point load _ocorresponding element represent) basis on increase again unit line displacement, unit damage or centre-point load unit change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable has unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is a centre-point load, just suppose that this centre-point load is at vectorial d _o(if this centre-point load is couple, centre-point load unit change can be got 1kNm, 2kNm, 3kNm etc. for unit change on the basis of the existing variable quantity of this centre-point load representing, to increase centre-point load unit change again; If this centre-point load is concentrated force, centre-point load unit change can be got 1kN, 2kN, 3kN etc. for unit change), if this evaluation object is the translational component of a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d _oon the basis of the existing displacement of the lines of this bearing representing, there is again unit line displacement (10mm etc. are unit line displacement for for example 2mm, 5mm), use D _ukrecord this unit line displacement, unit damage or centre-point load unit change, wherein k represents the numbering of the evaluation object that unit line displacement, unit damage or centre-point load unit change occur.With " evaluation object unit change vector D _u" (as the formula (4)) record all unit line displacements, unit damage or centre-point load unit change.In calculating each time, there is unit line displacement, the evaluation object of unit damage or centre-point load unit change is different from other calculating and occurs unit line displacement, the evaluation object of unit damage or centre-point load unit change, calculate each time the current calculated value that all utilizes mechanics method (for example finite element method) to calculate all monitored amounts of Cable Structure, monitored amount calculation current vector of current calculated value composition of all monitored amounts that calculate each time is (when k evaluation object of hypothesis has unit line displacement, when unit damage or centre-point load unit change, available formula (5) represents monitored amount calculation current vector C _t ^k), calculate each time monitored amount calculation current vector C _t ^kdeduct the current initial value vector of monitored amount C ^t _oafter calculate divided by this time the unit line displacement, unit damage or the centre-point load unit change numerical value D that suppose again _ukgained vector be exactly under this condition (with have unit line displacement, unit damage or centre-point load unit change evaluation object be numbered mark) monitored amount unit change vector (in the time that k evaluation object has unit line displacement, unit damage or centre-point load unit change, use δ C _krepresent monitored amount unit change vector, formula (6) is shown in definition), each element representation of monitored amount unit change vector supposition when calculating has the Unit alteration amount of the corresponding monitored amount of this element that unit line displacement, unit damage or the centre-point load unit change of that evaluation object of unit line displacement, unit damage or centre-point load unit change cause, there is N evaluation object just to have N monitored amount unit change vector, owing to there being M monitored amount, so each monitored amount unit change vector has M element, form successively by this N monitored amount unit change vector the monitored amount unit change matrix Δ C that has M × N element, the definition of Δ C as the formula (6).

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

Evaluation object unit change vector D in formula (4) _uelement D _uk(k=1,2,3 ...., N) represent unit line displacement, unit damage or the centre-point load unit change numerical value of k evaluation object.

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

Elements C in formula (5) _tj ^k(k=1,2,3 ...., N; J=1,2,3 ...., M) represent while having unit line displacement, unit damage or centre-point load unit change due to k evaluation object, according to the current calculated amount of the individual monitored amount of the corresponding j of coding rule.

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

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

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

4. monitored amount current value vector C(calculates or actual measurement) with the current initial value vector of monitored amount C ^t _o, linear approximate relationship between unit damage monitored numerical quantity transformation matrices Δ C and the vectorial d of the current name damage of evaluation object, shown in (8) or formula (9).The definition of the vectorial d of the current name damage of evaluation object is referring to formula (10).

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

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

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

D in formula (10) _k(k=1,2,3 ...., N) be the current health status of k evaluation object in Cable Structure, if this evaluation object is the rope (or pull bar) in cable system, d so _krepresent its current damage, d _kbe to represent not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion, if this evaluation object is translational component, so a d of a bearing _krepresent its present bit shift value, if this evaluation object is centre-point load, so a d _krepresent its variable quantity.

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

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

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

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

Because formula (8) or the represented linear relationship of formula (9) exist certain error, therefore can not be simply carry out direct solution according to formula (8) or formula (9) and the monitored amount current value vector of actual measurement C and obtain the vectorial d of the current name damage of evaluation object.If done like this, the element in the vectorial d of the current name damage of evaluation object obtaining even there will be larger negative value, namely negative damage, and this is obviously irrational.Therefore obtain the acceptable solution of the vectorial d of the current name damage of evaluation object (with reasonable error, but can be more exactly from cable system, determine the position of damaged cable and degree of injury thereof, also can determine more exactly centre-point load variation numerical value, also can determine more exactly support wire displacement numerical value) become a rational solution, available formula (12) is expressed this method.

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

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

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

G in formula (13) _j(j=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (8) or formula (9) described.Vector g can be selected according to the error vector e tentative calculation of formula (11) definition.

At the current initial value vector of monitored amount C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, survey monitored amount current value vector C when known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to solve formula (12), obtain the acceptable solution of the vectorial d of the current name damage of evaluation object.

The current actual damage vector of definition evaluation object d ^a(seeing formula (14)), can be by d ^adetermine the health status of evaluation object.

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

D in formula (14) ^a _k(k=1,2,3 ...., N) represent the current actual health status of k evaluation object, if this evaluation object is the rope (or pull bar) in cable system, d so ^a _krepresent its current actual damage, formula (15), d are shown in its definition ^a _kbe to represent not damaged at 0 o'clock, while being 100%, represent that this rope thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion; If this evaluation object is a translational component of a bearing, formula (15), d are so shown in its definition ^a _krepresent its current actual line displacement numerical value; If this evaluation object is a centre-point load, formula (15), d are so shown in its definition ^a _krepresent that it is with respect to setting up initial mechanical calculating benchmark model A _otime the structure corresponding centre-point load of bearing variable quantity.Vector d ^athe coding rule of element and formula (1) in vectorial d _othe coding rule of element identical.

D in formula (15) _ok(k=1,2,3 ...., N) be vectorial d _ok element, d _kk the element of vectorial d.

So far this method has realized three kinds of functions that existing method can not possess, be respectively, one, can reject support wire displacement, centre-point load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of damaged cable; Two, can reject the impact that support wire displacement, structure temperature variation and support cable health status change, realize the correct identification of centre-point load intensity of variation; Three, can reject the impact that centre-point load variation, structure temperature variation and support cable health status change, realize the correct identification of support wire displacement.

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

Hardware components comprises monitoring system (comprising monitored amount monitoring system, temperature monitoring system), signal picker and computing machine etc.Require Real-Time Monitoring to obtain temperature required measured data, require each monitored amount of Real-Time Monitoring simultaneously.

That software should complete is needed in this method, can be by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.

This method specifically comprises:

A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated support wire displacement component, support cable and centre-point load, if the quantity sum of the quantity of evaluated support wire displacement component, the quantity of support cable and centre-point load is N, the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be used for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; While determining hybrid monitoring, specify by the support cable of monitored Suo Li, establish in cable system total Q root support cable, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment support cable ₁individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment support cables; Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M ₁be one and be not less than 0 integer that is not more than Q; While determining hybrid monitoring, specify by the measured point of monitored strain, the monitored strain data of Cable Structure can be by K in Cable Structure ₂l individual specified point and each specified point ₂the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K ₂the variation of all tested strains of individual specified point; Each total M ₂individual strain measurement value or calculated value characterize Cable Structure strain, M ₂for K ₂and L ₂long-pending; M ₂to be not less than 0 integer; While determining hybrid monitoring, specify by the measured point of monitored angle, the monitored angle-data of Cable Structure is by K in Cable Structure ₃l individual specified point, that cross each specified point ₃the H of appointment straight line individual appointment straight line, each ₃individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments; Each total M ₃individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M ₃for K ₃, L ₃and H ₃long-pending; M ₃be one and be not less than 0 integer; While determining hybrid monitoring, specify by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure ₄l individual specified point and each specified point ₄the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K ₄the variation of all coordinate components of individual specified point; Each total M ₄individual measurement of coordinates value or calculated value characterize Cable Structure shape, M ₄for K ₄and L ₄long-pending; M ₄be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M ₁, M ₂, M ₃and M ₄sum, definition parameter K, K is M ₁, K ₂, K ₃and K ₄sum, K and M must not be less than the quantity N of evaluation object; For simplicity, in the method this is walked to listed M monitored amount referred to as " monitored amount "; In this method, must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring, the moment of survey record data is called the physical record data moment;

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

B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of 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 of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, do not represent necessarily can see the same day sun, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, 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 _r, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of narration _hunit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, 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 by Calculation of Heat Transfer, just claim that the temperature data calculating 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, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at least choose two points 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 ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, measure Cable Structure along each and be no less than three points along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, 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 by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual measurement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day, at the flat board of a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained to the temperature of the sunny slope of reference plate,

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

B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 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 Celsius, the b condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxafter deducting 2 degrees Celsius, be not more than Δ T _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax, only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition, Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 5 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 to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day, Section 6 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 Celsius, this method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to 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, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in 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 Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of 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 of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtaining, gets Δ T for convenience of narration _hunit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of 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 between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure,

C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " the Cable Structure steady temperature data that calculate under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T _o", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, obtain T in actual measurement _otime, namely obtaining initial Cable Structure steady temperature data vector T _othe synchronization in moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure bearing initial line displacement measurement data, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing spatial data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, in obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C _o, monitored amount initial value vector C _othe coding rule of coding rule and M monitored amount identical, utilize support cable initial health data, Cable Structure bearing initial line displacement measurement data and Cable Structure centre-point load measurement data to set up evaluation object initial damage vector d _o, vectorial d _orepresent with initial mechanical calculating benchmark model A _othe initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d _oelement number equal N, d _oelement and evaluation object be one-to-one relationship, vectorial d _othe coding rule of element identical with the coding rule of evaluation object, if d _oevaluation object corresponding to some elements be support cable, so a d in cable system _othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, represents this support cable, if d _oevaluation object corresponding to some elements be some translational components of some bearings, d so _othe numerical value of this element represent the initial value of this translational component of this bearing, if d _oevaluation object corresponding to some elements be some centre-point load, in this method, get d _othis element numerical value be 0, the initial value that represents the variation of this centre-point load is 0, if there is no Cable Structure bearing initial line displacement measurement data or can think that the displacement of Cable Structure bearing initial line is at 0 o'clock, vectorial d _oin the each element numerical value relevant to the displacement of Cable Structure support wire get 0, if while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d _oin the each element numerical value relevant to support cable get 0, initial Cable Structure bearing spatial data refers to the bearing spatial data under Cable Structure design point, and Cable Structure bearing initial line displacement measurement data refer to setting up initial mechanical calculating benchmark model A _otime, the displacement of the lines that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point,

The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, Cable Structure bearing initial line displacement measurement data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, initial Cable Structure steady temperature data vector T _owith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " _o, based on A _othe Cable Structure computational data calculating must approach its measured data very much, 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 _oevaluation object initial damage vector d for evaluation object health status _orepresent; Corresponding to A _omonitored amount initial value vector C for the initial value of all monitored amounts _orepresent; Set up for the first time the current initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " ^t _o, the current initial value of monitored amount vector C ^t _o" current initial Cable Structure steady temperature data vector T ^t _o"; Set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure ^t _owith the current initial value vector of monitored amount C ^t _otime, the current initial mechanical calculating benchmark model A of Cable Structure ^t _ojust equal the initial mechanical calculating benchmark model A of Cable Structure _o, the current initial value vector of monitored amount C ^t _ojust equal monitored amount initial value vector C _o; A ^t _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 ^t _o", set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure ^t _otime, T ^t _ojust equal T _o; A ^t _oinitial health and the A of evaluation object _othe health status of evaluation object identical, also use evaluation object initial damage vector d _orepresent A in cyclic process below ^t _othe initial health of evaluation object use all the time evaluation object initial damage vector d _orepresent; T _oand d _oa _oparameter, by A _oinitial value and the C of all monitored amounts of obtaining of Mechanics Calculation result _othe initial value of all monitored amounts that represent is identical, therefore also can say C _oby A _omechanics Calculation result composition; T ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition;

E. from entering the circulation that is walked m step by e here; In structure military service process, constantly according to " the temperature survey calculating method of the Cable Structure of this method " the constantly current data of Actual measurement acquisition " Cable Structure steady temperature data ", the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T ^t", vector T ^tdefinition mode and vector T _odefinition mode identical;

F. according to current cable structure steady temperature data vector T ^t, upgrade current initial mechanical calculating benchmark model A according to step f1 to f3 ^t _o, the current initial value of monitored amount vector C ^t _owith current initial Cable Structure steady temperature data vector T ^t _o;

F1. compare T ^twith T ^t _oif, T ^tequal T ^t _o, A ^t _o, C ^t _oand T ^t _oremain unchanged; Otherwise need to follow these steps to A ^t _o, U ^t _oand T ^t _oupgrade;

F2. calculate T ^twith T _opoor, T ^twith T _odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ^twith T _opoor represent with steady temperature change vector S, S equals T ^tdeduct T _o, S represents the variation of Cable Structure steady temperature data;

F3. to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A _oin the temperature variation that applies of Cable Structure after obtain the current initial mechanical calculating benchmark model A that upgrades ^t _o, upgrade A ^t _otime, T ^t _oall elements numerical value is also used T ^tcorresponding replacement of all elements numerical value, upgraded T ^t _o, so just obtained correctly corresponding to A ^t _ot ^t _o; Upgrade C ^t _omethod be: when upgrade A ^t _oafter, obtain A by Mechanics Calculation ^t _oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C ^t _o; A ^t _othe initial health of support cable use all the time evaluation object initial damage vector d _orepresent;

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

G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A ^t _o, the current initial value of monitored amount vector C ^t _owith current initial Cable Structure steady temperature data vector T ^t _oafterwards, must then upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D _u;

G2. at the current initial mechanical calculating benchmark model A of Cable Structure ^t _obasis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or displacement of the lines or centre-point load, to increase unit damage or unit line displacement or centre-point load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d _oon the basis of the existing damage of this support cable representing, increase again unit damage, if this evaluation object is the translational component of a direction of a bearing, just suppose that this bearing increases unit line displacement again at this sense of displacement, if this evaluation object is a centre-point load, just suppose that this centre-point load is at vectorial d _oon the basis of the existing variable quantity of this centre-point load representing, increase again centre-point load unit change, use D _ukthe unit damage or unit line displacement or the centre-point load unit change that record this increase, wherein k represents the numbering of the evaluation object that increases unit damage or unit line displacement or centre-point load unit change, D _ukevaluation object unit change vector D _uan element, evaluation object unit change vector D _ucoding rule and the vectorial d of element _othe coding rule of element identical; The evaluation object that increases unit damage or unit line displacement or centre-point load unit change in calculating is each time different from the evaluation object that increases unit damage or unit line displacement or centre-point load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, a monitored amount calculation current vector of current calculated value composition of all monitored amounts that calculate each time, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C _oelement coding rule identical;

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

G4. by this N monitored amount unit change vector according to the coding rule of N evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is listed as successively; Each of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is listed as corresponding to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is the different unit change amplitude in the time that different evaluation objects increase unit damage or unit line displacement or centre-point load unit change corresponding to same monitored amount; Coding rule and the vectorial d of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C _othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;

H. obtain current cable structure steady temperature data vector T in actual measurement ^ttime, actual measurement obtains obtaining current cable structure steady temperature data vector T ^tthe current measured value of all monitored amounts of Cable Structure of synchronization in moment, form monitored amount current value vector C; The current initial value vector of monitored amount current value vector C and monitored amount C ^t _owith monitored amount initial value vector C _odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at concrete numerical value in the same time not;

I. define the vectorial d of the current name damage of evaluation object, the element number of the vectorial d of the current name damage of evaluation object equals the quantity of evaluation object, between the element of the vectorial d of the current name damage of evaluation object and evaluation object, be one-to-one relationship, the element numerical value of the vectorial d of the current name damage of evaluation object represents nominal degree of injury or nominal displacement of the lines or the nominal centre-point load variable quantity of corresponding evaluation object; Coding rule and the vectorial d of the element of vector d _othe coding rule of element identical;

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

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

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

L. the current actual damage vector of the evaluation object d that utilizes formula 2 to express ^ak element d ^a _kwith evaluation object initial damage vector d _ok element d _okk the element d with the vectorial d of the current name damage of evaluation object _kbetween relation, calculate the current actual damage of evaluation object vector d ^aall elements;

K=1 in formula 2,2,3 ...., N, d ^a _krepresent the current actual health status of k evaluation object, if this evaluation object is support cable, so a d in cable system ^a _krepresent its current actual damage, d ^a _kbe to represent not damaged at 0 o'clock, while being 100%, represent that this support cable thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion; If this evaluation object is translational component, so a d of a bearing ^a _krepresent its current actual line displacement numerical value; If this evaluation object is centre-point load, so a d ^a _krepresent the actual change amount of this centre-point load; So according to the current actual damage vector of evaluation object d ^acan define the impaired and degree of injury of which support cable, can define which bearing displacement of the lines and numerical value thereof have occurred, can define which centre-point load the numerical value that changes and change has occurred; So far this method has realized and has rejected damaged cable identification impact, Cable Structure that support wire displacement, centre-point load variation and structure temperature change, realize and rejected centre-point load variation, structure temperature variation and identification support cable health status variable effect, support wire displacement, realized and rejected support wire displacement, structure temperature variation and identification support cable health status variable effect, centre-point load variable quantity;

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

Beneficial effect: in current published correlation technique, some only can in the time that other all conditions is constant, (load of only only having structure to bear changes, and structural health conditions etc. are all constant) variation of recognition structure bearing load, some only can (only only have structural health conditions to change in the time that other all conditions is constant, and the load that structure is born etc. are constant) variation of recognition structure health status, some only can (only only have structure temperature and structural health conditions to change in the time that other all conditions is constant, and the load that structure is born is constant) variation of recognition structure health status, the load of bearing in structure, when structure (environment) temperature and structural health conditions change simultaneously, and while there is support wire displacement in Cable Structure simultaneously, also do not have at present a kind of disclosed, the effective method load that recognition structure bears simultaneously, the variation of support wire displacement and structural health conditions, this method has realized three kinds of functions that existing method can not possess, respectively: one, reject the displacement of Cable Structure support wire, centre-point load variation and structure temperature change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of damaged cable, two, reject the impact that the displacement of Cable Structure support wire, structure temperature variation and support cable health status change, realize the correct identification of centre-point load intensity of variation, three, reject the impact that centre-point load variation, structure temperature variation and support cable health status change, realize the correct identification of Cable Structure support wire displacement.

Embodiment

For cable structure health monitoring problem, this method has realized three kinds of functions that existing method can not possess, respectively: one, reject the displacement of Cable Structure support wire, centre-point load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of damaged cable; Two, reject the impact that the displacement of Cable Structure support wire, structure temperature variation and support cable health status change, realize the correct identification of centre-point load intensity of variation; Three, reject the impact that centre-point load variation, structure temperature variation and support cable health status change, realize the correct identification of Cable Structure support wire displacement.The following describes of the embodiment of this method is in fact only exemplary, and object is never to limit application or the use of this method.

This method adopts a kind of algorithm, and this algorithm is for identifying the variation of support wire displacement, damaged cable and centre-point load.When concrete enforcement, the following step is the one in the various steps that can take.

The first step: the quantity of first confirming the centre-point load that may change that Cable Structure bears.The feature of the centre-point load of bearing according to Cable Structure, confirm wherein " centre-point load likely changing ", or all centre-point load is considered as " centre-point load likely changing ", establishes total JZW the centre-point load that may change.

Centre-point load is divided into two kinds of concentrated force and concentrated couples, in coordinate system, for example, in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, in the method a concentrated force component or a concentrated couple component is called to a centre-point load.

If the quantity sum of the quantity of the quantity of the support wire displacement component of Cable Structure, the support cable of Cable Structure and JZW " centre-point load likely changing " is N.For sake of convenience, it is " evaluation object " that this method unitedly calls evaluated support wire displacement component, support cable and " centre-point load likely changing ", total N evaluation object.Give evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.

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

If total Q root support cable in cable system, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment rope ₁individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment ropes.Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M ₁be one and be not less than 0 integer.

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

The monitored angle-data of Cable Structure is by K in Cable Structure ₃l individual specified point, that cross each specified point ₃the H of appointment straight line individual appointment straight line, each ₃individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments.Each total M ₃(M ₃=K ₃× L ₃× H ₃) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M ₃be one and be not less than 0 integer.

The monitored shape data of Cable Structure is by K in Cable Structure ₄l individual specified point and each specified point ₄the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K ₄the variation of all coordinate components of individual specified point.Each total M ₄(M ₄=K ₄× L ₄) individual measurement of coordinates value or calculated value characterize Cable Structure shape.M ₄be one and be not less than 0 integer.

Comprehensive above-mentioned monitored amount, whole Cable Structure has M(M=M ₁+ M ₂+ M ₃+ M ₄) individual monitored amount, definition parameter K(K=M ₁+ K ₂+ K ₃+ K ₄), K and M must not be less than N.

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

Determine " the temperature survey calculating method of the Cable Structure of this method ", the method concrete steps are as follows:

A step: inquiry or actual measurement (can be measured by conventional thermometry, 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 of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, 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 _r, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of narration _hunit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, the temperature of this R Cable Structure surface point will be obtained by actual observation record below, 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 by Calculation of Heat Transfer, just claim that the temperature data calculating 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, in Cable Structure, uniform choosing is no less than three different sea level elevations, if the sea level elevation of for example Cable Structure is between 0m to 200m, can choose so height above sea level 0m, 50m, 100m and height above sea level 200m, the sea level elevation place choosing at each is crossing with Cable Structure surface with imaginary surface level, obtain intersection, surface level obtain cross surface crossing with Cable Structure, intersection is the outer edge line of cross surface, choose 6 points 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 ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness.In the measurement Cable Structure of choosing along in 6 directions of the Temperature Distribution of wall thickness, first according to the meteorological data throughout the year in region, Cable Structure position and the physical dimension of Cable Structure, volume coordinate, the sunny slope of definite Cable Structure such as Cable Structure surrounding environment and in the shade, the sunny slope of Cable Structure and in the shade face are surperficial parts for Cable Structure, the sea level elevation place choosing at each, aforementioned intersection respectively has one section in sunny slope and in the shade, two sections of these of intersection respectively have a mid point, cross these two mid points and get the outer normal of Cable Structure, these two outer normals are called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure by this method, these two outer normal directions are called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure by this method, the outer normal of obvious sunny slope and the outer normal of in the shade are all crossing with aforementioned intersection, also just there are two intersection points, intersection is divided into two line segments by these two intersection points, on two line segments, get respectively 2 points, totally 4 points, each line segment in two line segments of intersection is divided into equal in length 3 sections by taken point, get the outer normal on Cable Structure surface at these 4 some places, just chosen altogether like this outer normal on 6 Cable Structure surfaces 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, these two intersection points are on Cable Structure outside surface, another is on inside surface, if Cable Structure is solid, these two intersection points are all on Cable Structure outside surface, connect these two intersection points and obtain a straight-line segment, on straight-line segment, choose again three points, 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, amount to the temperature of 5 points, concrete can first hole in Cable Structure, how temperature sensor is embedded in to this 5 some places, the temperature recording is called this place " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".If chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, 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 by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method.Measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual observation record, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day (as long as had sunrise the same day, this position just should be arrived by solar radiation), for example, for example, at the flat board (the square flat board that the wide 3mm of 30cm is thick) of carbon steel material of this position of sound production (No. 45 carbon steels), be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, reference plate can be placed in the top of the wooden thermometer screen that meets meteorology temperature measurement requirement, the one side of this reference plate on the sunny side, (be for example called sunny slope, on the Northern Hemisphere time, sunny slope faces up towards south, full daytime is all by sunshine, sunny slope should have the suitable gradient to make snow can not accumulate or clear up sunny slope after snow), the sunny slope of reference plate is coarse and (being conducive to accept solar irradiation) dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material (the thick calcium carbonate insulation material of for example 5mm), Real-Time Monitoring record is obtained to the temperature of the sunny slope of reference plate.

B step, Real-Time Monitoring (can be measured by conventional thermometry, for example use thermal resistance to measure, for example, every temperature data of 10 minutes survey records) record obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring (can be measured by conventional thermometry simultaneously, for example use thermal resistance to measure, for example, every temperature data of 10 minutes survey records) obtain the temperature profile data of previously defined Cable Structure along thickness, Real-Time Monitoring (can be measured by conventional thermometry simultaneously, for example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature, for example, every temperature data of 10 minutes survey records) record obtains meeting meteorology and measures the temperature record of the Cable Structure place environment of temperature requirement, (can measure by conventional thermometry by 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 temperature data of 10 minutes survey records) record obtain being carved at sunrise the sunrise moment next day same day after the temperature measured data sequence of Cable Structure place environment between 30 minutes, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains Cable Structure place environment by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T _emax, temperature measured data sequence by Cable Structure place environment (for example first carries out curve fitting to the temperature measured data sequence of Cable Structure place environment by conventional mathematical computations, then by ask the derivative of curve to the time or by ask the rate of change of each point corresponding to survey record data time to the time on curve by numerical method) obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes, (can measure by conventional thermometry by Real-Time Monitoring, for example use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, for example, every temperature data of 10 minutes survey records) obtain being carved at sunrise the sunrise moment next day same day after the measured data sequence of temperature of sunny slope of reference plate between 30 minutes, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data that was carved at sunrise the temperature of the sunny slope of the reference plate between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the measured data sequence of temperature of the sunny slope of reference plate, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of the sunny slope of reference plate by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T _pmax, (can measure by conventional thermometry by Real-Time Monitoring, for example use thermal resistance to measure Cable Structure surface point, for example, every temperature data of 10 minutes survey records) record obtain being carved at sunrise the sunrise moment next day same day after the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between 30 minutes, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by being carved at sunrise the Cable Structure surface temperature measured data between 30 minutes after the sunrise moment next day same day of a Cable Structure surface point, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of each Cable Structure surface point by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T _smax, (for example first each Cable Structure surface temperature measured data sequence is carried out curve fitting by conventional mathematical computations by each Cable Structure surface temperature measured data sequence, then by ask the derivative of curve to the time or by ask the rate of change of each point corresponding to survey record data time to the time on curve by numerical method) obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the sea level elevation place that chooses at each and amount to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", choose H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T _tmax.

C step, measures and calculates acquisition Cable Structure steady temperature data; First, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology; The a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, Δ T _pmaxwith Δ T _smaxall be not more than 5 degrees Celsius; Section 2 must be satisfied b condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the Δ T that calculates above _emaxbe not more than with reference to temperature difference per day Δ T _r, and measure the Δ T calculating above _pmaxdeduct 2 degrees Celsius and be not more than Δ T _emax, and measure the Δ T calculating above _smaxbe not more than Δ T _pmax; Only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition; Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; Section 5 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 to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day; Section 6 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 Celsius.This method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to 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, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly the moment of the acquisition Cable Structure steady temperature data of this method, then, according to Cable Structure heat transfer characteristic, utilize R Cable Structure surface temperature measured data and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model (for example finite element model) of Cable Structure, for example, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer (finite element method), now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in 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 Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of 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 of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".Get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%; Cable Structure surface comprises support cable surface; Second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtaining, gets Δ T for convenience of narration _hunit be ℃/m that the unit of getting Δ h for convenience of narration is m; " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of 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 between two; The 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure.

Second step: set up initial mechanical calculating benchmark model A _o.

In Cable Structure completion, or setting up before health monitoring systems, calculating " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " measurement (can measure by conventional thermometry, for example use thermal resistance to measure), " Cable Structure steady temperature data " now use vector T _orepresent, be called initial Cable Structure steady temperature data vector T _o.Obtain T in actual measurement _otime, namely at the synchronization in 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 obtaining the synchronization in moment of so-and-so Cable Structure steady temperature data vector such as (such as initial or current) according to following method, use so-and-so method measurement to calculate the data of the monitored amount of so-and-so measured amount (all monitored amount of for example Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) time, for example, every temperature of 10 minutes survey records, so simultaneously equally also every 10 minutes the monitored amount of so-and-so measured amount of survey record (all monitored amount of for example Cable Structure) data.Once determine the moment that obtains Cable Structure steady temperature data, for example, be just called and obtaining the synchronization in moment of Cable Structure steady temperature data with the data of the monitored amount of so-and-so measured amount (all monitored amount of Cable Structure) of moment synchronization that obtain 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 temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the various materials that Cable Structure uses.

Obtain initial Cable Structure steady temperature data vector T at Actual measurement _otime, namely, obtaining the synchronization in moment of Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise that Non-destructive Testing Data of support cable etc. can express the data of the health status of rope, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinate data (initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data), Cable Structure bearing initial line displacement measurement data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structure angle measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing spatial data refers to the bearing spatial data under Cable Structure design point, and Cable Structure bearing initial line displacement measurement data refer to setting up initial mechanical calculating benchmark model A _otime, the displacement of the lines that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is to determine according to these coordinate datas the geometric properties of Cable Structure.For 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 bridge two ends, so-called bridge type data that Here it is.Data, Cable Structure bearing initial line displacement measurement data and the Cable Structure centre-point load measurement data of utilizing the Non-destructive Testing Data etc. of support cable can express the health status of support cable are set up evaluation object initial damage vector d _o(as the formula (1)), use d _orepresent that Cable Structure is (with initial mechanical calculating benchmark model A _orepresent) the initial health of evaluation object.If while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d _oin the each element numerical value relevant to support cable get 0; If there is no Cable Structure bearing initial line displacement measurement data or can think that the displacement of Cable Structure bearing initial line is at 0 o'clock, vectorial d _oin the each element numerical value relevant to the displacement of Cable Structure support wire get 0; If d _oevaluation object corresponding to some elements be some centre-point load, in this method, get d _othis element numerical value be 0, the initial value that represents the variation of this centre-point load is 0.Utilize the Non-destructive Testing Data of the design drawing, as-constructed drawing of Cable Structure and the measured data of initial Cable Structure, support cable, temperature variant physical and mechanical properties parameter and the initial Cable Structure steady temperature data vector T of various materials that Cable Structure is used _o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical calculating benchmark model A _o.

No matter which kind of method to obtain initial mechanical calculating benchmark model A by _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 calculating must approach its measured data very much, and error generally must not be greater than 5%.Like this can utility A _ocalculate Suo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when approaching reliably institute's analog case and truly occurring.Model A _oevaluation object initial damage vector d for the health status of middle support cable _orepresent initial Cable Structure steady temperature data vector T for Cable Structure steady temperature data _orepresent.Due to based on A _othe evaluation that calculates all monitored amounts approaches the initial value (actual measurement obtains) of all monitored amounts very much, so also can be used in A _obasis on, carry out Mechanics Calculation obtains, A _othe evaluation of each monitored amount form monitored amount initial value vector C _o.Corresponding to A _o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T _o", corresponding to A _oevaluation object initial damage vector d for evaluation object health status _orepresent, corresponding to A _omonitored amount initial value vector C for the initial value of all monitored amounts _orepresent T _oand d _oa _oparameter, C _oby A _omechanics Calculation result composition.

The 3rd step: set up for the first time current initial mechanical calculating benchmark model A ^t _o, the current initial value of monitored amount vector C ^t _o" current initial Cable Structure steady temperature data vector T ^t _o", concrete grammar is: at initial time, set up for the first time current initial mechanical calculating benchmark model A ^t _owith the current initial value vector of monitored amount C ^t _otime, A ^t _ojust equal A _o, C ^t _ojust equal C _o, A ^t _ocorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T ^t _o", (namely set up for the first time A at initial time ^t _otime), T ^t _ojust equal T _o, vector T ^t _odefinition mode and vector T _odefinition mode identical.A ^t _ohealth status and the A of evaluation object _ohealth status (the evaluation object initial damage vector d of evaluation object _orepresent) identical, A in cyclic process ^t _othe health status of evaluation object use all the time evaluation object initial damage vector d _orepresent.T ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition.

The 4th step: in Cable Structure military service process, the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T ^t", vector T ^tdefinition mode and vector T _odefinition mode identical).Obtain current cable structure steady temperature data vector T in actual measurement ^ttime, namely obtaining current cable structure steady temperature data vector T ^tthe synchronization in moment, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, composition " monitored amount current value vector C ".

The 5th step: according to current cable structure steady temperature data vector T ^t, upgrade where necessary current initial mechanical calculating benchmark model A ^t _o, the current initial value of monitored amount vector C ^t _owith current initial Cable Structure steady temperature data vector T ^t _o.Obtain current cable structure steady temperature data vector T in the 4th step actual measurement ^tafter, relatively T ^tand T ^t _oif, T ^tequal T ^t _o, do not need A ^t _oand T ^t _oupgrade, otherwise need to be to A ^t _oand T ^t _oupgrade, update method is undertaken by following a step to c step:

A step is calculated T ^twith T _opoor, T ^twith T _odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ^twith T _opoor represent with steady temperature change vector S, S equals T ^tdeduct T _o, S represents the variation of Cable Structure steady temperature data.

B walks A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A _oin Cable Structure apply and obtain the current initial mechanical calculating benchmark model A that upgrades after temperature variation ^t _o.

C step is upgraded A ^t _otime, T ^t _oall elements numerical value is also used T ^tcorresponding replacement of all elements numerical value, upgraded T ^t _o, so just obtained correctly corresponding to A ^t _ot ^t _o; Upgrade C ^t _omethod be: when upgrade A ^t _oafter, obtain A by Mechanics Calculation ^t _oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C ^t _o.

The 6th step: at current initial mechanical calculating benchmark model A ^t _obasis on carry out several times Mechanics Calculation, obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D by calculating _u.Concrete grammar is: Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A ^t _otime, must upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C simultaneously; At the current initial mechanical calculating benchmark model A of Cable Structure ^t _obasis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity of all evaluation objects, there is N evaluation object just to have N calculating, calculating each time hypothesis only has an evaluation object on the basis of original damage or original translational component or centre-point load, to increase unit damage or unit line displacement or centre-point load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d _oon the basis of the existing damage of this support cable representing, increase again unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the translational component of a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d _oon the basis of the existing displacement of the lines of this bearing representing, there is again unit line displacement (10mm etc. are unit line displacement for for example 2mm, 5mm), if this evaluation object is a centre-point load, just suppose that this centre-point load is at vectorial d _o(if this centre-point load is couple, centre-point load unit change can be got 1kNm, 2kNm, 3kNm etc. for unit change on the basis of the existing variable quantity of this centre-point load representing, to increase centre-point load unit change again; If this centre-point load is concentrated force, centre-point load unit change can be got 1kN, 2kN, 3kN etc. for unit change), use D _ukrecord this unit damage or centre-point load unit change, the numbering that wherein k represents that unit damage or unit line displacement occur or the evaluation object of centre-point load unit change occurs; The evaluation object that occurs unit damage or unit line displacement or centre-point load unit change in calculating is each time different from the evaluation object that occurs unit damage or unit line displacement or centre-point load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, a monitored amount calculation current vector C of current calculated value composition of all monitored amounts that calculate each time, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C _oelement coding rule identical; The monitored amount calculation current vector C calculating each time deducts the current initial value vector of monitored amount C ^t _oafter calculate divided by this time unit damage or unit line displacement or the centre-point load unit change numerical value supposed again, obtain a monitored amount unit change vector, have N evaluation object just to have the individual monitored amount unit change vector of N; Form successively by this N monitored amount unit change vector the unit damage monitored numerical quantity transformation matrices Δ C that has N row; Each of unit damage monitored numerical quantity transformation matrices is listed as corresponding to a monitored amount unit change vector, and every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is the different unit change amplitude when different evaluation object generation unit damage or unit line displacement or the centre-point load unit change corresponding to same monitored amount; Coding rule and the vectorial d of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C _othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount.

The 7th step: set up linear relationship error vector e and vectorial g.Utilize data (the current initial value vector of monitored amount C above ^t _o, unit damage monitored numerical quantity transformation matrices Δ C), when the 6th step is calculated each time, calculating each time hypothesis and only have in evaluation object increase unit damage or unit line displacement or the centre-point load unit change D of an evaluation object _uk, the evaluation object that increases unit damage or unit line displacement or centre-point load unit change in calculating is each time different from the evaluation object that increases unit damage or unit line displacement or centre-point load unit change in other calculating, calculate each time the current value of all utilizing mechanics method (for example adopting finite element method) to calculate all monitored amounts in Cable Structure, when calculating each time a monitored amount calculation current vector C of composition, calculate each time vectorial d of damage of composition, originally walking out of existing damage vector d only uses in this step, damaging in all elements of vectorial d only has the numerical value of an element to get D _uk, the numerical value of other element gets 0, damages coding rule and the vectorial d of the element of vectorial d _othe coding rule of element identical, by C, C ^t _o, Δ C, D _u, d brings formula (12) into, obtains a linear relationship error vector e, calculates each time a linear relationship error vector e, there is N evaluation object just to have N calculating, just there is N linear relationship error vector e, will this N linear relationship error vector e obtain a vector after being added, 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 equals final error vector e.

The 8th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (for example measuring subsystem, signal conditioner etc. containing measurement of angle subsystem, cable force measurement subsystem, strain measurement subsystem, volume coordinate), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, computing machine and the panalarm of communicating by letter.Each monitored amount, each temperature must arrive by monitored system monitoring, and the signal monitoring is transferred to signal (data) collector by monitoring system; Signal is delivered to computing machine through signal picker; Computing machine is responsible for the health monitoring software of the evaluation object of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; In the time monitoring evaluation object health status and change, computer control communication panalarm is reported to the police to the personnel of monitor staff, owner and (or) appointment.

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

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

The 11 step: the monitored amount current value vector of foundation C is with the current initial value vector of monitored amount C ^t _o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D _uand the vectorial d(of the current name damage of evaluation object is made up of all Suo Dangqian name amount of damage) between the linear approximate relationship (formula (8)) that exists, calculate the noninferior solution of the vectorial d of the current name damage of evaluation object according to multi-objective optimization algorithm, namely with reasonable error but can reflect more exactly the solution of the variation of the health status of evaluation object.

The multi-objective optimization algorithm that can adopt has a variety of, for example: the multiple-objection optimization based on genetic algorithm, the multiple-objection optimization based on artificial neural network, the multi-objective optimization algorithm based on population, the multiple-objection optimization based on ant group algorithm, leash law (Constrain Method), weighted method (Weighted Sum Method), Objective Programming (Goal Attainment Method) etc.Because various multi-objective optimization algorithms are all conventional algorithms, can realize easily, this implementation step only solves the process of the vectorial d of current damage take Objective Programming as example provides, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.

According to Objective Programming, formula (8) can transform the multi-objective optimization question shown in an accepted way of doing sth (16) and formula (17), in formula (16), γ is a real number, R is real number field, area of space Ω has limited the span (each element of the present embodiment requirements vector d is not less than 0, is not more than 1) of each element of vectorial d.The meaning of formula (16) is to find a minimum real number γ, and formula (17) is met.In formula (17), G (d) is defined by formula (18), the deviation allowing between the middle G (d) of the product representation formula (17) of weighing vector W and γ and vectorial g in formula (17), and the definition of g is referring to formula (13), and its value calculates in the 7th step.When actual computation, vector W can be identical with vectorial g.The concrete programming of Objective Programming realizes has had universal program directly to adopt.Use Objective Programming just can damage vectorial d in the hope of the current name of evaluation object.

$\begin{array}{cc}\min \mathrm{imize}& \mathrm{\γ}\\ \mathrm{\γ}\∈R,& d\∈\mathrm{\Ω}\end{array}---\left(16\right)$

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

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

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

The 12 step: the current actual damage vector of definition evaluation object d ^a, the current actual damage vector of evaluation object d ^aelement number equal the quantity of evaluation object, the current actual damage vector of evaluation object d ^aelement and evaluation object between be one-to-one relationship, the current actual damage of evaluation object vector d ^aelement numerical value represent actual damage degree or actual line displacement or the actual centre-point load intensity of variation of corresponding evaluation object; Vector d ^acoding rule and the vectorial d of element _othe coding rule of element identical.The current actual damage vector of the evaluation object d that utilizes formula (15) to express ^ak element d ^a _kwith evaluation object initial damage vector d _ok element d _okk the element d with the vectorial d of the current name damage of evaluation object _kbetween relation, calculate the current actual damage of evaluation object vector d ^aall elements.

D ^a _krepresent the current actual health status of k evaluation object, if this evaluation object is support cable, so a d in cable system ^a _krepresent its current actual damage, d ^a _kbe to represent not damaged at 0 o'clock, while being 100%, represent that this support cable thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion.

D ^a _krepresent the current actual health status of k evaluation object, if this evaluation object is translational component, so a d of a bearing ^a _krepresent its current actual line displacement numerical value.

D ^a _krepresent the current actual health status of k evaluation object, if this evaluation object is a centre-point load, formula (15), d are so shown in its definition ^a _krepresent that it is with respect to setting up initial mechanical calculating benchmark model A _otime the structure corresponding centre-point load of bearing variable quantity; So according to the current actual damage vector of evaluation object d ^acan define the impaired and degree of injury of which support cable, define which bearing displacement of the lines and numerical value thereof have occurred, can define which centre-point load variation and numerical value thereof have occurred simultaneously.

So far, can say that this method has realized three kinds of functions that existing method can not possess, respectively: one, reject the displacement of Cable Structure support wire, centre-point load variation and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of damaged cable; Two, reject the impact that the displacement of Cable Structure support wire, structure temperature variation and support cable health status change, realized the correct identification of centre-point load intensity of variation; Three, reject the impact that centre-point load variation, structure temperature variation and support cable health status change, realized the correct identification of Cable Structure support wire displacement.

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

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

The 15 step: get back to the 4th step, start the circulation by the 4th step to the 15 steps.

Claims (1)

1. the damaged cable centre-point load support wire displacement recognition methods of hybrid monitoring, is characterized in that described method comprises:

A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated support wire displacement component, support cable and centre-point load, if the quantity sum of the quantity of evaluated support wire displacement component, the quantity of support cable and centre-point load is N, the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be used for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; While determining hybrid monitoring, specify by the support cable of monitored Suo Li, establish in cable system total Q root support cable, the monitored rope force data of Cable Structure is by M in Cable Structure ₁the M of individual appointment support cable ₁individual rope force data is described, and the variation of Cable Structure Suo Li is exactly the variation of the Suo Li of all appointment support cables; Each total M ₁individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M ₁be one and be not less than 0 integer that is not more than Q; While determining hybrid monitoring, specify by the measured point of monitored strain, the monitored strain data of Cable Structure can be by K in Cable Structure ₂l individual specified point and each specified point ₂the strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K ₂the variation of all tested strains of individual specified point; Each total M ₂individual strain measurement value or calculated value characterize Cable Structure strain, M ₂for K ₂and L ₂long-pending; M ₂to be not less than 0 integer; While determining hybrid monitoring, specify by the measured point of monitored angle, the monitored angle-data of Cable Structure is by K in Cable Structure ₃l individual specified point, that cross each specified point ₃the H of appointment straight line individual appointment straight line, each ₃individual angle coordinate component is described, and the variation of Cable Structure angle is exactly the variation of angle coordinate components appointment straight lines all specified points, all, all appointments; Each total M ₃individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M ₃for K ₃, L ₃and H ₃long-pending; M ₃be one and be not less than 0 integer; While determining hybrid monitoring, specify by monitored shape data, the monitored shape data of Cable Structure is by K in Cable Structure ₄l individual specified point and each specified point ₄the volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K ₄the variation of all coordinate components of individual specified point; Each total M ₄individual measurement of coordinates value or calculated value characterize Cable Structure shape, M ₄for K ₄and L ₄long-pending; M ₄be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M ₁, M ₂, M ₃and M ₄sum, definition parameter K, K is M ₁, K ₂, K ₃and K ₄sum, K and M must not be less than the quantity N of evaluation object; For simplicity, in the method this is walked to listed M monitored amount referred to as " monitored amount "; In this method, must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring, the moment of survey record data is called the physical record data moment;

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

B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of 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 of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, do not represent necessarily can see the same day sun, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, 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 _r, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation _h, get Δ T for convenience of narration _hunit be ℃/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, 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 by Calculation of Heat Transfer, just claim that the temperature data calculating 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, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at least choose two points 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 ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, measure Cable Structure along each and be no less than three points along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, 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 by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual measurement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day, at the flat board of a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained to the temperature of the sunny slope of reference plate,

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

B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 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 Celsius, the b condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates above _emaxbe not more than with reference to temperature difference per day Δ T _r, and reference plate maximum temperature difference Δ T _pmaxafter deducting 2 degrees Celsius, be not more than Δ T _emax, and Cable Structure surface maximum temperature difference Δ T _smaxbe not more than Δ T _pmax, only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition, Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 5 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 to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day, Section 6 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 Celsius, this method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to 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, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in 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 Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of 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 of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 ℃ divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation _hthe numerical value obtaining, gets Δ T for convenience of narration _hunit be ℃/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of 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 between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure,

C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " the Cable Structure steady temperature data that calculate under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T _o", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, obtain T in actual measurement _otime, namely obtaining initial Cable Structure steady temperature data vector T _othe synchronization in moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure bearing initial line displacement measurement data, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing spatial data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, in obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C _o, monitored amount initial value vector C _othe coding rule of coding rule and M monitored amount identical, utilize support cable initial health data, Cable Structure bearing initial line displacement measurement data and Cable Structure centre-point load measurement data to set up evaluation object initial damage vector d _o, vectorial d _orepresent with initial mechanical calculating benchmark model A _othe initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d _oelement number equal N, d _oelement and evaluation object be one-to-one relationship, vectorial d _othe coding rule of element identical with the coding rule of evaluation object, if d _oevaluation object corresponding to some elements be support cable, so a d in cable system _othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, represents this support cable, if d _oevaluation object corresponding to some elements be some translational components of some bearings, d so _othe numerical value of this element represent the initial value of this translational component of this bearing, if d _oevaluation object corresponding to some elements be some centre-point load, in this method, get d _othis element numerical value be 0, the initial value that represents the variation of this centre-point load is 0, if there is no Cable Structure bearing initial line displacement measurement data or can think that the displacement of Cable Structure bearing initial line is at 0 o'clock, vectorial d _oin the each element numerical value relevant to the displacement of Cable Structure support wire get 0, if while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d _oin the each element numerical value relevant to support cable get 0, initial Cable Structure bearing spatial data refers to the bearing spatial data under Cable Structure design point, and Cable Structure bearing initial line displacement measurement data refer to setting up initial mechanical calculating benchmark model A _otime, the displacement of the lines that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point,

The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, Cable Structure bearing initial line displacement measurement data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, initial Cable Structure steady temperature data vector T _owith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " _o, based on A _othe Cable Structure computational data calculating must approach its measured data very much, 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 _oevaluation object initial damage vector d for evaluation object health status _orepresent; Corresponding to A _omonitored amount initial value vector C for the initial value of all monitored amounts _orepresent; Set up for the first time the current initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " ^t _o, the current initial value of monitored amount vector C ^t _o" current initial Cable Structure steady temperature data vector T ^t _o"; Set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure ^t _owith the current initial value vector of monitored amount C ^t _otime, the current initial mechanical calculating benchmark model A of Cable Structure ^t _ojust equal the initial mechanical calculating benchmark model A of Cable Structure _o, the current initial value vector of monitored amount C ^t _ojust equal monitored amount initial value vector C _o; A ^t _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 ^t _o", set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure ^t _otime, T ^t _ojust equal T _o; A ^t _oinitial health and the A of evaluation object _othe health status of evaluation object identical, also use evaluation object initial damage vector d _orepresent A in cyclic process below ^t _othe initial health of evaluation object use all the time evaluation object initial damage vector d _orepresent; T _oand d _oa _oparameter, by A _oinitial value and the C of all monitored amounts of obtaining of Mechanics Calculation result _othe initial value of all monitored amounts that represent is identical, therefore also can say C _oby A _omechanics Calculation result composition; T ^t _oand d _oa ^t _oparameter, C ^t _oby A ^t _omechanics Calculation result composition;

E. from entering the circulation that is walked m step by e here; In structure military service process, constantly according to " the temperature survey calculating method of the Cable Structure of this method " the constantly current data of Actual measurement acquisition " Cable Structure steady temperature data ", the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T ^t", vector T ^tdefinition mode and vector T _odefinition mode identical;

F. according to current cable structure steady temperature data vector T ^t, upgrade current initial mechanical calculating benchmark model A according to step f1 to f3 ^t _o, the current initial value of monitored amount vector C ^t _owith current initial Cable Structure steady temperature data vector T ^t _o;

F1. compare T ^twith T ^t _oif, T ^tequal T ^t _o, A ^t _o, C ^t _oand T ^t _oremain unchanged; Otherwise need to follow these steps to A ^t _o, U ^t _oand T ^t _oupgrade;

F2. calculate T ^twith T _opoor, T ^twith T _odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T ^twith T _opoor represent with steady temperature change vector S, S equals T ^tdeduct T _o, S represents the variation of Cable Structure steady temperature data;

F3. to A _oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A _oin the temperature variation that applies of Cable Structure after obtain the current initial mechanical calculating benchmark model A that upgrades ^t _o, upgrade A ^t _otime, T ^t _oall elements numerical value is also used T ^tcorresponding replacement of all elements numerical value, upgraded T ^t _o, so just obtained correctly corresponding to A ^t _ot ^t _o; Upgrade C ^t _omethod be: when upgrade A ^t _oafter, obtain A by Mechanics Calculation ^t _oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C ^t _o; A ^t _othe initial health of support cable use all the time evaluation object initial damage vector d _orepresent;

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

G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A ^t _o, the current initial value of monitored amount vector C ^t _owith current initial Cable Structure steady temperature data vector T ^t _oafterwards, must then upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D _u;

G2. at the current initial mechanical calculating benchmark model A of Cable Structure ^t _obasis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or displacement of the lines or centre-point load, to increase unit damage or unit line displacement or centre-point load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d _oon the basis of the existing damage of this support cable representing, increase again unit damage, if this evaluation object is the translational component of a direction of a bearing, just suppose that this bearing increases unit line displacement again at this sense of displacement, if this evaluation object is a centre-point load, just suppose that this centre-point load is at vectorial d _oon the basis of the existing variable quantity of this centre-point load representing, increase again centre-point load unit change, use D _ukthe unit damage or unit line displacement or the centre-point load unit change that record this increase, wherein k represents the numbering of the evaluation object that increases unit damage or unit line displacement or centre-point load unit change, D _ukevaluation object unit change vector D _uan element, evaluation object unit change vector D _ucoding rule and the vectorial d of element _othe coding rule of element identical; The evaluation object that increases unit damage or unit line displacement or centre-point load unit change in calculating is each time different from the evaluation object that increases unit damage or unit line displacement or centre-point load unit change in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, a monitored amount calculation current vector of current calculated value composition of all monitored amounts that calculate each time, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C _oelement coding rule identical;

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

G4. by this N monitored amount unit change vector according to the coding rule of N evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is listed as successively; Each of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is listed as corresponding to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is the different unit change amplitude in the time that different evaluation objects increase unit damage or unit line displacement or centre-point load unit change corresponding to same monitored amount; Coding rule and the vectorial d of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C _othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;

H. obtain current cable structure steady temperature data vector T in actual measurement ^ttime, actual measurement obtains obtaining current cable structure steady temperature data vector T ^tthe current measured value of all monitored amounts of Cable Structure of synchronization in moment, form monitored amount current value vector C; The current initial value vector of monitored amount current value vector C and monitored amount C ^t _owith monitored amount initial value vector C _odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at concrete numerical value in the same time not;

I. define the vectorial d of the current name damage of evaluation object, the element number of the vectorial d of the current name damage of evaluation object equals the quantity of evaluation object, between the element of the vectorial d of the current name damage of evaluation object and evaluation object, be one-to-one relationship, the element numerical value of the vectorial d of the current name damage of evaluation object represents nominal degree of injury or nominal displacement of the lines or the nominal centre-point load variable quantity of corresponding evaluation object; Coding rule and the vectorial d of the element of vector d _othe coding rule of element identical;

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

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

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

L. the current actual damage vector of the evaluation object d that utilizes formula 2 to express ^ak element d ^a _kwith evaluation object initial damage vector d _ok element d _okk the element d with the vectorial d of the current name damage of evaluation object _kbetween relation, calculate the current actual damage of evaluation object vector d ^aall elements;

K=1 in formula 2,2,3 ...., N, d ^a _krepresent the current actual health status of k evaluation object, if this evaluation object is support cable, so a d in cable system ^a _krepresent its current actual damage, d ^a _kbe to represent not damaged at 0 o'clock, while being 100%, represent that this support cable thoroughly loses load-bearing capacity, between 0 and 100% time, represent to lose the load-bearing capacity of corresponding proportion; If this evaluation object is translational component, so a d of a bearing ^a _krepresent its current actual line displacement numerical value; If this evaluation object is centre-point load, so a d ^a _krepresent the actual change amount of this centre-point load; So according to the current actual damage vector of evaluation object d ^acan define the impaired and degree of injury of which support cable, can define which bearing displacement of the lines and numerical value thereof have occurred, can define which centre-point load the numerical value that changes and change has occurred; So far this method has realized and has rejected damaged cable identification impact, Cable Structure that support wire displacement, centre-point load variation and structure temperature change, realize and rejected centre-point load variation, structure temperature variation and identification support cable health status variable effect, support wire displacement, realized and rejected support wire displacement, structure temperature variation and identification support cable health status variable effect, centre-point load variable quantity;

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