CN102928243A - Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support - Google Patents

Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support Download PDF

Info

Publication number
CN102928243A
CN102928243A CN2012103793716A CN201210379371A CN102928243A CN 102928243 A CN102928243 A CN 102928243A CN 2012103793716 A CN2012103793716 A CN 2012103793716A CN 201210379371 A CN201210379371 A CN 201210379371A CN 102928243 A CN102928243 A CN 102928243A
Authority
CN
China
Prior art keywords
cable
temperature
vector
data
support
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2012103793716A
Other languages
Chinese (zh)
Inventor
韩玉林
王芳
韩佳邑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN2012103793716A priority Critical patent/CN102928243A/en
Publication of CN102928243A publication Critical patent/CN102928243A/en
Pending legal-status Critical Current

Links

Abstract

The invention relates to a slack cable progressing type identification method for general displacement and temperature variation strain monitoring of a support, which is based on strain monitoring and utilizes a reference mechanical calculation model of a cable structure for determining whether a cable structure needs to be updated or not by monitoring the general displacement of the support, and monitoring the cable structure temperature, the ambient temperature and the health degree of a bearing cable. According to the approximate linear relationship existing among the current numeric vector of a value to be monitored, the current initial numeric vector of the value to be monitored, a change matrix of the unit damage value to be monitored, and the current nominal damage vector to be calculated, a noninferior solution of the current nominal damage vector can be rapidly calculated by utilizing algorithms such as a multi-objective optimization algorithm, so that a virtual damaged cable can be recognized when the general displacement and the temperature of the support are changed. After the true damaged cable is identified from the nondestructive test method and the like, the rest damaged cable is the slack support cable, and the cable length to be adjusted of the slack support cable can be determined according to the mechanical equivalent relation.

Description

