CN105067354A - A progressive identification method for cable loads in simplified generalized displacement mixed monitoring problems - Google Patents

Abstract

Simplified generalized displacement mixed monitoring problem The progressive identification method of cable load is based on mixed monitoring. By monitoring the generalized displacement of the support, the temperature of the cable structure and the ambient temperature, it is determined whether it is necessary to update the mechanical calculation benchmark model of the cable structure, and the generalized displacement of the support is included. , cable structure temperature and ambient temperature mechanical calculation benchmark model of the cable structure, on the basis of this model, the numerical change matrix of the unit damage monitored quantity is calculated and obtained. According to the approximate linear relationship between the current numerical vector of the monitored quantity and the current initial numerical vector of the monitored quantity, the numerical change matrix of the monitored quantity per unit damage, and the current nominal damage vector of the assessed object to be calculated, the value of the current nominal damage vector of the assessed object is calculated. Non-inferior solutions, based on which the health status of the core evaluated objects can be identified.

Description

A progressive identification method for cable loads in simplified generalized displacement mixed monitoring problems

technical field

Cable-stayed bridges, suspension bridges, truss structures and other structures have one thing in common, that is, they have many parts that bear tensile loads, such as cable-stayed cables, main cables, slings, tie rods, etc. , cables, or rods that only bear tensile loads are supporting components. For the sake of convenience, this method expresses this type of structure as For the sake of convenience, the rods with tensile or axial compression loads (also known as two-force rods) are collectively referred to as "cable systems". In this method, the term "supporting cables" is used to refer to load-bearing cables, A member that bears axial tension or axial compression load is sometimes referred to as a "cable", so when the word "cable" is used later, it actually refers to a two-force member for a truss structure. The damage and slack of supporting cables is a major threat to the safety of cable structures. In this method, damaged cables and slack cables are collectively referred to as supporting cables with health problems, referred to as problem cables for short. During the service process of the structure, the correct identification of the health state of the supporting cables or the cable system is related to the safety of the entire cable structure. When the ambient temperature changes, the temperature of the cable structure will generally change accordingly. When the temperature of the cable structure changes, the support of the cable structure may undergo a generalized displacement, and the load on the cable structure may also change. The health status may also be changing. Under such complex conditions, this method is based on hybrid monitoring (this method judges the health status of the cable structure through the hybrid monitoring of changes in the measurable parameters of different types of cable structures mentioned in this section. In the method, all the monitored characteristic parameters of the cable structure are collectively referred to as "monitored quantity". Since the monitored quantity is composed of different types of measurable parameters of the cable structure, this method is called mixed monitoring) to identify problem cables. , which belongs to the field of engineering structure health monitoring.

Background technique

Eliminating the influence of load changes, generalized displacements of cable structure supports, and structural temperature changes on the identification results of the cable structure's health status, thereby accurately identifying changes in the health status of the structure is an urgent problem that needs to be solved. This method discloses the solution to this problem an effective and inexpensive method.

Contents of the invention

Technical problem: This method discloses a method, under the condition of lower construction cost, when the support has a generalized displacement, when the load on the structure and the structure (environment) temperature change, the generalized displacement and load change of the support can be eliminated And the impact of structural temperature changes on the identification results of the cable structure health status, so as to accurately identify the health status of the supporting cables.

Technical solution: In this method, "support space coordinates" refers to the coordinates of the support on the X, Y, and Z axes of the Cartesian coordinate system, and can also be said to be the space of the support on the X, Y, and Z axes. Coordinates, the specific value of the spatial coordinate of the support about a certain axis is called the spatial coordinate component of the support about the axis. In this method, a spatial coordinate component of the support is also used to express the specific value of the spatial coordinate of the support about a certain axis. Numerical value; use "support angular coordinates" to refer to the angular coordinates of the support with respect to the X, Y, and Z axes, and the specific value of the angular coordinates of the support with respect to a certain axis is called the angular coordinate component of the support with respect to the axis. In this method An angular coordinate component of the support is also used to express the specific value of the angular coordinate of the support with respect to a certain axis; the "generalized coordinate of the support" is used to refer to the angular coordinates of the support and the space coordinates of the support as a whole. In this method, the support's A generalized coordinate component expresses the specific value of the space coordinate or angular coordinate of the support about an axis; the change of the coordinate of the support about the X, Y, and Z axes is called the line displacement of the support, and it can also be said that the change of the space coordinate of the support is called is the linear displacement of the support. In this method, a linear displacement component of the support is also used to express the specific value of the linear displacement of the support about a certain axis; the change of the angular coordinates of the support about the X, Y, and Z axes is called the support Angular displacement, in this method, an angular displacement component of the support is also used to express the specific value of the angular displacement of the support with respect to a certain axis; the generalized displacement of the support refers to the overall displacement of the support line and the angular displacement of the support, and this method also uses A generalized displacement component of the support expresses the specific value of the linear displacement or angular displacement of the support about a certain axis; the linear displacement of the support can also be called translational displacement, and the settlement of the support is the linear displacement or translational displacement of the support in the direction of gravity portion.

The external force borne by objects and structures can be called load, and the load includes surface load and body load. Surface load, also known as surface load, is a load acting on the surface of an object, including concentrated load and distributed load. Body load is the load continuously distributed at various points inside the object, such as the object's own weight and inertial force.

Concentrated load is divided into concentrated force and concentrated force couple. In a coordinate system, such as a Cartesian rectangular coordinate system, a concentrated force can be decomposed into three components. Similarly, a concentrated force couple can also be decomposed into three components. If the load is actually a concentrated load, a concentrated force component or a concentrated force couple component is called a load in this method, and the change of the load at this time is embodied as a change of a concentrated force component or a concentrated force couple component.

Distributed load is divided into line distributed load and surface distributed load. The description of distributed load includes at least the action area of distributed load and the size of distributed load. , sine function and other distribution characteristics) and amplitude to express (for example, two distributed loads are uniformly distributed, but their amplitudes are different, uniform pressure can be used as an example to illustrate the concept of amplitude: the same structure bears two different loads Uniformly distributed pressure, both distributed loads are uniformly distributed loads, but the amplitude of one distributed load is 10MPa, and the amplitude of the other distributed load is 50MPa). If the load is actually a distributed load, when this method talks about the change of the load, it actually refers to the change of the amplitude of the distribution concentration of the distributed load, while the action area of the distributed load and the distribution characteristics of the distribution concentration remain unchanged. In the coordinate system, a distributed load can be decomposed into several components. If the amplitudes of the respective distribution concentration of several components of the distributed load change, and the ratios of the changes are not all the same, then in this method, the Several distributed load components are regarded as the same number of independent distributed loads. At this time, one load represents a distributed load component, and the components with the same amplitude change ratio of the distribution concentration can also be combined into a distributed load or called for a load.

Body load is the load continuously distributed at various points inside the object, such as the self-weight and inertial force of the object. The description of the body load includes at least the action area of the body load and the size of the body load. The size of the body load is expressed by the degree of distribution, and the distribution Concentration is expressed by distribution characteristics (such as uniform distribution, linear function and other distribution characteristics) and amplitude (for example, two body loads are uniformly distributed, but their amplitudes are different. The concept of amplitude can be illustrated by self-weight as an example: the same The two parts of a structure are made of different materials and therefore have different densities, so although the body loads on both parts are uniformly distributed, the magnitude of the body load on one part may be 10kN/m ³ and the other may be A part is subjected to a body load of magnitude 50kN/m ³ ). If the load is actually a body load, what is actually dealt with in this method is the change in the magnitude of the distribution concentration of the body load, while the action area of the body load and the distribution characteristics of the distribution concentration remain unchanged. At this time, in this method When referring to the change of the load in , it actually refers to the change of the magnitude of the distribution concentration of the body load. At this time, the changed load refers to the body load whose distribution concentration changes. In the coordinate system, a body load can be decomposed into several components (for example, in the Cartesian coordinate system, the body load can be decomposed into components about the three axes of the coordinate system, that is, in the Cartesian coordinate system body load can be decomposed into three components), if the magnitude of the respective distribution concentration of several components of the body load changes, and the ratios of changes are not all the same, then in this method, the The components are regarded as the same number of independent loads, and the body load components with the same amplitude change ratio of the distribution concentration can also be synthesized into a body load or called a load.

When the load is embodied as concentrated load, in this method, "load unit change" actually refers to "unit change of concentrated load", similarly, "load change" specifically refers to "change in magnitude of concentrated load", " "Load change" specifically refers to "the change in the size of the concentrated load", "the degree of load change" specifically refers to the "change degree of the concentrated load", and "the actual change in the load" refers to the "actual change in the size of the concentrated load". Amount", "changed load" refers to "concentrated load whose size changes". Simply put, at this time, "a certain change in a certain load" refers to "a certain change in the size of a certain concentrated load".

When the load is embodied as a distributed load, in this method, "the unit change of the load" actually refers to "the unit change of the amplitude of the distribution concentration of the distributed load", and the distribution characteristics of the distributed load are unchanged, similar to The "load change" specifically refers to "the change in the magnitude of the distribution concentration of the distributed load", while the distribution characteristics of the distributed load are unchanged, and the "load change" specifically refers to the "magnitude of the distribution concentration of the distributed load". The change amount of the load", "the degree of load change" specifically refers to the "change degree of the amplitude of the distribution concentration of the distributed load", and the "actual change of the load" specifically refers to the "actual change of the amplitude of the distribution concentration of the distributed load ", "changed load" refers to "distributed load whose amplitude of distribution intensity changes". A certain change in the amplitude of the load", while the distribution characteristics of the action area and distribution concentration of all distributed loads remain unchanged.

When the load is embodied as body load, in this method, "load unit change" actually refers to "unit change in the amplitude of distribution concentration of body load", similarly, "load change" means "body load The change of the magnitude of the distribution concentration of the body load", the "load change" refers to the "change of the magnitude of the distribution concentration of the body load", and the "load change degree" refers to the "magnitude of the distribution concentration of the body load The degree of change", "the actual change of the load" refers to the "actual change of the magnitude of the distribution concentration of the body load", and "the changed load" refers to the "body load whose distribution concentration changes. ", simply put, "a certain change in a certain load" refers to "a certain change in the magnitude of the distribution concentration of a certain body load", and the distribution characteristics of the action area and distribution concentration of all body loads are Changeless.

This method specifically includes:

