CN103852286A - Strain monitoring error cable load generalized displacement identification method - Google Patents

Strain monitoring error cable load generalized displacement identification method Download PDF

Info

Publication number
CN103852286A
CN103852286A CN201410085856.3A CN201410085856A CN103852286A CN 103852286 A CN103852286 A CN 103852286A CN 201410085856 A CN201410085856 A CN 201410085856A CN 103852286 A CN103852286 A CN 103852286A
Authority
CN
China
Prior art keywords
cable structure
temperature
data
vector
cable
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
Application number
CN201410085856.3A
Other languages
Chinese (zh)
Inventor
韩玉林
韩佳邑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201410085856.3A priority Critical patent/CN103852286A/en
Publication of CN103852286A publication Critical patent/CN103852286A/en
Pending legal-status Critical Current

Links

Landscapes

  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

The invention relates to a strain monitoring error cable load generalized displacement identification method. Based on strain monitoring, whether a mechanical calculation reference model of a cable structure is required to be updated or not is judged by monitoring the cable structure temperature and the environment temperature, the mechanical calculation reference model, including the cable structure temperature and the environment temperature, of the cable structure is obtained, and a unit damaged monitored quantity value change matrix is obtained through calculation on the basis of the mechanical calculation reference model. A noninferior solution of a current nominal damaged vector of an evaluated object is calculated according to the approximate linear relation among a monitored quantity current value vector, a monitored quantity current initial value vector, the unit damaged monitored quantity value change matrix and a to-be-obtained evaluated object current nominal damaged vector, and therefore support generalized displacement, the load fluctuation amount and an error cable can be simultaneously identified when the temperature is changed.

Description