The laddering recognition methods of slack line of generalized displacement of support temperature variation strain monitoring
Technical field
The 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, such as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of this class formation is that this method is " Cable Structure " with such structure representation for simplicity take rope, cable or the rod member that only bears tensile load as support unit.Along with the variation of environment temperature, the temperature of Cable Structure also can change, and when the Cable Structure temperature changed, (for example generalized displacement of support referred to that bearing is along the angular displacement around X, Y, Z axis of the displacement of the lines of X, Y, Z axis and bearing in that generalized displacement of support is arranged; Corresponding to generalized displacement of support, the bearing generalized coordinate refers to that bearing is about X, Y, the coordinate of Z axis and bearing are about X, Y, the angular coordinate of Z axis) time, the supporting system that this method is identified Cable Structure based on strain monitoring (refers 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 cable, also comprise the rod member that only bears tensile load, equally for convenient, censure all ropeway carrying-ropes and the rod member that only bears tensile load that all play supporting role with " support cable " this noun in this method, equally for convenient, censure all ropeway carrying-ropes and the rod member that only bears tensile load that all play supporting role with " support cable " this noun in this method) in damaged cable and need to adjust the support cable (truss-frame structure is just referred to the impaired rod member that only bears tensile load) of Suo Li, belong to the engineering structure health monitoring field.
Background technology
Cable system is Cable Structure (large-scale Cable Structure particularly normally; for example large-scale cable-stayed bridge, suspension bridge) key components; owing to the reason such as lax; the Suo Li of new construction completion a period of time rear support rope can change usually; the lax variation that also can cause the supporting cable force of its support cable behind the structure long service; these change the variation that all will cause structural internal force; safety to structure causes harmful effect; will cause the inefficacy of structure when serious, therefore accurately and timely identify the support cable that needs to adjust Suo Li and be very important.
Cable system is the key components of Cable Structure normally, its inefficacy usually brings the inefficacy of total, and the damaged cable (also referring to only bear as previously mentioned the rod member of tensile load) of identifying based on structural health monitoring technology in the cable system of Cable Structure is a kind of method that has potentiality.After the health status of cable system changes, can cause the variation of the measurable parameter of structure, for example the distortion of Cable Structure or strain meeting change, in fact the variation of strain has comprised the health status information of cable system, that is to say the health status that to utilize the structural strain data to judge structure, can (this method be called monitored strain " monitored amount " based on strain monitoring, the back is mentioned " monitored amount " and is just referred to monitored strain) (this method also is referred to as the support cable of unsoundness problem, refers to that support cable is impaired to identify damaged cable, lax or have both at the same time).Monitored amount is except the impact that is subjected to the cable system health status; also can be subjected to the impact of Cable Structure temperature variation (usually can occur) and Cable Structure generalized displacement of support (sedimentation is that generalized displacement is at the component of gravity direction); the Cable Structure temperature change and the condition of Cable Structure bearing generation generalized displacement under; if can realize based on the monitoring to monitored amount the identification to the support cable of unsoundness problem; safety to Cable Structure has important value, does not also have at present a kind of disclosed, effective health monitoring systems and method to solve this problem.
When Cable Structure has the generalized displacement of support 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 rationally effectively set up each monitored quantitative change with the method for the relation between the health status of all ropes in the cable system, can provides the health evaluating of more believable cable system based on the health monitoring systems of the method foundation.
Summary of the invention
Technical matters: the purpose of this method is when Cable Structure has the generalized displacement of support temperature variation, for the health monitoring problem of cable system in the Cable Structure, disclose a kind of based on health monitor method strain monitoring, that can monitor rationally and effectively cable system in the Cable Structure.
According to the Suo Li causes of change of support cable, the Suo Li of support cable can be changed the two kinds of situations that are divided into: the one, support cable has been subject to damage, and for example localized cracks and corrosion etc. have appearred in support cable; The 2nd, support cable and not damaged, but variation has also occured in Suo Li, the one of the main reasons that this variation occurs is that variation has occured the Suo Changdu (be called drift, this method specially refers to the drift of that section rope between support cable two supporting end points) under the support cable free state (this moment, Suo Zhangli claimed that also Suo Li is 0).One of fundamental purpose of this method will identify exactly drift the support cable that changes has occured, and identifies the change amount of their drift, and this change amount provides direct basis for the Suo Li adjustment of this rope.The reason that the support cable drift changes is not single, and for convenient, this method is referred to as slack line with the support cable that drift changes.Refer to the slack line recognition system with the cable system health monitoring systems in the method, refer to the slack line recognition methods with cable system health evaluating method, in other words this method " health monitoring " usually available " slack line identification " substitute.
Technical scheme: this method is comprised of three parts.Be respectively: one, " the temperature survey calculating method of the Cable Structure of this method "; Two, set up the required knowledge base of cable system health monitoring systems and parameter method, based on knowledge base (containing parameter) and the actual measurement monitored amount cable system health status appraisal procedure; Three, the software and hardware part of health monitoring systems.
The first of this method: " the temperature survey calculating method of the Cable Structure of this method ".
At first determine " the temperature survey calculating method of the Cable Structure of this method ".Because the temperature of Cable Structure may change, for example the temperature of the different parts of Cable Structure is to change along with the variation of intensity of sunshine, along with the variation of environment temperature changes, the surface of Cable Structure may be time dependent with the temperature of inside sometimes, the surface of Cable Structure may be different from the temperature of inside, the surface of Cable Structure is time dependent with the temperature difference of inside, this is just so that Mechanics Calculation and the monitoring very complex of the Cable Structure when considering temperature conditions, be the simplification problem, reduce calculated amount and reduce the measurement cost, especially in order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", and is specific as follows:
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 0 the highest temperature and the lowest temperature between rear 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to the base area revolutions constantly and the rule that revolves round the sun is determined constantly, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between rear 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 that the daily temperature at T cloudy day 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 rBe no less than between inquiry Cable Structure location and place Altitude Region temperature that the meteorological data in recent years in 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and place Altitude Region about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m." R Cable Structure surface point " got on surface in Cable Structure, the back will obtain by actual measurement the temperature of this R Cable Structure surface point, the temperature data that claims actual measurement to obtain 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 by Calculation of Heat Transfer, just claim the temperature data that calculates to be " R Cable Structure surface temperature computational data ".When " R Cable Structure surface point " got on the surface in Cable Structure, the quantity of " R Cable Structure surface point " was narrated in the back with the condition that must satisfy that distributes.From the residing minimum height above sea level of Cable Structure to the highest height above sea level, uniform choosing is no less than three different sea level elevations on Cable Structure, at each sea level elevation place that chooses, at least choose two points at the intersection place on surface level and Cable Structure surface, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of 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 the sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, especially, measure Cable Structure for support cable along each and only get a point along the direction of the Temperature Distribution of wall thickness, namely only measure the temperature of the surface point of support cable, measure all and be selected a little temperature, the temperature that records is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of 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, especially, E equals 1 for support cable, that " measures Cable Structure along the point of the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, the back will obtain by actual measurement the temperature of this HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", the temperature data that claims actual measurement to obtain is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim the temperature data that calculates to be " HBE Cable Structure is along thickness temperature computation data "; Will be 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 " in this method ".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; On-site spacious without choosing a position in the place of blocking in Cable Structure, this position should can both obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be 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 the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate.Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, the moment of survey record data is called the physical record data constantly.
Second step, Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously; Obtain being carved at sunrise the same day the constantly temperature measured data sequence of the Cable Structure place environment between rear 30 minutes of sunrise next day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged sequentially according to time order and function by the temperature measured data that was carved at sunrise the Cable Structure place environment between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T EmaxObtain the temperature of Cable Structure place environment about the rate of change of time by the temperature measured data sequence of Cable Structure place environment by conventional mathematical computations, this rate of change is also along with the time changes; Obtain being carved at sunrise the same day the constantly measured data sequence of the temperature of the sunny slope of the reference plate between rear 30 minutes of sunrise next day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged sequentially according to time order and function by the measured data of the temperature of the sunny slope that was carved at sunrise the reference plate between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day 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 the constantly Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between rear 30 minutes of sunrise next day by Real-Time Monitoring, there be R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged, the Cable Structure surface temperature measured data between rear 30 minutes is arranged sequentially according to time order and function constantly by being carved at sunrise sunrise next day on same day of a Cable Structure surface point for each Cable Structure surface temperature measured data sequence, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day 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 at sunrise the constantly maximum temperature difference numerical value between rear 30 minutes of sunrise next day R the same day, maximal value wherein is designated as Δ T SmaxObtain the temperature of each Cable Structure surface point about the rate of change of time by each Cable Structure surface temperature measured data sequence by conventional mathematical computations, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes.Obtain being carved at sunrise the same day between rear 30 minutes of the moment of sunrise next day by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculating amounts to maximum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature 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 sea level elevation 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 " 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, determine to obtain the moment of Cable Structure steady temperature data, the condition relevant with the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between rear 30 minutes between the same day, sunset constantly refer on the meteorology that base area revolutions and revolution rule determine sunset constantly, the sunset that can inquire about data or calculate each required day by conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between rear 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 rear 30 minutes during this period of time in, measure in front the environment maximum error Δ T that calculates 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 PmaxOnly needing to satisfy in second a condition and the b condition one just is called and satisfies second condition; The 3rd condition is that the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; 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 the rate of change of time; 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 rear 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 " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition of satisfying in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment of only satisfying the 6th condition in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during in constantly one of physical record data 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 be not any in constantly of physical record data in this method constantly, then get this method close to moment of mathematics those physical record data constantly that obtain Cable Structure steady temperature data 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, i.e. this Cable Structure temperature constantly temporal evolution not, 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 in the Cable Structure in the moment that obtains Cable Structure steady temperature data by conventional Calculation of Heat Transfer, 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 in front accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " of Cable Structure, the accounting temperature of HBE " measuring Cable Structure along the point of 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, 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 " R Cable Structure surface point " got on the surface in Cable Structure, 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 the temperature of any point on the Cable Structure surface be by " 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 the support cable 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 uniform along the Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " along sea level elevation all in twos the maximal value Δ h in the absolute value of the difference of the sea level elevation of adjacent Cable Structure surface point be 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; " R Cable Structure surface point " along sea level elevation the definition of adjacent Cable Structure surface point refers to only consider sea level elevation in twos the time, do not have a Cable Structure surface point in " R 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 inquiry or obtains rule at sunshine between Cable Structure location and place Altitude Region by the meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, find in Cable Structure to be subjected to the sunshine-duration position of sufficient those surface points the whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is an annual point that is subjected in the most sufficient those surface points of sunshine-duration on the Cable Structure.
The second portion of this method: set up the required knowledge base of cable system health monitoring systems and parameter method, based on knowledge base (containing parameter) and the actual measurement monitored amount cable system health status appraisal procedure.Can carry out successively as follows, to obtain the more accurately health status assessment of cable system.
The first step: establish total N root support cable, at first determine the coding rule of support cable, with support cable numberings all in the Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step by this rule." the whole monitored strain data of structure " can be described by the strain specified point of K on the structure, that reach L assigned direction of each specified point, and the variation of structural strain data is exactly the variation of all strains of K specified point.(the individual strain measurement value of M=K * L) or calculated value characterize structural strain information to each total M.K and M must not be less than the quantity N of support cable.
For simplicity, in the method with " the monitored strain data of structure " referred to as " monitored amount ".When mentioning in the back " so-and-so matrix of monitored amount or so-and-so vector ", also can be read as " strain so-and-so matrix or so-and-so vector ".
Set up initial mechanical calculating benchmark model A oThe time, 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 measure with conventional thermometry, for example use thermal resistance to measure), this moment " Cable Structure steady temperature data " use vector T oExpression is called initial Cable Structure steady temperature data vector T oObtain T in actual measurement oThe time, namely obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, use conventional method directly to measure the initial number of all monitored amounts that calculate Cable Structure.Obtain initial Cable Structure steady temperature data vector T at Actual measurement oThe time, use conventional method (consult reference materials or survey) to obtain temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T at Actual measurement oThe time, namely obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, use the conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise that Non-destructive Testing Data of support cable etc. can express the 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 bearing generalized coordinate data, Cable Structure modal data, structural strain data, structure angle measurement of coordinates 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 determine according to these coordinate datas the geometric properties of Cable Structure.For 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 data of utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope are set up cable system initial damage vector d o(as the formula (1)), use d oThe expression Cable Structure is (with initial mechanical calculating benchmark model A oThe initial health of cable system expression).If when not having the data of the Non-destructive Testing Data of rope and other health status that can express support cable, 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 support cable, temperature variant physical and mechanical properties parameter 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 mechanical calculating benchmark model A oCorresponding to A oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized 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 mechanical calculating benchmark model A of expression oIn the initial damage value of j root rope of cable system, d OjBe to represent j root rope not damaged at 0 o'clock, represent when being 100% that this rope thoroughly loses load-bearing capacity, represent the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, T represents the transposition (rear with) of vector.
Obtain T in actual measurement oThe time, namely obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, 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 in the narration of back, will represent this calculated value and measured value reliably close to the measured data of initial monitored amount with prosign.
C o=[C o1?C o2···C ok···C oM] T????????????????????(2)
C in the formula (2) Ok(k=1,2,3 ...., M) be k monitored amount in the Cable Structure.Vector C oBe by the monitored amount of M according to certain arranged sequentially forming, this is put in order there is no specific (special) requirements, only require all associated vector of back also in this order array data get final product.
No matter which kind of method to obtain initial mechanical calculating 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%.But such utility A 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 occuring near institute's analog case reliably.Model A oThe health status of middle support cable cable system initial damage vector d oExpression, Cable Structure 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 oCan say T o, U oAnd d oA oParameter, C oBy A oMechanics Calculation result form.
Second step: circulation beginning.During each time circulation beginning, the current initial damage vector of the cable system d when at first needing to set up or set up this circulation beginning i o(i=1,2,3 ...), set up the current initial mechanical calculating benchmark model A of Cable Structure i o(finite element benchmark model for example, A in each time circulation i oConstantly update), A i oTemperature Distribution with " current initial Cable Structure steady temperature data vector T i o" express.Letter i is except the place that represents significantly number of steps, and alphabetical i only represents cycle index in the method, i.e. the i time circulation.A oAnd A i oCount temperature parameter, can accounting temperature have changed the Effect on Mechanical Properties to Cable Structure.
The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d i o(as the formula (3)), use d i oCable Structure is (with current initial mechanical calculating benchmark model A when representing this time circulation beginning i oThe health status of cable system expression).
d o i = d o 1 i d o 2 i . . . d oj i . . . d oN i T - - - ( 3 )
D in the formula (3) i Oj(i=1,2,3, J=1,2,3 ...., when N) the i time circulation of expression begins, current initial mechanical calculating benchmark model A i oIn the initial damage value of j root rope of cable system, d i OjBe to represent j root rope not damaged at 0 o'clock, represent when being 100% that this rope thoroughly loses load-bearing capacity, represent the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%.
Current initial mechanical calculating benchmark model A corresponding to Cable Structure i oCable Structure bearing generalized coordinate data form current initial Cable Structure bearing generalized coordinate vector U i o, namely set up for the first time the current initial mechanical calculating benchmark model A of Cable Structure at initial time i oThe time, U i oJust equal U o
Set up and renewal d i oMethod as follows:
During for the first time circulation beginning, set up the current initial damage vector of cable system and (be designated as d according to formula (3) 1 o) time, d 1 oJust equal d oI (i=2,3,4,5,6 ...) the current initial damage vector of the cable system of needs d when inferior circulation begins i o, be front once (namely the i-1 time, i=2,3,4,5,6 ...) the front calculating acquisition of circulation end, concrete grammar is described below.
I (i=1,2,3,4,5,6 ...) the Mechanics Calculation benchmark model of the Mechanics Calculation benchmark model that need to set up when inferior circulation begins or the Cable Structure of having set up is designated as current initial mechanical calculating benchmark model A i oCorresponding to A i o" Cable Structure steady temperature data " use vector T i oExpression is called current initial Cable Structure steady temperature data vector T i oVector T i oDefinition mode and vector T oDefinition mode identical, when beginning circulation must be set up or set up and be called current initial Cable Structure steady temperature data vector T each time i o
Set up, upgrade A i o, U i oAnd T i oMethod as follows:
The Mechanics Calculation benchmark model of the Cable Structure of setting up during for the first time circulation beginning is designated as A 1 o, A 1 oEqual A o, T 1 oEqual T o, U 1 oEqual U oA in each time circulation i o, U i oAnd T i oConstantly update, concrete grammar is described below; When each time circulation end, upgrade A i o, U i oAnd T i oThe Mechanics Calculation benchmark model of required Cable Structure when next time being circulated beginning, concrete grammar is described below.
This method " the current initial value vector of monitored amount C i o" (i=1,2,3 ...) initial value (referring to formula (4)) of the monitored amount of all appointments when the i time (i=1,2,3,4,5,6 ...) circulation of expression begins, C i oAlso can be called " the i time current initial value of the monitored amount of circulation vector ".
C o i = C o 1 i C o 2 i . . . C ok i . . . C oM i T - - - ( 4 )
C in the formula (2) i Ok(i=1,2,3, K=1,2,3 ...., k monitored amount when M) being the i time circulation beginning, in the Cable Structure.Vector C i oBe by the monitored amount of previously defined M according to certain arranged sequentially forming, this is put in order there is no specific (special) requirements, only require all associated vector of back also in this order array data get final product.
Setting up model A i oThe time set up " the current initial value vector of monitored amount C i o", the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts, C i oElement and C oElement corresponding one by one, represent that respectively all monitored amounts are in A in Cable Structure i oAnd A oConcrete numerical value during two states.
Set up and renewal C i oConcrete grammar as follows:
During for the first time circulation beginning, C 1 o(i=1, C i oBe embodied as C 1 o) equal C oI (i=2,3,4,5,6 ...) " the current initial value vector of the monitored amount C of the i time circulation of needs when inferior circulation begins i o", be front once (namely the i-1 time, i=2,3,4,5,6 ...) circulation calculates before finishing and to obtain, concrete grammar is described below.The i time (i=1,2,3,4,5,6 ...) in the circulation, " the current initial value vector of monitored amount C i o" constantly update, concrete grammar is described below.Because according to model A i oCalculate the initial value of the monitored amount of gained reliably close to corresponding measured value, in the narration of back, will represent this calculated value composition of vector and measured value composition of vector with prosign.
Can say T i o, U i oAnd d i oA i oCharacterisitic parameter, C i oA i oAt T i o, U i oAnd d i oMechanics Calculation result under the condition forms.
The 3rd step: in Cable Structure military service process, in each time circulation, in other words in i (i=1,2,3,4,5,6 ...) in the inferior circulation, at known A i o, T i o, U i o, C i oAnd d i oAfter, according to the current data of " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement acquisition " Cable Structure steady temperature data ", the current data of all " Cable Structure steady temperature data " forms " current cable structure steady temperature data vector T i", vector T iDefinition mode and vector T oDefinition mode identical; In the actual measurement vector T iThe time, namely obtaining current cable structure steady temperature data vector T iThe synchronization in the moment, actual measurement obtains the currency of all monitored amounts in the Cable Structure, all these numerical value form monitored amount current value vector C iC iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.Obtaining vector T iThe time, actual measurement obtains Cable Structure bearing generalized coordinate current data, and all Cable Structure bearing generalized coordinate current datas form current cable structure actual measurement bearing generalized coordinate vector U i
Obtaining vector T iAfter, upgrade A according to following concrete grammar i o, T i o, U i o, C i oAnd d i o:
Compare respectively T iAnd T i o, U iAnd U i oIf, T iEqual T i oAnd U iEqual U i o, then do not need A i oUpgrade, otherwise need to be to A i o, U i oAnd T i oUpgrade, update method is: the first step is calculated U iWith U oPoor, U iWith U oDifference be exactly that the Cable Structure bearing is about the current generalized displacement of support of initial position, V represents generalized displacement of support with the generalized displacement of support vector, be one-to-one relationship between element among the generalized displacement of support vector V and the generalized displacement of support component, the numerical value of an element is corresponding to the displacement of an assigned direction of an appointment bearing among the generalized displacement of support vector V; Second step calculates T iWith T oPoor, T iWith T oDifference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T iWith T oPoor represent that with steady temperature change vector S S equals T iDeduct T o, S represents the variation of Cable Structure steady temperature data; The 3rd step is first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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 oMiddle Cable Structure bearing applies generalized displacement of support constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial mechanical calculating benchmark model A that obtains upgrading i o, upgrade A i oThe time, U i oAll elements numerical value is also used U iAll elements numerical value is corresponding to be replaced, and has namely upgraded U i o, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, namely upgraded T i o, so just obtained correctly corresponding to A i oT i oThis moment d i oRemain unchanged.When upgrading A i oAfter, A i oThe health status of rope with the current initial damage of cable system vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iExpression, A i oThe bearing generalized coordinate with current initial Cable Structure bearing generalized coordinate vector U i oExpression obtains A by Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, replace C with these concrete numerical value i oThe element of middle correspondence has so just been realized the current initial value vector of monitored amount C i oRenewal.
The 4th step: circulation time must be set up first " unit damage monitored numerical quantity transformation matrices " and " nominal unit damage vector " each time, and " unit damage monitored numerical quantity transformation matrices " that the i time circulation set up is designated as Δ C i(i=1,2,3 ...)." nominal unit damage vector " that the i time circulation set up is designated as D i uΔ C in each time circulation iAnd D i uNeed according to circumstances to constantly update, namely upgrading current initial mechanical calculating benchmark model A i o, current initial Cable Structure steady temperature data vector T i oWith the current initial value vector of monitored amount C i oAfter, upgrade unit damage monitored numerical quantity transformation matrices Δ C iWith nominal unit damage vector D i u
Set up in the steps below first unit damage monitored numerical quantity transformation matrices Δ C during each time circulation beginning iWith nominal unit damage vector D i uIf in the 3rd step, upgraded A i o, in this step, must rebulid so (namely upgrading) unit damage monitored numerical quantity transformation matrices Δ C iWith nominal unit damage vector D i uIf in the 3rd step, do not upgrade A i o, in this step, needn't rebulid so unit damage monitored numerical quantity transformation matrices Δ C iWith nominal unit damage vector D i uSet up and rebulid (namely upgrading) Δ C iAnd D i uDetailed process identical, be listed as follows:
Current initial mechanical calculating benchmark model A in Cable Structure i oThe basis on carry out several times and calculate, equal the quantity of all ropes on the calculation times numerical value.Calculating each time in the hypothesis cable system only has a rope to increase unit damage (for example getting 5%, 10%, 20% or 30% equivalent damage is unit damage) on the basis of original damage (original damage can be 0, can not be 0 also) again.Calculate for convenient, when setting unit damage in the circulation each time can all be structural health conditions during this time circulation beginning as being fully healthy, and set on this basis unit damage (in subsequent step, damage numerical value that calculate, rope---be called name damage d i c(i=1,2,3 ...), all with respect to this time when beginning circulation, with the health status of rope as being fully healthy speech, the name that the formula that therefore must foundation hereinafter provides will calculate is damaged and is converted into true damage).The rope of appearance damage is different from the rope that appearance damages in other time calculating in the each time calculating that once circulates together, and supposes that each time the unit damage value of the rope that damage is arranged can be different from the unit damage value of other ropes, with " the vectorial D of nominal unit damage i u" (as the formula (5)) record the unit damage of the supposition of all ropes in each time circulation, circulation time is designated as D for the first time 1 uCalculate each time all utilize mechanics method (for example finite element method) calculate Cable Structure, the current calculated value of the M of appointment monitored amount in front, the current calculated value that calculates each time gained M monitored amount forms one " monitored amount is calculated the current value vector ", and (when hypothesis j root rope had unit damage, available formula (6) represented that the monitored amount of M monitored amount of all appointments calculates current value vector C i Tj); The monitored amount that calculates is each time calculated the current value vector and is deducted the current initial value vector of monitored amount C i o, the gained vector is exactly that " the monitored numberical value of quantity change vector " of (take the position of rope that unit damage is arranged or numbering etc. as mark) (when j root rope has unit damage, used δ C under this condition i jRepresent monitored numberical value of quantity change vector, δ C i jDefinition see formula (7), formula (8) and formula (9), formula (7) deducts after the formula (4) again divided by vectorial D for formula (6) i uJ element D i UjGained), monitored numberical value of quantity change vector δ C i jEach element representation since when calculating supposition the unit damage (D for example of the Na Gensuo (for example j root rope) of unit damage is arranged i Uj), and the corresponding monitored numberical value of quantity change amount of this element that causes is with respect to the unit damage D of supposition i UjRate of change; There is N root rope that N " monitored numberical value of quantity change vector " just arranged, each monitored numberical value of quantity change vector has M element, forms successively " the unit damage monitored numerical quantity transformation matrices Δ C that M * N element arranged by this N " monitored numberical value of quantity change vector " i" (the capable N row of M), each vectorial δ C i j(j=1,2,3 ...., N) be matrix Δ C iOne row, Δ C iDefinition as the formula (10).
D u i = D u 1 i D u 2 i . . . D uj i . . . D uN i T - - - ( 5 )
Nominal unit damage vector D in the formula (5) i uElement D i Uj(i=1,2,3, J=1,2,3 ...., N) the unit damage numerical value of the j root rope of supposition in the i time circulation of expression, vectorial D i uIn the numerical value of each element can be the same or different.
C tj i = C tk 1 i C tk 2 i . . . C tjk i . . . C tjM i T - - - ( 6 )
Elements C in the formula (6) i Tjk(i=1,2,3, J=1,2,3 ...., N; K=1,2,3 ...., M) the i time circulation of expression be owing to j root rope when unit damage is arranged, according to the calculating current value of the monitored amount of corresponding k the appointment of coding rule.
δC j i = C tj i - C o i D uj i - - - ( 7 )
The subscript i(i=1 of each amount in the formula (7), 2,3 ...) the i time circulation of expression, subscript j (j=1,2,3 ...., N) expression j root rope has unit damage, D in the formula i UjVectorial D i uIn j element.Vector δ C i jDefinition as the formula (8), δ C i jK(k=1,2,3 ...., M) individual element δ C i JkRepresent to set up matrix Δ C in the i time circulation iThe time, suppose that the change amount of calculating a gained k monitored amount when j root rope has unit damage is with respect to the unit damage D that supposes i UjRate of change, it defines suc as formula shown in (9).
δ C j i = δC j 1 i δC j 2 i . . . δ C jk i . . . δC jM i T - - - ( 8 )
δC jk i = C tjk i - C ok i D uj i - - - ( 9 )
Each definition of quantity has been previously described in the formula (9).
δ C i = δC 1 i δC 2 i . . . δC j i . . . δC N i - - - ( 10 )
Vectorial δ C in the formula (10) i j(i=1,2,3 ....,, j=1,2,3 ...., N) in the i time circulation of expression, because j root rope has unit damage D i UjCause, the relative value of all monitored amounts changes.Matrix Δ C iThe coding rule of row (subscript j) and front vector d i oThe coding rule of subscript j of element identical.
The 5th step: the current health status of identification cable system.Detailed process is as follows.
I(i=1,2,3 ...) in the inferior circulation, utilize " the monitored amount current value vector C that obtains in the actual measurement of the 3rd step i" same " the current initial value of monitored amount vector C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between linear approximate relationship, shown in (11) or formula (12).
C i = C o i + ΔC i · d c i - - - ( 11 )
C i - C o i = ΔC i · d c i - - - ( 12 )
Monitored amount current value vector C in formula (11) and the formula (12) iDefinition be similar to the current initial value of monitored amount vector C i oDefinition, see formula (13); The vectorial d of the current name damage of cable system i cDefinition see formula (14).
C i = C 1 i C 2 i . . . C k i . . . C M i T - - - ( 13 )
Elements C in the formula (13) i k(i=1,2,3 ....; K=1,2,3 ...., M) be the i time circulation time Cable Structure, according to the current value of the monitored amount of the corresponding k of being numbered of coding rule.
d c i = d c 1 i d c 2 i . . . d cj i . . . d cN i T - - - ( 14 )
D in the formula (14) i Cj(i=1,2,3 ....; J=1,2,3 ...., N) be the current nominal impairment value of cable system j root rope in the i time circulation, vectorial d i cCoding rule and the matrix Δ C of subscript j of element iThe coding rule of row identical.
When the rope actual damage was not too large, because the Cable Structure material still is in the linear elasticity stage, the distortion of Cable Structure was also less, and the represented a kind of like this linear relationship of formula (11) or formula (12) is less with the error of actual conditions, and error can be used error vector e i(formula (15)) definition, the error of linear relationship shown in expression (11) or the formula (12).
e i = abs ( ΔC i · d c i - C i + C o i ) - - - ( 15 )
Abs () is the function that takes absolute value in the formula (15), and each element of the vector of trying to achieve in the bracket is taken absolute value.
Because there are certain error in formula (11) or the represented linear relationship of formula (12), therefore can not be simply according to formula (11) or formula (12) and " monitored amount current value vector C i" come direct solution to obtain the vectorial d of Suo Dangqian name damage i cIf done like this, the damage vector d that obtains i cIn element in addition larger negative value can appear, namely negative damage, this obviously is irrational.Therefore obtaining rope damages vectorial d i cAcceptable solution (namely with reasonable error, but can be more accurately from cable system, determining position and the degree of injury thereof of damaged cable) become a rational solution, available formula (16) is expressed this method.
abs ( ΔC i · d c i - C i + C o i ) ≤ g i - - - ( 16 )
Abs () is the function that takes absolute value in the formula (16), vectorial g iDescription departs from the legitimate skew of ideal linearity relation (formula (11) or formula (12)), is defined by formula (17).
g i = g 1 i g 2 i . . . g k i . . . g M i T - - - ( 17 )
G in the formula (17) i k(i=1,2,3 ....; K=1,2,3 ...., M) maximum allowable offset that departs from the ideal linearity relation shown in formula (11) or the formula (12) in the i time circulation has been described.Vector g iCan be according to the error vector e of formula (15) definition iTentative calculation is selected.
At the current initial value vector of monitored amount C i o, unit damage monitored numerical quantity transformation matrices Δ C iWith monitored amount current value vector C iWhen known, can utilize suitable algorithm (for example multi-objective optimization algorithm) to find the solution formula (16), obtain the vectorial d of the current name damage of cable system i cAcceptable solution, the current actual damage of cable system vector d iThe element of (formula (18) is seen in definition) can calculate according to formula (19), thereby can be by d iDetermine position and the degree of injury of damaged cable, namely realized the health monitoring of cable system, realized damaged cable identification.
d i = d 1 i d 2 i . . . d j i . . . d N i T - - - ( 18 )
D in the formula (18) i j(i=1,2,3, J=1,2,3 ...., N) the actual damage value of j root rope in the i time circulation of expression, formula (19), d are seen in its definition i jBe to represent j root rope not damaged at 0 o'clock, represent when being 100% that this rope thoroughly loses load-bearing capacity, represent the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, vectorial d iThe coding rule of element and formula (1) in vectorial d oThe coding rule of element identical.
d j i = 1 - ( 1 - d oj i ) ( 1 - d cj i ) - - - ( 19 )
D in the formula (19) i Oj(i=1,2,3,4, J=1,2,3 ...., N) be the current initial damage vector of cable system d i oJ element, d i CjThe vectorial d of the current name damage of cable system i cJ element.
The below's narration has obtained the current actual damage vector of cable system d iAfter, the position of how to confirm slack line and relax level.
If total N root support cable in the cable system, structure rope force data is described by the Suo Li of N root support cable.Available " initial rope force vector F o" the initial Suo Li (formula (20) is seen in definition) of all support cables in the expression Cable Structure.
F o=[F o1?F o2···F oj···F oN] T????????????????????(20)
F in the formula (20) o(j=1,2,3 ...., N) being the initial Suo Li of j root support cable in the Cable Structure, this element is according to the Suo Li of coding rule corresponding to the appointment support cable.Vector F oIt is constant.Obtain T in actual measurement oSynchronization, use conventional method directly to measure the rope force data that calculates all support cables, all these rope force datas form initial rope force vector F oSetting up the initial mechanical calculating benchmark model A of Cable Structure oThe time in fact used vectorial F o
The current cable power (formula (21) is seen in definition) of all support cables in the Cable Structure that obtains with " current cable force vector F " expression actual measurement in this method.
F=[F 1?F 2···F j···F N] T?????????????????????(21)
F in the formula (21) j(j=1,2,3 ...., N) be the current cable power of j root support cable in the Cable Structure.Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains the rope force data of all support cables in the Cable Structure, all these rope force datas form current cable force vector F.Element and the vectorial F of vector F oThe coding rule of element identical.According to the narration of front, vector T i oEqual vector T i
In this method, under the support cable original state, in the steady temperature data of Cable Structure with initial Cable Structure steady temperature data vector T oDuring expression, and support cable is when being in free state (free state refers to that Suo Li is 0, rear with), and the length of support cable is called initial drift, with " initial drift vector l o" the initial drift (formula (22) is seen in definition) of all support cables in the expression Cable Structure.Pass through vector T according to " the temperature survey calculating method of the Cable Structure of this method " oCan determine obtaining vector T oThe Temperature Distribution of all support cables constantly.
l o=[l o1?l o2···l oj···l oN] T???????????????????(22)
L in the formula (22) Oj(j=1,2,3 ...., N) be the initial drift of j root support cable in the Cable Structure.Vector l oBe constant, after when beginning, determining, just no longer change.
Similarly, under the support cable original state, in the steady temperature data of Cable Structure with initial Cable Structure steady temperature data vector T oDuring expression, and support cable is when being in free state, and the cross-sectional area of support cable is called initial free cross-sectional area, with " initial free cross-sectional area vector A o" the initial free cross-sectional area (formula (23) is seen in definition) of all support cables in the expression Cable Structure, the weight of the unit length of support cable is called the weight of initial free unit length, with " the weight vector ω of initial free unit length o" weight (formula (24) is seen in definition) of initial free unit length of all support cables in the expression Cable Structure.
A o=[A o1?A o2···A oj···A oN] T????????????????(23)
A in the formula (23) Oj(j=1,2,3 ...., N) be the initial free cross-sectional area of j root support cable in the Cable Structure.Vector A oBe constant, after when beginning, determining, just no longer change.
ω o=[ω o1o2···ω oj···ω oN] T????????????????(24)
ω in the formula (24) Oj(j=1,2,3 ...., N) be the weight of the free unit length of initial freedom of j root support cable in the Cable Structure.Vector ω oBe constant, after when beginning, determining, just no longer change.
In this method, in the steady temperature data of Cable Structure with current initial Cable Structure steady temperature data vector T t oDuring expression, with " current initial drift vector l t o" the current initial drift of all support cables in the expression Cable Structure (and formula (25) is seen in definition, refers to that hypothesis supporting cable force is at 0 o'clock, considered the impact on the support cable drift of thermal expansivity and temperature variation after, initial drift vector l oWith initial Cable Structure steady temperature data vector T oThe support cable of expression in temperature with current initial Cable Structure steady temperature data vector T i oSupport cable drift during expression).According to " the temperature survey calculating method of the Cable Structure of this method ", pass through vector T i oCan determine obtaining vector T i oThe Temperature Distribution of all support cables constantly.
l o t = l o 1 t l o 2 t . . . l oj t . . . l oN t T - - - ( 25 )
L in the formula (25) t Oj(j=1,2,3 ...., be with current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure N) i oDuring expression, the current initial drift of j root support cable in the Cable Structure can be utilized thermal expansivity, the l of support cable Oj, T oAnd T i oObtain l by conventional physical computing t Oj
Element, the vectorial l of vector F oElement, vectorial l t oElement, vectorial A oElement, vectorial ω oElement and vectorial F oThe coding rule of element identical, the different information of the same support cable of element representation of the identical numbering of these vectors.
In this method, in the steady temperature data of Cable Structure with current initial Cable Structure steady temperature data vector T i oDuring expression, with the current drift (formula (26) is seen in definition, and this moment, support cable may be intact, may be impaired, also may relax) of all support cables in " current drift vector l " expression Cable Structure.
l=[l 1?l 2···l j···l N] T??????????????????(26)
L in the formula (26) j(j=1,2,3 ...., N) be the current drift of j root support cable in the Cable Structure.
In this method, represent the change amount (formula (27) and formula (28) are seen in definition) of the drift of all support cables in the Cable Structure with " drift changes vectorial Δ l " (or claiming support cable current relax level vector).
Δl=[Δl 1?Δl 2···Δl j···Δl N] T???????????????(27)
Δ l in the formula (27) j(j=1,2,3 ...., N) being the change amount of the drift of j root support cable in the current cable structure, formula (28), Δ l are seen in its definition jBe not that 0 rope is slack line, Δ l jNumerical value be the slack of rope, and the current relax level of expression cable system j root support cable also is the long adjustment amount of rope of this rope when adjusting Suo Li.
Δl j = l j - l oj t - - - ( 28 )
By slack line is carried out the relax level identification that slack line is carried out in the mechanics equivalence with damaged cable, the mechanical condition of equivalence is in the method:
The mechanics parameters of initial drift, geometrical property parameter and material when one, the nothing of the rope of two equivalences relaxes with not damaged is identical;
Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.
When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in structure is exactly identical, if after namely replacing slack line with equivalent damaged cable, Cable Structure any variation can not occur, vice versa.
Obtained the current actual damage vector of cable system d iAfter, d iJ element d i j(j=1,2,3 ...., N) the actual damage value of expression j root rope, formula (19) is seen in its definition, although with d i jBe called the actual damage value of j root rope or the actual damage degree of j root rope, but since j root Suo Keneng be impaired also may be lax, so d iJ element d i jThe actual damage value of the j root rope of expression is actually the actual equivalent damage value of j root rope, when j root rope is actually impaired, and d i jThe actual damage value of the j root rope that just represents, when j root rope is actually lax, d i jJ root rope and actual damage value lax equivalence with regard to expression for sake of convenience, claims d in the method i jBe to represent j root rope not damaged at 0 o'clock, represent when being 100% that this rope thoroughly loses load-bearing capacity, represent the load-bearing capacity of j root rope forfeiture corresponding proportion in the time of between 0 and 100%, by the current actual damage vector of cable system d iJust can identify the support cable that health status goes wrong, but some is impaired in the support cable that these health status go wrong, some is to have relaxed, if j support cable is actually (its current relax level Δ l that relaxed jThe current relax level Δ l of j the support cable that definition), relaxes so j(Δ l jDefinition see formula (27)) with the current actual damage degree d of damaged cable of equivalence i jBetween relation determined by aforementioned two mechanics equivalent conditions.Δ l jSame d i jBetween physical relationship can adopt accomplished in many ways, for example can directly determine (referring to formula (29)) according to aforementioned equivalent condition, also can adopt based on the Ernst equivalent elastic modulus to replace the E in the formula (25) to revise rear determine (referring to formula (30)), also can adopt based on other methods such as trial and error procedure of finite element method and determine.
Δl j t = d j i 1 - d j i F j E j t A j t + F j l oj t - - - ( 29 )
Δl j t = d j i 1 - d j i F j [ E j t 1 + ( ω j t l jx t ) 2 A j t E j t 12 ( F j ) 3 ] A j t + F j l oj t - - - ( 30 )
E in formula (29) and the formula (30) t jWith current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure i oDuring expression, the elastic modulus of j support cable, A t jWith current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure i oDuring expression, the cross-sectional area of j support cable, F jWith current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure i oDuring expression, the current cable power of j support cable, d i jThe current actual damage degree of j support cable, ω t jWith current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure t oDuring expression, the weight of the unit length of j support cable, l t JxWith current initial Cable Structure steady temperature data vector T in the steady temperature data of Cable Structure i oDuring expression, the horizontal range of two supporting end points of j support cable.E t jA can be obtained according to the characteristic material data of looking into or survey j support cable t jAnd ω t jCan be according to thermal expansivity, the A of j support cable Oj, ω Oj, F j, T oAnd T i oObtain by conventional physics and Mechanics Calculation.Item in the formula (30) in [] is the Ernst equivalent elastic modulus of this support cable, can just can determine the current relax level vector of support cable Δ l by formula (29) or formula (30).Formula (30) is the correction to formula (29).
The 6th step: judge whether to finish this (the i time) circulation, if so, then finish the tailing in work before this circulation finishes, for next time (namely the i+1 time, i=1,2,3,4 ...) circulation preparation Mechanics Calculation benchmark model and necessary vector.Detailed process is as follows:
In this (the i time) circulation, try to achieve the vectorial d of current name damage i cAfter, at first, set up mark vector B according to formula (31) i, formula (32) has provided mark vector B iThe definition of j element; If mark vector B iElement be 0 entirely, then got back to for the 3rd step and proceed health monitoring and calculating to cable system; If mark vector B iElement be not 0 entirely, then finish subsequent step after, enter next time circulation.
So-called subsequent step is: at first, according to formula (33) calculate next time (namely the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation I+1 oEach element d I+1 OjThe second, at Mechanics Calculation benchmark model A oThe basis on, make A oIn the health status of rope be d I+1 oRather than be d oAfter, more further to A oIn Cable Structure apply temperature variation (as previously mentioned, the numerical value of the temperature variation that applies just taken from steady temperature change vector S, and steady temperature change vector S equals T iDeduct T o), so just obtained next time (namely the i+1 time, i=1,2,3,4 ...) the required current initial mechanical calculating benchmark mould A of circulation I+1 o, next time (namely the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual T i o, to A I+1 oCarrying out Mechanics Calculation obtains corresponding to A I+1 oConcrete numerical value all monitored amounts, current, these concrete numerical value form next time (namely the i+1 time, i=1,2,3,4 ...) the current initial value vector C of the required monitored amount of circulation I+1 o
B i = B 1 i B 2 i . . . B j i . . . B N i T - - - ( 31 )
Mark vector B in the formula (31) iSubscript i represent the i time the circulation, its element B i j(j=1,2,3 ..., subscript j N) represents the damage characteristic of j root rope, can only get 0 and 1 two amount, concrete value rule is seen formula (32).
B j i = 0 , if d cj i < D uj i 1 , if d cj i &GreaterEqual; D uj i - - - ( 32 )
Element B in the formula (32) i jMark vector B iJ element, D i UjNominal unit damage vector D i uJ element (seeing formula (3)), d i CjThe vectorial d of the current name damage of cable system i cJ element (seeing formula (14)), they all represent the relevant information of j root rope.
d oj i + 1 = 1 - ( 1 - d oj i ) ( 1 - D uj i B j i ) - - - ( 33 )
D in the formula (33) i UjNominal unit damage vector D i uJ element (seeing formula (5)), d i OjThe current initial damage vector of cable system d i oJ element (seeing formula (3)).
The third part of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (the horizontal range monitoring system, the Cable Structure bearing generalized coordinate monitoring system that comprise monitored amount monitoring system, temperature monitoring system, cable force monitoring system, support cable two supporting end points), signal picker and computing machine etc.The requirement Real-Time Monitoring obtains temperature required measured data, the Suo Li of requirement while each monitored amount of Real-Time Monitoring, while each support cable of Real-Time Monitoring, the horizontal range of each support cable two supporting end points of while Real-Time Monitoring, the data of while each Cable Structure bearing generalized coordinate of Real-Time Monitoring.
Software section should be finished the process that this method sets, namely finish needed in this method, can be with functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
This method specifically comprises:
A. establish total N root support cable, at first determine the coding rule of support cable, with support cable numberings all in the Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step by this rule; Determine the monitored point of appointment, monitored point namely characterizes all specified points of Cable Structure strain information, and gives all specified point numberings; Determine monitored should the changing direction of monitored point, and give the monitored strain numbering of all appointments; " monitored strain numbering " will be for generating the vector sum matrix in subsequent step; " the whole monitored strain data of Cable Structure " is comprised of above-mentioned all monitored strains; This method with " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity of monitored point must not be less than the quantity of support cable; The quantity sum of all monitored amounts must not be less than the quantity of support cable; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, 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 " is undertaken by step b1 to b3;
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, in the method can not be seen one of the sun daytime and be called all day the cloudy day, statistics obtains 0 the highest temperature and the lowest temperature between rear 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to the base area revolutions constantly and the rule that revolves round the sun is determined constantly, do not represent necessarily can see the sun same day, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between rear 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 that the daily temperature at T cloudy day 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 rBe no less than between inquiry Cable Structure location and place Altitude Region temperature that the meteorological data in recent years in 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and place Altitude Region about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; " R Cable Structure surface point " got on surface in Cable Structure, getting the Specific Principles of " R Cable Structure surface point " narrates in step b3, the back will obtain by actual measurement the temperature of this R Cable Structure surface point, the temperature data that claims actual measurement to obtain 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 by Calculation of Heat Transfer, just claim the temperature data that calculates to be " 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, uniform choosing is no less than three different sea level elevations on Cable Structure, at each sea level elevation place that chooses, at least choose two points at the intersection place on surface level and Cable Structure surface, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of 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 the sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, especially, measure Cable Structure for support cable along each and only get a point along the direction of the Temperature Distribution of wall thickness, namely only measure the temperature of the surface point of support cable, measure all and be selected a little temperature, the temperature that records is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of 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, especially, E equals 1 for support cable, that " measures Cable Structure along the point of the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, the back will obtain by actual measurement the temperature of this HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", the temperature data that claims actual measurement to obtain is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim the temperature data that calculates to be " HBE Cable Structure is along thickness temperature computation data "; Will be 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 " in this method "; 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; On-site spacious without choosing a position in the place of blocking in Cable Structure, this position should can both obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, 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 the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate;
B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously; Obtain being carved at sunrise the same day the constantly temperature measured data sequence of the Cable Structure place environment between rear 30 minutes of sunrise next day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged sequentially according to time order and function by the temperature measured data that was carved at sunrise the Cable Structure place environment between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day on same day that minimum temperature obtains Cable Structure place environment, be called the environment maximum temperature difference, be designated as Δ T EmaxObtain the temperature of Cable Structure place environment about the rate of change of time by the temperature measured data sequence of Cable Structure place environment by conventional mathematical computations, this rate of change is also along with the time changes; Obtain being carved at sunrise the same day the constantly measured data sequence of the temperature of the sunny slope of the reference plate between rear 30 minutes of sunrise next day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged sequentially according to time order and function by the measured data of the temperature of the sunny slope that was carved at sunrise the reference plate between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day 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 the constantly Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between rear 30 minutes of sunrise next day by Real-Time Monitoring, there be R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged, the Cable Structure surface temperature measured data between rear 30 minutes is arranged sequentially according to time order and function constantly by being carved at sunrise sunrise next day on same day of a Cable Structure surface point for each Cable Structure surface temperature measured data sequence, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there be R Cable Structure surface point just to have and be carved at sunrise the constantly maximum temperature difference numerical value between rear 30 minutes of sunrise next day R the same day, maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T SmaxObtain the temperature of each Cable Structure surface point about the rate of change of time by each Cable Structure surface temperature measured data sequence by conventional mathematical computations, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes; Obtain being carved at sunrise the same day between rear 30 minutes of the moment of sunrise next day by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculating amounts to maximum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature 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 sea level elevation 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 " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
B3: measure and calculate acquisition Cable Structure steady temperature data; At first, determine to obtain the moment of Cable Structure steady temperature data, the condition relevant with the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between rear 30 minutes between the same day, sunset constantly refer on the meteorology that base area revolutions and revolution rule determine sunset constantly, the sunset that can inquire about data or calculate each required day by conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between rear 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 rear 30 minutes during this period of time in, measure in front the environment maximum error Δ T that calculates 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 PmaxOnly needing to satisfy in second a condition and the b condition one just is called and satisfies second condition; The 3rd condition is that the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; 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 the rate of change of time; 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 rear 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 " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition of satisfying in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment of only satisfying the 6th condition in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during in constantly one of physical record data 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 be not any in constantly of physical record data in this method constantly, then get this method close to moment of mathematics those physical record data constantly that obtain Cable Structure steady temperature data 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, i.e. this Cable Structure temperature constantly temporal evolution not, 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 in the Cable Structure in the moment that obtains Cable Structure steady temperature data by conventional Calculation of Heat Transfer, 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 in front accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " of Cable Structure, the accounting temperature of HBE " measuring Cable Structure along the point of 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, 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 " R Cable Structure surface point " got on the surface in Cable Structure, 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 the temperature of any point on the Cable Structure surface be by " 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 the support cable 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 uniform along the Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " along sea level elevation all in twos the maximal value Δ h in the absolute value of the difference of the sea level elevation of adjacent Cable Structure surface point be 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; " R Cable Structure surface point " along sea level elevation the definition of adjacent Cable Structure surface point refers to only consider sea level elevation in twos the time, do not have a Cable Structure surface point in " R 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 inquiry or obtains rule at sunshine between Cable Structure location and place Altitude Region by the meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, find in Cable Structure to be subjected to the sunshine-duration position of sufficient those surface points the whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is an annual point that is subjected in the most sufficient those surface points of sunshine-duration on the Cable Structure;
C. directly measure the Cable Structure steady temperature data that calculate under the original state according to " the temperature survey calculating method of the Cable Structure of this method ", 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 physical and mechanical properties parameter of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is 0 o'clock length, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables when obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and initial rope force vector F of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum oThe coding rule of element identical; Obtain T in actual measurement oThe time, namely 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 initial value of the Non-destructive Testing Data of the health status of expressing support cable, all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, Cable Structure bearing generalized coordinate data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C oUtilization can be expressed the Non-destructive Testing Data of the health status of support cable and set up cable system initial damage vector d o, cable system initial damage vector d oElement number equal N, d oElement and support cable be one-to-one relationship, 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 support cable, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, no problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value between 0 and 100%, then represents this support cable and lost the load-bearing capacity of corresponding proportion, if when not having the data of the Non-destructive Testing Data of support cable and other health status that can express support cable, perhaps think when the Cable 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, represent that then the corresponding support cable of this element is problematic, in the method this support cable may be impaired also may be lax, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, the initial equivalent damage degree of this its corresponding support cable of number of elements value representation; Cable system initial damage vector d oThe coding rule and initial rope force vector F of element oThe coding rule of element identical; Corresponding to A oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U oThe bearing generalized coordinate comprises two kinds of line amount and angle amounts;
D. according to the measured data of the design drawing of Cable Structure, as-constructed drawing, initial Cable Structure, the Non-destructive Testing Data of support cable, the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure, initial Cable Structure bearing generalized coordinate vector U o, initial Cable Structure steady temperature data vector T oWith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating 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 oCable Structure bearing generalized coordinate data be exactly initial Cable Structure bearing generalized coordinate vector U oCorresponding to A oThe support cable 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; T o, U oAnd d oA oParameter, by A oInitial 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 can say C oBy A oMechanics Calculation result form A in the method o, U o, C o, d oAnd T oConstant;
E. in the method, alphabetical i is except the place that represents significantly number of steps, and alphabetical i only represents cycle index, i.e. the i time circulation; The current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A during i time circulation beginning i o, A oAnd A i oCount temperature parameter, can accounting temperature have changed the Effect on Mechanical Properties to Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; The current initial Cable Structure bearing generalized coordinate vector that needs during the i time circulation beginning is designated as U i o, U i oThe current initial mechanical calculating benchmark model A of data representation Cable Structure i oCable Structure bearing generalized coordinate; The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d i o, d i oCable Structure A when representing this time circulation beginning i oThe health status of cable system, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i o, U i oAnd d i oA i oCharacterisitic parameter, C i oBy A i oMechanics Calculation result form; During for the first time circulation beginning, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring for the first time circulation beginning, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring for the first time circulation beginning, U i oBe designated as U 1 o, set up U 1 oMethod for making U 1 oEqual U oDuring for the first time circulation beginning, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring for the first time circulation beginning, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
F. go on foot the s circulation in step from entering by f here; In Cable Structure military service process, according to the current data of " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement acquisition Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector T i, vector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement corresponding one by one; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains Cable Structure bearing generalized coordinate current data, all Cable Structure bearing generalized coordinate current datas form current cable structures actual measurement bearing generalized coordinates vector U iObtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains the rope force data of all support cables in the Cable Structure, all these rope force datas form current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current support cable two supporting end points horizontal range vectors oThe coding rule of element identical; Obtain vector T in actual measurement iThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T iThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value form monitored amount current value vector C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly;
G. according to current cable structure actual measurement bearing generalized coordinate vector U iWith current cable structure steady temperature data vector T i, upgrade current initial mechanical calculating benchmark model A according to step g 1 to g3 i o, current initial Cable Structure bearing generalized coordinate vector U i o, the current initial value of monitored amount vector C i oWith current initial Cable Structure steady temperature data vector T i o, and the current initial damage vector of cable system d i oRemain unchanged;
G1. compare respectively U iWith U i o, T iWith T i oIf, U iEqual U i oAnd T iEqual T i o, A then i o, U i o, C i oAnd T i oRemain unchanged; Otherwise need to follow these steps to A i o, U i oAnd T i oUpgrade;
G2. calculate U iWith U oPoor, U iWith U oDifference be exactly the Cable Structure bearing about the current generalized displacement of support of initial position, represent generalized displacement of support with generalized displacement of support vector V, V equals U iDeduct U o, be one-to-one relationship between the element among the generalized displacement of support vector V and the generalized displacement of support component, the numerical value of an element is corresponding to the generalized displacement of an assigned direction of an appointment bearing among the generalized displacement of support vector V; Calculate T iWith T oPoor, T iWith T oDifference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T iWith T oPoor represent that with steady temperature change vector S S equals T iDeduct T o, S represents the variation of Cable Structure steady temperature data;
G3. first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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 oMiddle Cable Structure bearing applies generalized displacement of support constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial mechanical calculating benchmark model A that obtains upgrading i o, upgrade A i oThe time, U i oAll elements numerical value is also used U iAll elements numerical value is corresponding to be replaced, and has namely upgraded U i o, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, namely upgraded T i o, so just obtained correctly corresponding to A i oT i oAnd U i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status of rope with the current initial damage of cable system vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iExpression, A i oThe bearing generalized coordinate with current initial Cable Structure bearing generalized coordinate vector U i oC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A by Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value form C i o
H. at current initial mechanical calculating benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity transformation matrices Δ C by calculating iWith nominal unit damage vector D i u
H1. when the i time circulation beginning, directly press step h2 to the listed method acquisition of step h4 Δ C iAnd D i uAt other constantly, when in step g to A i oAfter upgrading, must regain Δ C to the listed method of step h4 by step h2 iAnd D i uIf, in step g not to A i oUpgrade, then directly change herein step I over to and carry out follow-up work;
H2. at current initial mechanical calculating benchmark model A i oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all support cables on the calculation times numerical value, there is N root support cable that N calculating is just arranged, calculating each time in the hypothesis cable system only has a support cable to increase unit damage on the basis of original damage again, the support cable that occurs damage in calculating each time is different from the support cable that occurs damage in other time calculating, and each time supposition has the unit damage value of the support cable of damage can be different from the unit damage value of other support cables, uses " nominal unit damage vector D i u" record the unit damage of the supposition of all ropes, vectorial D i uElement coding rule and vectorial d oThe coding rule of element identical, calculate each time the current value of all monitored amounts in the Cable Structure, the current value of the monitored amount of all that calculate each time forms one " monitored amount is calculated the current value vector "; When hypothesis j root support cable has unit damage, available C i Tj" monitored amount is calculated the current value vector " that expression is corresponding; When giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target; C i TjDefinition mode and vectorial C oDefinition mode identical, C i TjElement and C oElement corresponding one by one;
H3. the vectorial C that calculates each time i TjDeduct vectorial C i oObtain a vector, should obtain " monitored numberical value of quantity change vector δ a C after the unit damage value of each element divided by supposition in this calculating of vector again i j"; There is N root support cable that N " monitored numberical value of quantity change vector " just arranged;
H4. form successively " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged by this N " monitored numberical value of quantity change vector " i"; " unit damage monitored numerical quantity transformation matrices Δ C i" each row corresponding to one " monitored numberical value of quantity change vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d oThe element coding rule identical;
I. define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal the quantity of support cable, d i cAnd d iElement and support cable between be one-to-one relationship, d i cAnd d iElement numerical value represent degree of injury or the health status of corresponding support cable, d i cAnd d iWith cable system initial damage vector d oThe element coding rule identical, d i cElement, d iElement and d oElement be one-to-one relationship;
J. according to monitored amount current value vector C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, in the formula 1 except d i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + &Delta;C i &CenterDot; d c i Formula 1
K. the current actual damage vector d that utilizes formula 2 to express iJ element d i jWith the current initial damage vector of cable system d i oJ element d i OjWith the vectorial d of current name damage i cJ element d i CjBetween relation, calculate current actual damage vector d iAll elements;
d j i = 1 - ( 1 - oj i ) ( 1 - d cj i ) Formula 2
J=1 in the formula 2,2,3 ...., N, current actual damage vector d iJ element d i jNumerical value be to represent j root support cable without health problem, d at 0 o'clock i jNumerical value is not to represent that j root support cable was the support cable of unsoundness problem at 0 o'clock, and the support cable of unsoundness problem may be slack line, also may be damaged cable, its numerical response the degree of lax or damage; The current actual damage vector of cable system d iElement numerical value be not less than 0, be not more than 100%, the current actual damage of cable system vector d iElement numerical value represent the degree of injury of corresponding support cable, if the current actual damage vector of cable system d iThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, without health problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, represent that then the corresponding support cable of this element is the unsoundness problem, the health problem of this support cable may be impaired also may be to have relaxed in the method, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, this its corresponding support cable of number of elements value representation and current actual equivalent damage degree its relax level mechanics equivalence;
L. identify damaged cable from the problematic support cable that k identified the step, remaining is exactly slack line;
M. utilize at current cable structure steady temperature data vector T iThe current actual damage vector of the cable system d that obtains in the k step under the condition iObtain slack line and current actual equivalent damage degree its relax level mechanics equivalence, utilize the f step obtain at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current support cable two supporting end points horizontal ranges vectors, utilize the c step obtain at initial Cable Structure steady temperature data vector T oThe initial drift vector of the support cable under the condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum utilize current cable structure steady temperature data vector T iThe current steady temperature data of support cable of expression, utilize the c step obtain at initial Cable Structure steady temperature data vector T oThe support cable initial steady state temperature data of expression, utilization is in the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure of c step acquisition, count temperature variation to the impact of support cable physics, mechanics and geometric parameter, by with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with equivalent slack line, Cable Structure any variation can not occur, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li; Lax identification and the damage identification of support cable have so just been realized; Institute's demand power is provided by current cable force vector F corresponding element during calculating;
N. try to achieve the vectorial d of current name damage i cAfter, set up mark vector B according to formula 3 i, formula 4 has provided mark vector B iThe definition of j element;
B i = B 1 i B 2 i . . . B j i . . . B N i T Formula 3
Formula 4
Element B in the formula 4 i jMark vector B iJ element, D i UjNominal unit damage vector D i uJ element, d i CjThe vectorial d of the current name damage of cable system i cJ element, they all represent the relevant information of j root support cable, j=1 in the formula 4,2,3 ..., N;
If mark vector B o. iElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B iElement be not 0 entirely, then enter next step, be step p;
P. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required cable system of circulation d I+1 oEach element;
d oj i + 1 = 1 - ( 1 - d oj i ) ( 1 - D uj i B j i ) Formula 5
D in the formula 5 I+1 OjThe current initial damage vector of the required cable system d that next time, namely circulates for the i+1 time I+1 oJ element, d i OjThe current initial damage vector of the cable system d that is this, namely circulates for the i time i oJ element, D i UjThe nominal unit damage vector D of the i time circulation i uJ element, B i jThe mark vector B of the i time circulation iJ element, j=1 in the formula 5,2,3 ..., N;
Q. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual the current initial Cable Structure steady temperature data vector T of the i time circulation i o
R. at initial mechanical calculating benchmark model A oThe basis on, first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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, the health status that makes again rope is d I+1 oAfter obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation I+1Obtain A I+1After, obtain A by Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that namely circulates for the i+1 time I+1 oNext time, the i.e. required current initial Cable Structure bearing generalized coordinate vector U of the i+1 time circulation I+1 oEqual the current initial Cable Structure bearing generalized coordinate vector U of the i time circulation i o
S. get back to step f, beginning is circulation next time.
Beneficial effect: when being subjected to the affecting of the 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 affect the monitored amount of Cable Structure, only have and partly rejected by the impact in temperature field monitored amount could to carry out rational cable structure health monitoring based on monitored amount, and the temperature field measurement of Cable Structure and calculating are very complicated, this method discloses and has comprised a kind of the simple of cable structure health monitoring that be suitable for, economical, feasible, the cable structure health monitoring method of efficient Cable Structure Calculation Method of Temperature Field, adopt this method to occur in the situation of generalized displacement at the Cable Structure bearing, when many ropes of Cable Structure are synchronously impaired and lax, and when the temperature of Cable Structure changed along with the time, very the monitor assessment identifies the health status (position that comprises all slack lines and damaged cable of cable system, and relax level or degree of injury), the disclosed system and method for this method is very useful to effective health monitoring of cable system.
Embodiment
Generalized displacement of support is being arranged and when temperature variation is arranged, for 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 describes of the embodiment of this method in fact only is exemplary, and purpose never is to limit application or the use of this method.
This method adopts a kind of algorithm, and this algorithm is used for the health status of the cable system of monitoring identification Cable Structure.During implementation, the following step is a kind of in the various steps that can take.
The first step: determine " the temperature survey calculating method of the Cable Structure of this method " that the method concrete steps are as follows:
The a step: inquiry or actual measurement (can be measured with conventional thermometry, 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 0 the highest temperature and the lowest temperature between rear 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to the base area revolutions constantly and the rule that revolves round the sun is determined constantly, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between rear 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 that the daily temperature at T cloudy day 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 rBe no less than between inquiry Cable Structure location and place Altitude Region temperature that the meteorological data in recent years in 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and place Altitude Region about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; " R Cable Structure surface point " got on surface in Cable Structure, getting the Specific Principles of " R Cable Structure surface point " narrates in step b3, the back will obtain by actual observation record the temperature of this R Cable Structure surface point, the temperature data that claims actual measurement to obtain 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 by Calculation of Heat Transfer, just claim the temperature data that calculates to be " 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, uniform choosing is no less than three different sea level elevations on Cable Structure, if for example the sea level elevation of Cable Structure is between 0m to 200m, can choose height above sea level 0m so, 50m, 100m and height above sea level 200m, intersect with imaginary surface level and Cable Structure surface at each sea level elevation place that chooses, obtain intersection, the crossing cross surface that obtains of surface level and Cable Structure, intersection is the outer edge line of cross surface, choose 6 points at the intersection place on surface level and Cable Structure surface, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of 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 meteorological data throughout the year in zone, Cable Structure position and the physical dimension of Cable Structure, volume coordinate, the sunny slope of definite Cable Structure such as 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, two sections of these of intersection respectively have a mid point, cross these two mid points and get the outer normal of Cable Structure, 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, two intersection points are also just arranged, these two intersection points are divided into two line segments with intersection, get 2 points at two line segments respectively, totally 4 points, taken point is divided into each line segment in two line segments of intersection 3 sections of equal in length, get the outer normal on Cable Structure surface at these 4 some places, just chosen altogether like this outer normal on 6 Cable Structure surfaces at each selected sea level elevation place, the direction of 6 outer normals is exactly " measuring Cable Structure along the direction of the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures Cable Structure along the direction of 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 again three points, 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 that amounts to 5 points, concrete can hole in Cable Structure first, temperature sensor is embedded in this 5 some places, especially, can not hole in support cable, support cable is only measured the temperature of support cable surface point, in any case, the temperature that records all is called this place " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ".If chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of 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, especially, E equals 1 for support cable, that " measures Cable Structure along the point of the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, the back will obtain by actual measurement the temperature of this HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", the temperature data that claims actual measurement to obtain is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim the temperature data that calculates to be " HBE Cable Structure is along thickness temperature computation data "; Will be 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 " in this method ".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; On-site spacious without choosing a position in the place of blocking in Cable Structure, this position should can both obtain in each day of the whole year this ground the most sufficient sunshine of getable this day (as long as the same day sunrise arranged, this position just should be arrived by solar radiation), flat board (square that for example the wide 3mm of 30cm is thick is dull and stereotyped) at 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, reference plate overhead distance is not less than 1.5 meters, reference plate can place and meet the top that meteorology temperature is measured the wooden thermometer screen that requires, the one side of this reference plate on the sunny side, (for example be called sunny slope, in the time of on the Northern Hemisphere, sunny slope faces up towards south, full daytime is all by sunshine, the suitable gradient should be arranged sunny slope so that snow can not accumulate or clear up sunny slope behind snow), the sunny slope of reference plate is coarse and (being conducive to accept 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 the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material (for example thick calcium carbonate insulation material of 5mm), the Real-Time Monitoring record is obtained the temperature of the sunny slope of reference plate.
The b step, Real-Time Monitoring (can be measured with conventional thermometry, for example use thermal resistance to measure, for example every temperature data of 10 minutes survey records) record obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring (can be measured with conventional thermometry simultaneously, for example use thermal resistance to measure, for example every temperature data of 10 minutes survey records) obtain previously defined Cable Structure along the temperature profile data of thickness, Real-Time Monitoring (can be measured with conventional thermometry simultaneously, for example in meeting the wooden thermometer screen that meteorology temperature measure to require, lay thermal resistance and measure temperature, for example every 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 measure with conventional thermometry by 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 temperature data of 10 minutes survey records) record obtains being carved at sunrise the constantly temperature measured data sequence of the Cable Structure place environment between rear 30 minutes of sunrise next day the same day, the temperature measured data sequence of Cable Structure place environment is arranged sequentially according to time order and function by the temperature measured data that was carved at sunrise the Cable Structure place environment between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day the same day that minimum temperature obtains Cable Structure place environment, be designated as Δ T Emax(for example first the temperature measured data sequence of Cable Structure place environment is carried out curve fitting by conventional mathematical computations by the temperature measured data sequence of Cable Structure place environment, then by asking curve to the derivative of time or by asking on the curve each with numerical method corresponding to the point of the survey record data time rate of change to the time) obtain the temperature of Cable Structure place environment about the rate of change of time, this rate of change is also along with the time changes; (can measure with conventional thermometry by Real-Time Monitoring, for example use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, for example every temperature data of 10 minutes survey records) obtain being carved at sunrise the same day the constantly measured data sequence of the temperature of the sunny slope of the reference plate between rear 30 minutes of sunrise next day, the measured data sequence of the temperature of the sunny slope of reference plate is arranged sequentially according to time order and function by the measured data of the temperature of the sunny slope that was carved at sunrise the reference plate between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day same day of temperature that minimum temperature obtains the sunny slope of reference plate, be designated as Δ T Pmax(can measure with conventional thermometry by Real-Time Monitoring, for example use thermal resistance to measure the Cable Structure surface point, for example every temperature data of 10 minutes survey records) record obtains being carved at sunrise the constantly Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between rear 30 minutes of sunrise next day the same day, there be R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged, the Cable Structure surface temperature measured data between rear 30 minutes is arranged sequentially according to time order and function constantly by being carved at sunrise sunrise next day on same day of a Cable Structure surface point for each Cable Structure surface temperature measured data sequence, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day 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 at sunrise the constantly maximum temperature difference numerical value between rear 30 minutes of sunrise next day R the same day, maximal value wherein is designated as Δ T Smax(for example first each Cable Structure surface temperature measured data sequence is carried out curve fitting by conventional mathematical computations by each Cable Structure surface temperature measured data sequence, then by asking curve to the derivative of time or by asking on the curve each with numerical method corresponding to the point of the survey record data time rate of change to the time) obtain the temperature of each Cable Structure surface point about the rate of change of time, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes.Obtain being carved at sunrise the same day between rear 30 minutes of the moment of sunrise next day by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculating amounts to maximum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature 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 sea level elevation 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 " 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, determine to obtain the moment of Cable Structure steady temperature data, the condition relevant with the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between rear 30 minutes between the same day, sunset constantly refer on the meteorology that base area revolutions and revolution rule determine sunset constantly, the sunset that can inquire about data or calculate each required day by conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between rear 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 rear 30 minutes during this period of time in, measure in front the Δ T that calculates EmaxBe not more than with reference to temperature difference per day Δ T r, and measure in front the Δ T that calculates PmaxDeduct 2 degrees centigrade and be not more than Δ T Emax, and measure in front the Δ T that calculates SmaxBe not more than Δ T PmaxOnly needing to satisfy in second a condition and the b condition one just is called and satisfies second condition; The 3rd condition is that the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; 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 the rate of change of time; 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 rear 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 " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition of satisfying in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment of only satisfying the 6th condition in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during in constantly one of physical record data 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 be not any in constantly of physical record data in this method constantly, then get this method close to moment of mathematics those physical record data constantly that obtain Cable Structure steady temperature data 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, i.e. this Cable Structure temperature constantly temporal evolution not, 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 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 in the Cable Structure in the moment that obtains Cable Structure steady temperature data by conventional Calculation of Heat Transfer (for example finite element method), 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 in front accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " of Cable Structure, the accounting temperature of HBE " measuring Cable Structure along the point of 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, 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 " R Cable Structure surface point " got on the surface in Cable Structure, 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 the temperature of any point on the Cable Structure surface be by " 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 the support cable 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 uniform along the Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " along sea level elevation all in twos the maximal value Δ h in the absolute value of the difference of the sea level elevation of adjacent Cable Structure surface point be 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; " R Cable Structure surface point " along sea level elevation the definition of adjacent Cable Structure surface point refers to only consider sea level elevation in twos the time, do not have a Cable Structure surface point in " R 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 inquiry or obtains rule at sunshine between Cable Structure location and place Altitude Region by the meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, find in Cable Structure to be subjected to the sunshine-duration position of sufficient those surface points the whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is an annual point that is subjected in the most sufficient those surface points of sunshine-duration on the Cable Structure.
Second step: set up initial mechanical calculating benchmark model A o
If total N root support cable, the coding rule of at first definite rope, with rope numberings all in the Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step by this rule.Determine measured point (namely all characterize the specified point of Cable Structure strain information, are provided with K specified point), give all specified point numberings; Determine that the measured strain of each specified point (establishes the strain of L assigned direction measuring each specified point, do not require that each specified point has the strain of the designated direction of same number, here the just strain that establish L assigned direction measuring each specified point in order to narrate convenient), and to all measured strains number; Above-mentioned numbering equally will be for generating the vector sum matrix in subsequent step.Each specified point can be exactly a near point the fixed endpoint (for example being that the drag-line of cable-stayed bridge is at the stiff end on the bridge floor) of each root rope, and this point generally should not be stress concentration point, to avoid occurring excessive strain measurement value; This numbering equally will be for generating the vector sum matrix in subsequent step.In the strain that each specified point can only be measured a direction, the strain that also can measure a plurality of directions." the whole monitored strain data of Cable Structure " described by strain K specified point, that cross L assigned direction of each specified point on top definite Cable Structure, and the variation of Cable Structure strain is exactly the variation of the strain of all assigned directions all specified points, all appointment straight lines.(the individual strain measurement value of M=K * L) or calculated value characterize the strain information of Cable Structure to each total M.K and M must not be less than the quantity N of support cable.For simplicity, in the method with " the monitored strain data of Cable Structure " referred to as " monitored amount ".
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 measure with conventional thermometry, for example use thermal resistance to measure), this moment " Cable Structure steady temperature data " use vector T oExpression is called initial Cable Structure steady temperature data vector T oObtain T in actual measurement oThe time, namely 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 (such as initial or current etc.) Cable Structure steady temperature data vector, 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): at the survey record temperature (temperature that comprises Cable Structure place environment, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) time, for example every temperature of 10 minutes survey records, so simultaneously equally also every 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 determined to obtain the moment of Cable Structure steady temperature data, just be called synchronization in the moment that obtains Cable Structure steady temperature data with the data of 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 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 temperature variant physical parameter (for example thermal expansivity) and the mechanical property parameters (for example elastic modulus, Poisson ratio) of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is 0 o'clock length, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables when obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and initial rope force vector F of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum oThe coding rule of element identical; Obtain initial Cable Structure steady temperature data vector T at Actual measurement oThe time, namely obtaining initial Cable Structure steady temperature data vector T oThe synchronization in the moment, use the conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise that the Non-destructive Testing Data of support cable etc. can express the 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 bearing generalized coordinate data, Cable Structure modal data, Cable Structure strain data, Cable Structure angular coordinate measurement data, Cable Structure volume coordinate measurement data.Corresponding to A oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized 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 Cable Structure, and purpose is to determine according to these coordinate datas the geometric properties of Cable Structure.For 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 data of utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope are set up cable system initial damage vector d oIf when not having the data of the Non-destructive Testing Data of rope and other health status that can express support cable, perhaps can think when the Cable 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 and support cable be one-to-one relationship, 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 support cable, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, no problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value between 0 and 100%, then represents this support cable and lost the load-bearing capacity of corresponding proportion, if when not having the data of the Non-destructive Testing Data of support cable and other health status that can express support cable, perhaps think when the Cable 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, represent that then the corresponding support cable of this element is problematic, in the method this support cable may be impaired also may be lax, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, the initial equivalent damage degree of this its corresponding support cable of number of elements value representation; Cable system initial damage vector d oThe coding rule and initial rope force vector F of element oThe coding rule of element identical.Utilize the measured data of design drawing, as-constructed drawing and the initial Cable Structure of Cable Structure, the Non-destructive Testing Data of support cable, the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure, initial Cable Structure bearing generalized 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 mechanical calculating benchmark model A o
No matter which kind of method to obtain initial mechanical calculating 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%.But such utility A 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 occuring near institute's analog case reliably.Model A oThe health status of middle support cable cable system initial damage vector d oExpression, Cable Structure 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 oThe support cable 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 bearing generalized coordinate data with initial Cable Structure bearing generalized coordinate vector U oExpression; T o, U oAnd d oA oParameter, C oBy A oMechanics Calculation result form.
The 3rd step: in the method, alphabetical i is except the place that represents significantly number of steps, and alphabetical i only represents cycle index, i.e. the i time circulation; The current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A during i time circulation beginning i o, A oAnd A i oCount temperature parameter, can accounting temperature have changed the Effect on Mechanical Properties to Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; The current initial Cable Structure bearing generalized coordinate vector that needs during the i time circulation beginning is designated as U i o, U i oThe current initial mechanical calculating benchmark model A of data representation Cable Structure i oCable Structure bearing generalized coordinate.The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d i o, d i oCable Structure A when representing this time circulation beginning i oThe health status of cable system, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i, U i oAnd d i oA i oCharacterisitic parameter; C i oBy A i oMechanics Calculation result form; During for the first time circulation beginning, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring for the first time circulation beginning, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring for the first time circulation beginning, U i oBe designated as U 1 o, set up U 1 oMethod for making U 1 oEqual U oDuring for the first time circulation beginning, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring for the first time circulation beginning, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
The 4th step: the hardware components of pass line structural healthy monitoring system.Hardware components comprises at least: monitored amount monitoring system (for example contains strain measurement system, signal conditioner etc.), Cable Structure bearing generalized coordinate monitoring system (for example contains total powerstation, angular transducer, signal conditioner etc.), 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.), the cable force monitoring system (for example contains acceleration transducer, signal conditioner etc.), the horizontal range monitoring system (for example monitoring with total powerstation) of each support cable two supporting end points, signal (data) collector, computing machine and the panalarm of communicating by letter.Horizontal range, each temperature of the bearing generalized coordinate of each monitored amount, each Cable Structure, the Suo Li of each root support cable, each root support cable two supporting end points 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 damage arranged, the computer control communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 5th step: establishment and the laddering recognition methods system software of slack line of installation and operation generalized displacement of support temperature variation strain monitoring on computers, this software will be finished the functions such as monitoring that this method " the laddering recognition methods of slack line of generalized displacement of support temperature variation strain monitoring " required by task wants, record, control, storage, calculating, notice, warning (being all work that can finish with computing machine in this specific implementation method).
The 6th step: thus step beginning circular flow, in Cable Structure military service process, according to the current data of " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement acquisition Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector T i, vector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement correspondingly one by one obtain current cable structure steady temperature data vector T in actual measurement iThe time, actual measurement obtains Cable Structure bearing generalized coordinate current data, and all data form current cable structure actual measurement bearing generalized coordinate vector U iObtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains the rope force data of all support cables in the Cable Structure, all these rope force datas form current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current support cable two supporting end points horizontal range vectors oThe coding rule of element identical.In the actual measurement vector T iThe time, namely obtaining current cable structure steady temperature data vector T iThe synchronization in the moment, actual measurement obtains the currency of all monitored amounts in the Cable Structure, all these numerical value form monitored amount current value vector C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly.
The 7th step: obtaining current cable structure actual measurement bearing generalized coordinate vector U iWith current cable structure steady temperature data vector T iAfter, compare respectively U iAnd U i o, T iAnd T i oIf, U iEqual U i oAnd T iEqual T i o, then do not need A i o, U i oAnd T i oUpgrade, otherwise need to be to current initial mechanical calculating benchmark model A i o, current initial Cable Structure bearing generalized coordinate vector U i o, current initial Cable Structure steady temperature data vector T i oWith the current initial value vector of monitored amount C i oUpgrade, and the current initial damage vector of cable system d i oRemain unchanged, update method follows these steps to a and carries out to step c:
A. calculate U iWith U oPoor, U iWith U oDifference be exactly the Cable Structure bearing about the current generalized displacement of support of initial position, represent generalized displacement of support with generalized displacement of support vector V, V equals U iDeduct U o, be one-to-one relationship between the element among the generalized displacement of support vector V and the generalized displacement of support component, the numerical value of an element is corresponding to the generalized displacement of an assigned direction of an appointment bearing among the generalized displacement of support vector V.
B. calculate T iWith T oPoor, T iWith T oDifference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T iWith T oPoor represent that with steady temperature change vector S S equals T iDeduct T o, S represents the variation of Cable Structure steady temperature data.
C. first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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 oMiddle Cable Structure bearing applies generalized displacement of support constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial mechanical calculating benchmark model A that obtains upgrading i o, upgrade A i oThe time, U i oAll elements numerical value is also used U iAll elements numerical value is corresponding to be replaced, and has namely upgraded U i o, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, namely upgraded T i o, so just obtained correctly corresponding to A i oT i oAnd U i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status of rope with the current initial damage of cable system vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iExpression, A i oThe bearing generalized coordinate with current initial Cable Structure bearing generalized coordinate vector U i oC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A by Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value form C i o
The 8th step: at current initial mechanical calculating benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step a to steps d, set up unit damage monitored numerical quantity transformation matrices Δ C by calculating iWith nominal unit damage vector D i u
A. when the i time circulation beginning, directly press step b to the listed method acquisition of steps d Δ C iAnd D i uAt other constantly, when in the 7th step to A i oAfter upgrading, must regain Δ C to the listed method of steps d by step b iAnd D i uIf, the 7th the step in not to A i oUpgrade, then directly changing herein for the 9th step over to carries out follow-up work.
B. at current initial mechanical calculating benchmark model A i oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all support cables on the calculation times numerical value, there is N root support cable that N calculating is just arranged, calculate each time and suppose that only having a support cable to increase unit damage on the basis of original damage in the cable system (for example gets 5% again, 10%, 20% or 30% equivalent damage is unit damage), the support cable that occurs damage in calculating each time is different from the support cable that occurs damage in other time calculating, and each time supposition has the unit damage value of the support cable of damage can be different from the unit damage value of other support cables, uses " nominal unit damage vector D i u" record the unit damage of the supposition of all ropes, vectorial D i uElement coding rule and vectorial d oThe coding rule of element identical, calculate each time the current value of all monitored amounts in the Cable Structure, the current value of the monitored amount of all that calculate each time forms one " monitored amount is calculated the current value vector "; When the hypothesis j (j=1,2,3 ..., when N) the root support cable has unit damage, available C i Tj" monitored amount is calculated the current value vector " that expression is corresponding; When giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target; C i TjDefinition mode and vectorial C oDefinition mode identical, C i TjElement and C oElement corresponding one by one.
C. the vectorial C that calculates each time i TjDeduct vectorial C i oObtain a vector, should obtain " monitored numberical value of quantity change vector δ a C after the unit damage value of each element divided by supposition in this calculating of vector again i j"; There is N root support cable that N " monitored numberical value of quantity change vector " just arranged.
D. form successively " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged by this N " monitored numberical value of quantity change vector " i"; " unit damage monitored numerical quantity transformation matrices Δ C i" each row corresponding to one " monitored numberical value of quantity change vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d oThe element coding rule identical.
The 9th step: set up linear relationship error vector e iWith vectorial g iUtilize data (" the current initial value vector of the monitored amount C of front i o", " unit damage monitored numerical quantity transformation matrices Δ C i"); when the 8th step calculated each time; namely in calculating each time, only have in the hypothesis cable system increase unit damage again on the basis of rope in original damage in; when hypothesis j (j=1,2,3; ..., when N) the root support cable has unit damage, calculate each time damage vector of composition, use d i TjRepresent this damage vector, it is C that corresponding monitored amount is calculated the current value vector i Tj(referring to the 8th step), damage vectorial d i TjElement number equal the quantity of rope, vectorial d i TjAll elements in only have the numerical value of an element to get to calculate each time in hypothesis increase the unit damage value of the rope of unit damage, d i TjThe numerical value of other element get 0, that is not numbering and the supposition of 0 the element corresponding relation that increases the rope of unit damage, be identical with the element of the same numbering of other vectors with the corresponding relation of this rope; d i TjWith cable system initial damage vector d oThe element coding rule identical, d i TjElement and d oElement be one-to-one relationship.With C i Tj, C i o, Δ C i, d i TjBring formula (23) into, obtain a linear relationship error vector e i j, calculate each time a linear relationship error vector e i je i jSubscript j represent j (j=1,2,3 ..., N) the root support cable has unit damage.There is N root rope that N calculating is just arranged, N linear relationship error vector e just arranged i j, with this N linear relationship error vector e i The j phaseObtain a vector after adding, the new vector that each element of this vector is obtained after divided by N is exactly final linear relationship error vector e iVector g iEqual final error vector e iWith vectorial g iBe kept on the hard disc of computer of operation health monitoring systems software, for the health monitoring systems software application.
e j i = abs ( &Delta;C i &CenterDot; d tj i - C tj i + C o i ) - - - ( 23 )
The tenth step: define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal the quantity of support cable, d i cAnd d iElement and support cable between be one-to-one relationship, d i cAnd d iElement numerical value represent degree of injury or the health status of corresponding support cable, d i cAnd d iWith cable system initial damage vector d oThe element coding rule identical, d i cElement, d iElement and d oElement be one-to-one relationship.
The 11 step: according to monitored amount current value vector C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula (11), calculates the vectorial d of current name damage according to multi-objective optimization algorithm i cNoninferior solution, namely with reasonable error but can be more exactly from all ropes, determine the position of damaged cable and the solution of nominal degree of injury thereof.
The multi-objective optimization algorithm that can adopt has a variety of, for example: based on the multiple-objection optimization of genetic algorithm, based on the multiple-objection optimization of artificial neural network, based on the multi-objective optimization algorithm of population, multiple-objection optimization, leash law (Constrain Method), weighted method (Weighted SUm Method), Objective Programming (Goal Attainment Method) etc. based on ant group algorithm.Because various multi-objective optimization algorithms all are conventional algorithms, can realize easily that this implementation step is only found the solution the vectorial d of current name damage take Objective Programming as example provides i cProcess, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.
According to Objective Programming, formula (11) can transform the multi-objective optimization question shown in an accepted way of doing sth (24) and the formula (25), and γ is a real number in the formula (24), and R is real number field, and area of space Ω has limited vectorial d i cSpan (the present embodiment requirements vector d of each element i cEach element be not less than 0, be not more than 1).The meaning of formula (24) is to seek the real number γ of a minimum, so that formula (25) is met.G (d in the formula (25) i c) by formula (25) definition, G (d in the product representation formula (25) of weighing vector W and γ in the formula (25) i c) and vectorial g iBetween the deviation that allows, g iDefinition referring to formula (17), its value the 9th the step calculate.During actual computation the vector W can with vectorial g iIdentical.The concrete programming of Objective Programming realizes having had universal program directly to adopt.Use Objective Programming just can damage vectorial d in the hope of current name i c
min mize &gamma; &gamma; &Element; R , d c i &Element; &Omega; - - - ( 24 )
G ( d c i ) - W&gamma; &le; g i - - - ( 25 )
G ( d c i ) = abs ( &Delta;C i &CenterDot; d c i - C i + C o i ) - - - ( 26 )
The 12 step: according to the current actual damage vector of cable system d iDefinition (seeing formula (18)) and the definition (seeing formula (19)) of its element calculate current actual damage vector d iEach element, thereby can be by d iDefine position and the degree of injury of the support cable of health problem.d i j(i=1,2,3, J=1,2,3 ...., the N) health status of j root rope in the i time circulation of expression, formula (19), d are seen in its definition i jBe to represent j root support cable without health problem, d at 0 o'clock i jNumerical value is not to represent that j root support cable was the support cable of unsoundness problem at 0 o'clock, and the support cable of unsoundness problem may be slack line, also may be damaged cable, its numerical response the degree of lax or damage; The current actual damage vector of cable system d iElement numerical value be not less than 0, be not more than 100%, the current actual damage of cable system vector d iElement numerical value represent the degree of injury of corresponding support cable, if the current actual damage vector of cable system d iThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, without health problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, represent that then the corresponding support cable of this element is the unsoundness problem, the health problem of this support cable may be impaired also may be to have relaxed in the method, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, this its corresponding support cable of number of elements value representation and current actual equivalent damage degree its relax level mechanics equivalence.
The 13 step: identify damaged cable in the problematic support cable that identifies from the 12 step, remaining is exactly slack line.The method of differentiating is varied; can pass through the protective seam of the support cable of removal unsoundness problem; support cable is carried out visual discriminating; perhaps carry out visual discriminating by optical imaging apparatus; also can be by lossless detection method to support cable impaired discriminating the whether, UT (Ultrasonic Testing) is exactly a kind of now widely used lossless detection method.Those support cables of not finding the unsoundness problem damaged are exactly that lax rope has occured after differentiating, need adjust exactly the rope of Suo Li.
The 14 step: utilize at current cable structure steady temperature data vector T iThe current actual damage vector of the cable system d that obtains in the 12 step under the condition iObtain slack line and current actual equivalent damage degree its relax level mechanics equivalence, utilize the 6th step obtain at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current support cable two supporting end points horizontal ranges vectors, utilize second step obtain at initial Cable Structure steady temperature data vector T oThe initial drift vector of the support cable under the condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum utilize current cable structure steady temperature data vector T tThe current steady temperature data of support cable of expression, utilize second step obtain at initial Cable Structure steady temperature data vector T oThe support cable initial steady state temperature data of expression, utilization is in the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure of second step acquisition, count temperature variation to the impact of support cable physics, mechanics and geometric parameter, by with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical.When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with equivalent slack line, Cable Structure any variation can not occur, vice versa.Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li.Particularly can be in the hope of the relax level (being the long adjustment amount of rope) of these ropes according to formula (29) or formula (30).Lax identification and the damage identification of support cable have so just been realized.Institute's demand power is provided by current cable force vector F corresponding element during calculating.
The 15 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.Under specified requirements, computing machine automatic operation in health monitoring systems communication panalarm to monitor staff, owner and (or) personnel of appointment report to the police.
The 16 step: set up mark vector B according to formula (31) i, formula (32) has provided mark vector B iThe definition of j element; If mark vector B iElement be 0 entirely, then got back to for the 6th step and proceed health monitoring and calculating to cable system; If mark vector B iElement be not 0 entirely, then finish subsequent step after, enter next time circulation.
The 17 the step: at first according to formula (33) calculate next time (namely the i+1 time, i=1,2,3,4 ...) the required initial damage vector d of circulation I+1 oEach element d I+1 Oj(j=1,2,3 ..., N); The second, at initial mechanical calculating benchmark model A oThe basis on, first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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, the health status that makes again rope is d I+1 oAfter obtain be exactly next time, namely the i+1 time (i=1,2,3,4 ...) the required Mechanics Calculation benchmark model A of circulation I+1Next time (namely the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual T i o, next time (namely the i+1 time, i=1,2,3,4 ...) the required current initial Cable Structure bearing generalized coordinate vector U of circulation I+1 oEqual U i oObtain A I+1, U I+1 o, d I+1 oAnd T I+1 oAfter, obtain A by Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that namely circulates for the i+1 time I+1 o
The 18 step: got back to for the 6th step, beginning is by the circulation in the 6th step to the 18 step.

Claims (1)

1. the laddering recognition methods of slack line of a generalized displacement of support temperature variation strain monitoring is characterized in that described method comprises:
A. establish total N root support cable, at first determine the coding rule of support cable, with support cable numberings all in the Cable Structure, this numbering will be for generating the vector sum matrix in subsequent step by this rule; Determine the monitored point of appointment, monitored point namely characterizes all specified points of Cable Structure strain information, and gives all specified point numberings; Determine monitored should the changing direction of monitored point, and give the monitored strain numbering of all appointments; " monitored strain numbering " will be for generating the vector sum matrix in subsequent step; " the whole monitored strain data of Cable Structure " is comprised of above-mentioned all monitored strains; This method with " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity of monitored point must not be less than the quantity of support cable; The quantity sum of all monitored amounts must not be less than the quantity of support cable; Must not be greater than 30 minutes to the time interval between any twice measurement of same amount Real-Time Monitoring in this method, 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 " is undertaken by step b1 to b3;
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, in the method can not be seen one of the sun daytime and be called all day the cloudy day, statistics obtains 0 the highest temperature and the lowest temperature between rear 30 minutes of the moment of sunrise next day at each cloudy day in T cloudy day, sunrise on the meteorology that sunrise refers to the base area revolutions constantly and the rule that revolves round the sun is determined constantly, do not represent necessarily can see the sun same day, the sunrise that can inquire about data or calculate each required day by conventional meteorology constantly, each cloudy day 0 up to next day sunrise constantly the highest temperature between rear 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 that the daily temperature at T cloudy day 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 rBe no less than between inquiry Cable Structure location and place Altitude Region temperature that the meteorological data in recent years in 2 years or actual measurement obtain Cable Structure environment of living in time with delta data and the Changing Pattern of sea level elevation, calculate the temperature of the Cable Structure environment of living in recent years that is no less than 2 years between Cable Structure location and place Altitude Region about the maximum rate of change Δ T of sea level elevation h, for Δ T is got in convenient narration hUnit be ℃/m; " R Cable Structure surface point " got on surface in Cable Structure, getting the Specific Principles of " R Cable Structure surface point " narrates in step b3, the back will obtain by actual measurement the temperature of this R Cable Structure surface point, the temperature data that claims actual measurement to obtain 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 by Calculation of Heat Transfer, just claim the temperature data that calculates to be " 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, uniform choosing is no less than three different sea level elevations on Cable Structure, at each sea level elevation place that chooses, at least choose two points at the intersection place on surface level and Cable Structure surface, outer normal from selected point straw line body structure surface, all outer normal directions of choosing are called " measuring Cable Structure along the direction of 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 the sunny slope outer normal direction that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness and Cable Structure, direction uniform choosing in Cable Structure along each measurement Cable Structure along the Temperature Distribution of wall thickness is no less than three points, especially, measure Cable Structure for support cable along each and only get a point along the direction of the Temperature Distribution of wall thickness, namely only measure the temperature of the surface point of support cable, measure all and be selected a little temperature, the temperature that records is called " Cable Structure is along the temperature profile data of thickness ", wherein along crossing with same " intersection on surface level and Cable Structure surface ", " measure Cable Structure along the direction of the Temperature Distribution of wall thickness " and measure " Cable Structure is along the temperature profile data of thickness " that obtain, be called in the method " identical sea level elevation Cable Structure is along the temperature profile data of thickness ", if chosen H different sea level elevation, at each sea level elevation place, choose B and measured Cable Structure along the direction of 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, especially, E equals 1 for support cable, that " measures Cable Structure along the point of the temperature profile data of thickness " on the meter Cable Structure adds up to HBE, the back will obtain by actual measurement the temperature of this HBE " measuring Cable Structure along the point of the temperature profile data of thickness ", the temperature data that claims actual measurement to obtain is " HBE Cable Structure is along thickness temperature measured data ", if utilize the thermal conduction study computation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure Cable Structure along the temperature of the point of the temperature profile data of thickness, just claim the temperature data that calculates to be " HBE Cable Structure is along thickness temperature computation data "; Will be 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 " in this method "; 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; On-site spacious without choosing a position in the place of blocking in Cable Structure, this position should can both obtain in each day of the whole year this ground the most sufficient sunshine of getable this day, flat board at a carbon steel material of this position of sound production, be called reference plate, reference plate can not contact with ground, reference plate overhead distance is not less than 1.5 meters, 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 the most sufficient sunshine of getable this day, the non-sunny slope of reference plate is covered with insulation material, Real-Time Monitoring is obtained the temperature of the sunny slope of reference plate;
B2: Real-Time Monitoring obtains R Cable Structure surface temperature measured data of above-mentioned R Cable Structure surface point, Real-Time Monitoring obtains previously defined Cable Structure along the temperature profile data of thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology is measured the Cable Structure place environment of temperature requirement simultaneously; Obtain being carved at sunrise the same day the constantly temperature measured data sequence of the Cable Structure place environment between rear 30 minutes of sunrise next day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged sequentially according to time order and function by the temperature measured data that was carved at sunrise the Cable Structure place environment between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day on same day that minimum temperature obtains Cable Structure place environment, be called the environment maximum temperature difference, be designated as Δ T EmaxObtain the temperature of Cable Structure place environment about the rate of change of time by the temperature measured data sequence of Cable Structure place environment by conventional mathematical computations, this rate of change is also along with the time changes; Obtain being carved at sunrise the same day the constantly measured data sequence of the temperature of the sunny slope of the reference plate between rear 30 minutes of sunrise next day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged sequentially according to time order and function by the measured data of the temperature of the sunny slope that was carved at sunrise the reference plate between rear 30 minutes of the moment of sunrise next day the same day, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day 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 the constantly Cable Structure surface temperature measured data sequence of all R Cable Structure surface points between rear 30 minutes of sunrise next day by Real-Time Monitoring, there be R Cable Structure surface point that R Cable Structure surface temperature measured data sequence just arranged, the Cable Structure surface temperature measured data between rear 30 minutes is arranged sequentially according to time order and function constantly by being carved at sunrise sunrise next day on same day of a Cable Structure surface point for each Cable Structure surface temperature measured data sequence, 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 at sunrise the constantly maximum temperature difference between rear 30 minutes of sunrise next day on same day that minimum temperature obtains the temperature of each Cable Structure surface point, there be R Cable Structure surface point just to have and be carved at sunrise the constantly maximum temperature difference numerical value between rear 30 minutes of sunrise next day R the same day, maximal value wherein is called Cable Structure surface maximum temperature difference, is designated as Δ T SmaxObtain the temperature of each Cable Structure surface point about the rate of change of time by each Cable Structure surface temperature measured data sequence by conventional mathematical computations, the temperature of each Cable Structure surface point about the rate of change of time also along with the time changes; Obtain being carved at sunrise the same day between rear 30 minutes of the moment of sunrise next day by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculating amounts to maximum temperature among the BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and the difference of minimum temperature 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 sea level elevation 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 " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T Tmax
B3: measure and calculate acquisition Cable Structure steady temperature data; At first, determine to obtain the moment of Cable Structure steady temperature data, the condition relevant with the moment that determines acquisition Cable Structure steady temperature data has six, first condition is the moment that obtains Cable Structure steady temperature data to be carved at sunset sunrise next day constantly between rear 30 minutes between the same day, sunset constantly refer on the meteorology that base area revolutions and revolution rule determine sunset constantly, the sunset that can inquire about data or calculate each required day by conventional meteorology is constantly; The a condition of second condition be the same day be carved at sunrise next day sunrise constantly between rear 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 rear 30 minutes during this period of time in, measure in front the environment maximum error Δ T that calculates 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 PmaxOnly needing to satisfy in second a condition and the b condition one just is called and satisfies second condition; The 3rd condition is that the temperature of Cable Structure place environment is not more than per hour 0.1 degree centigrade about the absolute value of the rate of change of time in the moment that obtains Cable Structure steady temperature data; 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 the rate of change of time; 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 rear 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 " obtaining the mathematics of Cable Structure steady temperature data constantly ", the first is first moment to the 5th condition of satisfying in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, the second is the moment of only satisfying the 6th condition in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly, and the third is first moment to the 6th condition of satisfying simultaneously in above-mentioned " condition relevant with the moment that determines acquisition Cable Structure steady temperature data " constantly; When the mathematics that obtains Cable Structure steady temperature data is exactly in this method during in constantly one of physical record data 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 be not any in constantly of physical record data in this method constantly, then get this method close to moment of mathematics those physical record data constantly that obtain Cable Structure steady temperature data 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, i.e. this Cable Structure temperature constantly temporal evolution not, 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 in the Cable Structure in the moment that obtains Cable Structure steady temperature data by conventional Calculation of Heat Transfer, 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 in front accounting temperature of selected HBE " measuring Cable Structure along the point of the temperature profile data of thickness " of Cable Structure, the accounting temperature of HBE " measuring Cable Structure along the point of 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, 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 " R Cable Structure surface point " got on the surface in Cable Structure, 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 the temperature of any point on the Cable Structure surface be by " 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 the support cable 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 uniform along the Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " along sea level elevation all in twos the maximal value Δ h in the absolute value of the difference of the sea level elevation of adjacent Cable Structure surface point be 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; " R Cable Structure surface point " along sea level elevation the definition of adjacent Cable Structure surface point refers to only consider sea level elevation in twos the time, do not have a Cable Structure surface point in " R 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 inquiry or obtains rule at sunshine between Cable Structure location and place Altitude Region by the meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, find in Cable Structure to be subjected to the sunshine-duration position of sufficient those surface points the whole year, having a Cable Structure surface point in " R Cable Structure surface point " at least is an annual point that is subjected in the most sufficient those surface points of sunshine-duration on the Cable Structure;
C. directly measure the Cable Structure steady temperature data that calculate under the original state according to " the temperature survey calculating method of the Cable Structure of this method ", 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 physical and mechanical properties parameter of the employed various materials of Cable Structure; Obtain initial Cable Structure steady temperature data vector T in actual measurement oSynchronization, directly measure the initial Suo Li that calculates all support cables, form initial rope force vector F oAccording to Cable Structure design data, completion data obtain that all support cables are in free state that Suo Li is 0 o'clock length, the weight of cross-sectional area during in free state and the unit length during in free state, and the temperature of all support cables when obtaining these three kinds of data, utilize on this basis temperature variant physical function parameter and the mechanical property parameters of all support cables, obtain all support cables at initial Cable Structure steady temperature data vector T according to conventional physical computing oSuo Li under the condition is that the length of 0 o'clock all support cable, cross-sectional area and the Suo Li that Suo Li is 0 o'clock all support cable are the weight of the unit length of 0 o'clock all support cable, form successively the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the coding rule and initial rope force vector F of the element of the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum oThe coding rule of element identical; Obtain T in actual measurement oThe time, namely 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 initial value of the Non-destructive Testing Data of the health status of expressing support cable, all monitored amounts, the initial rope force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, Cable Structure bearing generalized coordinate data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C oUtilization can be expressed the Non-destructive Testing Data of the health status of support cable and set up cable system initial damage vector d o, cable system initial damage vector d oElement number equal N, d oElement and support cable be one-to-one relationship, 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 support cable, if cable system initial damage vector d oThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, no problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value between 0 and 100%, then represents this support cable and lost the load-bearing capacity of corresponding proportion, if when not having the data of the Non-destructive Testing Data of support cable and other health status that can express support cable, perhaps think when the Cable 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, represent that then the corresponding support cable of this element is problematic, in the method this support cable may be impaired also may be lax, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, the initial equivalent damage degree of this its corresponding support cable of number of elements value representation; Cable system initial damage vector d oThe coding rule and initial rope force vector F of element oThe coding rule of element identical; Corresponding to A oCable Structure bearing generalized coordinate data form initial Cable Structure bearing generalized coordinate vector U oThe bearing generalized coordinate comprises two kinds of line amount and angle amounts;
D. according to the measured data of the design drawing of Cable Structure, as-constructed drawing, initial Cable Structure, the Non-destructive Testing Data of support cable, the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure, initial Cable Structure bearing generalized coordinate vector U o, initial Cable Structure steady temperature data vector T oWith all Cable Structure data that preceding step obtains, set up the initial mechanical calculating 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 oCable Structure bearing generalized coordinate data be exactly initial Cable Structure bearing generalized coordinate vector U oCorresponding to A oThe support cable 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; T o, U oAnd d oA oParameter, by A oInitial 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 can say C oBy A oMechanics Calculation result form A in the method o, U o, C o, d oAnd T oConstant;
E. in the method, alphabetical i is except the place that represents significantly number of steps, and alphabetical i only represents cycle index, i.e. the i time circulation; The current initial mechanical calculating benchmark model of Cable Structure that need to set up or that set up is designated as current initial mechanical calculating benchmark model A during i time circulation beginning i o, A oAnd A i oCount temperature parameter, can accounting temperature have changed the Effect on Mechanical Properties to Cable Structure; During the i time circulation beginning, corresponding to A i o" Cable Structure steady temperature data " with current initial Cable Structure steady temperature data vector T i oExpression, vector T i oDefinition mode and vector T oDefinition mode identical, T i oElement and T oElement corresponding one by one; The current initial Cable Structure bearing generalized coordinate vector that needs during the i time circulation beginning is designated as U i o, U i oThe current initial mechanical calculating benchmark model A of data representation Cable Structure i oCable Structure bearing generalized coordinate; The current initial damage vector of cable system that needs during the i time circulation beginning is designated as d i o, d i oCable Structure A when representing this time circulation beginning i oThe health status of cable system, d i oDefinition mode and d oDefinition mode identical, d i oElement and d oElement corresponding one by one; During the i time circulation beginning, the initial value of all monitored amounts is with the current initial value vector of monitored amount C i oExpression, vectorial C i oDefinition mode and vectorial C oDefinition mode identical, C i oElement and C oElement corresponding one by one, the current initial value vector of monitored amount C i oExpression is corresponding to A i oThe concrete numerical value of all monitored amounts; T i o, U i oAnd d i oA i oCharacterisitic parameter, C i oBy A i oMechanics Calculation result form; During for the first time circulation beginning, A i oBe designated as A 1 o, set up A 1 oMethod for making A 1 oEqual A oDuring for the first time circulation beginning, T i oBe designated as T 1 o, set up T 1 oMethod for making T 1 oEqual T oDuring for the first time circulation beginning, U i oBe designated as U 1 o, set up U 1 oMethod for making U 1 oEqual U oDuring for the first time circulation beginning, d i oBe designated as d 1 o, set up d 1 oMethod for making d 1 oEqual d oDuring for the first time circulation beginning, C i oBe designated as C 1 o, set up C 1 oMethod for making C 1 oEqual C o
F. go on foot the s circulation in step from entering by f here; In Cable Structure military service process, according to the current data of " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement acquisition Cable Structure steady temperature data, the current data of all " Cable Structure steady temperature data " forms current cable structure steady temperature data vector T iVector T iDefinition mode and vector T oDefinition mode identical, T iElement and T oElement corresponding one by one; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains Cable Structure bearing generalized coordinate current data, all Cable Structure bearing generalized coordinate current datas form current cable structures actual measurement bearing generalized coordinates vector U iObtain current cable structure steady temperature data vector T in actual measurement iSynchronization, actual measurement obtains the rope force data of all support cables in the Cable Structure, all these rope force datas form current cable force vector F, the element of vectorial F and vectorial F oThe coding rule of element identical; Obtain current cable structure steady temperature data vector T in actual measurement iSynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the volume coordinate of two the supporting end points in the horizontal direction difference of component is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two supporting end points horizontal range vectors, the coding rule and initial rope force vector F of the element of current support cable two supporting end points horizontal range vectors oThe coding rule of element identical; Obtain vector T in actual measurement iThe time, actual measurement obtains obtaining current cable structure steady temperature data vector T iThe Cable Structure of synchronization in the moment in the currency of all monitored amounts, all these numerical value form monitored amount current value vector C i, vectorial C iDefinition mode and vectorial C oDefinition mode identical, C iElement and C oElement corresponding one by one, represent that identical monitored amount is at difference numerical value constantly;
G. according to current cable structure actual measurement bearing generalized coordinate vector U iWith current cable structure steady temperature data vector T i, upgrade current initial mechanical calculating benchmark model A according to step g 1 to g3 i o, current initial Cable Structure bearing generalized coordinate vector U i o, the current initial value of monitored amount vector C i oWith current initial Cable Structure steady temperature data vector T i o, and the current initial damage vector of cable system d i oRemain unchanged;
G1. compare respectively U iWith U i o, T iWith T i oIf, U iEqual U i oAnd T iEqual T i o, A then i o, U i o, C i oAnd T i oRemain unchanged; Otherwise need to follow these steps to A i o, U i oAnd T i oUpgrade;
G2. calculate U iWith U oPoor, U iWith U oDifference be exactly the Cable Structure bearing about the current generalized displacement of support of initial position, represent generalized displacement of support with generalized displacement of support vector V, V equals U iDeduct U o, be one-to-one relationship between the element among the generalized displacement of support vector V and the generalized displacement of support component, the numerical value of an element is corresponding to the generalized displacement of an assigned direction of an appointment bearing among the generalized displacement of support vector V; Calculate T iWith T oPoor, T iWith T oDifference be exactly that current cable structure steady temperature data are about the variation of initial Cable Structure steady temperature data, T iWith T oPoor represent that with steady temperature change vector S S equals T iDeduct T o, S represents the variation of Cable Structure steady temperature data;
G3. first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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 oMiddle Cable Structure bearing applies generalized displacement of support constraint and to A oIn the temperature variation that applies of Cable Structure after the current initial mechanical calculating benchmark model A that obtains upgrading i o, upgrade A i oThe time, U i oAll elements numerical value is also used U iAll elements numerical value is corresponding to be replaced, and has namely upgraded U i o, T i oAll elements numerical value is also used T iCorresponding replacement of all elements numerical value, namely upgraded T i o, so just obtained correctly corresponding to A i oT i oAnd U i oThis moment d i oRemain unchanged; When upgrading A i oAfter, A i oThe health status of rope with the current initial damage of cable system vector d i oExpression, A i oThe Cable Structure steady temperature with current cable structure steady temperature data vector T iExpression, A i oThe bearing generalized coordinate with current initial Cable Structure bearing generalized coordinate vector U i oC is upgraded in expression i oMethod be: when upgrading A i oAfter, obtain A by Mechanics Calculation i oIn concrete numerical value all monitored amounts, current, these concrete numerical value form C i o
H. at current initial mechanical calculating benchmark model A i oThe basis on, carry out the several times Mechanics Calculation according to step h1 to step h4, set up unit damage monitored numerical quantity transformation matrices Δ C by calculating iWith nominal unit damage vector D i u
H1. when the i time circulation beginning, directly press step h2 to the listed method acquisition of step h4 Δ C iAnd D i uAt other constantly, when in step g to A i oAfter upgrading, must regain Δ C to the listed method of step h4 by step h2 iAnd D i uIf, in step g not to A i oUpgrade, then directly change herein step I over to and carry out follow-up work;
H2. at current initial mechanical calculating benchmark model A i oThe basis on carry out the several times Mechanics Calculation, equal the quantity of all support cables on the calculation times numerical value, there is N root support cable that N calculating is just arranged, calculating each time in the hypothesis cable system only has a support cable to increase unit damage on the basis of original damage again, the support cable that occurs damage in calculating each time is different from the support cable that occurs damage in other time calculating, and each time supposition has the unit damage value of the support cable of damage can be different from the unit damage value of other support cables, uses " nominal unit damage vector D i u" record the unit damage of the supposition of all ropes, vectorial D i uElement coding rule and vectorial d oThe coding rule of element identical, calculate each time the current value of all monitored amounts in the Cable Structure, the current value of the monitored amount of all that calculate each time forms one " monitored amount is calculated the current value vector "; When hypothesis j root support cable has unit damage, available C i Tj" monitored amount is calculated the current value vector " that expression is corresponding; When giving in this step each vectorial element numbering, should use same coding rule with other vector in this method, to guarantee any one element in each vector in this step, with element in other vector, that numbering is identical, expressed the relevant information of same monitored amount or same target; C i TjDefinition mode and vectorial C oDefinition mode identical, C i TjElement and C oElement corresponding one by one;
H3. the vectorial C that calculates each time i TjDeduct vectorial C i oObtain a vector, should obtain " monitored numberical value of quantity change vector δ a C after the unit damage value of each element divided by supposition in this calculating of vector again i j"; There is N root support cable that N " monitored numberical value of quantity change vector " just arranged;
H4. form successively " the unit damage monitored numerical quantity transformation matrices Δ C that the N row are arranged by this N " monitored numberical value of quantity change vector " i"; " unit damage monitored numerical quantity transformation matrices Δ C i" each row corresponding to one " monitored numberical value of quantity change vector "; The coding rule of the row of " unit damage monitored numerical quantity transformation matrices " and cable system initial damage vector d oThe element coding rule identical;
I. define the vectorial d of current name damage i cWith current actual damage vector d i, d i cAnd d iElement number equal the quantity of support cable, d i cAnd d iElement and support cable between be one-to-one relationship, d i cAnd d iElement numerical value represent degree of injury or the health status of corresponding support cable, d i cAnd d iIdentical with the element coding rule of cable system initial damage vector do, d i cElement, d iElement and d oElement be one-to-one relationship;
J. according to monitored amount current value vector C iWith " the current initial value vector of monitored amount C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " the vectorial d of current name damage i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, in the formula 1 except d i cOther outer amount is known, finds the solution formula 1 and just can calculate the vectorial d of current name damage i c
C i = C o i + &Delta;C i &CenterDot; d c i Formula 1
K. the current actual damage vector d that utilizes formula 2 to express iJ element d i jWith the current initial damage vector of cable system d i oJ element d i OjWith the vectorial d of current name damage i cJ element d i CjBetween relation, calculate current actual damage vector d iAll elements;
d j i = 1 - ( 1 - oj i ) ( 1 - d cj i ) Formula 2
J=1 in the formula 2,2,3 ...., N, current actual damage vector d iJ element d i jNumerical value be to represent j root support cable without health problem, d at 0 o'clock i jNumerical value is not to represent that j root support cable was the support cable of unsoundness problem at 0 o'clock, and the support cable of unsoundness problem may be slack line, also may be damaged cable, its numerical response the degree of lax or damage; The current actual damage vector of cable system d iElement numerical value be not less than 0, be not more than 100%, the current actual damage of cable system vector d iElement numerical value represent the degree of injury of corresponding support cable, if the current actual damage vector of cable system d iThe numerical value of a certain element be 0, represent that the corresponding support cable of this element is intact, without health problem, if its numerical value is 100%, represent that then the corresponding support cable of this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, represent that then the corresponding support cable of this element is the unsoundness problem, the health problem of this support cable may be impaired also may be to have relaxed in the method, when this support cable when being impaired, the degree of injury of this its corresponding support cable of number of elements value representation, when if this support cable is lax, this its corresponding support cable of number of elements value representation and current actual equivalent damage degree its relax level mechanics equivalence;
L. identify damaged cable from the problematic support cable that k identified the step, remaining is exactly slack line;
M. utilize at current cable structure steady temperature data vector T iThe current actual damage vector of the cable system d that obtains in the k step under the condition iObtain slack line and current actual equivalent damage degree its relax level mechanics equivalence, utilize the f step obtain at current cable structure steady temperature data vector T iCurrent cable force vector F under the condition and current support cable two supporting end points horizontal ranges vectors, utilize the c step obtain at initial Cable Structure steady temperature data vector T oThe initial drift vector of the support cable under the condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum utilize current cable structure steady temperature data vector T iThe current steady temperature data of support cable of expression, utilize the c step obtain at initial Cable Structure steady temperature data vector T oThe support cable initial steady state temperature data of expression, utilization is in the temperature variant physical and mechanical properties parameter of the employed various materials of Cable Structure of c step acquisition, count temperature variation to the impact of support cable physics, mechanics and geometric parameter, by with slack line with damaged cable carry out the mechanics equivalence calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, the mechanics parameters of lax initial drift, geometrical property parameter, density and the material during with not damaged of the nothing of the rope of two equivalences is identical; Two, after the lax or damage, the Suo Li of the slack line of two equivalences and damage rope be out of shape after overall length identical; When satisfying above-mentioned two equivalent conditions, the such mechanics function of two support cables in Cable Structure is exactly identical, if after namely replacing damaged cable with equivalent slack line, Cable Structure any variation can not occur, vice versa; Try to achieve the relax level that those are judged as slack line according to aforementioned mechanics equivalent condition, relax level is exactly the change amount of support cable drift, has namely determined the long adjustment amount of rope of the support cable that those need adjust Suo Li; Lax identification and the damage identification of support cable have so just been realized; Institute's demand power is provided by current cable force vector F corresponding element during calculating;
N. try to achieve the vectorial d of current name damage i cAfter, set up mark vector B according to formula 3 i, formula 4 has provided mark vector B iThe definition of j element;
B i = B 1 i B 2 i . . . B j i . . . B N i T Formula 3
Formula 4
Element B in the formula 4 i jMark vector B iJ element, D i UjNominal unit damage vector D i uJ element, d i CjThe vectorial d of the current name damage of cable system i cJ element, they all represent the relevant information of j root support cable, j=1 in the formula 4,2,3 ..., N;
If mark vector B o. iElement be 0 entirely, then get back to step f and continue this circulation; If mark vector B iElement be not 0 entirely, then enter next step, be step p;
P. according to formula 5 calculate next time, i.e. the i+1 time current initial damage vector of the required cable system of circulation d I+1 oEach element;
d oj i + 1 = 1 - ( 1 - d oj i ) ( 1 - D uj i B j i ) Formula 5
D in the formula 5 I+1 OjThe current initial damage vector of the required cable system d that next time, namely circulates for the i+1 time I+1 oJ element, d i OjThe current initial damage vector of the cable system d that is this, namely circulates for the i time i oJ element, D i UjThe nominal unit damage vector D of the i time circulation i uJ element, B i jThe mark vector B of the i time circulation iJ element, j=1 in the formula 5,2,3 ..., N;
Q. take off once, i.e. the i+1 time required current initial Cable Structure steady temperature data vector T of circulation I+1 oEqual the current initial Cable Structure steady temperature data vector T of the i time circulation i o
R. at initial mechanical calculating benchmark model A oThe basis on, first to A oIn the Cable Structure bearing apply the constraint of current generalized displacement of support, the numerical value of current generalized displacement of support constraint is just taken from the numerical value of corresponding element among the generalized displacement of support 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, the health status that makes again rope is d I+1 oAfter obtain be exactly next time, i.e. the i+1 time required Mechanics Calculation benchmark model A of circulation I+1Obtain A I+1After, obtain A by Mechanics Calculation I+1In concrete numerical value all monitored amounts, current, these concrete numerical value form next time, the vectorial C of the current initial value of required monitored amount that namely circulates for the i+1 time I+1 oNext time, the i.e. required current initial Cable Structure bearing generalized coordinate vector U of the i+1 time circulation I+1 oEqual the current initial Cable Structure bearing generalized coordinate vector U of the i time circulation i o
S. get back to step f, beginning is circulation next time.
CN2012103793716A 2012-10-08 2012-10-08 Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support Pending CN102928243A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012103793716A CN102928243A (en) 2012-10-08 2012-10-08 Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103793716A CN102928243A (en) 2012-10-08 2012-10-08 Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support

Publications (1)

Publication Number Publication Date
CN102928243A true CN102928243A (en) 2013-02-13

Family

ID=47643096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103793716A Pending CN102928243A (en) 2012-10-08 2012-10-08 Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support

Country Status (1)

Country Link
CN (1) CN102928243A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616241A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring progressive identification method for defective cable and concentrated loads
CN103616231A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102252863A (en) * 2011-05-31 2011-11-23 东南大学 Strain monitoring-based progressive method for identifying slack cables during generalized displacement of support saddle
JP4863806B2 (en) * 2006-03-24 2012-01-25 東京製綱株式会社 Corrosion degree evaluation method and apparatus in the vicinity of hanger rope end part such as suspension bridge
CN102706605A (en) * 2012-05-30 2012-10-03 东南大学 Problem cable and support generalized displacement progressive identification method on basis of strain monitoring during temperature variation

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4863806B2 (en) * 2006-03-24 2012-01-25 東京製綱株式会社 Corrosion degree evaluation method and apparatus in the vicinity of hanger rope end part such as suspension bridge
CN102252863A (en) * 2011-05-31 2011-11-23 东南大学 Strain monitoring-based progressive method for identifying slack cables during generalized displacement of support saddle
CN102706605A (en) * 2012-05-30 2012-10-03 东南大学 Problem cable and support generalized displacement progressive identification method on basis of strain monitoring during temperature variation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
葛耀君,翟东,张国泉: "混凝土斜拉桥温度场的试验研究", 《中国公路学报》 *
颜东煌等: "混凝土斜拉桥施工控制温度影响及其现场修正", 《中国公路学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616241A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring progressive identification method for defective cable and concentrated loads
CN103616231A (en) * 2013-12-09 2014-03-05 东南大学 Generalized displacement strain monitoring progressive identification method for damaged cable and concentrated loads

Similar Documents

Publication Publication Date Title
CN102735462B (en) The slack line progressive-type recognition method of angular displacement of support and temperature variation hybrid monitoring
CN102735465B (en) Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation
CN102928243A (en) Slack cable progressing type identification method for general displacement and temperature variation strain monitoring of support
CN102706622A (en) Loosened cable progressive-type identification method based on angle monitoring at moment of generalized displacement of supporting seat and temperature variation
CN102706615B (en) The slack line progressive-type recognition method of generalized displacement of support temperature variation cable force monitoring
CN102706584B (en) The slack line progressive-type recognition method of support settlement and temperature variation space coordinate monitoring
CN102706623B (en) The slack line progressive-type recognition method of generalized displacement of support temperature variation hybrid monitoring
CN102721552B (en) Slack cable identification method based on hybrid monitoring during temperature change
CN102735477B (en) The slack line progressive-type recognition method of strain monitoring during angular displacement of support temperature variation
CN102735478B (en) The slack line progressive-type recognition method of cable force monitoring during angular displacement of support temperature variation
CN102735463B (en) The slack line progressive-type recognition method of angle monitor during angular displacement of support temperature variation
CN102706618B (en) Slack cable approximant identification method on basis of strain monitoring during temperature variation
CN102706648B (en) The slack line progressive-type recognition method of cable force monitoring when support settlement and temperature variation
CN102706650B (en) The slack line progressive-type recognition method of strain monitoring when support settlement and temperature variation
CN102706586B (en) The slack line progressive-type recognition method of angle monitor when support settlement and temperature variation
CN102706675B (en) Angular displacement of support temperature variation space coordinate monitoring slack line progressive-type recognition method
CN102706617A (en) Loosened cable identification method based on angle monitoring at moment of generalized displacement of supporting seat and temperature variation
CN102735464A (en) Angle-monitoring-based slack cable identification method during angular support displacement and temperature change
CN102706656B (en) Based on the approximant recognition methods of the slack line of angle monitor during temperature variation
CN102706589B (en) Based on the approximant recognition methods of the slack line of space coordinate monitoring during temperature variation
CN102706613B (en) Generalized displacement of support temperature variation space coordinate monitoring slack line progressive-type recognition method
CN102706599B (en) The slack line progressive-type recognition method of hybrid monitoring when support settlement and temperature variation
CN102706616B (en) Generalized displacement of support temperature variation is based on the slack line recognition methods of hybrid monitoring
CN102706597B (en) Based on the approximant recognition methods of the slack line of hybrid monitoring during temperature variation
CN102721553B (en) Strain monitoring-based slack line identification method in existence of support settlement and 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
C05 Deemed withdrawal (patent law before 1993)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130213