a. When the load borne by the cable structure changes, but the load the cable structure is bearing does not exceed the initial allowable load of the cable structure, this method is applicable; the initial allowable load of the cable structure refers to the allowable load of the cable structure at the time of completion, It can be obtained through conventional mechanical calculations; this method collectively refers to the evaluated supporting cables and loads as "assessed objects", and the sum of the number of evaluated supporting cables and the number of loads is N, which is the number of "assessed objects" is N; this method uses the name "core assessed object" to specifically refer to the assessed support cables in the "assessed object", and this method uses the name "secondary assessed object" to specifically refer to the assessed support cable in the "assessed object". load; determine the numbering rule of the evaluated object, according to this rule, number all the evaluated objects in the cable structure, and the number will be used to generate vectors and matrices in subsequent steps; this method uses the variable k to represent this number, k = 1, 2, 3, ..., N; determine the specified supporting cables to be monitored in the mixed monitoring, assuming that there are Q supporting cables in the cable system, obviously the number of core objects to be evaluated is Q; The monitored cable force data is described by M ₁ cable force data of M ₁ designated supporting cables on the cable structure, and the change of the cable structure cable force is the change of the cable force of all designated supporting cables; each time there are M ₁ cable forces The measured value or calculated value is used to represent the cable force information of the cable structure; M ₁ is an integer not less than 0 and not greater than Q; the measured point to be monitored for the specified strain in the mixed monitoring is determined, and the monitored strain data of the cable structure It can be described by K ₂ designated points on the cable structure and the strains in L ₂ designated directions of each designated point. The change of the strain data of the cable structure is the change of all the measured strains of the K ₂ designated points; M ₂ strain measurement values or calculated values to characterize the strain of the cable structure, M ₂ is the product of K ₂ and L ₂ ; M ₂ is an integer not less than 0; the measured point to be monitored at the angle specified during the mixed monitoring is determined, The monitored angle data of the cable structure is described by K ₃ designated points on the cable structure, L ₃ designated straight lines passing through each designated point, and H ₃ angular coordinate components of each designated straight line. The change is the change of all specified points, all specified straight lines, and all specified angular coordinate components; each time there are M ₃ measured or calculated values of angular coordinate components to represent the angle information of the cable structure, M ₃ is K ₃ , L ₃ and H ₃ product; M ₃ is an integer not less than 0; the shape data to be monitored specified when determining the mixed monitoring, the monitored shape data of the cable structure is composed of K ₄ specified points on the cable structure, and The spatial coordinates of L ₄ designated directions of each designated point are used to describe the change of the cable structure shape data is the change of all coordinate components of the K ₄ designated points; each time there are M ₄ coordinate measurement values or calculated values to represent the cable structure. Structure shape, M ₄ is the product of K ₄ and L ₄ ; M ₄ is an integer not less than 0; based on the monitored quantities of the above mixed monitoring, the entire cable structure has a total of M monitored quantities, and M is M ₁ and M ₂ , the sum of M ₃ and M ₄ , Define the parameter K, K is the sum of M ₁ , K ₂ , K ₃ and K ₄ , M must be greater than the number of core evaluated objects, and M is smaller than the number of evaluated objects; for convenience, in this method, the The M monitored quantities listed are referred to as "monitored quantities" for short; the time interval between any two measurements of the same quantity monitored in real time in this method shall not be greater than 30 minutes, and the moment when the recorded data is measured is called the actual recorded data Time; the external force borne by objects and structures can be called load, which includes surface load and body load; surface load is also called surface load, which is a load acting on the surface of an object, including concentrated load and distributed load; body load is a continuous distribution The load on each point inside the object, including the self-weight and inertial force of the object; the concentrated load is divided into two types: concentrated force and concentrated force couple. In the coordinate system including the Cartesian coordinate system, a concentrated force can be decomposed into Three components. Similarly, a concentrated force couple can also be decomposed into three components. If the load is actually a concentrated load, in this method, a concentrated force component or a concentrated force couple component is counted or counted as a load. At this time The change of load is embodied as the change of a concentrated force component or a concentrated force couple component; distributed load is divided into line distributed load and surface distributed load, and the description of distributed load includes at least the action area of distributed load and the size of distributed load. The size is expressed by the distribution concentration, and the distribution concentration is expressed by the distribution characteristics and amplitude; if the load is actually a distributed load, when this method talks about the change of the load, it actually refers to the change of the amplitude of the distribution concentration of the distributed load , while the distribution characteristics of the action area and distribution concentration of all distributed loads remain unchanged; in a coordinate system including the Cartesian coordinate system, a distributed load can be decomposed into three components, if the three components of the distributed load The amplitudes of the respective distribution concentration of each component change, and the rate of change is not all the same, then in this method, the three components of the distributed load are counted or counted as three distributed loads, and one load is Represents a component of the distributed load; the body load is the load continuously distributed at each point inside the object, the description of the body load includes at least the action area of the body load and the size of the body load, the size of the body load is expressed by the degree of distribution, and the distribution set The degree is expressed by the distribution characteristics and amplitude; if the load is actually a body load, what is actually dealt with in this method is the change of the magnitude of the body load distribution concentration, while the action area of all body loads and the distribution of the distribution concentration The characteristics are unchanged. At this time, when the change of the load is mentioned in this method, it actually refers to the change of the magnitude of the distribution concentration of the body load. At this time, the changed load refers to the magnitude of the distribution concentration Body loads that change; in a coordinate system including Cartesian coordinates, a body load can be decomposed into three components, if the magnitude of the distribution of the three components of the body load changes, and the change ratios are not all the same, then in this method, the three components of the body load are counted or counted as three distributed loads;

b. This method defines "the temperature measurement and calculation method of the cable structure of this method" according to steps b1 to b3;

b1: Query or measure the temperature-dependent heat transfer parameters of the composition materials of the cable structure and the environment where the cable structure is located, use the design drawings, as-built drawings of the cable structure, and the geometrically measured data of the cable structure, and use these data and parameters to establish a cable structure. The heat transfer calculation model of the structure; query the meteorological data of not less than 2 years in recent years where the cable structure is located, and count the number of cloudy days during this period as T cloudy days. In this method, the number of cloudy days that cannot be seen during the day The whole day of the sun is called a cloudy day, and the highest and lowest temperature between 0:00 and 30 minutes after the sunrise time of the next day on each of the T cloudy days is counted. The meteorological sunrise time determined by the sun and the revolution law does not mean that the sun must be visible on that day, and the required sunrise time of each day can be obtained by querying the data or through conventional meteorological calculations. Every cloudy day from 0 o'clock to The maximum temperature difference between the highest temperature and the lowest temperature within 30 minutes after sunrise of the next day is called the maximum temperature difference of the daily temperature of the cloudy day. If there are T cloudy days, there will be the maximum temperature difference of the daily temperature of T cloudy days. Take T The maximum value of the maximum temperature difference of daily air temperature in a cloudy day is the reference daily temperature difference, which is recorded as ΔT _r ; query the location of the cable structure and the altitude interval of not less than 2 years of meteorological data in recent years or obtain the cable structure from actual measurements The change data and change law of the temperature of the environment with time and altitude, and the maximum change rate ΔT of the temperature of the environment where the cable structure is located with respect to the altitude in recent years for the location of the cable structure and the altitude interval of not less than 2 years in recent years _h , for the convenience of description, the unit of ΔT _h is ℃/m; take "R cable structure surface points" on the surface of the cable structure, and the specific principle of taking "R cable structure surface points" is described in step b3, and later The temperature of the surface points of the R cable structures will be obtained through actual measurement, and the temperature data obtained from the actual measurement is called "the measured data of the surface temperature of the R cable structures". The temperature of the surface points of the R cable structures is called the calculated temperature data as "calculation data of the surface temperature of the R cable structures"; For less than three different altitudes, at each selected altitude, at least two points are selected at the intersection of the horizontal plane and the surface of the cable structure, and the outer normal of the surface of the structure is indexed from the selected point, all selected The outer normal direction is called "the direction of measuring the temperature distribution of the cable structure along the wall thickness". The direction of temperature distribution along the wall thickness must include the outer normal direction of the sunny surface of the cable structure and the outer normal direction of the shaded surface of the cable structure, and the temperature distribution along the wall thickness of each measurement cable structure is uniform in the cable structure. Select no less than three points, for the support cable along the direction of each measuring cable structure temperature distribution along the wall thickness, only take one point, only measure the temperature of the surface point of the support cable, measure the temperature of all selected points, measure The obtained temperature is called "the temperature distribution data of the cable structure along the thickness", Among them, the "temperature distribution data of the cable structure along the thickness" measured along the same "intersection line between the horizontal plane and the surface of the cable structure" and "the direction of measuring the temperature distribution of the cable structure along the wall thickness" is called in this method "The temperature distribution data of the cable structure along the thickness at the same altitude", assume that H different altitudes are selected, and at each altitude, B are selected to measure the temperature distribution direction of the cable structure along the wall thickness, along each To measure the temperature distribution direction of the cable structure along the wall thickness, E points are selected in the cable structure, where H and E are not less than 3, and B is not less than 2. For the supporting cable, E is equal to 1. On the cable structure, "measurement cable structure The total number of "points of temperature distribution data along the thickness" is HBE, and the temperature of these HBE "points of temperature distribution data along the thickness of the cable structure" will be obtained through actual measurement later, and the temperature data obtained by actual measurement are called "HBE cable The measured data of the temperature along the thickness of the structure”, if the heat transfer calculation model of the cable structure is used to calculate the temperature of the HBE points that measure the temperature distribution data of the cable structure along the thickness, the calculated temperature data is called "HBE cable structure temperature calculation data along the thickness"; select a location at the location of the cable structure according to the meteorological temperature measurement requirements, and the temperature of the environment where the cable structure is located that meets the meteorological temperature measurement requirements will be measured at this location; at the location of the cable structure Choose a location in an open and unsheltered place, which should be able to get the fullest sunshine of the day every day of the year, and place a carbon steel plate at this location, called the reference The reference plate should not be in contact with the ground. The distance between the reference plate and the ground should not be less than 1.5 meters. One side of the reference plate faces the sun, called the sunny side. The sunny side of the reference plate is rough and dark. On every day of the year, the most sufficient sunshine of the day that a flat panel can get in this place can be obtained, and the non-sun-facing side of the reference flat-panel is covered with thermal insulation material, and the temperature of the sunny side of the reference flat-panel will be obtained through real-time monitoring;

