CN102706659A - Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change - Google Patents

Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change Download PDF

Info

Publication number
CN102706659A
CN102706659A CN201210173398XA CN201210173398A CN102706659A CN 102706659 A CN102706659 A CN 102706659A CN 201210173398X A CN201210173398X A CN 201210173398XA CN 201210173398 A CN201210173398 A CN 201210173398A CN 102706659 A CN102706659 A CN 102706659A
Authority
CN
China
Prior art keywords
cable
temperature
data
vector
initial
Prior art date
Application number
CN201210173398XA
Other languages
Chinese (zh)
Other versions
CN102706659B (en
Inventor
韩玉林
王芳
韩佳邑
Original Assignee
东南大学
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 东南大学 filed Critical 东南大学
Priority to CN201210173398.XA priority Critical patent/CN102706659B/en
Publication of CN102706659A publication Critical patent/CN102706659A/en
Application granted granted Critical
Publication of CN102706659B publication Critical patent/CN102706659B/en

Links

Abstract

The invention discloses a defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change. The method comprises the following steps of: on the basis of angular monitoring, judging whether a mechanical calculation reference model of a cable structure is required to be updated by monitoring the temperature of the cable structure, environment temperature and the health state of an evaluated object, and thus obtaining the mechanical calculation reference model of the cable structure, into which the temperature of the cable structure and the environment temperature are recorded; and on the basis of the mechanical calculation reference model, acquiring a unit damage monitored amount unit change matrix by calculation. A non-inferior solution of the current nominal damage vector of the evaluated object is obtained by calculation according to an approximately linear relationship between the current value vector of monitored amount, the current initial value vector of the monitored amount, the unit damage monitored amount unit change matrix, a unit damage or unit angular displacement vector and the to-be-obtained current nominal damage vector of the evaluated object, so a damaged cable, a loosened cable and support angular displacement can be quickly identified when the temperature change occurs.

Description

The problem rope of temperature variation angle monitor and the laddering recognition methods of bearing angular displacement
Technical field
Structures such as cable-stayed bridge, suspension bridge, truss-frame structure have a common ground; Be exactly that they have many parts that bear tensile load; Like suspension cable, main push-towing rope, hoist cable, pull bar or the like; The common ground of this class formation is to be support unit with rope, cable or the rod member that only bears tensile load, and this method is " Cable Structure " with such structure representation for simplicity.Variation along with environment temperature; The temperature of Cable Structure also can change; When the Cable Structure temperature changes; (for example bearing is around the rotation of coordinate axis X, Y, Z in that the bearing angular displacement is arranged; In fact be exactly the angular displacement of bearing around coordinate axis X, Y, Z) time, this method is discerned damaged cable, slack line (truss-frame structure just is meant the impaired or lax rod member that only bears tensile load) and the bearing angular displacement in the supporting system of Cable Structure based on angle monitor, belong to the engineering structure health monitoring field.The supporting system of Cable Structure refers to all ropeway carrying-ropes, reaches all rod members that only bear tensile load that play supporting role; For simplicity; This patent is called " cable system " with whole support unit unifications of this class formation; But in fact cable system not only refers to support rope, also comprises the rod member that only bears tensile load, censures all ropeway carrying-ropes and the rod member that only bears tensile load that all play supporting role with " supporting rope " this noun in this method.This method is called the supporting rope of unsoundness problem with damaged cable and slack line, abbreviates the problem rope as.
Background technology
Cable Structure (particularly large-scale Cable Structure; For example large-scale cable-stayed bridge, suspension bridge) damage, structural health problem such as lax can appear in its supporting rope behind the long service; Structural health problems such as angular displacement can appear in its bearing, and these structural health problems all will cause the variation of structural internal force, and the safety of structure is caused harmful effect; Will cause the inefficacy of structure when serious, therefore accurately and timely discerning needs slack line, damaged cable and bearing angular displacement to be very important.
The supporting cable system health status change with bearing generation angular displacement after; Can cause the variation of the measurable parameter of structure; For example can cause the variation of Suo Li; Can influence the distortion or the strain of Cable Structure, can influence the shape or the volume coordinate of Cable Structure, can cause variation (the for example variation of the angle coordinate of the straight line of any this point of mistake in any any the section of body structure surface of angle coordinate of any imaginary line of the every bit of Cable Structure; The perhaps body structure surface variation of the angle coordinate of the normal of any arbitrarily); All these change the health status information that has all comprised cable system, in fact the variation of these measurable parameters comprised cable system health status information, comprised the bearing angular displacement information, that is to say that the measurable parameter that can utilize structure discerns bearing angular displacement and problem rope.This method is come identification problem rope and bearing angular displacement based on angle monitor (this method is called monitored angle " monitored amount ").Monitored amount is except the influence that receives the Cable Structure health status; Also can receive the influence of Cable Structure temperature variation (usually can take place); Under the condition that the Cable Structure temperature changes; If can realize based on monitoring to the supporting rope of unsoundness problem and the identification of bearing angular displacement to monitored amount, the safety of Cable Structure is had significant values, also there are not a kind of disclosed, effective health monitoring systems and method to solve this problem at present.
Summary of the invention
Technical matters: the purpose of this method is when Cable Structure has temperature variation, discloses a kind of based on health monitor method angle monitor, that can discern damaged cable in the Cable Structure, slack line and bearing angular displacement rationally and effectively.
Reason according to the Suo Li that supports rope changes can change the four kinds of situation that be divided into the Suo Li of supporting rope: the one, and the supporting rope has received damage, for example supports rope and localized cracks and corrosion or the like have occurred; The 2nd, supporting rope and not damaged; But variation has also taken place in Suo Li; The one of the main reasons that this variation occurs is that the Suo Changdu under the supporting rope free state (this moment, Suo Zhangli claimed that also Suo Li is 0) (is called drift; This method specially refers to support the drift of that section rope between rope two supporting end points) variation has taken place, one of purpose of this method will identify drift exactly the supporting rope that changes has taken place, and identifies the change amount of their drift; This change amount provides direct foundation for the Suo Li adjustment of this rope; The reason that supporting rope drift changes is not single, and for ease, this method is referred to as slack line with the supporting rope that drift changes; The 3rd, the Suo Li that Cable Structure bearing generation angular displacement causes changes; The 4th, the Suo Li that the Cable Structure temperature variation causes changes.
Technical scheme: this method is made up of three parts.Be respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; Two, set up required knowledge base of cable structure health monitoring system and parameter method, based on knowledge base (containing parameter) and the actual measurement monitored amount Cable Structure health status appraisal procedure; Three, the software and hardware part of health monitoring systems.
If the quantity sum of the quantity of the supporting rope of Cable Structure and the bearing angular displacement component of Cable Structure is N.For the purpose of narrating conveniently, this method is unified claims that support cable and the bearing angular displacement assessed are " evaluation object ", total N evaluation object.Give by the evaluation object serial number, this numbering will be used to generate the vector sum matrix in subsequent step.
If total M in the cable system 1Root supporting rope, structure rope force data comprises this M 1The Suo Li of root supporting rope, obviously M 1Less than by the quantity N of evaluation object.
" the whole monitored angle-data of structure " described by H angle coordinate component L the appointment straight line specified point of K on the structure, that cross each specified point, each appointment straight line, and the variation of structure angle is exactly the variation of the angle coordinate component of all appointments all specified points, all appointment straight lines.(individual angle coordinate component measurement value of M=K * L * H) or calculated value characterize the angle information of structure to each total M.K and M must not be less than N.
Comprehensive above-mentioned monitored amount, total M the monitored amount of whole Cable Structure, M must not be less than by the quantity N of evaluation object.
For simplicity, in the method " all monitored parameters of Cable Structure " are abbreviated as " monitored amount ".Give M monitored amount serial number, this numbering will be used to generate the vector sum matrix in subsequent step.This method is with representing this numbering with variable j, j=1, and 2,3 ..., M.
The first of this method: " the temperature survey calculating method of the Cable Structure of this method ".
At first confirm " the temperature survey calculating method of the Cable Structure of this method ".Because the temperature of Cable Structure possibly change; For example the temperature of the different parts of Cable Structure change along with the variation of intensity of sunshine, along with the variation of environment temperature changes; The surface of Cable Structure and temperature inside possibly be time dependent sometimes; The surface of Cable Structure possibly be different with temperature inside, and the surface of Cable Structure and temperature inside difference are time dependent, and the Mechanics Calculation of the Cable Structure when this just makes the account temperature condition is quite complicated with monitoring; For cost is measured in simplification problem, minimizing calculated amount and reduction; In order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", and is specific as follows especially:
The first step; Inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; Statistics obtains in T cloudy day 0 the highest temperature and the lowest temperature between back 30 minutes of the moment of sunrise next day at each cloudy day; Sunrise be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunrise constantly, the sunrise that can inquire about data or calculate each required day through conventional meteorology constantly, each cloudy day 0 up to sunrise next day constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature; T cloudy day arranged; The maximum temperature difference of daily temperature that T cloudy day just arranged, the maximal value of getting in the maximum temperature difference of daily temperature at T cloudy day is with reference to temperature difference per day, is designated as Δ T with reference to temperature difference per day rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m.On the surface of Cable Structure, get " R Cable Structure surface point "; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data ".When on the surface of Cable Structure, getting " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " is narrated with the condition that must satisfy that distributes in the back.From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is called reference plate; The one side of this reference plate on the sunny side; Be called sunny slope, the sunny slope of reference plate is coarse and dark color, the sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day; The non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate.Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly.
Second step; Monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; Have R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise, maximal value wherein is designated as Δ T SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes.Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " at the sea level elevation place that each is chosen; The absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
In the 3rd step, measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, reference plate maximum temperature difference Δ T PmaxWith Cable Structure surface maximum temperature difference Δ T SmaxAll be not more than 5 degrees centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front EmaxBe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T PmaxBe not more than Δ T after deducting 2 degrees centigrade Emax, and Cable Structure surface maximum temperature difference Δ T SmaxBe not more than Δ T PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the Cable Structure heat transfer characteristic; Utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T hThe numerical value that obtains is for Δ T is got in convenient narration hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure.
The second portion of this method: set up required knowledge base of cable structure health monitoring system and parameter method, based on knowledge base (containing parameter) and the actual measurement monitored amount Cable Structure health status appraisal procedure.Can carry out successively as follows, assessed by the health status of evaluation object more accurately with acquisition.
The first step: set up initial Mechanics Calculation benchmark model A o
Set up initial Mechanics Calculation benchmark model A oIn Cable Structure completion; Perhaps before setting up the cable structure health monitoring system; Calculate " Cable Structure steady temperature data " (can use conventional thermometry to measure, for example use thermal resistance to measure) according to " the temperature survey calculating method of the Cable Structure of this method " measurement, this moment " Cable Structure steady temperature data " are used vector T oExpression is called initial Cable Structure steady temperature data vector T oObtain T in actual measurement oThe time, just obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.Calculate initial Cable Structure steady temperature data vector T in actual measurement oThe time, use conventional method (consult reference materials or survey) to obtain the temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Calculate initial Cable Structure steady temperature data vector T in actual measurement oThe time, just obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, use the conventional method actual measurement to calculate the actual measurement computational data of Cable Structure.The actual measurement computational data of Cable Structure comprises that Non-Destructive Testing data of supporting rope etc. can express measured datas such as the data of the health status of rope, Cable Structure bearing initial angle displacement measurement data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, Cable Structure support coordinate data, Cable Structure modal data, structural strain data, structure angle measurement data, structure space measurement of coordinates data.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point on the structure, and purpose is to confirm according to these coordinate datas the geometric properties of Cable Structure.Cable Structure bearing initial angle displacement measurement data refer to setting up initial Mechanics Calculation benchmark model A oThe time, the Cable Structure bearing is with respect to the angular displacement that bearing took place under the Cable Structure design point.As far as cable-stayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on the bridge two ends, so-called bridge type data that Here it is.The Non-Destructive Testing data of utilization supporting rope etc. can be expressed the data and the Cable Structure bearing initial angle displacement measurement data of the health status of supporting rope and set up by evaluation object initial damage vector d oD is used in (shown in (1)) oThe expression Cable Structure is (with initial Mechanics Calculation benchmark model A oThe expression) by the initial health of evaluation object.If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, can think that perhaps the structure original state is a not damaged when not having relaxed state, vectorial d oIn get 0 with each element numerical value of supporting Suo Xiangguan, if there are not Cable Structure bearing initial angle displacement measurement data or can think that the displacement of Cable Structure bearing initial angle is at 0 o'clock, vectorial d oIn each element numerical value relevant with the angular displacement of Cable Structure bearing get 0.Utilize temperature variant physics and the mechanical property parameters and the initial Cable Structure steady temperature data vector T of the measured data of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, the Non-Destructive Testing data that support rope, Cable Structure bearing initial angle displacement measurement data, the employed various materials of Cable Structure o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial Mechanics Calculation benchmark model A o
d o=[d o1?d o2···d ok···d oN] T (1)
D in the formula (1) Ok(k=1,2,3 ...., N) the initial Mechanics Calculation benchmark model A of expression oIn k by the original state of evaluation object, if should be the rope (or pull bar) in the cable system, d so by evaluation object OkRepresent its initial damage, d OkBeing to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be angular displacement component, a d so of a bearing by evaluation object OkRepresent its initial displacement numerical value, T representes the transposition (back together) of vector.
Obtain T in actual measurement oThe time, just at the synchronization in the moment that obtains Cable Structure steady temperature data, use conventional method directly to measure the initial value of all monitored amounts of the Cable Structure that calculates, form monitored amount initial value vector C o(seeing formula (2)).Requirement is obtaining A oThe time obtain C o, monitored amount initial value vector C oExpression is corresponding to A oThe concrete numerical value of " monitored amount ".Because of subject to the foregoing, the monitored amount of calculating gained based on the calculating benchmark model of Cable Structure approaches the measured data of initial monitored amount reliably, in the narration of back, will represent this calculated value and measured value with prosign.
C o=[C o1?C o2···C oj···C oM] T (2)
C in the formula (2) Oj(j=1,2,3 ...., M) be the original bulk of j monitored amount in the Cable Structure, this component according to coding rule corresponding to specific j monitored amount.Vector C oBe to be formed according to certain series arrangement by the monitored amount of M, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
No matter which kind of method to obtain initial Mechanics Calculation benchmark model A with o, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T o), based on A oThe Cable Structure computational data that calculates must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this oSuo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of calculating gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when truly taking place near institute's analog case reliably.Model A oIn used by evaluation object initial damage vector d by the health status of evaluation object oExpression, Cable Structure steady temperature data are with initial Cable Structure steady temperature data vector T oExpression.Because based on A oThe evaluation that calculates all monitored amounts is very near the initial value (actual measurement obtains) of all monitored amounts, so also can be used in A oThe basis on, carry out Mechanics Calculation obtains, A oThe evaluation of each monitored amount form monitored amount initial value vector C oT oAnd d oBe A oParameter, also we can say C oBy A oMechanics Calculation result form.
Second step: circulation beginning.When beginning circulation each time, in the time of at first need setting up or set up this circulation beginning by the current initial damage vector of evaluation object d i o(i=1,2,3 ...), set up the current initial Mechanics Calculation benchmark model A of Cable Structure i o(finite element benchmark model for example, A in circulation each time i oBring in constant renewal in), A i oTemperature Distribution with " current initial Cable Structure steady temperature data vector T i o" express.Letter i is except the place of representing number of steps significantly, and alphabetical in the method i only representes cycle index, i.e. the i time circulation.A oAnd A i oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure.
What need during the i time circulation beginning is designated as d by the current initial damage vector of evaluation object i oD is used in (shown in (3)) i oCable Structure is (with current initial Mechanics Calculation benchmark model A when representing this time circulation beginning i oThe expression) by the health status of evaluation object.
d o i = d o 1 i d o 2 i . . . d ok i . . . d oN i T - - - ( 3 )
D in the formula (3) i Ok(i=1,2,3, K=1,2,3 ...., when N) the i time circulation of expression begins, current initial Mechanics Calculation benchmark model A i oIn k by the original state of evaluation object, if should be the rope (or pull bar) in the cable system, d so by evaluation object i OkRepresent its initial damage, d i OkBeing to represent not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be angular displacement component, a d so of a bearing by evaluation object i OkRepresent its initial displacement numerical value.
Set up and renewal d i oMethod following:
During circulation beginning for the first time, set up by the current initial damage vector of evaluation object and (be designated as d according to formula (3) 1 o) time, d 1 oJust equal d oI (i=2,3,4,5,6 ...) need when inferior circulation begins by the current initial damage vector of evaluation object d i o, be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) the preceding calculating acquisition of loop ends, concrete grammar is civilian the narration in the back.
I (i=1,2,3,4,5,6 ...) the Mechanics Calculation benchmark model of the Mechanics Calculation benchmark model that need set up when inferior circulation begins or the Cable Structure of having set up is designated as current initial Mechanics Calculation benchmark model A i oCorresponding to A i o" Cable Structure steady temperature data " use vector T i oExpression is called current initial Cable Structure steady temperature data vector T i oVector T i oDefinition mode and vector T oDefinition mode identical, when beginning circulation must be set up or set up and be called current initial Cable Structure steady temperature data vector T each time i o
Set up, upgrade A i oAnd T i oMethod following:
The Mechanics Calculation benchmark model of the Cable Structure of setting up during circulation beginning for the first time is designated as A 1 o, A 1 oEqual A o, T 1 oEqual T oA in circulation each time i oAnd T i oBring in constant renewal in, concrete grammar is the literary composition narration in the back; When loop ends each time, upgrade A i oAnd T i oThe Mechanics Calculation benchmark model of required Cable Structure when next time being circulated beginning, concrete grammar is the literary composition narration in the back.
This method is with " the current initial value vector of monitored amount C i o" (i=1,2,3 ...) initial value (referring to formula (4)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation of expression begins, C i oAlso can be called " the i time current initial value of the monitored amount of circulation vector ".
C o i = C o 1 i C o 2 i . . . C oj i . . . C oM i T - - - ( 4 )
C in the formula (2) i Oj(i=1,2,3, J=1,2,3 ...., j monitored amount when M) being the i time circulation beginning, in the Cable Structure.Vector C i oBe that M monitored amount by front definition forms according to certain series arrangement, this is put in order does not have specific (special) requirements, only require all associated vector of back also in this order array data get final product.
Setting up model A i oThe time set up " the current initial value vector of monitored amount C i o", the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts, C i oElement and C oElement corresponding one by one, represent that respectively all monitored amounts are in A in Cable Structure i oAnd A oConcrete numerical value during two states.
Set up and renewal C i oConcrete grammar following:
During circulation beginning for the first time, C 1 o(i=1, C i oBe embodied as C 1 o) equal C oI (i=2,3,4,5,6 ...) " the current initial value vector of the monitored amount C of the i time circulation of needs when inferior circulation begins i o", be preceding once (promptly the i-1 time, i=2,3,4,5,6 ...) calculate before the loop ends and obtain, concrete grammar is in the back literary composition narration.The i time (i=1,2,3,4,5,6 ...) in the circulation, " the current initial value vector of monitored amount C i o" bring in constant renewal in, concrete grammar is the literary composition narration in the back.Because according to model A i oThe initial value of calculating the monitored amount of gained approaches corresponding measured value reliably, in the narration of back, will represent this calculated value composition of vector and measured value composition of vector with prosign.
T i oAnd d i oBe A i oCharacterisitic parameter, C i oBe A i oAt T i oAnd d i oMechanics Calculation result under the condition forms.
The 3rd step: in Cable Structure military service process, in circulation each time, in other words in i (i=1,2,3,4,5,6 ...) in the inferior circulation, at known A i o, T i o, C i oAnd d i oAfter; Constantly survey the current data of calculating acquisition " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed " current cable structure steady temperature data vector T i", vector T iDefinition mode and vector T oDefinition mode identical; In the actual measurement vector T iThe time, just obtaining current cable structure steady temperature data vector T iThe synchronization in the moment, actual measurement obtains the currency of all monitored amounts in the Cable Structure, all these numerical value are formed the current numerical value vector of monitored amount C iC iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.
Obtaining vector T iAfter, upgrade A according to following concrete grammar i o, T i o, C i oAnd d i o:
Compare T iAnd T i oIf, T iEqual T i o, then need be to A i oUpgrade, otherwise need be to A i oAnd T i oUpgrade, update method is: the first step is calculated T iWith T oPoor, T iWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T iWith T oThe difference represent that with steady temperature change vector S S equals T iDeduct T o, S representes the variation of Cable Structure steady temperature data; Second step is to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading i o, upgrade A i oThe time, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, promptly upgraded T i o, so just obtained correctly corresponding to A i oT i oThis moment d i oRemain unchanged.When upgrading A i oAfter, A i oThe health status by evaluation object use by the current initial damage of evaluation object vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iExpression obtains A through Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, with these concrete numerical value replacements C i oMiddle corresponding element has so just been realized the current initial value vector of monitored amount C i oRenewal.
The 4th step: circulation time must be set up " unit damage monitored numerical quantity unit change matrix " and " unit damage or unit angular displacement vector " earlier each time, and " unit damage monitored numerical quantity unit change matrix " that the i time circulation set up is designated as Δ C i(i=1,2,3 ...)." unit damage or unit angular displacement vector " that i time circulation set up is designated as D i uΔ C in circulation each time iAnd D i uNeed according to circumstances to bring in constant renewal in, promptly upgrading current initial Mechanics Calculation benchmark model A i o, current initial Cable Structure steady temperature data vector T i oWith the current initial value vector of monitored amount C i oAfter, upgrade unit damage monitored numerical quantity unit change matrix Δ C iWith unit damage or the angular displacement vector D of unit i u
Earlier set up unit damage monitored numerical quantity unit change matrix Δ C during circulation beginning each time by following step iWith unit damage or the angular displacement vector D of unit i uIf in the 3rd step, upgraded A i o, in this step, must rebulid (promptly upgrading) unit damage monitored numerical quantity unit change matrix Δ C so iWith unit damage or the angular displacement vector D of unit i uIf in the 3rd step, do not upgrade A i o, in this step, needn't rebulid unit damage monitored numerical quantity unit change matrix Δ C so iWith unit damage or the angular displacement vector D of unit i uSet up and rebulid (promptly upgrading) Δ C iAnd D i uDetailed process identical, row as follows:
Current initial Mechanics Calculation benchmark model A in Cable Structure i oThe basis on carry out several times and calculate, equal the quantity of all ropes on the calculation times numerical value.Calculate each time hypothesis have only one by evaluation object (original damage or bearing angular displacement can be 0 in original damage or bearing angular displacement; Can not be 0 also) the basis on increase unit damage or unit angular displacement again; Concrete; If should be a supporting rope in the cable system by evaluation object; So just this supporting rope of hypothesis increases unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) again; If should be the angular displacement component of a direction of a bearing, just suppose that this bearing increases unit angular displacement (for example get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. and be the unit angular displacement) again in this angular displacement direction by evaluation object.Calculate for convenient; When set increasing unit damage or unit angular displacement in the circulation each time can all be structural health conditions during this time circulation beginning as being healthy fully, and set on this basis unit damage or unit angular displacement (in subsequent step, calculate, by the damage numerical value or the angular displacement of evaluation object---be called name damage d i c(i=1,2,3 ...), all with respect to this time when beginning circulation, will be by the health status of evaluation object as being healthy fully speech, therefore must according to after the formula that provides of the literary composition name that will calculate damage be converted into true damage).With a round-robin occur in calculating each time unit damage or unit angular displacement be different from by evaluation object other time occur in calculating unit damage or unit angular displacement by evaluation object; And what supposition had unit damage or a unit angular displacement each time can be different from other by the unit damage value of evaluation object or unit value of angular displacement by the unit damage value of evaluation object or unit value of angular displacement, with " unit damage or the angular displacement vector D of unit i u" (shown in (5)) write down that all are by the unit damage of the supposition of evaluation object or unit angular displacement in each time circulation, circulation time is designated as D for the first time 1 uCalculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, the current calculated value of the M of appointment monitored amount in front; The current calculated value that calculates gained M monitored amount is each time formed one " monitored amount calculation current vector ", and (when k of hypothesis had unit damage or unit angular displacement by evaluation object, available formula (6) was represented the monitored amount calculation current vector C of M monitored amount of all appointments i Tk); The monitored amount calculation current vector that calculates each time deducts the current initial value vector of monitored amount C i oObtain a vector; Each element that again should vector obtains one " the unit change vector of monitored amount " after all calculating the unit damage supposed or unit angular displacement numerical value divided by this; " the unit change vector of monitored amount " that be exactly under this condition (is mark with the position of rope that unit damage is arranged or numbering etc.) (when k is had unit damage or unit angular displacement by evaluation object, used δ C i kThe unit change vector of representing monitored amount, δ C i kDefinition see formula (7), formula (8) and formula (9), formula (7) deducts after the formula (4) again divided by vectorial D for formula (6) i uK element D i UkGained), the unit change of monitored amount vector δ C i kEach element representation since that supposition has unit damage or a unit angular displacement when calculating by the unit damage of evaluation object (for example k by evaluation object) or unit angular displacement (D for example i Uk), and the numerical value change amount of the pairing monitored amount of this element that causes is with respect to the unit damage or the angular displacement numerical value D of unit of supposition i UkRate of change; There is N N " the unit change vector of monitored amount " just to be arranged by evaluation object; The unit change vector of each monitored amount has M element, forms " the unit damage monitored numerical quantity unit change matrix Δ C that M * N element arranged successively by this N " the unit change vector of monitored amount " i" (the capable N row of M), each vectorial δ C i k(k=1,2,3 ...., N) be matrix Δ C iOne row, Δ C iDefinition suc as formula shown in (10).
D u i = D u 1 i D u 2 i . . . D uk i . . . D uN i T - - - ( 5 )
Unit damage or the angular displacement vector D of unit in the formula (5) i uElement D i Uk(i=1,2,3, K=1,2,3 ...., N) k of supposition is individual by the unit damage of evaluation object or unit angular displacement numerical value, vectorial D in the i time circulation of expression i uIn the numerical value of each element can be the same or different.
C tk i = C tk 1 i C tk 2 i . . . C tkj i . . . C tkM i T - - - ( 6 )
Elements C in the formula (6) i Tkj(i=1,2,3, K=1,2,3 ...., N; J=1,2,3 ...., M) the i time circulation of expression be because k is individual when by evaluation object unit damage or unit angular displacement being arranged, according to the current numerical value of calculating of the monitored amount of pairing j the appointment of coding rule.
δ C k i = C tk i - C o i D uk i - - - ( 7 )
The subscript i of each amount in the formula (7) (i=1,2,3 ...) the i time circulation of expression, subscript k (k=1,2,3 ...., N) k unit damage or unit angular displacement that is increased by evaluation object of expression, D in the formula i UkBe vectorial D i uIn k element.Vector δ C i kDefinition suc as formula shown in (7) and the formula (8), δ C i kJ (j=1,2,3 ...., M) individual element δ C i KjIn the i time circulation of (definition is suc as formula shown in (9)) expression, set up matrix Δ C iThe time, suppose the change amount of calculating a gained j monitored amount when k has unit damage or unit angular displacement by evaluation object unit damage or the angular displacement D of unit with respect to supposition i UkRate of change.
δ C k i = δ C k 1 i δ C k 2 i . . . δ C kj i . . . δ C kM i T - - - ( 8 )
δ C kj i = C tkj i - C oj i D uk i - - - ( 9 )
Δ C i = δ C 1 i δ C 2 i . . . δ C k i . . . δ C N i - - - ( 10 )
Vectorial δ C in the formula (10) i k(i=1,2,3 ....,, k=1,2,3 ...., N) in the i time circulation of expression, increased unit damage or the angular displacement D of unit by evaluation object because k is individual i UkCause, the relative value of all monitored amounts changes.Matrix Δ C iCoding rule and the front vector d of row (subscript k) i oThe coding rule of subscript k of element identical.
The 5th step: the current health status of identification Cable Structure.Detailed process is following.
I (i=1,2,3 ...) in the inferior circulation, be utilized in " the current numerical value vector of the monitored amount C that the actual measurement of second step obtains i" " the current initial value of monitored amount vector C together i o", " unit damage monitored numerical quantity unit change matrix Δ C i" and " the vectorial d of current name damage i c" between linear approximate relationship, shown in (11) or formula (12).
C i = C o i + Δ C i · d c i - - - ( 11 )
C i - C o i = Δ C i · d c i - - - ( 12 )
The current numerical value vector of monitored amount C in formula (11) and the formula (12) iDefinition be similar to the current initial value of monitored amount vector C i oDefinition, see formula (13); By the vectorial d of the current name damage of evaluation object i cDefinition see formula (14).
C i = C 1 i C 2 i . . . C j i . . . C M i T - - - ( 13 )
Elements C in the formula (13) i j(i=1,2,3 ....; J=1,2,3 ...., M) be the i time circulation time Cable Structure, according to the current numerical value of the monitored amount of the pairing j of being numbered of coding rule.
d c i = d c 1 i d c 2 i . . . d ck i . . . d cN i T - - - ( 14 )
D in the formula (14) i Ck(i=1,2,3 ....; K=1,2,3 ...., be that k the current name by evaluation object damaged or current nominal value of angular displacement vectorial d in the i time circulation N) i cThe coding rule and the matrix Δ C of subscript k of element iThe coding rule of row identical.
When supporting rope actual damage is not too big, the bearing angular displacement hour; Because the Cable Structure material still is in the linear elasticity stage; The distortion of Cable Structure is also less, and the represented a kind of like this linear relationship of formula (11) or formula (12) is less with the error of actual conditions, and error can be used error vector e i(formula (15)) definition, the error of linear relationship shown in expression (11) or the formula (12).
e i = abs ( Δ C i · d c i - C i + C o i ) - - - ( 15 )
Abs () is the function that takes absolute value in the formula (15), and each element of the vector of trying to achieve in the bracket is taken absolute value.
Because there are certain error in formula (11) or the represented linear relationship of formula (12), therefore can not be simply according to formula (11) or formula (12) and " the current numerical value vector of monitored amount C i" come directly to find the solution and obtain by the vectorial d of the current name damage of evaluation object i cIf done like this, what obtain is damaged vectorial d by the current name of evaluation object i cIn element in addition bigger negative value can appear, just negative damage, this obviously is irrational.Therefore obtain by the vectorial d of the current name damage of evaluation object i cAcceptable separating (promptly have reasonable error; But can be more exactly from cable system, confirm damaged cable position and degree of injury thereof, also can confirm bearing angular displacement numerical value more exactly) become a rational solution, available formula (16) is expressed this method.
abs ( Δ C i · d c i - C i + C o i ) ≤ g i - - - ( 16 )
Abs () is the function that takes absolute value in the formula (16), vectorial g iDescription departs from the reasonable deviation of ideal linearity relation (formula (11) or formula (12)), is defined by formula (17).
g i = g 1 i g 2 i . . . g j i . . . g M i T - - - ( 17 )
G in the formula (17) i j(i=1,2,3 ....; J=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (11) or the formula (12) in the i time circulation has been described.Vector g iCan be according to the error vector e of formula (15) definition iTentative calculation is selected.
At the current initial value vector of monitored amount C i o, unit damage monitored numerical quantity unit change matrix Δ C iWith the current numerical value vector of monitored amount C iWhen known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (16), obtain by the vectorial d of the current name damage of evaluation object i cAcceptable separating, by the current actual damage of evaluation object vector d iThe element of (formula (18) is seen in definition) can calculate according to formula (19), thereby can be by d iConfirm by the health status of evaluation object.
d i = d 1 i d 2 i . . . d k i . . . d N i T - - - ( 18 )
D in the formula (18) i k(i=1,2,3, K=1,2,3 ...., N) k is by the current actual health status of evaluation object in the i time circulation of expression, and formula (19) is seen in its definition, if this was the supporting rope (or pull bar) in the cable system, d so by evaluation object i kRepresent its current actual damage, d i kBeing to represent not damaged at 0 o'clock, is to represent that this supporting rope thoroughly lost load-bearing capacity at 100% o'clock, representes to lose the load-bearing capacity of corresponding proportion in the time of between 0 and 100%, if should be angular displacement component, a d so of a bearing by evaluation object i kRepresent its current actual displacement angle numerical value, vectorial d iCoding rule and the formula (1) of element in vectorial d oThe coding rule of element identical.
D in the formula (19) i Ok(i=1,2,3,4, K=1,2,3 ...., N) by the current initial damage vector of evaluation object d i oK element, d i CkBy the vectorial d of the current name damage of evaluation object i cK element.
Narration has obtained by the current actual damage vector of evaluation object d below iAfter, how to confirm the position and the relax level of slack line.
By before can know in the Cable Structure total M 1Root supporting rope, Cable Structure rope force data is by M 1The Suo Li of root supporting rope describes.Available " initial rope force vector F o" represent that all support the initial Suo Li (formula (20) is seen in definition) of ropes in the Cable Structure.
F o = F o 1 F o 2 . . . F oh . . . F oM 1 T - - - ( 20 )
F in the formula (20) o(h=1,2,3 ...., M 1) be the initial Suo Li of h root supporting rope in the Cable Structure, this element is according to the Suo Li of coding rule corresponding to appointment supporting rope.Vector F oIt is constant.Obtain T in actual measurement oSynchronization, use conventional method directly to measure the rope force data that calculates all supporting ropes, all these rope force datas are formed initial rope force vector F oSetting up the initial Mechanics Calculation benchmark model A of Cable Structure oThe time in fact used vectorial F o
Will be by the current actual damage vector of evaluation object d iIn with the supporting Suo Xiangguan M 1Individual element takes out, and forms supporting Suo Dangqian actual damage vector d Ic, supporting Suo Dangqian actual damage vector d IcThe coding rule and initial rope force vector F of element oThe coding rule of element identical.Supporting Suo Dangqian actual damage vector d IcH element representation Cable Structure in the current actual damage amount of h root supporting rope, h=1,2,3 ...., M 1Current actual damage vector d IcIn numerical value be not 0 element supporting rope corresponding to the unsoundness problem, to the supporting Suo Jinhang Non-Destructive Testing of these unsoundness problems, after Non-Destructive Testing found out that this supporting rope does not damage, this element numerical value (was used d so Ic hExpression) expression should supporting rope and d Ic hRelaxing of impairment value mechanics equivalence just confirmed slack line thus, and the computing method of concrete slack are explained below.
The current cable power (formula (21) is seen in definition) of all supporting ropes in the Cable Structure that obtains with " current cable force vector F " expression actual measurement in this method.
F = F 1 F 2 . . . F h . . . F M 1 T - - - ( 21 )
F in the formula (21) h(h=1,2,3 ...., M 1) be the current cable power of h root supporting rope in the Cable Structure.Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains in the Cable Structure rope force data of all supporting ropes, all these rope force datas are formed current cable force vector F.Element and the vectorial F of vector F oThe coding rule of element identical.According to the narration of front, vector T i oEqual vector T i
In this method, under supporting rope original state, in the steady temperature data of Cable Structure with initial Cable Structure steady temperature data vector T oDuring expression, and the supporting rope is when being in free state (free state refers to that Suo Li is 0, back with), and the length of supporting rope is called initial drift, with " initial drift vector l o" represent that all support the initial drift (formula (22) is seen in definition) of ropes in the Cable Structure.Pass through vector T according to " the temperature survey calculating method of the Cable Structure of this method " oCan confirm obtaining vector T oThe Temperature Distribution of all supporting ropes constantly.
l o = l o 1 l o 2 . . . l oh . . . l oM 1 T - - - ( 22 )
L in the formula (22) Oh(h=1,2,3 ...., M 1) be the initial drift of h root supporting rope in the Cable Structure.Vector l oBe constant, after when beginning, confirming, just no longer change.
Similarly, under supporting rope original state, in the steady temperature data of Cable Structure with initial Cable Structure steady temperature data vector T oDuring expression, and the supporting rope is when being in free state, and the cross-sectional area of supporting rope is called initial free cross-sectional area, with " initial free cross-sectional area vector A o" represent that all support the initial free cross-sectional area (formula (23) is seen in definition) of ropes in the Cable Structure, the weight of the unit length of supporting rope is called the weight of initial free unit length, with " the weight vector ω of initial free unit length o" represent that all support the weight (formula (24) is seen in definition) of the initial free unit length of ropes in the Cable Structure.
A o = A o 1 A o 2 . . . A oh . . . A oM 1 T - - - ( 23 )
A in the formula (23) Oh(h=1,2,3 ...., M 1) be the initial free cross-sectional area of h root supporting rope in the Cable Structure.Vector A oBe constant, after when beginning, confirming, just no longer change.
ω o = ω o 1 ω o 2 . . . ω oh . . . ω o M 1 T - - - ( 24 )
ω in the formula (24) Oh(h=1,2,3 ...., M 1) be the weight of the free unit length of initial freedom of h root supporting rope in the Cable Structure.Vector ω oBe constant, after when beginning, confirming, just no longer change.
In this method, in the steady temperature data of Cable Structure with current initial Cable Structure steady temperature data vector T t oDuring expression, with " current initial drift vector l t o" the current initial drift of all supporting ropes in the expression Cable Structure (and formula (25) is seen in definition, refers to that hypothesis supporting cable force is at 0 o'clock, considered the influence of thermal expansivity and temperature variation to supporting rope drift after, the vectorial l of initial drift oWith initial Cable Structure steady temperature data vector T oThe expression the supporting rope in temperature with current initial Cable Structure steady temperature data vector T i oSupporting rope drift during expression).According to " the temperature survey calculating method of the Cable Structure of this method ", through vector T i oCan confirm obtaining vector T i oThe Temperature Distribution of all supporting ropes constantly.
l o t = l o 1 t l o 2 t . . . l oh t . . . l oM 1 t T - - - ( 25 )
L in the formula (25) t Oh(h=1,2,3 ...., M 1) be that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T i oDuring expression, the current initial drift of h root supporting rope can be utilized the thermal expansivity, the l that support rope in the Cable Structure Oh, T oAnd T i oObtain l through conventional physical computing t Oh
Element, the vectorial l of vector F oElement, vectorial l t oElement, vectorial A oElement, vectorial ω oElement and vectorial F oThe coding rule of element identical, the different information of the same supporting rope of element representation of the identical numbering of these vectors.
In this method, in the steady temperature data of Cable Structure with current initial Cable Structure steady temperature data vector T i oDuring expression, with the current drift (formula (26) is seen in definition, and support Suo Keneng this moment is intact, possibly be impaired, also possibly relax) of all supporting ropes in " current drift vector l " expression Cable Structure.
l = l 1 l 2 . . . l h . . . l M 1 T - - - ( 26 )
L in the formula (26) h(h=1,2,3 ...., M 1) be the current drift of h root supporting rope in the Cable Structure.
In this method, represent the change amount (formula (27) and formula (28) are seen in definition) of the drift of all supporting ropes in the Cable Structure with " drift changes vectorial Δ l " (or claiming supporting Suo Dangqian relax level vector).
Δl=[Δl 1?Δl 2···Δl j···Δl N] T (27)
Δ l in the formula (27) h(h=1,2,3 ...., M 1) be the change amount of the drift of h root supporting rope in the current cable structure, formula (28), Δ l are seen in its definition hBe not that 0 rope is a slack line, Δ l hNumerical value be the slack of rope, and the current relax level of expression cable system h root supporting rope, the long adjustment amount of rope of this rope when also being adjustment Suo Li.
Δ l h = l h - l oh t - - - ( 28 )
Through slack line is carried out the relax level identification that slack line is carried out in the mechanics equivalence with damaged cable, the mechanical condition of equivalence is in the method:
The mechanics parameters of initial drift, geometrical property parameter and material when one, the nothing of the rope of two equivalences relaxes with not damaged is identical;
Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.
When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after promptly replacing slack line with the damaged cable of equivalence, Cable Structure any variation can not take place, vice versa.
Obtained supporting Suo Dangqian actual damage vector d IcAfter, d IcH element d Ic h(h=1,2,3 ...., M 1) represent that the h root supports the actual damage value of rope, formula (19) is seen in its definition, though with d Ic hBe called the actual damage value of h root supporting rope or the actual damage degree of h root supporting rope, but since h root Suo Keneng be impaired also possibly be lax, so d IcH element d Ic hThe actual damage value of the h root supporting rope of expression is actually the actual equivalent damage value of h root supporting rope, when h root supporting rope is actually impaired, and d Ic hJust the h root of expression supports the actual damage value of rope, when h root supporting rope is actually lax, and d Ic hH root supporting rope and the lax equivalent actual damage value just represented for narrating conveniently, are claimed d in the method Ic hBeing to represent h root supporting rope not damaged at 0 o'clock, is to represent that this rope thoroughly lost load-bearing capacity at 100% o'clock, representes the load-bearing capacity of h root supporting rope forfeiture corresponding proportion in the time of between 0 and 100%, through by the vectorial d of the current actual damage of evaluation object IcJust can identify the supporting rope that health status goes wrong, but some is impaired in the supporting rope that these health status go wrong, some is to have relaxed, if h root supporting rope is actually and has relaxed that (its current relax level is used Δ l jDefinition), the current relax level Δ l of so lax h root supporting rope j(Δ l jDefinition see formula (27)) with the current actual damage degree d of damaged cable of equivalence Ic hBetween relation confirm by aforementioned two mechanics equivalent conditions.Δ l hSame d Ic hBetween physical relationship can adopt accomplished in many ways; For example can directly confirm (referring to formula (29)) according to aforementioned equivalent condition; Also can adopt based on the Ernst equivalent elastic modulus to replace the E in the formula (25) to revise back definite (referring to formula (30)), also can adopt and confirm based on other methods such as trial and error procedure of finite element method.
Δl h t = d h ci 1 - d h ci F h E h t A h t + F h l oh t - - - ( 29 )
Δl h t = d h ci 1 - d h ci F h [ E h t 1 + ( ω h t l xh t ) 2 A h t E h t 12 ( F h ) 3 ] A h t + F h l oh t - - - ( 30 )
E in formula (29) and the formula (30) t hBe that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T i oDuring expression, the elastic modulus of h root supporting rope, A t hBe that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T i oDuring expression, the cross-sectional area of h root supporting rope, F hBe that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T i oDuring expression, the current cable power of h root supporting rope, d i hBe the current actual damage degree of h root supporting rope, ω t hBe that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T t oDuring expression, the weight of the unit length of h root supporting rope, l t XhBe that steady temperature data in Cable Structure are with current initial Cable Structure steady temperature data vector T i oDuring expression, the horizontal range of two supporting end points of h root supporting rope, l t XhBe current supporting rope two supporting end points horizontal range vector l t xAn element, current supporting rope two supporting end points horizontal ranges vector l t xThe coding rule and initial drift vector l of element oThe coding rule of element identical, E t hCan obtain A according to the characteristic material data of looking into or survey h root supporting rope t hAnd ω t hCan be according to thermal expansivity, the A of h root supporting rope Oh, ω Oh, F h, T oAnd T i oObtain through conventional physics and Mechanics Calculation.Item in the formula (30) in [] is the Ernst equivalent elastic modulus of this supporting rope, can just can confirm to support Suo Dangqian relax level vector Δ l by formula (29) or formula (30).Formula (30) is the correction to formula (29).
The 6th step: judge whether to finish this (the i time) circulation, if, then accomplish the preceding tailing in work of this loop ends, for next time (promptly the i+1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is following:
In this (the i time) circulation, try to achieve the vectorial d of current name damage i cAfter, at first, set up mark vector B according to formula (31) i, formula (32) has provided mark vector B iThe definition of k element; If mark vector B iElement be 0 entirely, then got back to for the 3rd step and proceed health monitoring and calculating Cable Structure; If mark vector B iElement be not 0 entirely, then accomplish subsequent step after, get into circulation next time.
So-called subsequent step is: at first, according to formula (33) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation I+1 oEach element d I+1 OkThe second, at Mechanics Calculation benchmark model A oThe basis on, make A oIn be d by the health status of evaluation object I+1 oRather than be d oAfter, more further to A oIn Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation that applies just taken from steady temperature change vector S, and steady temperature change vector S equals T iDeduct T o), so just obtained next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Mechanics Calculation benchmark mould A of circulation I+1 o, next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual T i o, to A I+1 oCarrying out Mechanics Calculation obtains corresponding to A I+1 oConcrete numerical value all monitored amounts, current, these concrete numerical value form next time (promptly the i+1 time, i=1,2,3,4 ...) the current initial value vector C of the required monitored amount of circulation I+1 o
B i = B 1 i B 2 i . . . B k i . . . B N i T - - - ( 31 )
Mark vector B in the formula (31) iSubscript i represent the i time the circulation, its element B i k(k=1,2,3 ..., subscript k N) representes that k can only be got 0 and 1 two amount by the health status characteristic of evaluation object, concrete value rule is seen formula (32).
Element B in the formula (32) i kBe mark vector B iK element, D i UkBe unit damage or the angular displacement vector D of unit i uK element (seeing formula (3)), d i CkBy the vectorial d of the current name damage of evaluation object i cK element (seeing formula (14)), they represent that all k is by the relevant information of evaluation object.
D in the formula (33) i UkBe unit damage or the angular displacement vector D of unit i uK element (seeing formula (5)), d i OkBy the current initial damage vector of evaluation object d i oK element (seeing formula (3)).
The third part of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (comprising monitored amount monitoring system, temperature monitoring system), signal picker and computing machine etc.Require the temperature required measured data of monitoring acquisition in real time, require simultaneously each monitored amount of monitoring in real time.
Software section should be accomplished the process that this method sets, and promptly accomplishes functions such as needed in this method, as can to use computer realization monitoring, record, control, storage, calculating, notice, warning.
This method specifically comprises:
A. for for the purpose of narration is convenient, this method is unified claims that the supporting rope of being assessed is that establishing the quantity of the supporting rope of being assessed and the quantity sum of bearing angular displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with the bearing angular displacement component; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; This method is represented this numbering with variable k, k=1, and 2,3 ..., N; If total M in the cable system 1Root supporting rope, Cable Structure rope force data comprises this M 1The Suo Li of root supporting rope, obviously M 1Less than by the quantity N of evaluation object; Confirm the measured point of appointment, give all specified point numberings; Confirm the measured straight line of each measurement point, gave the measured straight line numbering of all appointments; Confirm the measured angle coordinate component of each measured straight line, give all measured angle coordinate component numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored angle-data of Cable Structure " is made up of above-mentioned all measured angle coordinate components; For simplicity, in the method " the monitored angle-data of Cable Structure " abbreviated as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " set by step b1 to b3 carry out;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure; Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; One be called the cloudy day all day with what can not see the sun daytime in the method; Statistics obtain in T cloudy day each cloudy day 0 up to the sunrise next day highest temperature and the lowest temperature between back 30 minutes constantly, sunrise is meant the sunrise moment on the meteorology that base area revolutions and revolution rule confirm constantly, does not represent necessarily can see the sun same day; The sunrise that can inquire about data or calculate each required day through conventional meteorology constantly; Each cloudy day 0 up to next day sunrise constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, T cloudy day arranged, the maximum temperature difference of T cloudy daily temperature is just arranged; Get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data "; From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness "; Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is called reference plate; Reference plate can not contact with ground; Reference plate overhead distance is not less than 1.5 meters, and the one side of this reference plate is called sunny slope on the sunny side; The sunny slope of reference plate is coarse and dark color; The sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day, the non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate;
B2: monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains Cable Structure place environment; Be called the environment maximum temperature difference, be designated as Δ T EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day of temperature that minimum temperature obtains the sunny slope of reference plate; Be called the reference plate maximum temperature difference, be designated as Δ T PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; There is R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise; Maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes; Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " at the sea level elevation place that each is chosen; The absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
B3: measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, reference plate maximum temperature difference Δ T PmaxWith Cable Structure surface maximum temperature difference Δ T SmaxAll be not more than 5 degrees centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front EmaxBe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T PmaxBe not more than Δ T after deducting 2 degrees centigrade Emax, and Cable Structure surface maximum temperature difference Δ T SmaxBe not more than Δ T PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the Cable Structure heat transfer characteristic; Utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T hThe numerical value that obtains is for Δ T is got in convenient narration hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure;
C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " and calculate the Cable Structure steady temperature data under the original state; Cable Structure steady temperature data under the original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T o"; Survey or consult reference materials and obtain the temperature variant physics and the mechanical property parameters of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all supporting ropes, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all supporting ropes are in free state that Suo Li is 0 o'clock a length, the weight of cross-sectional area during in free state and the unit length during in free state; And the temperature of all supporting ropes when obtaining these three kinds of data; Utilize the temperature variant physical function parameter and the mechanical property parameters of all supporting ropes on this basis, obtain all supporting ropes at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that length, the Suo Li of 0 o'clock all supporting rope is that 0 o'clock all cross-sectional area and Suo Li that supports ropes is the weight of the unit length of 0 o'clock all supporting rope; Form the initial drift vector of supporting rope, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule and initial rope force vector F of the element of the initial drift vector of supporting rope, the weight vector of the initial free unit length of initial free cross-sectional area vector sum oThe coding rule of element identical; Obtain T in actual measurement oThe time, just obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment; Directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure comprises the Non-Destructive Testing data, Cable Structure bearing initial angle displacement measurement data, the initial value of all monitored amounts, the initial rope force data of all supporting ropes, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure support coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data of the health status of expressing the supporting rope; The initial value of all monitored amounts is formed monitored amount initial value vector C o, monitored amount initial value vector C oThe coding rule of coding rule and M monitored amount identical; Utilize the Non-Destructive Testing data and the Cable Structure bearing initial angle displacement measurement data of the health status that can express the supporting rope to set up by evaluation object initial damage vector d o, vectorial d oExpression is with initial Mechanics Calculation benchmark model A oThe expression Cable Structure by the initial health of evaluation object; By evaluation object initial damage vector d oElement number equal N, d oElement be one-to-one relationship by evaluation object, vectorial d oThe coding rule of element with identical by the coding rule of evaluation object; If d oSome elements corresponding be supporting rope, a d so in the cable system by evaluation object oThe numerical value of this element represent the initial damage degree of corresponding supporting rope; If the numerical value of this element is 0, represent that the pairing supporting rope of this element is intact, do not damage; If its numerical value is 100%; Represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this supporting rope; If d oSome elements corresponding be some angular displacement components of some bearings, d so by evaluation object oThe numerical value of this element represent the initial value of this angular displacement component of this bearing; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, can think that perhaps the Cable Structure original state is a not damaged when not having relaxed state, vectorial d oIn get 0 with each element numerical value of supporting Suo Xiangguan, if there are not Cable Structure bearing initial angle displacement measurement data or can think that the displacement of Cable Structure bearing initial angle is at 0 o'clock, vectorial d oIn each element numerical value relevant with the angular displacement of Cable Structure bearing get 0; Initial Cable Structure support coordinate data refer to the support coordinate data under the Cable Structure design point, and Cable Structure bearing initial angle displacement measurement data refer to setting up initial Mechanics Calculation benchmark model A oThe time, the Cable Structure bearing is with respect to the angular displacement that bearing took place under the Cable Structure design point;
D. according to the temperature variant physics of the Non-Destructive Testing data of the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, supporting rope, Cable Structure bearing initial angle displacement measurement data, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure steady temperature data vector T oWith all Cable Structure data that obtain with preceding step, set up the initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " o, based on A oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A oUsed by the evaluation object health status by evaluation object initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression; T oAnd d oBe A oParameter, by A oThe initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result oThe initial value of all monitored amounts of expression is identical, therefore also we can say C oBy A oMechanics Calculation result form A in the method o, C o, d oAnd T oBe constant;
E. in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning i o, A oAnd A i oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; What need during the i time circulation beginning is designated as d by the current initial damage vector of evaluation object i o, d i oCable Structure A when representing this time circulation beginning i oBy the health status of evaluation object, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i oAnd d i oBe A i oCharacterisitic parameter, C i oBy A i oMechanics Calculation result form; During circulation beginning for the first time, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring circulation beginning for the first time, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring circulation beginning for the first time, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring circulation beginning for the first time, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
F. go on foot the s circulation in step from getting into here by f; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T i, vector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement corresponding one by one; Obtain vector T in actual measurement iThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T iThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value are formed the current numerical value vector of monitored amount C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains all M in the Cable Structure 1The rope force data of root supporting rope, all these rope force datas are formed current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement calculates all M 1The volume coordinate of two supporting end points of root supporting rope; The volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges; Two supporting end points horizontal range data of all supporting ropes are formed current supporting rope two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current supporting rope two supporting end points horizontal range vectors oThe coding rule of element identical;
G. according to current cable structure steady temperature data vector T i, upgrade current initial Mechanics Calculation benchmark model A according to step g 1 to g3 i o, the current initial value of monitored amount vector C i oWith current initial Cable Structure steady temperature data vector T i o, and by the current initial damage vector of evaluation object d i oRemain unchanged;
G1. compare T iWith T i oIf, T iEqual T i o, A then i o, C i oAnd T i oRemain unchanged; Otherwise need follow these steps to A i oAnd T i oUpgrade;
G2. calculate T iWith T oPoor, T iWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T iWith T oThe difference represent that with steady temperature change vector S S equals T iDeduct T o, S representes the variation of Cable Structure steady temperature data;
G3. to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading i o, upgrade A i oThe time, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, promptly upgraded T i o, so just obtained correctly corresponding to A i oT i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status of rope use by the current initial damage vector of evaluation object d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A through Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C i o
H. at current initial Mechanics Calculation benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity unit change matrix Δ C through calculating iWith unit damage or the angular displacement vector D of unit i u
H1. the i time when beginning circulation, directly h2 obtains Δ C to the listed method of step h4 set by step iAnd D i uAt other constantly, when in step g to A i oAfter upgrading, h2 regains Δ C to the listed method of step h4 set by step iAnd D i uIf, in step g not to A i oUpgrade, then directly change step I herein over to and carry out follow-up work;
H2. at the current initial Mechanics Calculation benchmark model A of Cable Structure i oThe basis on carry out the several times Mechanics Calculation, equal all on the calculation times numerical value by the quantity N of evaluation object, have N evaluation object that N calculating is just arranged; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or angular displacement, to be increased unit damage or unit angular displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object, so just this supporting rope of hypothesis increases unit damage again, if should be the angular displacement component of a direction of a bearing by evaluation object; Just suppose that this bearing increases the unit angular displacement again at this sense of displacement, uses D i UkWrite down the unit damage or the unit angular displacement of this increase, wherein k represent to increase unit damage or unit angular displacement by the numbering of evaluation object, D i UkBe unit damage or the angular displacement vector D of unit i uAn element, unit damage or unit angular displacement vector D i uThe coding rule and the vectorial d of element oThe coding rule of element identical; Increase again in calculating each time unit damage or unit angular displacement be different from by evaluation object other time increase again in calculating unit damage or unit angular displacement by evaluation object; Calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount calculation current vector; When hypothesis k increases unit damage or unit angular displacement again by evaluation object, use C i Tk" the monitored amount calculation current vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C i TkDefinition mode and vectorial C oDefinition mode identical, C i TkElement and C oElement corresponding one by one;
H3. the vectorial C that calculates each time i TkDeduct vectorial C i oObtain a vector, each element that again should vector obtains " unit change vector δ a C of monitored amount after all calculating the unit damage supposed or unit angular displacement numerical value divided by this i k"; There is N N " the unit change vector of monitored amount " just to be arranged by evaluation object;
H4. individual by the coding rule of evaluation object by this N " the unit change vector of monitored amount " according to N, form " the unit damage monitored numerical quantity unit change matrix Δ C that the N row are arranged successively i"; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to the different unit change amplitude of same monitored amount when difference is increased unit damage or unit angular displacement by evaluation object; Unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule and the vectorial d of row oThe coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule of coding rule and M monitored amount of row identical;
I. define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal to support the quantity of rope, d i cAnd d iElement and be one-to-one relationship between the evaluation object, d i cAnd d iElement numerical value represent corresponding by the degree of injury of evaluation object or health status, d i cAnd d iWith by evaluation object initial damage vector d oThe element coding rule identical, d i cElement, d iElement and d oElement be one-to-one relationship;
J. according to the current numerical value vector of monitored amount C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity unit change matrix Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + Δ C i · d c i Formula 1
K. the current actual damage vector d that utilizes formula 2 to express iK element d i kTogether by the current initial damage vector of evaluation object d i oK element d i OkWith the vectorial d of current name damage i cK element d i CkBetween relation, calculate current actual damage vector d iAll elements;
formula 2
K=1 in the formula 2,2,3 ..., N; d i kRepresent that k is individual by the current actual health status of evaluation object, d in the i time circulation i kBe to represent that k was not had health problem by evaluation object, d at 0 o'clock i kNumerical value be not represented in 0 o'clock k by evaluation object be the unsoundness problem by evaluation object, if should be supporting rope, a d so in the cable system by evaluation object i kThe order of severity of representing its current health problem, the supporting Suo Keneng of unsoundness problem is slack line, also possibly is damaged cable, d i kNumerical response the degree of lax or damage that should the supporting rope, if should be angular displacement component, a d so of a bearing by evaluation object i kRepresent its current actual displacement angle numerical value; So according to by the current actual damage of evaluation object vector d iCan define which the supporting rope unsoundness problem and order of severity thereof, define which bearing angular displacement and numerical value thereof have taken place;
L. will be by the current actual damage vector of evaluation object d iIn with M 1The M of root supporting Suo Xiangguan 1Individual element takes out, and forms supporting Suo Dangqian actual damage vector d Ca, supporting Suo Dangqian actual damage vector d CaThe coding rule and initial rope force vector F of element oThe coding rule of element identical; Supporting Suo Dangqian actual damage vector d CaH element representation Cable Structure in the current actual damage amount of h root supporting rope, h=1,2,3 ...., M 1Supporting Suo Dangqian actual damage vector d CaIn numerical value be not 0 element supporting rope corresponding to the unsoundness problem, from the supporting rope of these unsoundness problems, identify damaged cable, remaining is exactly slack line, supporting Suo Dangqian actual damage vector d CaIn with slack line corresponding to the element numerical expression be the current actual equivalent damage degree with slack line relax level mechanics equivalence;
M. be utilized in current cable structure steady temperature data vector T iSlack line under the condition, that identify in l step and with the vectorial d of supporting Suo Dangqian actual damage CaThese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, be utilized in that the f step obtains at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current supporting rope two supporting end points horizontal ranges vectors, be utilized in that the c step obtains at initial Cable Structure steady temperature data vector T oThe initial drift vector of the supporting rope under the condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, initial rope force vector F o, utilize current cable structure steady temperature data vector T tThe supporting Suo Dangqian steady temperature data of expression, be utilized in that c step obtains at initial Cable Structure steady temperature data vector T oThe supporting rope initial steady state temperature data of expression; Be utilized in the temperature variant physics and the mechanical property parameters of the employed various materials of Cable Structure of c step acquisition; Count the influence of temperature variation to supporting rope physics, mechanics and geometric parameter; Through with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two supporting ropes in Cable Structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of supporting rope drift, has just confirmed the long adjustment amount of rope of the supporting rope that those need adjust Suo Li; The lax identification and the damage identification of supporting rope have so just been realized; Institute's demand power is provided by current cable force vector F corresponding element during calculating;
N. try to achieve the vectorial d of current name damage i cAfter, set up mark vector B according to formula 3 i, formula 4 has provided mark vector B iThe definition of k element;
B i = B 1 i B 2 i . . . B k i . . . B N i T Formula 3
formula 4
Element B in the formula 4 i kBe mark vector B iK element, D i UkBe unit damage or the angular displacement vector D of unit i uK element, d i CkBy the vectorial d of the current name damage of evaluation object i cK element, they all represent k by the relevant information of evaluation object, k=1 in the formula 4,2,3 ..., N;
If mark vector B o. iElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B iElement be not 0 entirely, then get into next step, be step p;
P. according to formula 5 calculate next time, i.e. the i+1 time circulation be required by the current initial damage vector of evaluation object d I+1 oEach element;
formula 5
D in the formula 5 I+1 OkBe next time, i.e. the i+1 time circulation be required by the current initial damage vector of evaluation object d I+1 oK element, d i OkBe this, promptly the i time round-robin be by the current initial damage vector of evaluation object d i oK element, D i UkBe the i time round-robin unit damage or the angular displacement vector D of unit i uK element, B i kBe the i time round-robin mark vector B iK element, k=1 in the formula 5,2,3 ..., N;
Q. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual the current initial Cable Structure steady temperature data vector T of round-robin the i time i o
R. at initial Mechanics Calculation benchmark model A oThe basis on, to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, the health status that makes rope again is d I+1 oAfter obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation I+1Obtain A I+1After, obtain A through Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value are formed next time, the vectorial C of the current initial value of required monitored amount that promptly circulates for the i+1 time I+1 o
S. get back to step f, beginning is circulation next time.
Beneficial effect: when receiving the influencing of factors such as sunshine and environment temperature when the temperature field of Cable Structure; The temperature field of Cable Structure is constantly to change; The change of temperature field of Cable Structure must influence the monitored amount of Cable Structure; Have only and partly rejected by the influence in temperature field monitored amount could to carry out rational cable structure health monitoring based on monitored amount; And the temperature field measurement of Cable Structure is very complicated with calculating; This method discloses and has comprised a kind of simple, economic, feasible, cable structure health monitoring method of Cable Structure temperature field computing method efficiently that is suitable for cable structure health monitoring, adopts this method to occur at the Cable Structure bearing under the situation of angular displacement, when many ropes of Cable Structure are synchronously impaired or lax; And when the temperature of Cable Structure changed along with the time; Disclose and a kind ofly approached non-linear method with linear segmented, monitoring and evaluation identifies damaged cable, slack line and bearing angular displacement very exactly, and the disclosed system and method for this method is very useful to effective health monitoring of Cable Structure.
Embodiment
When temperature variation, to the damaged cable of Cable Structure and bearing angular displacement identification, this method discloses a kind ofly can monitor in the identification Cable Structure each by the system and method for the health status of evaluation object rationally and effectively.The following explanation of the embodiment of this method in fact only is exemplary, and purpose never is to limit the application or the use of this method.
This method adopts a kind of algorithm, and this algorithm is used to discern damaged cable, slack line and bearing angular displacement.During practical implementation, the following step is a kind of in the various steps that can take.
The first step: the quantity sum of bearing angular displacement component of quantity and Cable Structure of establishing the supporting rope of Cable Structure is N.For the purpose of narrating conveniently, this method is unified claims that support cable and the bearing angular displacement assessed are " evaluation object ", total N evaluation object.Give by the evaluation object serial number, this numbering will be used to generate the vector sum matrix in subsequent step, and this method is represented this numbering with variable k, k=1, and 2,3 ..., N.
If total M in the cable system 1Root supporting rope, Cable Structure rope force data comprises this M 1The Suo Li of root supporting rope, obviously M 1Less than by the quantity N of evaluation object.
Confirm that (promptly all characterize the specified point of Cable Structure angle displacements for the measured point of appointment; Be provided with K specified point); Each specified point can be exactly a near point the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) of each root rope; This specified point can also be a near point the Cable Structure bearing, gives all specified point numberings; Confirm the measured straight line (established each measurement point L appointment straight line arranged) of each measurement point, gave the measured straight line numbering of all appointments; Confirm the measured angle coordinate component (establish each measured straight line H angle coordinate component arranged) of each measured straight line, give all measured angle coordinate component numberings.Above-mentioned numbering will be used to generate the vector sum matrix equally in subsequent step." the whole monitored angle-data of Cable Structure " described by H angle coordinate component L appointment straight line K specified point, that cross each specified point, each appointment straight line on top definite Cable Structure, and the variation of Cable Structure angle is exactly the variation of the angle coordinate component of all appointments all specified points, all appointment straight lines.(individual angle coordinate component measurement value of M=K * L * H) or calculated value characterize the angle information of Cable Structure to each total M.K and M must not be less than by the quantity N of evaluation object.
For simplicity, in the method " all monitored parameters of Cable Structure " are abbreviated as " monitored amount ".Give M monitored amount serial number, this numbering will be used to generate the vector sum matrix in subsequent step.This method is with representing this numbering with variable j, j=1, and 2,3 ..., M.
Confirm " the temperature survey calculating method of the Cable Structure of this method " that these method concrete steps are following:
The a step: inquiry or actual measurement (can use conventional thermometry to measure; For example use thermal resistance to measure) obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure (for example finite element model).Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; Statistics obtains in T cloudy day 0 the highest temperature and the lowest temperature between back 30 minutes of the moment of sunrise next day at each cloudy day; Sunrise be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunrise constantly, the sunrise that can inquire about data or calculate each required day through conventional meteorology constantly, each cloudy day 0 up to sunrise next day constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature; T cloudy day arranged; The maximum temperature difference of daily temperature that T cloudy day just arranged, the maximal value of getting in the maximum temperature difference of daily temperature at T cloudy day is with reference to temperature difference per day, is designated as Δ T with reference to temperature difference per day rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual observation record; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data ".From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level,, can choose height above sea level 0m, 50m, 100m and height above sea level 200m so if for example the sea level elevation of Cable Structure is between 0m to 200m; Intersect with imaginary surface level and Cable Structure surface at each sea level elevation place that chooses; Obtain intersection, surface level and Cable Structure are crossing to obtain friendship face, and intersection is the outer edge line of friendship face; Intersection place on surface level and Cable Structure surface chooses 6 points; From the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", and it is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure.In the measurement Cable Structure of choosing along in 6 directions of the Temperature Distribution of wall thickness; At first according to the sunny slope of definite Cable Structure such as the physical dimension of the meteorological data throughout the year in Cable Structure position zone and Cable Structure, volume coordinate, Cable Structure surrounding environment and in the shade, the sunny slope of Cable Structure and in the shade face are the parts on the surface of Cable Structure, at each sea level elevation place that chooses; Aforementioned intersection respectively has one section in sunny slope and in the shade; This of intersection respectively has a mid point for two sections, crosses these two mid points and gets the outer normal of Cable Structure, and this method is called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure with these two outer normals; This method is called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure with these two outer normal directions; The outer normal of obvious sunny slope and in the shade outer normal all intersect with aforementioned intersection, and two intersection points are also just arranged, and these two intersection points are divided into two line segments with intersection; On two line segments, get 2 points respectively; Totally 4 points, each line segment is divided into 3 sections of equal in length in two line segments of the intersection of naming a person for a particular job of getting, and gets the outer normal on Cable Structure surface at these 4 some places; Just chosen the outer normal on 6 Cable Structure surfaces so altogether at each selected sea level elevation place, the direction of 6 outer normals is exactly " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures the direction of Cable Structure along the Temperature Distribution of wall thickness " line and Cable Structure, if Cable Structure is hollow, one at these two intersection points are on the Cable Structure outside surface; Another is on inside surface, if Cable Structure is solid, these two intersection points are all on the Cable Structure outside surface; Connect these two intersection points and obtain a straight-line segment, on straight-line segment, choose three points again, these three these straight-line segments of naming a person for a particular job are divided into four sections; Three points measuring that Cable Structure chooses at this and two end points of straight-line segment, the temperature of 5 points altogether; Concrete can hole on Cable Structure earlier, and temperature sensor is embedded in this 5 some places, special; Can not on the supporting rope, hole; The supporting rope is only measured the temperature that supports the rope surface point, and in any case, the temperature that records all is called this place " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " " Cable Structure is along temperature profile data of thickness " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measurement acquisition are called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method.If chosen H different altitude above sea level; At each sea level elevation place, chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness, measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness; Wherein H and E are not less than 3; B is not less than 2, and is special, and E equals 1 for the supporting rope; That " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE; The back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement, claim to survey the temperature data that obtains and will be " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure; Obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual observation record; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground getable this day fullest sunshine (as long as the same day sunrise arranged, this position just should by solar radiation to), at the flat board (the square flat board that for example the wide 3mm of 30cm is thick) of carbon steel material of this position of sound production (for example No. 45 carbon steels); Be called reference plate; Reference plate can not contact with ground, and reference plate overhead distance is not less than 1.5 meters, and reference plate can place the top of the wooden thermometer screen that meets meteorology temperature measurement requirement; The one side of this reference plate on the sunny side; (for example, in the time of on the Northern Hemisphere, sunny slope faces up towards south to be called sunny slope; Full daytime is all by sunshine; Sunny slope should have the suitable gradient to make snow not accumulate or behind snow, clear up sunny slope), the sunny slope of reference plate is coarse and (helping accepting solar irradiation) dark color, the sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day; The non-sunny slope of reference plate is covered with insulation material (the for example thick lime carbonate insulation material of 5mm), real-time monitoring record is obtained the temperature of the sunny slope of reference plate.
The b step; Monitoring in real time (can use conventional thermometry to measure; For example use thermal resistance to measure; For example every at a distance from temperature data of 10 minutes survey records) write down R the Cable Structure surface temperature measured data that obtains above-mentioned R Cable Structure surface point, monitoring in real time simultaneously (can be used conventional thermometry to measure, for example use the thermal resistance measurement; For example every at a distance from temperature data of 10 minutes survey records) obtain Cable Structure that the front defines temperature profile data along thickness; Monitoring in real time simultaneously (can use conventional thermometry to measure, for example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature, be for example every at a distance from temperature data of 10 minutes survey records) record obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; (can use conventional thermometry to measure through real-time monitoring; For example in the wooden thermometer screen that meets meteorology temperature measurement requirement, lay thermal resistance and measure temperature; For example every at a distance from temperature data of 10 minutes survey records) record obtain being carved at sunrise the same day sunrise next day constantly the Cable Structure between back 30 minutes belong to the temperature measured data sequence of environment; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T Emax(the temperature measured data sequence that for example earlier Cable Structure is belonged to environment carries out curve fitting through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment; Then through ask curve to the derivative of time or through ask with numerical method on the curve each corresponding to the point of survey record data time to the change of time rate) temperature that obtains Cable Structure place environment is about the change of time rate, this rate of change is also along with the time changes; (can use conventional thermometry to measure through real-time monitoring; For example use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark; For example every at a distance from temperature data of 10 minutes survey records) obtain being carved at sunrise the same day measured data sequence of temperature of the sunny slope of the reference plate between back 30 minutes of the moment of sunrise next day; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T Pmax(can use conventional thermometry to measure through real-time monitoring; For example use thermal resistance to measure the Cable Structure surface point; For example every at a distance from temperature data of 10 minutes survey records) record obtains being carved at sunrise the same day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes of the moment of sunrise next day; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; Have R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise, maximal value wherein is designated as Δ T Smax(for example earlier each Cable Structure surface temperature measured data sequence is carried out curve fitting through conventional mathematical computations by each Cable Structure surface temperature measured data sequence; Then through ask curve to the derivative of time or through ask with numerical method on the curve each corresponding to the point of survey record data time to the change of time rate) temperature that obtains each Cable Structure surface point is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes.Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " at the sea level elevation place that each is chosen; The absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
In the c step, measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, Δ T PmaxWith Δ T SmaX is not more than 5 degrees centigrade; Second b condition that must satisfy be the same day be carved at sunrise sunrise next day constantly between back 30 minutes during this period of time in, measure the Δ T that calculates in front EmaxBe not more than with reference to temperature difference per day Δ T r, and measure the Δ T that calculates in front PmaxDeduct 2 degrees centigrade and be not more than Δ T Emax, and measure the Δ T that calculates in front SmaxBe not more than Δ T PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T TmaxBe not more than 1 degree centigrade.This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly the moment of the acquisition Cable Structure steady temperature data of this method constantly; Then; According to the Cable Structure heat transfer characteristic; Utilize to obtain R Cable Structure surface temperature measured data and " HBE Cable Structure is along the thickness temperature measured data " in the moment of Cable Structure steady temperature data; Utilize the thermal conduction study computation model (for example finite element model) of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer (for example finite element method) in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; The accounting temperature that also comprises HBE " measuring the point of Cable Structure " that Cable Structure is selected in front along the temperature profile data of thickness; The accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data "; When R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal, the temperature profile data in the Cable Structure in moment of acquisition Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method; This moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T hThe numerical value that obtains is for Δ T is got in convenient narration hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure.
Second step: set up initial Mechanics Calculation benchmark model A o
In Cable Structure completion; Perhaps before setting up health monitoring (damaged cable identification) system; Calculating " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " measurement (can use conventional thermometry to measure; For example use thermal resistance to measure), this moment " Cable Structure steady temperature data " are used vector T oExpression is called initial Cable Structure steady temperature data vector T oObtain T in actual measurement oThe time, just at the synchronization in the moment that obtains initial Cable Structure steady temperature data vector, use conventional method directly to measure the initial value of all monitored amounts that calculate Cable Structure, form monitored amount initial value vector C o
Can be specifically in this method according to the synchronization of following method in the moment that obtains so-and-so Cable Structure steady temperature data vector such as (for example initial or current); Use so-and-so method measurement to calculate the data of the monitored amount of so-and-so measured amount (for example all monitored amounts of Cable Structure): in survey record temperature (comprising that Cable Structure belongs to the temperature and the Cable Structure surface temperature of the sunny slope of the temperature of environment, reference plate); For example every at a distance from temperature of 10 minutes survey records, equally also every so simultaneously at a distance from 10 minutes the monitored amount of so-and-so measured amount of survey record (for example all monitored amounts of Cable Structure) data.In case confirmed to obtain the moment of Cable Structure steady temperature data; Data with the monitored amount of so-and-so measured amount (for example all monitored amounts of Cable Structure) of the moment synchronization that obtains Cable Structure steady temperature data just are called the synchronization in the moment that obtains Cable Structure steady temperature data so, use so-and-so method to measure the data of the monitored amount of so-and-so measured amount that computing method obtain.
Use conventional method (consult reference materials or survey) to obtain the temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all supporting ropes, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all supporting ropes are in free state that Suo Li is 0 o'clock a length, the weight of cross-sectional area during in free state and the unit length during in free state; And the temperature of all supporting ropes when obtaining these three kinds of data; Utilize the temperature variant physical function parameter and the mechanical property parameters of all supporting ropes on this basis, obtain all supporting ropes at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that length, the Suo Li of 0 o'clock all supporting rope is that 0 o'clock all cross-sectional area and Suo Li that supports ropes is the weight of the unit length of 0 o'clock all supporting rope, forms the initial drift vector l of supporting rope successively o, initial free cross-sectional area vector A oWeight vector ω with initial free unit length o, the initial drift vector l of supporting rope o, initial free cross-sectional area vector A oWeight vector ω with initial free unit length oThe coding rule and initial rope force vector F of element oThe coding rule of element identical.Calculate initial Cable Structure steady temperature data vector T in actual measurement oThe time, just at the synchronization in the moment that obtains Cable Structure steady temperature data, use the conventional method actual measurement to calculate the actual measurement computational data of Cable Structure.The actual measurement computational data of Cable Structure comprises that the Non-Destructive Testing data that support rope etc. can express measured datas such as the initial geometric data of data, Cable Structure of the health status of rope, rope force data, draw-bar pull data, initial Cable Structure support coordinate data, Cable Structure bearing initial angle displacement measurement data, Cable Structure modal data, Cable Structure strain data, Cable Structure angle measurement data, Cable Structure volume coordinate measurement data.Initial Cable Structure support coordinate data refer to the support coordinate data under the Cable Structure design point, and Cable Structure bearing initial angle displacement measurement data refer to setting up initial Mechanics Calculation benchmark model A oThe time, the Cable Structure bearing is with respect to the angular displacement that bearing took place under the Cable Structure design point.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point on the Cable Structure, and purpose is to confirm according to these coordinate datas the geometric properties of Cable Structure.As far as cable-stayed bridge, the spatial data that initial geometric data can be the end points of all ropes adds the spatial data of some points on the bridge two ends, so-called bridge type data that Here it is.The Non-Destructive Testing data of utilization supporting rope etc. can be expressed the data and the Cable Structure bearing initial angle displacement measurement data of the health status of supporting rope and set up by evaluation object initial damage vector d o(shown in (1)) representes that with do Cable Structure is (with initial Mechanics Calculation benchmark model A oThe expression) by the initial health of evaluation object.If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, can think that perhaps the Cable Structure original state is a not damaged when not having relaxed state, vectorial d oIn get 0 with each element numerical value of supporting Suo Xiangguan, if there are not Cable Structure bearing initial angle displacement measurement data or can think that the displacement of Cable Structure bearing initial angle is at 0 o'clock, vectorial d oIn each element numerical value relevant with the angular displacement of Cable Structure bearing get 0.Utilize the measured data of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, the Non-Destructive Testing data of supporting rope, temperature variant physics and the mechanical property parameters and the initial Cable Structure steady temperature data vector T of the employed various materials of Cable Structure o, utilize mechanics method (for example finite element method) to count " Cable Structure steady temperature data " and set up initial Mechanics Calculation benchmark model A o
No matter which kind of method to obtain initial Mechanics Calculation benchmark model A with o, counting " Cable Structure steady temperature data " (is initial Cable Structure steady temperature data vector T o), based on A oThe Cable Structure computational data that calculates must be very near its measured data, and error generally must not be greater than 5%.Can guarantee to utilize A like this oSuo Li computational data, strain computational data, Cable Structure shape computational data and displacement computational data under the analog case of calculating gained, Cable Structure angle-data, Cable Structure spatial data etc., the measured data when truly taking place near institute's analog case reliably.Model A oIn used by evaluation object initial damage vector d by the health status of evaluation object oExpression, Cable Structure steady temperature data are with initial Cable Structure steady temperature data vector T oExpression.Because based on A oThe evaluation that calculates all monitored amounts is very near the initial value (actual measurement obtains) of all monitored amounts, so also can be used in A oThe basis on, carry out Mechanics Calculation obtains, A oThe evaluation of each monitored amount form monitored amount initial value vector C oCorresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A oUsed by the evaluation object health status by evaluation object initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression.T oAnd d oBe A oParameter, C oAlso can be by A oMechanics Calculation result form.
The 3rd step: in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning i o, A oAnd A i oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; What need during the i time circulation beginning is designated as d by the current initial damage vector of evaluation object i o, d i oCable Structure A when representing this time circulation beginning i oBy the health status of evaluation object, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i oAnd d i oBe A i oCharacterisitic parameter; C i oBy A i oMechanics Calculation result form; During circulation beginning for the first time, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring circulation beginning for the first time, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring circulation beginning for the first time, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring circulation beginning for the first time, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
The 4th step: the hardware components of pass line structural healthy monitoring system.Hardware components comprises at least: monitored amount monitoring system (for example containing angle measurement system, signal conditioner etc.), Cable Structure temperature monitoring system (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), signal (data) collector, the computing machine and the panalarm of communicating by letter.Each monitored amount, each temperature all must arrive by monitored system monitoring, and monitoring system is transferred to signal (data) collector with the signal that monitors; Signal is delivered to computing machine through signal picker; Computing machine then be responsible for the operation Cable Structure by the health monitoring software of evaluation object, comprise the signal that the transmission of tracer signal collector comes; When monitoring when being changed by the evaluation object health status, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 5th step: establishment and problem rope and the laddering recognition methods system software of bearing angular displacement of installation and operation temperature variation angle monitor on computers, this software will be accomplished functions such as monitoring that this method " the problem rope of temperature variation angle monitor and the laddering recognition methods of bearing angular displacement " required by task wants, record, control, storage, calculating, notice, warning (be in this practical implementation method all can with the work of computing machine completion).
The 6th step: the step begins the circulation running thus; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T i, vector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement corresponding one by one; In the actual measurement vector T iThe time, just obtaining current cable structure steady temperature data vector T iThe synchronization in the moment, actual measurement obtains the currency of all monitored amounts in the Cable Structure, all these numerical value are formed the current numerical value vector of monitored amount C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains all M in the Cable Structure 1The rope force data of root supporting rope, all these rope force datas are formed current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement calculates all M 1The volume coordinate of two supporting end points of root supporting rope, the volume coordinate of two supporting end points difference of component in the horizontal direction are exactly two supporting end points horizontal ranges, all M 1Two supporting end points horizontal range data of root supporting rope are formed current supporting rope two supporting end points horizontal range vector l t x, current supporting rope two supporting end points horizontal range vector l t xThe coding rule and initial rope force vector F of element oThe coding rule of element identical.
The 7th step: obtaining current cable structure steady temperature data vector T iAfter, compare T iAnd T i oIf, T iEqual T i o, then need be to A i oAnd T i oUpgrade, otherwise need be to current initial Mechanics Calculation benchmark model A i o, current initial Cable Structure steady temperature data vector T i oWith the current initial value vector of monitored amount C i oUpgrade, and by the current initial damage vector of evaluation object d i oRemain unchanged, update method follows these steps to a and carries out to step c:
A. calculate T iWith T oPoor, T iWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T iWith T oThe difference represent that with steady temperature change vector S S equals T iDeduct T o, S representes the variation of Cable Structure steady temperature data.
B. to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A oIn Cable Structure apply the current initial Mechanics Calculation benchmark model A that obtains upgrading after the temperature variation i o, upgrade A i oThe time, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, promptly upgraded T i o, so just obtained correctly corresponding to A i oT i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status by evaluation object use by the current initial damage of evaluation object vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A through Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C i o
The 8th step: at current initial Mechanics Calculation benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step a to steps d, set up unit damage monitored numerical quantity unit change matrix Δ C through calculating iWith unit damage or the angular displacement vector D of unit i u
A. the i time when beginning circulation, directly b obtains Δ C to the listed method of steps d set by step iAnd D i uAt other constantly, when in the 7th step to A i oAfter upgrading, b regains Δ C to the listed method of steps d set by step iAnd D i uIf, the 7th the step in not to A i oUpgrade, then directly changing for the 9th step herein over to carries out follow-up work.
B. at current initial Mechanics Calculation benchmark model A i oThe basis on carry out the several times Mechanics Calculation, vectorial d i oExpression A i oBy the health status of evaluation object, equal all on the calculation times numerical value by the quantity N of evaluation object, have N evaluation object that N calculating is just arranged; Calculate each time hypothesis have only one by evaluation object at vectorial d i oUnit damage or unit angular displacement are taken place in that representes on the basis of the health status of evaluation object, concrete, if should be a supporting rope in the cable system by evaluation object, so just hypothesis be somebody's turn to do the supporting rope at vectorial d i oOn the basis of the existing damage of this supporting rope of expression unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) is arranged again; If should be by evaluation object be the angular displacement component of a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d i oUnit angular displacement (for example get 100,000/radian, 2/100000ths radians, 3/100000ths radians etc. and be the unit angular displacement) takes place on the basis of the existing angular displacement of this bearing of expression again, uses D i UkWrite down this unit damage or unit angular displacement, wherein k represent to take place unit damage or unit angular displacement by the numbering of evaluation object, D i UkBe unit damage or the angular displacement vector D of unit i uAn element, unit damage or unit angular displacement vector D i uThe coding rule and the vectorial d of element oThe coding rule of element identical; Occur in calculating each time unit damage or unit angular displacement be different from by evaluation object other time occur in calculating unit damage or unit angular displacement by evaluation object; Calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount calculation current vector; When hypothesis k has unit damage or unit angular displacement by evaluation object, available C i Tk" the monitored amount calculation current vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C i TkDefinition mode and vectorial C oDefinition mode identical, C i TkElement and C oElement corresponding one by one.
C. the vectorial C that calculates each time i TkDeduct vectorial C i oObtain a vector, the unit damage or the angular displacement numerical value D of unit of supposition during each element that again should vector all calculates divided by this i UkAfter obtain " the unit change of a monitored amount vector δ C i k"; There is N N " the unit change vector of monitored amount " just to be arranged by evaluation object.
D. individual by the coding rule of evaluation object by this N " the unit change vector of monitored amount " according to N, form " the unit damage monitored numerical quantity unit change matrix Δ C that the N row are arranged successively i"; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to the different unit change amplitude of same monitored amount when difference is increased unit damage or unit angular displacement by evaluation object; Unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule and the vectorial d of row oThe coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule of coding rule and M monitored amount of row identical.
The 9th step: set up linear relationship error vector e iWith vectorial g iUtilize data (" the current initial value vector of the monitored amount C of front i o", " unit damage monitored numerical quantity unit change matrix Δ C i"), when the 8th step calculated each time, promptly calculate each time hypothesis have only in by evaluation object one by the increase unit damage of evaluation object or unit angular displacement in; when hypothesis k (k=1,2,3; ...; N) individual when being increased unit damage or unit angular displacement by evaluation object, calculate each time and form a damage vector, use d i TkExpression should damage vector, and corresponding monitored amount calculation current vector is C i Tk(referring to the 8th step), damage vectorial d i TkElement number equal by the quantity of evaluation object vectorial d i TkAll elements in have only the numerical value of an element to get to calculate each time in hypothesis increase unit damage or unit angular displacement by the unit damage of evaluation object or unit value of angular displacement, d i TkThe numerical value of other element get 0, what that was not that numbering and the supposition of 0 element increase unit damage or unit angular displacement is identical by the corresponding relation of evaluation object, with the element of the same numberings of other vectors with this corresponding relation by evaluation object; d i TkWith by evaluation object initial damage vector d oThe element coding rule identical, d i TkElement and d oElement be one-to-one relationship.With C i Tk, C i o, Δ C i, d i TkBring formula (34) into, obtain a linear relationship error vector e i k, calculate a linear relationship error vector e each time i ke i kSubscript k represent k (k=1,2,3 ..., N) individual by evaluation object increase unit damage or unit angular displacement.There is N N calculating just to be arranged, N linear relationship error vector e just arranged by evaluation object i k, with this N linear relationship error vector e i kObtain a vector after the addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e iVector g iEqual final error vector e iWith vectorial g iBe kept on the hard disc of computer of operation health monitoring systems software, supply health monitoring systems software to use.
e k i = abs ( Δ C i · d tk i - C tk i + C o i ) - - - ( 34 )
The tenth step: define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal by the quantity of evaluation object, d i cAnd d iElement and be one-to-one relationship between the evaluation object, d i cAnd d iElement numerical value represent corresponding by the degree of injury of evaluation object or bearing angular displacement, d i cAnd d iWith by evaluation object initial damage vector d oThe element coding rule identical, d i cElement, d iElement and d oElement be one-to-one relationship.
The 11 step: according to the current numerical value vector of monitored amount C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity unit change matrix Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula (11), calculates the vectorial d of current name damage according to multi-objective optimization algorithm i cNoninferior solution, just have reasonable error but can be more exactly from all ropes, confirm the position of damaged cable and separating of nominal degree of injury thereof.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted S Um Method), goal programming method (Goal Attainment Method) or the like based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is that example provides and finds the solution the vectorial d of current name damage with the goal programming method only i cProcess, the concrete implementation procedure of other algorithm can realize according to the requirement of its specific algorithm in a similar fashion.
According to the goal programming method, formula (11) can transform the multi-objective optimization question shown in an accepted way of doing sth (35) and the formula (36), and γ is a real number in the formula (35), and R is a real number field, and area of space Ω has limited vectorial d i cSpan (the present embodiment requirements vector d of each element i cEach element be not less than 0, be not more than 1).The meaning of formula (35) is to seek the real number γ of a minimum, makes formula (36) be met.G (d in the formula (36) i c) by formula (37) definition, G (d in the product representation formula (36) of weighing vector W and γ in the formula (36) i c) and vectorial g iBetween the deviation that allows, g iDefinition referring to formula (17), its value the 9th the step calculate.During actual computation the vector W can with vectorial g iIdentical.The concrete programming of goal programming method realizes having had universal program directly to adopt.Use the goal programming method just can damage vectorial d in the hope of current name i c
minimize γ
(35)
γ∈R, d c i ∈ Ω
G ( d c i ) - Wγ ≤ g i - - - ( 36 )
G ( d c i ) = abs ( Δ C i · d c i - C i + C o i ) - - - ( 37 )
The 12 step: foundation is by the current actual damage vector of evaluation object d iThe definition (seeing formula (19)) of definition (seeing formula (18)) and its element calculate current actual damage vector d iEach element, thereby can be by d iConfirm by the health status of evaluation object.By the current actual damage vector of evaluation object d iK element d i kRepresent that k is individual by the current actual health status of evaluation object in the i time circulation, if should be supporting rope, a d so in the cable system by evaluation object i kRepresent its current actual damage, d i kBe to represent that its corresponding supporting rope did not have health problem, d at 0 o'clock i kNumerical value is not to represent that the supporting rope that its corresponding supporting rope is the unsoundness problem, the supporting Suo Keneng of unsoundness problem were slack lines, also possibly are damaged cable at 0 o'clock, its numerical response the degree of lax or damage; If should be angular displacement component, a d so of a bearing by evaluation object i kRepresent its current actual displacement angle numerical value; So according to by the current actual damage of evaluation object vector d iCan define which the supporting rope unsoundness problem and order of severity thereof, define which bearing angular displacement and numerical value thereof have taken place; So far realized the bearing angular displacement identification of Cable Structure.
Will be by the current actual damage vector of evaluation object d iIn with the supporting Suo Xiangguan M 1Individual element takes out, and forms supporting Suo Dangqian actual damage vector d Ic, supporting Suo Dangqian actual damage vector d IcThe coding rule and initial rope force vector F of element oThe coding rule of element identical.Supporting Suo Dangqian actual damage vector d IcH element representation Cable Structure in the current actual damage amount of h root supporting rope, h=1,2,3 ...., M 1Current actual damage vector d IcIn numerical value be not 0 element supporting rope corresponding to the unsoundness problem, from the supporting rope of these unsoundness problems, identify damaged cable, remaining is exactly slack line.The method of differentiating is varied; Can pass through the protective seam of the supporting rope of removal unsoundness problem; To the visual discriminating of supporting Suo Jinhang; Perhaps carry out visual discriminating by optical imaging apparatus, also can be through lossless detection method to supporting rope impaired discriminating the whether, UT (Ultrasonic Testing) is exactly a kind of present widely used lossless detection method.Those do not find that supporting rope that damage, the unsoundness problem is exactly that lax rope has taken place to differentiate the back, need adjust the rope of Suo Li exactly, are exactly slack line, and these ropes that need adjust Suo Li are at supporting Suo Dangqian actual damage vector d IcIn pairing element numerical value (for example one of them element can be used d Ic hExpression) slack line has just been confirmed in expression and these degree of injury that supports the relax level mechanics equivalence of ropes thus, and the computing method of concrete slack are explained below.According to supporting Suo Dangqian actual damage vector d Ic, from the supporting rope of unsoundness problem, identify slack line after, remaining is exactly damaged cable, these damaged cables are at supporting Suo Dangqian actual damage vector d IcThe numerical value of the element of middle correspondence is just represented its degree of injury; The numerical value of corresponding element is to represent that this supporting rope thoroughly lost load-bearing capacity at 100% o'clock; Represent the load-bearing capacity of this supporting rope forfeiture corresponding proportion in the time of between 0 and 100%, so far just identified damaged cable and degree of injury thereof.
The 13 step: be utilized in current cable structure steady temperature data vector T iThe supporting Suo Dangqian actual damage vector d that obtains in the 12 step under the condition CaObtain slack line and the current actual equivalent damage degree equivalence of its relax level mechanics, be utilized in that the 6th step obtained at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current supporting rope two supporting end points horizontal range vector l t x, be utilized in that second step obtained at initial Cable Structure steady temperature data vector T oThe initial drift vector l of the supporting rope under the condition o, initial free cross-sectional area vector A oWeight vector ω with initial free unit length o, utilize current cable structure steady temperature data vector T iThe supporting Suo Dangqian steady temperature data of expression, be utilized in that second step obtained at initial Cable Structure steady temperature data vector T oThe supporting rope initial steady state temperature data of expression; Be utilized in the temperature variant physics and the mechanical property parameters of the employed various materials of Cable Structure of second step acquisition; Count the influence of temperature variation to supporting rope physics, mechanics and geometric parameter; Through with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical.When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa.Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of supporting rope drift, has just confirmed the long adjustment amount of rope of the supporting rope that those need adjust Suo Li.Particularly can be in the hope of the relax level (being the long adjustment amount of rope) of these ropes according to formula (29) or formula (30).So just realized the lax identification of supporting rope.So far damaged cable, slack line and bearing angular displacement have just all been discerned.
The 14 step: the computing machine in the health monitoring systems regularly generates cable system health condition form automatically or by the personnel operation health monitoring systems.
The 15 step: under specified requirements, the automatic operation communication panalarm of the computing machine in the health monitoring systems to monitor staff, owner and (or) personnel of appointment report to the police.
The 16 step: set up mark vector B according to formula (31) i, formula (32) has provided mark vector B iThe definition of k element; If mark vector B iElement be 0 entirely, then got back to for the 6th step and proceed health monitoring and calculating cable system; If mark vector B iElement be not 0 entirely, then accomplish subsequent step after, get into circulation next time.
The 17 the step: according to formula (33) calculate next time (promptly the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation I+1 oEach element d I+1 Ok(k=1,2,3 ..., N); The second, at initial Mechanics Calculation benchmark model A oThe basis on, to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, making by the health status of evaluation object again is d I+1 oAfter obtain be exactly next time, promptly the i+1 time (i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation I+1Next time (promptly the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual T i oObtain A I+1, d I+1 oAnd T I+1 oAfter, obtain A through Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value are formed next time, the vectorial C of the current initial value of required monitored amount that promptly circulates for the i+1 time I+1 o
The 18 step: got back to for the 6th step, beginning goes on foot the 17 circulation by the 6th.

Claims (1)

1. the problem rope of a temperature variation angle monitor and the laddering recognition methods of bearing angular displacement is characterized in that said method comprises:
A. for for the purpose of narration is convenient, this method is unified claims that the supporting rope of being assessed is that establishing the quantity of the supporting rope of being assessed and the quantity sum of bearing angular displacement component is N, is N by the quantity of evaluation object promptly by evaluation object with the bearing angular displacement component; Confirm that by the coding rule of evaluation object with all being numbered by evaluation object in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; This method is represented this numbering with variable k, k=1, and 2,3 ..., N; If total M in the cable system 1Root supporting rope, Cable Structure rope force data comprises this M 1The Suo Li of root supporting rope, obviously M 1Less than by the quantity N of evaluation object; Confirm the measured point of appointment, give all specified point numberings; Confirm the measured straight line of each measurement point, gave the measured straight line numbering of all appointments; Confirm the measured angle coordinate component of each measured straight line, give all measured angle coordinate component numberings; Above-mentioned numbering will be used to generate the vector sum matrix in subsequent step; " the whole monitored angle-data of Cable Structure " is made up of above-mentioned all measured angle coordinate components; For simplicity, in the method " the monitored angle-data of Cable Structure " abbreviated as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M must not be less than N; Must not be greater than 30 minutes in this method to the time interval between any twice measurement of same amount monitoring in real time, the moment of survey record data is called the physical record data constantly;
B. this method definition " the temperature survey calculating method of the Cable Structure of this method " set by step b1 to b3 carry out;
B1: inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of Cable Structure composition material and Cable Structure environment of living in; Utilize the geometry measured data of design drawing, as-constructed drawing and the Cable Structure of Cable Structure, utilize these data and parameter to set up the thermal conduction study computation model of Cable Structure; Inquiry Cable Structure location is no less than the meteorological data in recent years in 2 years; Cloudy quantity in statistics obtains during this period of time is designated as T cloudy day; One be called the cloudy day all day with what can not see the sun daytime in the method; Statistics obtain in T cloudy day each cloudy day 0 up to the sunrise next day highest temperature and the lowest temperature between back 30 minutes constantly, sunrise is meant the sunrise moment on the meteorology that base area revolutions and revolution rule confirm constantly, does not represent necessarily can see the sun same day; The sunrise that can inquire about data or calculate each required day through conventional meteorology constantly; Each cloudy day 0 up to next day sunrise constantly the highest temperature between back 30 minutes deduct the maximum temperature difference that the lowest temperature is called this cloudy daily temperature, T cloudy day arranged, the maximum temperature difference of T cloudy daily temperature is just arranged; Get maximal value in the maximum temperature difference of daily temperature at T cloudy day for reference to temperature difference per day, be designated as Δ T with reference to temperature difference per day rInquiry Cable Structure location and height above sea level interval, place be no less than temperature that meteorological data in recent years or the actual measurement in 2 years obtain Cable Structure environment of living in time with change of elevation data and Changing Pattern, calculate Cable Structure location and height above sea level interval, place and be no less than the temperature of Cable Structure environment of living in recent years in 2 years about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; On the surface of Cable Structure, get " R Cable Structure surface point "; Getting the concrete principle of " R Cable Structure surface point " narrates in step b3; The back will obtain the temperature of this R Cable Structure surface point through actual measurement; Claim that the temperature data that actual measurement obtains is " R Cable Structure surface temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain the temperature of this R Cable Structure surface point through Calculation of Heat Transfer, just claim that the temperature data that calculates is " R Cable Structure surface temperature computational data "; From the residing minimum height above sea level of Cable Structure to the highest height above sea level; On Cable Structure, be uniformly distributed with to choose and be no less than three different altitude above sea level; At each sea level elevation place that chooses, choose two points at least at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions of choosing are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness "; It is crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness to measure Cable Structure; In in the shade the outer normal direction of the measurement Cable Structure of choosing along sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, the direction along each measurement Cable Structure along the Temperature Distribution of wall thickness is uniformly distributed with to choose in Cable Structure and is no less than three points, and is special; Measure Cable Structure for the supporting rope along each and only get a point along the direction of the Temperature Distribution of wall thickness; Promptly only measure the temperature of the surface point of supporting rope, measure all and be selected temperature a little, the temperature that records is called " Cable Structure is along the temperature profile data of thickness "; Wherein edge and same " intersection on surface level and Cable Structure surface " crossing, " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness " measure " Cable Structure is along the temperature profile data of thickness " that obtain; Be called " identical sea level elevation Cable Structure is along the temperature profile data of thickness " in the method, establish and chosen H different altitude above sea level, at each sea level elevation place; Chosen B the direction of measuring Cable Structure along the Temperature Distribution of wall thickness; Measure Cable Structure along each and in Cable Structure, chosen E point along the direction of the Temperature Distribution of wall thickness, wherein H and E are not less than 3, and B is not less than 2; Special; E equals 1 for the supporting rope, and that " measures the point of Cable Structure along the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, and the back will obtain the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " through actual measurement; Claim that the temperature data that actual measurement obtains is " HBE Cable Structure is along thickness temperature measured data "; If utilize the thermal conduction study computation model of Cable Structure, obtain this HBE the temperature of measuring Cable Structure along the point of the temperature profile data of thickness through Calculation of Heat Transfer, just claim that the temperature data that calculates is " HBE Cable Structure is along thickness temperature computation data "; Will be in this method " at the number temperature profile data of each sea level elevation place that chooses " identical sea level elevation Cable Structure is along the temperature profile data of thickness "; Measure temperature in the Cable Structure location according to meteorology and require to choose a position, will obtain meeting the temperature that meteorology is measured the Cable Structure place environment of temperature requirement in this position actual measurement; The place of blocking chooses a position in the on-site spacious nothing of Cable Structure; This position should can both obtain in each day of the whole year this ground sunshine of fullest of getable this day, the flat board at a carbon steel material of this position of sound production is called reference plate; Reference plate can not contact with ground; Reference plate overhead distance is not less than 1.5 meters, and the one side of this reference plate is called sunny slope on the sunny side; The sunny slope of reference plate is coarse and dark color; The sunny slope of reference plate should can both obtain in each day of the whole year one flat plate on this ground sunshine of fullest of getable this day, the non-sunny slope of reference plate is covered with insulation material, monitoring is in real time obtained the temperature of the sunny slope of reference plate;
B2: monitoring in real time obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point; Monitoring in real time simultaneously obtains the temperature profile data of the Cable Structure of front definition along thickness, and monitoring in real time simultaneously obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement; Obtain being carved at sunrise the same day sunrise next day temperature measured data sequence of the place of the Cable Structure between back 30 minutes environment constantly through real-time monitoring; The temperature measured data sequence of Cable Structure place environment by be carved at sunrise the same day sunrise next day constantly the temperature measured data of the place of the Cable Structure between back 30 minutes environment according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the temperature measured data sequence of Cable Structure place environment; Deduct with the maximum temperature in the temperature measured data sequence of Cable Structure place environment and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains Cable Structure place environment; Be called the environment maximum temperature difference, be designated as Δ T EmaxThe temperature that obtains Cable Structure place environment through conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment is about the change of time rate, and this rate of change is also along with the time changes; Obtain being carved at sunrise the same day sunrise next day measured data sequence of the temperature of the sunny slope of the reference plate between back 30 minutes constantly through real-time monitoring; The measured data sequence of the temperature of the sunny slope of reference plate by be carved at sunrise the same day next day sunrise constantly the measured data of the temperature of the sunny slope of the reference plate between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in the measured data sequence of temperature of sunny slope of reference plate; Deduct with the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day of temperature that minimum temperature obtains the sunny slope of reference plate; Be called the reference plate maximum temperature difference, be designated as Δ T PmaxObtain being carved at sunrise the same day sunrise next day Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between back 30 minutes constantly through real-time monitoring; There is R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged; Each Cable Structure surface temperature measured data sequence by be carved at sunrise on same day of a Cable Structure surface point sunrise next day constantly the Cable Structure surface temperature measured data between back 30 minutes according to the time order and function series arrangement; Find maximum temperature and minimum temperature in each Cable Structure surface temperature measured data sequence; Deduct with the maximum temperature in each Cable Structure surface temperature measured data sequence and to be carved into the sunrise next day maximum temperature difference between back 30 minutes constantly at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point; There is R Cable Structure surface point just to have and be carved into the sunrise next day maximum temperature difference numerical value between back 30 minutes constantly R the same day at sunrise; Maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T SmaxThe temperature that obtains each Cable Structure surface point through conventional mathematical computations by each Cable Structure surface temperature measured data sequence is about the change of time rate, the temperature of each Cable Structure surface point about the change of time rate also along with the time changes; Through real-time monitoring obtain being carved at sunrise the same day sunrise next day constantly between back 30 minutes, behind synchronization, HBE " Cable Structure is along the temperature profile data of thickness "; Calculating amounts to maximum temperature and the difference of minimum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " at the sea level elevation place that each is chosen; The absolute value of this difference is called " identical sea level elevation place Cable Structure thickness direction maximum temperature difference "; Chosen H different altitude above sea level H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " just arranged; Claim that the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is " a Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
B3: measure and calculate acquisition Cable Structure steady temperature data; At first; Confirm to obtain the moment of Cable Structure steady temperature data; The relevant condition of the moment that obtains Cable Structure steady temperature data with decision has six; To be moment of obtaining Cable Structure steady temperature data be carved into sunrise next day constantly between back 30 minutes at sunset between the same day for first condition, sunset be meant constantly on the meteorology that base area revolutions and revolution rule confirm sunset constantly, the sunset that can inquire about data or calculate each required day through conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, reference plate maximum temperature difference Δ T PmaxWith Cable Structure surface maximum temperature difference Δ T SmaxAll be not more than 5 degrees centigrade; The b condition of second condition be the same day be carved at sunrise next day sunrise constantly between back 30 minutes during this period of time in, measure the environment maximum error Δ T that calculates in front EmaxBe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T PmaxBe not more than Δ T after deducting 2 degrees centigrade Emax, and Cable Structure surface maximum temperature difference Δ T SmaxBe not more than Δ T PmaxOne that only needs to satisfy in second a condition and the b condition just is called satisfied second condition; The 3rd condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 4th condition is in the moment that obtains Cable Structure steady temperature data, and the temperature of each the Cable Structure surface point in R Cable Structure surface point is not more than per hour 0.1 degree centigrade about the absolute value of change of time rate; The 5th condition is in the moment that obtains Cable Structure steady temperature data, and the Cable Structure surface temperature measured data of each the Cable Structure surface point in R Cable Structure surface point is the minimal value that was carved at sunrise the same day between back 30 minutes of the moment of sunrise next day; The 6th condition is at the moment that obtains Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T TmaxBe not more than 1 degree centigrade; This method is utilized above-mentioned six conditions; In following three kinds of moment any one is called the mathematics of Cable Structure steady temperature data " obtain constantly "; First kind of moment is first moment to the 5th condition of satisfying in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data "; Second kind of moment is the moment of only satisfying the 6th condition in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data ", and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " with the relevant condition of the moment of decision acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during physical record data in constantly constantly, the moment that obtains Cable Structure steady temperature data be exactly obtain Cable Structure steady temperature data mathematics constantly; If obtain the mathematics of Cable Structure steady temperature data and constantly is not any in constantly of physical record data in this method constantly, then get moment of mathematics those physical record data constantly that this method approaches to obtain Cable Structure steady temperature data most for obtaining the moment of Cable Structure steady temperature data; This method will be used in the amount of the moment survey record that obtains Cable Structure steady temperature data and carry out the relevant health monitoring analysis of Cable Structure; This method is approximate thinks that the Cable Structure temperature field in moment of obtaining Cable Structure steady temperature data is in stable state, and promptly this Cable Structure temperature does not constantly change in time, and this is exactly " obtaining the moment of Cable Structure steady temperature data " of this method constantly; Then; According to the Cable Structure heat transfer characteristic; Utilize " R the Cable Structure surface temperature measured data " and " HBE Cable Structure is along thickness temperature measured data " in the moment that obtains Cable Structure steady temperature data; Utilize the thermal conduction study computation model of Cable Structure; Obtain Temperature Distribution through conventional Calculation of Heat Transfer in the Cable Structure in the moment that obtains Cable Structure steady temperature data; This moment, calculated by stable state in the temperature field of Cable Structure; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates comprises the accounting temperature of R Cable Structure surface point on the Cable Structure; The accounting temperature of R Cable Structure surface point is called R Cable Structure stable state surface temperature computational data; Also comprise the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " that Cable Structure is selected in front, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure is along thickness temperature computation data ", when R Cable Structure surface temperature measured data and R Cable Structure stable state surface temperature computational data correspondent equal; And when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure is along thickness temperature computation data " correspondent equal; The temperature profile data in the Cable Structure in the moment that obtains Cable Structure steady temperature data that calculates is called " Cable Structure steady temperature data " in the method, and this moment " R Cable Structure surface temperature measured data " is called " R Cable Structure stable state surface temperature measured data ", and " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data "; When on the surface of Cable Structure, getting " R Cable Structure surface point "; The quantity of " R Cable Structure surface point " and necessary three conditions that satisfy that distribute; First condition is when the Cable Structure temperature field is in stable state; When on the Cable Structure surface arbitrarily the temperature of any be through " R Cable Structure surface point " in the Cable Structure surface on the observed temperature linear interpolation of the adjacent point in this arbitrfary point when obtaining, on the Cable Structure surface that linear interpolation obtains on the temperature of this arbitrfary point and the Cable Structure surface error of the actual temperature of this arbitrfary point be not more than 5%; The Cable Structure surface comprises supporting rope surface; Second condition is that the quantity at the point of same sea level elevation is not less than 4 in " R Cable Structure surface point ", and the point in same sea level elevation is uniformly distributed with along the Cable Structure surface in " R Cable Structure surface point "; Maximal value Δ h in the absolute value of the difference of the sea level elevation of all adjacent in twos Cable Structure surface points of " R Cable Structure surface point " coastal degree of lifting is not more than 0.2 ℃ divided by Δ T hThe numerical value that obtains is for Δ T is got in convenient narration hUnit be ℃/m that the unit of getting Δ h for convenient narration is m; The definition of the adjacent in twos Cable Structure surface point of " R Cable Structure surface point " coastal degree of lifting is meant when only considering sea level elevation; In " R Cable Structure surface point ", do not have a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point in twos; The 3rd condition is to inquire about or calculate the Cable Structure location and belong to the interval rule at sunshine of height above sea level by the meteorology routine; Again according to the geometric properties and the bearing data of Cable Structure; On Cable Structure, find the position of those surface points that receive the sunshine-duration fullest whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is a point in annual those surface points that receive the sunshine-duration fullest on the Cable Structure;
C. directly measure according to " the temperature survey calculating method of the Cable Structure of this method " and calculate the Cable Structure steady temperature data under the original state; Cable Structure steady temperature data under the original state are called initial Cable Structure steady temperature data, are designated as " initial Cable Structure steady temperature data vector T o"; Survey or consult reference materials and obtain the temperature variant physics and the mechanical property parameters of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all supporting ropes, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all supporting ropes are in free state that Suo Li is 0 o'clock a length, the weight of cross-sectional area during in free state and the unit length during in free state; And the temperature of all supporting ropes when obtaining these three kinds of data; Utilize the temperature variant physical function parameter and the mechanical property parameters of all supporting ropes on this basis, obtain all supporting ropes at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that length, the Suo Li of 0 o'clock all supporting rope is that 0 o'clock all cross-sectional area and Suo Li that supports ropes is the weight of the unit length of 0 o'clock all supporting rope; Form the initial drift vector of supporting rope, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule and initial rope force vector F of the element of the initial drift vector of supporting rope, the weight vector of the initial free unit length of initial free cross-sectional area vector sum oThe coding rule of element identical; Obtain T in actual measurement oThe time, just obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment; Directly measure the measured data that calculates initial Cable Structure, the measured data of initial Cable Structure comprises the Non-Destructive Testing data, Cable Structure bearing initial angle displacement measurement data, the initial value of all monitored amounts, the initial rope force data of all supporting ropes, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure support coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data of the health status of expressing the supporting rope; The initial value of all monitored amounts is formed monitored amount initial value vector C o, monitored amount initial value vector C oThe coding rule of coding rule and M monitored amount identical; Utilize the Non-Destructive Testing data and the Cable Structure bearing initial angle displacement measurement data of the health status that can express the supporting rope to set up by evaluation object initial damage vector d o, vectorial d oExpression is with initial Mechanics Calculation benchmark model A oThe expression Cable Structure by the initial health of evaluation object; By evaluation object initial damage vector d oElement number equal N, d oElement be one-to-one relationship by evaluation object, vectorial d oThe coding rule of element with identical by the coding rule of evaluation object; If d oSome elements corresponding be supporting rope, a d so in the cable system by evaluation object oThe numerical value of this element represent the initial damage degree of corresponding supporting rope; If the numerical value of this element is 0, represent that the pairing supporting rope of this element is intact, do not damage; If its numerical value is 100%; Represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, lost the load-bearing capacity of corresponding proportion if its numerical value between 0 and 100%, is then represented this supporting rope; If d oSome elements corresponding be some angular displacement components of some bearings, d so by evaluation object oThe numerical value of this element represent the initial value of this angular displacement component of this bearing; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, can think that perhaps the Cable Structure original state is a not damaged when not having relaxed state, vectorial d oIn get 0 with each element numerical value of supporting Suo Xiangguan, if there are not Cable Structure bearing initial angle displacement measurement data or can think that the displacement of Cable Structure bearing initial angle is at 0 o'clock, vectorial d oIn each element numerical value relevant with the angular displacement of Cable Structure bearing get 0; Initial Cable Structure support coordinate data refer to the support coordinate data under the Cable Structure design point, and Cable Structure bearing initial angle displacement measurement data refer to setting up initial Mechanics Calculation benchmark model A oThe time, the Cable Structure bearing is with respect to the angular displacement that bearing took place under the Cable Structure design point;
D. according to the temperature variant physics of the Non-Destructive Testing data of the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, supporting rope, Cable Structure bearing initial angle displacement measurement data, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure steady temperature data vector T oWith all Cable Structure data that obtain with preceding step, set up the initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " o, based on A oThe Cable Structure computational data that calculates must be very near its measured data, and difference therebetween must not be greater than 5%; Corresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A oUsed by the evaluation object health status by evaluation object initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression; T oAnd d oBe A oParameter, by A oThe initial value and the C of all monitored amounts of obtaining of Mechanics Calculation result oThe initial value of all monitored amounts of expression is identical, therefore also we can say C oBy A oMechanics Calculation result form A in the method o, C o, d oAnd T oBe constant;
E. in the method, alphabetical i is except the place of representing number of steps significantly, and alphabetical i only representes cycle index, i.e. the i time circulation; The current initial Mechanics Calculation benchmark model of Cable Structure that need set up or that set up is designated as current initial Mechanics Calculation benchmark model A during i time circulation beginning i o, A oAnd A i oCounted temperature parameter, can accounting temperature change mechanical property influence Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; What need during the i time circulation beginning is designated as d by the current initial damage vector of evaluation object i o, d i oCable Structure A when representing this time circulation beginning i oBy the health status of evaluation object, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i oAnd d i oBe A i oCharacterisitic parameter, C i oBy A i oMechanics Calculation result form; During circulation beginning for the first time, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring circulation beginning for the first time, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring circulation beginning for the first time, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring circulation beginning for the first time, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
F. go on foot the s circulation in step from getting into here by f; In Cable Structure military service process; Constantly survey the current data of calculating acquisition Cable Structure steady temperature data according to " the temperature survey calculating method of the Cable Structure of this method ", the current data of all " Cable Structure steady temperature data " is formed current cable structure steady temperature data vector T i, vector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement corresponding one by one; Obtain vector T in actual measurement iThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T iThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value are formed the current numerical value vector of monitored amount C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains all M in the Cable Structure 1The rope force data of root supporting rope, all these rope force datas are formed current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement calculates all M 1The volume coordinate of two supporting end points of root supporting rope; The volume coordinate of two the supporting end points difference of component in the horizontal direction is exactly two supporting end points horizontal ranges; Two supporting end points horizontal range data of all supporting ropes are formed current supporting rope two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current supporting rope two supporting end points horizontal range vectors oThe coding rule of element identical;
G. according to current cable structure steady temperature data vector T i, upgrade current initial Mechanics Calculation benchmark model A according to step g 1 to g3 i o, the current initial value of monitored amount vector C i oWith current initial Cable Structure steady temperature data vector T i o, and by the current initial damage vector of evaluation object d i oRemain unchanged;
G1. compare T iWith T i oIf, T iEqual T i o, A then i o, C i oAnd T i oRemain unchanged; Otherwise need follow these steps to A i oAnd T i oUpgrade;
G2. calculate T iWith T oPoor, T iWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T iWith T oThe difference represent that with steady temperature change vector S S equals T iDeduct T o, S representes the variation of Cable Structure steady temperature data;
G3. to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading i o, upgrade A i oThe time, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, promptly upgraded T i o, so just obtained correctly corresponding to A i oT i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status of rope use by the current initial damage vector of evaluation object d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A through Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C i o
H. at current initial Mechanics Calculation benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity unit change matrix Δ C through calculating iWith unit damage or the angular displacement vector D of unit i u
H1. the i time when beginning circulation, directly h2 obtains Δ C to the listed method of step h4 set by step iAnd D i uAt other constantly, when in step g to A i oAfter upgrading, h2 regains Δ C to the listed method of step h4 set by step iAnd D i uIf, in step g not to A i oUpgrade, then directly change step I herein over to and carry out follow-up work;
H2. at the current initial Mechanics Calculation benchmark model A of Cable Structure i oThe basis on carry out the several times Mechanics Calculation, equal all on the calculation times numerical value by the quantity N of evaluation object, have N evaluation object that N calculating is just arranged; According to by the coding rule of evaluation object, calculate successively; Calculating hypothesis each time has only one on the basis of original damage or angular displacement, to be increased unit damage or unit angular displacement again by evaluation object; Concrete; If should be a supporting rope in the cable system by evaluation object, so just this supporting rope of hypothesis increases unit damage again, if should be the angular displacement component of a direction of a bearing by evaluation object; Just suppose that this bearing increases the unit angular displacement again at this sense of displacement, uses D i UkWrite down the unit damage or the unit angular displacement of this increase, wherein k represent to increase unit damage or unit angular displacement by the numbering of evaluation object, D i UkBe unit damage or the angular displacement vector D of unit i uAn element, unit damage or unit angular displacement vector D i uThe coding rule and the vectorial d of element oThe coding rule of element identical; Increase again in calculating each time unit damage or unit angular displacement be different from by evaluation object other time increase again in calculating unit damage or unit angular displacement by evaluation object; Calculate the current calculated value that all utilizes mechanics method to calculate all monitored amounts of Cable Structure each time, the current calculated value of the monitored amount of all that calculate is each time formed a monitored amount calculation current vector; When hypothesis k increases unit damage or unit angular displacement again by evaluation object, use C i Tk" the monitored amount calculation current vector " that expression is corresponding; When in this step, giving each vectorial element numbering; Should use same coding rule with other vector in this method; To guarantee any element in each vector in this step,, expressed the relevant information of same monitored amount or same target with element in other vector, that numbering is identical; C i TkDefinition mode and vectorial C oDefinition mode identical, C i TkElement and C oElement corresponding one by one;
H3. the vectorial C that calculates each time i TkDeduct vectorial C i oObtain a vector, each element that again should vector obtains " unit change vector δ a C of monitored amount after all calculating the unit damage supposed or unit angular displacement numerical value divided by this i k"; There is N N " the unit change vector of monitored amount " just to be arranged by evaluation object;
H4. individual by the coding rule of evaluation object by this N " the unit change vector of monitored amount " according to N, form " the unit damage monitored numerical quantity unit change matrix Δ C that the N row are arranged successively i"; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to a monitored amount unit change vector; Unit damage monitored numerical quantity unit change matrix Δ C iEach row corresponding to the different unit change amplitude of same monitored amount when difference is increased unit damage or unit angular displacement by evaluation object; Unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule and the vectorial d of row oThe coding rule of element identical, unit damage monitored numerical quantity unit change matrix Δ C iThe coding rule of coding rule and M monitored amount of row identical;
I. define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal to support the quantity of rope, d i cAnd d iElement and be one-to-one relationship between the evaluation object, d i cAnd d iElement numerical value represent corresponding by the degree of injury of evaluation object or health status, d i cAnd d iWith by evaluation object initial damage vector d oThe element coding rule identical, d i cElement, d iElement and d oElement be one-to-one relationship;
J. according to the current numerical value vector of monitored amount C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity unit change matrix Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, removes d in the formula 1 i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + Δ C i · d c i Formula 1
K. the current actual damage vector d that utilizes formula 2 to express iK element d i kTogether by the current initial damage vector of evaluation object d i oK element d i OkWith the vectorial d of current name damage i cK element d i CkBetween relation, calculate current actual damage vector d iAll elements;
formula 2
K=1 in the formula 2,2,3 ..., N; d i kRepresent that k is individual by the current actual health status of evaluation object, d in the i time circulation i kBe to represent that k was not had health problem by evaluation object, d at 0 o'clock i kNumerical value be not represented in 0 o'clock k by evaluation object be the unsoundness problem by evaluation object, if should be supporting rope, a d so in the cable system by evaluation object i kThe order of severity of representing its current health problem, the supporting Suo Keneng of unsoundness problem is slack line, also possibly is damaged cable, d i kNumerical response the degree of lax or damage that should the supporting rope, if should be angular displacement component, a d so of a bearing by evaluation object i kRepresent its current actual displacement angle numerical value; So according to by the current actual damage of evaluation object vector d iCan define which the supporting rope unsoundness problem and order of severity thereof, define which bearing angular displacement and numerical value thereof have taken place;
L. will be by the current actual damage vector of evaluation object d iIn with M 1The M of root supporting Suo Xiangguan 1Individual element takes out, and forms supporting Suo Dangqian actual damage vector d Ca, supporting Suo Dangqian actual damage vector d CaThe coding rule and initial rope force vector F of element oThe coding rule of element identical; Supporting Suo Dangqian actual damage vector d CaH element representation Cable Structure in the current actual damage amount of h root supporting rope, h=1,2,3 ...., M 1Supporting Suo Dangqian actual damage vector d CaIn numerical value be not 0 element supporting rope corresponding to the unsoundness problem, from the supporting rope of these unsoundness problems, identify damaged cable, remaining is exactly slack line, supporting Suo Dangqian actual damage vector d CaIn with slack line corresponding to the element numerical expression be the current actual equivalent damage degree with slack line relax level mechanics equivalence;
M. be utilized in current cable structure steady temperature data vector T iSlack line under the condition, that identify in l step and with the vectorial d of supporting Suo Dangqian actual damage CaThese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanics equivalence, be utilized in that the f step obtains at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current supporting rope two supporting end points horizontal ranges vectors, be utilized in that the c step obtains at initial Cable Structure steady temperature data vector T oThe initial drift vector of the supporting rope under the condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, initial rope force vector F o, utilize current cable structure steady temperature data vector T tThe supporting Suo Dangqian steady temperature data of expression, be utilized in that c step obtains at initial Cable Structure steady temperature data vector T oThe supporting rope initial steady state temperature data of expression; Be utilized in the temperature variant physics and the mechanical property parameters of the employed various materials of Cable Structure of c step acquisition; Count the influence of temperature variation to supporting rope physics, mechanics and geometric parameter; Through with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after length overall identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two supporting ropes in Cable Structure is exactly identical, if after promptly replacing damaged cable with the slack line of equivalence, Cable Structure any variation can not take place, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of supporting rope drift, has just confirmed the long adjustment amount of rope of the supporting rope that those need adjust Suo Li; The lax identification and the damage identification of supporting rope have so just been realized; Institute's demand power is provided by current cable force vector F corresponding element during calculating;
N. try to achieve the vectorial d of current name damage i cAfter, set up mark vector B according to formula 3 i, formula 4 has provided mark vector B iThe definition of k element;
B i = B 1 i B 2 i . . . B k i . . . B N i T Formula 3
formula 4
Element B in the formula 4 i kBe mark vector B iK element, D i UkBe unit damage or the angular displacement vector D of unit i uK element, d i CkBy the vectorial d of the current name damage of evaluation object i cK element, they all represent k by the relevant information of evaluation object, k=1 in the formula 4,2,3 ..., N;
If mark vector B o. iElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B iElement be not 0 entirely, then get into next step, be step p;
P. according to formula 5 calculate next time, i.e. the i+1 time circulation be required by the current initial damage vector of evaluation object d I+1 oEach element;
formula 5
D in the formula 5 I+1 OkBe next time, i.e. the i+1 time circulation be required by the current initial damage vector of evaluation object d I+1 oK element, d i OkBe this, promptly the i time round-robin be by the current initial damage vector of evaluation object d i oK element, D i UkBe the i time round-robin unit damage or the angular displacement vector D of unit i uK element, B i kBe the i time round-robin mark vector B iK element, k=1 in the formula 5,2,3 ..., N;
Q. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual the current initial Cable Structure steady temperature data vector T of round-robin the i time i o
R. at initial Mechanics Calculation benchmark model A oThe basis on, to A oIn Cable Structure apply temperature variation, the numerical value of the temperature variation that applies is just taken from steady temperature change vector S, the health status that makes rope again is d I+1 oAfter obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation I+1Obtain A I+1After, obtain A through Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value are formed next time, the vectorial C of the current initial value of required monitored amount that promptly circulates for the i+1 time I+1 o
S. get back to step f, beginning is circulation next time.
CN201210173398.XA 2012-05-30 2012-05-30 Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change CN102706659B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210173398.XA CN102706659B (en) 2012-05-30 2012-05-30 Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210173398.XA CN102706659B (en) 2012-05-30 2012-05-30 Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change

Publications (2)

Publication Number Publication Date
CN102706659A true CN102706659A (en) 2012-10-03
CN102706659B CN102706659B (en) 2015-07-08

Family

ID=46899602

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210173398.XA CN102706659B (en) 2012-05-30 2012-05-30 Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change

Country Status (1)

Country Link
CN (1) CN102706659B (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604632A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load progressive recognition method on basis of generalized displacement angle monitoring
CN103616248A (en) * 2013-12-09 2014-03-05 东南大学 Progressive identifying method for problem cable concentrated loads and angular displacement based on angle monitoring
CN103630412A (en) * 2013-12-09 2014-03-12 东南大学 Method for progressively identifying concentrated loads of problematic cables by means of monitoring angles during angular displacement
CN103852299A (en) * 2014-03-10 2014-06-11 东南大学 Linear displacement space coordinate monitoring damaged cable load identification method
CN103852337A (en) * 2014-03-10 2014-06-11 东南大学 Strain monitoring error cable load angular displacement progressive identification method
CN103852332A (en) * 2014-03-10 2014-06-11 东南大学 Damaged cable load linear displacement progressive recognition method based on space coordinate monitoring
CN103852330A (en) * 2014-03-10 2014-06-11 东南大学 Method for recognizing damaged cables, loads and linear displacement based on cable force monitoring
CN103852300A (en) * 2014-03-10 2014-06-11 东南大学 Progressive recognition method for damaged cable loads through generalized displacement and angle monitoring
CN103884527A (en) * 2014-03-10 2014-06-25 东南大学 Method for recognizing load angular displacement and issue cables through angle monitoring
CN103913337A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement cable force monitoring defective cable load progressive identification method
CN103913322A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing fault cable and load based on angular displacement angle monitoring
CN103913329A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement hybrid monitoring damaged cable load identification method
CN103913332A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing damaged cable and load based on angular displacement strain monitoring
CN103913320A (en) * 2014-03-10 2014-07-09 东南大学 Method for progressively recognizing damaged cable, load and angular displacement based on angle monitoring
CN103913339A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing fault cable, load and generalized displacement based on angle monitoring

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006136884A1 (en) * 2005-04-04 2006-12-28 Czaloun Hans Guenter Device for measuring indirect cable tensions
JP2007286026A (en) * 2006-03-24 2007-11-01 Tokyo Seiko Co Ltd Method and apparatus for evaluating corrosion degree of area near hanger rope terminal part of hanging bridge or the like
US7992449B1 (en) * 2008-02-15 2011-08-09 Mahmoud Khaled M Method for assessment of cable strength and residual life
CN102288423A (en) * 2011-05-13 2011-12-21 东南大学 Progressive method for recognizing damaged cable and angular displacement of support based on angle monitoring
CN102288433A (en) * 2011-05-13 2011-12-21 东南大学 Progressive method for distinguishing slack cable based on angle monitoring in angular displacement of supporting seat
CN102297771A (en) * 2011-05-13 2011-12-28 东南大学 Progressive health monitoring method for cable system based on angle monitoring during supporting-base angular displacement

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006136884A1 (en) * 2005-04-04 2006-12-28 Czaloun Hans Guenter Device for measuring indirect cable tensions
JP2007286026A (en) * 2006-03-24 2007-11-01 Tokyo Seiko Co Ltd Method and apparatus for evaluating corrosion degree of area near hanger rope terminal part of hanging bridge or the like
US7992449B1 (en) * 2008-02-15 2011-08-09 Mahmoud Khaled M Method for assessment of cable strength and residual life
CN102288423A (en) * 2011-05-13 2011-12-21 东南大学 Progressive method for recognizing damaged cable and angular displacement of support based on angle monitoring
CN102288433A (en) * 2011-05-13 2011-12-21 东南大学 Progressive method for distinguishing slack cable based on angle monitoring in angular displacement of supporting seat
CN102297771A (en) * 2011-05-13 2011-12-28 东南大学 Progressive health monitoring method for cable system based on angle monitoring during supporting-base angular displacement

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
侯俊明 等: ""斜拉索索力的温度敏感性"", 《长安大学学报(自然科学版)》 *
冯伟文: ""多塔矮塔斜拉桥成桥状态温度影响分析"", 《广东土木与建筑》 *
吴志友: "基于索力监测的大跨索支承桥梁受损索识别的研究", 《东南大学硕士论文(万方数据库)》 *
缪长青 等: "大跨缆索支承型桥梁健康监测与评估系统的设计研究", 《特种结构》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103604632A (en) * 2013-12-09 2014-02-26 东南大学 Problematic cable/concentrated load progressive recognition method on basis of generalized displacement angle monitoring
CN103616248A (en) * 2013-12-09 2014-03-05 东南大学 Progressive identifying method for problem cable concentrated loads and angular displacement based on angle monitoring
CN103630412A (en) * 2013-12-09 2014-03-12 东南大学 Method for progressively identifying concentrated loads of problematic cables by means of monitoring angles during angular displacement
CN103852300A (en) * 2014-03-10 2014-06-11 东南大学 Progressive recognition method for damaged cable loads through generalized displacement and angle monitoring
CN103852337A (en) * 2014-03-10 2014-06-11 东南大学 Strain monitoring error cable load angular displacement progressive identification method
CN103852332A (en) * 2014-03-10 2014-06-11 东南大学 Damaged cable load linear displacement progressive recognition method based on space coordinate monitoring
CN103852330A (en) * 2014-03-10 2014-06-11 东南大学 Method for recognizing damaged cables, loads and linear displacement based on cable force monitoring
CN103852299A (en) * 2014-03-10 2014-06-11 东南大学 Linear displacement space coordinate monitoring damaged cable load identification method
CN103884527A (en) * 2014-03-10 2014-06-25 东南大学 Method for recognizing load angular displacement and issue cables through angle monitoring
CN103913337A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement cable force monitoring defective cable load progressive identification method
CN103913322A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing fault cable and load based on angular displacement angle monitoring
CN103913329A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement hybrid monitoring damaged cable load identification method
CN103913332A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing damaged cable and load based on angular displacement strain monitoring
CN103913320A (en) * 2014-03-10 2014-07-09 东南大学 Method for progressively recognizing damaged cable, load and angular displacement based on angle monitoring
CN103913339A (en) * 2014-03-10 2014-07-09 东南大学 Method for recognizing fault cable, load and generalized displacement based on angle monitoring

Also Published As

Publication number Publication date
CN102706659B (en) 2015-07-08

Similar Documents

Publication Publication Date Title
CN103913328A (en) Generalized displacement hybrid monitoring damaged cable load progressive identification method
CN102706659B (en) Defective cable and support angular displacement progressive identification method based on angular monitoring of temperature change
CN103913342A (en) Method for progressively recognizing fault cable, load and generalized displacement based on angle monitoring
CN102735459B (en) The problem cable generalized displacement of support progressive-type recognition method of temperature variation cable force monitoring
CN102706627B (en) The damaged cable of temperature variation hybrid monitoring and support angular displacement identification method
CN103604549A (en) Linear displacement progressive type identification method for strain monitoring of damaged cable and intensive load
CN103616249A (en) Strain-monitoring damaged cable centralized load generalized displacement recognition method
CN103604642A (en) Problematic cable/concentrated load progressive recognition method on basis of generalized displacement cable power monitoring
CN103616222A (en) Hybrid-monitoring linear displacement progressive identification method for defective cable and concentrated loads
CN102706605B (en) The problem cable generalized displacement of support progressive-type recognition method of temperature variation strain monitoring
CN103604634A (en) Linear displacement progressive type identification method for space coordinate monitoring of damaged cable and intensive load
CN103604623A (en) Progressive type damaged rope and concentrated load identification method based on angular displacement strain monitoring
CN102706575B (en) Damaged cable and supporting seat translation progressive-type identification method based on space coordinate monitoring at moment of temperature variation
CN103604628A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement angle monitoring
CN103852331A (en) Progressive recognition method for problem cable load generalized displacement based on mixing monitoring
CN103852304A (en) Method for recognizing damaged cables, loads and linear displacement based on cable force monitoring
CN102706595B (en) The damaged cable of angle monitor and support translation progressive identification method during temperature variation
CN103604644A (en) Problematic cable/concentrated load progressive recognition method on basis of linear displacement strain monitoring
CN103852309A (en) Progressive recognition method for problem cable load linear displacement based on mixing monitoring
CN103604649A (en) Problematic cable/concentrated load/linear displacement progressive recognition method on basis of strain monitoring
CN103616245A (en) Progressive concentrated load support angular displacement identification method for damaged cable in process of cable force monitoring
CN103604629A (en) Damaged cable/concentrated load/bracket angular displacement progressive recognition method on basis of strain monitoring
CN103604640A (en) Method for progressive recognition of concentrated loads of damaged cables based on generalized displacement angle monitoring
CN103604630A (en) Problematic cable/concentrated load/linear displacement recognition method on basis of hybrid monitoring
CN103913318A (en) Method for progressively recognizing damaged cable and load based on angular displacement space coordinates monitoring

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150708

Termination date: 20180530