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

Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process Download PDF

Info

Publication number
CN105067342A
CN105067342A CN201510439657.2A CN201510439657A CN105067342A CN 105067342 A CN105067342 A CN 105067342A CN 201510439657 A CN201510439657 A CN 201510439657A CN 105067342 A CN105067342 A CN 105067342A
Authority
CN
China
Prior art keywords
cable structure
data
temperature
cable
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201510439657.2A
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 CN201510439657.2A priority Critical patent/CN105067342A/en
Publication of CN105067342A publication Critical patent/CN105067342A/en
Pending legal-status Critical Current

Links

Abstract

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

Description

Simplify hybrid monitoring problem cable load generalized displacement recognition methods
Technical field
Cable-stayed bridge, suspension bridge, the structures such as truss structure have a common ground, it is exactly that they have many parts for bearing tensile load, such as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of the class formation is with rope, cable is subjected only to the rod member of tensile load for supporting member, for convenience, such structure representation is " Cable Structure " by this method, and by all carrying ropes of Cable Structure, carry cable, and all rod members (also known as two power rod members) for being subjected only to axial tension or axial compression load, it is collectively referred to as convenience " cable system ", in this method carrying rope is censured with " support cable " this noun, carrying cable and the rod member for being subjected only to axial tension or axial compression load, sometimes referred to simply as " rope ", so when " rope " this word is used below, two power rod members are just actually referred to truss structure.The impaired and lax pair Cable Structure of support cable is a significant threat safely, and damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method.During structure military service, the correct identification to support cable or the health status of cable system is related to the safety of whole Cable Structure.When environment temperature changes, the temperature of Cable Structure typically also can be with changing, when Cable Structure temperature changes, generalized displacement may occur for Cable Structure bearing, the load that Cable Structure is born is it can also happen that change, even if the temperature of actually Cable Structure does not change, the load that Cable Structure is born may also individually change, the health status of Cable Structure may also change simultaneously, in this complex condition, based on hybrid monitoring, (this method judges the health status of Cable Structure to this method by the hybrid monitoring of the change of the measurable parameter to the foregoing different types of Cable Structure of this section, all monitored Cable Structure characteristic parameters are referred to as " monitored amount " by this method, it is made up of due to being now monitored amount the different types of measurable parameter mixing of Cable Structure, this method is referred to as hybrid monitoring) recognize generalized displacement of support and problem cable, belong to engineering structure health monitoring field.
Background technology
Reject load change, Cable Structure generalized displacement of support and structure temperature and change influence to Cable Structure health status recognition result, so as to recognize the health status of support cable exactly, the problem of be in the urgent need to address at present;Reject load change, the change of Cable Structure health status and structure temperature and change influence to Cable Structure generalized displacement of support recognition result, so as to recognize Cable Structure generalized displacement of support exactly, the problem of be also in the urgent need to address at present;This method discloses a kind of effective, cheap method for solving the two problems.
The content of the invention
This method discloses a kind of method, under conditions of cost is lower, realizes two kinds of functions, is respectively:First, the influence of generalized displacement of support, load change and structure temperature change to Cable Structure health status recognition result is rejected, so as to identify the health status of support cable exactly;2nd, this method can also reject the influence of load change, the change of Cable Structure health status and structure temperature change to Cable Structure generalized displacement of support recognition result, so as to recognize Cable Structure generalized displacement of support exactly.
During Cable Structure military service, Suo Changdu under support cable free state (it is 0 that now rope tensility, which is also referred to as Suo Li) (is referred to as drift, this method refers exclusively to the drift of that section of rope between the supporting end points of support cable two) it can change, the first purpose of this method seeks to identify the support cable that drift is changed, and identify the knots modification of their drift, this knots modification provides direct basis for the cable force adjustment of the rope, for convenience, the support cable that drift changes is referred to as slack line by this method.
Technical scheme:In the method, coordinate of the bearing on the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing space coordinate ", space coordinate of the bearing on X, Y, Z axis can also be said to be, bearing is referred to as in space coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the space coordinate of some axle with a space coordinate subscale of bearing on the concrete numerical value of the space coordinate of some axle;Angular coordinate of the bearing on X, Y, Z axis is censured with " bearing angular coordinate ", bearing is referred to as in angular coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the angular coordinate of some axle with an angular coordinate subscale of bearing on the concrete numerical value of the angular coordinate of some axle;Censure that bearing angular coordinate and bearing space coordinate are all with " bearing generalized coordinates ", in this method also with a generalized coordinates subscale of bearing up to bearing on the space coordinate of axle or the concrete numerical value of angular coordinate;Bearing is referred to as bearing displacement of the lines on the change of the coordinate of X, Y, Z axis, it may also be said to which the change of bearing space coordinate is referred to as the concrete numerical value in bearing displacement of the lines, this method also with the translational component expression bearing of bearing on the displacement of the lines of some axle;Bearing is referred to as the concrete numerical value in angular displacement of support, this method also with the angular displacement component expression bearing of bearing on the angular displacement of some axle on the change of the angular coordinate of X, Y, Z axis;Generalized displacement of support censures bearing displacement of the lines and angular displacement of support is all, and also bearing is reached on the displacement of the lines of some axle or the concrete numerical value of angular displacement with a generalized displacement subscale of bearing in this method;Bearing displacement of the lines is alternatively referred to as translational displacement, and support settlement is the component of bearing displacement of the lines or translational displacement in gravity direction.
The external force that object, structure are born can be described as load, and load includes face load and volume load.Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load.Volume load is the continuously distributed load in interior of articles each point, the deadweight of such as object and inertia force.
Concentrfated load is divided into two kinds of concentrated force and concentrated couple, in a coordinate system, for example in Descartes's rectangular coordinate system, one concentrated force can resolve into three components, likewise, a concentrated couple can also resolve into three components, if load is actually concentrfated load, force component or a concentrated couple component is concentrated to be referred to as a load by one in the method, the now change of load is embodied as a change for concentrating force component or a concentrated couple component.
Distributed load is divided into line distributed load and EDS maps load, the size of the zone of action and distributed load of the description of distributed load at least including distributed load, the size of distributed load is expressed with distribution intensity, being distributed intensity, to express, (such as two distributed loads are all uniform with distribution characteristics (such as uniform, SIN function equal distribution feature) and amplitude, but its amplitude is different, and the concept of amplitude can be illustrated by taking well-distributed pressure as an example:Same structure bears two different well-distributed pressures, and two distributed loads are all uniform loads, but the amplitude of a distributed load is 10MPa, and the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, one distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is not all identical, the component of this several distributed load is so regarded as same amount of independent distributed load in the method, now a load just represents the component of a distributed load, the amplitude changing ratio identical component of wherein distribution intensity can also be synthesized into a distributed load or is load.
Volume load is the continuously distributed load in interior of articles each point, deadweight and inertia force such as object, the size of the zone of action and volume load of the description of volume load at least including volume load, the size of volume load is expressed with distribution intensity, with distribution characteristics (such as uniform, linear function equal distribution feature) and amplitude, to express, (it is all uniform that lotuses are storaged in such as two individuals to distribution intensity, but its amplitude is different, and the concept of amplitude is illustrated exemplified by can conducting oneself with dignity:The material of two parts of same structure is different, therefore density is different, so while the volume load suffered by the two parts is all uniform, but the amplitude of the volume load suffered by a part is probably 10kN/m3, the amplitude of the volume load suffered by another part is 50kN/m3).If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes.In a coordinate system, one individual stowage lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles on coordinate system, that is, volume load can resolve into three components in Descartes's rectangular coordinate system), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is not all identical, the component that lotus so is storaged in this several body in the method regards same amount of independent load as, the amplitude changing ratio identical volume sharing part of the load of wherein distribution intensity can also be synthesized to an individual stowage lotus or it is load.
When load is embodied as concentrfated load, in the method, " load unit change " actually refers to " unit change of concentrfated load ", similar, " load change " is referred specifically to " change of the size of concentrfated load ", " load change amount " is referred specifically to " variable quantity of the size of concentrfated load ", " load change degree " is referred specifically to " intensity of variation of the size of concentrfated load ", " the actual change amount of load " refers to " the actual change amount of the size of concentrfated load ", " load changed " refers to " concentrfated load that size changes ", briefly, now " so-and-so load so-and-so change " refers to " size of so-and-so concentrfated load so-and-so change ".
When load is embodied as distributed load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, similar, " load change " is referred specifically to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " is referred specifically to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " is referred specifically to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " is referred specifically to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that the amplitude of distribution intensity changes ", briefly, now " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so distributed load so-and-so change ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.
When load is embodied as volume load, in the method, " load unit change " actually refers to " unit change of the amplitude of the distribution intensity of volume load ", similar, " load change " refers to " change of the amplitude of the distribution intensity of volume load ", " load change amount " refers to " variable quantity of the amplitude of the distribution intensity of volume load ", " load change degree " refers to " intensity of variation of the amplitude of the distribution intensity of volume load ", " the actual change amount of load " refers to " the actual change amount of the amplitude of the distribution intensity of volume load ", " load changed " refers to " the volume load that the amplitude of distribution intensity changes ", briefly, " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.
This method is specifically included:
A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable;The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation;This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, i.e. evaluation object is N;The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step;This method represents this numbering, k=1,2,3 ..., N with variable k;This method refers exclusively to evaluated support cable and generalized displacement of support component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of generalized displacement of support component are P, i.e. the quantity of core evaluation object is P, and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object ";The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, the monitored rope force data of the Cable Structure M in Cable Structure1The M of individual specified support cable1Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change;M is had every time1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure;M1It is an integer for being not more than Q not less than 0;Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K2The L of individual specified point and each specified point2The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K2The change of all tested strains of individual specified point;M is had every time2Individual strain measurement value or calculated value strain to characterize Cable Structure, M2For K2And L2Product;M2It is no less than 0 integer;Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure3Individual specified point, excessively each specified point L3The H of individual specified straight line, each specified straight line3Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified;M is had every time3Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M3For K3、L3And H3Product;M3It is an integer not less than 0;The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure4The L of individual specified point and each specified point4The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K4The change of all coordinate components of individual specified point;M is had every time4Individual coordinates measurements or calculated value characterize Cable Structure shape, M4For K4And L4Product;M4It is an integer not less than 0;The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M1、M2、M3And M4Sum, it is M to define parameter K, K1、K2、K3And K4Sum, M should be more than the quantity of core evaluation object, and M is less than the quantity of evaluation object;For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method;Time interval between any measurement twice monitored in real time to same amount in this method cannot be greater than 30 minutes, be referred to as the physical record data moment at the time of measurement record data;The external force that object, structure are born can be described as load, and load includes face load and volume load;Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load;Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force;Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component;Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude;If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant;Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load;Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude;If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes;Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method;
B. this method definition " temperature survey of the Cable Structure of this method calculates method " is carried out by step b1 to b3;
b1:Inquiry or actual measurement obtain the thermal conduction study parameter varied with temperature of Cable Structure composition material and Cable Structure local environment, using the geometry measured data of the design drawing, as-built drawing and Cable Structure of Cable Structure, the Thermodynamic calculation model of Cable Structure is set up using these data and parameter;Inquire about the meteorological data in recent years that Cable Structure location is no less than 2 years, count the cloudy quantity obtained in this period and be designated as T cloudy day, it will can not see the one of the sun daytime in the method and be referred to as the cloudy day all day, obtain each cloudy day in T cloudy day 0 is counted up to the highest temperature and the lowest temperature after sunrise moment next day between 30 minutes, the sunrise moment referred to according to the sunrise moment on earth rotation and the meteorology of revolution rule determination, do not indicate that the same day necessarily can see that the sun, data can be inquired about or calculated by conventional meteorology and obtain the required sunrise moment of each day, the 0 of each cloudy day subtracts the maximum temperature difference that the lowest temperature is referred to as the cloudy daily temperature up to the highest temperature after sunrise moment next day between 30 minutes, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximum in the maximum temperature difference of T cloudy daily temperature is taken to refer to temperature difference per day, Δ T is designated as with reference to temperature difference per dayr;Inquiry Cable Structure location and meteorological data in recent years of the height above sea level interval in place no less than 2 years or actual measurement obtain the temperature of Cable Structure local environment with time and the delta data and changing rule of height above sea level, calculate maximum rate of change Δ T of the temperature for obtaining the Cable Structure local environment in recent years of Cable Structure location and place height above sea level interval no less than 2 years on height above sea levelh, Δ T is taken for convenience of narrationhUnit for DEG C/m;Taken on the surface of Cable Structure " R Cable Structure surface point ", the Specific Principles of " R Cable Structure surface point " are taken to be described in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, it is called " R Cable Structure surface temperature measured data " to survey obtained temperature data, if utilizing the Thermodynamic calculation model of Cable Structure, the temperature of this R Cable Structure surface point is obtained by Calculation of Heat Transfer, it is called " R Cable Structure land surface pyrometer count evidence " just to calculate obtained temperature data;From the minimum height above sea level residing for Cable Structure to highest height above sea level,It is uniform in Cable Structure to choose no less than three different height above sea levels,At the height above sea level of each selection,Two points are at least chosen at the intersection on horizontal plane Yu Cable Structure surface,The exterior normal of straw line body structure surface at selected point,The exterior normal direction of all selections is referred to as in " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness ",Direction of the Cable Structure along the Temperature Distribution of wall thickness is measured with " horizontal plane and the intersection on Cable Structure surface " to intersect,The sunny slope exterior normal direction of Cable Structure and the in the shade face exterior normal direction of Cable Structure must be included in direction of the measurement Cable Structure along the Temperature Distribution of wall thickness of selection,No less than three points are chosen along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness is uniform in Cable Structure,For support cable a point is only taken along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness,Only measure the temperature of the surface point of support cable,All temperature being selected a little of measurement,The temperature measured is referred to as " temperature profile data of the Cable Structure along thickness ",Wherein edge is intersected with same " horizontal plane and the intersection on Cable Structure surface "," temperature profile data of the Cable Structure along thickness " that " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness " measurement is obtained,It is referred to as in the method " temperature profile data of the identical height above sea level Cable Structure along thickness ",If have chosen H different height above sea levels,At each height above sea level,It has chosen direction of the B measurement Cable Structure along the Temperature Distribution of wall thickness,Direction along each measurement Cable Structure along the Temperature Distribution of wall thickness have chosen E point in Cable Structure,Wherein H and E are not less than 3,B is not less than 2,It is equal to 1 for support cable E,The sum for counting in Cable Structure " point of the measurement Cable Structure along the temperature profile data of thickness " is HBE,The temperature of this HBE " point of the measurement Cable Structure along the temperature profile data of thickness " will be obtained by actual measurement below,It is called " HBE Cable Structure along thickness temperature measured data " to survey obtained temperature data,If utilizing the Thermodynamic calculation model of Cable Structure,Temperature of this HBE measurement Cable Structure along the point of the temperature profile data of thickness is obtained by Calculation of Heat Transfer,It is called " HBE Cable Structure calculates data along thickness temperature " just to calculate obtained temperature data;Require to choose a position according to meteorology measurement temperature in Cable Structure location, the temperature of environment where meeting the Cable Structure of meteorology measurement temperature requirement will be obtained in the actual measurement of this position;A position is chosen at the spacious unobstructed place in Cable Structure location, the position should it is annual can obtain each day the ground this getable day most sufficient sunshine, in the flat board of one piece of carbon steel material of position of sound production, referred to as reference plate, reference plate not can contact with ground, distance is not less than 1.5 meters to reference plate from the ground, the one side of the reference plate faces south, referred to as sunny slope, the sunny slope of reference plate is coarse and dark, the sunny slope of reference plate should it is annual can obtain each day one flat plate the ground this getable day most sufficient sunshine, the non-sunny slope of reference plate is covered with insulation material, real-time monitoring is obtained to the temperature of the sunny slope of reference plate;
b2:Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point in real time, monitoring obtains temperature profile data of the previously defined Cable Structure along thickness in real time simultaneously, while monitoring in real time obtains the temperature record of environment where meeting the Cable Structure of meteorology measurement temperature requirement;By monitor in real time obtain the Cable Structure that the same day is carved into after sunrise moment next day between 30 minutes at sunrise where environment temperature measured data sequence, the temperature measured data of environment is arranged according to time order and function order where the temperature measured data sequence of environment was carved into after the sunrise moment next day Cable Structure between 30 minutes at sunrise by the same day where Cable Structure, maximum temperature and minimum temperature in the temperature measured data sequence of environment where finding Cable Structure, the same day for subtracting environment where minimum temperature obtains Cable Structure with the maximum temperature in the temperature measured data sequence of environment where Cable Structure is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as environment maximum temperature difference, it is designated as Δ Temax;Rate of change of the temperature of environment where obtaining Cable Structure on the time is calculated by Conventional mathematical by the temperature measured data sequence of environment where Cable Structure, the rate of change is also with time change;By the measured data sequence for monitoring the temperature for obtaining the sunny slope that the same day is carved into reference plate after sunrise moment next day between 30 minutes at sunrise in real time, the measured data that the measured data sequence of the temperature of the sunny slope of reference plate was carved into the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by the same day at sunrise is arranged according to time order and function order, maximum temperature and minimum temperature in the measured data sequence for the temperature for finding the sunny slope of reference plate, same day of temperature that the sunny slope that minimum temperature obtains reference plate is subtracted with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate is carved into maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as reference plate maximum temperature difference, it is designated as Δ Tpmax;The Cable Structure surface temperature measured data sequence for all R Cable Structure surface points that the same day is carved into after sunrise moment next day between 30 minutes at sunrise is obtained by monitoring in real time, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, the Cable Structure surface temperature measured data that each Cable Structure surface temperature measured data sequence was carved into after sunrise moment next day between 30 minutes by the same day of a Cable Structure surface point at sunrise is arranged according to time order and function order, find the maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, the same day for subtracting the temperature that minimum temperature obtains each Cable Structure surface point with the maximum temperature in each Cable Structure surface temperature measured data sequence is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, there is R Cable Structure surface point just to there is the R same day to be carved into the maximum temperature difference numerical value after sunrise moment next day between 30 minutes at sunrise, maximum therein is referred to as Cable Structure surface maximum temperature difference, it is designated as Δ Tsmax;Calculated by each Cable Structure surface temperature measured data sequence by Conventional mathematical and obtain rate of change of the temperature on the time of each Cable Structure surface point, the temperature of each Cable Structure surface point on the time rate of change also with time change;Obtain the same day by monitoring in real time and be carved at sunrise after sunrise moment next day between 30 minutes, in synchronization, after HBE " temperature profile data of the Cable Structure along thickness ", calculate the difference of the maximum temperature and minimum temperature that amount at the height above sea level of each selection in BE " temperature profile data of the identical height above sea level Cable Structure along thickness ", the absolute value of this difference is referred to as " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", have chosen H different height above sea levels just has H " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", maximum in this H " Cable Structure thickness direction maximum temperature difference at identical height above sea level " is called " Cable Structure thickness direction maximum temperature difference ", it is designated as Δ Ttmax
b3:Survey calculation obtains Cable Structure steady temperature data;First, it is determined that at the time of obtaining Cable Structure steady temperature data, the condition related at the time of Cable Structure steady temperature data to determining to obtain has six, Section 1 condition was carved into after sunrise moment next day between 30 minutes at sunset at the time of being and obtain Cable Structure steady temperature data between the same day, the sunset moment refers to that, according to the sunset moment on earth rotation and the meteorology of revolution rule determination, data can be inquired about or calculate by conventional meteorology obtaining the required sunset moment of each day;The a conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, reference plate maximum temperature difference Δ TpmaxWith Cable Structure surface maximum temperature difference Δ TsmaxAll it is not more than 5 degrees Celsius;The b conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, the environment maximum temperature difference Δ T obtained in above survey calculationemaxNo more than refer to temperature difference per day Δ Tr, and reference plate maximum temperature difference Δ TpmaxSubtract and be not more than Δ T after 2 degrees Celsiusemax, and Cable Structure surface maximum temperature difference Δ TsmaxNo more than Δ Tpmax;Need to only meet in a conditions and b conditions of Section 2 one is known as meeting Section 2 condition;Section 3 condition is that the temperature of environment is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time where Cable Structure at the time of Cable Structure steady temperature data are obtained;Section 4 condition is that at the time of Cable Structure steady temperature data are obtained, the temperature of each Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time;Section 5 condition is that at the time of Cable Structure steady temperature data are obtained, the Cable Structure surface temperature measured data of each Cable Structure surface point in R Cable Structure surface point was carved into the minimum after sunrise moment next day between 30 minutes for the same day at sunrise;Section 6 condition is " Cable Structure thickness direction maximum temperature difference " Δ T at the time of Cable Structure steady temperature data are obtainedtmaxNo more than 1 degree Celsius;This method utilizes above-mentioned six conditions, any one in the following three moment is referred to as " the mathematics moment for obtaining Cable Structure steady temperature data ", the first moment is at the time of meeting Section 1 to the Section 5 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", second of moment is at the time of only meeting the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", the third moment is while at the time of meeting Section 1 to the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ";It it is exactly the mathematics moment for obtaining Cable Structure steady temperature data at the time of acquisition Cable Structure steady temperature data when a moment in the physical record data moment during the mathematics moment for obtaining Cable Structure steady temperature data is exactly this method;If the mathematics moment for obtaining Cable Structure steady temperature data is not any one moment in the physical record data moment in this method, take this method closest to the mathematics moment for obtaining Cable Structure steady temperature data that physical record data at the time of to obtain at the time of Cable Structure steady temperature data;The amount that this method is used in measurement record at the time of obtaining Cable Structure steady temperature data carries out the monitoring analysis of Cable Structure relevant health;The Cable Structure temperature field that this method is approximately considered at the time of obtaining Cable Structure steady temperature data is in the Cable Structure temperature at stable state, i.e. this moment and not changed over time, and this moment is exactly " at the time of the obtaining Cable Structure steady temperature data " of this method;Then, according to Cable Structure heat-transfer character, utilize " the R Cable Structure surface temperature measured data " at the time of obtaining Cable Structure steady temperature data and " HBE Cable Structure along thickness temperature measured data ", utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution for obtaining Cable Structure at the time of Cable Structure steady temperature data are obtained is calculated by conventional heat transfer, now the temperature field of Cable Structure is calculated by stable state, calculating the temperature profile data of obtained Cable Structure at the time of Cable Structure steady temperature data are obtained includes the calculating temperature of R Cable Structure surface point in Cable Structure, the calculating temperature of R Cable Structure surface point is referred to as R Cable Structure steady-state surface temperature and calculates data, also include calculating temperature of the Cable Structure above selected HBE " point of the measurement Cable Structure along the temperature profile data of thickness ", the calculating temperature of HBE " points of the measurement Cable Structure along the temperature profile data of thickness " is referred to as " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data calculates data correspondent equal with R Cable Structure steady-state surface temperature, and " HBE Cable Structure along thickness temperature measured data " with " HBE Cable Structure along thickness temperature calculating data " correspondent equal when, the temperature profile data for calculating obtained Cable Structure at the time of Cable Structure steady temperature data are obtained is referred to as " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " now is referred to as " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure along thickness temperature measured data " is referred to as " HBE Cable Structure along thickness steady temperature measured data ";When " R Cable Structure surface point " is taken on the surface of Cable Structure, the quantity of " R Cable Structure surface point " must is fulfilled for three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the observed temperature linear interpolation of point adjacent with the arbitrfary point on Cable Structure surface during the temperature at any point on Cable Structure surface is by " R Cable Structure surface point " is obtained, the error of the temperature of the arbitrfary point and the actual temperature of the arbitrfary point on Cable Structure surface is not more than 5% on the Cable Structure surface that linear interpolation is obtained;Cable Structure surface includes support cable surface;Second condition is that the point being not less than in 4, and " R Cable Structure surface point " in same height above sea level in the quantity of the point of same height above sea level in " R Cable Structure surface point " is uniform along Cable Structure surface;Maximum Δ h in the absolute value of the difference of the height above sea level of " R Cable Structure surface point " along all Cable Structure surface points adjacent two-by-two of height above sea level is not more than 0.2 DEG C divided by Δ ThObtained numerical value, Δ T is taken for convenience of narrationhUnit for DEG C/m, for convenience of narration take Δ h unit be m;When the definition of " R Cable Structure surface point " along the Cable Structure surface point adjacent two-by-two of height above sea level refers to only consider height above sea level, a Cable Structure surface point is not present in " R Cable Structure surface point ", the height above sea level numerical value of the Cable Structure surface point is between the height above sea level numerical value of adjacent Cable Structure surface point two-by-two;3rd condition is inquiry or obtains Cable Structure location and the interval sunshine rule of place height above sea level by meteorology conventionally calculation, further according to the geometric properties and bearing data of Cable Structure, found in Cable Structure it is annual by the sunshine-duration most sufficient position of those surface points, in " R Cable Structure surface point " at least one Cable Structure surface point be in Cable Structure whole year by a point in those most sufficient surface points of sunshine-duration;
C. the Cable Structure steady temperature data obtained under original state are calculated according to " temperature survey of the Cable Structure of this method calculates method " direct measurement, Cable Structure steady temperature data under original state are referred to as initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector To”;Survey or consult reference materials and obtain the physical and mechanical properties parameter varied with temperature of various materials used in Cable Structure;Initial Cable Structure steady temperature data vector T is obtained in actual measurementoSynchronization, direct measurement, which is calculated, obtains the Initial cable forces of all support cables, composition Initial cable force vector Fo;According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector ToUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cableoElement coding rule it is identical;T is obtained in actual measurementoWhile, that is, obtaining initial Cable Structure steady temperature data vector ToAt the time of synchronization, direct measurement calculates the measured data for obtaining initial Cable Structure, the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable now are referred to as support cable initial health data;The monitored amount initial value vector C of initial value composition of all monitored amountso, it is monitored amount initial value vector CoCoding rule and M monitored amounts coding rules it is identical;Evaluation object initial damage vector d is set up using support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement datao, vectorial doRepresent with initial mechanical calculating benchmark model AoThe initial health of the evaluation object of the Cable Structure of expression;Evaluation object initial damage vector doElement number be equal to N, doElement and evaluation object be one-to-one relationship, vectorial doElement coding rule it is identical with the coding rule of evaluation object;If doThe corresponding evaluation object of some element be a support cable in cable system, then doThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion;;If doThe corresponding evaluation object of some element be some bearing some generalized displacement component, then doThe element numerical value represent this bearing the generalized displacement component initial value;If doThe corresponding evaluation object of some element be to take d in some load, this methodoThe element numerical value be 0, the initial value for representing the change of this load is 0;If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial doIn each element numerical value related to Cable Structure generalized displacement of support take 0;If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial doIn each element numerical value related to support cable take 0;Initial Cable Structure bearing generalized coordinates data refer to the bearing generalized coordinates data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model AoWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs;
D. the physical and mechanical properties parameter varied with temperature of various materials, initial Cable Structure steady temperature data vector T according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable StructureoAll Cable Structure data obtained with preceding step, set up the initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data "o, based on AoCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot be greater than 5%;Corresponding to Ao" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector To”;Corresponding to AoEvaluation object health status evaluation object initial damage vector doRepresent;Corresponding to AoAll monitored amounts initial value with monitored amount initial value vector CoRepresent;The current initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data " is set up for the first timet o, the monitored current initial value vector C of amountt o" current initial Cable Structure steady temperature data vector Tt o”;The current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timet oWith monitored amount current initial value vector Ct oWhen, the current initial mechanical calculating benchmark model A of Cable Structuret oIt is equal to the initial mechanical calculating benchmark model A of Cable Structureo, it is monitored the current initial value vector C of amountt oIt is equal to monitored amount initial value vector Co;At oCorresponding " Cable Structure steady temperature data " are referred to as " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector Tt o", the current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timet oWhen, Tt oIt is equal to To;At oEvaluation object initial health and AoEvaluation object health status it is identical, also with evaluation object initial damage vector doRepresent, the A in cyclic process belowt oEvaluation object initial health all the time use evaluation object initial damage vector doRepresent;ToAnd doIt is AoParameter, by AoThe obtained initial value of all monitored amounts of Mechanics Calculation result and CoThe initial value of all monitored amounts represented is identical, therefore alternatively CoBy AoMechanics Calculation result composition;Tt oAnd doIt is At oParameter, Ct oBy At oMechanics Calculation result composition;
E. the circulation walked by e to the n-th step is entered from here;During structure military service, constantly the current data of " Cable Structure steady temperature data " is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data of " Cable Structure steady temperature data " is referred to as " current Cable Structure steady temperature data ", is designated as " current Cable Structure steady temperature data vector Tt", vector TtDefinition mode and vector ToDefinition mode it is identical;Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, actual measurement obtain all M in Cable Structure1The rope force data of root support cable, all these rope force data composition current cable force vector F, vectorial F element and vector FoElement coding rule it is identical;Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, Actual measurement obtains all M1The space coordinate of two supporting end points of root support cable, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vectoroElement coding rule it is identical;
F. according to current Cable Structure steady temperature data vector Tt, current initial mechanical calculating benchmark model A is updated according to step f1 to f3t o, the monitored current initial value vector C of amountt oWith current initial Cable Structure steady temperature data vector Tt o
F1. T is comparedtWith Tt oIf, TtEqual to Tt o, then At o、Ct oAnd Tt oKeep constant;Otherwise need to follow these steps to At o、Ct oAnd Tt oIt is updated;
F2. T is calculatedtWith ToDifference, TtWith ToDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TtWith ToDifference represented with steady temperature change vector S, S be equal to TtSubtract To, S represents the change of Cable Structure steady temperature data;
F3. to AoIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, to AoIn the temperature change that applies of Cable Structure after the current initial mechanical calculating benchmark model A that is updatedt o, update At oWhile, Tt oAll elements numerical value also uses TtAll elements numerical value correspondence replace, that is, have updated Tt o, thus obtained properly corresponding to At oTt o;Update Ct oMethod be:As renewal At oAfterwards, A is obtained by Mechanics Calculationt oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Ct o;At oSupport cable initial health all the time use evaluation object initial damage vector doRepresent;
G. in current initial mechanical calculating benchmark model At oOn the basis of carry out Mechanics Calculation several times according to step g1 to g4, pass through to calculate and obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector Du
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is to constantly update, that is, is updating current initial mechanical calculating benchmark model At o, the monitored current initial value vector C of amountt oWith current initial Cable Structure steady temperature data vector Tt oAfterwards, it is necessary to then update Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector Du
G2. in the current initial mechanical calculating benchmark model A of Cable Structuret oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating;According to the coding rule of evaluation object, calculated successively;Calculate each time and assume that only one of which evaluation object is further added by unit damage or unit generalized displacement or load unit change on the basis of original damage or generalized displacement or load, specifically, if the evaluation object is a support cable in cable system, then it is assumed that the support cable is in vectorial doThe support cable represented is further added by unit damage on the basis of having damaged, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing is further added by unit generalized displacement in the direction of displacement, if the evaluation object is a load, it is assumed that the load is in vectorial doLoad unit change is further added by the basis of the existing variable quantity of the load represented, D is usedukThis increased unit damage or unit generalized displacement or load unit change are recorded, wherein k represents to increase the numbering of the evaluation object of unit damage or unit generalized displacement or load unit change, DukIt is evaluation object unit change vector DuAn element, evaluation object unit change vector DuElement coding rule and vector doElement coding rule it is identical;The evaluation object of increase unit damage or unit generalized displacement or load unit change during the evaluation object of increase unit damage or unit generalized displacement or load unit change is calculated different from other times in calculating each time, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, the current calculated value for calculating obtained all monitored amounts each time constitutes a monitored amount calculation current vector, is monitored the element number rule and monitored amount initial value vector C of amount calculation current vectoroElement number rule it is identical;
G3. obtained monitored amount calculation current vector is calculated each time subtracts the monitored current initial value vector C of amountt oObtain a vector, each element of the vector divided by the assumed unit damage of this time calculating or unit generalized displacement or load unit are changed into numerical value again, a monitored amount unit change vector is obtained, has N number of evaluation object just to have N number of monitored amount unit change vector;
G4. by coding rule of this N number of monitored amount unit change vector according to N number of evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is arranged successively;Cable Structure unit damage monitored numerical quantity transformation matrices Δ C each row correspond to a monitored amount unit change vector;Cable Structure unit damage monitored numerical quantity transformation matrices Δ C every a line corresponds to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit generalized displacement or load unit change;The coding rule of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C row and vector doElement coding rule it is identical, the coding rule of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C row is identical with the coding rules of M monitored amounts;
H. current Cable Structure steady temperature data vector T is obtained in actual measurementtWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TtAt the time of synchronization Cable Structure all monitored amounts current measured value, the monitored amount current value vector C of composition;The monitored amount current value vector C and monitored current initial value vector C of amountt oWith monitored amount initial value vector CoDefinition mode it is identical, the same monitored amount of element representation of three vectorial identical numberings is in concrete numerical value not in the same time;
I. the current nominal fatigue vector d of evaluation object is defined, the current nominal fatigue vector d of evaluation object element number is equal to the quantity of evaluation object, it is one-to-one relationship between the current nominal fatigue vector d of evaluation object element and evaluation object, the current nominal fatigue vector d of evaluation object element numerical value represents the nominal fatigue degree or nominal generalized displacement or nominal load variable quantity of correspondence evaluation object;The coding rule of vectorial d element and vector doElement coding rule it is identical;
J. according to monitored amount current value vector C with the monitored current initial value vector C of amountt o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and the current nominal fatigue vector d of evaluation object to be asked, the linear approximate relationship can be expressed as formula 1, other amounts in formula 1 in addition to d are, it is known that solution formula 1 can just calculate the current nominal fatigue vector d of evaluation object;
Formula 1
K. the currently practical injury vector d of evaluation object is defineda, the currently practical injury vector d of evaluation objectaElement number be equal to evaluation object quantity, the currently practical injury vector d of evaluation objectaElement and evaluation object between be one-to-one relationship, the currently practical injury vector d of evaluation objectaElement numerical value represent correspondence evaluation object actual damage degree or actual generalized displacement or real load variable quantity;Vectorial daElement coding rule and vector doElement coding rule it is identical;
L. the currently practical injury vector d of evaluation object expressed using formula 2aK-th of element da kWith evaluation object initial damage vector doK-th of element dokWith evaluation object current nominal fatigue vector d k-th of element dkBetween relation, calculating obtain the currently practical injury vector d of evaluation objectaAll elements;
K=1 in formula 2,2,3 ... ..., N, da kRepresent the currently practical health status of k-th of evaluation object, da kFor 0 when represent k-th of evaluation object without health problem, da kNumerical value represents that k-th of evaluation object is the evaluation object of unsoundness problem when not being 0, if the evaluation object is a support cable in cable system, then da kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, da kThe numerical response degree of relaxation or the damage of the support cable;Damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line, the currently practical injury vector d of evaluation objectaIn correspond to slack line element numerical expression be with slack line relax degree mechanic equivalent currently practical equivalent damage degree;If the evaluation object is a generalized displacement component of a bearing, then da kRepresent its currently practical generalized displacement numerical value;
M. utilize in current Cable Structure steady temperature data vector TtUnder the conditions of, in l walk the slack line that identifies and with the currently practical injury vector d of evaluation objectaExpression these slack lines, the degree that relaxed with it mechanic equivalent currently practical equivalent damage degree, using e walk acquisition in current Cable Structure steady temperature data vector TtUnder the conditions of current cable force vector F and the supporting end points horizontal range vector of current support cable two, using c walk obtain in initial Cable Structure steady temperature data vector ToUnder the conditions of support cable initial drift is vectorial, the weight of the initial free unit length of initial free cross-sectional area vector sum is vectorial, Initial cable force vector Fo, utilize current Cable Structure steady temperature data vector TtThe support cable current steady state temperature data of expression, using c walk obtain in initial Cable Structure steady temperature data vector ToThe support cable initial steady state temperature data of expression, utilize the physical and mechanical properties parameter varied with temperature that various materials used in the Cable Structure of acquisition are walked in c, it is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is:First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical;2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope;When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa;Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined;So it is achieved that relaxation identification and the non-destructive tests of support cable;Institute's demand power is provided by current cable force vector F corresponding elements during calculating;Damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method, so, this method is according to the currently practical injury vector d of evaluation objectaProblem cable can either be identified, also can determine which bearing there occurs generalized displacement and its numerical value;
N. e steps are returned to, start the circulation next time walked by e to the n-th step.
Beneficial effect:Structural healthy monitoring system is monitored on-line for a long time by using sensor to structural response first, (or offline) analysis online is carried out after acquisition Monitoring Data to it and obtains structural health conditions data, due to the complexity of structure, structural healthy monitoring system needs to use the equipment such as substantial amounts of sensor to carry out monitoring structural health conditions, therefore the usual suitable height of its cost, it may be said that cost problem is also a subject matter of limit structural health monitoring technique application.On the other hand, the correct identification of the health status of core evaluation object (such as suspension cable) is the indispensable part of the correct identification of structural health conditions, even its whole, and influence of the correct identification of the change of secondary evaluation object (load that such as structure is born) (such as by the change of the quality and quantity of the automobile of cable-stayed bridge) to the correct identification of the health status of Cable Structure be it is very little, it is even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is typically suitable, the quantity of secondary evaluation object is also frequently more than the quantity of core evaluation object, and the quantity of such evaluation object is often many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, the quantity of the monitored amount of conventional method requirement (being obtained using device measurings such as sensors) have to be larger than the quantity equal to evaluation object, when the secondary evaluation object changed quantity than it is larger when (practically always such), the quantity of the equipment such as the sensor required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, or even unacceptablely high.Fortunately inventor's research is found,In secondary evaluation object (such as normal load that structure is born,The normal load of structure refers to that the load that structure bearing is no more than the structure allowable load limited according to structure design book or structure completion book) change it is smaller when (be exactly that structure is only subjected only to normal load for load,Whether the load that structure is born is normal load,It can be determined by the observation of the methods such as naked eyes,If it find that the load that structure is born is not normal load,It is so artificial to remove,Remove after improper load,Structure is just solely subjected to normal load),Amplitude of variation of the amplitude of variation (this specification is called " secondary response ") of structural response caused by them much smaller than the structural response caused by the change (such as support cable is damaged) of core evaluation object (this specification is called " core response "),Secondary response responds total change that sum is structural response with core (this specification is called " global response "),Obvious core response occupies leading position in global response,Even if being chosen when now determining that monitored amount quantity based on this inventor research hair and being slightly larger than core evaluation object quantity,But much smaller than the numerical value of evaluation object quantity (this method is exactly so to do),Even if that is using equipment such as the relatively few many sensors of quantity,Still the state of health data of core evaluation object can accurately be obtained,Meet the core demand of structural health conditions monitoring,Therefore cost of the cost of the structural healthy monitoring system proposed by this method apparently than the structural healthy monitoring system required by conventional method is much lower,That is this method can realize the assessment of the health status to the core evaluation object of Cable Structure with the much lower condition of cost,This benefit is that can structural health monitoring technology be used is very important.
Embodiment
This method uses a kind of algorithm, and the algorithm is used to recognize generalized displacement of support and problem cable.When it is implemented, the following steps are one kind in the various steps that can be taken.
The first step:First confirm that the quantity for the load that the possibility that Cable Structure is born changes.The characteristics of load born according to Cable Structure, confirmation wherein " is possible to the load changed ", or all load is considered as " being possible to the load changed ", if the shared JZW load that may be changed, that is, have JZW secondary evaluation objects.
If evaluated support cable and the quantity sum of generalized displacement of support component are that P, the i.e. quantity of core evaluation object are P, if the quantity of evaluated generalized displacement of support component is Z.
If the quantity of the generalized displacement of support component of Cable Structure, the quantity of the support cable of Cable Structure and JZW " being possible to the load changed " quantity sum are N, that is, have N number of evaluation object.Evaluation object serial number is given, the numbering will be used to generate vector sum matrix in subsequent step.
Monitored multiclass parameter can include:Suo Li, strain, angle and space coordinate, are described below respectively:
If having Q root support cables in cable system, the monitored rope force data of the Cable Structure M in Cable Structure1The M of individual specified rope1Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified ropes change.M is had every time1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M1It is an integer not less than 0.
The monitored strain data of Cable Structure can in Cable Structure K2The L of individual specified point and each specified point2The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K2The change of all tested strains of individual specified point.M is had every time2(M2=K2×L2) individual strain measurement value or calculated value characterize Cable Structure strain.M2It is an integer not less than 0.
The monitored angle-data of Cable Structure K in Cable Structure3Individual specified point, excessively each specified point L3The H of individual specified straight line, each specified straight line3Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified.M is had every time3(M3=K3×L3×H3) individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure.M3It is an integer not less than 0.
The monitored shape data of Cable Structure K in Cable Structure4The L of individual specified point and each specified point4The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K4The change of all coordinate components of individual specified point.M is had every time4(M4=K4×L4) individual coordinates measurements or calculated value characterize Cable Structure shape.M4It is an integer not less than 0.
The monitored amount of summary, whole Cable Structure has M (M=M1+M2+M3+M4) individual monitored amount, define parameter K (K=M1+K2+K3+K4), M should add 4 not less than the quantity of core evaluation object, and M is less than the quantity of evaluation object.
For convenience, " monitored all parameters of Cable Structure " are referred to as in the method " monitored amount ".To M monitored amount serial numbers, the numbering will be used to generate vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 ..., M with variable j.
Provide that step is determined " temperature survey of the Cable Structure of this method calculates method " by technical scheme.
Second step:Set up initial mechanical calculating benchmark model Ao
When Cable Structure is completed, or before health monitoring systems are set up, obtain " Cable Structure steady temperature data " according to " temperature survey of the Cable Structure of this method calculates method " survey calculation (can be measured with ordinary temperature measuring method, for example measured using thermal resistance), " Cable Structure steady temperature data " now use vector ToRepresent, referred to as initial Cable Structure steady temperature data vector To.T is obtained in actual measurementoWhile, that is, synchronization at the time of initial Cable Structure steady temperature data vector is obtained, the initial value for all monitored amounts for obtaining Cable Structure, the monitored amount initial value vector C of composition are calculated using conventional method direct measuremento
Synchronization that can specifically in following manner at the time of so-and-so (such as initial or current) Cable Structure steady temperature data vector is obtained in this method, the data of so-and-so the measured monitored amount (all monitored amounts of such as Cable Structure) of amount are obtained using so-and-so method survey calculation:While measurement record temperature (temperature, the temperature of the sunny slope of reference plate and the Cable Structure surface temperature that include environment where Cable Structure), for example every 10 minutes measurement temperature of record, then while equally also recording the data of so-and-so the measured monitored amount (such as all monitored amounts of Cable Structure) of amount every measurement in 10 minutes.Once it is determined that at the time of acquisition Cable Structure steady temperature data, the data of so-and-so the measured monitored amount (all monitored amounts of such as Cable Structure) of amount of synchronization are known as synchronization at the time of Cable Structure steady temperature data are obtained at the time of so with obtaining Cable Structure steady temperature data, the data of so-and-so the measured monitored amount of amount obtained using so-and-so method survey calculation method.
The physical parameter (such as thermal coefficient of expansion) varied with temperature and mechanical property parameters (such as modulus of elasticity, Poisson's ratio) of various materials used in Cable Structure are obtained using conventional method (consult reference materials or survey).
Initial Cable Structure steady temperature data vector T is obtained in actual measurementoSynchronization, direct measurement, which is calculated, obtains the Initial cable forces of all support cables, composition Initial cable force vector Fo;Obtain length of all support cables when free state i.e. Suo Li is 0 according to Cable Structure design data, completion data, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector ToUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, the initial drift vector l of support cable is constituted successivelyo, initial free cross-sectional area vector AoWith the weight vector ω of initial free unit lengtho, the initial drift vector l of support cableo, initial free cross-sectional area vector AoWith the weight vector ω of initial free unit lengthoElement coding rule and Initial cable force vector FoElement coding rule it is identical.
The step of being provided by technical scheme, initial Cable Structure steady temperature data vector T is obtained in Actual measurementoWhile, that is, synchronization at the time of Cable Structure steady temperature data are obtained, the Actual measurement data of Cable Structure are obtained using conventional method Actual measurement.Utilize the physical and mechanical properties parameter varied with temperature and initial Cable Structure steady temperature data vector T of various materials used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structureo, it is included in " Cable Structure steady temperature data " using mechanics method (such as FInite Element) and sets up initial mechanical calculating benchmark model Ao。ToAnd doIt is AoParameter, CoBy AoMechanics Calculation result composition.
3rd step:Current initial mechanical calculating benchmark model A is set up for the first timet o, the monitored current initial value vector C of amountt o" current initial Cable Structure steady temperature data vector Tt o", specific method is:In initial time, i.e., current initial mechanical calculating benchmark model A is set up for the first timet oWith monitored amount current initial value vector Ct oWhen, At oIt is equal to Ao, Ct oIt is equal to Co, At oCorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector Tt o", in initial time (A is namely set up for the first timet oWhen), Tt oIt is equal to To, vector Tt oDefinition mode and vector ToDefinition mode it is identical.At oAssessment object health status and AoAssessment object health status (evaluation object initial damage vector doRepresent) identical, the A in cyclic processt oAssessment object health status all the time use evaluation object initial damage vector doRepresent.Tt oAnd doIt is At oParameter, Ct oBy At oMechanics Calculation result composition.
4th step:During Cable Structure military service, the current data for obtaining " Cable Structure steady temperature data " according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement (is referred to as " current Cable Structure steady temperature data vector Tt", vector TtDefinition mode and vector ToDefinition mode it is identical).Current Cable Structure steady temperature data vector T is obtained in actual measurementtWhile, that is, obtaining current Cable Structure steady temperature data vector TtAt the time of synchronization, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, composition " monitored amount current value vector C ".
Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, actual measurement obtain all M in Cable Structure1The rope force data of root support cable, all these rope force data composition current cable force vector F, vectorial F element and vector FoElement coding rule it is identical;Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, Actual measurement obtains all M1The space coordinate of two of root support cable supporting end points, the difference of the space coordinates of two supporting end points component in the horizontal direction is exactly two supporting end points horizontal ranges, all M1Two supporting end points horizontal range data of root support cable constitute the current supporting end points horizontal range vector of support cable two lt x, the current supporting end points horizontal range vector of support cable two lt xElement coding rule and Initial cable force vector FoElement coding rule it is identical.
5th step:According to current Cable Structure steady temperature data vector Tt, current initial mechanical calculating benchmark model A is updated when necessaryt o, the monitored current initial value vector C of amountt oWith current initial Cable Structure steady temperature data vector Tt o.Current Cable Structure steady temperature data vector T is obtained in the actual measurement of the 4th steptAfterwards, T is comparedtAnd Tt oIf, TtEqual to Tt o, then need not be to At o、Ct oAnd Tt oIt is updated, otherwise needs to At o、Ct oAnd Tt oIt is updated, the step of update method is provided by technical scheme is carried out.
6th step:The step of being provided by technical scheme, in current initial mechanical calculating benchmark model At oOn the basis of carry out Mechanics Calculation several times, pass through calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector Du.If specifically, the evaluation object is a support cable in cable system, then it is assumed that the support cable is in vectorial doThe support cable represented is further added by unit damage (for example taking 5%, 10%, 20% or 30% equivalent damage to be damaged for unit) on the basis of having damaged, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing in the direction of displacement in vectorial doUnit generalized displacement occurs again (if such as evaluation object is the translational component in the x directions of a bearing on the basis of the existing generalized displacement of the bearing represented, it is assumed that the bearing has unit displacement of the lines in x directions, for example take 1mm, if the evaluation object is the angular displacement component around x-axis of a bearing, it is assumed that the bearing has unit angular displacement around x-axis, for example take ten a ten thousandth radians), if the evaluation object is a load, it is assumed that the load is in vectorial doLoad unit is further added by the basis of the existing variable quantity of the load represented to change (if the load is distributed load, and the distributed load is line distributed load, load unit change can take 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. to change for unit;If the load is distributed load, and the distributed load is EDS maps load, and load unit change can take 1MPa, 2MPa, 3MPa or 1kNm/m2、2kNm/m2、3kNm/m2Change Deng for unit;If the load is concentrfated load, and the concentrfated load is couple, and load unit change can take 1kNm, 2kNm, 3kNm etc. to change for unit;If the load is concentrfated load, and the concentrfated load is concentrated force, and load unit change can take 1kN, 2kN, 3kN etc. to change for unit;If the load is volume load, load unit change can take 1kN/m3、2kN/m3、3kN/m3Change Deng for unit).
7th step:Set up linear relationship error vector e and vector g.Utilize data (the monitored current initial value vector C of amount abovet o, unit damage monitored numerical quantity transformation matrices Δ C), while the 6th step is calculated each time, i.e., each time calculate assume evaluation object in only one of which evaluation object increase unit damage unit generalized displacement or load unit change Duk, the evaluation object of increase unit damage or unit generalized displacement or load unit change during the evaluation object of increase unit damage or unit generalized displacement or load unit change is calculated different from other times in calculating each time, the current value that all monitored amounts in Cable Structure are all calculated using mechanics method (for example with FInite Element) is calculated each time, while calculating composition one is monitored amount calculation current vector C each time, one injury vector d of composition is calculated each time, originally existing injury vector d is walked out of only to use in this step, the numerical value of only one of which element takes D in injury vector d all elementsuk, the numerical value of other elements takes 0, coding rule and the vector d of injury vector d elementoElement coding rule it is identical;By C, Ct o、ΔC、Du, d bring formula (1) into, obtain a linear relationship error vector e, each time calculate obtain a linear relationship error vector e;There is N number of evaluation object just to have n times calculating, just there is N number of linear relationship error vector e, a vector is obtained after this N number of linear relationship error vector e is added, is exactly final linear relationship error vector e by the new vector obtained after each element divided by N of this vector.Vectorial g is equal to final error vector e.
Abs () is the function that takes absolute value in formula (1), and each the vectorial element tried to achieve in bracket is taken absolute value.
8th step:The hardware components of cable structure health monitoring system are installed.Hardware components at least include:Monitored amount monitoring system (such as subsystem containing angular surveying, cable force measurement subsystem, strain measurement subsystem, space coordinate measures subsystem, signal conditioner etc.), Cable Structure temperature monitoring system (contains temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (contain temperature sensor, signal conditioner etc.), support cable cable force monitoring system, the space coordinate monitoring system of the supporting end points of support cable, signal (data) collector, computer and communication alert equipment.Each monitored amount, each temperature, the Suo Li of each support cable, the space coordinate of the supporting end points of each support cable must be monitored system monitoring and arrive, and the signal monitored is transferred to signal (data) collector by monitoring system;Signal is delivered to computer through signal picker;Computer is then responsible for the health monitoring software of the evaluation object of operation Cable Structure, including the signal that the transmission of tracer signal collector comes;When monitoring that evaluation object health status is changed, computer control communication warning device is alarmed to monitoring personnel, owner and (or) the personnel specified.
9th step:By the current initial value vector C of monitored amountt o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector DuParameter is stored in the way of data file on the hard disc of computer of operation health monitoring systems software.
Tenth step:Establishment and on computers installation and operation this method system software, the software will complete the functions (all work that can be completed with computer i.e. in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, alarm that this method required by task is wanted
11st step:According to monitored amount current value vector C with the monitored current initial value vector C of amountt o, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector DuThe linear approximate relationship (formula (2)) existed between the current nominal fatigue vector d (being made up of all Suo Dangqian nominal fatigues amounts) of evaluation object, calculate evaluation object current nominal fatigue vector d noninferior solution according to multi-objective optimization algorithm, that is, with reasonable error but can relatively accurately reflect evaluation object health status change solution.The Objective Programming (GoalAttainmentMethod) in multi-objective optimization algorithm can be used to solve formula (2) and obtain current injury vector d.
12nd step:Utilize the currently practical injury vector d of evaluation objectaK-th of element da kWith evaluation object initial damage vector doK-th of element dokWith evaluation object current nominal fatigue vector d k-th of element dkBetween relation, calculating obtain the currently practical injury vector d of evaluation objectaAll elements.da kThe currently practical health status of k-th of evaluation object is represented, if the evaluation object is a support cable in cable system, then da kRepresent its currently practical damage, da kFor 0 when represent its corresponding support cable without health problem, da kNumerical value represents that its corresponding support cable is the support cable of unsoundness problem when not being 0, and the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, the degree of the relaxation of its numerical response or damage;If the evaluation object is a generalized displacement component of a bearing, then da kRepresent its currently practical generalized displacement numerical value.
13rd step:By the currently practical injury vector d of evaluation objectaIn the Q element related to support cable take out, the currently practical injury vector d of composition support cableca, the currently practical injury vector d of support cablecaElement coding rule and Initial cable force vector FoElement coding rule it is identical.The currently practical injury vector d of support cablecaH-th of element representation Cable Structure in h root support cables currently practical amount of damage, h=1,2,3 ... ..., Q;The currently practical injury vector d of support cablecaMiddle numerical value does not correspond to the support cable of unsoundness problem for 0 element, damaged cable is identified from the support cable of these unsoundness problems by lossless detection method, remaining is exactly slack line, that is, needs to adjust Suo Li rope, and these need to adjust Suo Li rope in the currently practical injury vector d of support cablecaIn corresponding element numerical value (such as one of element can use dca hRepresent) represent degree of injury with the relaxation degree mechanic equivalents of these support cables.Damaged cable is in the currently practical injury vector d of support cablecaIn the numerical value of corresponding element mean that its degree of injury, the numerical value of corresponding element represents that the support cable thoroughly loses bearing capacity when being 100%, represent that the support cable loses the bearing capacity of corresponding proportion when between 0 and 100%, so far just have identified damaged cable and its degree of injury.
14th step:It is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent come calculate slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, specifically can according to formula (3) can in the hope of these ropes relaxation degree (i.e. the long adjustment amount of rope).So it is achieved that the relaxation identification of support cable.So far damaged cable and slack line all just are identified.
E in formula (3)t hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structuret oDuring expression, the modulus of elasticity of h root support cables, At hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structuret oDuring expression, the cross-sectional area of h root support cables, FhIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structuret oDuring expression, the current cable power of h root support cables, dca hIt is the currently practical degree of injury of h root support cables, ωt hIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structuret oDuring expression, the weight of the unit length of h root support cables, lt xhIt is the steady temperature data currently initial Cable Structure steady temperature data vector T in Cable Structuret oDuring expression, the horizontal range of two supporting end points of h root support cables, lt xhIt is the current supporting end points horizontal range vector of support cable two lt xAn element, current support cable two supporting end points horizontal range vector lt xElement coding rule and initial drift vector loElement coding rule it is identical, Et hIt can be obtained according to the characteristic material data for looking into or surveying h root support cables, At hAnd ωt hThermal coefficient of expansion, A that can be according to h root support cablesoh、ωoh、Fh、ToAnd Tt oObtained by Typical physical and Mechanics Calculation.
15th step:Computer in health monitoring systems is periodically automatic or generates cable system health condition form by human users' health monitoring systems.
16th step:Under specified requirements, the computer in health monitoring systems is automatically brought into operation communication alert equipment and alarmed to monitoring personnel, owner and (or) the personnel specified.
17th step:The 4th step is returned to, starts the circulation by the 4th step to the 17th step.

Claims (1)

1. simplify hybrid monitoring problem cable load generalized displacement recognition methods, it is characterised in that methods described includes:
This method is specifically included:
A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable;The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation;This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, i.e. evaluation object is N;The coding rule of evaluation object is determined, is numbered evaluation object all in Cable Structure by this rule, the numbering will be used to generate vector sum matrix in subsequent step;This method represents this numbering, k=1,2,3 ..., N with variable k;This method refers exclusively to evaluated support cable and generalized displacement of support component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of generalized displacement of support component are P, i.e. the quantity of core evaluation object is P, and this method refers exclusively to the evaluated load in " evaluation object " with title " secondary evaluation object ";The support cable by monitored Suo Li specified during hybrid monitoring is determined, if having Q root support cables in cable system, the monitored rope force data of the Cable Structure M in Cable Structure1The M of individual specified support cable1Individual rope force data is described, and Cable Structure Suo Li change is exactly the Suo Li of all specified support cables change;M is had every time1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure;M1It is an integer for being not more than Q not less than 0;Determine the measured point by monitored strain specified during hybrid monitoring, the monitored strain data of Cable Structure can in Cable Structure K2The L of individual specified point and each specified point2The strain of individual assigned direction is described, and the change of Cable Structure strain data is exactly K2The change of all tested strains of individual specified point;M is had every time2Individual strain measurement value or calculated value strain to characterize Cable Structure, M2For K2And L2Product;M2It is no less than 0 integer;Determine the measured point by monitored angle specified during hybrid monitoring, the monitored angle-data of the Cable Structure K in Cable Structure3Individual specified point, excessively each specified point L3The H of individual specified straight line, each specified straight line3Individual angle coordinate component is described, and the change of Cable Structure angle is exactly change of all specified points, all specified straight lines, all angle coordinate components specified;M is had every time3Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M3For K3、L3And H3Product;M3It is an integer not less than 0;The shape data that will be monitored for determining to specify during hybrid monitoring, the monitored shape data of the Cable Structure K in Cable Structure4The L of individual specified point and each specified point4The space coordinate of individual assigned direction is described, and the change of Cable Structure shape data is exactly K4The change of all coordinate components of individual specified point;M is had every time4Individual coordinates measurements or calculated value characterize Cable Structure shape, M4For K4And L4Product;M4It is an integer not less than 0;The monitored amount of summary hybrid monitoring, whole Cable Structure has M monitored amounts, and M is M1、M2、M3And M4Sum, it is M to define parameter K, K1、K2、K3And K4Sum, M should be more than the quantity of core evaluation object, and M is less than the quantity of evaluation object;For convenience, listed M monitored amount of this step is referred to as " monitored amount " in the method;Time interval between any measurement twice monitored in real time to same amount in this method cannot be greater than 30 minutes, be referred to as the physical record data moment at the time of measurement record data;The external force that object, structure are born can be described as load, and load includes face load and volume load;Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load;Volume load is the continuously distributed load in interior of articles each point, including the deadweight of object and inertia force;Concentrfated load is divided into two kinds of concentrated force and concentrated couple, including in the coordinate system including Descartes's rectangular coordinate system, one concentrated force can resolve into three components, same, one concentrated couple can also resolve into three components, if load is actually concentrfated load, it is a load to concentrate force component or a concentrated couple component to be calculated as or count by one in the method, and the now change of load is embodied as a change for concentrating force component or a concentrated couple component;Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least includes the zone of action of distributed load and the size of distributed load, and the size of distributed load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude;If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant;Including in the coordinate system including Descartes's rectangular coordinate system, one distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three components of this distributed load changes, and the ratio of change is not all identical, it is three distributed loads that so three components of this distributed load, which are calculated as or counted, in the method, and now a load just represents the one-component of distributed load;Volume load is the continuously distributed load in interior of articles each point, and the description of volume load at least includes the zone of action of volume load and the size of volume load, and the size of volume load is expressed with distribution intensity, and distribution intensity is expressed with distribution characteristics and amplitude;If load is actually volume load, in the method actual treatment be volume load distribution intensity amplitude change, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, actually refer to the change of the amplitude of the distribution intensity of volume load during the change for now mentioning load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes;Including in the coordinate system including Descartes's rectangular coordinate system, one individual stowage lotus can resolve into three components, if the amplitude of the respective distribution intensity of three components of this volume load changes, and the ratio of change is not all identical, then three components of this volume load are calculated as or counted as three distributed loads in the method;
B. this method definition " temperature survey of the Cable Structure of this method calculates method " is carried out by step b1 to b3;
b1:Inquiry or actual measurement obtain the thermal conduction study parameter varied with temperature of Cable Structure composition material and Cable Structure local environment, using the geometry measured data of the design drawing, as-built drawing and Cable Structure of Cable Structure, the Thermodynamic calculation model of Cable Structure is set up using these data and parameter;Inquire about the meteorological data in recent years that Cable Structure location is no less than 2 years, count the cloudy quantity obtained in this period and be designated as T cloudy day, it will can not see the one of the sun daytime in the method and be referred to as the cloudy day all day, obtain each cloudy day in T cloudy day 0 is counted up to the highest temperature and the lowest temperature after sunrise moment next day between 30 minutes, the sunrise moment referred to according to the sunrise moment on earth rotation and the meteorology of revolution rule determination, do not indicate that the same day necessarily can see that the sun, data can be inquired about or calculated by conventional meteorology and obtain the required sunrise moment of each day, the 0 of each cloudy day subtracts the maximum temperature difference that the lowest temperature is referred to as the cloudy daily temperature up to the highest temperature after sunrise moment next day between 30 minutes, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximum in the maximum temperature difference of T cloudy daily temperature is taken to refer to temperature difference per day, Δ T is designated as with reference to temperature difference per dayr;Inquiry Cable Structure location and meteorological data in recent years of the height above sea level interval in place no less than 2 years or actual measurement obtain the temperature of Cable Structure local environment with time and the delta data and changing rule of height above sea level, calculate maximum rate of change Δ T of the temperature for obtaining the Cable Structure local environment in recent years of Cable Structure location and place height above sea level interval no less than 2 years on height above sea levelh, Δ T is taken for convenience of narrationhUnit for DEG C/m;Taken on the surface of Cable Structure " R Cable Structure surface point ", the Specific Principles of " R Cable Structure surface point " are taken to be described in step b3, the temperature of this R Cable Structure surface point will be obtained by actual measurement below, it is called " R Cable Structure surface temperature measured data " to survey obtained temperature data, if utilizing the Thermodynamic calculation model of Cable Structure, the temperature of this R Cable Structure surface point is obtained by Calculation of Heat Transfer, it is called " R Cable Structure land surface pyrometer count evidence " just to calculate obtained temperature data;From the minimum height above sea level residing for Cable Structure to highest height above sea level,It is uniform in Cable Structure to choose no less than three different height above sea levels,At the height above sea level of each selection,Two points are at least chosen at the intersection on horizontal plane Yu Cable Structure surface,The exterior normal of straw line body structure surface at selected point,The exterior normal direction of all selections is referred to as in " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness ",Direction of the Cable Structure along the Temperature Distribution of wall thickness is measured with " horizontal plane and the intersection on Cable Structure surface " to intersect,The sunny slope exterior normal direction of Cable Structure and the in the shade face exterior normal direction of Cable Structure must be included in direction of the measurement Cable Structure along the Temperature Distribution of wall thickness of selection,No less than three points are chosen along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness is uniform in Cable Structure,For support cable a point is only taken along each direction of measurement Cable Structure along the Temperature Distribution of wall thickness,Only measure the temperature of the surface point of support cable,All temperature being selected a little of measurement,The temperature measured is referred to as " temperature profile data of the Cable Structure along thickness ",Wherein edge is intersected with same " horizontal plane and the intersection on Cable Structure surface "," temperature profile data of the Cable Structure along thickness " that " direction of the measurement Cable Structure along the Temperature Distribution of wall thickness " measurement is obtained,It is referred to as in the method " temperature profile data of the identical height above sea level Cable Structure along thickness ",If have chosen H different height above sea levels,At each height above sea level,It has chosen direction of the B measurement Cable Structure along the Temperature Distribution of wall thickness,Direction along each measurement Cable Structure along the Temperature Distribution of wall thickness have chosen E point in Cable Structure,Wherein H and E are not less than 3,B is not less than 2,It is equal to 1 for support cable E,The sum for counting in Cable Structure " point of the measurement Cable Structure along the temperature profile data of thickness " is HBE,The temperature of this HBE " point of the measurement Cable Structure along the temperature profile data of thickness " will be obtained by actual measurement below,It is called " HBE Cable Structure along thickness temperature measured data " to survey obtained temperature data,If utilizing the Thermodynamic calculation model of Cable Structure,Temperature of this HBE measurement Cable Structure along the point of the temperature profile data of thickness is obtained by Calculation of Heat Transfer,It is called " HBE Cable Structure calculates data along thickness temperature " just to calculate obtained temperature data;Require to choose a position according to meteorology measurement temperature in Cable Structure location, the temperature of environment where meeting the Cable Structure of meteorology measurement temperature requirement will be obtained in the actual measurement of this position;A position is chosen at the spacious unobstructed place in Cable Structure location, the position should it is annual can obtain each day the ground this getable day most sufficient sunshine, in the flat board of one piece of carbon steel material of position of sound production, referred to as reference plate, reference plate not can contact with ground, distance is not less than 1.5 meters to reference plate from the ground, the one side of the reference plate faces south, referred to as sunny slope, the sunny slope of reference plate is coarse and dark, the sunny slope of reference plate should it is annual can obtain each day one flat plate the ground this getable day most sufficient sunshine, the non-sunny slope of reference plate is covered with insulation material, real-time monitoring is obtained to the temperature of the sunny slope of reference plate;
b2:Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point in real time, monitoring obtains temperature profile data of the previously defined Cable Structure along thickness in real time simultaneously, while monitoring in real time obtains the temperature record of environment where meeting the Cable Structure of meteorology measurement temperature requirement;By monitor in real time obtain the Cable Structure that the same day is carved into after sunrise moment next day between 30 minutes at sunrise where environment temperature measured data sequence, the temperature measured data of environment is arranged according to time order and function order where the temperature measured data sequence of environment was carved into after the sunrise moment next day Cable Structure between 30 minutes at sunrise by the same day where Cable Structure, maximum temperature and minimum temperature in the temperature measured data sequence of environment where finding Cable Structure, the same day for subtracting environment where minimum temperature obtains Cable Structure with the maximum temperature in the temperature measured data sequence of environment where Cable Structure is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as environment maximum temperature difference, it is designated as Δ Temax;Rate of change of the temperature of environment where obtaining Cable Structure on the time is calculated by Conventional mathematical by the temperature measured data sequence of environment where Cable Structure, the rate of change is also with time change;By the measured data sequence for monitoring the temperature for obtaining the sunny slope that the same day is carved into reference plate after sunrise moment next day between 30 minutes at sunrise in real time, the measured data that the measured data sequence of the temperature of the sunny slope of reference plate was carved into the temperature of the sunny slope of the reference plate between 30 minutes after sunrise moment next day by the same day at sunrise is arranged according to time order and function order, maximum temperature and minimum temperature in the measured data sequence for the temperature for finding the sunny slope of reference plate, same day of temperature that the sunny slope that minimum temperature obtains reference plate is subtracted with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate is carved into maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, referred to as reference plate maximum temperature difference, it is designated as Δ Tpmax;The Cable Structure surface temperature measured data sequence for all R Cable Structure surface points that the same day is carved into after sunrise moment next day between 30 minutes at sunrise is obtained by monitoring in real time, there is R Cable Structure surface point just to have R Cable Structure surface temperature measured data sequence, the Cable Structure surface temperature measured data that each Cable Structure surface temperature measured data sequence was carved into after sunrise moment next day between 30 minutes by the same day of a Cable Structure surface point at sunrise is arranged according to time order and function order, find the maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence, the same day for subtracting the temperature that minimum temperature obtains each Cable Structure surface point with the maximum temperature in each Cable Structure surface temperature measured data sequence is carved into the maximum temperature difference after sunrise moment next day between 30 minutes at sunrise, there is R Cable Structure surface point just to there is the R same day to be carved into the maximum temperature difference numerical value after sunrise moment next day between 30 minutes at sunrise, maximum therein is referred to as Cable Structure surface maximum temperature difference, it is designated as Δ Tsmax;Calculated by each Cable Structure surface temperature measured data sequence by Conventional mathematical and obtain rate of change of the temperature on the time of each Cable Structure surface point, the temperature of each Cable Structure surface point on the time rate of change also with time change;Obtain the same day by monitoring in real time and be carved at sunrise after sunrise moment next day between 30 minutes, in synchronization, after HBE " temperature profile data of the Cable Structure along thickness ", calculate the difference of the maximum temperature and minimum temperature that amount at the height above sea level of each selection in BE " temperature profile data of the identical height above sea level Cable Structure along thickness ", the absolute value of this difference is referred to as " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", have chosen H different height above sea levels just has H " Cable Structure thickness direction maximum temperature difference at identical height above sea level ", maximum in this H " Cable Structure thickness direction maximum temperature difference at identical height above sea level " is called " Cable Structure thickness direction maximum temperature difference ", it is designated as Δ Ttmax
b3:Survey calculation obtains Cable Structure steady temperature data;First, it is determined that at the time of obtaining Cable Structure steady temperature data, the condition related at the time of Cable Structure steady temperature data to determining to obtain has six, Section 1 condition was carved into after sunrise moment next day between 30 minutes at sunset at the time of being and obtain Cable Structure steady temperature data between the same day, the sunset moment refers to that, according to the sunset moment on earth rotation and the meteorology of revolution rule determination, data can be inquired about or calculate by conventional meteorology obtaining the required sunset moment of each day;The a conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, reference plate maximum temperature difference Δ TpmaxWith Cable Structure surface maximum temperature difference Δ TsmaxAll it is not more than 5 degrees Celsius;The b conditions of Section 2 condition were carved at sunrise on the same day in this period after sunrise moment next day between 30 minutes, the environment maximum temperature difference Δ T obtained in above survey calculationemaxNo more than refer to temperature difference per day Δ Tr, and reference plate maximum temperature difference Δ TpmaxSubtract and be not more than Δ T after 2 degrees Celsiusemax, and Cable Structure surface maximum temperature difference Δ TsmaxNo more than Δ Tpmax;Need to only meet in a conditions and b conditions of Section 2 one is known as meeting Section 2 condition;Section 3 condition is that the temperature of environment is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time where Cable Structure at the time of Cable Structure steady temperature data are obtained;Section 4 condition is that at the time of Cable Structure steady temperature data are obtained, the temperature of each Cable Structure surface point in R Cable Structure surface point is not more than 0.1 degree Celsius per hour on the absolute value of the rate of change of time;Section 5 condition is that at the time of Cable Structure steady temperature data are obtained, the Cable Structure surface temperature measured data of each Cable Structure surface point in R Cable Structure surface point was carved into the minimum after sunrise moment next day between 30 minutes for the same day at sunrise;Section 6 condition is " Cable Structure thickness direction maximum temperature difference " Δ T at the time of Cable Structure steady temperature data are obtainedtmaxNo more than 1 degree Celsius;This method utilizes above-mentioned six conditions, any one in the following three moment is referred to as " the mathematics moment for obtaining Cable Structure steady temperature data ", the first moment is at the time of meeting Section 1 to the Section 5 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", second of moment is at the time of only meeting the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ", the third moment is while at the time of meeting Section 1 to the Section 6 condition in above-mentioned " condition related at the time of Cable Structure steady temperature data to determining to obtain ";It it is exactly the mathematics moment for obtaining Cable Structure steady temperature data at the time of acquisition Cable Structure steady temperature data when a moment in the physical record data moment during the mathematics moment for obtaining Cable Structure steady temperature data is exactly this method;If the mathematics moment for obtaining Cable Structure steady temperature data is not any one moment in the physical record data moment in this method, take this method closest to the mathematics moment for obtaining Cable Structure steady temperature data that physical record data at the time of to obtain at the time of Cable Structure steady temperature data;The amount that this method is used in measurement record at the time of obtaining Cable Structure steady temperature data carries out the monitoring analysis of Cable Structure relevant health;The Cable Structure temperature field that this method is approximately considered at the time of obtaining Cable Structure steady temperature data is in the Cable Structure temperature at stable state, i.e. this moment and not changed over time, and this moment is exactly " at the time of the obtaining Cable Structure steady temperature data " of this method;Then, according to Cable Structure heat-transfer character, utilize " the R Cable Structure surface temperature measured data " at the time of obtaining Cable Structure steady temperature data and " HBE Cable Structure along thickness temperature measured data ", utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution for obtaining Cable Structure at the time of Cable Structure steady temperature data are obtained is calculated by conventional heat transfer, now the temperature field of Cable Structure is calculated by stable state, calculating the temperature profile data of obtained Cable Structure at the time of Cable Structure steady temperature data are obtained includes the calculating temperature of R Cable Structure surface point in Cable Structure, the calculating temperature of R Cable Structure surface point is referred to as R Cable Structure steady-state surface temperature and calculates data, also include calculating temperature of the Cable Structure above selected HBE " point of the measurement Cable Structure along the temperature profile data of thickness ", the calculating temperature of HBE " points of the measurement Cable Structure along the temperature profile data of thickness " is referred to as " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data calculates data correspondent equal with R Cable Structure steady-state surface temperature, and " HBE Cable Structure along thickness temperature measured data " with " HBE Cable Structure along thickness temperature calculating data " correspondent equal when, the temperature profile data for calculating obtained Cable Structure at the time of Cable Structure steady temperature data are obtained is referred to as " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " now is referred to as " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure along thickness temperature measured data " is referred to as " HBE Cable Structure along thickness steady temperature measured data ";When " R Cable Structure surface point " is taken on the surface of Cable Structure, the quantity of " R Cable Structure surface point " must is fulfilled for three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the observed temperature linear interpolation of point adjacent with the arbitrfary point on Cable Structure surface during the temperature at any point on Cable Structure surface is by " R Cable Structure surface point " is obtained, the error of the temperature of the arbitrfary point and the actual temperature of the arbitrfary point on Cable Structure surface is not more than 5% on the Cable Structure surface that linear interpolation is obtained;Cable Structure surface includes support cable surface;Second condition is that the point being not less than in 4, and " R Cable Structure surface point " in same height above sea level in the quantity of the point of same height above sea level in " R Cable Structure surface point " is uniform along Cable Structure surface;Maximum Δ h in the absolute value of the difference of the height above sea level of " R Cable Structure surface point " along all Cable Structure surface points adjacent two-by-two of height above sea level is not more than 0.2 DEG C divided by Δ ThObtained numerical value, Δ T is taken for convenience of narrationhUnit for DEG C/m, for convenience of narration take Δ h unit be m;When the definition of " R Cable Structure surface point " along the Cable Structure surface point adjacent two-by-two of height above sea level refers to only consider height above sea level, a Cable Structure surface point is not present in " R Cable Structure surface point ", the height above sea level numerical value of the Cable Structure surface point is between the height above sea level numerical value of adjacent Cable Structure surface point two-by-two;3rd condition is inquiry or obtains Cable Structure location and the interval sunshine rule of place height above sea level by meteorology conventionally calculation, further according to the geometric properties and bearing data of Cable Structure, found in Cable Structure it is annual by the sunshine-duration most sufficient position of those surface points, in " R Cable Structure surface point " at least one Cable Structure surface point be in Cable Structure whole year by a point in those most sufficient surface points of sunshine-duration;
C. the Cable Structure steady temperature data obtained under original state are calculated according to " temperature survey of the Cable Structure of this method calculates method " direct measurement, Cable Structure steady temperature data under original state are referred to as initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector To”;Survey or consult reference materials and obtain the physical and mechanical properties parameter varied with temperature of various materials used in Cable Structure;Initial Cable Structure steady temperature data vector T is obtained in actual measurementoSynchronization, direct measurement, which is calculated, obtains the Initial cable forces of all support cables, composition Initial cable force vector Fo;According to include the data including Cable Structure design data, completion data obtain all support cables free state i.e. Suo Li for 0 when length, in free state when cross-sectional area and the unit length in free state weight, and obtain the temperature of all support cables during these three data, on this basis using the physical function parameter varied with temperature and mechanical property parameters of all support cables, calculated according to Typical physical and obtain all support cables in initial Cable Structure steady temperature data vector ToUnder the conditions of Suo Li unit lengths of all support cables when the cross-sectional area and Suo Li of all support cables are 0 when the length of all support cables, Suo Li are 0 when being 0 weight, successively composition support cable initial drift is vectorial, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and Initial cable force vector F of the element that initial drift is vectorial, the initial free unit length of initial free cross-sectional area vector sum weight is vectorial of support cableoElement coding rule it is identical;T is obtained in actual measurementoWhile, that is, obtaining initial Cable Structure steady temperature data vector ToAt the time of synchronization, direct measurement calculates the measured data for obtaining initial Cable Structure, the measured data of initial Cable Structure is to include Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinates data, initial Cable Structure angle-data, measured data including initial Cable Structure spatial data, while the measured data of initial Cable Structure is obtained, survey calculation obtains the data of the health status that can express support cable including the Non-destructive Testing Data of support cable, the data of the health status that can express support cable now are referred to as support cable initial health data;The monitored amount initial value vector C of initial value composition of all monitored amountso, it is monitored amount initial value vector CoCoding rule and M monitored amounts coding rules it is identical;Evaluation object initial damage vector d is set up using support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement datao, vectorial doRepresent with initial mechanical calculating benchmark model AoThe initial health of the evaluation object of the Cable Structure of expression;Evaluation object initial damage vector doElement number be equal to N, doElement and evaluation object be one-to-one relationship, vectorial doElement coding rule it is identical with the coding rule of evaluation object;If doThe corresponding evaluation object of some element be a support cable in cable system, then doThe element numerical value represent correspondence support cable initial damage degree, if the numerical value of the element is 0, it is intact to represent the support cable corresponding to the element, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to the element has completely lost bearing capacity, if its numerical value is between 0 and 100%, then it represents that the support cable loses the bearing capacity of corresponding proportion;;If doThe corresponding evaluation object of some element be some bearing some generalized displacement component, then doThe element numerical value represent this bearing the generalized displacement component initial value;If doThe corresponding evaluation object of some element be to take d in some load, this methodoThe element numerical value be 0, the initial value for representing the change of this load is 0;If without the initial generalized displacement measurement data of Cable Structure bearing or can consider the initial generalized displacement of Cable Structure bearing be 0 when, vectorial doIn each element numerical value related to Cable Structure generalized displacement of support take 0;If during data without the Non-destructive Testing Data of support cable and other health status that can express support cable, or can consider structure original state be not damaged without relaxed state when, vectorial doIn each element numerical value related to support cable take 0;Initial Cable Structure bearing generalized coordinates data refer to the bearing generalized coordinates data under Cable Structure design point, and the initial generalized displacement measurement data of Cable Structure bearing, which refers to, is setting up initial mechanical calculating benchmark model AoWhen, Cable Structure bearing is relative to the generalized displacement that the bearing under Cable Structure design point occurs;
D. the physical and mechanical properties parameter varied with temperature of various materials, initial Cable Structure steady temperature data vector T according to used in the design drawing of Cable Structure, the measured data of as-built drawing and initial Cable Structure, support cable initial health data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable StructureoAll Cable Structure data obtained with preceding step, set up the initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data "o, based on AoCalculate obtained Cable Structure calculate data must closely its measured data, difference therebetween cannot be greater than 5%;Corresponding to Ao" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector To”;Corresponding to AoEvaluation object health status evaluation object initial damage vector doRepresent;Corresponding to AoAll monitored amounts initial value with monitored amount initial value vector CoRepresent;The current initial mechanical calculating benchmark model A for the Cable Structure for being included in " Cable Structure steady temperature data " is set up for the first timet o, the monitored current initial value vector C of amountt o" current initial Cable Structure steady temperature data vector Tt o”;The current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timet oWith monitored amount current initial value vector Ct oWhen, the current initial mechanical calculating benchmark model A of Cable Structuret oIt is equal to the initial mechanical calculating benchmark model A of Cable Structureo, it is monitored the current initial value vector C of amountt oIt is equal to monitored amount initial value vector Co;At oCorresponding " Cable Structure steady temperature data " are referred to as " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector Tt o", the current initial mechanical calculating benchmark model A of Cable Structure is set up for the first timet oWhen, Tt oIt is equal to To;At oEvaluation object initial health and AoEvaluation object health status it is identical, also with evaluation object initial damage vector doRepresent, the A in cyclic process belowt oEvaluation object initial health all the time use evaluation object initial damage vector doRepresent;ToAnd doIt is AoParameter, by AoThe obtained initial value of all monitored amounts of Mechanics Calculation result and CoThe initial value of all monitored amounts represented is identical, therefore alternatively CoBy AoMechanics Calculation result composition;Tt oAnd doIt is At oParameter, Ct oBy At oMechanics Calculation result composition;
E. the circulation walked by e to the n-th step is entered from here;During structure military service, constantly the current data of " Cable Structure steady temperature data " is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data of " Cable Structure steady temperature data " is referred to as " current Cable Structure steady temperature data ", is designated as " current Cable Structure steady temperature data vector Tt", vector TtDefinition mode and vector ToDefinition mode it is identical;Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, actual measurement obtain all M in Cable Structure1The rope force data of root support cable, all these rope force data composition current cable force vector F, vectorial F element and vector FoElement coding rule it is identical;Current Cable Structure steady temperature data vector T is obtained in actual measurementtSynchronization, Actual measurement obtains all M1The space coordinate of two supporting end points of root support cable, the difference of the space coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables constitute the current supporting end points of support cable two horizontal range vector, the coding rule and Initial cable force vector F of the element of the current supporting end points of support cable two horizontal range vectoroElement coding rule it is identical;
F. according to current Cable Structure steady temperature data vector Tt, current initial mechanical calculating benchmark model A is updated according to step f1 to f3t o, the monitored current initial value vector C of amountt oWith current initial Cable Structure steady temperature data vector Tt o
F1. T is comparedtWith Tt oIf, TtEqual to Tt o, then At o、Ct oAnd Tt oKeep constant;Otherwise need to follow these steps to At o、Ct oAnd Tt oIt is updated;
F2. T is calculatedtWith ToDifference, TtWith ToDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TtWith ToDifference represented with steady temperature change vector S, S be equal to TtSubtract To, S represents the change of Cable Structure steady temperature data;
F3. to AoIn Cable Structure apply temperature change, the numerical value of the temperature change of application is just derived from steady temperature change vector S, to AoIn the temperature change that applies of Cable Structure after the current initial mechanical calculating benchmark model A that is updatedt o, update At oWhile, Tt oAll elements numerical value also uses TtAll elements numerical value correspondence replace, that is, have updated Tt o, thus obtained properly corresponding to At oTt o;Update Ct oMethod be:As renewal At oAfterwards, A is obtained by Mechanics Calculationt oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Ct o;At oSupport cable initial health all the time use evaluation object initial damage vector doRepresent;
G. in current initial mechanical calculating benchmark model At oOn the basis of carry out Mechanics Calculation several times according to step g1 to g4, pass through to calculate and obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector Du
G1. Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is to constantly update, that is, is updating current initial mechanical calculating benchmark model At o, the monitored current initial value vector C of amountt oWith current initial Cable Structure steady temperature data vector Tt oAfterwards, it is necessary to then update Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector Du
G2. in the current initial mechanical calculating benchmark model A of Cable Structuret oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating;According to the coding rule of evaluation object, calculated successively;Calculate each time and assume that only one of which evaluation object is further added by unit damage or unit generalized displacement or load unit change on the basis of original damage or generalized displacement or load, specifically, if the evaluation object is a support cable in cable system, then it is assumed that the support cable is in vectorial doThe support cable represented is further added by unit damage on the basis of having damaged, if the evaluation object is the generalized displacement component in a direction of a bearing, it is assumed that the bearing is further added by unit generalized displacement in the direction of displacement, if the evaluation object is a load, it is assumed that the load is in vectorial doLoad unit change is further added by the basis of the existing variable quantity of the load represented, D is usedukThis increased unit damage or unit generalized displacement or load unit change are recorded, wherein k represents to increase the numbering of the evaluation object of unit damage or unit generalized displacement or load unit change, DukIt is evaluation object unit change vector DuAn element, evaluation object unit change vector DuElement coding rule and vector doElement coding rule it is identical;The evaluation object of increase unit damage or unit generalized displacement or load unit change during the evaluation object of increase unit damage or unit generalized displacement or load unit change is calculated different from other times in calculating each time, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, the current calculated value for calculating obtained all monitored amounts each time constitutes a monitored amount calculation current vector, is monitored the element number rule and monitored amount initial value vector C of amount calculation current vectoroElement number rule it is identical;
G3. obtained monitored amount calculation current vector is calculated each time subtracts the monitored current initial value vector C of amountt oObtain a vector, each element of the vector divided by the assumed unit damage of this time calculating or unit generalized displacement or load unit are changed into numerical value again, a monitored amount unit change vector is obtained, has N number of evaluation object just to have N number of monitored amount unit change vector;
G4. by coding rule of this N number of monitored amount unit change vector according to N number of evaluation object, composition has the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C that N is arranged successively;Cable Structure unit damage monitored numerical quantity transformation matrices Δ C each row correspond to a monitored amount unit change vector;Cable Structure unit damage monitored numerical quantity transformation matrices Δ C every a line corresponds to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit generalized displacement or load unit change;The coding rule of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C row and vector doElement coding rule it is identical, the coding rule of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C row is identical with the coding rules of M monitored amounts;
H. current Cable Structure steady temperature data vector T is obtained in actual measurementtWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TtAt the time of synchronization Cable Structure all monitored amounts current measured value, the monitored amount current value vector C of composition;The monitored amount current value vector C and monitored current initial value vector C of amountt oWith monitored amount initial value vector CoDefinition mode it is identical, the same monitored amount of element representation of three vectorial identical numberings is in concrete numerical value not in the same time;
I. the current nominal fatigue vector d of evaluation object is defined, the current nominal fatigue vector d of evaluation object element number is equal to the quantity of evaluation object, it is one-to-one relationship between the current nominal fatigue vector d of evaluation object element and evaluation object, the current nominal fatigue vector d of evaluation object element numerical value represents the nominal fatigue degree or nominal generalized displacement or nominal load variable quantity of correspondence evaluation object;The coding rule of vectorial d element and vector doElement coding rule it is identical;
J. according to monitored amount current value vector C with the monitored current initial value vector C of amountt o, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and the current nominal fatigue vector d of evaluation object to be asked, the linear approximate relationship can be expressed as formula 1, other amounts in formula 1 in addition to d are, it is known that solution formula 1 can just calculate the current nominal fatigue vector d of evaluation object;
C = C o t + Δ C · d Formula 1
K. the currently practical injury vector d of evaluation object is defineda, the currently practical injury vector d of evaluation objectaElement number be equal to evaluation object quantity, the currently practical injury vector d of evaluation objectaElement and evaluation object between be one-to-one relationship, the currently practical injury vector d of evaluation objectaElement numerical value represent correspondence evaluation object actual damage degree or actual generalized displacement or real load variable quantity;Vectorial daElement coding rule and vector doElement coding rule it is identical;
L. the currently practical injury vector d of evaluation object expressed using formula 2aK-th of element da kWith evaluation object initial damage vector doK-th of element dokWith evaluation object current nominal fatigue vector d k-th of element dkBetween relation, calculating obtain the currently practical injury vector d of evaluation objectaAll elements;
K=1,2 in formula 2,3 ... ..., N, da kRepresent the currently practical health status of k-th of evaluation object, da kFor 0 when represent k-th of evaluation object without health problem, da kNumerical value represents that k-th of evaluation object is the evaluation object of unsoundness problem when not being 0, if the evaluation object is a support cable in cable system, then da kThe order of severity of its current health problem is represented, the support cable of unsoundness problem is probably slack line, is also likely to be damaged cable, da kThe numerical response degree of relaxation or the damage of the support cable;Damaged cable is identified from the support cable of these unsoundness problems, remaining is exactly slack line, the currently practical injury vector d of evaluation objectaIn correspond to slack line element numerical expression be with slack line relax degree mechanic equivalent currently practical equivalent damage degree;If the evaluation object is a generalized displacement component of a bearing, then da kRepresent its currently practical generalized displacement numerical value;
M. utilize in current Cable Structure steady temperature data vector TtUnder the conditions of, in l walk the slack line that identifies and with the currently practical injury vector d of evaluation objectaExpression these slack lines, the degree that relaxed with it mechanic equivalent currently practical equivalent damage degree, using e walk acquisition in current Cable Structure steady temperature data vector TtUnder the conditions of current cable force vector F and the supporting end points horizontal range vector of current support cable two, using c walk obtain in initial Cable Structure steady temperature data vector ToUnder the conditions of support cable initial drift is vectorial, the weight of the initial free unit length of initial free cross-sectional area vector sum is vectorial, Initial cable force vector Fo, utilize current Cable Structure steady temperature data vector TtThe support cable current steady state temperature data of expression, using c walk obtain in initial Cable Structure steady temperature data vector ToThe support cable initial steady state temperature data of expression, utilize the physical and mechanical properties parameter varied with temperature that various materials used in the Cable Structure of acquisition are walked in c, it is included in influence of the temperature change to support cable physics, mechanics and geometric parameter, calculated by the way that slack line is carried out into mechanic equivalent with damaged cable slack line, with the equivalent relaxation degree of currently practical equivalent damage degree, mechanic equivalent condition is:First, two equivalent ropes without relaxation with not damaged when initial drift, geometrical property parameter, the mechanics parameters of density and material it is identical;2nd, after relaxing or damage, two equivalent slack lines are identical with the overall length after deformation with the Suo Li of damage rope;When meeting above-mentioned two mechanic equivalent condition, mechanics function of such two support cables in Cable Structure is exactly identical, if replaced with equivalent slack line after damaged cable, any change will not occur for Cable Structure, and vice versa;Those relaxation degree for being judged as slack line are tried to achieve according to foregoing mechanic equivalent condition, relaxation degree is exactly the knots modification of support cable drift, that is, the long adjustment amount of rope of those support cables that need to adjust Suo Li is determined;So it is achieved that relaxation identification and the non-destructive tests of support cable;Institute's demand power is provided by current cable force vector F corresponding elements during calculating;Damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method, so, this method is according to the currently practical injury vector d of evaluation objectaProblem cable can either be identified, also can determine which bearing there occurs generalized displacement and its numerical value;
N. e steps are returned to, start the circulation next time walked by e to the n-th step.
CN201510439657.2A 2015-07-23 2015-07-23 Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process Pending CN105067342A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510439657.2A CN105067342A (en) 2015-07-23 2015-07-23 Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510439657.2A CN105067342A (en) 2015-07-23 2015-07-23 Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process

Publications (1)

Publication Number Publication Date
CN105067342A true CN105067342A (en) 2015-11-18

Family

ID=54496766

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510439657.2A Pending CN105067342A (en) 2015-07-23 2015-07-23 Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process

Country Status (1)

Country Link
CN (1) CN105067342A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3902359A (en) * 1973-12-14 1975-09-02 Rucker Co Ton cycle indicator
DE10312052A1 (en) * 2003-03-18 2004-09-30 Christa Reiners Arrangement for testing the stability under load and bending strength of a mast, whereby a force application unit is provided and displacement measurement sensors for both the mast itself and a neighboring mast
CN103852331A (en) * 2014-03-10 2014-06-11 东南大学 Progressive recognition method for problem cable load generalized displacement based on mixing monitoring
CN103913338A (en) * 2014-03-10 2014-07-09 东南大学 Hybrid monitoring defective cable load generalized displacement identification method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3902359A (en) * 1973-12-14 1975-09-02 Rucker Co Ton cycle indicator
DE10312052A1 (en) * 2003-03-18 2004-09-30 Christa Reiners Arrangement for testing the stability under load and bending strength of a mast, whereby a force application unit is provided and displacement measurement sensors for both the mast itself and a neighboring mast
CN103852331A (en) * 2014-03-10 2014-06-11 东南大学 Progressive recognition method for problem cable load generalized displacement based on mixing monitoring
CN103913338A (en) * 2014-03-10 2014-07-09 东南大学 Hybrid monitoring defective cable load generalized displacement identification method

Similar Documents

Publication Publication Date Title
CN105067032A (en) Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined linear displacement
CN105067033A (en) Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined linear displacement
CN105067300A (en) Simplified cable-force-monitoring defective cable load generalized displacement progressive identifying method
CN105067342A (en) Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process
CN105136497A (en) Simplified cable force monitoring damaged cable load generalized displacement progressive identification method
CN105067343A (en) Method for recognizing load angular displacement of damaged cable based on streamlined hybrid monitoring process
CN105136494A (en) Simplified space coordinate monitoring faulty cable load generalized displacement identification method
CN105067320A (en) Method for recognizing load of damaged cable based on angle monitoring process of streamlined and generalized displacement
CN105067332A (en) Method for recognizing load linear displacement of damaged cable based on streamlined angle monitoring process
CN105136491A (en) Simplified strain monitoring faulty cable load identification method
CN105067335A (en) Method for recognizing load generalized displacement of damaged cable based on streamlined space coordinate monitoring process
CN105067313A (en) Method for recognizing load of damaged cable based on angle monitoring process of streamlined angular displacement
CN105067356A (en) Method for recognizing load of damaged cable based on streamlined hybrid monitoring process
CN105136490A (en) Simplified cable force monitoring faulty cable load linear displacement identification method
CN105136498A (en) Faulty cable load generalized displacement progressive identification method based on angle monitoring
CN105067347A (en) Method for progressively recognizing load linear displacement of damaged cable based on streamlined hybrid monitoring process
CN105067354A (en) Method for progressively recognizing load of damaged cable based on streamlined hybrid monitoring process of generalized displacement
CN105067305A (en) Simplified generalized displacement strain-monitoring damaged cable load identifying method
CN105067298A (en) Simplified generalized displacement strain-monitoring defective cable load progressive identifying method
CN105136492A (en) Simplified linear displacement angle monitoring faulty cable load progressive identification method
CN105067318A (en) Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined linear displacement
CN105067297A (en) Simplified cable-force-monitoring damaged cable load support angular displacement identifying method
CN105067337A (en) Method for progressively recognizing load generalized displacement of damaged cable based on hybrid monitoring process
CN105136499A (en) Simplified space coordinate monitoring faulty cable load progressive identification method
CN105067314A (en) Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined annular displacement

Legal Events

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

Application publication date: 20151118