b2: Obtain the measured surface temperature data of the R cable structures at the surface points of the above R cable structures through real-time monitoring, and at the same time obtain the temperature distribution data of the previously defined cable structures along the thickness through real-time monitoring. The temperature data of the environment where the structure is located; through real-time monitoring, the temperature measured data series of the environment where the cable structure is located is obtained from the sunrise time of the day to 30 minutes after the sunrise time of the next day. The measured temperature data of the environment where the cable structure is located between the time and 30 minutes after sunrise of the next day are arranged in chronological order, and the highest temperature and the lowest temperature in the air temperature measurement data sequence of the environment where the cable structure is located are found. The maximum temperature difference in the air temperature measured data sequence minus the minimum temperature is the maximum temperature difference between the sunrise time of the day and 30 minutes after the sunrise time of the next day in the environment where the cable structure is located, which is called the maximum temperature difference of the environment, and is recorded as ΔT _emax ; The temperature measured data sequence of the environment where the cable structure is located is obtained through conventional mathematical calculations. The measured data series of the temperature of the sunny side of the reference plate between 30 minutes, the measured data series of the temperature of the sunny side of the reference plate from the sunrise time of the current day to the sunny side of the reference plate 30 minutes after the sunrise time of the next day The measured data of the temperature of the reference plate are arranged in chronological order, find the highest temperature and the lowest temperature in the measured data sequence of the temperature of the sunny side of the reference plate, and subtract the lowest temperature from the highest temperature in the measured data sequence of the temperature of the sunny side of the reference plate Temperature obtains the maximum temperature difference between the sunrise time of the day and 30 minutes after the sunrise time of the next day from the temperature of the sunny side of the reference plate, which is called the maximum temperature difference of the reference plate, and is recorded as ΔT _pmax ; The measured data series of cable structure surface temperature of all R cable structure surface points between 30 minutes after the sunrise time of the next day, there are R cable structure surface temperature measured data sequences of R cable structure surface points, each cable structure The measured data sequence of the surface temperature is arranged in chronological order by the measured data of the surface temperature of the cable structure between the sunrise time of the day and 30 minutes after the sunrise time of the next day at a surface point of the cable structure, and the measured data sequence of the surface temperature of each cable structure is found From the highest temperature and the lowest temperature in the data series, subtract the lowest temperature from the highest temperature in the measured data series of the surface temperature of each cable structure to obtain the temperature of each cable structure surface point from the sunrise time of the day to 30 minutes after the sunrise time of the next day If there are R surface points of the cable structure, there are R maximum temperature difference values between the sunrise time of the day and 30 minutes after the sunrise time of the next day, and the maximum value is called the maximum temperature difference of the cable structure surface, which is denoted as ΔT _smax ; From the measured data sequence of the surface temperature of each cable structure, the rate of change of the temperature of each cable structure surface point with respect to time is obtained through conventional mathematical calculations, and the temperature of each cable structure surface point with respect to time The rate of change also changes with time; after obtaining HBE "temperature distribution data along the thickness of the cable structure" at the same time between the sunrise time of the day and 30 minutes after the sunrise time of the next day through real-time monitoring, calculate the The difference between the maximum temperature and the minimum temperature in a total of BE "temperature distribution data of the cable structure along the thickness at the same altitude" at a selected altitude, the absolute value of this difference is called "thickness direction of the cable structure at the same altitude". If H different altitudes are selected, there will be H "maximum temperature differences in the thickness direction of the cable structure at the same altitude", and the maximum value of the H "maximum temperature differences in the thickness direction of the cable structure at the same altitude" is "The maximum temperature difference in the thickness direction of the cable structure", denoted as ΔT _tmax ;

b3: Measurement and calculation to obtain the steady-state temperature data of the cable structure; firstly, determine the moment to obtain the steady-state temperature data of the cable The time of the structural steady-state temperature data is between the sunset time of the current day and 30 minutes after the sunrise time of the next day. The sunset time refers to the meteorological sunset time determined according to the laws of the earth's rotation and revolution. The data can be queried or through conventional weather The required sunset time of each day can be obtained through scientific calculation; the second condition a condition is that during the period between the sunrise time of the current day and 30 minutes after the sunrise time of the next day, the maximum temperature difference ΔT _pmax and the maximum temperature difference of the reference plate The maximum temperature difference ΔT _smax on the surface of the cable structure is not greater than 5 degrees Celsius; the b condition of the second condition is that during the period between the sunrise of the day and 30 minutes after the sunrise of the next day, in the environment obtained by the previous measurement and calculation The maximum temperature difference ΔT _emax is not greater than the reference daily temperature difference ΔT _r , and the maximum temperature difference ΔT _pmax of the reference plate minus 2 degrees Celsius is not greater than ΔT _emax , and the maximum temperature difference ΔT _smax on the surface of the cable structure is not greater than ΔT _pmax ; only the second item a One of the condition and b condition is said to meet the second condition; the third condition is that at the moment when the steady-state temperature data of the cable structure is obtained, the absolute value of the temperature change rate of the environment where the cable structure is located with respect to time is not greater than each 0.1 degrees Celsius per hour; the fourth condition is that at the moment when the steady-state temperature data of the cable structure is obtained, the absolute value of the temperature change rate of each cable structure surface point in the R cable structure surface points with respect to time is not greater than 0.1 degrees Celsius per hour The fifth condition is that at the moment when the steady-state temperature data of the cable structure is obtained, the measured data of the cable structure surface temperature of each cable structure surface point in the R cable structure surface points is from the sunrise time of the day to the day after the sunrise time of the next day The minimum value between 30 minutes; the sixth condition is that at the moment when the steady-state temperature data of the cable structure is obtained, the "maximum temperature difference in the thickness direction of the cable structure" ΔT _tmax is not greater than 1 degree Celsius; this method utilizes the above six conditions, and the following Any one of the three kinds of time is called "mathematical time for obtaining the steady-state temperature data of the cable structure". Item 1 to item 5 conditions, the second type of time is the time that only meets the sixth condition in the above "conditions related to the time to determine the time to obtain the steady-state temperature data of the cable structure", and the third type of time is when the above-mentioned conditions are met at the same time The moment of the first item to the sixth condition in "Conditions related to the moment of determining the time to obtain the steady-state temperature data of the cable structure"; At one moment, the moment of obtaining the steady-state temperature data of the cable structure is the mathematical moment of obtaining the steady-state temperature data of the cable structure; if the mathematical moment of obtaining the steady-state temperature data of the cable structure is not any moment in the actual recording data moment in this method, Then take this method closest to the mathematical time of obtaining the steady-state temperature data of the cable structure The actual recorded data moment is the moment when the steady-state temperature data of the cable structure is obtained; this method will use the measured and recorded quantity at the moment when the steady-state temperature data of the cable structure is obtained to perform cable structure-related health monitoring analysis; this method approximately considers that the obtained The temperature field of the cable structure at the moment of the steady-state temperature data of the cable structure is in a steady state, that is, the temperature of the cable structure at this moment does not change with time, and this moment is the "moment of obtaining the steady-state temperature data of the cable structure" of this method; Heat transfer characteristics of the structure, using the “measured surface temperature data of R cable structures” and “measured temperature data of HBE cable structures along the thickness” at the moment when the steady-state temperature data of the cable structure is obtained, and using the heat transfer calculation model of the cable structure, through The temperature distribution of the cable structure at the moment of obtaining the steady-state temperature data of the cable structure is obtained by conventional heat transfer calculation. At this time, the temperature field of the cable structure is calculated according to the steady state. The temperature distribution data of the structure includes the calculated temperature of the R cable structure surface points on the cable structure, and the calculated temperature of the R cable structure surface points is called the R cable structure steady-state surface temperature calculation data, and also includes the previously selected cable structure The calculated temperature of the HBE "points for measuring the temperature distribution data of the cable structure along the thickness", and the calculated temperature of the HBE "points for measuring the temperature distribution data of the cable structure along the thickness" are called "calculated data of the temperature of the HBE cable structure along the thickness" , when the measured surface temperature data of the R cable structures is equal to the calculated data of the steady-state surface temperature of the R cable structures, and the "measured data of the temperature along the thickness of the HBE cable structures" corresponds to the "calculated data of the temperature along the thickness of the HBE cable structures". When they are equal, the calculated temperature distribution data of the cable structure at the moment when the cable structure steady-state temperature data is obtained is called "cable structure steady-state temperature data" in this method, and the "R cable structure surface temperature measured data " is called "the measured data of the steady-state surface temperature of R cable structures", and "the measured data of the temperature of the HBE cable structures along the thickness" is called "the measured data of the steady-state temperature of the HBE cable structures along the thickness"; When "R cable structure surface points", the number and distribution of "R cable structure surface points" must meet three conditions. The first condition is that when the cable structure temperature field is in a steady state, when any point on the cable structure surface When the temperature of is obtained by the linear interpolation of the measured temperature of the point adjacent to the arbitrary point on the surface of the cable structure among the "R cable structure surface points", the temperature of the arbitrary point on the surface of the cable structure obtained by linear interpolation is the same as that of the surface of the cable structure The error of the actual temperature at this arbitrary point is not greater than 5%; the cable structure surface includes the supporting cable surface; the second condition is that the number of points at the same altitude in the "R cable structure surface points" is not less than 4, and " The points at the same altitude in the "R cable structure surface points" are uniformly distributed along the cable structure surface; The maximum value Δh among the values is not greater than the value obtained by dividing ΔT _h by 0.2°C. For the convenience of description, the unit of ΔT _h is ℃/m, and for the convenience of description, the unit of Δh is m;" The definition of two adjacent cable structure surface points along the altitude of "R cable structure surface points" means that when only the altitude is considered, there is no cable structure surface point in the "R cable structure surface points". The altitude value of the surface point is between the altitude values of two adjacent cable structure surface points; the third condition is to query or calculate according to meteorological routines to obtain the sunshine law of the cable structure location and the altitude interval, and then according to the cable structure According to the geometric characteristics and orientation data of the structure, find the positions of those surface points on the cable structure that receive the most sunshine time throughout the year, and at least one cable structure surface point in the "R cable structure surface points" is on the cable structure that is exposed to sunlight throughout the year. one of those surface points at which the hours of sunshine are greatest;