Strain monitoring problem rope load generalized displacement recognition methods
Technical field
Cable-stayed bridge, suspension bridge, the structures such as truss-frame structure have a common ground, be exactly that they have many parts that bear tensile load, as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or the rod member that only bears tensile load are support unit, for simplicity, this method is " Cable Structure " by such structure representation, and by all ropeway carrying-ropes of Cable Structure, carrying cable, and all rod members (being called again two power rod members) that only bear axial tension or axial compression load, unified be called " cable system " for simplicity, in this method, censure ropeway carrying-rope with " support cable " this noun, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word in the back, truss-frame structure reality is just referred to two power rod members.Impaired and the lax pair Cable Structure of support cable is safely a significant threat, and this method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope.In structure military service process, the correct identification of the health status to support cable or cable system is related to the safety of whole Cable Structure.In the time that environment temperature changes, the temperature of Cable Structure generally also can be along with changing, in the time that Cable Structure temperature changes, may there is generalized displacement in Cable Structure bearing, the load that Cable Structure is born also may change, even if in fact the temperature of Cable Structure does not change, the load that Cable Structure is born also may change separately, the health status of Cable Structure also may change simultaneously, at this complex condition, this method is identified generalized displacement of support based on strain monitoring (this method is called monitored strain " monitored amount "), the variable quantity of the load that problem rope and Cable Structure are born, belong to engineering structure health monitoring field.
Background technology
Reject load change, Cable Structure generalized displacement of support and structure temperature and change the impact on Cable Structure health status recognition result, thereby the variation of the health status of recognition structure is exactly current problem in the urgent need to address; Rejecting load change, the variation of Cable Structure health status and structure temperature and change the impact on Cable Structure generalized displacement of support recognition result, thereby identify exactly Cable Structure generalized displacement of support, is also current problem in the urgent need to address; Same, the variation of rejecting structure temperature, Cable Structure generalized displacement of support and structural health conditions change the impact of the recognition result of the variable quantity of the load that structure is born, significant equally to structural safety.Based on structural health monitoring technology, this method discloses a kind of effective ways that solve these three problems.
Summary of the invention
Technical matters: this method discloses a kind of method, three kinds of functions that existing method can not possess are realized, be respectively, one, reject generalized displacement of support, load change and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the health status of support cable; Two, this method, in identifying problem rope, can also identify the variation of load simultaneously, and this method can be rejected the impact that generalized displacement of support, structure temperature variation and support cable health status change, and realizes the correct identification of load change degree; Three, this method can also be rejected load change, Cable Structure health status changes and structure temperature changes the impact on Cable Structure generalized displacement of support recognition result, thereby identifies exactly Cable Structure generalized displacement of support.
In Cable Structure military service process, Suo Changdu under support cable free state (now Suo Zhangli also claims that Suo Li is 0) (is called drift, this method specially refers to that support cable two supports the drift of that section of rope between end points) can change, one of object of this method will identify exactly drift the support cable changing has occurred, and identify the change amount of their drift, the Suo Li that this change amount is this rope adjusts provides direct basis, for convenient, the support cable that this method changes drift is referred to as slack line.
Technical scheme: in the method, censure the coordinate of bearing about the X, Y, Z axis of Descartes's rectangular coordinate system with " bearing volume coordinate ", also can be said to is the volume coordinate of bearing about X, Y, Z axis, bearing is called the volume coordinate component of bearing about this axle about the concrete numerical value of the volume coordinate of some axles, and in this method, also a volume coordinate component with bearing is expressed the concrete numerical value of bearing about the volume coordinate of some axles; Censure the angular coordinate of bearing about X, Y, Z axis with " bearing angular coordinate ", bearing is called the angular coordinate component of bearing about this axle about the concrete numerical value of the angular coordinate of some axles, and in this method, also an angular coordinate component with bearing is expressed the concrete numerical value of bearing about the angular coordinate of some axles; All by " bearing generalized coordinate " denotion bearing angular coordinate and bearing volume coordinate, in this method, also a generalized coordinate component with bearing is expressed the concrete numerical value of bearing about volume coordinate or the angular coordinate of an axle; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, also can say that the change of bearing volume coordinate is called support wire displacement, and in this method, also a translational component with bearing is expressed the concrete numerical value of bearing about the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, and in this method, also an angular displacement component with bearing is expressed the concrete numerical value of bearing about the angular displacement of some axles; Generalized displacement of support denotion support wire displacement and angular displacement of support are all, and in this method, also a generalized displacement component with bearing is expressed bearing about the displacement of the lines of some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is support wire displacement or the translational displacement component at gravity direction.
The external force that object, structure are born can be described as load, and load comprises face load and volume load.Face load claims again surface load, is the load that acts on body surface, comprises two kinds of centre-point load and distributed loads.Volume load is that continuous distribution is in the load of interior of articles each point, as the deadweight of object and inertial force.
Centre-point load is divided into two kinds of concentrated force and concentrated couples, in coordinate system, for example, in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component are called to a load, the now variation of load is embodied as the variation of a concentrated force component or a concentrated couple component.
Distributed load is divided into line distributed load and face distributed load, the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, the size of distributed load is expressed by distribution intensity, distribution for intensity distribution characteristics (for example uniform, the distribution characteristicss such as sine function) and amplitude is expressed, and (for example two distributed loads are all uniform, but its amplitude difference, can well-distributed pressure be the concept that example illustrates amplitude: same structure is born two different well-distributed 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 is talked about the variation 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 coordinate system, a distributed load can be resolved into several components, if the amplitude of the distribution intensity separately of several components of this distributed load changes, and the ratio changing is all not identical, so in the method the component of these several distributed loads is regarded as to the independently distributed load of same quantity, now load just represents the component of a distributed load, also component identical the amplitude changing ratio of the intensity that wherein distributes can be synthesized to 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 the deadweight of object and inertial force, the description of volume load at least comprises the zone of action of volume load and the size of volume load, the size of volume load is expressed by distribution intensity, distribution for intensity distribution characteristics (for example uniform, the distribution characteristicss such as linear function) and amplitude is expressed, and (for example two individual stow lotuses are all uniform, but its amplitude difference, can conduct oneself with dignity the concept of amplitude is described for example: the material difference of two parts of same structure, therefore density difference, so although the suffered volume load of these two parts is all uniform, but the amplitude of the suffered volume load of part may be 10kN/m 3, the amplitude of the suffered volume load of another part is 50kN/m 3).If load is actually volume load, actual treatment is the change of the amplitude of volume load distribution intensity in the method, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, while now mentioning in the method the change of load, in fact refer to the change of the amplitude of the distribution intensity of volume load, now, the load changing refers to the volume load that the amplitude of those distribution intensities changes.In coordinate system, one individual stow lotus can be resolved into several components (for example, in Descartes's rectangular coordinate system, volume load can resolve into the component about three axles of coordinate system, that is to say, in Descartes's rectangular coordinate system, volume load can resolve into three components), if the amplitude of the distribution intensity separately of several components of this volume load changes, and the ratio changing is all not identical, so in the method the component of this several body stow lotus is regarded as to the independently load of same quantity, also the volume sharing part of the load identical the amplitude changing ratio of the intensity that wherein distributes can be synthesized to an individual stow lotus or be called a load.
In the time that load is embodied as centre-point load, in the method, " load unit variation " in fact refers to " unit change of centre-point load ", similarly, " load change " specifically refers to " the big or small variation of centre-point load ", " load change amount " specifically refers to " the big or small variable quantity of centre-point load ", " load change degree " specifically refers to " the big or small intensity of variation of centre-point load ", " the actual change amount of load " refers to " the big or small actual change amount of centre-point load ", " load changing " refers to " centre-point load that size changes ", briefly, now " so-and-so load so-and-so variation " refers to " so-and-so centre-point load big or small so-and-so variation ".
In the time that load is embodied as distributed load, in the method, " load unit variation " 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, similarly, " load change " specifically refers to " variation 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 changing " refers to " distributed load that the amplitude of distribution intensity changes ", briefly, now " so-and-so load so-and-so variation " refers to " amplitude of the distribution intensity of so-and-so distributed load so-and-so variation ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.
In the time that load is embodied as volume load, in the method, " load unit variation " in fact refers to " unit change of the amplitude of the distribution intensity of volume load ", similarly, " load change " refers to " variation 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 changing " refers to " the volume load that the amplitude of distribution intensity changes ", briefly, " so-and-so load so-and-so variation " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so variation ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.
This method specifically comprises:
A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated generalized displacement of support component, support cable and load, if the quantity sum of the quantity of evaluated generalized displacement of support component, the quantity of support cable and load is N, the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be used for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; Determine the point being monitored of specifying, point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified point numberings; Determine monitored should the changing direction of point being monitored, and number to the monitored strain of all appointments, " monitored strain numbering " will be used for generating vector sum matrix in subsequent step, and " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; In this method, must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring, the moment of survey record data is called the physical record data moment; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load claims again surface load, is the load that acts on body surface, comprises two kinds of centre-point load and distributed loads; Volume load is that continuous distribution is in the load of interior of articles each point, including the deadweight and inertial force of object; Centre-point load is divided into two kinds of concentrated force and concentrated couples, comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now variation of load is embodied as the variation of a concentrated force component or a concentrated couple component; Distributed load is divided into line distributed load and face distributed 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 of distributed load is expressed by distribution intensity, and distribution intensity is expressed by distribution characteristics and amplitude; If load is actually distributed load, when this method is talked about the variation 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; Comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, a distributed load can be resolved into three components, if the amplitude of the distribution intensity separately of three components of this distributed load changes, and the ratio changing is all not identical, so in the method three components of this distributed load being counted or added up is three distributed loads, and now load just represents the one-component 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, the size of volume load is expressed by distribution intensity, distribution intensity is expressed by distribution characteristics and amplitude; If load is actually volume load, actual treatment is the change of the amplitude of volume load distribution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, while now mentioning in the method the change of load, in fact refer to the change of the amplitude of the distribution intensity of volume load, now, the load changing refers to the volume load that the amplitude of those distribution intensities changes; Comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, one individual stow lotus can be resolved into three components, if the amplitude of the distribution intensity separately of three components of this volume load changes, and the ratio changing is all not identical, and so in the method three components of this volume load 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 Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, do not represent necessarily can see the same day sun, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day r, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of narration hunit be DEG C/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ", from the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at least choose two points at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, measure Cable Structure along each and be no less than three points along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual measurement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day, at the flat board of a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained to the temperature of the sunny slope of reference plate,
B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains the temperature record of the Cable Structure place environment that meets the requirement of meteorology measurement temperature simultaneously, obtain being carved at sunrise the same day temperature measured data sequence of the Cable Structure place environment between 30 minutes after sunrise moment next day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains Cable Structure place environment by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T emax, temperature measured data sequence by Cable Structure place environment obtains the temperature of Cable Structure place environment about the rate of change of time by conventional mathematical computations, and this rate of change is also along with the time changes, obtain being carved at sunrise the same day measured data sequence of the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data that was carved at sunrise the temperature of the sunny slope of the reference plate between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the measured data sequence of temperature of the sunny slope of reference plate, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of the sunny slope of reference plate by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T pmax, obtain being carved at sunrise the same day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between 30 minutes after sunrise moment next day by Real-Time Monitoring, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by being carved at sunrise the Cable Structure surface temperature measured data between 30 minutes after the sunrise moment next day same day of a Cable Structure surface point, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of each Cable Structure surface point by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T smax, obtain the temperature of each Cable Structure surface point about the rate of change of time by each Cable Structure surface temperature measured data sequence by conventional mathematical computations, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the sea level elevation place that chooses at each and amount to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", choose H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax,
B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, reference plate maximum temperature difference Δ T pmaxwith Cable Structure surface maximum temperature difference Δ T smaxall be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates above emaxbe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T pmaxafter deducting 2 degrees Celsius, be not more than Δ T emax, and Cable Structure surface maximum temperature difference Δ T smaxbe not more than Δ T pmax, only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition, Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 5 condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day, Section 6 condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius, this method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to carry out the relevant health monitoring analysis of Cable Structure, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 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 hthe numerical value obtaining, gets Δ T for convenience of narration hunit be DEG C/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure,
C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " the Cable Structure steady temperature data that calculate under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T o", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, obtain initial Cable Structure steady temperature data vector T in actual measurement osynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F o, according to data including Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing osuo Li under condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the weight vector of the initial drift vector of support cable, the initial free unit length of initial free cross-sectional area vector sum, the coding rule of the element of the weight vector of the initial drift vector of support cable, the initial free unit length of initial free cross-sectional area vector sum and initial rope force vector F othe coding rule of element identical, obtain T in actual measurement otime, namely obtaining initial Cable Structure steady temperature data vector T othe synchronization in moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, in obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C o, monitored amount initial value vector C othe coding rule of coding rule and M monitored amount identical, utilize support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement data to set up evaluation object initial damage vector d o, vectorial d orepresent with initial mechanical calculating benchmark model A othe initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d oelement number equal N, d oelement and evaluation object be one-to-one relationship, vectorial d othe coding rule of element identical with the coding rule of evaluation object, if d oevaluation object corresponding to some elements be support cable, so a d in cable system othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, represents this support cable, , if d oevaluation object corresponding to some elements be some generalized displacement components of some bearings, d so othe numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d oevaluation object corresponding to some elements be some load, in this method, get d othis element numerical value be 0, the initial value that represents the variation of this load is 0, if there is no the initial generalized displacement measurement data of Cable Structure bearing or can think that the initial generalized displacement of Cable Structure bearing is at 0 o'clock, vectorial d oin the each element numerical value relevant to Cable Structure generalized displacement of support get 0, if while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d oin the each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data refer to the bearing generalized coordinate data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing refers to setting up initial mechanical calculating benchmark model A otime, the generalized displacement that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point,
The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, initial Cable Structure steady temperature data vector T owith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " o, based on A othe Cable Structure computational data calculating must approach its measured data very much, and difference therebetween must not be greater than 5%; Corresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A oevaluation object initial damage vector d for evaluation object health status orepresent; Corresponding to A omonitored amount initial value vector C for the initial value of all monitored amounts orepresent; Set up for the first time the current initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure t owith the current initial value vector of monitored amount C t otime, the current initial mechanical calculating benchmark model A of Cable Structure t ojust equal the initial mechanical calculating benchmark model A of Cable Structure o, the current initial value vector of monitored amount C t ojust equal monitored amount initial value vector C o; A t ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure t otime, T t ojust equal T o; A t oinitial health and the A of evaluation object othe health status of evaluation object identical, also use evaluation object initial damage vector d orepresent A in cyclic process below t othe initial health of evaluation object use all the time evaluation object initial damage vector d orepresent; T oand d oa oparameter, by A oinitial value and the C of all monitored amounts of obtaining of Mechanics Calculation result othe initial value of all monitored amounts that represent is identical, therefore also can say C oby A omechanics Calculation result composition; T t oand d oa t oparameter, C t oby A t omechanics Calculation result composition;
E. from entering the circulation that is walked n step by e here; In structure military service process, constantly according to " the temperature survey calculating method of the Cable Structure of this method " the constantly current data of Actual measurement acquisition " Cable Structure steady temperature data ", the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical; Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, actual measurement obtains all M in Cable Structure 1the rope force data of root support cable, all these rope force data composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, Actual measurement obtains all M 1the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports the coding rule and initial rope force vector F of the element of end points horizontal range vector othe coding rule of element identical;
F. according to current cable structure steady temperature data vector T t, upgrade current initial mechanical calculating benchmark model A according to step f1 to f3 t o, the current initial value of monitored amount vector C t owith current initial Cable Structure steady temperature data vector T t o;
F1. compare T twith T t oif, T tequal T t o, A t o, C t oand T t oremain unchanged; Otherwise need to follow these steps to A t oand T t oupgrade;
F2. calculate T twith T opoor, T twith T odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T opoor represent with steady temperature change vector S, S equals T tdeduct T o, S represents the variation of Cable Structure steady temperature data;
F3. to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A oin the temperature variation that applies of Cable Structure after obtain the current initial mechanical calculating benchmark model A that upgrades t o, upgrade A t otime, T t oall elements numerical value is also used T tcorresponding replacement of all elements numerical value, upgraded T t o, so just obtained correctly corresponding to A t ot t o; Upgrade C t omethod be: when upgrade A t oafter, obtain A by Mechanics Calculation t oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C t o; A t othe initial health of support cable use all the time evaluation object initial damage vector d orepresent;
G. at current initial mechanical calculating benchmark model A t obasis on carry out several times Mechanics Calculation according to step g 1 to g4, obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D by calculating u;
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A t o, the current initial value of monitored amount vector C t owith current initial Cable Structure steady temperature data vector T t oafterwards, must then upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D u;
g2. at the current initial mechanical calculating benchmark model A of Cable Structure t obasis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or generalized displacement or load, to increase unit damage or unit generalized displacement or load unit variation again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d oon the basis of the existing damage of this support cable representing, increase again unit damage, if this evaluation object is the generalized displacement component of a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, if this evaluation object is a load, just suppose that this load is at vectorial d oon the basis of the existing variable quantity of this load representing, increase again load unit and change, use D ukunit damage or unit generalized displacement or the load unit variation of recording this increase, wherein k represents the numbering of the evaluation object that increases unit damage or unit generalized displacement or load unit variation, D ukevaluation object unit change vector D uan element, evaluation object unit change vector D ucoding rule and the vectorial d of element othe coding rule of element identical; The evaluation object that increases unit damage or unit generalized displacement or load unit variation in calculating is each time different from the evaluation object that increases unit damage or unit generalized displacement or load unit variation in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, a monitored amount calculation current vector of current calculated value composition of all monitored amounts that calculate each time, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C oelement coding rule identical;
G3. the monitored amount calculation current vector calculating each time deducts the current initial value vector of monitored amount C t oobtain a vector, again each element of this vector is calculated to unit damage or unit generalized displacement or the load unit variation numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N evaluation object just to have N monitored amount unit change vector;
G4. by this N monitored amount unit change vector according to the coding rule of N evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is listed as successively; Each of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is listed as corresponding to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is the different unit change amplitude in the time that different evaluation objects increase unit damage or unit generalized displacement or load unit variation corresponding to same monitored amount; Coding rule and the vectorial d of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
H. obtain current cable structure steady temperature data vector T in actual measurement ttime, actual measurement obtains obtaining current cable structure steady temperature data vector T tthe current measured value of all monitored amounts of Cable Structure of synchronization in moment, form monitored amount current value vector C; The current initial value vector of monitored amount current value vector C and monitored amount C t owith monitored amount initial value vector C odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at concrete numerical value in the same time not;
I. define the vectorial d of the current name damage of evaluation object, the element number of the vectorial d of the current name damage of evaluation object equals the quantity of evaluation object, between the element of the vectorial d of the current name damage of evaluation object and evaluation object, be one-to-one relationship, the element numerical value of the vectorial d of the current name damage of evaluation object represents nominal degree of injury or nominal generalized displacement or the nominal load variable quantity of corresponding evaluation object; Coding rule and the vectorial d of the element of vector d othe coding rule of element identical;
J. the monitored amount current value vector of foundation C is with the current initial value vector of monitored amount C t o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and the vectorial d of the current name damage of evaluation object to be asked, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate the vectorial d of the current name damage of evaluation object;
C = C o t + ΔC · d Formula 1
K. define the current actual damage vector of evaluation object d a, the current actual damage vector of evaluation object d aelement number equal the quantity of evaluation object, the current actual damage vector of evaluation object d aelement and evaluation object between be one-to-one relationship, the current actual damage of evaluation object vector d aelement numerical value represent actual damage degree or actual generalized displacement or the real load variable quantity of corresponding evaluation object; Vector d acoding rule and the vectorial d of element othe coding rule of element identical;
L. the current actual damage vector of the evaluation object d that utilizes formula 2 to express ak element d a kwith evaluation object initial damage vector d ok element d okk the element d with the vectorial d of the current name damage of evaluation object kbetween relation, calculate the current actual damage of evaluation object vector d aall elements;
Figure BDA0000474929560000152
formula 2
K=1 in formula 2,2,3 ...., N, d a krepresent the current actual health status of k evaluation object, d a kbe to represent that k evaluation object is without health problem at 0 o'clock, d a knumerical value is not to represent that k evaluation object is the evaluation object of unsoundness problem at 0 o'clock, if this evaluation object is support cable, so a d in cable system a krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d a knumerical response the degree of lax or damage of this support cable; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, the current actual damage vector of evaluation object d ain with slack line corresponding to element numerical expression be and the current actual equivalent damage degree of slack line relax level mechanics equivalence; If this evaluation object is generalized displacement component, so a d of a bearing a krepresent its current actual generalized displacement numerical value; If this evaluation object is load, so a d a krepresent the actual change amount of this load;
M. utilize at current cable structure steady temperature data vector T tunder condition, walk the slack line that identifies and with the current actual damage vector of evaluation object d at l athese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, utilize obtain in e step at current cable structure steady temperature data vector T tcurrent cable force vector F under condition and current support cable two support end points horizontal range vector, utilize obtain in c step at initial Cable Structure steady temperature data vector T othe weight vector of the initial drift vector of the support cable under condition, the initial free unit length of initial free cross-sectional area vector sum, initial rope force vector F o, utilize current cable structure steady temperature data vector T tthe current steady temperature data of support cable that represent, utilize obtain in c step at initial Cable Structure steady temperature data vector T othe support cable initial steady state temperature data representing, utilize the temperature variant physical and mechanical properties parameter of the various materials that use in the Cable Structure of c step acquisition, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanics equivalent condition is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; While meeting above-mentioned two mechanics equivalent conditions, the mechanics function of two such support cables in Cable Structure is exactly identical, if replaced after damaged cable with equivalent slack line, any variation can not occur Cable Structure, and vice versa; The relax level of trying to achieve those and be judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of those support cables that need adjust Suo Li; Lax identification and the damage identification of support cable are so just realized; When calculating, institute's demand power is provided by current cable force vector F corresponding element; This method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope, so this method is according to the current actual damage vector of evaluation object d acan either identify problem rope, can define which bearing generalized displacement and numerical value thereof have occurred, also can define which load the numerical value that changes and change has occurred; So far this method has realized and has rejected problem rope identification impact, Cable Structure that generalized displacement of support, load change and structure temperature change, realize simultaneously and rejected generalized displacement of support, structure temperature variation and identification support cable health status variable effect, load change amount, also realized and rejected load change, structure temperature variation and identification variable effect, generalized displacement of support of support cable health status;
N. get back to e step, start to be walked by e the circulation next time of n step.
Beneficial effect: this method has realized three kinds of functions that existing method can not possess, respectively: one, reject Cable Structure generalized displacement of support, load change and structure temperature and change the impact on Cable Structure health status recognition result, thereby identify exactly the structure health monitoring method of problem rope; Two, reject the impact that Cable Structure generalized displacement of support, structure temperature variation and support cable health status change, realize the correct identification of load change degree; Three, reject the impact that load change, structure temperature variation and support cable health status change, realize the correct identification of Cable Structure generalized displacement of support.
Embodiment
This method adopts a kind of algorithm, and this algorithm is for identifying the variation of generalized displacement of support, problem rope and load.When concrete enforcement, the following step is the one in the various steps that can take.
The first step: the quantity of first confirming the load that may change that Cable Structure bears.The feature of the load of bearing according to Cable Structure, confirm wherein " load likely changing ", or all load is considered as " load likely changing ", if total JZW the load that may change, the variable quantity that this method is expressed " load likely changing " by the intensity of variation of identification this JZW " load likely changing ".
If the quantity sum of the quantity of the quantity of the generalized displacement of support component of Cable Structure, the support cable of Cable Structure and JZW " load likely changing " is N.For sake of convenience, it is " evaluation object " that this method unitedly calls evaluated generalized displacement of support component, support cable and " load likely changing ", total N evaluation object.Give evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.
" the whole monitored strain data of structure " can be described by the strain of a L assigned direction specified point and each specified point of K in structure, and the variation of structural strain data is exactly the variation of all strains of K specified point.Each total M(M=K × L) individual strain measurement value or calculated value characterize structural strain information.K and M generally must not be less than N.
Comprehensive above-mentioned monitored amount, total M the monitored amount of whole Cable Structure, M must not be less than the quantity N of evaluation object.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".Give M monitored amount serial number, this numbering will be used for generating vector sum matrix in subsequent step.This method is with representing this numbering with variable j, j=1, and 2,3 ..., M.
The step specifying by technical scheme is determined " the temperature survey calculating method of the Cable Structure of this method ".
Second step: set up initial mechanical calculating benchmark model A o.
In Cable Structure completion, or setting up before health monitoring systems, calculating " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " measurement (can measure by conventional thermometry, for example use thermal resistance to measure), " Cable Structure steady temperature data " now use vector T orepresent, be called initial Cable Structure steady temperature data vector T o.Obtain T in actual measurement otime, namely at the synchronization in moment that obtains initial Cable Structure steady temperature data vector, use conventional method directly to measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial value vector C o.
Can be specifically in this method obtaining the synchronization in moment of so-and-so Cable Structure steady temperature data vector such as (such as initial or current) according to following method, use so-and-so method measurement to calculate the data of the monitored amount of so-and-so measured amount (all monitored amount of for example Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) time, for example, every temperature of 10 minutes survey records, so simultaneously equally also every 10 minutes the monitored amount of so-and-so measured amount of survey record (all monitored amount of for example Cable Structure) data.Once determine the moment that obtains Cable Structure steady temperature data, for example, be just called and obtaining the synchronization in moment of Cable Structure steady temperature data with the data of the monitored amount of so-and-so measured amount (all monitored amount of Cable Structure) of moment synchronization that obtain Cable Structure steady temperature data so, use so-and-so method to measure the data of the monitored amount of so-and-so measured amount that computing method obtain.
Use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the various materials that Cable Structure uses.
Obtain initial Cable Structure steady temperature data vector T in actual measurement osynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F o; According to Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing osuo Li under condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, forms successively the initial drift vector l of support cable o, initial free cross-sectional area vector A oweight vector ω with initial free unit length o, the initial drift vector l of support cable o, initial free cross-sectional area vector A oweight vector ω with initial free unit length othe coding rule and initial rope force vector F of element othe coding rule of element identical.
The step specifying by technical scheme, obtains initial Cable Structure steady temperature data vector T at Actual measurement otime, namely, obtaining the synchronization in moment of Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.Utilize the Non-destructive Testing Data of the design drawing, as-constructed drawing of Cable Structure and the measured data of initial Cable Structure, support cable, temperature variant physical and mechanical properties parameter and the initial Cable Structure steady temperature data vector T of various materials that Cable Structure is used o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical calculating benchmark model A o.T oand d oa oparameter, C oby A omechanics Calculation result composition.
The 3rd step: set up for the first time current initial mechanical calculating benchmark model A t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o", concrete grammar is: at initial time, set up for the first time current initial mechanical calculating benchmark model A t owith the current initial value vector of monitored amount C t otime, A t ojust equal A o, C t ojust equal C o, A t ocorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T t o", (namely set up for the first time A at initial time t otime), T t ojust equal T o, vector T t odefinition mode and vector T odefinition mode identical.A t ohealth status and the A of evaluation object ohealth status (the evaluation object initial damage vector d of evaluation object orepresent) identical, A in cyclic process t othe health status of evaluation object use all the time evaluation object initial damage vector d orepresent.T t oand d oa t oparameter, C t oby A t omechanics Calculation result composition.
The 4th step: in Cable Structure military service process, the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical).Obtain current cable structure steady temperature data vector T in actual measurement ttime, namely obtaining current cable structure steady temperature data vector T tthe synchronization in moment, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, composition " monitored amount current value vector C ".
Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, actual measurement obtains all M in Cable Structure 1the rope force data of root support cable, all these rope force data composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, Actual measurement obtains all M 1the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, all M 1two supporting end points horizontal range data of root support cable form current support cable two and support end points horizontal range vector l t x, current support cable two supports end points horizontal range vector l t xthe coding rule and initial rope force vector F of element othe coding rule of element identical.
The 5th step: according to current cable structure steady temperature data vector T t, upgrade where necessary current initial mechanical calculating benchmark model A t o, the current initial value of monitored amount vector C t owith current initial Cable Structure steady temperature data vector T t o.Obtain current cable structure steady temperature data vector T in the 4th step actual measurement tafter, relatively T tand T t oif, T tequal T t o, do not need A t oand T t oupgrade, otherwise need to be to A t oand T t oupgrade, the step that update method specifies by technical scheme is carried out.
The 6th step: the step specifying by technical scheme, at current initial mechanical calculating benchmark model A t obasis on carry out several times Mechanics Calculation, obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D by calculating u.Concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d oon the basis of the existing damage of this support cable representing, increase again unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the generalized displacement component of a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d oon the basis of the existing generalized displacement of this bearing representing, there is again unit generalized displacement (if for example this evaluation object is the translational component of the x direction of a bearing, just suppose that this bearing has unit line displacement in x direction, for example get 1mm, if this evaluation object is the angular displacement component around x axle of a bearing, just suppose that this bearing is around the angular displacement of x Zhou You unit, for example get 100,000/radian), if this evaluation object is a load, just suppose that this load is at vectorial d oon the basis of the existing variable quantity of this load representing, increasing load unit changes (if this load is distributed load again, and this distributed load is line distributed load, load unit changes 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 face distributed load, and load unit changes can get 1MPa, 2MPa, 3MPa or 1kNm/m 2, 2kNm/m 2, 3kNm/m 2deng being unit change; If this load is centre-point load, and this centre-point load is couple, and load unit changes can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centre-point load, and this centre-point load is concentrated force, and load unit changes can get 1kN, 2kN, 3kN etc. for unit change; If this load is volume load, load unit changes can get 1kN/m 3, 2kN/m 3, 3kN/m 3deng being unit change).
The 7th step: set up linear relationship error vector e and vectorial g.Utilize data (the current initial value vector of monitored amount C above t o, unit damage monitored numerical quantity transformation matrices Δ C), when the 6th step is calculated each time, only have in evaluation object the increase unit damage of an evaluation object or unit generalized displacement or load unit to change D calculating each time hypothesis uk, the evaluation object that increases unit damage or unit generalized displacement or load unit variation in calculating is each time different from the evaluation object that increases unit damage or unit generalized displacement or load unit variation in other calculating, calculate each time the current value of all utilizing mechanics method (for example adopting finite element method) to calculate all monitored amounts in Cable Structure, when calculating each time a monitored amount calculation current vector C of composition, calculate each time vectorial d of damage of composition, originally walking out of existing damage vector d only uses in this step, damaging in all elements of vectorial d only has the numerical value of an element to get D uk, the numerical value of other element gets 0, damages coding rule and the vectorial d of the element of vectorial d othe coding rule of element identical, by C, C t o, Δ C, D u, d brings formula (1) into, obtains a linear relationship error vector e, calculates each time a linear relationship error vector e, there is N evaluation object just to have N calculating, just there is N linear relationship error vector e, will this N linear relationship error vector e obtain a vector after being added, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e.Vector g equals final error vector e.
e = abs ( ΔC · d - C + C o t ) - - - ( 1 )
In formula (1), abs () is the function that takes absolute value, and each vectorial element of trying to achieve in bracket is taken absolute value.
The 8th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: space coordinate monitoring system, signal (data) collector, the computing machine and the panalarm of communicating by letter of monitored amount monitoring system (for example, containing strain measurement system, signal conditioner etc.), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and the supporting end points of Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), support cable cable force monitoring system, support cable.The volume coordinate of the Suo Li of each monitored amount, each temperature, each root support cable, the supporting end points of each root support cable must arrive by monitored system monitoring, and the signal monitoring is transferred to signal (data) collector by monitoring system; Signal is delivered to computing machine through signal picker; Computing machine is responsible for the health monitoring software of the evaluation object of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; In the time monitoring evaluation object health status and change, computer control communication panalarm is reported to the police to the personnel of monitor staff, owner and (or) appointment.
The 9th step: by current monitored amount initial value vector C t o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D uparameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
The tenth step: establishment installation and operation this method system software on computers, this software will complete the functions (being all work that can complete with computing machine in this specific implementation method) such as monitoring that this method required by task wants, record, control, storage, calculating, notice, warning
The 11 step: the monitored amount current value vector of foundation C is with the current initial value vector of monitored amount C t o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D uand the vectorial d(of the current name damage of evaluation object is made up of all Suo Dangqian name amount of damage) between the linear approximate relationship (formula (2)) that exists, calculate the noninferior solution of the vectorial d of the current name damage of evaluation object according to multi-objective optimization algorithm, namely with reasonable error but can reflect more exactly the solution of the variation of the health status of evaluation object.Can adopt Objective Programming (Goal Attainment Method) in multi-objective optimization algorithm to solve formula (2) and obtain the vectorial d of current damage.
C = C o t + ΔC · d - - - ( 8 )
The 12 step: utilize the current actual damage vector of evaluation object d ak element d a kwith evaluation object initial damage vector d ok element d okk the element d with the vectorial d of the current name damage of evaluation object kbetween relation, calculate the current actual damage of evaluation object vector d aall elements.D a krepresent the current actual health status of k evaluation object, if this evaluation object is support cable, so a d in cable system a krepresent its current actual damage, d a kbe to represent that its corresponding support cable is without health problem at 0 o'clock, d a knumerical value is not the support cable that represents that its corresponding support cable is unsoundness problem at 0 o'clock, and the support cable of unsoundness problem may be slack line, also may be damaged cable, its numerical response the degree of lax or damage; If this evaluation object is generalized displacement component, so a d of a bearing a krepresent its current actual generalized displacement numerical value; If this evaluation object is load, so a d a krepresent that it is with respect to setting up initial mechanical calculating benchmark model A otime the structure corresponding load of bearing variable quantity.
The 13 step: by current evaluation object actual damage vector d ain the M relevant to support cable 1individual element takes out, the current actual damage vector of composition support cable d ca, the current actual damage vector of support cable d cathe coding rule and initial rope force vector F of element othe coding rule of element identical.The current actual damage vector of support cable d cah element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M 1; The current actual damage vector of support cable d camiddle numerical value is not 0 the element support cable corresponding to unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable by lossless detection method, remaining is exactly slack line, namely needs to adjust the rope of Suo Li, and these ropes that need adjust Suo Li are at the current actual damage vector of support cable d cain corresponding element numerical value (for example one of them element can be used d ca hrepresent) represent and the degree of injury of the relax level mechanics equivalence of these support cables.Damaged cable is at the current actual damage vector of support cable d cathe numerical value of the element of middle correspondence just represents its degree of injury, the numerical value of corresponding element represents while being 100% that this support cable thoroughly loses load-bearing capacity, the load-bearing capacity that represents this support cable forfeiture corresponding proportion between 0 and 100% time, has so far just identified damaged cable and degree of injury thereof.
The 14 step: count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, particularly can be in the hope of the relax level of these ropes (being the long adjustment amount of rope) according to formula (3).So just realize the lax identification of support cable.So far damaged cable and slack line have just all been identified.
Δl h t = d h ca 1 - d h ca F h [ E h t 1 + ( ω h t l xh t ) 2 A h t E h t 12 ( F h ) 3 ] A h t + F h l oh t - - - ( 3 )
E in formula (3) t hthe current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure t owhen expression, the elastic modulus of h root support cable, A t hthe current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure t owhen expression, the cross-sectional area of h root support cable, F hthe current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure t owhen expression, the current cable power of h root support cable, d ca hthe current actual damage degree of h root support cable, ω t hthe current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure t owhen expression, the weight of the unit length of h root support cable, l t xhthe current initial Cable Structure steady temperature data vector T for steady temperature data in Cable Structure t owhen expression, the horizontal range of two supporting end points of h root support cable, l t xhthat current support cable two supports end points horizontal range vector l t xan element, current support cable two supports end points horizontal range vector l t xthe coding rule and initial drift vector l of element othe coding rule of element identical, E t hcan obtain according to the characteristic material data of looking into or survey h root support cable A t hand ω t hcan be according to the thermal expansivity of h root support cable, A oh, ω oh, F h, T oand T t oobtain by conventional physics and Mechanics Calculation.
The 15 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by personnel's operational health monitoring system.
The 16 step: under specified requirements, the computing machine automatic operation communication panalarm in health monitoring systems is reported to the police to the personnel of monitor staff, owner and (or) appointment.
The 17 step: get back to the 4th step, start the circulation by the 4th step to the 17 steps.

