CN105067328A - Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process - Google Patents
Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process Download PDFInfo
- Publication number
- CN105067328A CN105067328A CN201510439492.9A CN201510439492A CN105067328A CN 105067328 A CN105067328 A CN 105067328A CN 201510439492 A CN201510439492 A CN 201510439492A CN 105067328 A CN105067328 A CN 105067328A
- Authority
- CN
- China
- Prior art keywords
- cable structure
- temperature
- data
- vector
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention provides a method for progressively recognizing the load linear displacement of a damaged cable based on the cable force monitoring process. According to the method, based on the cable force 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
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.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, displacement of the lines may occur for Cable Structure bearing, the load that Cable Structure is born is it can also happen that change, the health status of Cable Structure may also change simultaneously, in this complex condition, this method recognizes bearing displacement of the lines and damaged cable (this method is referred to as the health status of core evaluation object), category engineering structure health monitoring field based on cable force monitoring (monitored Suo Li is referred to as " monitored amount " by this method).
Background technology
Reject load change, Cable Structure bearing displacement of the lines and structure temperature and change influence to Cable Structure health status recognition result, so that the change of the health status of identification structure 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 support wire displacement identification result, so as to recognize Cable Structure bearing displacement of the lines 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
Technical problem:This method discloses a kind of method, under conditions of cost is lower, two kinds of functions are realized, are respectively, first, the influence of bearing displacement of the lines, 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 support wire displacement identification result, so as to recognize Cable Structure bearing displacement of the lines exactly.
Technical scheme:In the method, coordinate of the bearing on the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing space coordinate ", space coordinate of the bearing on X, Y, Z axis can also be said to be, bearing is referred to as in space coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the space coordinate of some axle with a space coordinate subscale of bearing on the concrete numerical value of the space coordinate of some axle;Angular coordinate of the bearing on X, Y, Z axis is censured with " bearing angular coordinate ", bearing is referred to as in angular coordinate component of the bearing on the axle, this method also reaching concrete numerical value of the bearing on the angular coordinate of some axle with an angular coordinate subscale of bearing on the concrete numerical value of the angular coordinate of some axle;Censure that bearing angular coordinate and bearing space coordinate are all with " bearing generalized coordinates ", in this method also with a generalized coordinates subscale of bearing up to bearing on the space coordinate of axle or the concrete numerical value of angular coordinate;Bearing is referred to as bearing displacement of the lines on the change of the coordinate of X, Y, Z axis, it may also be said to which the change of bearing space coordinate is referred to as the concrete numerical value in bearing displacement of the lines, this method also with the translational component expression bearing of bearing on the displacement of the lines of some axle;Bearing is referred to as the concrete numerical value in angular displacement of support, this method also with the angular displacement component expression bearing of bearing on the angular displacement of some axle on the change of the angular coordinate of X, Y, Z axis;Generalized displacement of support censures bearing displacement of the lines and angular displacement of support is all, and also bearing is reached on the displacement of the lines of some axle or the concrete numerical value of angular displacement with a generalized displacement subscale of bearing in this method;Bearing displacement of the lines is alternatively referred to as translational displacement, and support settlement is the component of bearing displacement of the lines or translational displacement in gravity direction.
The external force that object, structure are born can be described as load, and load includes face load and volume load.Face load is also known as surface load, is the load for acting on body surface, including two kinds of concentrfated load and distributed load.Volume load is the continuously distributed load in interior of articles each point, the deadweight of such as object and inertia force.
Concentrfated load is divided into two kinds of concentrated force and concentrated couple, in a coordinate system, for example in Descartes's rectangular coordinate system, one concentrated force can resolve into three components, likewise, a concentrated couple can also resolve into three components, if load is actually concentrfated load, force component or a concentrated couple component is concentrated to be referred to as a load by one in the method, the now change of load is embodied as a change for concentrating force component or a concentrated couple component.
Distributed load is divided into line distributed load and EDS maps load, the size of the zone of action and distributed load of the description of distributed load at least including distributed load, the size of distributed load is expressed with distribution intensity, being distributed intensity, to express, (such as two distributed loads are all uniform with distribution characteristics (such as uniform, SIN function equal distribution feature) and amplitude, but its amplitude is different, and the concept of amplitude can be illustrated by taking well-distributed pressure as an example:Same structure bears two different well-distributed pressures, and two distributed loads are all uniform loads, but the amplitude of a distributed load is 10MPa, and the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, actually refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, one distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is not all identical, the component of this several distributed load is so regarded as same amount of independent distributed load in the method, now a load just represents the component of a distributed load, the amplitude changing ratio identical component of wherein distribution intensity can also be synthesized into a distributed load or is load.
Volume load is the continuously distributed load in interior of articles each point, deadweight and inertia force such as object, the size of the zone of action and volume load of the description of volume load at least including volume load, the size of volume load is expressed with distribution intensity, with distribution characteristics (such as uniform, linear function equal distribution feature) and amplitude, to express, (it is all uniform that lotuses are storaged in such as two individuals to distribution intensity, but its amplitude is different, and the concept of amplitude is illustrated exemplified by can conducting oneself with dignity:The material of two parts of same structure is different, therefore density is different, so while the volume load suffered by the two parts is all uniform, but the amplitude of the volume load suffered by a part is probably 10kN/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 bearing translational component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of bearing translational 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 bearing translational component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of bearing translational 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 ";If having M in cable system1Root support cable, Cable Structure rope force data includes this M1The Suo Li of root support cable, this method is in monitoring whole M1On the basis of root supporting cable force, the artificial increase M in Cable Structure2Root rope, referred to as senses rope, and this M newly increased will be monitored during cable structure health monitoring2Root senses the Suo Li of rope;The monitored amount of summary, the M Suo Li that whole Cable Structure has M root ropes is monitored, that is, has M monitored amounts, and wherein M is M1With M2Sum;M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object;The M newly increased2The rigidity of root sensing rope, should be much smaller compared with the rigidity of any one support cable of Cable Structure;The M newly increased2The Suo Li of each sensing rope of root sensing rope should be more much smaller than the Suo Li of any one support cable of Cable Structure, even if can so ensure this M newly increased2Root sensing rope occurs in that damage or relaxation, and the influence to the stress of other components of Cable Structure, strain, deformation is very little;The M newly increased2Direct stress should be less than its fatigue limit on the cross section of root sensing rope, and these requirements can ensure the M newly increased2Fatigue damage occurs for root sensing Suo Buhui;The M newly increased2The two ends of root sensing rope should fully be anchored, it is ensured that be not in relaxation;The M newly increased2Root sensing rope should obtain sufficient anti-corrosion protection, it is ensured that the M newly increased2Damage and relaxation occur for root sensing Suo Buhui;For convenience, " monitored all parameters of Cable Structure " are referred to as in the method " monitored amount ";To M monitored amount serial numbers, this method is with this numbering, j=1,2,3 ..., M is represented with variable j, and the numbering will be used to generate vector sum matrix in subsequent step;The M newly increased in the method2Root senses a part for Suo Zuowei Cable Structures, when mentioning Cable Structure again hereinafter, and Cable Structure includes increase M2Cable Structure and the M newly increased before root sensing rope2Root senses rope, that is to say, that refer to the M including newly increasing when mentioning Cable Structure hereinafter2Root senses the Cable Structure of rope;Therefore mention hereinafter when obtaining " Cable Structure steady temperature data " according to the temperature survey of the Cable Structure of this method " calculate method " survey calculation, Cable Structure therein includes the M newly increased2Root senses rope, and obtained " Cable Structure steady temperature data " include the M newly increased2Root senses the steady temperature data of rope, obtains the M newly increased2The method of the steady temperature data of root sensing rope is same as the M of Cable Structure1The preparation method of the steady temperature data of root support cable, is no longer handed over one by one later;Measure the M newly increased2The Suo Li of root sensing rope method is same as the M of Cable Structure1The Suo Li of root support cable measuring method, is no longer handed over one by one later;When carrying out any measurement to the support cable of Cable Structure, while to the M newly increased2The measurements of root sensing Suo Jinhang equally, are no longer handed over one by one later;The M newly increased2Root senses rope in addition to not damaging and relaxing, to the M newly increased2Requirement and preparation method of the requirement and preparation method of the information content of root sensing rope with the information content of the support cable of Cable Structure are identical, no longer hand over one by one later;When setting up the various mechanical models of Cable Structure later, by the M newly increased2The support cable that root sensing rope treats as Cable Structure is treated;Below, in addition to the occasion of damage and relaxation for mentioning support cable, the M that described support cable includes the support cable of Cable Structure and newly increased when mentioning support cable2Root senses rope;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;T is obtained in actual measurementoWhile, 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 linear movement measuring 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;Using the initial linear movement measuring data of support cable initial health data, Cable Structure bearing can be expressed and Cable Structure load measurement data set up evaluation object initial damage vector do, 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 translational component, then doThe element numerical value represent this bearing the translational 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 linear movement measuring data of Cable Structure bearing or can consider the initial displacement of the lines of Cable Structure bearing be 0 when, vectorial doIn each element numerical value related to Cable Structure bearing displacement of the lines 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 spatial data refers to the bearing spatial data under Cable Structure design point, and the initial linear movement measuring data of Cable Structure bearing, which refer to, is setting up initial mechanical calculating benchmark model AoWhen, Cable Structure bearing is relative to the displacement of the lines 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 linear movement measuring 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;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, A in the methodo、Co、doAnd ToIt is constant;
E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated;Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingi o, AoAnd Ai oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated;When ith circulation starts, corresponding to Ai o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Ti oRepresent, vector Ti oDefinition mode and vector ToDefinition mode it is identical, Ti oElement and ToElement correspond;The current initial damage vector of evaluation object that ith circulation needs when starting is designated as di o, di oRepresent Cable Structure A during this time circulation beginningi oEvaluation object health status, di oDefinition mode and doDefinition mode it is identical, di oElement and doElement correspond;When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amounti oRepresent, vectorial Ci oDefinition mode and vector CoDefinition mode it is identical, Ci oElement and CoElement correspond, be monitored the current initial value vector C of amounti oRepresent to correspond to Ai oAll monitored amounts concrete numerical value;Ti oAnd di oIt is Ai oCharacterisitic parameter, Ci oBy Ai oMechanics Calculation result composition;When circulation starts for the first time, Ai oIt is designated as A1 o, set up A1 oMethod to make A1 oEqual to Ao;When circulation starts for the first time, Ti oIt is designated as T1 o, set up T1 oMethod to make T1 oEqual to To;When circulation starts for the first time, di oIt is designated as d1 o, set up d1 oMethod to make d1 oEqual to do;When circulation starts for the first time, Ci oIt is designated as C1 o, set up C1 oMethod to make C1 oEqual to Co;
F. enter from here and the circulation walked to q is walked by f;During structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Ti, vector TiDefinition mode and vector ToDefinition mode it is identical, TiElement and ToElement correspond;Vector T is obtained in actual measurementiWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TiAt the time of synchronization Cable Structure in all monitored amounts currency, the monitored amount current value vector C of all these numerical value compositioni, vectorial CiDefinition mode and vector CoDefinition mode it is identical, CiElement and CoElement correspond, represent identical monitored amount in numerical value not in the same time;Current Cable Structure steady temperature data vector T is obtained in actual measurementiWhile, to the M newly increased2Root senses Suo Jinhang Non-Destructive Testings, therefrom identify the sensing rope for occurring damaging or relaxing, according to monitored amount coding rule, appearances that removes with identifying is damaged or the corresponding element of loose sensing rope in each vector numbered according to monitored amount coding rule occurred before this method, also no longer occur in each vector sum matrix occurred after this method the appearance with identifying damage or relaxation the corresponding element of sensing rope, no longer include being identified the sensing rope for occurring damaging or relaxing here when mentioning sensing rope after this method, no longer include the Suo Li for being identified the sensing rope for occurring damaging or relaxing here when mentioning monitored amount after this method;Several sensing ropes for occurring damaging or relaxing are identified from Cable Structure, just by M2Reduce same quantity with M;
G. according to current Cable Structure steady temperature data vector Ti, current initial mechanical calculating benchmark model A is updated according to step g1 to g3i o, the monitored current initial value vector C of amounti oWith current initial Cable Structure steady temperature data vector Ti o, and the current initial damage vector d of evaluation objecti oKeep constant;
G1. T is comparediWith Ti oIf, TiEqual to Ti o, then Ai o、Ci oAnd Ti oKeep constant;Otherwise need to follow these steps to Ai o、Ci oAnd Ti oIt is updated;
G2. T is calculatediWith ToDifference, TiWith ToDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TiWith ToDifference represented with steady temperature change vector S, S be equal to TiSubtract To, S represents the change of Cable Structure steady temperature data;
G3. 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 updatedi o, update Ai oWhile, Ti oAll elements numerical value also uses TiAll elements numerical value correspondence replace, that is, have updated Ti o, thus obtained properly corresponding to Ai oTi o;Now di oKeep constant;As renewal Ai oAfterwards, Ai oRope the health status current initial damage vector d of evaluation objecti oRepresent, Ai oCable Structure steady temperature with current Cable Structure steady temperature data vector TiRepresent, update Ci oMethod be:As renewal Ai oAfterwards, A is obtained by Mechanics Calculationi oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Ci o;
H. in current initial mechanical calculating benchmark model Ai oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingiWith evaluation object unit change vector Di u;
H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4iAnd Di u;At other moment, when in step g to Ai oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CiAnd Di uIf, not to A in step gi oIt is updated, then is directly transferred to step i here and carries out follow-up work;
H2. in current initial mechanical calculating benchmark model Ai oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating;According to the coding rule of evaluation object, calculated successively;Calculate each time and assume that only one of which evaluation object is further added by unit damage or load unit change on the basis of original damage or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is the translational component in a direction of a bearing, it is assumed that the bearing is further added by unit displacement of the lines in the direction of displacement, if the evaluation object is a load, it is assumed that the load is further added by load unit change, D is usedi ukThis increased unit damage or unit displacement of the lines or load unit change is recorded, wherein k represents the numbering for increasing unit damage or unit displacement of the lines or the evaluation object of load unit change, Di ukIt is evaluation object unit change vector Di uAn element, evaluation object unit change vector Di uElement coding rule and vector doElement coding rule it is identical;The evaluation object that unit damage or unit displacement of the lines or load unit change are further added by calculating each time is different from the evaluation object that unit damage or unit displacement of the lines or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector;When assuming that k-th of evaluation object is further added by unit damage or unit displacement of the lines or load unit change, C is usedi tkRepresent corresponding " monitored amount calculation current vector ";When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed;Ci tkDefinition mode and vector CoDefinition mode it is identical, Ci tkElement and CoElement correspond;
H3. obtained vectorial C is calculated each timei tkSubtract vectorial Ci oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or unit displacement of the lines or load unit change numerical valuei k”;There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount ";
H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "i”;Unit damage monitored numerical quantity transformation matrices Δ CiEach row correspond to a monitored amount unit change vector;Unit damage monitored numerical quantity transformation matrices Δ CiEvery a line correspond to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit displacement of the lines or load unit change;Unit damage monitored numerical quantity transformation matrices Δ CiRow coding rule and vector doElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CiRow coding rule and M monitored amounts coding rules it is identical;
I. current nominal fatigue vector d is definedi cWith currently practical injury vector di, di cAnd diElement number be equal to evaluation object quantity, di cAnd diElement and evaluation object between be one-to-one relationship, di cElement numerical value represent the nominal fatigue degree or nominal displacement of the lines or nominal load variable quantity of correspondence evaluation object, di cAnd diWith evaluation object initial damage vector doElement number rule it is identical, di cElement, diElement and doElement be one-to-one relationship;
J. according to monitored amount current value vector CiWith " the monitored current initial value vector C of amounti o", " unit damage monitored numerical quantity transformation matrices Δ Ci" and " current nominal fatigue vector di c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1i cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector di c;
Formula 1
K. the currently practical injury vector d expressed using formula 2iK-th of element di kWith the current initial damage vector d of evaluation objecti oK-th of element di okWith current nominal fatigue vector di cK-th of element di ckBetween relation, calculating obtain currently practical injury vector diAll elements;
K=1,2,3 in formula 2 ..., N;di kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then di kRepresent its currently practical damage, di kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity;If the evaluation object is a translational component of a bearing, then di kRepresent its currently practical displacement of the lines numerical value;So far this method realizes the identification of the health status of core evaluation object;
L. current nominal fatigue vector d is tried to achievei cAfterwards, mark vector B is set up according to formula 3i, formula 4 gives mark vector BiK-th of element definition;
Formula 3
Element B in formula 4i kIt is mark vector BiK-th of element, Di ukIt is evaluation object unit change vector Di uK-th of element, di ckIt is the current nominal fatigue vector d of evaluation objecti cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4;
If m. mark vector BiElement be all 0, then return to step f continue this circulation;If mark vector BiElement be not all 0, then enter next step, i.e. step n;
N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5i+1 oEach element;
D in formula 5i+1 okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationi+1 oK-th of element, di okThis, i.e. the current initial damage vector d of evaluation object of ith circulationi oK-th of element, Di ukIt is the evaluation object unit change vector D of ith circulationi uK-th of element, Bi kIt is the mark vector B of ith circulationiK-th of element, k=1,2,3 ... ..., N in formula 5;
O. remove once, i.e. i+1 time circulates required current initial Cable Structure steady temperature data vector Ti+1 oThe current initial Cable Structure steady temperature data vector T circulated equal to ithi o;
P. in initial mechanical calculating benchmark model AoOn the basis of, 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, then it is d to make the health status of ropei+1 oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationi+1;Obtain Ai+1Afterwards, A is obtained by Mechanics Calculationi+1In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationi+1 o;
Q. step f is returned to, starts to circulate next time.
Beneficial effect:Structural healthy monitoring system is monitored on-line for a long time by using sensor to structural response first, (or offline) analysis online is carried out after acquisition Monitoring Data to it and obtains structural health conditions data, due to the complexity of structure, structural healthy monitoring system needs to use the equipment such as substantial amounts of sensor to carry out monitoring structural health conditions, therefore 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 for the health status for recognizing core evaluation object.When it is implemented, the following steps are one kind in the various steps that can be taken.
The first step:First confirm that the quantity for the load that the possibility that Cable Structure is born changes.The characteristics of load born according to Cable Structure, confirmation wherein " is possible to the load changed ", or all load is considered as " being possible to the load changed ", if the shared JZW load that may be changed, that is, have JZW secondary evaluation objects.
If the quantity of the bearing translational 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.
If evaluated support cable and the quantity sum of bearing translational component are that P, the i.e. quantity of core evaluation object are P, if the quantity of evaluated bearing translational component is Z, if the quantity of evaluated support cable is M1, structure rope force data include this M1The Suo Li of root support cable, this method is in monitoring whole M1On the basis of root supporting cable force, increase M2Other individual monitored amounts.
Increased M2Other individual monitored amounts are still Suo Li, are described below:
It is first artificial in Cable Structure to increase M before the start-up operation of structural health detecting system2(M2Not less than Z+4) root rope, rope is referred to as sensed, the M newly increased2The rigidity of root sensing rope, should be much smaller compared with the rigidity of any one support cable of Cable Structure, such as the small M 20 times, newly increased2The Suo Li of root sensing rope should be smaller, and for example its cross section direct stress should be less than its fatigue limit, and these requirements can ensure the M newly increased2Fatigue damage, the M newly increased occur for root sensing Suo Buhui2The two ends of root sensing rope should fully be anchored, it is ensured that be not in relaxation, the M newly increased2Root sensing rope should obtain sufficient anti-corrosion protection, it is ensured that the M newly increased2Damage and relaxation occur for root sensing Suo Buhui, and this M newly increased will be monitored during monitoring structural health conditions2Root senses the Suo Li of rope.
The reliability of health monitoring can also be ensured by the way of many increase sensing ropes, for example, make M2Not less than Z+8, the rope force data that intact sensing rope therein is only selected during monitoring structural health conditions (is referred to as the monitored amount that can actually use, its quantity is recorded for K, K cannot be less than Z+4) and the monitored amount unit change matrix Δ C progress health state evaluation of corresponding Cable Structure.This M newly increased will be monitored during monitoring structural health conditions2Root senses the Suo Li of rope.The M newly increased2Root sensing rope should be arranged on structure on, personnel be easy to reach position, be easy to personnel to carry out Non-Destructive Testing to it.
The M newly increased in the method2Root senses a part for Suo Zuowei Cable Structures, when mentioning Cable Structure again hereinafter, and Cable Structure includes increase M2Cable Structure and the M newly increased before root sensing rope2Root senses rope, that is to say, that refer to the M including newly increasing when mentioning Cable Structure hereinafter2Root senses the Cable Structure of rope.Therefore mention hereinafter when obtaining " Cable Structure steady temperature data " according to the temperature survey of the Cable Structure of this method " calculate method " survey calculation, Cable Structure therein includes the M newly increased2Root senses rope, and obtained " Cable Structure steady temperature data " include the M newly increased2Root senses the steady temperature data of rope, obtains the M newly increased2The method of the steady temperature data of root sensing rope is same as the M of Cable Structure1The preparation method of the steady temperature data of root support cable, is no longer handed over one by one later;Measure the M newly increased2The Suo Li of root sensing rope method is same as the M of Cable Structure1The Suo Li of root support cable measuring method, is no longer handed over one by one later;When carrying out any measurement to the support cable of Cable Structure, while to the M newly increased2The measurements of root sensing Suo Jinhang equally, are no longer handed over one by one later;The M newly increased2Root senses rope in addition to not damaging and relaxing, the M newly increased2The information content of root rope is identical with the information content of the support cable of Cable Structure, no longer hands over one by one later;The M newly increased2The Suo Li of root sensing rope is exactly increased M2Other individual monitored amounts.When setting up the various mechanical models of Cable Structure later, by the M newly increased2Root sensing rope treats as the M of Cable Structure1Root support cable is treated, the occasion of damage and relaxation except mentioning support cable, the M including newly increasing when other occasions mention support cable2Root rope.
The monitored amount of summary, whole Cable Structure has M (M=M1+M2) individual monitored amount, M adds 4 not less than P.
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.
According to the method provided in technical scheme and claims, using the conventional method including ordinary temperature survey calculation method, it is determined that " temperature survey of the Cable Structure of this method calculates method ".
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, 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).
T is obtained in actual measurementoWhile, the Actual measurement data of Cable Structure are obtained using conventional method Actual measurement.The Actual measurement data of Cable Structure include the data that Non-destructive Testing Data of support cable etc. can express the health status of rope, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinates data (initial Cable Structure bearing generalized coordinates data include initial Cable Structure bearing spatial data and initial Cable Structure bearing angular data), the initial linear movement measuring data of Cable Structure bearing, Cable Structure concentrfated load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structural point measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing spatial data refers to the bearing spatial data under Cable Structure design point, and the initial linear movement measuring data of Cable Structure bearing, which refer to, is setting up initial mechanical calculating benchmark model AoWhen, Cable Structure bearing is relative to the displacement of the lines that the bearing under Cable Structure design point occurs.The initial geometric data of Cable Structure can be spatial data of the spatial data plus a series of point in structure of the end points of all ropes, it is therefore intended that the geometric properties of Cable Structure are determined according to these coordinate datas.For cable-stayed bridge, if initial geometric data, which can be the spatial data of the end points of all ropes, adds the spatial data done on bridge two ends, here it is so-called bridge type data.Non-destructive Testing Data using support cable etc. can express the initial linear movement measuring data of the data of the health status of support cable, Cable Structure bearing and Cable Structure load measurement data set up evaluation object initial damage vector do, use doRepresent Cable Structure (with initial mechanical calculating benchmark model AoRepresent) evaluation object initial health.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;If without the initial linear movement measuring data of Cable Structure bearing or can consider the initial displacement of the lines of Cable Structure bearing be 0 when, vectorial doIn each element numerical value related to Cable Structure bearing displacement of the lines take 0;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.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。
No matter which kind of method to obtain initial mechanical calculating benchmark model A witho, it is included in " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector To), based on AoCalculate obtained Cable Structure calculate data must closely its measured data, error typically cannot be greater than 5%.So can utility AoSuo Li under analog case obtained by calculating calculates that data, strain calculation data, Cable Structure shapometer count evidence and displacement meter counts evidence, Cable Structure angle-data, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model AoThe health status of middle support cable evaluation object initial damage vector doRepresent, Cable Structure steady temperature data are with initial Cable Structure steady temperature data vector ToRepresent.Due to based on AoThe initial value (actual measurement is obtained) for obtaining the evaluations of all monitored amounts closely all monitored amounts is calculated, so A can also be used inoOn the basis of, carry out Mechanics Calculation obtain, AoEach monitored amount the monitored amount initial value vector C of evaluation compositiono.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.ToAnd doIt is AoParameter, CoBy AoMechanics Calculation result composition.
3rd step:In the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated;Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingi o, AoAnd Ai oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated;When ith circulation starts, corresponding to Ai o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Ti oRepresent, vector Ti oDefinition mode and vector ToDefinition mode it is identical, Ti oElement and ToElement correspond;The current initial damage vector of evaluation object that ith circulation needs when starting is designated as di o, di oRepresent Cable Structure A during this time circulation beginningi oEvaluation object health status, di oDefinition mode and doDefinition mode it is identical, di oElement and doElement correspond;When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amounti oRepresent, vectorial Ci oDefinition mode and vector CoDefinition mode it is identical, Ci oElement and CoElement correspond, be monitored the current initial value vector C of amounti oRepresent to correspond to Ai oAll monitored amounts concrete numerical value;Ti oAnd di oIt is Ai oCharacterisitic parameter;Ci oBy Ai oMechanics Calculation result composition;When circulation starts for the first time, Ai oIt is designated as A1 o, set up A1 oMethod to make A1 oEqual to Ao;When circulation starts for the first time, Ti oIt is designated as T1 o, set up T1 oMethod to make T1 oEqual to To;When circulation starts for the first time, di oIt is designated as d1 o, set up d1 oMethod to make d1 oEqual to do;When circulation starts for the first time, Ci oIt is designated as C1 o, set up C1 oMethod to make C1 oEqual to Co。
4th step:The hardware components of cable structure health monitoring system are installed.Hardware components at least include:Monitored amount monitoring system (such as system containing cable force measurement, signal conditioner), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, computer and communication alert equipment.Each monitored amount, each temperature must be monitored system monitoring and arrive, and the signal monitored is transferred to signal (data) collector by monitoring system;Signal is delivered to computer through signal picker;Computer is then responsible for the health monitoring software of the evaluation object of operation Cable Structure, including the signal that the transmission of tracer signal collector comes;When monitoring that evaluation object health status is changed, computer control communication warning device is alarmed to monitoring personnel, owner and (or) the personnel specified.
5th step:Establishment and the on computers system software of installation and operation this method, the software will complete the functions (all work that can be completed with computer i.e. in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, alarm that this method required by task is wanted.
6th step:Thus step starts the cycle over running, during structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Ti, vector TiDefinition mode and vector ToDefinition mode it is identical, TiElement and ToElement correspond;In actual measurement vector TiWhile, actual measurement obtains the currency of all monitored amounts in Cable Structure, all these monitored amount current value vector C of numerical value compositioni, vectorial CiDefinition mode and vector CoDefinition mode it is identical, CiElement and CoElement correspond, represent identical monitored amount in numerical value not in the same time.
Current Cable Structure steady temperature data vector T is obtained in actual measurementiWhile, to the M newly increased2Root senses Suo Jinhang Non-Destructive Testings, such as ultrasonic examination, visual inspection, infrared imaging is checked, therefrom identify the sensing rope for occurring damaging or relaxing, according to monitored amount coding rule, appearances that removes with identifying is damaged or the corresponding element of loose sensing rope in each vector numbered according to monitored amount coding rule occurred before this method, also no longer occur in each vector sum matrix occurred after this method the appearance with identifying damage or relaxation the corresponding element of sensing rope, no longer include being identified the sensing rope for occurring damaging or relaxing here when mentioning sensing rope after this method, no longer include the Suo Li for being identified the sensing rope for occurring damaging or relaxing here when mentioning monitored amount after this method;Several sensing ropes for occurring damaging or relaxing are identified from Cable Structure, just by M2Reduce same quantity with M.
7th step:Obtaining current Cable Structure steady temperature data vector TiAfterwards, T is comparediAnd Ti oIf, TiEqual to Ti o, then need not be to Ai o、Ci oAnd Ti oIt is updated, otherwise needs to current initial mechanical calculating benchmark model Ai o, current initial Cable Structure steady temperature data vector Ti oWith monitored amount current initial value vector Ci oIt is updated, and the current initial damage vector d of evaluation objecti oKeep constant, update method is carried out as the step of being provided in technical scheme and claims.
8th step:In current initial mechanical calculating benchmark model Ai oOn the basis of, according to Mechanics Calculation several times is carried out the step of being provided in technical scheme and claims, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingiWith evaluation object unit change vector Di u, wherein can take 5%, 10%, 20% or 30% equivalent damage be unit damage, if the evaluation object is the translational component in a direction of a bearing, it is assumed that the bearing in the direction of displacement in vectorial di oOccurs unit displacement of the lines (such as 2mm, 5mm, 10mm etc. are unit displacement of the lines) on the basis of the existing displacement of the lines of the bearing represented again, if the evaluation object is a load, it is assumed that the load is in vectorial di oLoad unit is further added by the basis of the existing variable quantity of the load represented to change (if the load is distributed load, and the distributed load is line distributed load, load unit change can take 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. to change for unit;If the load is distributed load, and the distributed load is EDS maps load, and load unit change can take 1MPa, 2MPa, 3MPa or 1kNm/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), use Di ukThis unit damage or unit displacement of the lines or load unit change are recorded, wherein k represents the numbering for occurring unit damage or unit displacement of the lines or the evaluation object of load unit change, Di ukIt is evaluation object unit change vector Di uAn element, evaluation object unit change vector Di uElement coding rule and vector doElement coding rule it is identical.
9th step:Set up linear relationship error vector eiWith vectorial gi.Utilize data (" the monitored current initial value vector C of amount abovei o", " unit damage monitored numerical quantity transformation matrices Δ Ci"); while the 8th step is calculated each time; i.e. while the increase unit damage or unit displacement of the lines of assuming only one of which evaluation object in evaluation object or load unit change is calculated each time; when assuming that kth (k=1,2,3; ...; N) individual evaluation object increase unit damage or unit displacement of the lines or during load unit change, one injury vector of composition is calculated each time, d is usedi tkThe injury vector is represented, corresponding monitored amount calculation current vector is Ci tk(referring to the 8th step), injury vector di tkElement number be equal to evaluation object quantity, vectorial di tkAll elements in only one of which element numerical value take each time calculate in assume increase unit damage unit displacement of the lines or load unit change evaluation object unit damage or unit displacement of the lines or load unit changing value, di tkThe numerical value of other elements take 0, that for 0 element evaluation object of the numbering with assuming increase unit damage or unit displacement of the lines or load unit change corresponding relation, with the elements of the same numberings of other vectors with the corresponding relation of the evaluation object be identical;di tkWith evaluation object initial damage vector doElement number rule it is identical, di tkElement and doElement be one-to-one relationship.By Ci tk、Ci o、ΔCi、di tkBring formula (1) into, obtain a linear relationship error vector ei k, calculate obtain a linear relationship error vector e each timei k;ei kSubscript k represent the individual evaluation object increase unit damage of kth (k=1,2,3 ..., N) or unit displacement of the lines or load unit change.There is N number of evaluation object just to have n times calculating, just there is N number of linear relationship error vector ei k, by this N number of linear relationship error vector ei kA vector is obtained after addition, is exactly final linear relationship error vector e by the new vector obtained after each element divided by N of this vectori.Vectorial giEqual to final error vector ei.By vectorial giOn the hard disc of computer for being stored in operation health monitoring systems software, used for health monitoring systems software.
Tenth step:Define current nominal fatigue vector di cWith currently practical injury vector di, di cAnd diElement number be equal to evaluation object quantity, di cAnd diElement and evaluation object between be one-to-one relationship, di cAnd diElement numerical value represent correspondence evaluation object degree of injury or displacement of the lines amount or load change degree, di cAnd diWith evaluation object initial damage vector doElement number rule it is identical, di cElement, diElement and doElement be one-to-one relationship.
11st step:According to monitored amount current value vector CiWith " the monitored current initial value vector C of amounti o", " unit damage monitored numerical quantity transformation matrices Δ Ci" and " current nominal fatigue vector di c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as formula (2), current nominal fatigue vector d calculated according to multi-objective optimization algorithmi cNoninferior solution, 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 nominal fatigue vector di c, the specific programming realization of Objective Programming had general program directly to use.
12nd step:According to the currently practical injury vector d of cable systemiDefinition and its element definition calculate obtain currently practical injury vector diEach element, so as to by diDetermine the health status of evaluation object.Currently practical injury vector diK-th of element di kRepresent the currently practical health status of k-th of evaluation object in ith circulation.
di kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then di kRepresent its currently practical damage, di kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity.
di kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a translational component of a bearing, then di kRepresent its currently practical displacement of the lines numerical value.
13rd step:Computer in health monitoring systems is periodically automatic or generates cable system health condition form by human users' health monitoring systems.
14th step:Under specified requirements, the computer in health monitoring systems is automatically brought into operation communication alert equipment and alarmed to monitoring personnel, owner and (or) the personnel specified.
15th step:Set up mark vector BiIf, mark vector BiElement be all 0, then return to the 6th step and proceed health monitoring and calculating to cable system;If mark vector BiElement be not all 0, then complete after subsequent step, into circulating next time.
16th step:Calculating obtains initial damage vector d next time needed for (i.e. i+1 time, i=1,2,3,4 ...) circulationi+1 oEach element di+1 ok(k=1,2,3 ..., N);In initial mechanical calculating benchmark model AoOn the basis of, 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, then it is d to make the health status of ropei+1 oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time (i=1,2,3,4 ...) circulationi+1;Current initial Cable Structure steady temperature data vector T needed for (i.e. i+1 time, i=1,2,3,4 ...) is circulated next timei+1 oEqual to Ti o.Obtain Ai+1、di+1 oAnd Ti+1 oAfterwards, A is obtained by Mechanics Calculationi+1In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationi+1 o。
17th step:The 6th step is returned to, starts the circulation by the 6th step to the 17th step.
Claims (1)
1. cable force monitoring damaged cable load load/linear displacement progressive recognition method, it is characterised in that methods described includes:
A. when though the load that Cable Structure is born is changed, during initial without departing from the Cable Structure allowable load of the load that Cable Structure is being born, this method is applicable;The initial allowable load of Cable Structure refers to allowable load of the Cable Structure in completion, can be obtained by conventional Mechanics Calculation;This method unitedly calls evaluated bearing translational component, support cable and load to be evaluation object, if the quantity that the quantity sum of the evaluated quantity of bearing translational 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 bearing translational component in " evaluation object " with title " core evaluation object ", if evaluated support cable and the quantity sum of bearing translational 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 ";If having M in cable system1Root support cable, Cable Structure rope force data includes this M1The Suo Li of root support cable, this method is in monitoring whole M1On the basis of root supporting cable force, the artificial increase M in Cable Structure2Root rope, referred to as senses rope, and this M newly increased will be monitored during cable structure health monitoring2Root senses the Suo Li of rope;The monitored amount of summary, the M Suo Li that whole Cable Structure has M root ropes is monitored, that is, has M monitored amounts, and wherein M is M1With M2Sum;M have to be larger than the quantity of core evaluation object, and M is less than the quantity of evaluation object;The M newly increased2The rigidity of root sensing rope, should be much smaller compared with the rigidity of any one support cable of Cable Structure;The M newly increased2The Suo Li of each sensing rope of root sensing rope should be more much smaller than the Suo Li of any one support cable of Cable Structure, even if can so ensure this M newly increased2Root sensing rope occurs in that damage or relaxation, and the influence to the stress of other components of Cable Structure, strain, deformation is very little;The M newly increased2Direct stress should be less than its fatigue limit on the cross section of root sensing rope, and these requirements can ensure the M newly increased2Fatigue damage occurs for root sensing Suo Buhui;The M newly increased2The two ends of root sensing rope should fully be anchored, it is ensured that be not in relaxation;The M newly increased2Root sensing rope should obtain sufficient anti-corrosion protection, it is ensured that the M newly increased2Damage and relaxation occur for root sensing Suo Buhui;For convenience, " monitored all parameters of Cable Structure " are referred to as in the method " monitored amount ";To M monitored amount serial numbers, this method is with this numbering, j=1,2,3 ..., M is represented with variable j, and the numbering will be used to generate vector sum matrix in subsequent step;The M newly increased in the method2Root senses a part for Suo Zuowei Cable Structures, when mentioning Cable Structure again hereinafter, and Cable Structure includes increase M2Cable Structure and the M newly increased before root sensing rope2Root senses rope, that is to say, that refer to the M including newly increasing when mentioning Cable Structure hereinafter2Root senses the Cable Structure of rope;Therefore mention hereinafter when obtaining " Cable Structure steady temperature data " according to the temperature survey of the Cable Structure of this method " calculate method " survey calculation, Cable Structure therein includes the M newly increased2Root senses rope, and obtained " Cable Structure steady temperature data " include the M newly increased2Root senses the steady temperature data of rope, obtains the M newly increased2The method of the steady temperature data of root sensing rope is same as the M of Cable Structure1The preparation method of the steady temperature data of root support cable, is no longer handed over one by one later;Measure the M newly increased2The Suo Li of root sensing rope method is same as the M of Cable Structure1The Suo Li of root support cable measuring method, is no longer handed over one by one later;When carrying out any measurement to the support cable of Cable Structure, while to the M newly increased2The measurements of root sensing Suo Jinhang equally, are no longer handed over one by one later;The M newly increased2Root senses rope in addition to not damaging and relaxing, to the M newly increased2Requirement and preparation method of the requirement and preparation method of the information content of root sensing rope with the information content of the support cable of Cable Structure are identical, no longer hand over one by one later;When setting up the various mechanical models of Cable Structure later, by the M newly increased2The support cable that root sensing rope treats as Cable Structure is treated;Below, in addition to the occasion of damage and relaxation for mentioning support cable, the M that described support cable includes the support cable of Cable Structure and newly increased when mentioning support cable2Root senses rope;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;T is obtained in actual measurementoWhile, 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 linear movement measuring 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;Using the initial linear movement measuring data of support cable initial health data, Cable Structure bearing can be expressed and Cable Structure load measurement data set up evaluation object initial damage vector do, 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 translational component, then doThe element numerical value represent this bearing the translational 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 linear movement measuring data of Cable Structure bearing or can consider the initial displacement of the lines of Cable Structure bearing be 0 when, vectorial doIn each element numerical value related to Cable Structure bearing displacement of the lines 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 spatial data refers to the bearing spatial data under Cable Structure design point, and the initial linear movement measuring data of Cable Structure bearing, which refer to, is setting up initial mechanical calculating benchmark model AoWhen, Cable Structure bearing is relative to the displacement of the lines 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 linear movement measuring 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;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, A in the methodo、Co、doAnd ToIt is constant;
E. in the method, alphabetical i is in addition to the place for clearly indicating that number of steps, and alphabetical i only represents that cycle-index, i.e. ith are circulated;Ith circulation needs the current initial mechanical calculating benchmark model of Cable Structure setting up or having set up to be designated as current initial mechanical calculating benchmark model A when startingi o, AoAnd Ai oTemperature parameter has been included in, Effect on Mechanical Properties of the temperature change to Cable Structure can be calculated;When ith circulation starts, corresponding to Ai o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector Ti oRepresent, vector Ti oDefinition mode and vector ToDefinition mode it is identical, Ti oElement and ToElement correspond;The current initial damage vector of evaluation object that ith circulation needs when starting is designated as di o, di oRepresent Cable Structure A during this time circulation beginningi oEvaluation object health status, di oDefinition mode and doDefinition mode it is identical, di oElement and doElement correspond;When ith circulation starts, the initial value of all monitored amounts, with the current initial value vector C of monitored amounti oRepresent, vectorial Ci oDefinition mode and vector CoDefinition mode it is identical, Ci oElement and CoElement correspond, be monitored the current initial value vector C of amounti oRepresent to correspond to Ai oAll monitored amounts concrete numerical value;Ti oAnd di oIt is Ai oCharacterisitic parameter, Ci oBy Ai oMechanics Calculation result composition;When circulation starts for the first time, Ai oIt is designated as A1 o, set up A1 oMethod to make A1 oEqual to Ao;When circulation starts for the first time, Ti oIt is designated as T1 o, set up T1 oMethod to make T1 oEqual to To;When circulation starts for the first time, di oIt is designated as d1 o, set up d1 oMethod to make d1 oEqual to do;When circulation starts for the first time, Ci oIt is designated as C1 o, set up C1 oMethod to make C1 oEqual to Co;
F. enter from here and the circulation walked to q is walked by f;During structure military service, the current data of Cable Structure steady temperature data is obtained according to " temperature survey of the Cable Structure of this method calculates method " constantly Actual measurement, the current data for owning " Cable Structure steady temperature data " constitutes current Cable Structure steady temperature data vector Ti, vector TiDefinition mode and vector ToDefinition mode it is identical, TiElement and ToElement correspond;Vector T is obtained in actual measurementiWhile, actual measurement obtains obtaining current Cable Structure steady temperature data vector TiAt the time of synchronization Cable Structure in all monitored amounts currency, the monitored amount current value vector C of all these numerical value compositioni, vectorial CiDefinition mode and vector CoDefinition mode it is identical, CiElement and CoElement correspond, represent identical monitored amount in numerical value not in the same time;Current Cable Structure steady temperature data vector T is obtained in actual measurementiWhile, to the M newly increased2Root senses Suo Jinhang Non-Destructive Testings, therefrom identify the sensing rope for occurring damaging or relaxing, according to monitored amount coding rule, appearances that removes with identifying is damaged or the corresponding element of loose sensing rope in each vector numbered according to monitored amount coding rule occurred before this method, also no longer occur in each vector sum matrix occurred after this method the appearance with identifying damage or relaxation the corresponding element of sensing rope, no longer include being identified the sensing rope for occurring damaging or relaxing here when mentioning sensing rope after this method, no longer include the Suo Li for being identified the sensing rope for occurring damaging or relaxing here when mentioning monitored amount after this method;Several sensing ropes for occurring damaging or relaxing are identified from Cable Structure, just by M2Reduce same quantity with M;
G. according to current Cable Structure steady temperature data vector Ti, current initial mechanical calculating benchmark model A is updated according to step g1 to g3i o, the monitored current initial value vector C of amounti oWith current initial Cable Structure steady temperature data vector Ti o, and the current initial damage vector d of evaluation objecti oKeep constant;
G1. T is comparediWith Ti oIf, TiEqual to Ti o, then Ai o、Ci oAnd Ti oKeep constant;Otherwise need to follow these steps to Ai o、Ci oAnd Ti oIt is updated;
G2. T is calculatediWith ToDifference, TiWith ToDifference be exactly change of the current Cable Structure steady temperature data on initial Cable Structure steady temperature data, TiWith ToDifference represented with steady temperature change vector S, S be equal to TiSubtract To, S represents the change of Cable Structure steady temperature data;
G3. 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 updatedi o, update Ai oWhile, Ti oAll elements numerical value also uses TiAll elements numerical value correspondence replace, that is, have updated Ti o, thus obtained properly corresponding to Ai oTi o;Now di oKeep constant;As renewal Ai oAfterwards, Ai oRope the health status current initial damage vector d of evaluation objecti oRepresent, Ai oCable Structure steady temperature with current Cable Structure steady temperature data vector TiRepresent, update Ci oMethod be:As renewal Ai oAfterwards, A is obtained by Mechanics Calculationi oIn all monitored amounts, current concrete numerical value, these concrete numerical values composition Ci o;
H. in current initial mechanical calculating benchmark model Ai oOn the basis of, Mechanics Calculation several times is carried out according to step h1 to step h4, unit damage monitored numerical quantity transformation matrices Δ C is set up by calculatingiWith evaluation object unit change vector Di u;
H1. when ith circulates beginning, method obtains Δ C directly as listed by step h2 to step h4iAnd Di u;At other moment, when in step g to Ai oAfter being updated, it is necessary to which the method as listed by step h2 to step h4 regains Δ CiAnd Di uIf, not to A in step gi oIt is updated, then is directly transferred to step i here and carries out follow-up work;
H2. in current initial mechanical calculating benchmark model Ai oOn the basis of carry out Mechanics Calculation several times, calculation times are numerically equal to the quantity N of all evaluation objects, have it is N number of assessment object just have n times calculating;According to the coding rule of evaluation object, calculated successively;Calculate each time and assume that only one of which evaluation object is further added by unit damage or load unit change on the basis of original damage or load, specifically, if the evaluation object is a support cable in cable system, it is assumed that the support cable is further added by unit damage, if the evaluation object is the translational component in a direction of a bearing, it is assumed that the bearing is further added by unit displacement of the lines in the direction of displacement, if the evaluation object is a load, it is assumed that the load is further added by load unit change, D is usedi ukThis increased unit damage or unit displacement of the lines or load unit change is recorded, wherein k represents the numbering for increasing unit damage or unit displacement of the lines or the evaluation object of load unit change, Di ukIt is evaluation object unit change vector Di uAn element, evaluation object unit change vector Di uElement coding rule and vector doElement coding rule it is identical;The evaluation object that unit damage or unit displacement of the lines or load unit change are further added by calculating each time is different from the evaluation object that unit damage or unit displacement of the lines or load unit change are further added by other calculating, the current calculated value for all monitored amounts that Cable Structure is all calculated using mechanics method is calculated each time, and the current calculated value that obtained all monitored amounts are calculated each time constitutes a monitored amount calculation current vector;When assuming that k-th of evaluation object is further added by unit damage or unit displacement of the lines or load unit change, C is usedi tkRepresent corresponding " monitored amount calculation current vector ";When in this step to each vectorial element number, same coding rule should be used with other vectors in this method, to ensure any one element in this step in each vector, with other vectors, numbering identical element, same monitored amount or the relevant information of same target are expressed;Ci tkDefinition mode and vector CoDefinition mode it is identical, Ci tkElement and CoElement correspond;
H3. obtained vectorial C is calculated each timei tkSubtract vectorial Ci oA vector is obtained, then " numerical value change vector δ a C for monitored amount will be obtained after each element of the vector divided by the assumed unit damage of this calculating or unit displacement of the lines or load unit change numerical valuei k”;There is N number of evaluation object just to have N number of " the numerical value change vector of monitored amount ";
H4. " the unit damage monitored numerical quantity transformation matrices Δ C for having N to arrange is constituted successively according to the coding rule of N number of evaluation object by this N number of " numerical value change vector of monitored amount "i”;Unit damage monitored numerical quantity transformation matrices Δ CiEach row correspond to a monitored amount unit change vector;Unit damage monitored numerical quantity transformation matrices Δ CiEvery a line correspond to different unit change amplitudes of the same monitored amount when different evaluation objects increase unit damage or unit displacement of the lines or load unit change;Unit damage monitored numerical quantity transformation matrices Δ CiRow coding rule and vector doElement coding rule it is identical, unit damage monitored numerical quantity transformation matrices Δ CiRow coding rule and M monitored amounts coding rules it is identical;
I. current nominal fatigue vector d is definedi cWith currently practical injury vector di, di cAnd diElement number be equal to evaluation object quantity, di cAnd diElement and evaluation object between be one-to-one relationship, di cElement numerical value represent the nominal fatigue degree or nominal displacement of the lines or nominal load variable quantity of correspondence evaluation object, di cAnd diWith evaluation object initial damage vector doElement number rule it is identical, di cElement, diElement and doElement be one-to-one relationship;
J. according to monitored amount current value vector CiWith " the monitored current initial value vector C of amounti o", " unit damage monitored numerical quantity transformation matrices Δ Ci" and " current nominal fatigue vector di c" between the linear approximate relationship that exists, the linear approximate relationship can be expressed as removing d in formula 1, formula 1i cOuter other amounts are, it is known that solution formula 1 can just calculate current nominal fatigue vector di c;
K. the currently practical injury vector d expressed using formula 2iK-th of element di kWith the current initial damage vector d of evaluation objecti oK-th of element di okWith current nominal fatigue vector di cK-th of element di ckBetween relation, calculating obtain currently practical injury vector diAll elements;
Formula 2
K=1,2,3 in formula 2 ..., N;di kThe currently practical health status of k-th of evaluation object in ith circulation is represented, if the evaluation object is a support cable in cable system, then di kRepresent its currently practical damage, di kFor 0 when represent not damaged, represent that the support cable thoroughly loses bearing capacity when being 100%, represented when between 0 and 100% lose corresponding proportion bearing capacity;If the evaluation object is a translational component of a bearing, then di kRepresent its currently practical displacement of the lines numerical value;So far this method realizes influence, Cable Structure the damaged cable identification for rejecting bearing displacement of the lines, load change and structure temperature change, and load change, structure temperature change and the change influence of support cable health status, bearing displacement of the lines identification are rejected while realizing;
L. current nominal fatigue vector d is tried to achievei cAfterwards, mark vector B is set up according to formula 3i, formula 4 gives mark vector BiK-th of element definition;
Formula 4
Element B in formula 4i kIt is mark vector BiK-th of element, Di ukIt is evaluation object unit change vector Di uK-th of element, di ckIt is the current nominal fatigue vector d of evaluation objecti cK-th of element, they all represent k=1,2,3 ... ..., N in the relevant information of k-th of evaluation object, formula 4;
If m. mark vector BiElement be all 0, then return to step f continue this circulation;If mark vector BiElement be not all 0, then enter next step, i.e. step n;
N. the current initial damage vector d of evaluation object obtained next time, i.e. needed for i+1 time circulation is calculated according to formula 5i+1 oEach element;
Formula 5
D in formula 5i+1 okIt is the current initial damage vector d of evaluation object next time, i.e. needed for i+1 time circulationi+1 oK-th of element, di okThis, i.e. the current initial damage vector d of evaluation object of ith circulationi oK-th of element, Di ukIt is the evaluation object unit change vector D of ith circulationi uK-th of element, Bi kIt is the mark vector B of ith circulationiK-th of element, k=1,2,3 ... ..., N in formula 5;
O. remove once, i.e. i+1 time circulates required current initial Cable Structure steady temperature data vector Ti+1 oThe current initial Cable Structure steady temperature data vector T circulated equal to ithi o;
P. in initial mechanical calculating benchmark model AoOn the basis of, 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, then it is d to make the health status of ropei+1 oThat obtain afterwards is exactly Mechanics Calculation benchmark model A next time, i.e. needed for i+1 time circulationi+1;Obtain Ai+1Afterwards, A is obtained by Mechanics Calculationi+1In all monitored amounts, current concrete numerical value, these concrete numerical values constitute next time, i.e. the required current initial value vector C of monitored amount of i+1 time circulationi+1 o;
Q. step f is returned to, starts to circulate next time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510439492.9A CN105067328A (en) | 2015-07-23 | 2015-07-23 | Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510439492.9A CN105067328A (en) | 2015-07-23 | 2015-07-23 | Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105067328A true CN105067328A (en) | 2015-11-18 |
Family
ID=54496752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510439492.9A Pending CN105067328A (en) | 2015-07-23 | 2015-07-23 | Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105067328A (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103604656A (en) * | 2013-12-09 | 2014-02-26 | 东南大学 | Damaged cable/concentrated load/linear displacement progressive recognition method on basis of cable power monitoring |
CN103604636A (en) * | 2013-12-09 | 2014-02-26 | 东南大学 | Linear displacement progressive type identification method for cable force monitoring of damaged cable and intensive load |
CN103852319A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Method for progressively recognizing damaged cables, loads and linear displacement based on cable force monitoring |
CN103852309A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Progressive recognition method for problem cable load linear displacement based on mixing monitoring |
CN103852323A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Method for progressively recognizing damaged cables, loads and linear displacement based on cable force monitoring |
CN103868736A (en) * | 2014-03-10 | 2014-06-18 | 东南大学 | Linear displacement cable force monitoring-based progressive load identification method for damaged cable |
CN103868714A (en) * | 2014-03-10 | 2014-06-18 | 东南大学 | Linear displacement strain monitoring-based progressive load identification method for damaged cable |
-
2015
- 2015-07-23 CN CN201510439492.9A patent/CN105067328A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103604656A (en) * | 2013-12-09 | 2014-02-26 | 东南大学 | Damaged cable/concentrated load/linear displacement progressive recognition method on basis of cable power monitoring |
CN103604636A (en) * | 2013-12-09 | 2014-02-26 | 东南大学 | Linear displacement progressive type identification method for cable force monitoring of damaged cable and intensive load |
CN103852319A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Method for progressively recognizing damaged cables, loads and linear displacement based on cable force monitoring |
CN103852309A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Progressive recognition method for problem cable load linear displacement based on mixing monitoring |
CN103852323A (en) * | 2014-03-10 | 2014-06-11 | 东南大学 | Method for progressively recognizing damaged cables, loads and linear displacement based on cable force monitoring |
CN103868736A (en) * | 2014-03-10 | 2014-06-18 | 东南大学 | Linear displacement cable force monitoring-based progressive load identification method for damaged cable |
CN103868714A (en) * | 2014-03-10 | 2014-06-18 | 东南大学 | Linear displacement strain monitoring-based progressive load identification method for damaged cable |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105136497A (en) | Simplified cable force monitoring damaged cable load generalized displacement progressive identification method | |
CN105067300A (en) | Simplified cable-force-monitoring defective cable load generalized displacement progressive identifying method | |
CN105067033A (en) | Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined linear displacement | |
CN105067328A (en) | Method for progressively recognizing load linear displacement of damaged cable based on cable force monitoring process | |
CN105067032A (en) | Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined linear displacement | |
CN105067314A (en) | Method for progressively recognizing load of damaged cable based on strain monitoring process of streamlined annular displacement | |
CN105067308A (en) | Simplified angular displacement cable-force-monitoring defective cable load progressive identifying method | |
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 | |
CN105067354A (en) | Method for progressively recognizing load of damaged cable based on streamlined hybrid monitoring process of generalized displacement | |
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 | |
CN105067347A (en) | Method for progressively recognizing load linear displacement of damaged cable based on streamlined hybrid monitoring process | |
CN105136499A (en) | Simplified space coordinate monitoring faulty cable load progressive identification method | |
CN105136498A (en) | Faulty cable load generalized displacement progressive identification method based on angle monitoring | |
CN105067356A (en) | Method for recognizing load of damaged cable based on streamlined hybrid monitoring process | |
CN105067313A (en) | Method for recognizing load of damaged cable based on angle monitoring process of streamlined angular displacement | |
CN105067318A (en) | Method for recognizing load of damaged cable based on space coordinate monitoring process of streamlined linear displacement | |
CN105067332A (en) | Method for recognizing load linear displacement of damaged cable based on streamlined angle monitoring process | |
CN105067335A (en) | Method for recognizing load generalized displacement of damaged cable based on streamlined space coordinate monitoring process | |
CN105067305A (en) | Simplified generalized displacement strain-monitoring damaged cable load identifying method | |
CN105067320A (en) | Method for recognizing load of damaged cable based on angle monitoring process of streamlined and generalized displacement | |
CN105067343A (en) | Method for recognizing load angular displacement of damaged cable based on streamlined hybrid monitoring process | |
CN105067342A (en) | Method for recognizing load generalized displacement of damaged cable based on streamlined hybrid monitoring process | |
CN105136491A (en) | Simplified strain monitoring faulty cable load identification method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151118 |
|
RJ01 | Rejection of invention patent application after publication |