c. According to the "temperature measurement and calculation method of the cable structure of this method", the steady-state temperature data of the cable structure in the initial state is directly measured and calculated, and the steady-state temperature data of the cable structure in the initial state is called the initial steady-state temperature data of the cable structure. It is recorded as "initial cable structure steady-state temperature data vector T _o "; the physical and mechanical performance parameters of various materials used in the cable structure that vary with temperature are obtained from actual measurement or information; the initial cable structure steady-state temperature data vector is obtained from the actual measurement At the same moment T _o , the initial cable forces of all supporting cables are directly measured and calculated to form the initial cable force vector F _o ; according to the data including cable structure design data and completed data, all supporting cables are in the free state, that is, the cable force is The length at 0, the cross-sectional area in the free state and the weight per unit length in the free state, as well as the temperature of all supporting cables when these three data are obtained, on this basis, the temperature-dependent Physical performance parameters and mechanical performance parameters, according to conventional physical calculations, the lengths of all supporting cables when the cable force is 0 and the lengths of all supporting cables when the cable force is 0 under the condition of the initial cable structure steady-state temperature data vector T _o The cross-sectional area and the weight per unit length of all supporting cables when the cable force is 0 form the initial free length vector of the supporting cables, the initial free cross-sectional area vector and the weight vector of the initial free unit length, and the initial free length vector of the supporting cables , the initial free cross-sectional area vector and the initial free unit length weight vector have the same _numbering rules as the elements of the initial cable force vector F _o ; At the same moment of the state temperature data vector T _o , the measured data of the initial cable structure are directly measured and calculated, and the measured data of the initial cable structure include the concentrated load measurement data of the cable structure, the distributed load measurement data of the cable structure, and the body load measurement data of the cable structure data, initial values of all monitored quantities, initial cable force data of all supporting cables, initial cable structure modal data, initial cable structure strain data, initial cable structure geometry data, initial cable structure support generalized coordinate data, initial cable structure The measured data including the angle data and the initial cable structure space coordinate data, the generalized coordinate data of the initial cable structure support include the space coordinate data of the initial cable structure support and the initial cable structure support angular coordinate data, after obtaining the measured data of the initial cable structure At the same time, the data that can express the health state of the support cable including the non-destructive testing data of the support cable are obtained through measurement and calculation. At this time, the data that can express the health state of the support cable is called the initial health state data of the support cable; all monitored The initial value of the quantity constitutes the initial value vector C _o of the monitored quantity, and the numbering rule of the initial value vector C _o of the monitored quantity is the same as that of the M monitored quantities; it is established by using the initial health state data of the supporting cable and the load measurement data of the cable structure The initial damage vector d _o of the evaluated object, the vector d _o represents the initial health state of the evaluated object of the cable structure represented by the initial mechanical calculation benchmark model A _o ; The number of elements of the damage vector d _o is equal to N, and the elements of d _o have a one-to-one correspondence with the evaluated object. The numbering rule of the elements of the vector d _o is the same as that of the evaluated object; if a certain element of d _o corresponds to The object to be evaluated is a supporting cable in the cable system, then the value of this element of d _o represents the initial damage degree of the corresponding supporting cable, if the value of this element is 0, it means that the supporting cable corresponding to this element is intact , without damage, if its value is 100%, it means that the supporting cable corresponding to this element has completely lost its bearing capacity; if its value is between 0 and 100%, it means that the supporting cable has lost its corresponding proportion of bearing capacity capacity; if a certain element of d _o corresponds to a certain load, the value of this element of d _o is taken as 0 in this method, which means that the initial value of the change of this load is 0; if there is no non-destructive testing of supporting cables data and other data that can express the health state of the supporting cable, or when the initial state of the structure can be considered as a state of no damage and no relaxation, the value of each element related to the supporting cable in the vector d _o is set to 0; the generalized coordinates of the initial cable structure support The data constitute the generalized coordinate vector U _o of the initial cable structure support;

d. According to the design drawing of the cable structure, the as-built drawing and the actual measurement data of the initial cable structure, the initial health status data of the supporting cable, the measurement data of the concentrated load of the cable structure, the measurement data of the distributed load of the cable structure, the measurement data of the body load of the cable structure, and the data of the cable structure. The physical and mechanical performance parameters of the various materials used vary with temperature, the generalized coordinate vector U _o of the initial cable structure support, the initial cable structure steady-state temperature data vector T _o and all the cable structure data obtained in the previous steps. The initial mechanical calculation benchmark model A _o of the cable structure is entered into the "steady-state temperature data of the cable structure". The calculated data of the cable structure based on A _o must be very close to the measured data, and the difference between them shall not be greater than 5%; corresponding to A _o The "cable structure steady-state temperature data" is the "initial cable structure steady-state temperature data vector T _o "; the cable structure support generalized coordinate data corresponding to A _o is the initial cable structure support generalized coordinate vector U _o ; corresponding to A The health status of the evaluated object of _o is represented by the initial damage vector d _o of the evaluated object; the initial values of all monitored quantities corresponding to A _o are represented by the initial value vector C _o of the monitored quantities; U _o , T _o and d _o are The parameters of A _o , the initial values of all the monitored quantities obtained from the mechanical calculation results of A _o are the same as the initial values of all the monitored quantities represented by C _o , so it can also be said that C _o is composed of the mechanical calculation results of A _o , In this method, A _o , C _o , d _o , U _o and T _o are unchanged;

e. In this method, except where the letter i clearly indicates the step number, the letter i only indicates the number of cycles, that is, the i-th cycle; the current initial value of the cable structure that needs to be established or has been established at the beginning of the i-th cycle The mechanical calculation benchmark model is recorded as the current initial mechanical calculation benchmark model A ⁱ _o , A _o and A ⁱ _o include temperature parameters, and the influence of temperature changes on the mechanical properties of the cable structure can be calculated; at the beginning of the i-th cycle, corresponding to A The "cable structure steady-state temperature data" of ⁱ _o is represented by the current initial cable structure steady-state temperature data vector T ⁱ _o , the definition of vector T ⁱ _o is the same as that of vector T _o , and the elements of T ⁱ _o are the same as T _o The elements of are in one-to-one correspondence; at the beginning of the i-th cycle, the "cable structure support generalized coordinate data" corresponding to A ⁱ _o is represented by the current initial cable structure support generalized coordinate vector U ⁱ _o , and the definition method of vector U ⁱ _o In the same way as the definition of the vector U _o , the elements of U ⁱ _o correspond to the elements of U _o one by one; the current initial damage vector of the evaluated object required at the beginning of the i-th cycle is recorded as d ⁱ _o , and d ⁱ _o represents the At the beginning of the cycle, the health status of the evaluated object of the search structure A ⁱ _o is defined. The definition of d ⁱ _o is the same as that of d _o , and the elements of d ⁱ _o correspond to the elements of d _o one by one; when the i-th cycle starts , the initial values of all the monitored quantities are represented by the current initial value vector C ⁱ _o of the monitored quantity. The definition method of the vector C ⁱ _o is the same as that of the vector C _o , and the elements of C ⁱ _o are one by one with the elements of C _o Correspondingly, the current initial value vector C ⁱ _o of the monitored quantity represents the specific values of all the monitored quantities corresponding to A ⁱ _o ; U ⁱ _o , T ⁱ _o and d ⁱ _o are the characteristic parameters of A ⁱ _o , and C ⁱ _o is determined by A ⁱ _o is composed of mechanical calculation results; at the beginning of the first cycle, A ⁱ _o is recorded as A ¹ _o , and the method of establishing A ¹ _o is to make A ¹ _o equal to A _o ; at the beginning of the first cycle, T ⁱ _o Denote it as T ¹ _o , the method of establishing T ¹ _o is to make T ¹ _o equal to T _o ; at the beginning of the first cycle, U ⁱ _o is recorded as U ¹ _o , and the method of establishing U ¹ _o is to make U ¹ _o equal to U _o ; at the beginning of the first cycle, d ⁱ _o is recorded as d ¹ _o , and the method of establishing d ¹ _o is to make d ¹ _o equal to d _o ; at the beginning of the first cycle, C ⁱ _o is recorded as C ¹ _o , and the establishment The method of C ¹ _o is to make C ¹ _o equal to C _o ;

f. From here, enter the cycle from step f to step q; during the service process of the structure, according to "the temperature measurement and calculation method of the cable structure of this method", the current data of the steady-state temperature data of the cable structure are continuously measured and calculated, and all The current data of "cable structure steady-state temperature data" constitutes the current cable structure steady-state temperature data vector T ⁱ , the definition method of vector T ⁱ is the same as that of vector T _o , and the elements of T ⁱ correspond to the elements of T _o one by one ; At the same moment when the steady-state temperature data vector T ⁱ of the current cable structure is obtained from the actual measurement, the current data of the generalized coordinates of the cable structure supports are obtained from the actual measurement, and all the current data of the generalized coordinates of the cable structure supports form the generalized coordinate vector U ⁱ , the vector U ⁱ is defined in the same way as the vector U _o ; when the vector T ⁱ is obtained from the actual measurement, all the cable structures in the cable structure at the same moment when the current cable structure steady-state temperature data vector T ⁱ is obtained The current value of the monitored quantity. All these values form the current value vector C ⁱ of the monitored quantity. The definition method of the vector C ⁱ is the same as that of the vector C _o . The elements of C ⁱ correspond to the elements of C _o one by one, indicating that the same The value of the monitoring quantity at different times; at the same moment when the current cable structure steady-state temperature data vector T ⁱ is measured, the measured cable force data of all M ₁ supporting cables in the cable structure, all these cable force data constitute the current cable force Vector F ⁱ , the elements of vector F ⁱ have the same numbering rules as the elements of vector F _o ; at the same moment when the current steady-state temperature data vector T ⁱ of the cable structure is obtained from the actual measurement, the two supporting cables of all M ₁ supporting cables are obtained from the actual measurement and calculation. The spatial coordinates of the end points, the difference between the spatial coordinates of the two support end points in the horizontal direction is the horizontal distance between the two support end points, the horizontal distance data of the two support end points of all the support cables form the horizontal distance vector between the two support end points of the current support cable, and the current support The numbering rule of the elements of the horizontal distance vector between the two supporting ends of the cable is the same as the numbering rule of the elements of the initial cable force vector F _o ;

g. According to the generalized coordinate vector U ⁱ of the measured support of the current cable structure and the current steady-state temperature data vector T ⁱ of the cable structure, follow steps g1 to g3 to update the current initial mechanical calculation benchmark model A ⁱ _o and the current initial value vector C of the monitored quantity ⁱ _o , the current initial cable structure steady-state temperature data vector T ⁱ _o and the current initial cable structure support generalized coordinate vector U ⁱ _o , while the current initial damage vector d ⁱ _o of the evaluated object remains unchanged;

g1. Compare T ⁱ and T ⁱ _o , U ⁱ and U ⁱ _o respectively, if T ⁱ is equal to T ⁱ _o and U ⁱ is equal to U ⁱ _o , there is no need to compare A ⁱ _o , U ⁱ _o , C ⁱ _o and T ⁱ _o to be updated, otherwise A ⁱ _o , U ⁱ _o , C ⁱ _o and T ⁱ _o need to be updated according to the following steps;

