CN105004561A  Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement  Google Patents
Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement Download PDFInfo
 Publication number
 CN105004561A CN105004561A CN201510439869.0A CN201510439869A CN105004561A CN 105004561 A CN105004561 A CN 105004561A CN 201510439869 A CN201510439869 A CN 201510439869A CN 105004561 A CN105004561 A CN 105004561A
 Authority
 CN
 China
 Prior art keywords
 cable
 temperature
 data
 load
 vector
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Pending
Links
 238000006073 displacement reaction Methods 0.000 title claims abstract description 49
 238000000034 method Methods 0.000 title claims abstract description 11
 238000004364 calculation method Methods 0.000 claims abstract description 104
 230000003862 health status Effects 0.000 claims abstract description 34
 239000011159 matrix material Substances 0.000 claims abstract description 8
 238000011156 evaluation Methods 0.000 claims description 209
 238000009826 distribution Methods 0.000 claims description 93
 238000005259 measurement Methods 0.000 claims description 71
 230000000875 corresponding Effects 0.000 claims description 40
 230000001131 transforming Effects 0.000 claims description 29
 239000000203 mixture Substances 0.000 claims description 17
 239000000463 material Substances 0.000 claims description 14
 238000009659 nondestructive testing Methods 0.000 claims description 8
 206010022114 Injury Diseases 0.000 claims description 7
 229910052799 carbon Inorganic materials 0.000 claims description 7
 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 6
 240000002853 Nelumbo nucifera Species 0.000 claims description 6
 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 6
 230000035852 Tmax Effects 0.000 claims description 4
 241001269238 Data Species 0.000 claims description 3
 238000004458 analytical method Methods 0.000 claims description 3
 125000004122 cyclic group Chemical group 0.000 claims description 3
 230000001771 impaired Effects 0.000 claims description 3
 229910000975 Carbon steel Inorganic materials 0.000 claims description 2
 239000010962 carbon steel Substances 0.000 claims description 2
 239000002131 composite material Substances 0.000 claims description 2
 238000009413 insulation Methods 0.000 claims description 2
 238000004519 manufacturing process Methods 0.000 claims description 2
 239000010902 straw Substances 0.000 claims description 2
 238000000547 structure data Methods 0.000 claims description 2
 239000002965 rope Substances 0.000 description 11
 238000007796 conventional method Methods 0.000 description 5
 238000004891 communication Methods 0.000 description 3
 239000000725 suspension Substances 0.000 description 3
 230000005540 biological transmission Effects 0.000 description 2
 238000007906 compression Methods 0.000 description 2
 238000005516 engineering process Methods 0.000 description 2
 238000005457 optimization Methods 0.000 description 2
 230000015572 biosynthetic process Effects 0.000 description 1
 238000009529 body temperature measurement Methods 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000005755 formation reaction Methods 0.000 description 1
 230000005484 gravity Effects 0.000 description 1
 239000004615 ingredient Substances 0.000 description 1
 238000009434 installation Methods 0.000 description 1
 238000003860 storage Methods 0.000 description 1
 239000000700 tracer Substances 0.000 description 1