Claims (1)

1. strain monitoring problem rope load generalized displacement recognition methods, is characterized in that described method comprises:
A. for sake of convenience, it is evaluation object that this method unitedly calls evaluated generalized displacement of support component, support cable and load, if the quantity sum of the quantity of evaluated generalized displacement of support component, the quantity of support cable and load is N, the quantity of evaluation object is N; Determine the coding rule of evaluation object, by this rule, by evaluation object numberings all in Cable Structure, this numbering will be used for generating vector sum matrix in subsequent step; This method represents this numbering with variable k, k=1, and 2,3 ..., N; Determine the point being monitored of specifying, point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified point numberings; Determine monitored should the changing direction of point being monitored, and number to the monitored strain of all appointments, " monitored strain numbering " will be used for generating vector sum matrix in subsequent step, and " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; In this method, must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring, the moment of survey record data is called the physical record data moment; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load claims again surface load, is the load that acts on body surface, comprises two kinds of centre-point load and distributed loads; Volume load is that continuous distribution is in the load of interior of articles each point, including the deadweight and inertial force of object; Centre-point load is divided into two kinds of concentrated force and concentrated couples, comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now variation of load is embodied as the variation of a concentrated force component or a concentrated couple component; Distributed load is divided into line distributed load and face distributed 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 of distributed load is expressed by distribution intensity, and distribution intensity is expressed by distribution characteristics and amplitude; If load is actually distributed load, when this method is talked about the variation 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; Comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, a distributed load can be resolved into three components, if the amplitude of the distribution intensity separately of three components of this distributed load changes, and the ratio changing is all not identical, so in the method three components of this distributed load being counted or added up is three distributed loads, and now load just represents the one-component 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, the size of volume load is expressed by distribution intensity, distribution intensity is expressed by distribution characteristics and amplitude; If load is actually volume load, actual treatment is the change of the amplitude of volume load distribution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, while now mentioning in the method the change of load, in fact refer to the change of the amplitude of the distribution intensity of volume load, now, the load changing refers to the volume load that the amplitude of those distribution intensities changes; Comprising in the coordinate system of Descartes's rectangular coordinate in tying up to, one individual stow lotus can be resolved into three components, if the amplitude of the distribution intensity separately of three components of this volume load changes, and the ratio changing is all not identical, and so in the method three components of this volume load 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 Cable Structure composition material and Cable Structure environment of living in, utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, statistics obtains interior during this period of time cloudy quantity and is designated as T cloudy day, in the method can not be seen to one of the sun daytime and be called all day the cloudy day, statistics obtain each cloudy day in T cloudy day 0 after the sunrise moment next day highest temperature and the lowest temperature between 30 minutes, the sunrise moment on the meteorology that the sunrise moment refers to base area revolutions and the rule that revolves round the sun is definite, do not represent necessarily can see the same day sun, can inquire about data or calculate sunrise moment of each required day by conventional meteorology, each cloudy day 0 after the sunrise moment next day highest temperature between 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day r, between inquiry Cable Structure location and Altitude Region, place, be no less than temperature that the meteorological data in recent years of 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and Altitude Region, place about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of narration hunit be DEG C/m, on the surface of Cable Structure, get " R Cable Structure surface point ", get the Specific Principles of " R Cable Structure surface point " narrates in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim that the temperature data calculating is " R Cable Structure surface temperature computational data ", from the residing minimum height above sea level of Cable Structure to the highest height above sea level, in Cable Structure, uniform choosing is no less than three different sea level elevations, the sea level elevation place choosing at each, at least choose two points at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in the measurement Cable Structure of choosing along comprising the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure in the direction of the Temperature Distribution of wall thickness, measure Cable Structure along each and be no less than three points along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, measure all temperature that are selected a little, the temperature recording is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure the direction of Cable Structure along the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured the direction of Cable Structure along the Temperature Distribution of wall thickness, direction along each measurement Cable Structure along the Temperature Distribution of wall thickness has been chosen E point in Cable Structure, wherein H and E are not less than 3, B is not less than 2, if HBE is the product of H and B and E, corresponding total HBE " measuring the point of Cable Structure along the temperature profile data of thickness ", to obtain by actual measurement the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " below, claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, just claim that the temperature data calculating is " HBE Cable Structure is along thickness temperature computation data ", if BE is the product of B and E, total BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in sea level elevation place of choosing at each in this method, measure temperature in Cable Structure location according to meteorology and require to choose a position, will obtain the temperature of the Cable Structure place environment that meets the requirement of meteorology measurement temperature in this position actual measurement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of getable this day, at the flat board of a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, the one side of this reference plate on the sunny side, be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained to the temperature of the sunny slope of reference plate,
B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains the temperature record of the Cable Structure place environment that meets the requirement of meteorology measurement temperature simultaneously, obtain being carved at sunrise the same day temperature measured data sequence of the Cable Structure place environment between 30 minutes after sunrise moment next day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data that was carved at sunrise the Cable Structure place environment between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains Cable Structure place environment by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T emax, temperature measured data sequence by Cable Structure place environment obtains the temperature of Cable Structure place environment about the rate of change of time by conventional mathematical computations, and this rate of change is also along with the time changes, obtain being carved at sunrise the same day measured data sequence of the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data that was carved at sunrise the temperature of the sunny slope of the reference plate between 30 minutes after the sunrise moment next day same day, find maximum temperature and minimum temperature in the measured data sequence of temperature of the sunny slope of reference plate, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of the sunny slope of reference plate by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T pmax, obtain being carved at sunrise the same day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between 30 minutes after sunrise moment next day by Real-Time Monitoring, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by being carved at sunrise the Cable Structure surface temperature measured data between 30 minutes after the sunrise moment next day same day of a Cable Structure surface point, find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, deduct and be carved at sunrise the maximum temperature difference between 30 minutes after sunrise moment next day on same day that minimum temperature obtains the temperature of each Cable Structure surface point by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T smax, obtain the temperature of each Cable Structure surface point about the rate of change of time by each Cable Structure surface temperature measured data sequence by conventional mathematical computations, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the sea level elevation place that chooses at each and amount to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature, the absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", choose H different sea level elevation and just had H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax,
B3: measure and calculate acquisition Cable Structure steady temperature data, first, determine the moment that obtains Cable Structure steady temperature data, the condition relevant to the moment that determines acquisition Cable Structure steady temperature data has six, Section 1 condition be obtain Cable Structure steady temperature data moment after being carved at sunset sunrise moment next day between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the definite meteorology of revolution rule, can inquire about data or calculate sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, reference plate maximum temperature difference Δ T pmaxwith Cable Structure surface maximum temperature difference Δ T smaxall be not more than 5 degrees Celsius, the b condition of Section 2 condition be after being carved at sunrise sunrise moment next day on same day between 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates above emaxbe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T pmaxafter deducting 2 degrees Celsius, be not more than Δ T emax, and Cable Structure surface maximum temperature difference Δ T smaxbe not more than Δ T pmax, only need meet in a condition of Section 2 and b condition one is just called and meets Section 2 condition, Section 3 condition is that the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 4 condition is that the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data, Section 5 condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is to be carved at sunrise the minimal value between 30 minutes after the sunrise moment next day same day, Section 6 condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius, this method is utilized above-mentioned six conditions, any one in following three kinds of moment is called to " the mathematics moment that obtain Cable Structure steady temperature data ", the first moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 5 condition, the second moment is the moment that only meets the Section 6 condition in above-mentioned " condition relevant to the moment that determines acquisition Cable Structure steady temperature data ", the third moment is to meet Section 1 in above-mentioned " condition relevant to the moment that determines the to obtain Cable Structure steady temperature data " moment to Section 6 condition simultaneously, be exactly in this method when one in the physical record data moment when obtaining the mathematics moment of Cable Structure steady temperature data, the moment that obtains Cable Structure steady temperature data is exactly the mathematics moment that obtains Cable Structure steady temperature data, be not any moment in the physical record data moment in this method if obtain the mathematics moment of Cable Structure steady temperature data, getting this method is the moment that obtains Cable Structure steady temperature data close to moment of those physical record data in the mathematics moment that obtains Cable Structure steady temperature data, this method will be used the amount at the moment survey record that obtains Cable Structure steady temperature data to carry out the relevant health monitoring analysis of Cable Structure, this method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, i.e. the not temporal evolution of Cable Structure temperature in this moment, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to Cable Structure heat transfer characteristic, utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data, utilize the thermal conduction study computation model of Cable Structure, obtain obtaining the Temperature Distribution of Cable Structure in moment of Cable Structure steady temperature data by conventional Calculation of Heat Transfer, now calculate by stable state in the temperature field of Cable Structure, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", in the time of R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculating is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure stable state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", get " R Cable Structure surface point " on the surface of Cable Structure time, the quantity of " R Cable Structure surface point " and necessary three conditions that meet that distribute, first condition is when Cable Structure temperature field is during in stable state, when the temperature of any point on Cable Structure surface be by " R Cable Structure surface point " in Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point while obtaining, on the Cable Structure surface that linear interpolation obtains, on the temperature of this arbitrfary point and Cable Structure surface, the error of the actual temperature of this arbitrfary point is not more than 5%, Cable Structure surface comprises support cable surface, second condition is that in " R Cable Structure surface point ", the quantity at the point of same sea level elevation is not less than 4, and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 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 hthe numerical value obtaining, gets Δ T for convenience of narration hunit be DEG C/m that the unit of getting Δ h for convenience of narration is m, " R Cable Structure surface point " refers to while only considering sea level elevation along the definition of adjacent Cable Structure surface point between two of sea level elevation, in " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, the 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to the geometric properties of Cable Structure and bearing data, in Cable Structure, find and be subject to the sunshine-duration position of those surface points the most fully the whole year, in " R Cable Structure surface point ", having a Cable Structure surface point at least is the annual point being subject in the most sufficient those surface points of sunshine-duration in Cable Structure,
C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " the Cable Structure steady temperature data that calculate under original state, Cable Structure steady temperature data under original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T o", actual measurement or consult reference materials and obtain the temperature variant physical and mechanical properties parameter of the various materials that Cable Structure uses, obtain initial Cable Structure steady temperature data vector T in actual measurement osynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F o, according to data including Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is the length of 0 o'clock, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables while obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing osuo Li under condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the weight vector of the initial drift vector of support cable, the initial free unit length of initial free cross-sectional area vector sum, the coding rule of the element of the weight vector of the initial drift vector of support cable, the initial free unit length of initial free cross-sectional area vector sum and initial rope force vector F othe coding rule of element identical, obtain T in actual measurement otime, namely obtaining initial Cable Structure steady temperature data vector T othe synchronization in moment, directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure is to comprise Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, in obtaining the measured data of initial Cable Structure, measurement calculates the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C o, monitored amount initial value vector C othe coding rule of coding rule and M monitored amount identical, utilize support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement data to set up evaluation object initial damage vector d o, vectorial d orepresent with initial mechanical calculating benchmark model A othe initial health of the evaluation object of the Cable Structure representing, evaluation object initial damage vector d oelement number equal N, d oelement and evaluation object be one-to-one relationship, vectorial d othe coding rule of element identical with the coding rule of evaluation object, if d oevaluation object corresponding to some elements be support cable, so a d in cable system othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the corresponding support cable of this element is intact, do not damage, if its numerical value is 100%, represent that the corresponding support cable of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, represents this support cable, , if d oevaluation object corresponding to some elements be some generalized displacement components of some bearings, d so othe numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d oevaluation object corresponding to some elements be some load, in this method, get d othis element numerical value be 0, the initial value that represents the variation of this load is 0, if there is no the initial generalized displacement measurement data of Cable Structure bearing or can think that the initial generalized displacement of Cable Structure bearing is at 0 o'clock, vectorial d oin the each element numerical value relevant to Cable Structure generalized displacement of support get 0, if while not having the Non-destructive Testing Data of support cable and other can express the data of health status of support cable, or can think that structure original state is not damaged during without relaxed state, vectorial d oin the each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data refer to the bearing generalized coordinate data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing refers to setting up initial mechanical calculating benchmark model A otime, the generalized displacement that Cable Structure bearing occurs with respect to the bearing under Cable Structure design point,
The temperature variant physical and mechanical properties parameter of the various materials that d. use according to measured data, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the Cable Structure of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, initial Cable Structure steady temperature data vector T owith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " o, based on A othe Cable Structure computational data calculating must approach its measured data very much, and difference therebetween must not be greater than 5%; Corresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A oevaluation object initial damage vector d for evaluation object health status orepresent; Corresponding to A omonitored amount initial value vector C for the initial value of all monitored amounts orepresent; Set up for the first time the current initial mechanical calculating benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure t owith the current initial value vector of monitored amount C t otime, the current initial mechanical calculating benchmark model A of Cable Structure t ojust equal the initial mechanical calculating benchmark model A of Cable Structure o, the current initial value vector of monitored amount C t ojust equal monitored amount initial value vector C o; A t ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure t otime, T t ojust equal T o; A t oinitial health and the A of evaluation object othe health status of evaluation object identical, also use evaluation object initial damage vector d orepresent A in cyclic process below t othe initial health of evaluation object use all the time evaluation object initial damage vector d orepresent; T oand d oa oparameter, by A oinitial value and the C of all monitored amounts of obtaining of Mechanics Calculation result othe initial value of all monitored amounts that represent is identical, therefore also can say C oby A omechanics Calculation result composition; T t oand d oa t oparameter, C t oby A t omechanics Calculation result composition;
E. from entering the circulation that is walked n step by e here; In structure military service process, constantly according to " the temperature survey calculating method of the Cable Structure of this method " the constantly current data of Actual measurement acquisition " Cable Structure steady temperature data ", the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical; Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, actual measurement obtains all M in Cable Structure 1the rope force data of root support cable, all these rope force data composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement tsynchronization, Actual measurement obtains all M 1the volume coordinate of two supporting end points of root support cable, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports the coding rule and initial rope force vector F of the element of end points horizontal range vector othe coding rule of element identical;
F. according to current cable structure steady temperature data vector T t, upgrade current initial mechanical calculating benchmark model A according to step f1 to f3 t o, the current initial value of monitored amount vector C t owith current initial Cable Structure steady temperature data vector T t o;
F1. compare T twith T t oif, T tequal T t o, A t o, C t oand T t oremain unchanged; Otherwise need to follow these steps to A t oand T t oupgrade;
F2. calculate T twith T opoor, T twith T odifference be exactly the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T opoor represent with steady temperature change vector S, S equals T tdeduct T o, S represents the variation of Cable Structure steady temperature data;
F3. to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation applying is just taken from steady temperature change vector S, to A oin the temperature variation that applies of Cable Structure after obtain the current initial mechanical calculating benchmark model A that upgrades t o, upgrade A t otime, T t oall elements numerical value is also used T tcorresponding replacement of all elements numerical value, upgraded T t o, so just obtained correctly corresponding to A t ot t o; Upgrade C t omethod be: when upgrade A t oafter, obtain A by Mechanics Calculation t oin concrete numerical value all monitored amounts, current, these concrete numerical value compositions C t o; A t othe initial health of support cable use all the time evaluation object initial damage vector d orepresent;
G. at current initial mechanical calculating benchmark model A t obasis on carry out several times Mechanics Calculation according to step g 1 to g4, obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D by calculating u;
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C constantly updates, and is upgrading current initial mechanical calculating benchmark model A t o, the current initial value of monitored amount vector C t owith current initial Cable Structure steady temperature data vector T t oafterwards, must then upgrade Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D u;
G2. at the current initial mechanical calculating benchmark model A of Cable Structure t obasis on carry out several times Mechanics Calculation, on calculation times numerical value, equal the quantity N of all evaluation objects, have N evaluation object just to have N calculating; According to the coding rule of evaluation object, calculate successively; Calculating each time hypothesis only has an evaluation object on the basis of original damage or generalized displacement or load, to increase unit damage or unit generalized displacement or load unit variation again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d oon the basis of the existing damage of this support cable representing, increase again unit damage, if this evaluation object is the generalized displacement component of a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, if this evaluation object is a load, just suppose that this load is at vectorial d oon the basis of the existing variable quantity of this load representing, increase again load unit and change, use D ukunit damage or unit generalized displacement or the load unit variation of recording this increase, wherein k represents the numbering of the evaluation object that increases unit damage or unit generalized displacement or load unit variation, D ukevaluation object unit change vector D uan element, evaluation object unit change vector D ucoding rule and the vectorial d of element othe coding rule of element identical; The evaluation object that increases unit damage or unit generalized displacement or load unit variation in calculating is each time different from the evaluation object that increases unit damage or unit generalized displacement or load unit variation in other calculating, calculate each time the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure, a monitored amount calculation current vector of current calculated value composition of all monitored amounts that calculate each time, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C oelement coding rule identical;
G3. the monitored amount calculation current vector calculating each time deducts the current initial value vector of monitored amount C t oobtain a vector, again each element of this vector is calculated to unit damage or unit generalized displacement or the load unit variation numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N evaluation object just to have N monitored amount unit change vector;
G4. by this N monitored amount unit change vector according to the coding rule of N evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is listed as successively; Each of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is listed as corresponding to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is the different unit change amplitude in the time that different evaluation objects increase unit damage or unit generalized displacement or load unit variation corresponding to same monitored amount; Coding rule and the vectorial d of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C othe coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
H. obtain current cable structure steady temperature data vector T in actual measurement ttime, actual measurement obtains obtaining current cable structure steady temperature data vector T tthe current measured value of all monitored amounts of Cable Structure of synchronization in moment, form monitored amount current value vector C; The current initial value vector of monitored amount current value vector C and monitored amount C t owith monitored amount initial value vector C odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at concrete numerical value in the same time not;
I. define the vectorial d of the current name damage of evaluation object, the element number of the vectorial d of the current name damage of evaluation object equals the quantity of evaluation object, between the element of the vectorial d of the current name damage of evaluation object and evaluation object, be one-to-one relationship, the element numerical value of the vectorial d of the current name damage of evaluation object represents nominal degree of injury or nominal generalized displacement or the nominal load variable quantity of corresponding evaluation object; Coding rule and the vectorial d of the element of vector d othe coding rule of element identical;
J. the monitored amount current value vector of foundation C is with the current initial value vector of monitored amount C t o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and the vectorial d of the current name damage of evaluation object to be asked, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate the vectorial d of the current name damage of evaluation object;
C = C o t + ΔC · d Formula 1
K. define the current actual damage vector of evaluation object d a, the current actual damage vector of evaluation object d aelement number equal the quantity of evaluation object, the current actual damage vector of evaluation object d aelement and evaluation object between be one-to-one relationship, the current actual damage of evaluation object vector d aelement numerical value represent actual damage degree or actual generalized displacement or the real load variable quantity of corresponding evaluation object; Vector d acoding rule and the vectorial d of element othe coding rule of element identical;
L. the current actual damage vector of the evaluation object d that utilizes formula 2 to express ak element d a kwith evaluation object initial damage vector d ok element d okk the element d with the vectorial d of the current name damage of evaluation object kbetween relation, calculate the current actual damage of evaluation object vector d aall elements;
Figure FDA0000474929550000112
formula 2
K=1 in formula 2,2,3 ...., N, d a krepresent the current actual health status of k evaluation object, d a kbe to represent that k evaluation object is without health problem at 0 o'clock, d a knumerical value is not to represent that k evaluation object is the evaluation object of unsoundness problem at 0 o'clock, if this evaluation object is support cable, so a d in cable system a krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d a knumerical response the degree of lax or damage of this support cable; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, the current actual damage vector of evaluation object d ain with slack line corresponding to element numerical expression be and the current actual equivalent damage degree of slack line relax level mechanics equivalence; If this evaluation object is generalized displacement component, so a d of a bearing a krepresent its current actual generalized displacement numerical value; If this evaluation object is load, so a d a krepresent the actual change amount of this load;
M. utilize at current cable structure steady temperature data vector T tunder condition, walk the slack line that identifies and with the current actual damage vector of evaluation object d at l athese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, utilize obtain in e step at current cable structure steady temperature data vector T tcurrent cable force vector F under condition and current support cable two support end points horizontal range vector, utilize obtain in c step at initial Cable Structure steady temperature data vector T othe weight vector of the initial drift vector of the support cable under condition, the initial free unit length of initial free cross-sectional area vector sum, initial rope force vector F o, utilize current cable structure steady temperature data vector T tthe current steady temperature data of support cable that represent, utilize obtain in c step at initial Cable Structure steady temperature data vector T othe support cable initial steady state temperature data representing, utilize the temperature variant physical and mechanical properties parameter of the various materials that use in the Cable Structure of c step acquisition, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, mechanics equivalent condition is: one, the mechanics parameters without lax initial drift, geometrical property parameter, density and material during with not damaged of the rope of two equivalences is identical; Two, after lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; While meeting above-mentioned two mechanics equivalent conditions, the mechanics function of two such support cables in Cable Structure is exactly identical, if replaced after damaged cable with equivalent slack line, any variation can not occur Cable Structure, and vice versa; The relax level of trying to achieve those and be judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of those support cables that need adjust Suo Li; Lax identification and the damage identification of support cable are so just realized; When calculating, institute's demand power is provided by current cable force vector F corresponding element; This method is referred to as damaged cable and slack line the support cable of unsoundness problem, referred to as problem rope, so this method is according to the current actual damage vector of evaluation object d acan either identify problem rope, can define which bearing generalized displacement and numerical value thereof have occurred, also can define which load the numerical value that changes and change has occurred; So far this method has realized and has rejected problem rope identification impact, Cable Structure that generalized displacement of support, load change and structure temperature change, realize simultaneously and rejected generalized displacement of support, structure temperature variation and identification support cable health status variable effect, load change amount, also realized and rejected load change, structure temperature variation and identification variable effect, generalized displacement of support of support cable health status;
N. get back to e step, start to be walked by e the circulation next time of n step.
CN201410085856.3A 2014-03-10 2014-03-10 Strain monitoring error cable load generalized displacement identification method Pending CN103852286A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410085856.3A CN103852286A (en) 2014-03-10 2014-03-10 Strain monitoring error cable load generalized displacement identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410085856.3A CN103852286A (en) 2014-03-10 2014-03-10 Strain monitoring error cable load generalized displacement identification method