g2. Calculate the difference between U ⁱ and U _o . The difference between U ⁱ and U _o is the generalized displacement of the cable structure support relative to the initial position. Use the generalized displacement vector V of the support to represent the generalized displacement of the support. V is equal to U ⁱ minus Remove U _o ; calculate the difference between T ⁱ and T _o , the difference between T ⁱ and T _o is the change of the current cable structure steady-state temperature data with respect to the initial cable structure steady-state temperature data, the difference between T ⁱ and T _o is the steady-state temperature The change vector S indicates that S is equal to T ⁱ minus T _o , and S indicates the change of the steady-state temperature data of the cable structure;

g3. First apply the support generalized displacement constraint to the cable structure support in A _o , the value of the support generalized displacement constraint is taken from the value of the corresponding element in the support generalized displacement vector V, and then apply the support to the cable structure in A _o The temperature change, the value of the applied temperature change is taken from the steady-state temperature change vector S, and the current initial mechanical calculation benchmark model is updated after the generalized displacement constraint of the support is applied to the support of the cable structure in A _o and the temperature change is applied to the cable structure A ⁱ _o , while updating A ⁱ _o , all element values of U ⁱ _o are also replaced by corresponding values of all elements of U ⁱ , that is, U ⁱ _o is updated, and all element values of T ⁱ _o are also replaced by corresponding values of all elements of T ⁱ , that is, T ⁱ _o is updated, so that U ⁱ _o and T ⁱ _o corresponding to A ⁱ _o are obtained correctly, and d ⁱ _o remains unchanged at this time; when A ⁱ _o is updated, the index of A ⁱ _o The health status is expressed by the current initial damage vector d ⁱ _o of the evaluated object, the steady-state temperature of the cable structure of A ⁱ _o is expressed by the current data vector T ⁱ _o of the steady-state temperature of the cable structure, and the generalized coordinates of the support of A ⁱ _o are expressed by the current initial cable The generalized coordinate vector U ⁱ _o of the structural support is represented; the method of updating C ⁱ _o is: after updating A ⁱ _o , the current specific values of all monitored quantities in A ⁱ _o are obtained through mechanical calculations, and these specific values form C ⁱ _o ;

h. On the basis of the current initial mechanical calculation benchmark model A ⁱ _o , perform several mechanical calculations according to steps h1 to h4, and establish the numerical change matrix ΔC ⁱ of the unit damage monitored quantity and the unit change vector D ⁱ of the evaluated object through calculation _u ;

h1. At the beginning of the ith cycle, directly follow the methods listed in step h2 to step h4 to obtain ΔC ⁱ and D ⁱ _u ; at other times, after updating A ⁱ _o in step g, you must follow steps h2 to h4 The method listed in step h4 regains ΔC ⁱ and D ⁱ _u , if A ⁱ _o is not updated in step g, then directly transfer to step i for follow-up work here;

h2. Carry out several mechanical calculations on the basis of the current initial mechanical calculation benchmark model A ⁱ _o . The number of calculations is numerically equal to the number N of all evaluated objects. There are N evaluation objects and there are N calculations; Numbering rules, calculations are performed sequentially; each calculation assumes that there is only one evaluated object, and the unit damage or load unit change is added on the basis of the original damage or load. Specifically, if the evaluated object is a supporting cable in the cable system , then it is assumed that the support cable will increase the unit damage, if the evaluated object is a load, it is assumed that the load will increase the load unit change, and record the increased unit damage or load unit change with D ⁱ _uk , where k represents Increment the number of the assessed object of the unit damage or load unit change, D ⁱ _uk is an element of the assessed object unit change vector D ⁱ _u , the numbering rule of the elements of the assessed object unit change vector D ⁱ _u is the same as that of the vector d _o The numbering rules of the elements are the same; the evaluated object with additional unit damage or load unit change in each calculation is different from the evaluated object with additional unit damage or load unit change in other calculations, and each calculation uses the mechanical method to calculate the index The current calculated values of all the monitored quantities of the structure, the current calculated values of all the monitored quantities obtained by each calculation form a current vector of the monitored quantity calculation; when it is assumed that the kth assessed object is added with unit damage or load unit change , use C ⁱ _tk to represent the corresponding "monitored quantity calculation current vector"; when numbering the elements of each vector in this step, the same numbering rule should be used with other vectors in this method to ensure that the elements in each vector in this step Any element, with the elements of the same number in other vectors, expresses the relevant information of the same monitored quantity or the same object; the definition of C ⁱ _tk is the same as that of the vector C _o , and the elements of C ⁱ _tk are the same as C The elements of _o correspond one by one;

h3. The vector C ⁱ _tk obtained by each calculation is subtracted from the vector C ⁱ _o to obtain a vector, and then each element of the vector is divided by the value of unit damage or load unit change assumed in this calculation to obtain a "being The numerical change vector of the monitored quantity δC ⁱ _k "; if there are N evaluated objects, there are N "numerical change vectors of the monitored quantity";

h4. From these N "value change vectors of monitored quantities" according to the numbering rules of N evaluated objects, a "value change matrix ΔC ⁱ of unit damage monitored quantities" with N columns is formed in turn; the value of unit damage monitored quantities Each column of the change matrix ΔC ⁱ corresponds to a unit change vector of the monitored quantity; each row of the value change matrix ΔC ⁱ of the unit damage monitored quantity corresponds to the same monitored quantity when the unit damage or load unit changes for different evaluated objects Different unit change ranges; the numbering rules of the columns of the unit damage monitored quantity numerical change matrix ΔC ⁱ are the same as the numbering rules of the elements of the vector d _o , and the numbering rules of the rows of the unit damage monitored quantity numerical change matrix ΔC ⁱ are the same as M The numbering rules of the monitored quantities are the same;

i. Define the current nominal damage vector d ⁱ _c and the current actual damage vector d ⁱ , the number of elements of d ⁱ _c and d ⁱ is equal to the number of evaluated objects, and the distance between the elements of d ⁱ _c and d ⁱ and the evaluated object is One-to-one correspondence, the element value of d ⁱ _c represents the nominal damage degree or nominal load change of the corresponding evaluated object, d ⁱ _c and d ⁱ are the same as the element numbering rules of the initial damage vector d _o of the evaluated object, d ⁱ _c The elements of , the elements of d ⁱ and the elements of d _o are in one-to-one correspondence;

j. According to the relationship between the current numerical vector C ⁱ of the monitored quantity and the "current initial numerical vector C ⁱ _o of the monitored quantity", "the numerical change matrix ΔC ⁱ of the monitored quantity with unit damage" and "the current nominal damage vector d ⁱ _c " Approximate linear relationship, the approximate linear relationship can be expressed as formula 1. In formula 1, other quantities except d ⁱ _c are known, and the current nominal damage vector d ⁱ _c can be calculated by solving formula 1;

$C^i=C_o^i+{\mathrm{\ΔC}}^i\&Center\; Dot;d_c^i$ Formula 1

k. The kth element d ⁱ _k of the current actual damage vector d ⁱ expressed by formula 2 is the same as the kth element d i ok of the current initial damage vector d ⁱ _o of the evaluated object and the kth element d ⁱ _ok of the current nominal damage vector d ⁱ _c The relationship between k elements d ⁱ _ck is calculated to obtain all the elements of the current actual damage vector d ⁱ ;

In formula 2, ^k =1,2,3,...,N; the numbering rule of the elements of the _vector d ⁱ is the same as the numbering rule of the elements of the vector d _o in the formula (1); The current actual health status of k evaluated objects, if the evaluated object is a supporting cable in the cable system, then d ⁱ _k represents the severity of its current health problems, and the supporting cables with health problems may be slack cables, It may also be a damaged cable, and the d ⁱ _k value reflects the degree of relaxation or damage of the supporting cable; the M ₁ elements related to the M ₁ supporting cables in the current actual damage vector d ⁱ of the evaluated object are taken out to form a supporting The current actual damage vector d ^ci of the cable, the numbering rule of the elements of the current actual damage vector d ^ci of the supporting cable is the same as the numbering rule of the elements of the initial cable force vector F _o ; the hth element of the current actual damage vector d ^ci of the supporting cable represents the The current actual damage amount of the hth support cable in the structure, h=1, 2, 3,..., M ₁ ; the elements whose value is not 0 in the current actual damage vector d ^ci of the support cable correspond to supports with health problems The damaged cables are identified from these supporting cables with health problems, and the rest are slack cables; the value of the elements in the current actual damage vector d ^ci of the supporting cables corresponding to the damaged cables expresses When the element value is 100%, it means that the supporting cable completely loses its bearing capacity, and when it is between 0 and 100%, it means that the corresponding proportion of bearing capacity is lost; using the current cable structure steady-state temperature data vector T ⁱ condition The following, the slack cables identified in step l and the current actual equivalent damage degree of these slack cables expressed by the current actual damage vector d ^ci of the support cables, which are mechanically equivalent to their slack degrees, use the obtained in step f Under the current cable structure steady-state temperature data vector T ⁱ condition, the current cable force vector F ⁱ and the horizontal distance vector between the two supporting ends of the current supporting cable, use the initial cable structure steady-state temperature data vector T _o condition obtained in step c The initial free length vector of the supporting cable, the initial free cross-sectional area vector, the initial free weight vector per unit length, the initial cable force vector F _o , and the current steady state of the supporting cable represented by the current cable structure steady-state temperature data vector T ⁱ Temperature data, using the initial steady-state temperature data of the supporting cable represented by the initial cable structure steady-state temperature data vector T _o obtained in step c, using the temperature changes of various materials used in the cable structure obtained in step c Taking into account the influence of temperature changes on the physical, mechanical and geometric parameters of the supporting cable, the mechanical equivalent of the slack cable and the damaged cable is used to calculate the equivalent damage degree of the slack cable and the current actual equivalent damage degree. The mechanical equivalent conditions are: 1. The initial free length, geometric characteristic parameters, density and material mechanical characteristic parameters of the two equivalent cables are the same when there is no relaxation and no damage; 2. After relaxation or damage, the two The cable force and the total length after deformation of the equivalent slack cable and the damaged cable are the same; when the above two mechanical equivalence conditions are satisfied, such two supports The mechanical function of the cable in the cable structure is exactly the same, that is, if the damaged cable is replaced by an equivalent slack cable, the cable structure will not change, and vice versa; is the slack degree of the slack cable, the slack degree is the change of the free length of the support cable, that is, the adjustment amount of the cable length of the support cable whose force needs to be adjusted is determined; in this way, the slack identification and damage identification of the support cable are realized; the calculation The required cable force is given by the corresponding elements of the current cable force vector F ⁱ ; in this method, the damaged cable and the slack cable are collectively referred to as the support cable with health problems, referred to as the problem cable. , the influence of load change and structural temperature change, and the identification of problem cables in the cable structure, and at the same time realize the identification of the load change that eliminates the influence of generalized displacement of the support, structural temperature change and health state change of the support cable;

