CN102706638A - Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation - Google Patents

Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation Download PDF

Info

Publication number
CN102706638A
CN102706638A CN2012101726524A CN201210172652A CN102706638A CN 102706638 A CN102706638 A CN 102706638A CN 2012101726524 A CN2012101726524 A CN 2012101726524A CN 201210172652 A CN201210172652 A CN 201210172652A CN 102706638 A CN102706638 A CN 102706638A
Authority
CN
China
Prior art keywords
cable structure
temperature
data
vector
cable
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101726524A
Other languages
Chinese (zh)
Other versions
CN102706638B (en
Inventor
韩玉林
万江
韩佳邑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201210172652.4A priority Critical patent/CN102706638B/en
Publication of CN102706638A publication Critical patent/CN102706638A/en
Application granted granted Critical
Publication of CN102706638B publication Critical patent/CN102706638B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

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

Abstract

The invention relates to a damaged cable identification method on the basis of the hybrid monitoring during support settlement and temperature variation. According to the damaged cable identification method, on the basis of the hybrid monitoring, the translation displacement of a support, the temperatures of cable structures and the environment temperature are monitored to determine whether a mechanical calculation reference model of the cable structures needs to be updated so as to obtain a mechanical calculation reference model of the cable structures, in which the translation displacement, the temperatures of the cable structures and the environment temperature are counted; on the basis of the model, a change matrix of a unit damage monitored quantity is obtained by calculation; and according to approximate linear relations between a current value vector of the monitored quantity and a current initial value vector of the monitored quantity, the change matrix of the unit damage monitored quantity, a unit damage scalar and a current nominal damage vector to be obtained of a cable system, a noninferior solution of the current nominal damage vector of the cable system is calculated, so that when the support has displacement and the temperature is changed, the positions and the damage degrees of damaged cables can be accurately determined.

Description

When support settlement and temperature variation are arranged based on the damaged cable recognition methods of hybrid monitoring
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; This method is discerned the supporting system of Cable Structure based on hybrid monitoring, and (refer to all ropeway carrying-ropes, and 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 comprise the rod member that only bears tensile load, censure 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) in damaged cable (truss-frame structure just is meant the impaired rod member that only bears tensile load), belong to the engineering structure health monitoring field.
Background technology
The health status of supporting cable system changes after (for example damaging); 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, therefore can be through the hybrid monitoring of the variation of the characteristic parameter of these dissimilar structures being judged the health status of structure, and this method architectural feature parameter that all are monitored is referred to as " monitored amount "; Because this moment, monitored amount was made up of the dissimilar measurable parameter mixing of structure, this method claims that this is a hybrid monitoring.Monitored amount is except the influence that receives the cable system health status; Also can receive the influence of Cable Structure temperature variation (usually can take place) and Cable Structure bearing translational displacement (sedimentation is the component of translational displacement at gravity direction); Under the Cable Structure temperature changes condition with Cable Structure bearing generation translational displacement; If can realize identification based on monitoring to monitored amount to the supporting rope of unsoundness problem; Safety to Cable Structure has significant values, does not also have a kind of disclosed, effective health monitoring systems and method to solve this problem at present.
When Cable Structure has support settlement and temperature variation; In order reliable monitoring and judgement to be arranged to the health status of the cable system of Cable Structure; Must there be one can rational and effective set up the method for each monitored quantitative changeization, can provides the health evaluating of more believable cable system based on the health monitoring systems of this method foundation with the relation between the health status of all ropes in the cable system.
Summary of the invention
Technical matters: the purpose of this method is when Cable Structure has support settlement and temperature variation; To the health monitoring problem of cable system in the Cable Structure, a kind of health monitor method hybrid monitoring, that can monitor cable system in the Cable Structure rationally and effectively to the multiclass parameter is disclosed.
Technical scheme: this method is made up of three parts.Be respectively the method for setting up required knowledge base of cable system health monitoring systems and parameter, based on knowledge base (containing parameter) and the cable system health status appraisal procedure of the monitored amount of actual measurement, the software and hardware part of health monitoring systems.
The first of this method: the method for setting up required knowledge base of cable system health monitoring systems and parameter.Specific as follows:
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.
2. set up the initial Mechanics Calculation benchmark model A of Cable Structure o(for example finite element benchmark model) and current initial Mechanics Calculation benchmark model A t oThe method of (for example finite element benchmark model) is set up and A oCorresponding monitored amount initial value vector C oMethod, set up and A t oThe current initial value vector of corresponding monitored amount C t oMethod.A in the method oAnd C oBe constant.A t oAnd C t oBring in constant renewal in.Set up A oAnd C o, set up and upgrade A t oAnd C t oMethod following.
If total N root supporting rope, the coding rule of at first definite supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule.
Monitored multiclass parameter can comprise: Suo Li, strain, angle and volume coordinate, and narration is as follows respectively:
If total N root rope in the cable system, the monitored rope force data of structure is by M on the structure 1The M of individual appointment rope 1Individual rope force data is described, and the variation of structure Suo Li is exactly all variations of specifying the Suo Li of rope.Each total M 1Individual cable force measurement value or calculated value characterize the rope force information of structure.M 1Be one and be not less than 0 integer.
The monitored strain data of structure can be by K on the structure 2L individual specified point, that reach each specified point 2The strain of individual assigned direction is described, and the variation of structural strain data is exactly K 2The variation of the tested strain of all of individual specified point.Each total M 2(M 2=K 2* L 2) individual strain measurement value or calculated value characterize structural strain.M 2Be one and be not less than 0 integer.
The monitored angle-data of structure is by K on the structure 3L individual specified point, that cross each specified point 3H individual appointment straight line, each appointment straight line 3Individual angle coordinate component is described, and the variation of structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M 3(M 3=K 3* L 3* H 3) individual angle coordinate component measurement value or calculated value characterize the angle information of structure.M 3Be one and be not less than 0 integer.
The monitored shape data of structure is by K on the structure 4L individual specified point, that reach each specified point 4The volume coordinate of individual assigned direction is described, and the variation of planform data is exactly K 4The variation of all coordinate components of individual specified point.Each total M 4(M 4=K 4* L 4) individual measurement of coordinates value or calculated value characterize planform.M 4Be one and be not less than 0 integer.
Comprehensive above-mentioned monitored amount, total has M (M=M 1+ M 2+ M 3+ M 4) individual monitored amount, definition parameter K (K=M 1+ K 2+ K 3+ K 4), K and M must not be less than the quantity N of rope.Because M monitored amount is dissimilar, so this method is called hybrid monitoring.
For simplicity, in the method " all monitored parameters of structure " are abbreviated as " monitored amount ".
Set up initial Mechanics Calculation benchmark model A oThe time; In Cable Structure completion; Perhaps before setting up health monitoring (damaged cable identification) 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 at the synchronization in the moment that obtains Cable Structure steady temperature data, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.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 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 Non-Destructive Testing data of supporting 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, 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.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 of the health status of rope and set up cable system initial damage vector d oD is used in (shown in (1)) oThe expression Cable Structure is (with initial Mechanics Calculation benchmark model A oThe initial health of cable system expression).If when not having Non-Destructive Testing data and other of rope can express the data of health status of supporting rope, perhaps can think when the structure original state is the not damaged state vectorial d oEach element numerical value 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 oCorresponding to A oCable Structure support coordinate data form initial Cable Structure support coordinate vector U o
d o=[d o1?d o2···d oj···d oN] T (1)
D in the formula (1) Oj(j=1,2,3 ...., N) the initial Mechanics Calculation benchmark model A of expression oIn the initial damage value of j root rope of cable system, d OjBeing to represent j root 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 j root rope forfeiture corresponding proportion in the time of between 0 and 100%, and T representes the transposition (back with) 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.
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 oThe health status of middle supporting rope is with cable system initial damage vector d 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 oWe can say T o, U oAnd d oBe A oParameter, C oBy A oMechanics Calculation result form.
Set up and upgrade current initial Mechanics Calculation benchmark model A t oMethod be: (just set up for the first time A at initial time t oThe time), A t oJust equal A o, A t oCorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T t o", at initial time, T t oJust equal T o, vector T t oDefinition mode and vector T oDefinition mode identical.Current initial Mechanics Calculation benchmark model A corresponding to Cable Structure t oCable Structure support coordinate data form current initial Cable Structure support coordinate vector U t o, just set up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure at initial time t oThe time, U t oJust equal U oA t oThe initial health and the A of supporting rope oThe health status of supporting rope identical, also use cable system initial damage vector d oExpression, A in the cyclic process of back t oThe initial health of supporting rope use cable system initial damage vector d all the time oExpression; Cable Structure is in A t oDuring state, this method is with the current initial value vector of monitored amount C t oThe concrete numerical value of representing all monitored amounts, C t oElement and C oElement corresponding one by one, represent that respectively all monitored amounts are in A in Cable Structure t oAnd A oConcrete numerical value during two states.At initial time, C t oJust equal C o, T t o, U t oAnd d oBe A t oParameter, C t oBy A t oMechanics Calculation result form; 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 " and (be called " current cable structure steady temperature data vector T t", vector T tDefinition mode and vector T oDefinition mode identical); Obtaining vector T tThe time, actual measurement obtains Cable Structure support coordinate current data, and all Cable Structure support coordinate current datas are formed current cable structure actual measurement support coordinate vector U tIf T tEqual T t oAnd U tEqual U t o, then need be to A t oUpgrade, otherwise need be to A t o, U t oAnd T t oUpgrade, update method is: the first step is calculated U tWith U oPoor, U tWith U oDifference be exactly the current bearing translational displacement of Cable Structure bearing about initial position; V representes the bearing translational displacement with bearing translational displacement vector; Between element among the bearing translational displacement vector V and the bearing translational displacement component is one-to-one relationship; The numerical value of an element is corresponding to the displacement of an assigned direction of an appointment bearing among the bearing translational displacement vector V, and wherein the bearing translational displacement is exactly the support settlement amount at the component of gravity direction; Second step was calculated T tWith T oPoor, T tWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T tWith T oThe difference represent that with steady temperature change vector S S equals T tDeduct T o, S representes the variation of Cable Structure steady temperature data; The 3rd step is earlier to A oIn the Cable Structure bearing apply the constraint of current bearing translational displacement, the numerical value of current bearing translational displacement constraint is just taken from the numerical value of corresponding element among the bearing translational displacement vector V, again 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 Cable Structure bearing apply bearing translational displacement constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading t o, upgrade A t oThe time, U t oAll elements numerical value is also used U tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U t o, T t oAll elements numerical value is also used T tCorresponding replacement of all elements numerical value, promptly upgraded T t o, so just obtained correctly corresponding to A t oT t oUpgrade C t oMethod be: when upgrading A t oAfter, obtain A through Mechanics Calculation t oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C t o
The currency of all monitored amounts is formed the current numerical value vector of monitored amount C (formula (3) is seen in definition) in the Cable Structure.
C=[C 1?C 2···C j···C M] T (3)
C in the formula (3) j(j=1,2,3 ...., M) be the currency of j monitored amount in the Cable Structure, this component C jAccording to coding rule and C OjCorresponding to same " monitored amount ".Obtain current cable structure steady temperature data vector T in actual measurement tSynchronization, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, forms the current numerical value vector of monitored amount C.
3. set up and upgrade the method for Cable Structure unit damage monitored quantitative change matrix Δ C.
Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current initial Mechanics Calculation benchmark model A t oWith the current initial value vector of monitored amount C t oThe time, upgrade Cable Structure unit damage monitored quantitative change matrix Δ C.Concrete grammar is following:
Current initial Mechanics Calculation benchmark model A in Cable Structure t oThe basis on carry out several times and calculate, equal the quantity of all supporting ropes on the calculation times numerical value.Calculating each time in the hypothesis cable system has only a supporting rope (to use vectorial d at initial damage oCorresponding element represent) the basis on increase unit damage D again u(for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage); The rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating; Calculate the current calculated value that all utilizes mechanics method (for example finite element 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, and (when hypothesis i root rope had unit damage, available formula (4) was represented monitored amount calculation current vector ); Calculate monitored amount calculation current vector each time and deduct the current initial value vector of monitored amount C t o, the gained vector is exactly that the monitored quantitative changeization vector of (is mark with the position of supporting rope that unit damage is arranged or numbering etc.) (when i root rope has unit damage, is used δ C under this condition iRepresent monitored quantitative changeization vector, formula (5) is seen in definition), each element representation of monitored quantitative changeization vector is owing to suppose the change amount of the pairing monitored amount of this element that the unit damage of the Na Gensuo that unit damage is arranged causes when calculating; There is N root rope that N monitored quantitative changeization vector just arranged; Because N monitored amount arranged; So each monitored quantitative change vector has N element; Being made up of successively this N monitored quantitative change vector has the unit damage monitored of M * N element quantitative change matrix Δ C, and the definition of Δ C is suc as formula shown in (6).
C t i = C t 1 i C t 2 i · · · C tj i · · · C tM i T - - - ( 4 )
Element
Figure BDA00001694030700123
(i=1 in the formula (4); 2; 3 ...., N; J=1,2,3 ...., when M) expression has unit damage owing to i root rope, according to the current calculated amount of the pairing j of coding rule monitored amount.
δC i = C t i - C o t - - - ( 5 )
ΔC = Δ C 1,1 Δ C 1,2 · Δ C 1 , i · Δ C 1 , N Δ C 2,1 Δ C 2,2 · Δ C 2 , i · Δ C 2 , N · · · · · · Δ C j , i Δ C j , 2 · Δ C j , i · Δ C j , N · · · · · · Δ C M , 1 Δ C M , 2 . Δ C M , i . Δ C M , N - - - ( 6 )
Formula (Δ C in 60 J, i(i=1,2,3 ...., N; J=1,2,3 ...., M) expression only since i root rope have that unit damage causes, according to the variation (algebraic value) of the current numerical value of calculating of the individual monitored amount of the pairing j of coding rule.Monitored quantitative changeization vector δ C iBe actually the row among the matrix Δ C, that is to say that formula (6) also can write an accepted way of doing sth (7).
ΔC=[δC 1?δC 2···δC i···δC N] (7)
4. the current numerical value vector of monitored amount C (calculating or actual measurement) is with the current initial value vector of monitored amount C t o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D UAnd the linear approximate relationship between the vectorial d of the current name damage of cable system, shown in (8) or formula (9).The definition of the vectorial d of the current name damage of cable system is referring to formula (10).
C = C o t + 1 D u ΔC · d - - - ( 8 )
C - C o t = 1 D u ΔC · d - - - ( 9 )
d=[d 1?d 2···d i···d N] T (10)
D in the formula (10) i(i=1,2,3 ...., N) be the current name damage of i root rope (or pull bar) in the cable system.
Represented that rope thoroughly lost load-bearing capacity at 100% o'clock if establish rope damage and be; When actual damage is not too big, (for example be not more than 30% damage) so; 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 (8) or formula (9) is less with the error of actual conditions.Error with linear relationship shown in the linear relationship error vector e expression (8) of formula (11) definition or the formula (9).
e = abs ( 1 D u ΔC · d - C + C o t ) - - - ( 11 )
Abs () is the function that takes absolute value in the formula (11), and each element of the vector of trying to achieve in the bracket is taken absolute value.
The second portion of this method: based on the cable system health status appraisal procedure of knowledge base (containing parameter) and the monitored amount of actual measurement.
Because there are certain error in formula (8) or the represented linear relationship of formula (9), therefore can not be simply come directly to find the solution to obtain the current name of cable system with the current numerical value vector of the monitored amount of actual measurement C and damage vectorial d according to formula (8) or formula (9).If done like this, the current name of the cable system that obtains is damaged the element among the vectorial d even bigger negative value can be occurred, just negative damage, and this obviously is irrational.Therefore acceptable the separating that obtains the vectorial d of the current name damage of cable system (promptly has reasonable error; But can from cable system, confirm the position and the degree of injury thereof of damaged cable more accurately) become a rational solution, available formula (12) is expressed this method.
abs ( 1 D u ΔC · d - C + C o t ) ≤ g - - - ( 12 )
Abs () is the function that takes absolute value in the formula (12), and vectorial g describes the reasonable deviation that departs from ideal linearity relation (formula (8) or formula (9)), is defined by formula (13).
g=[g 1?g 2···g j···g M] T (13)
G in the formula (13) j(j=1,2,3 ...., M) maximum allowable offset of the ideal linearity relation that departs from shown in formula (8) or the formula (9) has been described.Vector g can be selected according to the error vector e tentative calculation of formula (11) definition.
At the current initial value vector of monitored amount C t o, unit damage monitored quantitative change matrix Δ C, the current numerical value of the monitored amount of actual measurement vector C and unit damage D uWhen (setting before calculating Δ C, is scalar) is known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (12), obtain the current name damage of cable system the acceptable of vectorial d and separate.
The current actual damage vector of definition cable system d a(seeing formula (14)), the current actual damage vector of cable system d aElement can calculate according to formula (15), just obtained Suo Dangqian actual damage vector d aThereby, can be by d aConfirm the position and the degree of injury of damaged cable, just realized the health monitoring of cable system, realized damaged cable identification.
d a = d 1 a d 2 a · · · d j a · · · d N a T - - - ( 14 )
D in the formula (14) a j(j=1,2,3 ...., N) the actual damage value of expression j root rope, formula (15), d are seen in its definition a jBeing to represent j root 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 j root rope forfeiture corresponding proportion in the time of between 0 and 100%, vectorial d aCoding rule and the formula (1) of element in vectorial d oThe coding rule of element identical.
d j a = 1 - ( 1 - d oj ) ( 1 - d j ) - - - ( 15 )
D in the formula (15) Oj(i=1,2,3,4, J=1,2,3 ...., N) be vectorial d oJ element, d jBe j the element of vectorial d.Thereby confirm the position and the degree of injury of damaged cable.
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, Cable Structure support coordinate monitoring system), signal picker and computing machine etc.Require monitoring in real time to obtain measured data temperature required and the Cable Structure support coordinate, require simultaneously each monitored amount of monitoring in real time.
Functions such as the monitoring that software should be accomplished is needed in this method, can use computer realization, record, control, storage, calculating, notice, warning.
This method specifically comprises:
A. establish total N root supporting rope, at first confirm the coding rule of supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The supporting rope with monitored Suo Li of appointment when confirming hybrid monitoring is established total N root supporting rope in the cable system, and the monitored rope force data of Cable Structure is by M on the Cable Structure 1The M of individual appointment supporting rope 1Individual rope force data is described, and the variation of Cable Structure Suo Li is exactly all variations of specifying the Suo Li of supporting rope; Each total M 1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M 1Be one and be not less than 0 integer; The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of Cable Structure can be by K on the Cable Structure 2L individual specified point, that reach each specified point 2The strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K 2The variation of the tested strain of all of individual specified point; Each total M 2Individual strain measurement value or calculated value characterize Cable Structure strain, M 2Be K 2And L 2Long-pending; M 2Be to be not less than 0 integer; The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of Cable Structure is by K on the Cable Structure 3L individual specified point, that cross each specified point 3H individual appointment straight line, each appointment straight line 3Individual angle coordinate component is described, and the variation of Cable Structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total M 3Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M 3Be K 3, L 3And H 3Long-pending; M 3Be one and be not less than 0 integer; When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of Cable Structure is by K on the Cable Structure 4L individual specified point, that reach each specified point 4The volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K 4The variation of all coordinate components of individual specified point; Each total M 4Individual measurement of coordinates value or calculated value characterize Cable Structure shape, M 4Be K 4And L 4Long-pending; M 4Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M 1, M 2, M 3And M 4Sum, definition parameter K, K is M 1, K 2, K 3And K 4Sum, K and M must not be less than the quantity N of rope; For simplicity, in the method " all monitored parameters of Cable Structure " are abbreviated as " monitored amount "; 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 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 of the health status of expressing the supporting rope, 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; The initial value of all monitored amounts is formed monitored amount initial value vector C oUtilize the Non-Destructive Testing data of the health status that can express the supporting rope to set up cable system initial damage vector d o, cable system initial damage vector d oElement number equal N, d oElement be one-to-one relationship with the supporting rope, cable system initial damage vector d oElement numerical value represent the degree of injury of corresponding supporting rope, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the pairing supporting rope of this element is intact, not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%; Represent that then this supporting rope lost the load-bearing capacity of corresponding proportion; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, perhaps think when the Cable Structure original state is the not damaged state vectorial d oEach element numerical value get 0; Initial Cable Structure support coordinate data are formed initial Cable Structure support coordinate vector U o
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, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure support coordinate vector U o, 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 oSupporting rope health status with cable system initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression; Set up the current initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " for the first time t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure t oWith the current initial value vector of monitored amount C t oThe time, the current initial Mechanics Calculation benchmark model A of Cable Structure t oJust equal the initial Mechanics Calculation benchmark model A of Cable Structure o, the current initial value vector of monitored amount C t oJust equal monitored amount initial value vector C oA t oCorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up the current initial Mechanics Calculation benchmark model A of Cable Structure the first time t oThe time, T t oJust equal T oCurrent initial Mechanics Calculation benchmark model A corresponding to Cable Structure t oCable Structure support coordinate data form current initial Cable Structure support coordinate vector U t o, set up the current initial Mechanics Calculation benchmark model A of Cable Structure the first time t oThe time, U t oJust equal U oA t oThe initial health and the A of supporting rope oThe health status of supporting rope identical, also use cable system initial damage vector d oExpression, A in the cyclic process of back t oThe initial health of supporting rope use cable system initial damage vector d all the time oExpression; T o, U 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 T t o, U t oAnd d oBe A t oParameter, C t oBy A t oMechanics Calculation result form; A in the method o, U o, C o, d oAnd T oBe constant;
E. go on foot the m circulation in step from getting into here by e; In Cable Structure military service process; Constantly constantly survey and calculate the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method "; The current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tDefinition mode and vector T oDefinition mode identical; Obtain current cable structure steady temperature data vector T in actual measurement tSynchronization, actual measurement obtains Cable Structure support coordinate current data, all Cable Structure support coordinate current datas are formed current cable structures actual measurement support coordinate vector U t
F. according to current cable structure actual measurement support coordinate vector U tWith current cable structure steady temperature data vector T t, upgrade current initial Mechanics Calculation benchmark model A according to step f1 to f3 t o, current initial Cable Structure support coordinate vector U t o, the current initial value of monitored amount vector C t oWith current initial Cable Structure steady temperature data vector T t o
F1 compares U respectively tWith U t o, T tWith T t oIf, U tEqual U t oAnd T tEqual T t o, A then t o, U t o, C t oAnd T t oRemain unchanged; Otherwise need follow these steps to A t o, U t oAnd T t oUpgrade;
F2. calculate U tWith U oPoor, U tWith U oDifference be exactly the current bearing translational displacement of Cable Structure bearing about initial position, represent the bearing translational displacement with bearing translational displacement vector V, V equals U tDeduct U oBetween element among the bearing translational displacement vector V and the bearing translational displacement component is one-to-one relationship; The numerical value of an element is corresponding to the translational displacement of an assigned direction of an appointment bearing among the bearing translational displacement vector V, and wherein the bearing translational displacement is exactly the support settlement amount at the component of gravity direction; Calculate T tWith T oPoor, T tWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T tWith T oThe difference represent that with steady temperature change vector S S equals T tDeduct T o, S representes the variation of Cable Structure steady temperature data;
F3. earlier to A oIn the Cable Structure bearing apply the constraint of current bearing translational displacement, the numerical value of current bearing translational displacement constraint is just taken from the numerical value of corresponding element among the bearing translational displacement vector V, again 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 Cable Structure bearing apply bearing translational displacement constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading t o, upgrade A t oThe time, U t oAll elements numerical value is also used U tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U t o, T t oAll elements numerical value is also used T tCorresponding replacement of all elements numerical value, promptly upgraded T t o, so just obtained correctly corresponding to A t oT t oAnd U t oUpgrade C t oMethod be: when upgrading A t oAfter, obtain A through Mechanics Calculation t oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C t oA t oThe initial health of supporting rope use cable system initial damage vector d all the time oExpression;
G. at current initial Mechanics Calculation benchmark model A t oThe basis on carry out the several times Mechanics Calculation according to step g 1 to g4, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D u
G1. Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current initial Mechanics Calculation benchmark model A t o, current initial Cable Structure support coordinate vector U t o, the current initial value of monitored amount vector C t oWith current initial Cable Structure steady temperature data vector T t oAfterwards, must then upgrade Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D u
G2. at the current initial Mechanics Calculation benchmark model A of Cable Structure t oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have N root supporting rope that N calculatings is just arranged, calculating each time supposes to have only a supporting rope that unit damage scalar D is arranged in the cable system uThe rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating; Calculate the current calculated value of all monitored amounts in the 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, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C oThe element coding rule identical;
G3. the monitored amount calculation current vector that calculates each time deducts the current initial value vector of monitored amount C t oObtain a monitored quantitative changeization vector; There is N root supporting rope that N monitored quantitative changeization vector just arranged;
G4. form the Cable Structure unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N monitored quantitative change vector; Each row of Cable Structure unit damage monitored quantitative change matrix Δ C are corresponding to a monitored quantitative changeization vector;
H. obtain current cable structure steady temperature data vector T in actual measurement tThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T tThe current measured value of all monitored amounts of Cable Structure of synchronization in the moment, form the current numerical value vector of monitored amount C; Current numerical value vector C of monitored amount and the current initial value vector of monitored amount C t oWith monitored amount initial value vector C oDefinition mode identical, the same monitored amount of the element representation of the identical numbering of three vectors is at difference concrete numerical value constantly;
I. define the vectorial d of the current name damage of cable system; The element number of the vectorial d of the current name damage of cable system equals to support the quantity of rope; Between the element of the vectorial d of the current name damage of cable system and the supporting rope is one-to-one relationship, and the element numerical value of the vectorial d of the current name damage of cable system is represented the nominal degree of injury or the nominal health status of corresponding supporting rope; The coding rule and the vectorial d of the element of vector d oThe coding rule of element identical;
J. the current numerical value vector of the monitored amount of foundation C is with the current initial value vector of monitored amount C t o, Cable Structure unit damage monitored quantitative change matrix Δ C, unit damage scalar D uAnd the linear approximate relationship that exists between the vectorial d of the current name damage of cable system to be asked, this linear approximate relationship can be expressed as formula 1, and other amount in the formula 1 except that d is known, finds the solution formula 1 and just can calculate the vectorial d of the current name damage of cable system;
C = C o t + 1 D u Δ C · d Formula 1
K. define the current actual damage vector of cable system d a, the current actual damage vector of cable system d aElement number equal to support the quantity of rope, the current actual damage vector of cable system d aElement and supporting be one-to-one relationship between the rope, the current actual damage vector of cable system d aElement numerical value represent the actual damage degree or the actual health status of corresponding supporting rope; Vector d aThe coding rule and the vectorial d of element oThe coding rule of element identical;
1. the current actual damage vector of the cable system that utilizes formula 2 to express d aJ element d a jWith cable system initial damage vector d oJ element d OjJ element d with the vectorial d of the current name damage of cable system jBetween relation, calculate the current actual damage of cable system vector d aAll elements;
d j a = 1 - ( 1 - d Oj ) ( 1 - d j ) Formula 2
J=1 in the formula 2,2,3 ...., N, d a jBeing to represent j 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 in the time of between 0 and 100% that j root supporting rope loses the load-bearing capacity of corresponding proportion; The current actual damage vector of cable system d aElement numerical value represent the degree of injury of corresponding supporting rope, so according to the current actual damage vector of cable system d aCan define the impaired and degree of injury of which rope, promptly realize the damaged cable identification or the health monitoring of cable system in the Cable Structure;
M. get back to the e step, beginning goes on foot the m circulation next time in step by e.
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 the rational structure health monitoring based on monitored amount; And the temperature field measurement of Cable Structure with calculate to be very complicated, this method discloses and has comprised and adopt this method to occur under the translational displacement situation of (comprising sedimentation) at the Cable Structure bearing by a kind of simple, economic, feasible, cable structure health monitoring method of structure temperature field computing method efficiently that is suitable for monitoring structural health conditions; Synchronous when impaired at many ropes of Cable Structure; And when the temperature of Cable Structure changed along with the time, monitoring and evaluation identified the health status (position and the degree of injury thereof that comprise all damaged cables) of cable system very exactly, and the disclosed system and method for this method is very useful to effective health monitoring of cable system.
Embodiment
When bearing translational displacement (support settlement is the component of bearing translational displacement) and temperature variation are arranged; To the health monitoring of the cable system of Cable Structure, this method discloses a kind of system and method for the health status of each root rope in the cable system of can monitoring rationally and effectively in the identification Cable Structure.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 for monitoring the health status of the cable system of Cable Structure.During practical implementation, the following step is a kind of in the various steps that can take.
The first step: 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 SmaxAll be 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
If total N root supporting rope, the coding rule of at first definite rope, with rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule.The supporting rope with monitored Suo Li of appointment when confirming hybrid monitoring is established total N root rope in the cable system, and the monitored rope force data of Cable Structure is by M on the structure 1The M of individual appointment supporting rope 1Individual rope force data is described, and the variation of Cable Structure Suo Li is exactly all variations of specifying the Suo Li of supporting rope.Each total M 1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure.M 1Be one and be not less than 0 integer.When reality was selected the rope of monitored Suo Li, the rope that can select those Suo Li to be easy to measure was monitored rope.
The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of Cable Structure can be by K on the Cable Structure 2L individual specified point, that reach each specified point 2The strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K 2The variation of the tested strain of all of individual specified point.Each total M 2Individual strain measurement value or calculated value characterize structural strain, M 2Be K 2And L 2Long-pending.M 2Be one and be not less than 0 integer.A near point each fixed endpoint that can be exactly each root rope (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge) with the measured point of monitored strain; This point generally should not be a stress concentration point; Avoiding occurring excessive strain measurement value, near the fixed endpoint of the rope of the monitored Suo Li of appointment or it when these points generally should all be not hybrid monitoring yet.
The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of Cable Structure is by K on the Cable Structure 3L individual specified point, that cross each specified point 3H individual appointment straight line, each appointment straight line 3Individual angle coordinate component is described, and the variation of Cable Structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments.Each total M 3Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M 3Be K 3, L 3And H 3Long-pending.M 3Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (drag-line that for example is cable-stayed bridge is at the stiff end on the bridge floor) or a near point it of each root rope, and the point of monitored angle-data generally should all not be chosen as " fixed endpoint of the rope of the monitored Suo Li of appointment or near the point it in the hybrid monitoring " and " point of the monitored strain of appointment or near the point it in the hybrid monitoring "; Can only measure an angle coordinate of specifying straight line at each specified point, for example measure body structure surface normal or the tangent line of the specified point angle coordinate with respect to the acceleration of gravity direction, in fact be exactly measurement of dip angle here.
When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of Cable Structure is by K on the structure 4L individual specified point, that reach each specified point 4The volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K 4The variation of all coordinate components of individual specified point.Each total M 4Individual measurement of coordinates value or calculated value characterize Cable Structure shape, M 4Be K 4And L 4Long-pending.M 4Be one and be not less than 0 integer.Each specified point can be exactly the fixed endpoint (for example being the stiff end of drag-line on bridge of cable-stayed bridge) of each root rope; Here selected monitored point should all not selected " fixed endpoint of the rope of the monitored Suo Li of appointment or near the point it in the hybrid monitoring ", " point of the monitored strain of appointment or near the point it in the hybrid monitoring " and " point of the monitored angle-data of appointment or near the point it in the hybrid monitoring " for use.
Comprehensive above-mentioned monitored amount, whole Cable Structure be total M monitored amount with regard to hybrid monitoring, and M is M 1, M 2, M 3And M 4Sum, definition parameter K, K is M 1, K 2, K 3And K 4Sum, K and M must not be less than the quantity N of rope.Because M monitored amount is dissimilar, so the present invention is called hybrid monitoring.For simplicity, in the present invention that this step is listed " all monitored parameters of Cable Structure during hybrid monitoring " abbreviate " monitored amount " as.
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; 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 modal data, structural strain data, structure angle measurement data, structure space measurement of coordinates data.Initial Cable Structure support coordinate data are formed initial Cable Structure support coordinate vector U oThe 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.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 of the health status of rope and set up cable system initial damage vector d oIf when not having Non-Destructive Testing data and other of rope can express the data of health status of supporting rope, perhaps can think when the structure original state is the not damaged state vectorial d oEach element numerical value get 0.Cable system initial damage vector d oElement number equal N, d oElement be one-to-one relationship with the supporting rope, cable system initial damage vector d oElement numerical value be not less than 0, be not more than 100%, d oElement numerical value represent the degree of injury of corresponding supporting rope, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the pairing supporting rope of this element is intact, no problem, if its numerical value is 100%; Represent that then the pairing supporting rope of this element has completely lost load-bearing capacity; If its numerical value is between 0 and 100%, represent that then this supporting rope lost the load-bearing capacity of corresponding proportion, if when not supporting Non-Destructive Testing data and other of rope and can express the data of the health status that support rope; Perhaps think when the structure original state is the not damaged state vectorial d oEach element numerical value get 0; If d oThe numerical value of a certain element be not 0, then represent the degree of injury of the pairing supporting rope of this element.Utilize temperature variant physics and the mechanical property parameters of Non-Destructive Testing data, the employed various materials of Cable Structure of measured data, the supporting rope of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, initial Cable Structure support coordinate vector U oWith initial Cable Structure steady temperature data vector T 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 oThe health status of middle supporting rope is with cable system initial damage vector d 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 oSupporting rope health status with cable system initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression.Corresponding to A oCable Structure support coordinate data with initial Cable Structure support coordinate vector U oExpression; T o, U oAnd d oBe A oParameter, C oBy A oMechanics Calculation result form.
The 3rd step: set up current initial Mechanics Calculation benchmark model A for the first time t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o", concrete grammar is: at initial time, promptly set up for the first time current initial Mechanics Calculation benchmark model A t oWith the current initial value vector of monitored amount C t oThe time, A t oJust equal A o, C t oJust equal C o, A t oCorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T t o", (just set up for the first time A at initial time t oThe time), T t oJust equal T o, vector T t oDefinition mode and vector T oDefinition mode identical.Current initial Mechanics Calculation benchmark model A corresponding to Cable Structure t oCable Structure support coordinate data form current initial Cable Structure support coordinate vector U t oSet up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure t oThe time, U t oJust equal U oA t oThe health status and the A of supporting rope oHealth status (the cable system initial damage vector d of supporting rope oExpression) identical, A in cyclic process t oThe health status of supporting rope use cable system initial damage vector d all the time oExpression.T t o, U t oAnd d oBe A t oParameter, C t oBy A t oMechanics Calculation result form.
The 4th step: 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 " and (be called " current cable structure steady temperature data vector T t", vector T tDefinition mode and vector T oDefinition mode identical).Obtain current cable structure steady temperature data vector T in actual measurement tThe time, just obtaining current cable structure steady temperature data vector T tThe synchronization in the moment, actual measurement obtains the current measured value of all monitored amounts of Cable Structure, forms " the current numerical value vector of monitored amount C ".Obtain current cable structure steady temperature data vector T in actual measurement tThe time, actual measurement obtains Cable Structure support coordinate current data, and all data are formed current cable structure actual measurement support coordinate vector U t
The 5th step: according to current cable structure actual measurement support coordinate vector U tWith current cable structure steady temperature data vector T t, upgrade current initial Mechanics Calculation benchmark model A where necessary t o, current initial Cable Structure support coordinate vector U t o, the current initial value of monitored amount vector C t oWith current initial Cable Structure steady temperature data vector T t oObtain current cable structure actual measurement support coordinate vector U in the actual measurement of the 4th step tWith current cable structure steady temperature data vector T tAfter, compare U respectively tAnd U t o, T tAnd T t oIf, U tEqual U t oAnd T tEqual T t o, then need be to A t o, U t oAnd T t oUpgrade, otherwise need be to A t o, U t oAnd T t oUpgrade, update method was undertaken by following a step to the c step:
The a step is calculated U tWith U oPoor, U tWith U oDifference be exactly the current bearing translational displacement of Cable Structure bearing about initial position, represent the bearing translational displacement with bearing translational displacement vector V, V equals U tDeduct U oBetween element among the bearing translational displacement vector V and the bearing translational displacement component is one-to-one relationship; The numerical value of an element is corresponding to the displacement of an assigned direction of an appointment bearing among the bearing translational displacement vector V, and wherein the bearing translational displacement is exactly the support settlement amount at the component of gravity direction.
The b step is calculated T tWith T oPoor, T tWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T tWith T oThe difference represent that with steady temperature change vector S S equals T tDeduct T o, S representes the variation of Cable Structure steady temperature data.
The c step is earlier to A oIn the Cable Structure bearing apply the constraint of current bearing translational displacement, the numerical value of current bearing translational displacement constraint is just taken from the numerical value of corresponding element among the bearing translational displacement vector V, again 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 Cable Structure bearing apply bearing translational displacement constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading t o, upgrade A t oThe time, U t oAll elements numerical value is also used U tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U t o, T t oAll elements numerical value is also used T tCorresponding replacement of all elements numerical value, promptly upgraded T t o, so just obtained correctly corresponding to A t oT t oAnd U t oUpgrade C t oMethod be: when upgrading A t oAfter, obtain A through Mechanics Calculation t oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C t o
The 6th step: at current initial Mechanics Calculation benchmark model A t oThe basis on carry out the several times Mechanics Calculation, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D uConcrete grammar is: Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current initial Mechanics Calculation benchmark model A t oWith current cable structural bearings angular coordinate vector U t oThe time, must upgrade Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D simultaneously uCurrent initial Mechanics Calculation benchmark model A in Cable Structure t oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have N root rope that N calculating is just arranged, calculate each time in the hypothesis cable system and have only a rope that unit damage D is arranged u(for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage); The rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating; Calculate the current calculated value of all monitored amounts in the 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 C; Calculate monitored amount calculation current vector C each time and deduct the current initial value vector of monitored amount C t oObtain a monitored quantitative changeization vector; There is N root rope that N monitored quantitative changeization vector just arranged; Form the unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N monitored quantitative change vector; Each row of unit damage monitored quantitative change matrix are corresponding to a monitored quantitative changeization vector.
The 7th step: set up linear relationship error vector e and vectorial g.Utilize data (the current initial value vector of the monitored amount C of front t o, unit damage monitored quantitative change matrix Δ C), when the 6th step calculated each time, promptly in calculating the hypothesis cable system each time, have only a rope that unit damage D is arranged uThe rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating; Calculate each time and all utilize mechanics method (for example adopting finite element method) to calculate in the Cable Structure the current numerical value of all monitored amounts in the cable system, calculate each time when forming a monitored amount calculation current vector C, calculate each time and form a vectorial d of damage; This is walked out of existing damage vector d and only uses in this step, and this damages in all elements of vectorial d has only the numerical value of an element to get D u, the numerical value of other element gets 0, damages that numerical value is D among the vectorial d uThe unit damage degree D of element unique damaged cable when calculating corresponding to this time uWith C, C t o, Δ C, D u, d brings formula (12) into, obtains a linear relationship error vector e, calculates a linear relationship error vector e each time; Have N root rope that N calculating is just arranged, N linear relationship error vector e just arranged, with obtaining a vector after this N the linear relationship error vector e addition, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e.Vector g equals final error vector e.
The 8th step: the hardware components of pass line structural healthy monitoring system.Hardware components comprises at least: monitored amount monitoring system (for example contains the measurement of angle subsystem; The cable force measurement subsystem; The strain measurement subsystem; Volume coordinate is measured subsystem; Signal conditioner etc.); Cable Structure support coordinate monitoring system (for example measuring) with total powerstation; The Cable Structure temperature monitoring system (contains temperature sensor; Signal conditioner etc.) and Cable Structure ambient temperature measurement system (contain temperature sensor; Signal conditioner etc.); Signal (data) collector; The computing machine and the panalarm of communicating by letter.The support coordinate of each monitored amount, each Cable Structure, 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 is responsible for the health monitoring software of the cable system of operation Cable Structure, comprises the signal that the transmission of tracer signal collector comes; When monitoring rope when damage is arranged, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 9th step: with the current initial value vector of monitored amount C t o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D uParameter is kept on the hard disc of computer of operation health monitoring systems software with the mode of data file.
The tenth step: establishment and when installation and operation has support settlement and temperature variation on computers based on the damaged cable recognition methods system software of hybrid monitoring, this software will be accomplished functions such as monitoring that this method damaged cable recognition methods of hybrid monitoring " when support settlement and temperature variation are arranged based on " required by task wants, record, control, storage, calculating, notice, warning (in this practical implementation method all can with the work of computing machine completion)
The 11 step: the current numerical value vector of the monitored amount of foundation C is with the current initial value vector of monitored amount C t o, unit damage monitored quantitative change matrix Δ C, unit damage scalar D uAnd the linear approximate relationship (formula (8)) that exists between the vectorial d of the current name damage of cable system (forming) by all Suo Dangqian name amount of damage; Calculate the noninferior solution of the vectorial d of the current name damage of cable system according to multi-objective optimization algorithm, just have reasonable error but can from all ropes, confirm the position of damaged cable and separating of nominal degree of injury thereof more exactly.
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 SUm 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; This implementation step is that example provides the process of finding the solution the vectorial d of current damage with the goal programming method only, and 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 (8) can transform the multi-objective optimization question shown in an accepted way of doing sth (16) and the formula (17); γ is a real number in the formula (16); R is a real number field, and area of space Ω has limited the span (each element of present embodiment requirements vector d is not less than 0, is not more than 1) of each element of vectorial d.The meaning of formula (16) is to seek the real number γ of a minimum, makes formula (17) be met.G (d) is defined by formula (18) in the formula (17), the deviation that allows between middle G (d) of the product representation formula (17) of weighing vector W and γ and the vectorial g in the formula (17), and the definition of g is referring to formula (13), and its value calculates in the 7th step.Vector W can be identical with vectorial g during actual computation.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 the current name of cable system.
minimize γ (16)
γ∈R,d?∈Ω
G(d)-Wγ≤g (17)
G ( d ) = abs ( 1 D u ΔC · d - C + C o t ) - - - ( 18 )
The element number of the vectorial d of the current name damage of cable system equals the quantity of rope; Between the element of the vectorial d of the current name damage of cable system and the rope is one-to-one relationship, and the element numerical value of the vectorial d of the current name damage of cable system is represented the nominal degree of injury or the nominal health status of corresponding rope.Vector d the coding rule and the vectorial d of element oThe coding rule of element identical.
The 12 step: the current actual damage vector of definition cable system d a, the current actual damage vector of cable system d aElement number equal to support the quantity of rope, the current actual damage vector of cable system d aElement and supporting be one-to-one relationship between the rope, the current actual damage vector of cable system d aElement numerical value represent the actual damage degree or the actual health status of corresponding supporting rope; Vector d aThe coding rule and the vectorial d of element oThe coding rule of element identical.The current actual damage vector of the cable system that utilizes formula (15) to express d aJ element d a jWith cable system initial damage vector d oJ element d OjJ element d with the vectorial d of the current name damage of cable system jBetween relation, calculate the current actual damage of cable system vector d aAll elements; d a jBeing to represent j root supporting rope not damaged at 0 o'clock, is to represent that this supporting rope thoroughly lost load-bearing capacity at 100% o'clock, representes in the time of between 0 and 100% that the j root supports the load-bearing capacity that rope is lost corresponding proportion; That is to say the current actual damage vector of cable system d aElement numerical value represent the degree of injury of corresponding supporting rope, so according to the current actual damage vector of cable system d aCan define the impaired and degree of injury of which rope, promptly realize the damaged cable identification or the health monitoring of cable system in the Cable Structure.
The 13 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 14 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 15 the step: got back to for the 4th step, the beginning by the 4th go on foot the 15 the step circulation.

Claims (1)

  1. One kind when support settlement and temperature variation are arranged based on the damaged cable recognition methods of hybrid monitoring, it is characterized in that said method comprises:
    A. establish total N root supporting rope, at first confirm the coding rule of supporting rope, with supporting rope numberings all in the Cable Structure, this numbering will be used to generate the vector sum matrix in subsequent step by this rule; The supporting rope with monitored Suo Li of appointment when confirming hybrid monitoring is established total N root supporting rope in the cable system, and the monitored rope force data of Cable Structure is by M on the Cable Structure 1The M of individual appointment supporting rope 1Individual rope force data is described, and the variation of Cable Structure Suo Li is exactly all variations of specifying the Suo Li of supporting rope; Each total M 1Individual cable force measurement value or calculated value characterize the rope force information of Cable Structure; M 1Be one and be not less than 0 integer; The measured point with monitored strain of appointment when confirming hybrid monitoring, the monitored strain data of Cable Structure can be by K on the Cable Structure 2L individual specified point, that reach each specified point 2The strain of individual assigned direction is described, and the variation of Cable Structure strain data is exactly K 2The variation of the tested strain of all of individual specified point; Each total M 2Individual strain measurement value or calculated value characterize Cable Structure strain, M 2Be K 2And L 2Long-pending; M 2Be to be not less than 0 integer; The measured point with monitored angle of appointment when confirming hybrid monitoring, the monitored angle-data of Cable Structure is by K on the Cable Structure 3L individual specified point, that cross each specified point 3H individual appointment straight line, each appointment straight line 3Individual angle coordinate component is described, and the variation of Cable Structure angle is exactly variations all specified points, all appointments angle coordinate components straight line, all appointments; Each total M 3Individual angle coordinate component measurement value or calculated value characterize the angle information of Cable Structure, M 3Be K 3, L 3And H 3Long-pending; M 3Be one and be not less than 0 integer; When confirming hybrid monitoring appointment with monitored shape data, the monitored shape data of Cable Structure is by K on the Cable Structure 4L individual specified point, that reach each specified point 4The volume coordinate of individual assigned direction is described, and the variation of Cable Structure shape data is exactly K 4The variation of all coordinate components of individual specified point; Each total M 4Individual measurement of coordinates value or calculated value characterize Cable Structure shape, M 4Be K 4And L 4Long-pending; M 4Be one and be not less than 0 integer; The monitored amount of comprehensive above-mentioned hybrid monitoring, total M the monitored amount of whole Cable Structure, M is M 1, M 2, M 3And M 4Sum, definition parameter K, K is M 1, K 2, K 3And K 4Sum, K and M must not be less than the quantity N of rope; For simplicity, in the method " all monitored parameters of Cable Structure " are abbreviated as " monitored amount "; 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 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 of the health status of expressing the supporting rope, 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; The initial value of all monitored amounts is formed monitored amount initial value vector C oUtilize the Non-Destructive Testing data of the health status that can express the supporting rope to set up cable system initial damage vector d o, cable system initial damage vector d oElement number equal N, d oElement be one-to-one relationship with the supporting rope, cable system initial damage vector d oElement numerical value represent the degree of injury of corresponding supporting rope, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the pairing supporting rope of this element is intact, not damage; If its numerical value is 100%, represent that then the pairing supporting rope of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%; Represent that then this supporting rope lost the load-bearing capacity of corresponding proportion; If when not supporting Non-Destructive Testing data and other of rope and can express the data of health status of supporting rope, perhaps think when the Cable Structure original state is the not damaged state vectorial d oEach element numerical value get 0; Initial Cable Structure support coordinate data are formed initial Cable Structure support coordinate vector U o
    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, the employed various materials of Cable Structure and mechanical property parameters, initial Cable Structure support coordinate vector U o, 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 oSupporting rope health status with cable system initial damage vector d oExpression; Corresponding to A oThe initial value of all monitored amounts with monitored amount initial value vector C oExpression; Set up the current initial Mechanics Calculation benchmark model A of the Cable Structure that counts " Cable Structure steady temperature data " for the first time t o, the current initial value of monitored amount vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up for the first time the current initial Mechanics Calculation benchmark model A of Cable Structure t oWith the current initial value vector of monitored amount C t oThe time, the current initial Mechanics Calculation benchmark model A of Cable Structure t oJust equal the initial Mechanics Calculation benchmark model A of Cable Structure o, the current initial value vector of monitored amount C t oJust equal monitored amount initial value vector C oA t oCorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up the current initial Mechanics Calculation benchmark model A of Cable Structure the first time t oThe time, T t oJust equal T oCurrent initial Mechanics Calculation benchmark model A corresponding to Cable Structure t oCable Structure support coordinate data form current initial Cable Structure support coordinate vector U t o, set up the current initial Mechanics Calculation benchmark model A of Cable Structure the first time t oThe time, U t oJust equal U oA t oThe initial health and the A of supporting rope oThe health status of supporting rope identical, also use cable system initial damage vector d oExpression, A in the cyclic process of back t oThe initial health of supporting rope use cable system initial damage vector d all the time oExpression; T o, U 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 T t o, U t oAnd d oBe A t oParameter, C t oBy A t oMechanics Calculation result form; A in the method o, U o, C o, d oAnd T oBe constant;
    E. go on foot the m circulation in step from getting into here by e; In Cable Structure military service process; Constantly constantly survey and calculate the current data that obtains " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method "; The current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tDefinition mode and vector T oDefinition mode identical; Obtain current cable structure steady temperature data vector T in actual measurement tSynchronization, actual measurement obtains Cable Structure support coordinate current data, all Cable Structure support coordinate current datas are formed current cable structures actual measurement support coordinate vector U t
    F. according to current cable structure actual measurement support coordinate vector U tWith current cable structure steady temperature data vector T t, upgrade current initial Mechanics Calculation benchmark model A according to step f1 to f3 t o, current initial Cable Structure support coordinate vector U t o, the current initial value of monitored amount vector C t oWith current initial Cable Structure steady temperature data vector T t o
    F1 compares U respectively tWith U t o, T tWith T t oIf, U tEqual U t oAnd T tEqual T t o, A then t o, U t o, C t oAnd T t oRemain unchanged; Otherwise need follow these steps to A t o, U t oAnd T t oUpgrade;
    F2. calculate U tWith U oPoor, U tWith U oDifference be exactly the current bearing translational displacement of Cable Structure bearing about initial position, represent the bearing translational displacement with bearing translational displacement vector V, V equals U tDeduct U oBetween element among the bearing translational displacement vector V and the bearing translational displacement component is one-to-one relationship; The numerical value of an element is corresponding to the translational displacement of an assigned direction of an appointment bearing among the bearing translational displacement vector V, and wherein the bearing translational displacement is exactly the support settlement amount at the component of gravity direction; Calculate T tWith T oPoor, T tWith T oDifference be exactly of the variations of current cable structure steady temperature data about initial Cable Structure steady temperature data, T tWith T oThe difference represent that with steady temperature change vector S S equals T tDeduct T o, S representes the variation of Cable Structure steady temperature data;
    F3. earlier to A oIn the Cable Structure bearing apply the constraint of current bearing translational displacement, the numerical value of current bearing translational displacement constraint is just taken from the numerical value of corresponding element among the bearing translational displacement vector V, again 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 Cable Structure bearing apply bearing translational displacement constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial Mechanics Calculation benchmark model A that obtains upgrading t o, upgrade A t oThe time, U t oAll elements numerical value is also used U tAll elements numerical value is corresponding to be replaced, and has promptly upgraded U t o, T t oAll elements numerical value is also used T tCorresponding replacement of all elements numerical value, promptly upgraded T t o, so just obtained correctly corresponding to A t oT t oAnd U t oUpgrade C t oMethod be: when upgrading A t oAfter, obtain A through Mechanics Calculation t oIn concrete numerical value all monitored amounts, current, these concrete numerical value are formed C t oA t oThe initial health of supporting rope use cable system initial damage vector d all the time oExpression;
    G. at current initial Mechanics Calculation benchmark model A t oThe basis on carry out the several times Mechanics Calculation according to step g 1 to g4, through calculate obtaining Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D u
    G1. Cable Structure unit damage monitored quantitative change matrix Δ C brings in constant renewal in, and is promptly upgrading current initial Mechanics Calculation benchmark model A t o, current initial Cable Structure support coordinate vector U t o, the current initial value of monitored amount vector C t oWith current initial Cable Structure steady temperature data vector T t oAfterwards, must then upgrade Cable Structure unit damage monitored quantitative change matrix Δ C and unit damage scalar D u
    G2. at the current initial Mechanics Calculation benchmark model A of Cable Structure t oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all ropes on the calculation times numerical value, have N root supporting rope that N calculatings is just arranged, calculating each time supposes to have only a supporting rope that unit damage scalar D is arranged in the cable system uThe rope that occurs damage in calculating each time is different from the rope that occurs damage in other time calculating; Calculate the current calculated value of all monitored amounts in the 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, the element coding rule of monitored amount calculation current vector and monitored amount initial value vector C oThe element coding rule identical;
    G3. the monitored amount calculation current vector that calculates each time deducts the current initial value vector of monitored amount C t oObtain a monitored quantitative changeization vector; There is N root supporting rope that N monitored quantitative changeization vector just arranged;
    G4. form the Cable Structure unit damage monitored quantitative change matrix Δ C that the N row are arranged successively by this N monitored quantitative change vector; Each row of Cable Structure unit damage monitored quantitative change matrix Δ C are corresponding to a monitored quantitative changeization vector;
    H. obtain current cable structure steady temperature data vector T in actual measurement tThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T tThe current measured value of all monitored amounts of Cable Structure of synchronization in the moment, form the current numerical value vector of monitored amount C; Current numerical value vector C of monitored amount and the current initial value vector of monitored amount C t oWith monitored amount initial value vector C oDefinition mode identical, the same monitored amount of the element representation of the identical numbering of three vectors is at difference concrete numerical value constantly;
    I. define the vectorial d of the current name damage of cable system; The element number of the vectorial d of the current name damage of cable system equals to support the quantity of rope; Between the element of the vectorial d of the current name damage of cable system and the supporting rope is one-to-one relationship, and the element numerical value of the vectorial d of the current name damage of cable system is represented the nominal degree of injury or the nominal health status of corresponding supporting rope; The coding rule and the vectorial d of the element of vector d oThe coding rule of element identical;
    J. the current numerical value vector of the monitored amount of foundation C is with the current initial value vector of monitored amount C t o, Cable Structure unit damage monitored quantitative change matrix Δ C, unit damage scalar D uAnd the linear approximate relationship that exists between the vectorial d of the current name damage of cable system to be asked, this linear approximate relationship can be expressed as formula 1, and other amount in the formula 1 except that d is known, finds the solution formula 1 and just can calculate the vectorial d of the current name damage of cable system;
    C = C o t + 1 D u Δ C · d Formula 1
    K. define the current actual damage vector of cable system d a, the current actual damage vector of cable system d aElement number equal to support the quantity of rope, the current actual damage vector of cable system d aElement and supporting be one-to-one relationship between the rope, the current actual damage vector of cable system d aElement numerical value represent the actual damage degree or the actual health status of corresponding supporting rope; Vector d aThe coding rule and the vectorial d of element oThe coding rule of element identical;
    L. the current actual damage vector of the cable system that utilizes formula 2 to express d aJ element d a jWith cable system initial damage vector d oJ element d OjJ element d with the vectorial d of the current name damage of cable system jBetween relation, calculate the current actual damage of cable system vector d aAll elements;
    d j a = 1 - ( 1 - d Oj ) ( 1 - d j ) Formula 2
    J=1 in the formula 2,2,3 ...., N, d a jBeing to represent j 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 in the time of between 0 and 100% that j root supporting rope loses the load-bearing capacity of corresponding proportion; The current actual damage vector of cable system d aElement numerical value represent the degree of injury of corresponding supporting rope, so according to the current actual damage vector of cable system d aCan define the impaired and degree of injury of which rope, promptly realize the damaged cable identification or the health monitoring of cable system in the Cable Structure;
    M. get back to the e step, beginning goes on foot the m circulation next time in step by e.
CN201210172652.4A 2012-05-29 2012-05-29 Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation Expired - Fee Related CN102706638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210172652.4A CN102706638B (en) 2012-05-29 2012-05-29 Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210172652.4A CN102706638B (en) 2012-05-29 2012-05-29 Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation

Publications (2)

Publication Number Publication Date
CN102706638A true CN102706638A (en) 2012-10-03
CN102706638B CN102706638B (en) 2015-04-08

Family

ID=46899581

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210172652.4A Expired - Fee Related CN102706638B (en) 2012-05-29 2012-05-29 Damaged cable identification method on basis of hybrid monitoring during support settlement and temperature variation

Country Status (1)

Country Link
CN (1) CN102706638B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793624A (en) * 2010-03-17 2010-08-04 东南大学 Progressive method for distinguishing slack supporting cable based on mixed monitoring during support settlement
CN101806663A (en) * 2010-03-17 2010-08-18 东南大学 Health monitoring method of cable system based on mixed monitoring in presence of support seat settlement
CN101806665A (en) * 2010-03-17 2010-08-18 东南大学 Method for identifying untensioned support cables based on mixed monitoring in presence of support seat settlement
CN102288428A (en) * 2011-05-13 2011-12-21 东南大学 Method for distinguishing slack cable based on hybrid monitoring in angular displacement of supporting seat
CN102297769A (en) * 2011-05-13 2011-12-28 东南大学 Health monitoring method for identifying damaged cable and supporting-base angular displacement based on hybrid monitoring
CN102305723A (en) * 2011-05-31 2012-01-04 东南大学 Health monitoring method based on mixed monitoring and identification of damaged cable and generalized displacement of support

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101793624A (en) * 2010-03-17 2010-08-04 东南大学 Progressive method for distinguishing slack supporting cable based on mixed monitoring during support settlement
CN101806663A (en) * 2010-03-17 2010-08-18 东南大学 Health monitoring method of cable system based on mixed monitoring in presence of support seat settlement
CN101806665A (en) * 2010-03-17 2010-08-18 东南大学 Method for identifying untensioned support cables based on mixed monitoring in presence of support seat settlement
CN102288428A (en) * 2011-05-13 2011-12-21 东南大学 Method for distinguishing slack cable based on hybrid monitoring in angular displacement of supporting seat
CN102297769A (en) * 2011-05-13 2011-12-28 东南大学 Health monitoring method for identifying damaged cable and supporting-base angular displacement based on hybrid monitoring
CN102305723A (en) * 2011-05-31 2012-01-04 东南大学 Health monitoring method based on mixed monitoring and identification of damaged cable and generalized displacement of support

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭良友等: "武汉长江二桥的索力、温度和应力测量", 《桥梁建设》, no. 3, 31 December 1995 (1995-12-31), pages 44 - 53 *

Also Published As

Publication number Publication date
CN102706638B (en) 2015-04-08

Similar Documents

Publication Publication Date Title
CN102706627A (en) Damaged cable and support angular displacement identification method on basis of hybrid monitoring during temperature variation
CN102735468A (en) Identification method for damaged cables based on mixed monitoring of generalized displacement and temperature change of support
CN102721560A (en) Damaged cable identification method used in case of angular displacement of support and temperature variation on basis of space coordinate monitoring
CN102749212A (en) Generalized displacement identification method for temperature variation hybrid monitored problem cables and supports
CN102706628A (en) Damaged cable and support angular displacement identification method on basis of space coordinate monitoring during temperature variation
CN102706625A (en) Damaged cable and support angular displacement identification method on basis of strain monitoring during temperature variation
CN102706666A (en) Method for identifying damaged cable and support generalized displacement based on space coordinate monitoring during temperature variation
CN102706669A (en) Damaged cable and support generalized displacement identification method based on strain monitoring of temperature change
CN102706663A (en) Method for identifying damaged cable based on angle monitoring during support generalized displacement and temperature variation
CN102706631A (en) Damaged cable and support angular displacement identification method on basis of angle monitoring during temperature variation
CN102706667A (en) Damaged cable and support generalized displacement identification method based on angular monitoring of temperature change
CN102735466A (en) Identification method for damaged cables based on mixed monitoring in change of angular support displacement and temperatures
CN102735467A (en) Method for identifying damaged rope based on angle monitoring in the presence of supporting seat angle displacement and temperature change
CN102735469A (en) Strain monitoring-based damaged cable recognition method in support angular displacement and temperature change
CN102735471A (en) Strain monitoring-based damaged cable recognition method in support generalized displacement and temperature change
CN102735470A (en) Damaged cable identification method based on cable force monitoring of support angular displacement and temperature change
CN102706577A (en) Method for identifying problem cable and support translation based on hybrid monitoring during temperature change
CN102706661A (en) Defective cable and support angular displacement identification method based on mixed monitoring of temperature change
CN102706664A (en) Damaged cable approximation identification method based on mixed monitoring of support generalized displacement and temperature change
CN102706637A (en) Damaged cable and support translation progressive identification method on basis of hybrid monitoring during temperature variation
CN102721551A (en) Problematic cable and support translation identification method based on space coordinate monitoring during temperature change
CN102706585A (en) Damaged cable identification method based on combined monitoring at moment of temperature variation
CN102735472A (en) Damaged cable identification method based on spatial coordinate monitoring of support generalized displacement temperature change
CN102706609A (en) Damaged cable identification method on the basis of angle monitoring during temperature variation
CN102706580A (en) Method for identifying damaged cable based on cable force monitoring during temperature change

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

Granted publication date: 20150408

Termination date: 20180529

CF01 Termination of patent right due to non-payment of annual fee