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

Abstract

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

Description

Simplify generalized displacement hybrid monitoring problem cable load progressive-type recognition method

Technical field

Cable-stayed bridge, suspension bridge, the structures such as truss structure have a common ground, it is exactly that they have many parts for bearing tensile load, such as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of the class formation is with rope, cable is subjected only to the rod member of tensile load for supporting member, for convenience, such structure representation is " Cable Structure " by this method, and by all carrying ropes of Cable Structure, carry cable, and all rod members (also known as two power rod members) for being subjected only to axial tension or axial compression load, it is collectively referred to as convenience " cable system ", in this method carrying rope is censured with " support cable " this noun, carrying cable and the rod member for being subjected only to axial tension or axial compression load, sometimes referred to simply as " rope ", so when " rope " this word is used below, two power rod members are just actually referred to truss structure.The impaired and lax pair Cable Structure of support cable is a significant threat safely, and damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method.During structure military service, the correct identification to support cable or the health status of cable system is related to the safety of whole Cable Structure.When environment temperature changes, the temperature of Cable Structure typically also can be with changing, when Cable Structure temperature changes, generalized displacement may occur for Cable Structure bearing, the load that Cable Structure is born is it can also happen that change, the health status of Cable Structure may also change simultaneously, in 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 change of the measurable parameter to the foregoing different types of Cable Structure of this section, all monitored Cable Structure characteristic parameters are referred to as " monitored amount " by this method, it is made up of due to being now monitored amount the different types of measurable parameter mixing of Cable Structure, this method is referred to as hybrid monitoring) recognize problem cable, belong to engineering structure health monitoring field.

Background technology

Reject the influence of load change, Cable Structure generalized displacement of support and structure temperature change to Cable Structure health status recognition result, so that the change of the health status of identification structure exactly, the problem of being in the urgent need to address at present, this method discloses a kind of effective, the cheap method for solving this problem.

The content of the invention

Technical problem：This method discloses a kind of method, under conditions of cost is lower, when bearing has generalized displacement, during load and structure (environment) temperature change born in structure, the influence of generalized displacement of support, load change and structure temperature change to Cable Structure health status recognition result can be rejected, so as to identify the health status of support cable exactly.

Technical scheme：In the method, coordinate of the bearing on the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing space coordinate ", space coordinate of the bearing on X, Y, Z axis can also be said to be, bearing is referred to as in space coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the space coordinate of some axle with a space coordinate subscale of bearing on the concrete numerical value of the space coordinate of some axle；Angular coordinate of the bearing on X, Y, Z axis is censured with " bearing angular coordinate ", bearing is referred to as in angular coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the angular coordinate of some axle with an angular coordinate subscale of bearing on the concrete numerical value of the angular coordinate of some axle；Censure that bearing angular coordinate and bearing space coordinate are all with " bearing generalized coordinates ", in this method also with a generalized coordinates subscale of bearing up to bearing on the space coordinate of axle or the concrete numerical value of angular coordinate；Bearing is referred to as bearing displacement of the lines on the change of the coordinate of X, Y, Z axis, it may also be said to which the change of bearing space coordinate is referred to as the concrete numerical value in bearing displacement of the lines, this method also with the translational component expression bearing of bearing on the displacement of the lines of some axle；Bearing is referred to as the concrete numerical value in angular displacement of support, this method also with the angular displacement component expression bearing of bearing on the angular displacement of some axle on the change of the angular coordinate of X, Y, Z axis；Generalized displacement of support censures bearing displacement of the lines and angular displacement of support is all, and also bearing is reached on the displacement of the lines of some axle or the concrete numerical value of angular displacement with a generalized displacement subscale of bearing in this method；Bearing displacement of the lines is alternatively referred to as translational displacement, and support settlement is the component of bearing displacement of the lines or translational displacement in gravity direction.

The external force that object, structure are born can be described as load, and load includes face load and volume load.Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load.Volume load is the continuously distributed load in interior of articles each point, the deadweight of such as object and inertia force.

Concentrfated load is divided into two kinds of concentrated force and concentrated couple, in a coordinate system, for example in Descartes's rectangular coordinate system, one concentrated force can resolve into three components, likewise, a concentrated couple can also resolve into three components, if load is actually concentrfated load, force component or a concentrated couple component is concentrated to be referred to as a load by one in the method, the now change of load is embodied as a change for concentrating force component or a concentrated couple component.