l. After obtaining the current nominal damage vector d ⁱ _c , establish the identification vector B ⁱ according to formula 3, and formula 4 gives the definition of the kth element of the identification vector B ⁱ ;

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

In formula 4, the element B ⁱ _k is the kth element of the identification vector B ⁱ , D ⁱ _uk is the kth element of the unit change vector D ⁱ _u of the evaluated object, d ⁱ _ck is the current nominal damage vector d ⁱ of the evaluated object The kth element of _c , they all represent the relevant information of the kth evaluated object, k=1,2,3,...,N in formula 4;

m. If the elements of the identification vector B ⁱ are all 0, then get back to step f to continue this cycle; if the elements of the identification vector B ⁱ are not all 0, then enter the next step, namely step n;

n. Calculate each element of the current initial damage vector d ⁱ⁺¹ _o of the evaluated object required for the next, ie i+1th cycle, according to formula 5;

In formula 5, d ⁱ⁺¹ _ok is the kth element of the current initial damage vector d ⁱ⁺¹ _o of the evaluated object required for the next cycle, that is, the i+1th cycle, and d ⁱ _ok is this time, that is, the ith The k-th element of the current initial damage vector d ⁱ _o of the evaluated object in the second cycle, D ⁱ _uk is the k-th element of the unit change vector D ⁱ _u of the evaluated object in the i-th cycle, and B ⁱ _k is the i-th The kth element of the cyclic identification vector B ⁱ , k=1,2,3,...,N in formula 5;

o. On the basis of the initial mechanical calculation benchmark model A _o , the support generalized displacement constraint is first imposed on the cable structure support in A _o , and the value of the support generalized displacement constraint is taken from the corresponding element in the support generalized displacement vector V value, and then apply a temperature change to the cable structure in A _o , the value of the applied temperature change is taken from the steady-state temperature change vector S, and then let the health status of the cable be d ⁱ⁺¹ _o to get the next time, That is, the mechanical calculation benchmark model A ⁱ⁺¹ required for the i+1th cycle; after obtaining A ⁱ⁺¹ , obtain the current specific values of all monitored quantities in A ⁱ⁺¹ through mechanical calculations, and these specific values consist of The current initial value vector C ⁱ⁺¹ _o of the monitored quantity required for the next time, that is, the i+1th cycle;

p. Take the current initial cable structure steady-state temperature data vector T ⁱ⁺¹ _o required for the next cycle, that is, the i+1th cycle, which is equal to the current initial cable structure steady-state temperature data vector T ⁱ _o of the i-th cycle; The current initial cable structure support generalized coordinate vector U ⁱ⁺¹ _o required for the i+1th cycle is equal to the current initial cable structure support generalized coordinate vector U ⁱ _o of the ith cycle;

q. Go back to step f and start the next cycle.

Beneficial effects: the structural health monitoring system first uses sensors to monitor the structural response on-line for a long time, and then analyzes the monitoring data online (or offline) to obtain the structural health status data. Due to the complexity of the structure, the structural health monitoring system needs to use A large number of sensors and other equipment are used for structural health monitoring, so the cost is usually quite high, so the cost problem is a major problem restricting the application of structural health monitoring technology. On the other hand, the correct identification of the health status of the core assessed objects (such as stay cables) is an integral part of, or even the whole of, the correct identification of the structural health status, while the secondary assessed objects (such as structural bearing The correct identification of changes in the load (such as changes in the number and mass of cars passing a cable-stayed bridge) has little or no effect on the correct identification of the state of health of the cable structure. However, the number of secondary evaluated objects is usually equal to the number of core evaluated objects, and the number of secondary evaluated objects is often greater than the number of core evaluated objects, so that the number of evaluated objects is often the number of core evaluated objects multiple times the number. When the secondary evaluated object (load) changes, in order to accurately identify the core evaluated object, the conventional method requires that the quantity of the monitored quantity (measured using sensors and other equipment) must be greater than or equal to the number of evaluated objects. When the number of secondary assessment objects is relatively large (in fact, it is often the case), the number of sensors and other equipment required by the structural health monitoring system is very large, so the cost of the structural health monitoring system will become very high, even high unacceptable. The inventors have found that the changes in the secondary evaluated objects (such as the normal load of the structure, the normal load of the structure means that the load that the structure is bearing does not exceed the structural allowable load defined in the structural design document or the structural completion document) are relatively small. Hours (For loads, the structure only bears normal loads. Whether the loads borne by the structure are normal loads can be determined by visual inspection. , the structure is only subjected to normal loads), the magnitude of the change in the structural response caused by them (this specification is called "secondary response") is much smaller than that caused by the change of the core assessed object (such as damage to the support cable) The magnitude of change in the response (this specification calls it "core response"), the sum of the secondary response and the core response is the total change in the structural response (this specification calls it the "overall response"), obviously the core response occupies the overall response Based on this, the inventor found that when determining the number of monitored quantities, even if a value slightly larger than the number of core evaluated objects is selected, but far smaller than the number of evaluated objects (this method is done in this way), that is to say, even if A relatively small number of sensors and other equipment can still accurately obtain the health status data of the core assessed objects and meet the core needs of structural health status monitoring. Therefore, the cost of the structural health monitoring system proposed by this method is obviously higher than that required by conventional methods. The cost of the structural health monitoring system is much lower, that is to say, this method can realize the evaluation of the health status of the core evaluated object of the cable structure at a much lower cost. Adoption is pivotal.

Detailed ways

The method employs an algorithm for identifying the health status of core assessed subjects. During specific implementation, the following steps are one of various steps that may be taken.

Step 1: First identify the amount of possible varying loads that the cable structure will be subjected to. According to the characteristics of the loads borne by the cable structure, confirm "all the loads that may change", or regard all the loads as "all the loads that may change", let there be a total of JZW loads that may change, that is, a total of JZW A secondary object to be evaluated. Let the sum of the number of supporting cables of the cable structure and the number of JZW "all possible changing loads" be N, that is, there are N evaluated objects. Consecutively number the evaluated objects, which will be used to generate vectors and matrices in subsequent steps.

The multi-type parameters to be monitored may include: cable force, strain, angle and spatial coordinates, which are described as follows:

Assuming that there are Q supporting cables in the cable system, that is, there are Q core objects to be evaluated, the monitored cable force data of the cable structure is described by M ₁ cable force data of M ₁ specified cables on the cable structure, and the cable structure cable The force change is the change in cable force for all specified cables. Each time there are M1 cable force measurements or calculations to represent the cable force information _of the cable structure. M ₁ is an integer not less than 0.

The monitored strain data of the cable structure can be described by the strains of K ₂ designated points on the cable structure and L ₂ designated directions of each designated point. The change of the strain data of the cable structure is all the K ₂ designated points Measure the change in strain. Each time there are M ₂ (M ₂ =K ₂ ×L ₂ ) measured or calculated strain values to characterize the strain of the cable structure. M ₂ is an integer not less than 0.

The monitored angle data of the cable structure is described by K ₃ designated points on the cable structure, L ₃ designated straight lines passing through each designated point, and H ₃ angular coordinate components of each designated straight line. The change is the change of all specified points, all specified lines, and all specified angular coordinate components. Each time there are M ₃ (M ₃ =K ₃ ×L ₃ ×H ₃ ) measured or calculated values of the angle coordinate components to represent the angle information of the cable structure. M ₃ is an integer not less than 0.

The monitored shape data of the cable structure is described by the spatial coordinates of K ₄ specified points on the cable structure and L ₄ specified directions of each specified point. The change of the shape data of the cable structure is all of the K ₄ specified points. Changes in coordinate components. Each time there are M ₄ (M ₄ =K ₄ ×L ₄ ) coordinate measurement values or calculation values to characterize the shape of the cable structure. M ₄ is an integer not less than 0.

Based on the above-mentioned monitored quantities, the entire cable structure has M (M=M ₁ +M ₂ +M ₃ +M ₄ ) monitored quantities, and the defined parameter K (K=M ₁ +K ₂ +K ₃ +K ₄ ), M must not be less than the number of core evaluated objects plus 4, and M is smaller than the number N of evaluated objects.

For convenience, in this method, "all monitored parameters of the cable structure" are referred to as "monitored quantities" for short. Number the M monitored quantities consecutively, and this number will be used to generate vectors and matrices in subsequent steps. The method uses variable j to denote this number, j=1,2,3,...,M.

Determine the "temperature measurement and calculation method of the cable structure of this method" according to the steps specified in the technical plan.

The second step: establish the initial mechanical calculation benchmark model A _o .

When the cable structure is completed, or before the health monitoring system is established, the "cable structure steady-state temperature data" is measured and calculated according to the "temperature measurement and calculation method of the cable structure of this method" (it can be measured by conventional temperature measurement methods, such as using Thermal resistance measurement), the "cable structure steady-state temperature data" at this time is represented by a vector T _o , which is called the initial cable structure steady-state temperature data vector T _o . When T _o is obtained from the actual measurement, the initial values of all monitored quantities of the cable structure are directly measured and calculated by conventional methods to form the initial value vector C _o of the monitored quantities.

In this method, the following method can be used to measure and calculate a certain measured quantity and monitored quantity (for example, a cable structure) at the same moment when a certain (such as initial or current, etc.) cable structure steady-state temperature data vector is obtained. All monitored quantities of the structure) data: while measuring and recording the temperature (including the air temperature of the environment where the cable structure is located, the temperature of the sunny side of the reference plate and the surface temperature of the cable structure), for example, measuring and recording the temperature every 10 minutes, then At the same time, the data of a certain measured quantity and monitored quantity (such as all monitored quantities of the cable structure) are also measured and recorded every 10 minutes. Once the time to obtain the steady-state temperature data of the cable structure is determined, the data of a certain measured quantity and monitored quantity (for example, all monitored quantities of the cable structure) at the same time as the time when the steady-state temperature data of the cable structure is obtained is called At the same moment when the steady-state temperature data of the cable structure is obtained, the data of a certain measured quantity and a certain monitored quantity obtained by a certain method are measured by a certain method.