Abstract
The invention provides a method for recognizing the load of a damaged cable based on the space coordinate monitoring process of the streamlined annular displacement. According to the method, based on the space coordinate monitoring process, whether a mechanical calculation benchmark model of a cable structure needs to be updated or not is determined through monitoring the annular displacement of a bearer, the temperature of the cable structure and the ambient temperature. After that, the mechanical calculation benchmark model of the cable structure, with the annular displacement of the bearer, the temperature of the cable structure and the ambient temperature being taken into account, is obtained. On the basis of the above model, a numerical value variation matrix of a monitored quantity with unit damage is obtained through calculation. Finally, according to the approximately linear relationships between the current numeric vector of the monitored quantity and the current initial numeric vector of the monitored quantity, the numerical value variation matrix of the monitored quantity with unit damage and the current nominal damage vector of an unknown tobeevaluated object, the noninferior solution of the current nominal damage vector of the tobeevaluated object is figured out. Therefore, the health status of a core tobeevaluated object can be recognized.
Description
Technical field
Cablestayed bridge, suspension bridge, the structures such as trussframe structure have a common ground, be exactly that they have many parts bearing tensile load, as suspension cable, main pushtowing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or the rod member only bearing tensile load are support unit, for simplicity, such structure representation is " Cable Structure " by this method, and by all ropeway carryingropes of Cable Structure, carrying cable, and all rod members (being also called two power rod members) only bearing axial tension or axial compression load, be collectively referred to as simplicity " cable system ", ropeway carryingrope is censured with " support cable " this noun in this method, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word below, two power rod members are just referred to trussframe structure reality.In structure military service process, the correct identification of the health status of support cable or cable system is related to the safety of whole Cable Structure.When environment temperature changes, the temperature of Cable Structure generally also can along with changing, when Cable Structure temperature changes, may angular displacement be there is in Cable Structure bearing, the load that Cable Structure is born also may change, the health status of Cable Structure also may change simultaneously, at this complex condition, this method identifies damaged cable (this method is referred to as the health status of core evaluation object) based on space coordinate monitoring (monitored volume coordinate is called " monitored amount " by this method), belongs to engineering structure health monitoring field.
Background technology
Reject the impact of load change, Cable Structure angular displacement of support and structure temperature change on Cable Structure health status recognition result, thus identify the change of the health status of structure exactly, be problem in the urgent need to address at present, this method discloses a kind of effective, the cheap method addressed this problem.
Summary of the invention
Technical matters: this method discloses a kind of method, under the condition that cost is lower, when bearing has angular displacement, when the load that structure is born and structure (environment) temperature variation, the impact of angular displacement of support, load change and structure temperature change on Cable Structure health status recognition result can be rejected, thus identify the health status of support cable exactly.
Technical scheme: in the method, the coordinate of bearing about the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing volume coordinate ", alternatively becoming is the volume coordinate of bearing about X, Y, Z axis, bearing is called the volume coordinate component of bearing about this axle about the concrete numerical value of the volume coordinate of some axles, and a volume coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the volume coordinate of some axles; The angular coordinate of bearing about X, Y, Z axis is censured with " bearing angular coordinate ", bearing is called the angular coordinate component of bearing about this axle about the concrete numerical value of the angular coordinate of some axles, and an angular coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the angular coordinate of some axles; Censure bearing angular coordinate and bearing volume coordinate entirety with " bearing generalized coordinate ", a generalized coordinate component also with bearing in this method expresses the volume coordinate of bearing about an axle or the concrete numerical value of angular coordinate; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, and alternatively the change of bearing volume coordinate is called support wire displacement, and a translational component also with bearing in this method expresses the concrete numerical value of bearing about the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, and an angular displacement component also with bearing in this method expresses the concrete numerical value of bearing about the angular displacement of some axles; Generalized displacement of support censures support wire displacement and angular displacement of support is all, and a generalized displacement component also with bearing in this method expresses the displacement of the lines of bearing about some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is support wire displacement or the translational displacement component at gravity direction.
The external force that object, structure are born can be described as load, and load comprises face load and volume load.Face load, also known as surface load, is the load acting on body surface, comprises centrepoint load and distributed load two kinds.Volume load be continuous distribution in the load of interior of articles each point, as deadweight and the inertial force of object.
Centrepoint load is divided into concentrated force and concentrated couple two kinds, in a coordinate system, such as in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centrepoint load, in the method a concentrated force component or a concentrated couple component are called a load, the now change of load is embodied as the change of a concentrated force component or a concentrated couple component.
Distributed load is divided into line distributed load and EDS maps load, the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, the size distribution intensity of distributed load is expressed, distribution intensity distribution characteristics is (such as uniform, sine function equal distribution feature) and amplitude is expressed, and (such as two distributed loads are all uniform, but its amplitude is different, can welldistributed pressure be example so that the concept of amplitude to be described: same structure bears two different welldistributed pressures, two distributed loads are all uniformly distributed loads, but the amplitude of a distributed load is 10MPa, the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, in fact refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, a distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is all not identical, so in the method the component of these several distributed loads is regarded as the independently distributed load of same quantity, now load just represents the component of a distributed load, also component identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized a distributed load or be called a load.
Volume load is that continuous distribution is in the load of interior of articles each point, as deadweight and the inertial force of object, the description of volume load at least comprises the zone of action of volume load and the size of volume load, the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics is (such as uniform, linear function equal distribution feature) and amplitude is expressed, and (such as two individual stow lotuses are all uniform, but its amplitude is different, can conduct oneself with dignity for example is to illustrate the concept of amplitude: the material of two parts of same structure is different, therefore density is different, so although this volume load suffered by two parts is all uniform, but the amplitude of the volume load suffered by a part may be 10kN/m
^{3}, the amplitude of the volume load suffered by another part is 50kN/m
^{3}).If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes.In a coordinate system, one individual stow lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles about coordinate system, that is, in Descartes's rectangular coordinate system, volume load can resolve into three components), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is all not identical, so in the method the component of this several body stow lotus is regarded as the independently load of same quantity, also the volume sharing part of the load identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized an individual stow lotus or be called a load.
When load is embodied as centrepoint load, in the method, " load unit change " in fact refers to " unit change of centrepoint load ", similar, " load change " specifically refers to " change of the size of centrepoint load ", " load change amount " specifically refers to " variable quantity of the size of centrepoint load ", " load change degree " specifically refers to " intensity of variation of the size of centrepoint load ", " the actual change amount of load " refers to " the actual change amount of the size of centrepoint load ", " load changed " refers to " centrepoint load that size changes ", briefly, now " soandso load soandso change " refers to " size of soandso centrepoint load soandso change ".
When load is embodied as distributed load, in the method, " load unit change " in fact refers to " unit change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, similar, " load change " specifically refers to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " specifically refers to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " specifically refers to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " specifically refers to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that changes of amplitude of distribution intensity ", briefly, now " soandso load soandso change " refers to " amplitude of the distribution intensity of soandso distributed load soandso change ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.
When load is embodied as volume load, in the method, " load unit change " in fact refers to " unit change of the amplitude of the distribution intensity of volume load ", similar, " load change " refers to " change of the amplitude of the distribution intensity of volume load ", " load change amount " refers to " variable quantity of the amplitude of the distribution intensity of volume load ", " load change degree " refers to " intensity of variation of the amplitude of the distribution intensity of volume load ", " the actual change amount of load " refers to " the actual change amount of the amplitude of the distribution intensity of volume load ", " load changed " refers to " the volume load that changes of amplitude of distribution intensity ", briefly, " soandso load soandso change " refers to " amplitude of the distribution intensity of soandso volume load soandso change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.
This method specifically comprises:
Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated support cable and load to be " evaluation object ", if the quantity sum of the quantity of evaluated support cable and load is N, namely the quantity of " evaluation object " is N; This method title " core evaluation object " specially refers to the evaluated support cable in " evaluation object ", and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; If total M in cable system
_{1}root support cable, the quantity of obvious core evaluation object is exactly M
_{1}; Determine to specify by the measured point of monitored volume coordinate, number to all specified points; Determined each measurement point by monitored volume coordinate component, to all measured volume coordinate component numberings; Abovementioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored spatial data of Cable Structure " is made up of abovementioned all measured volume coordinate components; For simplicity, in the method by " the monitored spatial data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M is greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount RealTime Monitoring in this method, the moment of survey record data is called the physical record data moment; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centrepoint load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centrepoint load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centrepoint load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated 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 comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact 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; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the onecomponent of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, only measure the temperature of the surface point of support cable, measure all temperature be selected a little, the temperature recorded is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained the temperature of the sunny slope of reference plate,
B2: RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, reference plate maximum temperature difference Δ T
_{pmax}with Cable Structure surface maximum temperature difference Δ T
_{smax}all be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum temperature difference Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}Δ T is not more than after deducting 2 degrees Celsius
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly a moment in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T
_{o}", survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, T is obtained in actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises Cable Structure centrepoint 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 coordinate data, initial Cable Structure angledata, initial Cable Structure spatial data is in interior measured data, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Nondestructive Testing Data comprising support cable, and the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C
_{o}, monitored amount initial value vector C
_{o}the coding rule of coding rule and M monitored amount identical, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this element numerical value be 0, the initial value representing the change of this load is 0, if there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o},
Temperature variant physical and mechanical properties parameter, the initial Cable Structure bearing angular coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, support cable initial health data, Cable Structure centrepoint 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
_{o}with all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}"; Corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}" current initial Cable Structure steady temperature data vector T
^{t} _{o}"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}with monitored amount current initial value vector C
^{t} _{o}time, the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}just equal the initial mechanical Calculation Basis model A of Cable Structure
_{o}, monitored amount current initial value vector C
^{t} _{o}just equal monitored amount initial value vector C
_{o}; A
^{t} _{o}corresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T
^{t} _{o}", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}time, T
^{t} _{o}just equal T
_{o}; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}cable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}time, U
^{t} _{o}just equal U
_{o}; A
^{t} _{o}the initial health of evaluation object and A
_{o}the health status of evaluation object identical, also use evaluation object initial damage vector d
_{o}represent, A in cyclic process below
^{t} _{o}the initial health of evaluation object use evaluation object initial damage vector d all the time
_{o}represent; T
_{o}, U
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition; T
^{t} _{o}, U
^{t} _{o}and d
_{o}a
^{t} _{o}parameter, C
^{t} _{o}by A
^{t} _{o}mechanics Calculation result composition;
E. from entering the circulation being walked to m step by e here; In structure military service process, the current data of " Cable Structure steady temperature data " is constantly obtained according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T
^{t}", vector T
^{t}definition mode and vector T
_{o}definition mode identical; Current cable structure steady temperature data vector T is obtained in actual measurement
^{t}synchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas composition current cable structure actual measurement bearing angular coordinate vector U
^{t}, vectorial U
^{t}definition mode and vectorial U
_{o}definition mode identical;
F. according to current cable structure actual measurement bearing angular coordinate vector U
^{t}with current cable structure steady temperature data vector T
^{t}, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3
^{t} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}with current initial Cable Structure steady temperature data vector T
^{t} _{o};
F1. U is compared respectively
^{t}with U
^{t} _{o}, T
^{t}with T
^{t} _{o}if, U
^{t}equal U
^{t} _{o}and T
^{t}equal T
^{t} _{o}, then A
^{t} _{o}, U
^{t} _{o}, C
^{t} _{o}and T
^{t} _{o}remain unchanged, otherwise need to follow these steps to A
^{t} _{o}, U
^{t} _{o}, C
^{t} _{o}and T
^{t} _{o}upgrade;
F2. U is calculated
^{t}with U
_{o}difference, U
^{t}with U
_{o}difference be exactly the angular displacement of support of Cable Structure bearing about initial position, with angular displacement of support vector V represent angular displacement of support, V equals U
^{t}deduct U
_{o}, be onetoone relationship between the element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing; Calculate T
^{t}with T
_{o}difference, T
^{t}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{t}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{t}deduct T
_{o}, S represents the change of Cable Structure steady temperature data;
F3. first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}middle Cable Structure bearing applies angular displacement of support constraint and to A
_{o}in Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A
^{t} _{o}, upgrade A
^{t} _{o}while, U
^{t} _{o}all elements numerical value also uses U
^{t}all elements numerical value correspondence replaces, and namely have updated U
^{t} _{o}, T
^{t} _{o}all elements numerical value also uses T
^{t}all elements numerical value correspondence replace, namely have updated T
^{t} _{o}, so just obtain and correctly correspond to A
^{t} _{o}t
^{t} _{o}and U
^{t} _{o}; Upgrade C
^{t} _{o}method be: when renewal A
^{t} _{o}after, obtain A by Mechanics Calculation
^{t} _{o}in all monitored amounts, current concrete numerical value, these concrete numerical value composition C
^{t} _{o}; A
^{t} _{o}the initial health of support cable use evaluation object initial damage vector d all the time
_{o}represent;
G. at current initial mechanical Calculation Basis model A
^{t} _{o}basis on carry out several times Mechanics Calculation according to step g 1 to g4, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u};
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal
^{t} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}with current initial Cable Structure steady temperature data vector T
^{t} _{o}afterwards, Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D must then be upgraded
_{u};
G2. at the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object on the basis of original damage or load, increase unit damage or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just supposes that this support cable is at vectorial d
_{o}the basis that this support cable represented has a damage increases unit damage again, if this evaluation object is a load, just supposes that this load is at vectorial d
_{o}the basis that this load represented has a variable quantity increases load unit change again, use D
_{uk}record unit damage or the load unit change of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or load unit change, D
_{uk}evaluation object unit change vector D
_{u}an element, evaluation object unit change vector D
_{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical;
G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C
^{t} _{o}obtain a vector, then each element of this vector is calculated the unit damage or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;
G4. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
H. current cable structure steady temperature data vector T is obtained in actual measurement
^{t}while, actual measurement obtains at acquisition current cable structure steady temperature data vector T
^{t}the current measured value of all monitored amount of Cable Structure of synchronization in moment, form monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C
^{t} _{o}with monitored amount initial value vector C
_{o}definition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time;
I. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical;
J. according to monitored amount current value vector C with monitored amount current initial value vector C
^{t} _{o}, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;
formula 1
K. evaluation object current actual damage vector d is defined
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent actual damage degree or the real load variable quantity of corresponding evaluation object; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical;
L. the evaluation object utilizing formula 2 to express current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements;
K=1 in formula 2,2,3 ...., N, d
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; So according to evaluation object current actual damage vector d
^{a}impaired and the degree of injury of which support cable can be defined;
M. get back to e step, start the circulation next time being walked to m step by e.
Beneficial effect: structural healthy monitoring system is first by using sensor to carry out longterm online monitoring to structural response, after obtaining Monitoring Data, (or offline) analysis is online carried out to it and obtain structural health conditions data, due to the complicacy of structure, structural healthy monitoring system needs to use a large amount of sensor equipment to carry out monitoring structural health conditions, therefore its cost is usually quite high, and therefore cost problem is a subject matter of limit structural health monitoring technique application.On the other hand, the correct identification of the health status of core evaluation object (such as suspension cable) is the indispensable ingredient of the correct identification of structural health conditions, or even they are whole, and the impact of correct identification on the correct identification of the health status of Cable Structure of the change (such as by the change of the quality and quantity of the automobile of cablestayed bridge) of secondary evaluation object (load that such as structure is born) is very little, or even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is normally suitable, the quantity of secondary evaluation object is also usually greater than the quantity of core evaluation object, and the quantity of such evaluation object is usually many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, conventional method requires that the quantity of monitored amount (using sensor device measuring to obtain) must be more than or equal to the quantity of evaluation object, when the number ratio of the secondary evaluation object changed is larger (in fact often so), the quantity of the sensor equipment required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, unacceptablely even high.Inventor studies discovery, in the secondary evaluation object (normal load that such as structure is born, the normal load of structure refers to that the load that structure is being born is no more than the structure allowable load limited according to structural design book or structure completion book) change less time (be exactly that structure only bears normal load for load, whether the load that structure is born is normal load, can be observed by methods such as naked eyes and determine, if find that the load that structure is born is not normal load, so artificially remove, after removing improper load, structure just only bears normal load), the amplitude of variation (this instructions is called " secondary response ") of the structural response caused by them much smaller than core evaluation object change (such as support cable is impaired) caused by the amplitude of variation (this instructions is called " core response ") of structural response, secondary response and core respond total change (this instructions is called " global response ") that sum is structural response, obvious core response dominate in global response, based on this, find to choose when determining monitored amount quantity to be a bit larger tham core evaluation object quantity even if inventor studies, but much smaller than the numerical value (this method is exactly do like this) of evaluation object quantity, even if that is adopt the relatively few a lot of sensor equipment of quantity, still the state of health data of core evaluation object can accurately be obtained, meet the core demand of structural health conditions monitoring, therefore this method cost of structural healthy monitoring system of advising is more much lower than the cost of the structural healthy monitoring system required by conventional method apparently, that is this method can realize to the health status of the core evaluation object of Cable Structure with the much lower condition of cost assessment, can this benefit be used structural health monitoring technology is very important.
Embodiment
This method adopts a kind of algorithm, and this algorithm is for identifying the health status of core evaluation object.During concrete enforcement, the following step is the one in the various steps that can take.
The first step: the quantity first confirming the load that may change that Cable Structure is born.According to the feature of the load that Cable Structure is born, confirm wherein " load likely changed ", or all load is considered as " load likely changed ", if total JZW the load that may change, i.e. total JZW secondary evaluation object.If total M in cable system
_{1}root support cable, i.e. total M
_{1}individual core evaluation object.
If the quantity sum of the quantity of the support cable of Cable Structure and JZW " load likely changed " is N, i.e. total N number of evaluation object.To evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.
" the whole monitored spatial data of structure " by the specified point of K in structure and the volume coordinate of L assigned direction of each specified point describe, the change of structure space coordinate data is exactly the change of all volume coordinate components of K specified point.Each total individual volume coordinate measured value of M (M=K × L) or calculated value carry out characterisation of structures spatial coordinated information.Comprehensive abovementioned monitored amount, whole Cable Structure has M monitored amount, and the quantity that M must not be less than core evaluation object adds 4.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".To M monitored amount serial number, this numbering will be used for generating vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 with variable j ..., M.
Determine " the temperature survey calculating method of the Cable Structure of this method " according to the step specified in technical scheme and claims.
Second step: set up initial mechanical Calculation Basis model A
_{o}.
When Cable Structure is completed, or before setting up health monitoring systems, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T
_{o}represent, be called initial Cable Structure steady temperature data vector T
_{o}.T is obtained in actual measurement
_{o}while, use the direct survey calculation of conventional method to obtain the initial value of all monitored amount of Cable Structure, form monitored amount initial value vector C
_{o}.
Specifically the synchronization in moment of soandso Cable Structure steady temperature data vector such as (such as initial or current) can be being obtained according to following method in this method, soandso method survey calculation is used to obtain the data of the monitored amount of soandso measured amount (all monitored amount of such as Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) while, such as every 10 minutes survey records temperature, so simultaneously equally also every 10 minutes the monitored amount of soandso measured amount of survey record (all monitored amount of such as Cable Structure) data.Once determine the moment obtaining Cable Structure steady temperature data, so just be called the synchronization in the moment obtaining Cable Structure steady temperature data, the data of the monitored amount of soandso measured amount using soandso method survey calculation method to obtain with the data of the monitored amount of soandso measured amount (all monitored amount of such as Cable Structure) of the moment synchronization of acquisition Cable Structure steady temperature data.
Conventional method (consult reference materials or survey) is used to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses.
T is obtained in actual measurement
_{o}while, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Nondestructive Testing Data etc. that the Actual measurement data of Cable Structure comprise support cable can express the data of the health status of rope, the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data (comprise bearing about Descartes rectangular coordinate system X, Y, the volume coordinate of Z axis and angular coordinate and initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data), Cable Structure centrepoint load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structural point measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o}.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is the geometric properties according to these coordinate data determination Cable Structure.For cablestayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, socalled bridge type data that Here it is.Utilize the Nondestructive Testing Data etc. of support cable can express the data of the health status of support cable and Cable Structure load measurement data set up evaluation object initial damage vector d
_{o}(such as formula (1) Suo Shi), uses d
_{o}represent that Cable Structure is (with initial mechanical Calculation Basis model A
_{o}represent) the initial health of evaluation object.If there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0; If d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this element numerical value be 0, the initial value representing the change of this load is 0.The temperature variant physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, the Nondestructive Testing Data of support cable, Cable Structure to use, initial Cable Structure bearing angular coordinate vector U
_{o}with initial Cable Structure steady temperature data vector T
_{o}, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A
_{o}.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by
_{o}, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T
_{o}), based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A
_{o}the Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, Cable Structure angledata, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A
_{o}the health status evaluation object initial damage vector d of middle support cable
_{o}represent, the initial Cable Structure steady temperature data vector T of Cable Structure steady temperature data
_{o}represent.Due to based on A
_{o}the initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A
_{o}basis on, carry out Mechanics Calculation obtains, A
_{o}the evaluation of each monitored amount form monitored amount initial value vector C
_{o}.Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}", corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent, corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent, corresponding to A
_{o}cable Structure bearing angular data initial Cable Structure bearing angular coordinate vector U
_{o}represent, T
_{o}, U
_{o}and d
_{o}a
_{o}parameter, C
_{o}by A
_{o}mechanics Calculation result composition.
3rd step: first time sets up current initial mechanical Calculation Basis model A
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}" current initial Cable Structure steady temperature data vector T
^{t} _{o}", concrete grammar is: at initial time, and namely first time sets up current initial mechanical Calculation Basis model A
^{t} _{o}with monitored amount current initial value vector C
^{t} _{o}time, A
^{t} _{o}just equal A
_{o}, C
^{t} _{o}just equal C
_{o}, A
^{t} _{o}corresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T
^{t} _{o}", at initial time, (namely first time sets up A
^{t} _{o}time), T
^{t} _{o}just equal T
_{o}, vector T
^{t} _{o}definition mode and vector T
_{o}definition mode identical.Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}cable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}time, U
^{t} _{o}just equal U
_{o}.A
^{t} _{o}the health status of evaluation object and A
_{o}evaluation object health status (evaluation object initial damage vector d
_{o}represent) identical, A in cyclic process
^{t} _{o}the health status of evaluation object use evaluation object initial damage vector d all the time
_{o}represent.T
^{t} _{o}, U
^{t} _{o}and d
_{o}a
^{t} _{o}parameter, C
^{t} _{o}by A
^{t} _{o}mechanics Calculation result composition.
4th step: in Cable Structure military service process, the current data obtaining " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T
^{t}", vector T
^{t}definition mode and vector T
_{o}definition mode identical).T is obtained in actual measurement
^{t}while, actual measurement obtains the current measured value of all monitored amount of Cable Structure, composition " monitored amount current value vector C ".Current cable structure steady temperature data vector T is obtained in actual measurement
^{t}while, actual measurement obtains Cable Structure bearing angular coordinate current data, all data composition current cable structure actual measurement bearing angular coordinate vector U
^{t}.
5th step: according to current cable structure actual measurement bearing angular coordinate vector U
^{t}with current cable structure steady temperature data vector T
^{t}, upgrade current initial mechanical Calculation Basis model A where necessary
^{t} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}with current initial Cable Structure steady temperature data vector T
^{t} _{o}.Current cable structure actual measurement bearing angular coordinate vector U is obtained in the 4th step actual measurement
^{t}with current cable structure steady temperature data vector T
^{t}after, compare U respectively
^{t}and U
^{t} _{o}, T
^{t}and T
^{t} _{o}if, U
^{t}equal U
^{t} _{o}and T
^{t}equal T
^{t} _{o}, then do not need A
^{t} _{o}, U
^{t} _{o}and T
^{t} _{o}upgrade, otherwise need A
^{t} _{o}, U
^{t} _{o}and T
^{t} _{o}upgrade, update method is undertaken by the step specified in technical scheme and claims.
6th step: at current initial mechanical Calculation Basis model A
^{t} _{o}basis on carry out several times Mechanics Calculation, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u}.Concrete grammar is: Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal
^{t} _{o}while, Cable Structure unit damage monitored numerical quantity transformation matrices Δ C must be upgraded simultaneously; At the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity of all evaluation objects, N number of evaluation object is had just to have N calculating, calculating hypothesis each time only has an evaluation object on the basis of original damage or load, increase unit damage or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d
_{o}the basis that this support cable represented has a damage increases unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again, if this evaluation object is a load, just suppose that this load is at vectorial d
_{o}the basis that this load represented has a variable quantity increases again load unit change (if this load is distributed load, and this distributed load is line distributed load, load unit change can get 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. for unit change; If this load is distributed load, and this distributed load is EDS maps load, and load unit change can get 1MPa, 2MPa, 3MPa or 1kNm/m
^{2}, 2kNm/m
^{2}, 3kNm/m
^{2}deng be unit change; If this load is centrepoint load, and this centrepoint load is couple, and load unit change can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centrepoint load, and this centrepoint load is concentrated force, and load unit change can get 1kN, 2kN, 3kN etc. for unit change; If this load is volume load, load unit change can get 1kN/m
^{3}, 2kN/m
^{3}, 3kN/m
^{3}deng be unit change), use D
_{uk}record this unit damage or load unit change, wherein k represents the numbering of evaluation object unit damage occurring or load unit change occurs; Occur in calculating each time that the evaluation object of unit damage or load unit change is different from during other time calculates the evaluation object occurring unit damage or load unit change, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector C, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical; The monitored amount calculation current vector C calculated each time deducts monitored amount current initial value vector C
^{t} _{o}after calculate the unit damage supposed or load unit change numerical value divided by this time again, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector; The unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively by this N number of monitored amount unit change vector; Each row of unit damage monitored numerical quantity transformation matrices correspond to a monitored amount unit change vector, and every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object generation unit damage or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount.
7th step: set up linear relationship error vector e and vectorial g.Utilize (the monitored amount current initial value vector C of data above
^{t} _{o}, unit damage monitored numerical quantity transformation matrices Δ C), while the 6th step calculates each time, namely only have the increase unit damage of an evaluation object or load unit change D calculating each time in hypothesis evaluation object
_{uk}the evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current value all utilizing mechanics method (such as adopting finite element method) to calculate all monitored amounts in Cable Structure each time, while calculating the monitored amount calculation current vector C of composition one each time, calculate composition injury vector d each time, originally walk out of existing injury vector d only to use in this step, in all elements of injury vector d, only have the numerical value of an element to get D
_{uk}, the numerical value of other element gets 0, the coding rule of the element of injury vector d and vectorial d
_{o}the coding rule of element identical; By C, C
^{t} _{o}, Δ C, D
_{u}, d brings formula (1) into, obtain a linear relationship error vector e, calculate a linear relationship error vector e each time; N number of evaluation object is had just to have N calculating, just there is N number of linear relationship error vector e, obtaining a vector after being added by this N number of linear relationship error vector e, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N.Vector g equals final error vector e.
In formula (1), abs () is the function that takes absolute value, and takes absolute value to each element of the vector of trying to achieve in bracket.
8th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (such as containing spatial coordinate measuring system, signal conditioner etc.), Cable Structure bearing angular coordinate monitoring system (containing angle measuring sensor, signal conditioner etc.), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, computing machine and communication alert equipment.Each bearing angular coordinate, each temperature of each monitored amount, Cable Structure must arrive by monitored system monitoring, monitoring system by the Signal transmissions that monitors to signal (data) collector; Signal is delivered to computing machine through signal picker; The health monitoring software of the evaluation object running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring evaluation object health status and changing, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.
9th step: by current for monitored amount initial value vector C
^{t} _{o}, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D
_{u}parameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
Tenth step: establishment installation and operation this method system software on computers, the function (i.e. all work that can complete with computing machine in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method required by task and wants
11 step: according to monitored amount current value vector C with monitored amount current initial value vector C
^{t} _{o}, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D
_{u}and the linear approximate relationship (formula (2)) existed between evaluation object current nominal fatigue vector d (being made up of all Suo Dangqian nominal fatigue amounts), calculate the noninferior solution of evaluation object current nominal fatigue vector d according to multiobjective optimization algorithm, namely can determine the position of damaged cable and the solution of nominal fatigue degree thereof more exactly with reasonable error from all ropes.
The Objective Programming in multiobjective optimization algorithm (Goal Attainment Method) can be adopted to solve current injury vector d, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load intensity of variation of corresponding evaluation object; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical.
12 step: definition evaluation object current actual damage vector d
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent actual damage degree or the real load intensity of variation of corresponding evaluation object; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical.Utilize evaluation object current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements.
D
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion.
13 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.
14 step: under specified requirements, the computing machine automatic operation communication alert equipment in health monitoring systems is reported to the police to monitor staff, owner and (or) the personnel that specify.
15 step: get back to the 4th step, starts by the circulation of the 4th step to the 15 step.
Claims (1)
1. simplify angular displacement space coordinate monitoring damaged cable load recognition method, it is characterized in that described method comprises:
Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated support cable and load to be " evaluation object ", if the quantity sum of the quantity of evaluated support cable and load is N, namely the quantity of " evaluation object " is N; This method title " core evaluation object " specially refers to the evaluated support cable in " evaluation object ", and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; If total M in cable system
_{1}root support cable, the quantity of obvious core evaluation object is exactly M
_{1}; Determine to specify by the measured point of monitored volume coordinate, number to all specified points; Determined each measurement point by monitored volume coordinate component, to all measured volume coordinate component numberings; Abovementioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored spatial data of Cable Structure " is made up of abovementioned all measured volume coordinate components; For simplicity, in the method by " the monitored spatial data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M is greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount RealTime Monitoring in this method, the moment of survey record data is called the physical record data moment; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centrepoint load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centrepoint load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centrepoint load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated 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 comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact 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; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the onecomponent of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " is undertaken by step b1 to b3;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, asconstructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day
_{r}, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation
_{h}, get Δ T for convenience of describing
_{h}unit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, only measure the temperature of the surface point of support cable, measure all temperature be selected a little, the temperature recorded is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " to measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the onsite spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the nonsunny slope of reference plate is covered with insulation material, RealTime Monitoring is obtained the temperature of the sunny slope of reference plate,
B2: RealTime Monitoring obtains R Cable Structure surface temperature measured data of abovementioned R Cable Structure surface point, RealTime Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and RealTime Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T
_{emax}, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T
_{pmax}, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by RealTime Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T
_{smax}, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by RealTime Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and minimum temperature at the sea level elevation place that each is chosen, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T
_{tmax},
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, reference plate maximum temperature difference Δ T
_{pmax}with Cable Structure surface maximum temperature difference Δ T
_{smax}all be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum temperature difference Δ T that survey calculation obtains above
_{emax}be not more than with reference to temperature difference per day Δ T
_{r}, and reference plate maximum temperature difference Δ T
_{pmax}Δ T is not more than after deducting 2 degrees Celsius
_{emax}, and Cable Structure surface maximum temperature difference Δ T
_{smax}be not more than Δ T
_{pmax}, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T
_{tmax}be not more than 1 degree Celsius, this method utilizes abovementioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in abovementioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in abovementioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly a moment in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steadystate surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steadystate surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steadystate surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation
_{h}the numerical value obtained, gets Δ T for convenience of describing
_{h}unit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshineduration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshineduration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T
_{o}", survey or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, T is obtained in actual measurement
_{o}while, namely at the initial Cable Structure steady temperature data vector T of acquisition
_{o}the synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises Cable Structure centrepoint 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 coordinate data, initial Cable Structure angledata, initial Cable Structure spatial data is in interior measured data, initial Cable Structure bearing generalized coordinate data comprise initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Nondestructive Testing Data comprising support cable, and the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C
_{o}, monitored amount initial value vector C
_{o}the coding rule of coding rule and M monitored amount identical, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this element numerical value be 0, the initial value representing the change of this load is 0, if there is no the Nondestructive Testing Data of support cable and other are when can express the data of the health status of support cable, or can think structure original state be not damaged without relaxed state time, vectorial d
_{o}in each element numerical value relevant to support cable get 0, initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U
_{o},
Temperature variant physical and mechanical properties parameter, the initial Cable Structure bearing angular coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, support cable initial health data, Cable Structure centrepoint 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
_{o}with all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
_{o}, based on A
_{o}the Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T
_{o}"; Corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data "
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}" current initial Cable Structure steady temperature data vector T
^{t} _{o}"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}with monitored amount current initial value vector C
^{t} _{o}time, the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}just equal the initial mechanical Calculation Basis model A of Cable Structure
_{o}, monitored amount current initial value vector C
^{t} _{o}just equal monitored amount initial value vector C
_{o}; A
^{t} _{o}corresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T
^{t} _{o}", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}time, T
^{t} _{o}just equal T
_{o}; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}cable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time
^{t} _{o}time, U
^{t} _{o}just equal U
_{o}; A
^{t} _{o}the initial health of evaluation object and A
_{o}the health status of evaluation object identical, also use evaluation object initial damage vector d
_{o}represent, A in cyclic process below
^{t} _{o}the initial health of evaluation object use evaluation object initial damage vector d all the time
_{o}represent; T
_{o}, U
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition; T
^{t} _{o}, U
^{t} _{o}and d
_{o}a
^{t} _{o}parameter, C
^{t} _{o}by A
^{t} _{o}mechanics Calculation result composition;
E. from entering the circulation being walked to m step by e here; In structure military service process, the current data of " Cable Structure steady temperature data " is constantly obtained according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T
^{t}", vector T
^{t}definition mode and vector T
_{o}definition mode identical; Current cable structure steady temperature data vector T is obtained in actual measurement
^{t}synchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas composition current cable structure actual measurement bearing angular coordinate vector U
^{t}, vectorial U
^{t}definition mode and vectorial U
_{o}definition mode identical;
F. according to current cable structure actual measurement bearing angular coordinate vector U
^{t}with current cable structure steady temperature data vector T
^{t}, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3
^{t} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}with current initial Cable Structure steady temperature data vector T
^{t} _{o};
F1. U is compared respectively
^{t}with U
^{t} _{o}, T
^{t}with T
^{t} _{o}if, U
^{t}equal U
^{t} _{o}and T
^{t}equal T
^{t} _{o}, then A
^{t} _{o}, U
^{t} _{o}, C
^{t} _{o}and T
^{t} _{o}remain unchanged, otherwise need to follow these steps to A
^{t} _{o}, U
^{t} _{o}, C
^{t} _{o}and T
^{t} _{o}upgrade;
F2. U is calculated
^{t}with U
_{o}difference, U
^{t}with U
_{o}difference be exactly the angular displacement of support of Cable Structure bearing about initial position, with angular displacement of support vector V represent angular displacement of support, V equals U
^{t}deduct U
_{o}, be onetoone relationship between the element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing; Calculate T
^{t}with T
_{o}difference, T
^{t}with T
_{o}difference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T
^{t}with T
_{o}difference represent with steady temperature change vector S, S equals T
^{t}deduct T
_{o}, S represents the change of Cable Structure steady temperature data;
F3. first to A
_{o}in Cable Structure bearing apply angular displacement of support constraint, the numerical value of angular displacement of support constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A
_{o}in Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A
_{o}middle Cable Structure bearing applies angular displacement of support constraint and to A
_{o}in Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A
^{t} _{o}, upgrade A
^{t} _{o}while, U
^{t} _{o}all elements numerical value also uses U
^{t}all elements numerical value correspondence replaces, and namely have updated U
^{t} _{o}, T
^{t} _{o}all elements numerical value also uses T
^{t}all elements numerical value correspondence replace, namely have updated T
^{t} _{o}, so just obtain and correctly correspond to A
^{t} _{o}t
^{t} _{o}and U
^{t} _{o}; Upgrade C
^{t} _{o}method be: when renewal A
^{t} _{o}after, obtain A by Mechanics Calculation
^{t} _{o}in all monitored amounts, current concrete numerical value, these concrete numerical value composition C
^{t} _{o}; A
^{t} _{o}the initial health of support cable use evaluation object initial damage vector d all the time
_{o}represent;
G. at current initial mechanical Calculation Basis model A
^{t} _{o}basis on carry out several times Mechanics Calculation according to step g 1 to g4, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u};
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal
^{t} _{o}, current initial Cable Structure bearing angular coordinate vector U
^{t} _{o}, monitored amount current initial value vector C
^{t} _{o}with current initial Cable Structure steady temperature data vector T
^{t} _{o}afterwards, Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D must then be upgraded
_{u};
G2. at the current initial mechanical Calculation Basis model A of Cable Structure
^{t} _{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating hypothesis each time only has an evaluation object on the basis of original damage or load, increase unit damage or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just supposes that this support cable is at vectorial d
_{o}the basis that this support cable represented has a damage increases unit damage again, if this evaluation object is a load, just supposes that this load is at vectorial d
_{o}the basis that this load represented has a variable quantity increases load unit change again, use D
_{uk}record unit damage or the load unit change of this increase, wherein k represents the numbering of the evaluation object increasing unit damage or load unit change, D
_{uk}evaluation object unit change vector D
_{u}an element, evaluation object unit change vector D
_{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or load unit change, calculate the current calculated value all utilizing mechanics method to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical;
G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C
^{t} _{o}obtain a vector, then each element of this vector is calculated the unit damage or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;
G4. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
H. current cable structure steady temperature data vector T is obtained in actual measurement
^{t}while, actual measurement obtains at acquisition current cable structure steady temperature data vector T
^{t}the current measured value of all monitored amount of Cable Structure of synchronization in moment, form monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C
^{t} _{o}with monitored amount initial value vector C
_{o}definition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time;
I. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents nominal fatigue degree or the nominal load variable quantity of corresponding evaluation object; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical;
J. according to monitored amount current value vector C with monitored amount current initial value vector C
^{t} _{o}, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;
formula 1
K. evaluation object current actual damage vector d is defined
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent actual damage degree or the real load variable quantity of corresponding evaluation object; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical;
L. the evaluation object utilizing formula 2 to express current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements;
K=1 in formula 2,2,3 ...., N, d
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; So according to evaluation object current actual damage vector d
^{a}impaired and the degree of injury of which support cable can be defined;
M. get back to e step, start the circulation next time being walked to m step by e.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201510439869.0A CN105004561A (en)  20150723  20150723  Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201510439869.0A CN105004561A (en)  20150723  20150723  Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement 
Publications (1)
Publication Number  Publication Date 

CN105004561A true CN105004561A (en)  20151028 
Family
ID=54377317
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201510439869.0A Pending CN105004561A (en)  20150723  20150723  Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement 
Country Status (1)
Country  Link 

CN (1)  CN105004561A (en) 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

JP2006337144A (en) *  20050601  20061214  Kawasaki Heavy Ind Ltd  Fatigue life diagnostic method and diagnostic support device of bridge 
KR100806367B1 (en) *  20070119  20080227  한국유지관리 주식회사  Apparatus for exciting the bridge cable 
CN103616198A (en) *  20131209  20140305  东南大学  Identifying method for damaged cable and concentrated loads based on space coordinate monitoring in time of angular displacement 
CN103852315A (en) *  20140310  20140611  东南大学  Angular displacement space coordinates monitoring load identification method for problematic cable 
CN103884522A (en) *  20140310  20140625  东南大学  Method for recognizing damaged cables and loads by monitoring angular displacement space coordinates 
CN103913318A (en) *  20140310  20140709  东南大学  Method for progressively recognizing damaged cable and load based on angular displacement space coordinates monitoring 

2015
 20150723 CN CN201510439869.0A patent/CN105004561A/en active Pending
Patent Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

JP2006337144A (en) *  20050601  20061214  Kawasaki Heavy Ind Ltd  Fatigue life diagnostic method and diagnostic support device of bridge 
KR100806367B1 (en) *  20070119  20080227  한국유지관리 주식회사  Apparatus for exciting the bridge cable 
CN103616198A (en) *  20131209  20140305  东南大学  Identifying method for damaged cable and concentrated loads based on space coordinate monitoring in time of angular displacement 
CN103852315A (en) *  20140310  20140611  东南大学  Angular displacement space coordinates monitoring load identification method for problematic cable 
CN103884522A (en) *  20140310  20140625  东南大学  Method for recognizing damaged cables and loads by monitoring angular displacement space coordinates 
CN103913318A (en) *  20140310  20140709  东南大学  Method for progressively recognizing damaged cable and load based on angular displacement space coordinates monitoring 
Similar Documents
Publication  Publication Date  Title 

CN105136495A (en)  Simplified angular displacement strain monitoring damaged cable load identification method  
CN104990747A (en)  Simplified method of recognizing problematic cable load through generalized displacement space coordinate monitoring  
CN105004561A (en)  Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined annular displacement  
CN105004549A (en)  Method for recognizing load of damaged cable based on streamlined spacecoordinate monitoring process of generalized displacement  
CN105004376A (en)  Recognition method for damaged cable load support angular displacement based on streamlined angle monitoring process  
CN104990750A (en)  Simplified method of recognizing damaged cable load through linear displacement hybrid monitoring  
CN105067338A (en)  Method for recognizing load of damaged cable based on streamlined angle monitoring process  
CN104990582A (en)  Simplified damaged cable load identification method through angle monitoring under linear displacement  
CN105067355A (en)  Method for recognizing load of damaged cable based on streamlined strain monitoring process  
CN105067334A (en)  Method for recognizing load generalized displacement of damaged cable based on streamlined strain monitoring process  
CN105115765A (en)  Identification method for damaged cable, load and support angular displacement through simplified strain monitoring  
CN105067315A (en)  Method for recognizing load of damaged cable based on hybrid monitoring process of streamlined angular displacement  
CN105004543A (en)  Method for recognizing load linear displacement of damaged cable based on streamlined space coordinate monitoring process  
CN105115766A (en)  Identification method for damaged cable and load through simplified monitoring of linear displacement and strain  
CN104990739A (en)  Simplified damaged cable load identification method through angle monitoring under generalized displacement  
CN104990581A (en)  Simplified damaged cable load support angular displacement identification method through space coordinate monitoring  
CN105115760A (en)  Identification method for load of damaged cable through simplified space coordinate monitoring  
CN104977178A (en)  Simplified hybrid monitoring damaged cable load bearing angular displacement recognition method  
CN105115759A (en)  Identification method for damaged cable load and generalized displacement through simplified angle monitoring  
CN105067329A (en)  Method for recognizing load of damaged cable based on streamlined hybrid monitoring process of generalized displacement  
CN105067336A (en)  Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process  
CN105115763A (en)  Identification method for damaged cable and load through simplified monitoring of angular displacement and cable force  
CN105004548A (en)  Method for recognizing load of damaged cable based on angle monitoring process of streamlined linear displacement  
CN105115750A (en)  Identification method for damaged cable load and linear displacement through simplified angle monitoring  
CN105004560A (en)  Method for recognizing load of damaged cable based on cable force monitoring process of streamlined linear displacement 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
RJ01  Rejection of invention patent application after publication 
Application publication date: 20151028 

RJ01  Rejection of invention patent application after publication 