Distributed load is divided into line distributed load and EDS maps load, the size of the zone of action and distributed load of the description of distributed load at least including distributed load, the size of distributed load is expressed with distribution intensity, being distributed intensity, to express, (such as two distributed loads are all uniform with distribution characteristics (such as uniform, SIN function equal distribution feature) and amplitude, but its amplitude is different, and the concept of amplitude can be illustrated by taking well-distributed pressure as an example：Same structure bears two different well-distributed pressures, and two distributed loads are all uniform loads, but the amplitude of a distributed load is 10MPa, and the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, one distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is not all identical, the component of this several distributed load is so regarded as same amount of independent distributed load in the method, now a load just represents the component of a distributed load, the amplitude changing ratio identical component of wherein distribution intensity can also be synthesized into a distributed load or is load.

Volume load is the continuously distributed load in interior of articles each point, deadweight and inertia force such as object, the size of the zone of action and volume load of the description of volume load at least including volume load, the size of volume load is expressed with distribution intensity, with distribution characteristics (such as uniform, linear function equal distribution feature) and amplitude, to express, (it is all uniform that lotuses are storaged in such as two individuals to distribution intensity, but its amplitude is different, and the concept of amplitude is illustrated exemplified by can conducting oneself with dignity：The material of two parts of same structure is different, therefore density is different, so while the volume load suffered by the two parts is all uniform, but the amplitude of the volume load suffered by a part is probably 10kN/m³, the amplitude of the volume load suffered by another part is 50kN/m³).If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes.In a coordinate system, one individual stowage lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles on coordinate system, that is, volume load can resolve into three components in Descartes's rectangular coordinate system), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is not all identical, the component that lotus so is storaged in this several body in the method regards same amount of independent load as, the amplitude changing ratio identical volume sharing part of the load of wherein distribution intensity can also be synthesized to an individual stowage lotus or it is load.

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

When load is embodied as distributed load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, similar, " load change " is referred specifically to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " is referred specifically to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " is referred specifically to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " is referred specifically to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that the amplitude of distribution intensity changes ", briefly, now " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so distributed load so-and-so change ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.

When load is embodied as volume load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of volume load ", similar, " load change " refers to " change of the amplitude of the distribution intensity of volume load ", " load change amount " refers to " variable quantity of the amplitude of the distribution intensity of volume load ", " load change degree " refers to " intensity of variation of the amplitude of the distribution intensity of volume load ", " the actual change amount of load " refers to " the actual change amount of the amplitude of the distribution intensity of volume load ", " load changed " refers to " the volume load that the amplitude of distribution intensity changes ", briefly, " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.

This method is specifically included：

A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable；The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation；This method unitedly calls evaluated support cable and load is " evaluation object ", if the evaluated quantity of support cable and the quantity sum of load is that N, the i.e. quantity of " evaluation object " are N；This method refers exclusively to the evaluated support cable in " evaluation object " with title " core evaluation object ", and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step；This method represents this numbering, k=1,2,3 ..., N with variable k；The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, it is clear that the quantity of core evaluation object is exactly Q；The monitored rope force data of Cable Structure M in Cable Structure₁The M of individual specified support cable₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change；M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure；M₁It is an integer for being not more than Q not less than 0；Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point；M is had every time₂Individual strain measurement value or calculated value strain to characterize Cable Structure, M₂For K₂And L₂Product；M₂It is no less than 0 integer；Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified；M is had every time₃Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M₃For K₃、L₃And H₃Product；M₃It is an integer not less than 0；The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point；M is had every time₄Individual coordinates measurements or calculated value characterize Cable Structure shape, M₄For K₄And L₄Product；M₄It is an integer not less than 0；The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M₁、M₂、M₃And M₄Sum, it is M to define parameter K, K₁、K₂、K₃And K₄Sum, M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object；For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method；Time interval between any measurement twice monitored in real time to same amount in this method cannot be greater than 30 minutes, be referred to as the physical record data moment at the time of measurement record data；The external force that object, structure are born can be described as load, and load includes face load and volume load；Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load；Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force；Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component；Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant；Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load；Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes；Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method；

B. this method definition " temperature survey of the Cable Structure of this method calculates method " is carried out by step b1 to b3；

b1：Inquiry or actual measurement obtain the thermal conduction study parameter varied with temperature of Cable Structure composition material and Cable Structure local environment, using the geometry measured data of the design drawing, as-built drawing and Cable Structure of Cable Structure, the Thermodynamic calculation model of Cable Structure is set up using these data and parameter；Inquire about the meteorological data in recent years that Cable Structure location is no less than 2 years, count the cloudy quantity obtained in this period and be designated as T cloudy day, it will can not see the one of the sun daytime in the method and be referred to as the cloudy day all day, obtain each cloudy day in T cloudy day 0 is counted up to the highest temperature and the lowest temperature after sunrise moment next day between 30 minutes, the sunrise moment referred to according to the sunrise moment on earth rotation and the meteorology of revolution rule determination, do not indicate that the same day necessarily can see that the sun, data can be inquired about or calculated by conventional meteorology and obtain the required sunrise moment of each day, the 0 of each cloudy day subtracts the maximum temperature difference that the lowest temperature is referred to as the cloudy daily temperature up to the highest temperature after sunrise moment next day between 30 minutes, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximum in the maximum temperature difference of T cloudy daily temperature is taken to refer to temperature difference per day, Δ T is designated as with reference to temperature difference per day_r；Inquiry Cable Structure location and meteorological data in recent years of the height above sea level interval in place no less than 2 years or actual measurement obtain the temperature of Cable Structure local environment with time and the delta data and changing rule of height above sea level, calculate maximum rate of change Δ T of the temperature for obtaining the Cable Structure local environment in recent years of Cable Structure location and place height above sea level interval no less than 2 years on height above sea level_h, Δ T is taken for convenience of narration_hUnit for DEG C/m；Taken on the surface of Cable Structure " R Cable Structure surface point ", the Specific Principles of " R Cable Structure surface point " are taken to be described in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, it is called " R Cable Structure surface temperature measured data " to survey obtained temperature data, if utilizing the Thermodynamic calculation model of Cable Structure, the temperature of this R Cable Structure surface point is obtained by Calculation of Heat Transfer, it is called " R Cable Structure land surface pyrometer count evidence " just to calculate obtained temperature data；From the minimum height above sea level residing for Cable Structure to highest height above sea level,It is uniform in Cable Structure to choose no less than three different height above sea levels,At the height above sea level of each selection,Two points are at least chosen at the intersection on horizontal plane Yu Cable Structure surface,The exterior normal of straw line body structure surface at selected point,The exterior normal direction of all selections is referred to as in " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness ",Direction of the Cable Structure along the Temperature Distribution of wall thickness is measured with " horizontal plane and the intersection on Cable Structure surface " to intersect,The sunny slope exterior normal direction of Cable Structure and the in the shade face exterior normal direction of Cable Structure must be included in direction of the measurement Cable Structure along the Temperature Distribution of wall thickness of selection,No less than three points are chosen along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness is uniform in Cable Structure,For support cable a point is only taken along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness,Only measure the temperature of the surface point of support cable,All temperature being selected a little of measurement,The temperature measured is referred to as " temperature profile data of the Cable Structure along thickness ",Wherein edge is intersected with same " horizontal plane and the intersection on Cable Structure surface "," temperature profile data of the Cable Structure along thickness " that " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness " measurement is obtained,It is referred to as in the method " temperature profile data of the identical height above sea level Cable Structure along thickness ",If have chosen H different height above sea levels,At each height above sea level,It has chosen direction of the B measurement Cable Structure along the Temperature Distribution of wall thickness,Direction along each measurement Cable Structure along the Temperature Distribution of wall thickness have chosen E point in Cable Structure,Wherein H and E are not less than 3,B is not less than 2,It is equal to 1 for support cable E,The sum for counting in Cable Structure " point of the measurement Cable Structure along the temperature profile data of thickness " is HBE,The temperature of this HBE " point of the measurement Cable Structure along the temperature profile data of thickness " will be obtained by actual measurement below,It is called " HBE Cable Structure along thickness temperature measured data " to survey obtained temperature data,If utilizing the Thermodynamic calculation model of Cable Structure,Temperature of this HBE measurement Cable Structure along the point of the temperature profile data of thickness is obtained by Calculation of Heat Transfer,It is called " HBE Cable Structure calculates data along thickness temperature " just to calculate obtained temperature data；Require to choose a position according to meteorology measurement temperature in Cable Structure location, the temperature of environment where meeting the Cable Structure of meteorology measurement temperature requirement will be obtained in the actual measurement of this position；A position is chosen at the spacious unobstructed place in Cable Structure location, the position should it is annual can obtain each day the ground this getable day most sufficient sunshine, in the flat board of one piece of carbon steel material of position of sound production, referred to as reference plate, reference plate not can contact with ground, distance is not less than 1.5 meters to reference plate from the ground, the one side of the reference plate faces south, referred to as sunny slope, the sunny slope of reference plate is coarse and dark, the sunny slope of reference plate should it is annual can obtain each day one flat plate the ground this getable day most sufficient sunshine, 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：Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point in real time, monitoring obtains temperature profile data of the previously defined Cable Structure along thickness in real time simultaneously, while monitoring in real time obtains the temperature record of environment where meeting the Cable Structure of meteorology measurement temperature requirement；By monitor in real time obtain the Cable Structure that the same day is carved into after sunrise moment next day between 30 minutes at sunrise where environment temperature measured data sequence, the temperature measured data of environment is arranged according to time order and function order where the temperature measured data sequence of environment was carved into after the sunrise moment next day Cable Structure between 30 minutes at sunrise by the same day where Cable Structure, maximum temperature and minimum temperature in the temperature measured data sequence of environment where finding Cable Structure, the same day for subtracting environment where minimum temperature obtains Cable Structure with the maximum temperature in the temperature measured data sequence of environment where Cable Structure is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as environment maximum temperature difference, it is designated as Δ T_emax；Rate of change of the temperature of environment where obtaining Cable Structure on the time is calculated by Conventional mathematical by the temperature measured data sequence of environment where Cable Structure, the rate of change is also with time change；By the measured data sequence for monitoring the temperature for obtaining the sunny slope that the same day is carved into reference plate after sunrise moment next day between 30 minutes at sunrise in real time, the measured data that the measured data sequence of the temperature of the sunny slope of reference plate was carved into the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by the same day at sunrise is arranged according to time order and function order, maximum temperature and minimum temperature in the measured data sequence for the temperature for finding the sunny slope of reference plate, same day of temperature that the sunny slope that minimum temperature obtains reference plate is subtracted with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate is carved into maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as reference plate maximum temperature difference, it is designated as Δ T_pmax；The Cable Structure surface temperature measured data sequence for all R Cable Structure surface points that the same day is carved into after sunrise moment next day between 30 minutes at sunrise is obtained by monitoring in real time, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, the Cable Structure surface temperature measured data that each Cable Structure surface temperature measured data sequence was carved into after sunrise moment next day between 30 minutes by the same day of a Cable Structure surface point at sunrise is arranged according to time order and function order, find the maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, the same day for subtracting the temperature that minimum temperature obtains each Cable Structure surface point with the maximum temperature in each Cable Structure surface temperature measured data sequence is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, there is R Cable Structure surface point just to there is the R same day to be carved into the maximum temperature difference numerical value after sunrise moment next day between 30 minutes at sunrise, maximum therein is referred to as Cable Structure surface maximum temperature difference, it is designated as Δ T_smax；Calculated by each Cable Structure surface temperature measured data sequence by Conventional mathematical and obtain rate of change of the temperature on the time of each Cable Structure surface point, the temperature of each Cable Structure surface point on the time rate of change also with time change；Obtain the same day by monitoring in real time and be carved at sunrise after sunrise moment next day between 30 minutes, in synchronization, after HBE " temperature profile data of the Cable Structure along thickness ", calculate the difference of the maximum temperature and minimum temperature that amount at the height above sea level of each selection in BE " temperature profile data of the identical height above sea level Cable Structure along thickness ", the absolute value of this difference is referred to as " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", have chosen H different height above sea levels just has H " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", maximum in this H " Cable Structure thickness direction maximum temperature difference at identical height above sea level " is called " Cable Structure thickness direction maximum temperature difference ", it is designated as Δ T_tmax；

b3：Survey calculation obtains Cable Structure steady temperature data；First, it is determined that at the time of obtaining Cable Structure steady temperature data, the condition related at the time of Cable Structure steady temperature data to determining to obtain has six, Section 1 condition was carved into after sunrise moment next day between 30 minutes at sunset at the time of being and obtain Cable Structure steady temperature data between the same day, the sunset moment refers to that, according to the sunset moment on earth rotation and the meteorology of revolution rule determination, data can be inquired about or calculate by conventional meteorology obtaining the required sunset moment of each day；The a conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, reference plate maximum temperature difference Δ T_pmaxWith Cable Structure surface maximum temperature difference Δ T_smaxAll it is not more than 5 degrees Celsius；The b conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, the environment maximum temperature difference Δ T obtained in above survey calculation_emaxNo more than refer to temperature difference per day Δ T_r, and reference plate maximum temperature difference Δ T_pmaxSubtract and be not more than Δ T after 2 degrees Celsius_emax, and Cable Structure surface maximum temperature difference Δ T_smaxNo more than Δ T_pmax；Need to only meet in a conditions and b conditions of Section 2 one is known as meeting Section 2 condition；Section 3 condition is that the temperature of environment is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time where Cable Structure at the time of Cable Structure steady temperature data are obtained；Section 4 condition is that at the time of Cable Structure steady temperature data are obtained, the temperature of each Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time；Section 5 condition is that at the time of Cable Structure steady temperature data are obtained, the Cable Structure surface temperature measured data of each Cable Structure surface point in R Cable Structure surface point was carved into the minimum after sunrise moment next day between 30 minutes for the same day at sunrise；Section 6 condition is " Cable Structure thickness direction maximum temperature difference " Δ T at the time of Cable Structure steady temperature data are obtained_tmaxNo more than 1 degree Celsius；This method utilizes above-mentioned six conditions, any one in the following three moment is referred to as " the mathematics moment for obtaining Cable Structure steady temperature data ", the first moment is at the time of meeting Section 1 to the Section 5 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", second of moment is at the time of only meeting the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", the third moment is while at the time of meeting Section 1 to the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain "；It it is exactly the mathematics moment for obtaining Cable Structure steady temperature data at the time of acquisition Cable Structure steady temperature data when a moment in the physical record data moment during the mathematics moment for obtaining Cable Structure steady temperature data is exactly this method；If the mathematics moment for obtaining Cable Structure steady temperature data is not any one moment in the physical record data moment in this method, take this method closest to the mathematics moment for obtaining Cable Structure steady temperature data that physical record data at the time of to obtain at the time of Cable Structure steady temperature data；The amount that this method is used in measurement record at the time of obtaining Cable Structure steady temperature data carries out the monitoring analysis of Cable Structure relevant health；The Cable Structure temperature field that this method is approximately considered at the time of obtaining Cable Structure steady temperature data is in the Cable Structure temperature at stable state, i.e. this moment and not changed over time, and this moment is exactly " at the time of the obtaining Cable Structure steady temperature data " of this method；Then, according to Cable Structure heat-transfer character, utilize " the R Cable Structure surface temperature measured data " at the time of obtaining Cable Structure steady temperature data and " HBE Cable Structure along thickness temperature measured data ", utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution for obtaining Cable Structure at the time of Cable Structure steady temperature data are obtained is calculated by conventional heat transfer, now the temperature field of Cable Structure is calculated by stable state, calculating the temperature profile data of obtained Cable Structure at the time of Cable Structure steady temperature data are obtained includes the calculating temperature of R Cable Structure surface point in Cable Structure, the calculating temperature of R Cable Structure surface point is referred to as R Cable Structure steady-state surface temperature and calculates data, also include calculating temperature of the Cable Structure above selected HBE " point of the measurement Cable Structure along the temperature profile data of thickness ", the calculating temperature of HBE " points of the measurement Cable Structure along the temperature profile data of thickness " is referred to as " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data calculates data correspondent equal with R Cable Structure steady-state surface temperature, and " HBE Cable Structure along thickness temperature measured data " with " HBE Cable Structure along thickness temperature calculating data " correspondent equal when, the temperature profile data for calculating obtained Cable Structure at the time of Cable Structure steady temperature data are obtained is referred to as " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " now is referred to as " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure along thickness temperature measured data " is referred to as " HBE Cable Structure along thickness steady temperature measured data "；When " R Cable Structure surface point " is taken on the surface of Cable Structure, the quantity of " R Cable Structure surface point " must is fulfilled for three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the observed temperature linear interpolation of point adjacent with the arbitrfary point on Cable Structure surface during the temperature at any point on Cable Structure surface is by " R Cable Structure surface point " is obtained, the error of the temperature of the arbitrfary point and the actual temperature of the arbitrfary point on Cable Structure surface is not more than 5% on the Cable Structure surface that linear interpolation is obtained；Cable Structure surface includes support cable surface；Second condition is that the point being not less than in 4, and " R Cable Structure surface point " in same height above sea level in the quantity of the point of same height above sea level in " R Cable Structure surface point " is uniform along Cable Structure surface；Maximum Δ h in the absolute value of the difference of the height above sea level of " R Cable Structure surface point " along all Cable Structure surface points adjacent two-by-two of height above sea level is not more than 0.2 DEG C divided by Δ T_hObtained numerical value, Δ T is taken for convenience of narration_hUnit for DEG C/m, for convenience of narration take Δ h unit be m；When the definition of " R Cable Structure surface point " along the Cable Structure surface point adjacent two-by-two of height above sea level refers to only consider height above sea level, a Cable Structure surface point is not present in " R Cable Structure surface point ", the height above sea level numerical value of the Cable Structure surface point is between the height above sea level numerical value of adjacent Cable Structure surface point two-by-two；3rd condition is inquiry or obtains Cable Structure location and the interval sunshine rule of place height above sea level by meteorology conventionally calculation, further according to the geometric properties and bearing data of Cable Structure, found in Cable Structure it is annual by the sunshine-duration most sufficient position of those surface points, in " R Cable Structure surface point " at least one Cable Structure surface point be in Cable Structure whole year by a point in those most sufficient surface points of sunshine-duration；

C. the Cable Structure steady temperature data obtained under original state are calculated according to " temperature survey of the Cable Structure of this method calculates method " direct measurement, Cable Structure steady temperature data under original state are referred to as initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T_o”；Survey or consult reference materials and obtain the physical and mechanical properties parameter varied with temperature of various materials used in Cable Structure；Initial Cable Structure steady temperature data vector T is obtained in actual measurement_oSynchronization, direct measurement, which is calculated, obtains the Initial cable forces of all support cables, composition Initial cable force vector F_o；According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector T_oUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cable_oElement coding rule it is identical；T is obtained in actual measurement_oWhile, that is, obtaining initial Cable Structure steady temperature data vector T_oAt the time of synchronization, direct measurement, which is calculated, obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, initial Cable Structure bearing generalized coordinates data include initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, and the data of the health status that can express support cable now are referred to as support cable initial health data；The monitored amount initial value vector C of initial value composition of all monitored amounts_o, it is monitored amount initial value vector C_oCoding rule and M monitored amounts coding rules it is identical；Evaluation object initial damage vector d is set up using support cable initial health data and Cable Structure load measurement data_o, vectorial d_oRepresent with initial mechanical calculating benchmark model A_oThe initial health of the evaluation object of the Cable Structure of expression；Evaluation object initial damage vector d_oElement number be equal to N, d_oElement and evaluation object be one-to-one relationship, vectorial d_oElement coding rule it is identical with the coding rule of evaluation object；If d_oThe corresponding evaluation object of some element be a support cable in cable system, then d_oThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion；If d_oThe corresponding evaluation object of some element be to take d in some load, this method_oThe element numerical value be 0, the initial value for representing the change of this load is 0；If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0；Initial Cable Structure bearing generalized coordinates data constitute initial Cable Structure bearing generalized coordinates vector U_o；

D. the physical and mechanical properties parameter varied with temperature of various materials, initial Cable Structure bearing generalized coordinates vector U according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure_o, initial Cable Structure steady temperature data vector T_oAll Cable Structure data obtained with preceding step, set up the initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data "_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot 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_oCable Structure bearing generalized coordinates data be exactly initial Cable Structure bearing generalized coordinates vector U_o；Corresponding to A_oEvaluation object health status evaluation object initial damage vector d_oRepresent；Corresponding to A_oAll monitored amounts initial value with monitored amount initial value vector C_oRepresent；U_o、T_oAnd d_oIt is A_oParameter, by A_oThe obtained initial value of all monitored amounts of Mechanics Calculation result and C_oThe initial value of all monitored amounts represented is identical, therefore alternatively C_oBy A_oMechanics Calculation result composition, A in the method_o、C_o、d_o、U_oAnd T_oIt is constant；

E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o, A_oAnd Aⁱ _oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Tⁱ _oRepresent, vector Tⁱ _oDefinition mode and vector T_oDefinition mode it is identical, Tⁱ _oElement and T_oElement correspond；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure bearing generalized coordinates data " with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent, vectorial Uⁱ _oDefinition mode and vector U_oDefinition mode it is identical, Uⁱ _oElement and U_oElement correspond；The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；Uⁱ _o、Tⁱ _oAnd dⁱ _oIt is Aⁱ _oCharacterisitic parameter, Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, Tⁱ _oIt is designated as T¹ _o, set up T¹ _oMethod to make T¹ _oEqual to T_o；When circulation starts for the first time, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o；

F. enter from here and the circulation walked to q is walked by f；During structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Tⁱ, vector TⁱDefinition mode and vector T_oDefinition mode it is identical, TⁱElement and T_oElement correspond；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, actual measurement obtains Cable Structure bearing generalized coordinates current data, and all Cable Structure bearing generalized coordinates current datas constitute current Cable Structure actual measurement bearing generalized coordinates vector Uⁱ, vectorial UⁱDefinition mode and vector U_oDefinition mode it is identical；Vector T is obtained in actual measurementⁱWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TⁱAt the time of synchronization Cable Structure in all monitored amounts currency, the monitored amount current value vector C of all these numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, actual measurement obtain all M in Cable Structure₁The rope force data of root support cable, all these rope force data composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, Actual measurement obtains all M₁The space coordinate of two supporting end points of root support cable, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vector_oElement coding rule it is identical；

G. bearing generalized coordinates vector U is surveyed according to current Cable StructureⁱWith current Cable Structure steady temperature data vector Tⁱ, current initial mechanical calculating benchmark model A is updated according to step g1 to g3ⁱ _o, the monitored current initial value vector C of amountⁱ _o, current initial Cable Structure steady temperature data vector Tⁱ _oWith current initial Cable Structure bearing generalized coordinates vector Uⁱ _o, and the current initial damage vector d of evaluation objectⁱ _oKeep constant；

G1. it is respectively compared TⁱAnd Tⁱ _o、UⁱAnd Uⁱ _oIf, TⁱEqual to Tⁱ _oAnd UⁱEqual to Uⁱ _o, then need not be to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated, otherwise needs to follow these steps to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated；

G2. U is calculatedⁱWith U_oDifference, UⁱWith U_oDifference be exactly generalized displacement of support of the Cable Structure bearing on initial position, represent generalized displacement of support with generalized displacement of support vector V, V is equal to UⁱSubtract U_o；Calculate TⁱWith T_oDifference, TⁱWith T_oDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TⁱWith T_oDifference represented with steady temperature change vector S, S be equal to TⁱSubtract T_o, S represents the change of Cable Structure steady temperature data；

G3. first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then to A_oIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, to A_oMiddle Cable Structure bearing applies generalized displacement of support constraint and applies the current initial mechanical calculating benchmark model A updated after temperature change to Cable Structureⁱ _o, update Aⁱ _oWhile, Uⁱ _oAll elements numerical value also uses UⁱAll elements numerical value correspondence is replaced, that is, have updated Uⁱ _o, Tⁱ _oAll elements numerical value also uses TⁱAll elements numerical value correspondence replace, that is, have updated Tⁱ _o, thus obtained properly corresponding to Aⁱ _oUⁱ _oAnd Tⁱ _o, now dⁱ _oKeep constant；As renewal Aⁱ _oAfterwards, Aⁱ _oRope the health status current initial damage vector d of evaluation objectⁱ _oRepresent, Aⁱ _oCable Structure steady temperature with current Cable Structure steady temperature data vector Tⁱ _oRepresent, Aⁱ _oBearing generalized coordinates with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent；Update Cⁱ _oMethod be：As renewal Aⁱ _oAfterwards, A is obtained by Mechanics Calculationⁱ _oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Cⁱ _o；

H. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u；

H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4ⁱAnd Dⁱ _u；At other moment, when in step g to Aⁱ _oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CⁱAnd Dⁱ _uIf, not to A in step gⁱ _oIt is updated, then is directly transferred to step i here and carries out follow-up work；

H2. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating；According to the coding rule of evaluation object, calculated successively；Calculate each time and assume that only one of which evaluation object is further added by unit damage or load unit change on the basis of original damage or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is a load, it is assumed that the load is further added by load unit change, D is usedⁱ _ukThis increased unit damage or load unit change is recorded, wherein k represents the numbering for increasing unit damage or the evaluation object of load unit change, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uAn element, evaluation object unit change vector Dⁱ _uElement coding rule and vector d_oElement coding rule it is identical；The evaluation object that unit damage or load unit change are further added by calculating each time is different from the evaluation object that unit damage or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector；When assuming that k-th of evaluation object is further added by unit damage or load unit change, C is usedⁱ _tkRepresent corresponding " monitored amount calculation current vector "；When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed；Cⁱ _tkDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _tkElement and C_oElement correspond；

H3. obtained vectorial C is calculated each timeⁱ _tkSubtract vectorial Cⁱ _oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or load unit change numerical valueⁱ _k”；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount "；

H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "ⁱ”；Unit damage monitored numerical quantity transformation matrices Δ CⁱEach row correspond to a monitored amount unit change vector；Unit damage monitored numerical quantity transformation matrices Δ CⁱEvery a line correspond to same monitored amount different evaluation objects increase unit damage or load unit change when different unit change amplitudes；Unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and vector d_oElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and M monitored amounts coding rules it is identical；

I. current nominal fatigue vector d is definedⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cElement numerical value represent correspondence evaluation object nominal fatigue degree or nominal load variable quantity, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship；

J. according to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1ⁱ _cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector dⁱ _c；

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

K. the currently practical injury vector d expressed using formula 2ⁱK-th of element dⁱ _kWith the current initial damage vector d of evaluation objectⁱ _oK-th of element dⁱ _okWith current nominal fatigue vector dⁱ _cK-th of element dⁱ _ckBetween relation, calculating obtain currently practical injury vector dⁱAll elements；

K=1,2,3 ... ..., N in formula 2；Vectorial dⁱElement coding rule and formula (1) in vector d_oElement coding rule it is identical；dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, dⁱ _kThe numerical response degree of relaxation or the damage of the support cable；By the currently practical injury vector d of evaluation objectⁱIn with M₁The related M of root support cable₁Individual element takes out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical；The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., M₁；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line；Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity；Using in current Cable Structure steady temperature data vector TⁱUnder the conditions of, in l walk the slack line that identifies and with the currently practical injury vector d of support cable^ciExpression these slack lines, the degree that relaxed with it mechanic equivalent currently practical equivalent damage degree, using f walk acquisition in current Cable Structure steady temperature data vector TⁱUnder the conditions of current cable force vector FⁱWith current support cable two supporting end points horizontal range vector, using c walk obtain in initial Cable Structure steady temperature data vector T_oUnder the conditions of support cable initial drift is vectorial, the weight of the initial free unit length of initial free cross-sectional area vector sum is vectorial, Initial cable force vector F_o, utilize current Cable Structure steady temperature data vector TⁱThe support cable current steady state temperature data of expression, using c walk obtain in initial Cable Structure steady temperature data vector T_oThe support cable initial steady state temperature data of expression, utilize the physical and mechanical properties parameter varied with temperature that various materials used in the Cable Structure of acquisition are walked in c, it is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is：First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical；2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope；When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa；Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined；So it is achieved that relaxation identification and the non-destructive tests of support cable；Institute's demand power is by current cable force vector F during calculatingⁱCorresponding element is provided；Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far this method realizes the problem of rejecting generalized displacement of support, load change and the influence of structure temperature change, Cable Structure rope and recognized, generalized displacement of support, structure temperature change and the change influence of support cable health status, load change amount identification are rejected while realizing；

L. current nominal fatigue vector d is tried to achieveⁱ _cAfterwards, mark vector B is set up according to formula 3ⁱ, formula 4 gives mark vector BⁱK-th of element definition；

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

Element B in formula 4ⁱ _kIt is mark vector BⁱK-th of element, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element, dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4；

If m. mark vector BⁱElement be all 0, then return to step f continue this circulation；If mark vector BⁱElement be not all 0, then enter next step, i.e. step n；

N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5ⁱ⁺¹ _oEach element；

D in formula 5ⁱ⁺¹ _okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationⁱ⁺¹ _oK-th of element, dⁱ _okThis, i.e. the current initial damage vector d of evaluation object of ith circulationⁱ _oK-th of element, Dⁱ _ukIt is the evaluation object unit change vector D of ith circulationⁱ _uK-th of element, Bⁱ _kIt is the mark vector B of ith circulationⁱK-th of element, k=1,2,3 ... ..., N in formula 5；

O. in initial mechanical calculating benchmark model A_oOn the basis of, first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then to A_oIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, then it is d to make the health status of ropeⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationⁱ⁺¹；Obtain Aⁱ⁺¹Afterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o；

P. remove once, i.e. i+1 time circulates required current initial Cable Structure steady temperature data vector Tⁱ⁺¹ _oThe current initial Cable Structure steady temperature data vector T circulated equal to ithⁱ _o；Next time, the current initial Cable Structure bearing generalized coordinates vector U i.e. needed for i+1 time circulationⁱ⁺¹ _oThe current initial Cable Structure bearing generalized coordinates vector U circulated equal to ithⁱ _o；

Q. step f is returned to, starts to circulate next time.

Beneficial effect：Structural healthy monitoring system is monitored on-line for a long time by using sensor to structural response first, (or offline) analysis online is carried out after acquisition Monitoring Data to it and obtains structural health conditions data, due to the complexity of structure, structural healthy monitoring system needs to use the equipment such as substantial amounts of sensor to carry out monitoring structural health conditions, therefore its cost is generally at a relatively high, therefore cost problem is a subject matter of limit structural health monitoring technique application.On the other hand, the correct identification of the health status of core evaluation object (such as suspension cable) is the indispensable part of the correct identification of structural health conditions, even its whole, and influence of the correct identification of the change of secondary evaluation object (load that such as structure is born) (such as by the change of the quality and quantity of the automobile of cable-stayed bridge) to the correct identification of the health status of Cable Structure be it is very little, it is even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is typically suitable, the quantity of secondary evaluation object is also frequently more than the quantity of core evaluation object, and the quantity of such evaluation object is often many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, the quantity of the monitored amount of conventional method requirement (being obtained using device measurings such as sensors) have to be larger than the quantity equal to evaluation object, when the secondary evaluation object changed quantity than it is larger when (practically always such), the quantity of the equipment such as the sensor required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, or even unacceptablely high.Inventor's research is found,In secondary evaluation object (such as normal load that structure is born,The normal load of structure refers to that the load that structure bearing is no more than the structure allowable load limited according to structure design book or structure completion book) change it is smaller when (be exactly that structure is only subjected only to normal load for load,Whether the load that structure is born is normal load,It can be determined by the observation of the methods such as naked eyes,If it find that the load that structure is born is not normal load,It is so artificial to remove,Remove after improper load,Structure is just solely subjected to normal load),Amplitude of variation of the amplitude of variation (this specification is called " secondary response ") of structural response caused by them much smaller than the structural response caused by the change (such as support cable is damaged) of core evaluation object (this specification is called " core response "),Secondary response responds total change that sum is structural response with core (this specification is called " global response "),Obvious core response occupies leading position in global response,Based on this,Even if inventor's research hair is chosen when now determining that monitored amount quantity is slightly larger than core evaluation object quantity,But much smaller than the numerical value of evaluation object quantity (this method is exactly so to do),Even if that is using equipment such as the relatively few many sensors of quantity,Still the state of health data of core evaluation object can accurately be obtained,Meet the core demand of structural health conditions monitoring,Therefore cost of the cost of the structural healthy monitoring system proposed by this method apparently than the structural healthy monitoring system required by conventional method is much lower,That is this method can realize the assessment of the health status to the core evaluation object of Cable Structure with the much lower condition of cost,This benefit is that can structural health monitoring technology be used is very important.

Embodiment

This method uses a kind of algorithm, and the algorithm is used for the health status for recognizing core evaluation object.When it is implemented, the following steps are one kind in the various steps that can be taken.

The first step：First confirm that the quantity for the load that the possibility that Cable Structure is born changes.The characteristics of load born according to Cable Structure, confirmation wherein " is possible to the load changed ", or all load is considered as " being possible to the load changed ", if the shared JZW load that may be changed, that is, have JZW secondary evaluation objects.If the quantity sum of the quantity and JZW " being possible to the load changed " of the support cable of Cable Structure is N, that is, have N number of evaluation object.Evaluation object serial number is given, the numbering will be used to generate vector sum matrix in subsequent step.

Monitored multiclass parameter can include：Suo Li, strain, angle and space coordinate, are described below respectively：

If having Q root support cables in cable system, that is, have a Q core evaluation object, the monitored rope force data of the Cable Structure M in Cable Structure₁The M of individual specified rope₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified ropes change.M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M₁It is an integer not less than 0.

The monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point.M is had every time₂(M₂=K₂×L₂) individual strain measurement value or calculated value characterize Cable Structure strain.M₂It is an integer not less than 0.

The monitored angle-data of Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified.M is had every time₃(M₃=K₃×L₃×H₃) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M₃It is an integer not less than 0.

The monitored shape data of Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point.M is had every time₄(M₄=K₄×L₄) individual coordinates measurements or calculated value characterize Cable Structure shape.M₄It is an integer not less than 0.

The monitored amount of summary, whole Cable Structure has M (M=M₁+M₂+M₃+M₄) individual monitored amount, define parameter K (K=M₁+K₂+K₃+K₄), the quantity that M cannot be less than core evaluation object plus 4, and M is less than the quantity N of evaluation object.

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

Provide that step is determined " temperature survey of the Cable Structure of this method calculates method " by technical scheme.

Second step：Set up initial mechanical calculating benchmark model A_o。

When Cable Structure is completed, or before health monitoring systems are set up, obtain " Cable Structure steady temperature data " according to " temperature survey of the Cable Structure of this method calculates method " survey calculation (can be measured with ordinary temperature measuring method, for example measured using thermal resistance), " Cable Structure steady temperature data " now use vector T_oRepresent, referred to as initial Cable Structure steady temperature data vector T_o.T is obtained in actual measurement_oWhile, the initial value for all monitored amounts for obtaining Cable Structure, the monitored amount initial value vector C of composition are calculated using conventional method direct measurement_o。

Synchronization that can specifically in following manner at the time of so-and-so (such as initial or current) Cable Structure steady temperature data vector is obtained in this method, the data of so-and-so the measured monitored amount (all monitored amounts of such as Cable Structure) of amount are obtained using so-and-so method survey calculation：While measurement record temperature (temperature, the temperature of the sunny slope of reference plate and the Cable Structure surface temperature that include environment where Cable Structure), for example every 10 minutes measurement temperature of record, then while equally also recording the data of so-and-so the measured monitored amount (such as all monitored amounts of Cable Structure) of amount every measurement in 10 minutes.Once it is determined that at the time of acquisition Cable Structure steady temperature data, the data of so-and-so the measured monitored amount (all monitored amounts of such as Cable Structure) of amount of synchronization are known as synchronization at the time of Cable Structure steady temperature data are obtained at the time of so with obtaining Cable Structure steady temperature data, the data of so-and-so the measured monitored amount of amount obtained using so-and-so method survey calculation method.

The physical parameter (such as thermal coefficient of expansion) varied with temperature and mechanical property parameters (such as modulus of elasticity, Poisson's ratio) of various materials used in Cable Structure are obtained using conventional method (consult reference materials or survey).

T is obtained in actual measurement_oWhile, direct measurement calculates the Initial cable force for obtaining all support cables, composition Initial cable force vector F_o；Obtain length of all support cables when free state i.e. Suo Li is 0 according to Cable Structure design data, completion data, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector T_oUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, the initial drift vector l of support cable is constituted successively_o, initial free cross-sectional area vector A_oWith the weight vector ω of initial free unit length_o, the initial drift vector l of support cable_o, initial free cross-sectional area vector A_oWith the weight vector ω of initial free unit length_oElement coding rule and Initial cable force vector F_oElement coding rule it is identical.

Step is provided by technical scheme, T is obtained in actual measurement_oWhile, the Actual measurement data of Cable Structure are obtained using conventional method Actual measurement.Utilize the physical and mechanical properties parameter varied with temperature of various materials used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure, initial Cable Structure bearing generalized coordinates vector U_oWith initial Cable Structure steady temperature data vector T_o, it is included in " Cable Structure steady temperature data " using mechanics method (such as FInite Element) and sets up initial mechanical calculating benchmark model A_o。T_o、U_oAnd d_oIt is A_oParameter, C_oBy A_oMechanics Calculation result composition.

3rd step：In the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o, A_oAnd Aⁱ _oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Tⁱ _oRepresent, vector Tⁱ _oDefinition mode and vector T_oDefinition mode it is identical, Tⁱ _oElement and T_oElement correspond；Current initial mechanical calculating benchmark model A that ith circulation needs when starting, corresponding to Cable Structureⁱ _oThe current initial Cable Structure bearing generalized coordinates vector U of Cable Structure bearing generalized coordinates data compositionⁱ _o, the current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timeⁱ _oWhen, Uⁱ _oIt is equal to U_o.The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；Tⁱ _oAnd dⁱ _oIt is Aⁱ _oCharacterisitic parameter；Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, Tⁱ _oIt is designated as T¹ _o, set up T¹ _oMethod to make T¹ _oEqual to T_o；When circulation starts for the first time, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o。

4th step：The hardware components of cable structure health monitoring system are installed.Hardware components at least include：Monitored amount monitoring system (such as subsystem containing angular surveying, cable force measurement subsystem, strain measurement subsystem, space coordinate measures subsystem, signal conditioner etc.), Cable Structure bearing generalized coordinates monitoring system (contains total powerstation, angle measuring sensor, signal conditioner etc.), Cable Structure temperature monitoring system (contains temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (contain temperature sensor, signal conditioner etc.), support cable cable force monitoring system, the space coordinate monitoring system of the supporting end points of support cable, signal (data) collector, computer and communication alert equipment.Each monitored amount, each bearing generalized coordinates of Cable Structure, each temperature, the Suo Li of each support cable, the space coordinate of the supporting end points of each support cable must be monitored system monitoring and arrive, and the signal monitored is transferred to signal (data) collector by monitoring system；Signal is delivered to computer through signal picker；Computer is then responsible for the health monitoring software of the evaluation object of operation Cable Structure, including the signal that the transmission of tracer signal collector comes；When monitoring that evaluation object health status is changed, computer control communication warning device is alarmed to monitoring personnel, owner and (or) the personnel specified.

5th step：Establishment and the on computers system software of installation and operation this method, the software will complete the functions (all work that can be completed with computer i.e. in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, alarm that this method required by task is wanted.

6th step：Thus step starts the cycle over running, during structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Tⁱ, vector TⁱDefinition mode and vector T_oDefinition mode it is identical, TⁱElement and T_oElement correspond；In actual measurement vector TⁱWhile, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time.

In actual measurement vector TⁱWhile, actual measurement obtains Cable Structure bearing generalized coordinates current data, all data composition current cable structure actual measurement bearing generalized coordinates vector Uⁱ。

In actual measurement vector TⁱWhile, actual measurement obtains all M in Cable Structure₁The rope force data of root support cable, all these rope force data composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；In actual measurement vector TⁱWhile, Actual measurement obtains all M₁The space coordinate of two of root support cable supporting end points, the difference of the space coordinates of two supporting end points component in the horizontal direction is exactly two supporting end points horizontal ranges, all M₁Two supporting end points horizontal range data of root support cable constitute the current supporting end points horizontal range vector of support cable two lⁱ _x, the current supporting end points horizontal range vector of support cable two lⁱ _xElement coding rule and Initial cable force vector F_oElement coding rule it is identical.

7th step：Obtaining current Cable Structure actual measurement bearing generalized coordinates vector UⁱWith current Cable Structure steady temperature data vector TⁱAfterwards, it is respectively compared UⁱAnd Uⁱ _o、TⁱAnd Tⁱ _oIf, UⁱEqual to Uⁱ _oAnd TⁱEqual to Tⁱ _o, then need not be to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated, otherwise needs to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated, and the current initial damage vector d of evaluation objectⁱ _oKeep constant, update method is carried out by technical scheme regulation step.

8th step：Step is provided by technical scheme, in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u.If specifically, the evaluation object is a support cable in cable system, then it is assumed that the support cable is in vectorial dⁱ _oThe support cable represented has unit damage (for example taking 5%, 10%, 20% or 30% equivalent damage to be damaged for unit) again on the basis of having damaged, if the evaluation object is a load, it is assumed that the load is in vectorial dⁱ _oLoad unit is further added by the basis of the existing variable quantity of the load represented to change (if the load is distributed load, and the distributed load is line distributed load, load unit change can take 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. to change for unit；If the load is distributed load, and the distributed load is EDS maps load, and load unit change can take 1MPa, 2MPa, 3MPa or 1kNm/m²、2kNm/m²、3kNm/m²Change Deng for unit；If the load is concentrfated load, and the concentrfated load is couple, and load unit change can take 1kNm, 2kNm, 3kNm etc. to change for unit；If the load is concentrfated load, and the concentrfated load is concentrated force, and load unit change can take 1kN, 2kN, 3kN etc. to change for unit；If the load is volume load, load unit change can take 1kN/m³、2kN/m³、3kN/m³Change Deng for unit).

9th step：Set up linear relationship error vector eⁱWith vectorial gⁱ.Utilize data (" the monitored current initial value vector C of amount aboveⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ"); while the 8th step is calculated each time; i.e. while the increase unit damage for assuming only one of which evaluation object in evaluation object or load unit change is calculated each time; when assuming that kth (k=1,2,3; ...; N) individual evaluation object increase unit damage or during load unit change, one injury vector of composition is calculated each time, d is usedⁱ _tkThe injury vector is represented, corresponding monitored amount calculation current vector is Cⁱ _tk(referring to the 8th step), injury vector dⁱ _tkElement number be equal to evaluation object quantity, vectorial dⁱ _tkAll elements in only one of which element numerical value take each time calculate in assume increase unit damage or load unit change evaluation object unit damage or load unit changing value, dⁱ _tkThe numerical value of other elements take 0, that for 0 element evaluation object of the numbering with assuming increase unit damage or load unit change corresponding relation, with the elements of the same numberings of other vectors with the corresponding relation of the evaluation object be identical；dⁱ _tkWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _tkElement and d_oElement be one-to-one relationship.By Cⁱ _tk、Cⁱ _o、ΔCⁱ、dⁱ _tkBring formula (1) into, obtain a linear relationship error vector eⁱ _k, calculate obtain a linear relationship error vector e each timeⁱ _k；eⁱ _kSubscript k represent the individual evaluation object increase unit damage of kth (k=1,2,3 ..., N) or load unit change.There is N number of evaluation object just to have n times calculating, just there is N number of linear relationship error vector eⁱ _k, by this N number of linear relationship error vector eⁱ _kA vector is obtained after addition, is exactly final linear relationship error vector e by the new vector obtained after each element divided by N of this vectorⁱ.Vectorial gⁱEqual to final error vector eⁱ.By vectorial gⁱOn the hard disc of computer for being stored in operation health monitoring systems software, used for health monitoring systems software.

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

Tenth step：Define current nominal fatigue vector dⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cAnd dⁱElement numerical value represent correspondence evaluation object degree of injury or load change degree, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship.

11st step：According to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as formula (2), current nominal fatigue vector d calculated according to multi-objective optimization algorithmⁱ _cNoninferior solution, that is, position and its solution of nominal fatigue degree of damaged cable can be determined with reasonable error but relatively accurately from all ropes.

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

The Objective Programming (GoalAttainmentMethod) in multi-objective optimization algorithm can be used to solve formula (2) and obtain current nominal fatigue vector dⁱ _c。

12nd step：According to the currently practical injury vector d of cable systemⁱDefinition and its element definition calculate obtain currently practical injury vector dⁱEach element, so as to by dⁱDetermine the health status of evaluation object.Currently practical injury vector dⁱK-th of element dⁱ _kRepresent the currently practical health status of k-th of evaluation object in ith circulation.

dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kRepresent its currently practical damage, dⁱ _kFor 0 when represent its corresponding support cable without health problem, dⁱ _kNumerical value represents that its corresponding support cable is the support cable of unsoundness problem when not being 0, and the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, the degree of the relaxation of its numerical response or damage.

By the currently practical injury vector d of evaluation objectⁱIn the Q element related to support cable take out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical.The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., Q；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems by lossless detection method, remaining is exactly slack line, exactly needs to adjust Suo Li rope, these need to adjust Suo Li rope in the currently practical injury vector d of support cable^ciIn corresponding element numerical value (such as one of element can use d^ci _hRepresent) degree of injury with the relaxation degree mechanic equivalents of these support cables is represented, thus determine that slack line.Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity.

It is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent come calculate slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, the relaxation degree (i.e. the long adjustment amount of rope) of these ropes can be specifically tried to achieve according to formula (3).So it is achieved that the relaxation identification of support cable.So far damaged cable and slack line all just are identified.

${\mathrm{\Δl}}_h^i=\frac{d_h^{ci}}{1-d_h^{ci}}\frac{F_h^i}{\[\frac{E_h^i}{1+\frac{{\left({\mathrm{\ω}}_h^i{l}_{xh}^i\right)}^2{A}_h^i{E}_h^i}{12{\left(F_h^i\right)}^3}}\]A_h^i+F_h^i}{l}_{oh}^i---\left(3\right)$

E in formula (3)ⁱ _hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structureⁱ _oDuring expression, the modulus of elasticity of h root support cables, Aⁱ _hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structureⁱ _oDuring expression, the cross-sectional area of h root support cables, Fⁱ _hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structureⁱ _oDuring expression, the current cable power of h root support cables, d^ci _hIt is the currently practical degree of injury of h root support cables, ωⁱ _hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structureⁱ _oDuring expression, the weight of the unit length of h root support cables, lⁱ _xhIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structureⁱ _oDuring expression, the horizontal range of two supporting end points of h root support cables, lⁱ _xhIt is the current supporting end points horizontal range vector of support cable two lⁱ _xAn element, current support cable two supporting end points horizontal range vector lⁱ _xElement coding rule and initial drift vector l_oElement coding rule it is identical, Eⁱ _hIt can be obtained according to the characteristic material data for looking into or surveying h root support cables, Aⁱ _hAnd ωⁱ _hThermal coefficient of expansion, A that can be according to h root support cables_oh、ω_oh、Fⁱ _h、T_oAnd Tⁱ _oObtained by Typical physical and Mechanics Calculation.

13rd step：Computer in health monitoring systems is periodically automatic or generates cable system health condition form by human users' health monitoring systems.

14th step：Under specified requirements, the computer in health monitoring systems is automatically brought into operation communication alert equipment and alarmed to monitoring personnel, owner and (or) the personnel specified.

15th step：Set up mark vector BⁱIf, mark vector BⁱElement be all 0, then return to the 6th step and proceed health monitoring and calculating to cable system；If mark vector BⁱElement be not all 0, then complete after subsequent step, into circulating next time.

16th step：Calculating obtains initial damage vector d next time needed for (i.e. i+1 time, i=1,2,3,4 ...) circulationⁱ⁺¹ _oEach element dⁱ⁺¹ _ok(k=1,2,3 ..., N)；Second, in initial mechanical calculating benchmark model A_oOn the basis of, first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then to A_oIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, then it is d to make the health status of ropeⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time (i=1,2,3,4 ...) circulationⁱ⁺¹；Current initial Cable Structure steady temperature data vector T needed for (i.e. i+1 time, i=1,2,3,4 ...) is circulated next timeⁱ⁺¹ _oEqual to Tⁱ _o, the current initial Cable Structure bearing generalized coordinates vector U needed for (i.e. i+1 time, i=1,2,3,4 ...) is circulated next timeⁱ⁺¹ _oEqual to Uⁱ _o.Obtain Aⁱ⁺¹、dⁱ⁺¹ _o、Uⁱ⁺¹ _oAnd Tⁱ⁺¹ _oAfterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o。

17th step：The 6th step is returned to, starts the circulation by the 6th step to the 17th step.

Claims (1)

1. simplify generalized displacement hybrid monitoring problem cable load progressive-type recognition method, it is characterised in that methods described includes：

A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable；The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation；This method unitedly calls evaluated support cable and load is " evaluation object ", if the evaluated quantity of support cable and the quantity sum of load is that N, the i.e. quantity of " evaluation object " are N；This method refers exclusively to the evaluated support cable in " evaluation object " with title " core evaluation object ", and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object "；The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step；This method represents this numbering, k=1,2,3 ..., N with variable k；The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, it is clear that the quantity of core evaluation object is exactly Q；The monitored rope force data of Cable Structure M in Cable Structure₁The M of individual specified support cable₁Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change；M is had every time₁Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure；M₁It is an integer for being not more than Q not less than 0；Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K₂The L of individual specified point and each specified point₂The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K₂The change of all tested strains of individual specified point；M is had every time₂Individual strain measurement value or calculated value strain to characterize Cable Structure, M₂For K₂And L₂Product；M₂It is no less than 0 integer；Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure₃Individual specified point, excessively each specified point L₃The H of individual specified straight line, each specified straight line₃Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified；M is had every time₃Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M₃For K₃、L₃And H₃Product；M₃It is an integer not less than 0；The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure₄The L of individual specified point and each specified point₄The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K₄The change of all coordinate components of individual specified point；M is had every time₄Individual coordinates measurements or calculated value characterize Cable Structure shape, M₄For K₄And L₄Product；M₄It is an integer not less than 0；The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M₁、M₂、M₃And M₄Sum, it is M to define parameter K, K₁、K₂、K₃And K₄Sum, M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object；For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method；Time interval between any measurement twice monitored in real time to same amount in this method cannot be greater than 30 minutes, be referred to as the physical record data moment at the time of measurement record data；The external force that object, structure are born can be described as load, and load includes face load and volume load；Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load；Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force；Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component；Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant；Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load；Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude；If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes；Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method；

B. this method definition " temperature survey of the Cable Structure of this method calculates method " is carried out by step b1 to b3；

b1：Inquiry or actual measurement obtain the thermal conduction study parameter varied with temperature of Cable Structure composition material and Cable Structure local environment, using the geometry measured data of the design drawing, as-built drawing and Cable Structure of Cable Structure, the Thermodynamic calculation model of Cable Structure is set up using these data and parameter；Inquire about the meteorological data in recent years that Cable Structure location is no less than 2 years, count the cloudy quantity obtained in this period and be designated as T cloudy day, it will can not see the one of the sun daytime in the method and be referred to as the cloudy day all day, obtain each cloudy day in T cloudy day 0 is counted up to the highest temperature and the lowest temperature after sunrise moment next day between 30 minutes, the sunrise moment referred to according to the sunrise moment on earth rotation and the meteorology of revolution rule determination, do not indicate that the same day necessarily can see that the sun, data can be inquired about or calculated by conventional meteorology and obtain the required sunrise moment of each day, the 0 of each cloudy day subtracts the maximum temperature difference that the lowest temperature is referred to as the cloudy daily temperature up to the highest temperature after sunrise moment next day between 30 minutes, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximum in the maximum temperature difference of T cloudy daily temperature is taken to refer to temperature difference per day, Δ T is designated as with reference to temperature difference per day_r；Inquiry Cable Structure location and meteorological data in recent years of the height above sea level interval in place no less than 2 years or actual measurement obtain the temperature of Cable Structure local environment with time and the delta data and changing rule of height above sea level, calculate maximum rate of change Δ T of the temperature for obtaining the Cable Structure local environment in recent years of Cable Structure location and place height above sea level interval no less than 2 years on height above sea level_h, Δ T is taken for convenience of narration_hUnit for DEG C/m；Taken on the surface of Cable Structure " R Cable Structure surface point ", the Specific Principles of " R Cable Structure surface point " are taken to be described in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, it is called " R Cable Structure surface temperature measured data " to survey obtained temperature data, if utilizing the Thermodynamic calculation model of Cable Structure, the temperature of this R Cable Structure surface point is obtained by Calculation of Heat Transfer, it is called " R Cable Structure land surface pyrometer count evidence " just to calculate obtained temperature data；From the minimum height above sea level residing for Cable Structure to highest height above sea level,It is uniform in Cable Structure to choose no less than three different height above sea levels,At the height above sea level of each selection,Two points are at least chosen at the intersection on horizontal plane Yu Cable Structure surface,The exterior normal of straw line body structure surface at selected point,The exterior normal direction of all selections is referred to as in " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness ",Direction of the Cable Structure along the Temperature Distribution of wall thickness is measured with " horizontal plane and the intersection on Cable Structure surface " to intersect,The sunny slope exterior normal direction of Cable Structure and the in the shade face exterior normal direction of Cable Structure must be included in direction of the measurement Cable Structure along the Temperature Distribution of wall thickness of selection,No less than three points are chosen along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness is uniform in Cable Structure,For support cable a point is only taken along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness,Only measure the temperature of the surface point of support cable,All temperature being selected a little of measurement,The temperature measured is referred to as " temperature profile data of the Cable Structure along thickness ",Wherein edge is intersected with same " horizontal plane and the intersection on Cable Structure surface "," temperature profile data of the Cable Structure along thickness " that " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness " measurement is obtained,It is referred to as in the method " temperature profile data of the identical height above sea level Cable Structure along thickness ",If have chosen H different height above sea levels,At each height above sea level,It has chosen direction of the B measurement Cable Structure along the Temperature Distribution of wall thickness,Direction along each measurement Cable Structure along the Temperature Distribution of wall thickness have chosen E point in Cable Structure,Wherein H and E are not less than 3,B is not less than 2,It is equal to 1 for support cable E,The sum for counting in Cable Structure " point of the measurement Cable Structure along the temperature profile data of thickness " is HBE,The temperature of this HBE " point of the measurement Cable Structure along the temperature profile data of thickness " will be obtained by actual measurement below,It is called " HBE Cable Structure along thickness temperature measured data " to survey obtained temperature data,If utilizing the Thermodynamic calculation model of Cable Structure,Temperature of this HBE measurement Cable Structure along the point of the temperature profile data of thickness is obtained by Calculation of Heat Transfer,It is called " HBE Cable Structure calculates data along thickness temperature " just to calculate obtained temperature data；Require to choose a position according to meteorology measurement temperature in Cable Structure location, the temperature of environment where meeting the Cable Structure of meteorology measurement temperature requirement will be obtained in the actual measurement of this position；A position is chosen at the spacious unobstructed place in Cable Structure location, the position should it is annual can obtain each day the ground this getable day most sufficient sunshine, in the flat board of one piece of carbon steel material of position of sound production, referred to as reference plate, reference plate not can contact with ground, distance is not less than 1.5 meters to reference plate from the ground, the one side of the reference plate faces south, referred to as sunny slope, the sunny slope of reference plate is coarse and dark, the sunny slope of reference plate should it is annual can obtain each day one flat plate the ground this getable day most sufficient sunshine, 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：Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point in real time, monitoring obtains temperature profile data of the previously defined Cable Structure along thickness in real time simultaneously, while monitoring in real time obtains the temperature record of environment where meeting the Cable Structure of meteorology measurement temperature requirement；By monitor in real time obtain the Cable Structure that the same day is carved into after sunrise moment next day between 30 minutes at sunrise where environment temperature measured data sequence, the temperature measured data of environment is arranged according to time order and function order where the temperature measured data sequence of environment was carved into after the sunrise moment next day Cable Structure between 30 minutes at sunrise by the same day where Cable Structure, maximum temperature and minimum temperature in the temperature measured data sequence of environment where finding Cable Structure, the same day for subtracting environment where minimum temperature obtains Cable Structure with the maximum temperature in the temperature measured data sequence of environment where Cable Structure is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as environment maximum temperature difference, it is designated as Δ T_emax；Rate of change of the temperature of environment where obtaining Cable Structure on the time is calculated by Conventional mathematical by the temperature measured data sequence of environment where Cable Structure, the rate of change is also with time change；By the measured data sequence for monitoring the temperature for obtaining the sunny slope that the same day is carved into reference plate after sunrise moment next day between 30 minutes at sunrise in real time, the measured data that the measured data sequence of the temperature of the sunny slope of reference plate was carved into the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by the same day at sunrise is arranged according to time order and function order, maximum temperature and minimum temperature in the measured data sequence for the temperature for finding the sunny slope of reference plate, same day of temperature that the sunny slope that minimum temperature obtains reference plate is subtracted with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate is carved into maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as reference plate maximum temperature difference, it is designated as Δ T_pmax；The Cable Structure surface temperature measured data sequence for all R Cable Structure surface points that the same day is carved into after sunrise moment next day between 30 minutes at sunrise is obtained by monitoring in real time, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, the Cable Structure surface temperature measured data that each Cable Structure surface temperature measured data sequence was carved into after sunrise moment next day between 30 minutes by the same day of a Cable Structure surface point at sunrise is arranged according to time order and function order, find the maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, the same day for subtracting the temperature that minimum temperature obtains each Cable Structure surface point with the maximum temperature in each Cable Structure surface temperature measured data sequence is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, there is R Cable Structure surface point just to there is the R same day to be carved into the maximum temperature difference numerical value after sunrise moment next day between 30 minutes at sunrise, maximum therein is referred to as Cable Structure surface maximum temperature difference, it is designated as Δ T_smax；Calculated by each Cable Structure surface temperature measured data sequence by Conventional mathematical and obtain rate of change of the temperature on the time of each Cable Structure surface point, the temperature of each Cable Structure surface point on the time rate of change also with time change；Obtain the same day by monitoring in real time and be carved at sunrise after sunrise moment next day between 30 minutes, in synchronization, after HBE " temperature profile data of the Cable Structure along thickness ", calculate the difference of the maximum temperature and minimum temperature that amount at the height above sea level of each selection in BE " temperature profile data of the identical height above sea level Cable Structure along thickness ", the absolute value of this difference is referred to as " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", have chosen H different height above sea levels just has H " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", maximum in this H " Cable Structure thickness direction maximum temperature difference at identical height above sea level " is called " Cable Structure thickness direction maximum temperature difference ", it is designated as Δ T_tmax；

b3：Survey calculation obtains Cable Structure steady temperature data；First, it is determined that at the time of obtaining Cable Structure steady temperature data, the condition related at the time of Cable Structure steady temperature data to determining to obtain has six, Section 1 condition was carved into after sunrise moment next day between 30 minutes at sunset at the time of being and obtain Cable Structure steady temperature data between the same day, the sunset moment refers to that, according to the sunset moment on earth rotation and the meteorology of revolution rule determination, data can be inquired about or calculate by conventional meteorology obtaining the required sunset moment of each day；The a conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, reference plate maximum temperature difference Δ T_pmaxWith Cable Structure surface maximum temperature difference Δ T_smaxAll it is not more than 5 degrees Celsius；The b conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, the environment maximum temperature difference Δ T obtained in above survey calculation_emaxNo more than refer to temperature difference per day Δ T_r, and reference plate maximum temperature difference Δ T_pmaxSubtract and be not more than Δ T after 2 degrees Celsius_emax, and Cable Structure surface maximum temperature difference Δ T_smaxNo more than Δ T_pmax；Need to only meet in a conditions and b conditions of Section 2 one is known as meeting Section 2 condition；Section 3 condition is that the temperature of environment is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time where Cable Structure at the time of Cable Structure steady temperature data are obtained；Section 4 condition is that at the time of Cable Structure steady temperature data are obtained, the temperature of each Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time；Section 5 condition is that at the time of Cable Structure steady temperature data are obtained, the Cable Structure surface temperature measured data of each Cable Structure surface point in R Cable Structure surface point was carved into the minimum after sunrise moment next day between 30 minutes for the same day at sunrise；Section 6 condition is " Cable Structure thickness direction maximum temperature difference " Δ T at the time of Cable Structure steady temperature data are obtained_tmaxNo more than 1 degree Celsius；This method utilizes above-mentioned six conditions, any one in the following three moment is referred to as " the mathematics moment for obtaining Cable Structure steady temperature data ", the first moment is at the time of meeting Section 1 to the Section 5 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", second of moment is at the time of only meeting the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", the third moment is while at the time of meeting Section 1 to the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain "；It it is exactly the mathematics moment for obtaining Cable Structure steady temperature data at the time of acquisition Cable Structure steady temperature data when a moment in the physical record data moment during the mathematics moment for obtaining Cable Structure steady temperature data is exactly this method；If the mathematics moment for obtaining Cable Structure steady temperature data is not any one moment in the physical record data moment in this method, take this method closest to the mathematics moment for obtaining Cable Structure steady temperature data that physical record data at the time of to obtain at the time of Cable Structure steady temperature data；The amount that this method is used in measurement record at the time of obtaining Cable Structure steady temperature data carries out the monitoring analysis of Cable Structure relevant health；The Cable Structure temperature field that this method is approximately considered at the time of obtaining Cable Structure steady temperature data is in the Cable Structure temperature at stable state, i.e. this moment and not changed over time, and this moment is exactly " at the time of the obtaining Cable Structure steady temperature data " of this method；Then, according to Cable Structure heat-transfer character, utilize " the R Cable Structure surface temperature measured data " at the time of obtaining Cable Structure steady temperature data and " HBE Cable Structure along thickness temperature measured data ", utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution for obtaining Cable Structure at the time of Cable Structure steady temperature data are obtained is calculated by conventional heat transfer, now the temperature field of Cable Structure is calculated by stable state, calculating the temperature profile data of obtained Cable Structure at the time of Cable Structure steady temperature data are obtained includes the calculating temperature of R Cable Structure surface point in Cable Structure, the calculating temperature of R Cable Structure surface point is referred to as R Cable Structure steady-state surface temperature and calculates data, also include calculating temperature of the Cable Structure above selected HBE " point of the measurement Cable Structure along the temperature profile data of thickness ", the calculating temperature of HBE " points of the measurement Cable Structure along the temperature profile data of thickness " is referred to as " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data calculates data correspondent equal with R Cable Structure steady-state surface temperature, and " HBE Cable Structure along thickness temperature measured data " with " HBE Cable Structure along thickness temperature calculating data " correspondent equal when, the temperature profile data for calculating obtained Cable Structure at the time of Cable Structure steady temperature data are obtained is referred to as " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " now is referred to as " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure along thickness temperature measured data " is referred to as " HBE Cable Structure along thickness steady temperature measured data "；When " R Cable Structure surface point " is taken on the surface of Cable Structure, the quantity of " R Cable Structure surface point " must is fulfilled for three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the observed temperature linear interpolation of point adjacent with the arbitrfary point on Cable Structure surface during the temperature at any point on Cable Structure surface is by " R Cable Structure surface point " is obtained, the error of the temperature of the arbitrfary point and the actual temperature of the arbitrfary point on Cable Structure surface is not more than 5% on the Cable Structure surface that linear interpolation is obtained；Cable Structure surface includes support cable surface；Second condition is that the point being not less than in 4, and " R Cable Structure surface point " in same height above sea level in the quantity of the point of same height above sea level in " R Cable Structure surface point " is uniform along Cable Structure surface；Maximum Δ h in the absolute value of the difference of the height above sea level of " R Cable Structure surface point " along all Cable Structure surface points adjacent two-by-two of height above sea level is not more than 0.2 DEG C divided by Δ T_hObtained numerical value, Δ T is taken for convenience of narration_hUnit for DEG C/m, for convenience of narration take Δ h unit be m；When the definition of " R Cable Structure surface point " along the Cable Structure surface point adjacent two-by-two of height above sea level refers to only consider height above sea level, a Cable Structure surface point is not present in " R Cable Structure surface point ", the height above sea level numerical value of the Cable Structure surface point is between the height above sea level numerical value of adjacent Cable Structure surface point two-by-two；3rd condition is inquiry or obtains Cable Structure location and the interval sunshine rule of place height above sea level by meteorology conventionally calculation, further according to the geometric properties and bearing data of Cable Structure, found in Cable Structure it is annual by the sunshine-duration most sufficient position of those surface points, in " R Cable Structure surface point " at least one Cable Structure surface point be in Cable Structure whole year by a point in those most sufficient surface points of sunshine-duration；

C. the Cable Structure steady temperature data obtained under original state are calculated according to " temperature survey of the Cable Structure of this method calculates method " direct measurement, Cable Structure steady temperature data under original state are referred to as initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T_o”；Survey or consult reference materials and obtain the physical and mechanical properties parameter varied with temperature of various materials used in Cable Structure；Initial Cable Structure steady temperature data vector T is obtained in actual measurement_oSynchronization, direct measurement, which is calculated, obtains the Initial cable forces of all support cables, composition Initial cable force vector F_o；According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector T_oUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cable_oElement coding rule it is identical；T is obtained in actual measurement_oWhile, that is, obtaining initial Cable Structure steady temperature data vector T_oAt the time of synchronization, direct measurement, which is calculated, obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, initial Cable Structure bearing generalized coordinates data include initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, and the data of the health status that can express support cable now are referred to as support cable initial health data；The monitored amount initial value vector C of initial value composition of all monitored amounts_o, it is monitored amount initial value vector C_oCoding rule and M monitored amounts coding rules it is identical；Evaluation object initial damage vector d is set up using support cable initial health data and Cable Structure load measurement data_o, vectorial d_oRepresent with initial mechanical calculating benchmark model A_oThe initial health of the evaluation object of the Cable Structure of expression；Evaluation object initial damage vector d_oElement number be equal to N, d_oElement and evaluation object be one-to-one relationship, vectorial d_oElement coding rule it is identical with the coding rule of evaluation object；If d_oThe corresponding evaluation object of some element be a support cable in cable system, then d_oThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion；If d_oThe corresponding evaluation object of some element be to take d in some load, this method_oThe element numerical value be 0, the initial value for representing the change of this load is 0；If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial d_oIn each element numerical value related to support cable take 0；Initial Cable Structure bearing generalized coordinates data constitute initial Cable Structure bearing generalized coordinates vector U_o；

D. the physical and mechanical properties parameter varied with temperature of various materials, initial Cable Structure bearing generalized coordinates vector U according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure_o, initial Cable Structure steady temperature data vector T_oAll Cable Structure data obtained with preceding step, set up the initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data "_o, based on A_oCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot 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_oCable Structure bearing generalized coordinates data be exactly initial Cable Structure bearing generalized coordinates vector U_o；Corresponding to A_oEvaluation object health status evaluation object initial damage vector d_oRepresent；Corresponding to A_oAll monitored amounts initial value with monitored amount initial value vector C_oRepresent；U_o、T_oAnd d_oIt is A_oParameter, by A_oThe obtained initial value of all monitored amounts of Mechanics Calculation result and C_oThe initial value of all monitored amounts represented is identical, therefore alternatively C_oBy A_oMechanics Calculation result composition, A in the method_o、C_o、d_o、U_oAnd T_oIt is constant；

E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated；Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingⁱ _o, A_oAnd Aⁱ _oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Tⁱ _oRepresent, vector Tⁱ _oDefinition mode and vector T_oDefinition mode it is identical, Tⁱ _oElement and T_oElement correspond；When ith circulation starts, corresponding to Aⁱ _o" Cable Structure bearing generalized coordinates data " with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent, vectorial Uⁱ _oDefinition mode and vector U_oDefinition mode it is identical, Uⁱ _oElement and U_oElement correspond；The current initial damage vector of evaluation object that ith circulation needs when starting is designated as dⁱ _o, dⁱ _oRepresent Cable Structure A during this time circulation beginningⁱ _oEvaluation object health status, dⁱ _oDefinition mode and d_oDefinition mode it is identical, dⁱ _oElement and d_oElement correspond；When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amountⁱ _oRepresent, vectorial Cⁱ _oDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _oElement and C_oElement correspond, be monitored the current initial value vector C of amountⁱ _oRepresent to correspond to Aⁱ _oAll monitored amounts concrete numerical value；Uⁱ _o、Tⁱ _oAnd dⁱ _oIt is Aⁱ _oCharacterisitic parameter, Cⁱ _oBy Aⁱ _oMechanics Calculation result composition；When circulation starts for the first time, Aⁱ _oIt is designated as A¹ _o, set up A¹ _oMethod to make A¹ _oEqual to A_o；When circulation starts for the first time, Tⁱ _oIt is designated as T¹ _o, set up T¹ _oMethod to make T¹ _oEqual to T_o；When circulation starts for the first time, Uⁱ _oIt is designated as U¹ _o, set up U¹ _oMethod to make U¹ _oEqual to U_o；When circulation starts for the first time, dⁱ _oIt is designated as d¹ _o, set up d¹ _oMethod to make d¹ _oEqual to d_o；When circulation starts for the first time, Cⁱ _oIt is designated as C¹ _o, set up C¹ _oMethod to make C¹ _oEqual to C_o；

F. enter from here and the circulation walked to q is walked by f；During structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Tⁱ, vector TⁱDefinition mode and vector T_oDefinition mode it is identical, TⁱElement and T_oElement correspond；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, actual measurement obtains Cable Structure bearing generalized coordinates current data, and all Cable Structure bearing generalized coordinates current datas constitute current Cable Structure actual measurement bearing generalized coordinates vector Uⁱ, vectorial UⁱDefinition mode and vector U_oDefinition mode it is identical；Vector T is obtained in actual measurementⁱWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TⁱAt the time of synchronization Cable Structure in all monitored amounts currency, the monitored amount current value vector C of all these numerical value compositionⁱ, vectorial CⁱDefinition mode and vector C_oDefinition mode it is identical, CⁱElement and C_oElement correspond, represent identical monitored amount in numerical value not in the same time；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, actual measurement obtain all M in Cable Structure₁The rope force data of root support cable, all these rope force data composition current cable force vector Fⁱ, vectorial FⁱElement and vector F_oElement coding rule it is identical；Current Cable Structure steady temperature data vector T is obtained in actual measurementⁱSynchronization, Actual measurement obtains all M₁The space coordinate of two supporting end points of root support cable, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vector_oElement coding rule it is identical；

G. bearing generalized coordinates vector U is surveyed according to current Cable StructureⁱWith current Cable Structure steady temperature data vector Tⁱ, current initial mechanical calculating benchmark model A is updated according to step g1 to g3ⁱ _o, the monitored current initial value vector C of amountⁱ _o, current initial Cable Structure steady temperature data vector Tⁱ _oWith current initial Cable Structure bearing generalized coordinates vector Uⁱ _o, and the current initial damage vector d of evaluation objectⁱ _oKeep constant；

G1. it is respectively compared TⁱAnd Tⁱ _o、UⁱAnd Uⁱ _oIf, TⁱEqual to Tⁱ _oAnd UⁱEqual to Uⁱ _o, then need not be to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated, otherwise needs to follow these steps to Aⁱ _o、Uⁱ _o、Cⁱ _oAnd Tⁱ _oIt is updated；

G2. U is calculatedⁱWith U_oDifference, UⁱWith U_oDifference be exactly generalized displacement of support of the Cable Structure bearing on initial position, represent generalized displacement of support with generalized displacement of support vector V, V is equal to UⁱSubtract U_o；Calculate TⁱWith T_oDifference, TⁱWith T_oDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TⁱWith T_oDifference represented with steady temperature change vector S, S be equal to TⁱSubtract T_o, S represents the change of Cable Structure steady temperature data；

G3. first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then to A_oIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, to A_oMiddle Cable Structure bearing applies generalized displacement of support constraint and applies the current initial mechanical calculating benchmark model A updated after temperature change to Cable Structureⁱ _o, update Aⁱ _oWhile, Uⁱ _oAll elements numerical value also uses UⁱAll elements numerical value correspondence is replaced, that is, have updated Uⁱ _o, Tⁱ _oAll elements numerical value also uses TⁱAll elements numerical value correspondence replace, that is, have updated Tⁱ _o, thus obtained properly corresponding to Aⁱ _oUⁱ _oAnd Tⁱ _o, now dⁱ _oKeep constant；As renewal Aⁱ _oAfterwards, Aⁱ _oRope the health status current initial damage vector d of evaluation objectⁱ _oRepresent, Aⁱ _oCable Structure steady temperature with current Cable Structure steady temperature data vector Tⁱ _oRepresent, Aⁱ _oBearing generalized coordinates with current initial Cable Structure bearing generalized coordinates vector Uⁱ _oRepresent；Update Cⁱ _oMethod be：As renewal Aⁱ _oAfterwards, A is obtained by Mechanics Calculationⁱ _oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Cⁱ _o；

H. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingⁱWith evaluation object unit change vector Dⁱ _u；

H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4ⁱAnd Dⁱ _u；At other moment, when in step g to Aⁱ _oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CⁱAnd Dⁱ _uIf, not to A in step gⁱ _oIt is updated, then is directly transferred to step i here and carries out follow-up work；

H2. in current initial mechanical calculating benchmark model Aⁱ _oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating；According to the coding rule of evaluation object, calculated successively；Calculate each time and assume that only one of which evaluation object is further added by unit damage or load unit change on the basis of original damage or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is a load, it is assumed that the load is further added by load unit change, D is usedⁱ _ukThis increased unit damage or load unit change is recorded, wherein k represents the numbering for increasing unit damage or the evaluation object of load unit change, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uAn element, evaluation object unit change vector Dⁱ _uElement coding rule and vector d_oElement coding rule it is identical；The evaluation object that unit damage or load unit change are further added by calculating each time is different from the evaluation object that unit damage or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector；When assuming that k-th of evaluation object is further added by unit damage or load unit change, C is usedⁱ _tkRepresent corresponding " monitored amount calculation current vector "；When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed；Cⁱ _tkDefinition mode and vector C_oDefinition mode it is identical, Cⁱ _tkElement and C_oElement correspond；

H3. obtained vectorial C is calculated each timeⁱ _tkSubtract vectorial Cⁱ _oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or load unit change numerical valueⁱ _k”；There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount "；

H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "ⁱ”；Unit damage monitored numerical quantity transformation matrices Δ CⁱEach row correspond to a monitored amount unit change vector；Unit damage monitored numerical quantity transformation matrices Δ CⁱEvery a line correspond to same monitored amount different evaluation objects increase unit damage or load unit change when different unit change amplitudes；Unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and vector d_oElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CⁱRow coding rule and M monitored amounts coding rules it is identical；

I. current nominal fatigue vector d is definedⁱ _cWith currently practical injury vector dⁱ, dⁱ _cAnd dⁱElement number be equal to evaluation object quantity, dⁱ _cAnd dⁱElement and evaluation object between be one-to-one relationship, dⁱ _cElement numerical value represent correspondence evaluation object nominal fatigue degree or nominal load variable quantity, dⁱ _cAnd dⁱWith evaluation object initial damage vector d_oElement number rule it is identical, dⁱ _cElement, dⁱElement and d_oElement be one-to-one relationship；

J. according to monitored amount current value vector CⁱWith " the monitored current initial value vector C of amountⁱ _o", " unit damage monitored numerical quantity transformation matrices Δ Cⁱ" and " current nominal fatigue vector dⁱ _c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1ⁱ _cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector dⁱ _c；

Formula 1

K. the currently practical injury vector d expressed using formula 2ⁱK-th of element dⁱ _kWith the current initial damage vector d of evaluation objectⁱ _oK-th of element dⁱ _okWith current nominal fatigue vector dⁱ _cK-th of element dⁱ _ckBetween relation, calculating obtain currently practical injury vector dⁱAll elements；

Formula 2

K=1,2,3 in formula 2 ..., N；Vectorial dⁱElement coding rule and formula (1) in vector d_oElement coding rule it is identical；dⁱ _kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then dⁱ _kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, dⁱ _kThe numerical response degree of relaxation or the damage of the support cable；By the currently practical injury vector d of evaluation objectⁱIn with M₁The related M of root support cable₁Individual element takes out, the currently practical injury vector d of composition support cable^ci, the currently practical injury vector d of support cable^ciElement coding rule and Initial cable force vector F_oElement coding rule it is identical；The currently practical injury vector d of support cable^ciH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., M₁；The currently practical injury vector d of support cable^ciMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line；Support cable currently practical injury vector d corresponding with damaged cable^ciIn element numerical expression be the damaged cable currently practical damage, element numerical value be 100% when represent that the support cable thoroughly loses bearing capacity, represented when between 0 and 100% lose corresponding proportion bearing capacity；Using in current Cable Structure steady temperature data vector TⁱUnder the conditions of, in l walk the slack line that identifies and with the currently practical injury vector d of support cable^ciExpression these slack lines, the degree that relaxed with it mechanic equivalent currently practical equivalent damage degree, using f walk acquisition in current Cable Structure steady temperature data vector TⁱUnder the conditions of current cable force vector FⁱWith current support cable two supporting end points horizontal range vector, using c walk obtain in initial Cable Structure steady temperature data vector T_oUnder the conditions of support cable initial drift is vectorial, the weight of the initial free unit length of initial free cross-sectional area vector sum is vectorial, Initial cable force vector F_o, utilize current Cable Structure steady temperature data vector TⁱThe support cable current steady state temperature data of expression, using c walk obtain in initial Cable Structure steady temperature data vector T_oThe support cable initial steady state temperature data of expression, utilize the physical and mechanical properties parameter varied with temperature that various materials used in the Cable Structure of acquisition are walked in c, it is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is：First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical；2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope；When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa；Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined；So it is achieved that relaxation identification and the non-destructive tests of support cable；Institute's demand power is by current cable force vector F during calculatingⁱCorresponding element is provided；Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far this method realizes the problem of rejecting generalized displacement of support, load change and the influence of structure temperature change, Cable Structure rope and recognized, generalized displacement of support, structure temperature change and the change influence of support cable health status, load change amount identification are rejected while realizing；

L. current nominal fatigue vector d is tried to achieveⁱ _cAfterwards, mark vector B is set up according to formula 3ⁱ, formula 4 gives mark vector BⁱK-th of element definition；

Formula 3

Formula 4

Element B in formula 4ⁱ _kIt is mark vector BⁱK-th of element, Dⁱ _ukIt is evaluation object unit change vector Dⁱ _uK-th of element, dⁱ _ckIt is the current nominal fatigue vector d of evaluation objectⁱ _cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4；

If m. mark vector BⁱElement be all 0, then return to step f continue this circulation；If mark vector BⁱElement be not all 0, then enter next step, i.e. step n；

N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5ⁱ⁺¹ _oEach element；

Formula 5

D in formula 5ⁱ⁺¹ _okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationⁱ⁺¹ _oK-th of element, dⁱ _okThis, i.e. the current initial damage vector d of evaluation object of ith circulationⁱ _oK-th of element, Dⁱ _ukIt is the evaluation object unit change vector D of ith circulationⁱ _uK-th of element, Bⁱ _kIt is the mark vector B of ith circulationⁱK-th of element, k=1,2,3 ... ..., N in formula 5；

O. in initial mechanical calculating benchmark model A_oOn the basis of, first to A_oIn Cable Structure bearing apply generalized displacement of support constraint, the numerical value of generalized displacement of support constraint is just derived from the numerical value of corresponding element in generalized displacement of support vector V, then to A_oIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, then it is d to make the health status of ropeⁱ⁺¹ _oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationⁱ⁺¹；Obtain Aⁱ⁺¹Afterwards, A is obtained by Mechanics Calculationⁱ⁺¹In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationⁱ⁺¹ _o；

P. remove once, i.e. i+1 time circulates required current initial Cable Structure steady temperature data vector Tⁱ⁺¹ _oThe current initial Cable Structure steady temperature data vector T circulated equal to ithⁱ _o；Next time, the current initial Cable Structure bearing generalized coordinates vector U i.e. needed for i+1 time circulationⁱ⁺¹ _oThe current initial Cable Structure bearing generalized coordinates vector U circulated equal to ithⁱ _o；

Q. step f is returned to, starts to circulate next time.