The temperature-dependent physical parameters (such as thermal expansion coefficient) and mechanical performance parameters (such as elastic modulus and Poisson's ratio) of various materials used in the cable structure are obtained by conventional methods (research information or actual measurement).

While T _o is obtained from the actual measurement, the initial cable forces of all supporting cables are directly measured and calculated to form the initial cable force vector F _o ; according to the cable structure design data and completion data, the values of all supporting cables in the free state, that is, when the cable force is 0, are obtained. Length, cross-sectional area in the free state and weight per unit length in the free state, as well as the temperature of all supporting cables when these three data are obtained, based on which the physical performance parameters of all supporting cables that vary with temperature are used and mechanical performance parameters, according to conventional physical calculations, the length of all supporting cables when the cable force is 0 and the cross-sectional area of all supporting cables when the cable force is 0 under the condition of the initial cable structure steady-state temperature data vector T _o And the weight per unit length of all supporting cables when the cable force is 0, constitute the initial free length vector l _o of the supporting cables, the initial free cross-sectional area vector A _o and the weight vector ω _o of the initial free unit length of the supporting cables in turn, the initial The numbering rules of the elements of the free length vector l _o , the initial free cross-sectional area vector A _o and the initial free unit length weight vector ω _o are the same as the numbering rules of the elements of the initial cable force vector F _o .

According to the steps specified in the technical plan, while T _o is obtained through actual measurement, the actual measurement and calculation data of the cable structure are obtained through actual measurement and calculation using conventional methods. Utilize the design drawing, as-built drawing of the cable structure, the measured data of the initial cable structure, the non-destructive testing data of the supporting cable, the physical and mechanical performance parameters of various materials used in the cable structure as a function of temperature, and the generalized coordinates of the initial cable structure support The vector U _o and the initial cable structure steady-state temperature data vector T _o are incorporated into the "cable structure steady-state temperature data" using mechanical methods (such as finite element method) to establish the initial mechanical calculation benchmark model A _o . T _o , U _o and d _o are parameters of A _o , and C _o is composed of mechanical calculation results of A _o .

The third step: in this method, except the place where the letter i clearly represents the step number, the letter i only represents the number of cycles, i.e. the ith cycle; The current initial mechanical calculation benchmark model is denoted as the current initial mechanical calculation benchmark model A ⁱ _o , A _o and A ⁱ _o include temperature parameters, and the influence of temperature changes on the mechanical properties of the cable structure can be calculated; at the beginning of the i-th cycle, the corresponding The "cable structure steady-state temperature data" in A ⁱ _o is represented by the current initial cable structure steady-state temperature data vector T ⁱ _o , the definition of vector T ⁱ _o is the same as that of vector T _o , and the elements of T ⁱ _o are the same as The elements of T _o are in one-to-one correspondence; the generalized coordinate data of the cable structure support required at the beginning of the i-th cycle and corresponding to the current initial mechanical calculation benchmark model A ⁱ _o of the cable structure constitute the current initial cable structure support generalized coordinate vector U ⁱ _o , when the current initial mechanical calculation benchmark model A ⁱ _o of the cable structure is established for the first time, U ⁱ _o is equal to U _o . The current initial damage vector of the evaluated object required at the beginning of the i-th cycle is denoted as d ⁱ _o , d ⁱ _o represents the health status of the evaluated object of the cable structure A ⁱ _o at the beginning of the cycle, and the definition of d ⁱ _o is the same as d _o is defined in the same way, and the elements of d ⁱ _o correspond to the elements of d _o one by one; at the beginning of the i-th cycle, the initial values of all monitored quantities are represented by the current initial value vector C ⁱ _o of the monitored quantity, and the vector The definition method of C ⁱ _o is the same as the definition method of vector C _o , the elements of C ⁱ _o correspond to the elements of C _o one by one, and the current initial value vector C ⁱ _o of the monitored quantity represents all the monitored quantities corresponding to A ⁱ _o The specific value of ; T ⁱ _o and d ⁱ _o are the characteristic parameters of A ⁱ _o ; C ⁱ _o is composed of the mechanical calculation results of A ⁱ _o ; at the beginning of the first cycle, A ⁱ _o is recorded as A ¹ _o , and A The method of ¹ _o is to make A ¹ _o equal to A _o ; at the beginning of the first cycle, T ⁱ _o is recorded as T ¹ _o , and the method of establishing T ¹ _o is to make T ¹ _o equal to T _o ; at the beginning of the first cycle , U ⁱ _o is recorded as U ¹ _o , the method of establishing U ¹ _o is to make U ¹ _o equal to U _o ; at the beginning of the first cycle, d ⁱ _o is recorded as d ¹ _o , and the method of establishing d ¹ _o is to make d ¹ _o is equal to d _o ; at the beginning of the first cycle, C ⁱ _o is recorded as C ¹ _o , and the method of establishing C ¹ _o is to make C ¹ _o equal to C _o .

Step 4: Install the hardware part of the cable structure health monitoring system. The hardware part includes at least: the monitored quantity monitoring system (such as angle measurement subsystem, cable force measurement subsystem, strain measurement subsystem, space coordinate measurement subsystem, signal conditioner, etc.), cable structure support generalized coordinate monitoring system ( Including total station, angle measurement sensor, signal conditioner, etc.), cable structure temperature monitoring system (including temperature sensor, signal conditioner, etc.) and cable structure environment temperature measurement system (including temperature sensor, signal conditioner, etc.), supporting cable A force monitoring system, a space coordinate monitoring system of the supporting end point of the supporting cable, a signal (data) collector, a computer and a communication alarm device. Each monitored quantity, the generalized coordinates of each support of the cable structure, each temperature, the cable force of each supporting cable, and the spatial coordinates of the supporting end points of each supporting cable must be monitored by the monitoring system. The monitored signal is transmitted to the signal (data) collector; the signal is transmitted to the computer through the signal collector; the computer is responsible for running the health monitoring software of the evaluated object of the cable structure, including recording the signal transmitted by the signal collector; when monitoring When there is a change in the health status of the evaluated object, the computer controls the communication alarm device to alarm the monitoring personnel, the owner and (or) the designated personnel.

The 5th step: compile and install and run the system software of this method on the computer, this software will finish the functions such as monitoring, record, control, storage, computing, notification, alarm that this method task needs (that is all in this specific implementation method work that can be done with a computer).

Step 6: This step starts the cycle operation. During the service process of the structure, the current data of the steady-state temperature data of the cable structure are obtained through continuous actual measurement and calculation according to the "temperature measurement and calculation method of the cable structure of this method". All "cable structure steady-state The current data of "Temperature Data" form the current cable structure steady-state temperature data vector T ⁱ , the definition of vector T ⁱ is the same as that of vector T _o , and the elements of T ⁱ correspond to the elements of T _o one by one; in the measured vector T At the same time as ⁱ , the current values of all the monitored quantities in the cable structure are actually measured, and all these values form the current value vector C ⁱ of the monitored quantity. The definition of the vector C ⁱ is the same as that of the vector C _o . The elements of C ⁱ are the same as The elements of C _o correspond to each other, indicating the value of the same monitored quantity at different times.

At the same time as the measured vector T ⁱ , the current data of the generalized coordinates of the cable structure support are obtained from the actual measurement, and all the data form the current generalized coordinate vector U ⁱ of the measured support of the cable structure.

At the same time as the actual measurement of the vector T ⁱ , the cable force data of all M ₁ supporting cables in the cable structure are obtained by actual measurement, all these cable force data form the current cable force vector F ⁱ , the numbers of the elements of the vector F ⁱ and the elements of the vector F _o The rules are the same; at the same time as the measured vector T ⁱ , the measured and calculated spatial coordinates of the two supporting end points _of all M1 supporting cables, the difference in the horizontal component of the spatial coordinates of the two supporting end points is the horizontal distance between the two supporting end points, The horizontal distance data of the two supporting ends of all M ₁ supporting cables constitute the horizontal distance vector l ⁱ _x of the two supporting ends of the current supporting cables, the numbering rules of the elements of the horizontal distance vector l ⁱ _x between the two supporting ends of the current supporting cables and the initial cable force vector The numbering rules for the elements of F _o are the same.

Step 7: After obtaining the generalized coordinate vector U ⁱ of the measured support of the current cable structure and the steady-state temperature data vector T ⁱ of the current cable structure, compare U ⁱ and U ⁱ _o , T ⁱ and T ⁱ _o respectively, if U ⁱ is equal to U ⁱ _o and T ⁱ is equal to T ⁱ _o , then there is no need to update A ⁱ _o , U ⁱ _o , C ⁱ _o and T ⁱ _o , otherwise it is necessary to update A ⁱ _o , U ⁱ _o , C ⁱ _o and T ⁱ _o is updated, while the current initial damage vector d ⁱ _o of the evaluated object remains unchanged, and the update method is carried out according to the steps specified in the technical plan.

Step 8: According to the steps specified in the technical plan, on the basis of the current initial mechanical calculation benchmark model A ⁱ _o , perform several mechanical calculations, and establish the numerical change matrix ΔC ⁱ of the monitored quantity of unit damage and the unit change vector of the evaluated object through calculation D ⁱ _u . Specifically, if the evaluated object is a supporting cable in the cable system, then it is assumed that the supporting cable has unit damage on the basis of the existing damage of the supporting cable represented by the vector d ⁱ _o (such as taking 5%, 10%, 20% or 30% damage is the unit damage), if the object to be evaluated is a load, it is assumed that the load will increase the load unit change on the basis of the existing change of the load represented by the vector d ⁱ _o ( If the load is a distributed load, and the distributed load is a linear distributed load, the load unit can be changed in units of 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m, etc.; if The load is a distributed load, and the distributed load is a surface distributed load, and the load unit can be changed in units of 1MPa, 2MPa, 3MPa or 1kNm/m ² , 2kNm/m ² , 3kNm/m ² ; if the load is Concentrated load, and the concentrated load is a force couple, the load unit can be changed in units of 1kNm, 2kNm, 3kNm, etc.; if the load is a concentrated load, and the concentrated load is a concentrated force, the load unit can be changed in 1kN, 2kN, 3kN etc.; if the load is body load, the load unit change can take 1kN/m ³ , 2kN/m ³ , 3kN/m ³ etc. as the unit change).

Step 9: Establish linear relationship error vector e ⁱ and vector g ⁱ . Using the previous data ("the current initial value vector C ⁱ _o of the monitored quantity", "the numerical change matrix of the monitored quantity per unit damage ΔC ⁱ "), each calculation is performed in the eighth step, that is, each calculation assumes that When only one of the assessed objects has increased unit damage or load unit change, when it is assumed that the kth (k=1,2,3,...,N)th assessed object has increased unit damage or load unit change, Each calculation forms a damage vector, which is represented by d ⁱ _tk , and the current vector of the corresponding monitored quantity calculation is C ⁱ _tk (see step 8), and the number of elements of the damage vector d ⁱ _tk is equal to the value of the evaluated object Quantity, the value of only one element among all the elements of vector d ⁱ _tk is the unit damage or load unit change value of the evaluated object assuming an increase in unit damage or load unit change in each calculation, and the values of other elements of d ⁱ _tk are 0, the corresponding relationship between the number of the element that is not 0 and the evaluated object assuming the increase of unit damage or the change of the load unit, and the corresponding relationship between the elements with the same number of other vectors and the evaluated object are the same; d ⁱ _tk The element numbering rule is the same as the element numbering rule of the initial damage vector d _o of the evaluated object, and the elements of d ⁱ _tk are in one-to-one correspondence with the elements of d _o . Bring C ⁱ _tk , C ⁱ _o , ΔC ⁱ , d ⁱ _tk into formula (1), and obtain a linear relationship error vector e ⁱ _k , and each calculation obtains a linear relationship error vector e ⁱ _k ; the next step of e ⁱ _k The mark k indicates that the kth (k=1,2,3,...,N) evaluated object increases unit damage or load unit change. If there are N evaluated objects, there will be N calculations, and there will be N linear relationship error vectors e ⁱ _k , and the N linear relationship error vectors e ⁱ _k will be added to obtain a vector, and each element of this vector will be divided by The new vector obtained after N is the final linear relationship error vector e ⁱ . The vector g ⁱ is equal to the final error vector e ⁱ . Save the vector g ⁱ on the computer hard disk running the health monitoring system software for use by the health monitoring system software.

${\mathrm{<span\; class="notranslate">\; e</span>}}_{\mathrm{<span\; class="notranslate">\; k</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; b</span>}\mathrm{<span\; class="notranslate">\; the\; s</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&Delta;C</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&CenterDot;</span>{\mathrm{<span\; class="notranslate">\; d</span>}}_{\mathrm{<span\; class="notranslate">\; t</span>}\mathrm{<span\; class="notranslate">\; k</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; C</span>}}_{\mathrm{<span\; class="notranslate">\; t</span>}\mathrm{<span\; class="notranslate">\; k</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; C</span>}}_{\mathrm{<span\; class="notranslate">\; o</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; )</span>$

Step 10: Define the current nominal damage vector d ⁱ _c and the current actual damage vector d ⁱ , the number of elements of d ⁱ _c and d ⁱ is equal to the number of evaluated objects, and the number of elements of d ⁱ _c and d ⁱ is equal to the number of evaluated objects There is a one-to-one correspondence between d ⁱ _c and d ⁱ , the element values of d i c and d i represent the damage degree or load change degree of the corresponding evaluated object, d ⁱ _c and d ⁱ are the same as the element numbering rules of the initial damage vector d _o of the evaluated object, The elements of d ⁱ _c , the elements of d ⁱ and the elements of d _o are in one-to-one correspondence.

Step 11: According to the current numerical vector C ⁱ of the monitored quantity, the "current initial numerical vector C ⁱ _o of the monitored quantity", "the numerical change matrix of the monitored quantity with unit damage ΔC ⁱ ", and "the current nominal damage vector d ⁱ _c " The approximate linear relationship exists between , the approximate linear relationship can be expressed as formula (2), and the non-inferior solution of the current nominal damage vector d ⁱ _c is calculated according to the multi-objective optimization algorithm, that is, there is a reasonable error, but it can be obtained from A solution that determines the location of damaged cables and their nominal damage levels in all cables.

${\mathrm{<span\; class="notranslate">\; C</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; C</span>}}_{\mathrm{<span\; class="notranslate">\; o</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; \&Delta;C</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&CenterDot;</span>{\mathrm{<span\; class="notranslate">\; d</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; )</span>$

The current nominal damage vector d ⁱ _c can be obtained by solving Equation (2) by using the Goal Attainment Method in the multi-objective optimization algorithm.

Step 12: Calculate each element of the current actual damage vector d ⁱ according to the definition of the current actual damage vector d ⁱ of the cable system and the definition of its elements, so that the health status of the evaluated object can be determined by d ⁱ . The k-th element d ⁱ _k of the current actual damage vector d ⁱ represents the current actual health status of the k-th evaluated object in the i-th cycle.

d ⁱ _k represents the current actual health status of the kth evaluated object in the i-th cycle, if the evaluated object is a supporting cable in the cable system, then d ⁱ _k represents its current actual damage, d ⁱ _k is When 0, it means that the corresponding supporting cable has no health problems. When the d ⁱ _k value is not 0, it means that the corresponding supporting cable has health problems. The supporting cables with health problems may be loose or damaged. cable, whose value reflects the degree of laxity or damage.

Take out the Q elements related to the supporting cable from the current actual damage vector d ⁱ of the evaluated object to form the current actual damage vector d ^ci of the supporting cable, and the numbering rules of the elements of the current actual damage vector d ^ci of the supporting cable are the same as the initial cable force vector F The numbering rules for the elements of _o are the same. The hth element of the current actual damage vector d ^ci of the support cable represents the current actual damage amount of the h-th support cable in the cable structure, h=1,2,3,...,Q; the current actual damage vector d ^ci of the support cable The elements whose values are not 0 correspond to the supporting cables with health problems, and the damaged cables are identified from these supporting cables with health problems by non-destructive testing methods, and the rest are slack cables, which are the cables whose force needs to be adjusted. The values of the elements corresponding to the current actual damage vector d ^ci of the cables whose forces need to be adjusted (for example, one of the elements can be represented by d ^ci _h ) represent the mechanically equivalent damage degrees to the slack of these support cables, thus The slackline is identified. The value of the element in the current actual damage vector d ^ci of the supporting cable corresponding to the damaged cable expresses the current actual damage of the damaged cable. When the value of the element is 100%, it means that the supporting cable completely loses its bearing capacity, between 0 and Between 100% means that the corresponding proportion of bearing capacity is lost.

Taking into account the influence of temperature changes on the physical, mechanical and geometric parameters of the supporting cable, the relaxation degree of the slack cable, which is equivalent to the current actual equivalent damage degree, can be calculated by mechanically equivalenting the slack cable to the damaged cable. Specifically, it can be The degree of relaxation of these cables (that is, the adjustment amount of the cable length) is obtained according to formula (3). This enables slack detection of the support cables. Damaged and slack cables have all been identified so far.

${\mathrm{<span\; class="notranslate">\; \&Delta;l</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; d</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; c</span>}\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; d</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; c</span>}\mathrm{<span\; class="notranslate">\; i</span>}}}\frac{{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}}{<span\; class="notranslate">\; \&lsqb;</span>\frac{{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; +</span>\frac{{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&omega;</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; l</span>}}_{\mathrm{<span\; class="notranslate">\; x</span>}\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 2</span>}}{\mathrm{<span\; class="notranslate">\; A</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}{\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; 12</span>}{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 3</span>}}}}<span\; class="notranslate">\; \&rsqb;</span>{\mathrm{<span\; class="notranslate">\; A</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; l</span>}}_{\mathrm{<span\; class="notranslate">\; o</span>}\mathrm{<span\; class="notranslate">\; h</span>}}^{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 3</span>}<span\; class="notranslate">\; )</span>$

In formula (3), E ⁱ _h is the elastic modulus of the hth supporting cable when the steady-state temperature data of the cable structure is represented by the current initial cable structure steady-state temperature data vector T ⁱ _o , A ⁱ _h is the elastic modulus of the cable structure When the steady-state temperature data of the current initial cable structure steady-state temperature data vector T ⁱ _o is represented, the cross-sectional area of the hth supporting cable, F ⁱ _h is the steady-state temperature data of the cable structure and the current initial cable structure steady-state When the temperature data vector T ⁱ _o represents the current cable force of the h-th supporting cable, d ^ci _h is the current actual damage degree of the h-th supporting cable, ω ⁱ _h is the steady-state temperature data of the cable structure using the current initial cable When the structural steady-state temperature data vector T ⁱ _o represents the weight per unit length of the hth supporting cable, l ⁱ _xh is when the steady-state temperature data of the cable structure is represented by the current initial cable structure steady-state temperature data vector T ⁱ _o , the horizontal distance between the two supporting endpoints of the h-th supporting cable, l ⁱ _xh is an element of the horizontal distance vector l ⁱ _x between the two supporting endpoints of the current supporting cable, and the element of the horizontal distance vector l ⁱ _x between the two supporting endpoints of the current supporting cable The numbering rules are the same as the numbering rules of the elements of the initial free length vector l _o , E ⁱ _h can be obtained according to the material property data of the h-th support cable, A ⁱ _h and ω ⁱ _h can be obtained according to the h-th support cable The coefficients of thermal expansion, A _oh , ω _oh , F ⁱ _h , T _o and T ⁱ _o are calculated by conventional physics and mechanics.

Step 13: The computer in the health monitoring system generates reports on the health status of the cable system automatically or by personnel operating the health monitoring system on a regular basis.

Step 14: Under the specified conditions, the computer in the health monitoring system automatically operates the communication alarm equipment to alarm the monitoring personnel, the owner and (or) the designated personnel.

Step 15: Establish the identification vector B ⁱ , if the elements of the identification vector B ⁱ are all 0, go back to the sixth step to continue the health monitoring and calculation of the cable system; if the elements of the identification vector B ⁱ are not all 0, After completing the subsequent steps, enter the next cycle.

Step 16: Calculate each element d ⁱ +1 o of the initial damage vector d ⁱ⁺¹ _o required for the next (i+ _1th , i=1,2,3,4,…) cycle (k=1,2,3,……,N); Second, on the basis of the initial mechanical calculation benchmark model A _o , the support generalized displacement constraint is imposed on the cable structure support in A _o , the support generalized The value of the displacement constraint is obtained from the value of the corresponding element in the generalized displacement vector V of the support, and then the temperature change is applied to the cable structure in A _o , and the value of the temperature change is obtained from the steady-state temperature change vector S, and then the cable The health status of d ⁱ⁺¹ _o is the mechanical calculation benchmark model A ⁱ⁺¹ required for the next cycle, that is, the i+1th (i=1,2,3,4,...) cycle; the next time (i.e. the i+1th time, i=1,2,3,4,...) The current initial cable structure steady-state temperature data vector T ⁱ⁺¹ _o required for the cycle is equal to T ⁱ _o , the next time (i.e. the i+th 1 time, i=1, 2, 3, 4,...) The current initial cable structure support generalized coordinate vector U ⁱ⁺¹ _o required by the cycle is equal to U ⁱ _o . After obtaining A ⁱ⁺¹ , d ⁱ⁺¹ _o , U ⁱ⁺¹ _o and T ⁱ⁺¹ _o , the current specific values of all the monitored quantities in A ⁱ⁺¹ are obtained through mechanical calculations, and these specific values are composed of the following The current initial value vector C ⁱ⁺¹ _o of the monitored quantity required for one cycle, that is, the i+1th cycle.

Step 17: Go back to Step 6 and start the cycle from Step 6 to Step 17.

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.