CN105136491A - Simplified strain monitoring faulty cable load identification method - Google Patents

Simplified strain monitoring faulty cable load identification method Download PDF

Info

Publication number
CN105136491A
CN105136491A CN201510437882.2A CN201510437882A CN105136491A CN 105136491 A CN105136491 A CN 105136491A CN 201510437882 A CN201510437882 A CN 201510437882A CN 105136491 A CN105136491 A CN 105136491A
Authority
CN
China
Prior art keywords
cable structure
temperature
data
vector
cable
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201510437882.2A
Other languages
Chinese (zh)
Inventor
韩玉林
韩佳邑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201510437882.2A priority Critical patent/CN105136491A/en
Publication of CN105136491A publication Critical patent/CN105136491A/en
Pending legal-status Critical Current

Links

Abstract

A simplified strain monitoring faulty cable load identification method, based on strain monitoring, includes the steps: monitoring the temperature of a cable structure and an ambient temperature to determine whether a mechanical calculating reference model of the cable structure needs to be updated, obtaining a mechanical calculating reference model, of the cable structure, that takes the temperature of the cable structure and the ambient temperature into consideration, and obtaining a unit damage monitored quantity value change matrix through calculation on the basis of the model. A non-inferior solution of the current nominal damage vector of an evaluated object is calculated according to an approximation linear relation between a monitored quantity current value vector and a monitored quantity current initial value vector as well as an approximation linear relation between the unit damage monitored quantity value change matrix and the to-be-solved current nominal damage vector of the evaluated object. In this way, the influence of interference factors can be removed and faulty cables can be identified if the linear displacement of the bearer exists and the temperate changes.

Description

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

Claims (1)

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

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510437882.2A CN105136491A (en) 2015-07-23 2015-07-23 Simplified strain monitoring faulty cable load identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510437882.2A CN105136491A (en) 2015-07-23 2015-07-23 Simplified strain monitoring faulty cable load identification method

Publications (1)

Publication Number Publication Date
CN105136491A true CN105136491A (en) 2015-12-09

Family

ID=54721921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510437882.2A Pending CN105136491A (en) 2015-07-23 2015-07-23 Simplified strain monitoring faulty cable load identification method

Country Status (1)

Country Link
CN (1) CN105136491A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604649A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring
CN103604638A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load/linear displacement recognition method on basis of strain monitoring
CN103616243A (en) * 2013-12-09 2014-03-05 东南大学 Strain-monitoring damaged cable centralized load generalized displacement progressive recognition method
CN103616206A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring identification method for defective cable and concentrated loads
CN103616223A (en) * 2013-12-09 2014-03-05 东南大学 Method for identifying concentrated loads of damaged cables in linear displacement process based on strain monitoring
CN103616202A (en) * 2013-12-09 2014-03-05 东南大学 Angular displacement strain monitoring identification method for defective cable and concentrated loads
CN103852279A (en) * 2014-03-10 2014-06-11 东南大学 Strain monitoring error cable load identification method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604649A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring
CN103604638A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load/linear displacement recognition method on basis of strain monitoring
CN103616243A (en) * 2013-12-09 2014-03-05 东南大学 Strain-monitoring damaged cable centralized load generalized displacement progressive recognition method
CN103616206A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring identification method for defective cable and concentrated loads
CN103616223A (en) * 2013-12-09 2014-03-05 东南大学 Method for identifying concentrated loads of damaged cables in linear displacement process based on strain monitoring
CN103616202A (en) * 2013-12-09 2014-03-05 东南大学 Angular displacement strain monitoring identification method for defective cable and concentrated loads
CN103852279A (en) * 2014-03-10 2014-06-11 东南大学 Strain monitoring error cable load identification method

Similar Documents

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

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: 20151209

RJ01 Rejection of invention patent application after publication