Publications (1)

Publication Number Publication Date
CN103852286A true CN103852286A (en) 2014-06-11

Family

ID=50860204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410085856.3A Pending CN103852286A (en) 2014-03-10 2014-03-10 Strain monitoring error cable load generalized displacement identification method

Country Status (1)

Country Link
CN (1) CN103852286A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104990587A (en) * 2015-07-23 2015-10-21 东南大学 Simplified method of recognizing problematic load cable generalized displacement through strain monitoring
CN105067334A (en) * 2015-07-23 2015-11-18 东南大学 Method for recognizing load generalized displacement of damaged cable based on streamlined strain monitoring process
CN105115768A (en) * 2015-07-23 2015-12-02 东南大学 Identification method for damaged cable, load and generalized displacement through simplified strain monitoring

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006337144A (en) * 2005-06-01 2006-12-14 Kawasaki Heavy Ind Ltd Fatigue life diagnostic method and diagnostic support device of bridge
CN102128725A (en) * 2010-12-02 2011-07-20 李惠 Method for monitoring health and pre-warning safety of large-span space structure
CN103268371A (en) * 2013-04-26 2013-08-28 重庆交通大学 Real-time bridge load identification method based on influence matrix
CN103604629A (en) * 2013-12-09 2014-02-26 东南大学 Damaged cable/concentrated load/bracket angular displacement progressive recognition method on basis of strain monitoring
CN103616206A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring identification method for defective cable and concentrated loads
CN103616244A (en) * 2013-12-09 2014-03-05 东南大学 Strain-monitoring damaged cable centralized load generalized displacement recognition method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006337144A (en) * 2005-06-01 2006-12-14 Kawasaki Heavy Ind Ltd Fatigue life diagnostic method and diagnostic support device of bridge
CN102128725A (en) * 2010-12-02 2011-07-20 李惠 Method for monitoring health and pre-warning safety of large-span space structure
CN103268371A (en) * 2013-04-26 2013-08-28 重庆交通大学 Real-time bridge load identification method based on influence matrix
CN103604629A (en) * 2013-12-09 2014-02-26 东南大学 Damaged cable/concentrated load/bracket angular displacement progressive recognition method on basis of strain monitoring
CN103616206A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring identification method for defective cable and concentrated loads
CN103616244A (en) * 2013-12-09 2014-03-05 东南大学 Strain-monitoring damaged cable centralized load generalized displacement recognition method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭良友等: "武汉长江二桥的索力、温度和应力测量", 《桥梁建设》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104990587A (en) * 2015-07-23 2015-10-21 东南大学 Simplified method of recognizing problematic load cable generalized displacement through strain monitoring
CN105067334A (en) * 2015-07-23 2015-11-18 东南大学 Method for recognizing load generalized displacement of damaged cable based on streamlined strain monitoring process
CN105115768A (en) * 2015-07-23 2015-12-02 东南大学 Identification method for damaged cable, load and generalized displacement through simplified strain monitoring

Similar Documents

Publication Publication Date Title
CN103913342A (en) Method for progressively recognizing fault cable, load and generalized displacement based on angle monitoring
CN103868731A (en) Generalized displacement strain monitoring damaged cable load identification method
CN103884524A (en) Method for recognizing fault rod loads based on generalized displacement angle monitoring
CN103913338A (en) Hybrid monitoring defective cable load generalized displacement identification method
CN103852285A (en) Recognition method for problematic cable loads through generalized displacement and cable force monitoring
CN103852331A (en) Progressive recognition method for problem cable load generalized displacement based on mixing monitoring
CN103852289A (en) Problem cable load generalized displacement recognition method based on space coordinate monitoring
CN103913330A (en) Generalized displacement space coordinate monitoring damaged cable load identification method
CN103852291A (en) Progressive recognition method for problematic cable loads through generalized displacement and space coordinate monitoring
CN103852327A (en) Method for recognizing damaged cables, loads and generalized displacement based on cable force monitoring
CN103852329A (en) Damaged cable load generalized displacement recognition method through mixed monitoring
CN103852307A (en) Damaged cable load generalized displacement recognition method based on space coordinate monitoring
CN103868743A (en) Generalized displacement strain monitoring problematical cable load progressive identification method
CN103852333A (en) Strain monitoring damaged cable load linear displacement identification method
CN103852286A (en) Strain monitoring error cable load generalized displacement identification method
CN103913339A (en) Method for recognizing fault cable, load and generalized displacement based on angle monitoring
CN103852290A (en) Recognition method for problematic cable loads through generalized displacement and strain monitoring
CN103852316A (en) Recognition method for problematic cable loads through generalized displacement and space coordinate monitoring
CN103868717A (en) Strain monitoring-based identification method for defective cable, load and angular displacement
CN103852296A (en) Strain monitoring error cable load linear displacement identification method
CN103913321A (en) Generalized displacement hybrid monitoring defective cable load identification method
CN103868737A (en) Angle monitoring problem cable load linear displacement identification method
CN103852292A (en) Strain monitoring error cable load generalized displacement progressive identification method
CN103868733A (en) Linear displacement angle monitoring-based load identification method for defective cable
CN103853926A (en) Space coordinate monitoring problem cable load generalized displacement progressive identification method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20140611