CN102735465B - Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation - Google Patents

Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation Download PDF

Info

Publication number
CN102735465B
CN102735465B CN201210173656.4A CN201210173656A CN102735465B CN 102735465 B CN102735465 B CN 102735465B CN 201210173656 A CN201210173656 A CN 201210173656A CN 102735465 B CN102735465 B CN 102735465B
Authority
CN
China
Prior art keywords
cable
temperature
data
vector
initial
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.)
Expired - Fee Related
Application number
CN201210173656.4A
Other languages
Chinese (zh)
Other versions
CN102735465A (en
Inventor
韩玉林
郑可
韩佳邑
王芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201210173656.4A priority Critical patent/CN102735465B/en
Publication of CN102735465A publication Critical patent/CN102735465A/en
Application granted granted Critical
Publication of CN102735465B publication Critical patent/CN102735465B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

During angular displacement of support temperature variation based on the slack line recognition methods of strain monitoring based on strain monitoring, the Mechanics Calculation benchmark model needing to upgrade Cable Structure is determined whether by monitoring angular displacement of support, Cable Structure temperature and environment temperature, obtain the Mechanics Calculation benchmark model of the Cable Structure counting angular displacement of support, Cable Structure temperature and environment temperature, the basis of this model calculates and obtains unit damage monitored amount transformation matrices.Calculate the noninferior solution of the current nominal fatigue vector of cable system with the linear approximate relationship existed between monitored amount current initial value vector, unit damage monitored amount transformation matrices, unit damage scalar sum cable system current nominal fatigue vector to be asked according to monitored amount current value vector, after use lossless detection method identifies true damaged cable, the rope of remaining unsoundness problem is exactly slack line.

Description

Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation
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 bearing tensile load, as suspension cable, main push-towing rope, hoist cable, pull bar etc., the common ground of this class formation is with rope, cable or only bears the rod member of tensile load for support unit, and such structure representation is " Cable Structure " by this method for simplicity.Along with the change of environment temperature, the temperature of Cable Structure also can change, when Cable Structure temperature changes, there being angular displacement of support, (such as bearing is around coordinate axis X, Y, the rotation of Z, in fact be exactly that bearing is around coordinate axis X, Y, the angular displacement of Z) time, based on strain monitoring, this method identifies that the supporting system of Cable Structure (refers to all ropeway carrying-ropes, and all rod members only bearing tensile load play supporting role, for simplicity, whole support units of this class formation are collectively referred to as " cable system " by this patent, but in fact cable system not only refers to support cable, also the rod member only bearing tensile load is comprised, censure all ropeway carrying-ropes and all rod members only bearing tensile load play supporting role with " support cable " this noun in this method) in the support cable (the impaired rod member only bearing tensile load is just referred to truss-frame structure) of need adjustment Suo Li, belong to engineering structure health monitoring field.
Background technology
Cable system is Cable Structure (particularly large-scale Cable Structure normally; such as large-scale cable-stayed bridge, suspension bridge) key components; due to the reason such as lax; the be completed Suo Li of a period of time rear support rope of new construction can change usually; the lax change that also can cause supporting cable force of its support cable after structure long service; these changes all will cause the change of structural internal force; harmful effect is caused to the safety of structure; the inefficacy of structure will be caused time serious, therefore identify that the support cable that need adjust Suo Li is very important accurately and timely.
The health status of support cable system changes and (such as occurs lax, damage etc.) after, the change of the measurable parameter of structure can be caused, the distortion of such as Cable Structure or strain can change, in fact the change strained contains the health status information of cable system, that is structural strain data can be utilized to judge the health status of structure, can (monitored strain be called " monitored amount " by this method based on strain monitoring, mention " monitored amount " below and just refer to monitored strain) identify that (this method is also referred to as the support cable of unsoundness problem to damaged cable, refer to that support cable is impaired, relax or have both at the same time).Monitored amount is except being subject to the impact of cable system health status; also can be subject to the impact of Cable Structure temperature variation (usually can occur) and Cable Structure angular displacement of support (usually can occur); change with under the condition of Cable Structure bearing generation angular displacement in Cable Structure temperature; if can based on the identification monitoring of monitored amount being realized to the support cable to unsoundness problem; to the safety of Cable Structure, there is important value, also do not have a kind of disclosed, effective health monitoring systems and method to solve this problem at present.
When Cable Structure has angular displacement of support and temperature variation, in order to can monitoring reliably be had to the health status of the cable system of Cable Structure and judge, must have one can rationally effectively set up each monitored amount change with ropes all in cable system health status between the method for relation, the health monitoring systems set up based on the method can provide the health evaluating of more believable cable system.
Summary of the invention
Technical matters: the object of this method is when Cable Structure has support displacement and temperature variation, for the health monitoring problem of cable system in Cable Structure, disclose a kind of based on strain monitoring, health monitor method that cable system in Cable Structure can be monitored rationally and effectively.
According to the reason that the Suo Li of support cable changes, the Suo Li of support cable change can be divided into two kinds of situations: one is that support cable receives damage, and localized cracks and corrosion etc. have appearred in such as support cable; Two are support cable and not damaged, but Suo Li also there occurs change, occur that the one of the main reasons of this change is that Suo Changdu (be called drift, this method specially refers to that support cable two supports the drift of that section of rope between end points) under support cable free state (now rope tensility also claims Suo Li to be 0) there occurs change.One of fundamental purpose of this method will identify the support cable that drift there occurs change exactly, and identifies the knots modification of their drift, and this knots modification is that the cable force adjustment of this rope provides direct basis.The reason that support cable drift changes is not single, and conveniently, the support cable that drift changes by this method is referred to as slack line.Refer to slack line recognition system with cable system health monitoring systems in the method, refer to slack line recognition methods by cable system health evaluating method, in other words this method " health monitoring " usually available " slack line identification " substitute.
Technical scheme: this method is made up of three parts.Method, knowledge based storehouse (containing parameter) and actual measurement the cable system health state evaluation method of monitored amount, the software and hardware part of health monitoring systems of setting up knowledge base needed for cable system health monitoring systems and parameter respectively.
The Part I of this method: the method setting up knowledge base needed for cable system health monitoring systems and parameter.Specific as follows:
1. first determine " the temperature survey calculating method of the Cable Structure of this method ".Temperature due to Cable Structure may be change, the temperature of the different parts of such as Cable Structure changes along with the change of intensity of sunshine, along with the change of environment temperature changes, the surface of Cable Structure may be time dependent with inner temperature sometimes, the surface of Cable Structure may be different from inner temperature, the surface of Cable Structure is time dependent with inner temperature difference, this just makes the Mechanics Calculation of Cable Structure when considering temperature conditions and monitoring quite complicated, for simplifying problem, reduce calculated amount and reduce and measure cost, especially in order to improve computational accuracy, this method proposes " the temperature survey calculating method of the Cable Structure of this method ", specific as follows:
The first step, inquiry or actual measurement obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day r.Be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of describing hunit be DEG C/m.The surface of Cable Structure is got " R Cable Structure surface point ", the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " describes later with the condition that must meet that distributes.From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, 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 moment.
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 the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and 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 ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax.
3rd step, survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above emaxbe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T pmaxΔ T is not more than after deducting 2 degrees Celsius emax, and Cable Structure surface maximum temperature difference Δ T smaxbe not more than Δ T pmax, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining 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 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius, this method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation hthe numerical value obtained, gets Δ T for convenience of describing hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure.
2. set up the initial mechanical Calculation Basis model A of Cable Structure o(such as finite element benchmark model) and current initial mechanical Calculation Basis model A t othe method of (such as finite element benchmark model), sets up and A ocorresponding monitored amount initial value vector C omethod, set up and A t ocorresponding monitored amount current initial value vector C t omethod.A in the method oand C oconstant.A t oand C t oconstantly update.Set up A oand C o, set up and upgrade A t oand C t omethod as follows.
If total N root support cable, first determines the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step." the whole monitored strain data of structure " can by the specified point of K in structure and the strain of L assigned direction of each specified point describe, the change of structural strain data is exactly the change of all strains of K specified point.Each total individual strain measurement value of M (M=K × L) or calculated value carry out characterisation of structures strain information.K and M must not be less than the quantity N of support cable.
For simplicity, in the method by " the monitored strain data of structure " referred to as " monitored amount ".When mentioning " so-and-so matrix of monitored amount or so-and-so vector " later, also can be read as " strain so-and-so matrix or so-and-so vector ".
Set up initial mechanical Calculation Basis model A otime, when Cable Structure is completed, or before setting up health monitoring (damaged cable identification) system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T orepresent, be called initial Cable Structure steady temperature data vector T o.T is obtained in actual measurement owhile, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, use the direct survey calculation of conventional method to obtain the initial number of all monitored amount of Cable Structure.Conventional method (consult reference materials or survey) is used to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained at Actual measurement owhile, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing spatial data, initial Cable Structure bearing angular data, Cable Structure modal data, structural strain data, structure angle measurement of coordinates data, structure space measurement of coordinates data that the Non-destructive Testing Data of support cable etc. can express the health status of rope.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 in structure, and object is the geometric properties according to these coordinate data determination Cable Structure.For cable-stayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, so-called bridge type data that Here it is.The data utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope set up cable system initial damage vector d o(as the formula (1)), d is used orepresent that Cable Structure is (with initial mechanical Calculation Basis model A orepresent) the initial health of cable system.If there is no the Non-destructive Testing Data of rope and other are when can express the data of the health status of support cable, or when can think that structure original state is not damaged state, vectorial d oeach element numerical value get 0.The temperature variant physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure to use and initial Cable Structure steady temperature data vector T o, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A o.Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U o.
d o=[d o1d o2···d oj···d oN] T(1)
D in formula (1) oj(j=1,2,3 ...., N) represent initial mechanical Calculation Basis model A oin the initial damage value of jth root rope of cable system, d ojrepresent jth root rope not damaged when being 0, represent that this rope thoroughly loses load-bearing capacity when being 100%, represent that the load-bearing capacity of corresponding proportion lost by jth root rope time between 0 and 100%, T represents the transposition (same afterwards) of vector.
T is obtained in actual measurement owhile, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, the initial value of all monitored amount of the Cable Structure using the direct survey calculation of conventional method to obtain, forms monitored amount initial value vector C o(see formula (2)).Require at acquisition A owhile obtain C o, monitored amount initial value vector C orepresent and correspond to A othe concrete numerical value of " monitored amount ".Because of subject to the foregoing, the Calculation Basis model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in describing below, will represent this calculated value and measured value with prosign.
C o=[C o1C o2···C oj···C oM] T(2)
C in formula (2) oj(j=1,2,3 ...., M) be the original bulk of jth monitored amount in Cable Structure, this component corresponds to a specific jth monitored amount according to coding rule.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by o, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T o), based on A othe Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A othe Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, Cable Structure angle-data, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A othe health status of middle support cable represents with cable system initial damage vector do, the initial Cable Structure steady temperature data vector T of Cable Structure steady temperature data orepresent.Due to based on A othe initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A obasis on, carry out Mechanics Calculation obtains, A othe evaluation of each monitored amount form monitored amount initial value vector C o.T can be said o, U oand d oa oparameter, C oby A omechanics Calculation result composition.
Set up and upgrade current initial mechanical Calculation Basis model A t omethod be: initial time (namely first time set up A t otime), A t ojust equal A o, A t ocorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T t o", at initial time, T t ojust equal T o, vector T t odefinition mode and vector T odefinition mode identical.Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure t ocable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U t o, the current initial mechanical Calculation Basis model A of Cable Structure is namely set up for the first time at initial time t otime, U t ojust equal U o.A t othe initial health of support cable and A othe health status of support cable identical, also use cable system initial damage vector d orepresent, A in cyclic process below t othe initial health of support cable use cable system initial damage vector d all the time orepresent; Cable Structure is in A t oduring state, this method monitored amount current initial value vector C t orepresent the concrete numerical value of all monitored amounts, C t oelement and C oelement one_to_one corresponding, represent that all monitored amounts are in A in Cable Structure respectively t oand A oconcrete numerical value during two states.At initial time, C t ojust equal C o, T t o, U t oand d oa t oparameter, C t oby A t omechanics Calculation result composition; In Cable Structure military service process, the current data obtaining " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical); Obtaining vector T twhile, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current data composition current cable structure actual measurement bearing angular coordinate vector U t; If T tequal T t oand U tequal U t o, then do not need A t oupgrade, otherwise need A t o, U t oand T t oupgrade, update method is: the first step calculates U twith U odifference, U twith U odifference be exactly the front holder angular displacement of Cable Structure bearing about initial position, angular displacement of support is represented with angular displacement of support vector V, be one-to-one relationship between element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing; Second step calculates T twith T odifference, T twith T odifference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T odifference represent with steady temperature change vector S, S equals T tdeduct T o, S represents the change of Cable Structure steady temperature data; 3rd step is first to A oin Cable Structure bearing apply front holder angular displacement constraint, the numerical value of front holder angular displacement constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A omiddle Cable Structure bearing applies angular displacement of support constraint and to A oin Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A t o, upgrade A t owhile, U t oall elements numerical value also uses U tall elements numerical value correspondence replaces, and namely have updated U t o, T t oall elements numerical value also uses T tall elements numerical value correspondence replace, namely have updated T t o, so just obtain and correctly correspond to A t ot t o; Upgrade C t omethod be: when renewal A t oafter, obtain A by Mechanics Calculation t oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C t o.
In Cable Structure the currency of all monitored amounts form monitored amount current value vector C(definition see formula (3)).
C=[C 1C 2···C j···C M] T(3)
C in formula (3) j(j=1,2,3 ...., M) be the currency of jth monitored amount in Cable Structure, this component C jaccording to coding rule and C ojcorresponding to same " monitored amount ".Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C.
3. set up and upgrade the method for Cable Structure unit damage monitored amount transformation matrices Δ C.
Cable Structure unit damage monitored amount transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal t owith monitored amount current initial value vector C t owhile, upgrade Cable Structure unit damage monitored amount transformation matrices Δ C.Concrete grammar is as follows:
At the current initial mechanical Calculation Basis model A of Cable Structure t obasis on carry out several times calculating, calculation times numerically equals the quantity of all support cables.Calculating each time in hypothesis cable system only has a support cable (to use vectorial d at initial damage ocorresponding element represent) basis on increase unit damage D again u(such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), occur in calculating each time that the rope damaged is different from during other time calculates the rope occurring damaging, calculate the current calculated value all utilizing mechanics method (such as finite element method) to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, and (when hypothesis i-th rope has unit damage, available formula (4) represents monitored amount calculation current vector ); Calculate monitored amount calculation current vector each time and deduct monitored amount current initial value vector C t o, gained vector is exactly that the monitored amount change vector of (to have the position of the support cable of unit damage or numbering etc. for mark) (when i-th rope has unit damage, uses δ C under this condition irepresent monitored amount change vector, formula (5) is shown in definition), each element representation of monitored amount change vector is owing to suppose there is the knots modification of the monitored amount corresponding to the unit damage of the Na Gensuo of unit damage and this element of causing when calculating; N root rope is had just to have N number of monitored amount change vector, owing to there being N number of monitored amount, so each monitored amount change vector has N number of element, be made up of the unit damage monitored amount transformation matrices Δ C having M × N number of element successively this N number of monitored amount change vector, the definition of Δ C as the formula (6).
C t i = C t 1 i C t 2 i · · · C tj i · · · C tM i T - - - ( 4 )
Elements C in formula (4) tj i(i=1,2,3 ...., N; J=1,2,3 ...., M) represent due to i-th rope have a unit damage time, according to the current calculated amount of the monitored amount of the jth corresponding to coding rule.
δC i = C t i - C o i - - - ( 5 )
ΔC = ΔC 1,1 ΔC 1,2 · ΔC 1 , i · ΔC 1 , N ΔC 2,1 ΔC 2,2 · ΔC 2 , i · ΔC 2 , N · · · · · · ΔC j , 1 ΔC j , 2 · ΔC j , i · ΔC j , N · · · · · · ΔC M , 1 ΔC M , 2 · ΔC M , i · ΔC M , N - - - ( 6 )
Δ C in formula (6) j,i(i=1,2,3 ...., N; J=1,2,3 ...., M) represent only because i-th rope has unit damage to cause, according to the change (algebraic value) of the calculating current value of the monitored amount of the jth corresponding to coding rule.Monitored amount change vector δ C ibe actually the row in matrix Δ C, that is formula (6) also can be write as formula (7).
ΔC=[δC 1δC 2···δC i···δC N] (7)
4. monitored amount current value vector C(calculates or actual measurement) with monitored amount current initial value vector C t o, unit damage monitored amount transformation matrices Δ C, unit damage scalar D uand the linear approximate relationship between cable system current nominal fatigue vector d, shown in (8) or formula (9).The definition of cable system current nominal fatigue vector d is see formula (10).
C = C o t + 1 D u ΔC · d - - - ( 8 )
C - C o t = 1 D u ΔC · d - - - ( 9 )
d=[d 1d 2···d i···d N] T(10)
D in formula (10) i(i=1,2,3 ...., N) be the current nominal fatigue of i-th rope (or pull bar) in cable system.
If set rope damage as 100% time represent that rope thoroughly loses load-bearing capacity, so (be such as not more than the damage of 30%) when actual damage is not too large, because Cable Structure material is still in the linear elasticity stage, the distortion of Cable Structure is also less, formula (8) or a kind of like this linear relationship represented by formula (9) less with the error of actual conditions.The error of the linear relationship error vector e expression (8) defined by formula (11) or the shown linear relationship of formula (9).
e = abs ( 1 D u ΔC · d - C + C o t ) - - - ( 11 )
In formula (11), abs () is the function that takes absolute value, and takes absolute value to each element of the vector of trying to achieve in bracket.
The Part II of this method: the cable system health state evaluation method of knowledge based storehouse (containing parameter) and the monitored amount of actual measurement.
There is certain error in the linear relationship represented by formula (8) or formula (9), therefore simply can not carry out direct solution according to formula (8) or formula (9) and actual measurement monitored amount current value vector C and obtain cable system current nominal fatigue vector d.If this has been doned, the element in the cable system obtained current nominal fatigue vector d even there will be larger negative value, namely negative damage, and this is obviously irrational.Therefore the acceptable solution of cable system current nominal fatigue vector d is obtained (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more accurately from cable system) become a rational solution, available formula (12) expresses this method.
abs ( 1 D u ΔC · d - C + C o t ) ≤ g - - - ( 12 )
In formula (12), abs () is the function that takes absolute value, and vectorial g describes the legitimate skew departing from ideal linearity relation (formula (8) or formula (9)), is defined by formula (13).
g=[g 1g 2···g j···g M] T(13)
G in formula (13) j(j=1,2,3 ...., M) describe the maximum allowable offset departing from formula (8) or the ideal linearity relation shown in formula (9).The error vector e tentative calculation that vector g can define according to formula (11) is selected.
At monitored amount current initial value vector C t o, unit damage monitored amount transformation matrices Δ C, survey monitored amount current value vector C and unit damage D uwhen (setting before calculating Δ C, is scalar) is known, suitable algorithm (such as multi-objective optimization algorithm) can be utilized to solve formula (12), obtain the acceptable solution of cable system current nominal fatigue vector d.
Definition cable system current actual damage vector d a(see formula (14)), cable system current actual damage vector d aelement can calculate according to formula (15), namely obtain Suo Dangqian actual damage vector d a, thus can by d adetermine position and the degree of injury of damaged cable, namely achieve the health monitoring of cable system, achieve damaged cable identification.
d a = d 1 a d 2 a · · · d j a · · · d N a T - - - ( 14 )
D in formula (14) a j(j=1,2,3 ...., N) represent the actual damage value of jth root rope, its definition is shown in formula (15), d a jrepresent jth root rope not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, represent time between 0 and 100% that the load-bearing capacity of corresponding proportion lost by jth root rope, vectorial d athe coding rule of element and formula (1) in vectorial d othe coding rule of element identical.
d j a = 1 - ( 1 - d oj ) ( 1 - d j ) - - - ( 15 )
D in formula (15) oj(i=1,2,3,4, J=1,2,3 ...., N) be vectorial d oa jth element, d jit is a jth element of vectorial d.Thus determine position and the degree of injury of damaged cable.
Describe below and obtain cable system current actual damage vector d aafter, how to determine position and the relax level of slack line.
If total N root support cable in cable system, structure rope force data is described by the Suo Li of N root support cable.Available " Initial cable force vector F o" represent the Initial cable force (formula (16) is shown in definition) of all support cables in Cable Structure.
F o=[F o1F o2···F oj···F oN] T(16)
F in formula (16) o(j=1,2,3 ...., N) be the Initial cable force of jth root support cable in Cable Structure, this element corresponds to the Suo Li specifying support cable according to coding rule.Vector F oit is constant.T is obtained in actual measurement osynchronization, use conventional method direct survey calculation to obtain the rope force data of all support cables, all these rope force datas composition Initial cable force vector F o.Setting up the initial mechanical Calculation Basis model A of Cable Structure otime in fact employ vectorial F o.
Represent the current cable power (formula (17) is shown in definition) of all support cables in the Cable Structure of surveying and obtaining with " current cable force vector F " in this method.
F=[F 1F 2···F j···F N] T(17)
F in formula (17) j(j=1,2,3 ...., N) be the current cable power of jth root support cable in Cable Structure.Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F.The element of vector F and vectorial F othe coding rule of element identical.According to describing above, vector T t oequal vector T t.
In this method, under support cable original state, at the initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure oduring expression, and when support cable is in free state (free state refers to that Suo Li is 0, rear same), the length of support cable is called initial drift, with " initial drift vector l o" represent the initial drift (formula (18) is shown in definition) of all support cables in Cable Structure.Vector T is passed through 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 in moment.
l o=[l o1l o2···l oj···l oN] T(18)
L in formula (18) oj(j=1,2,3 ...., N) be the initial drift of jth root support cable in Cable Structure.Vector l obeing constant, after determining when starting, just no longer changing.
Similar, under support cable original state, at the initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure oduring expression, and when support cable is in free state, the cross-sectional area of support cable is called initial free cross-sectional area, with " initial free cross-sectional area vector A o" representing the initial free cross-sectional area (formula (19) is shown in definition) of all support cables in 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" represent the weight (formula (20) is shown in definition) of the initial free unit length of all support cables in Cable Structure.
A o=[A o1A o2···A oj···A oN] T(19)
A in formula (19) oj(j=1,2,3 ...., N) be the initial free cross-sectional area of jth root support cable in Cable Structure.Vector A obeing constant, after determining when starting, just no longer changing.
ω o=[ω o1ω o2···ω oj···ω oN] T(20)
ω in formula (20) oj(j=1,2,3 ...., N) be the weight of the free unit length of initial freedom of jth root support cable in Cable Structure.Vector ω obeing constant, after determining when starting, just no longer changing.
In this method, at the current initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure t oduring expression, with " current initial drift vector l t o" (formula (21) is shown in definition, when referring to that hypothesis supporting cable force is 0, considers thermal expansivity and temperature variation to after the impact of support cable drift, initial drift vector l to represent the current initial drift of all support cables in Cable Structure owith initial Cable Structure steady temperature data vector T othe support cable represented is at the current initial Cable Structure steady temperature data vector T of temperature t osupport cable drift during expression).According to " the temperature survey calculating method of the Cable Structure of this method ", pass through vector T t ocan determine obtaining vector T t othe Temperature Distribution of all support cables in moment.
l o t = l o 1 t l o 2 t · · · l oj t · · · l oN t T - - - ( 21 )
L in formula (21) t oj(j=1,2,3 ...., N) be the current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, in Cable Structure, the current initial drift of jth root support cable, can utilize the thermal expansivity of support cable, l oj, T oand T t ol is calculated by Typical physical 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, at the current initial Cable Structure steady temperature data vector T of steady temperature data of Cable Structure t oduring expression, represent the current drift (formula (22) is shown in definition, and now support cable may be intact, also may be impaired, also may be lax) of all support cables in Cable Structure with " current drift vector l ".
l=[l 1l 2···l j···l N] T(22)
L in formula (22) j(j=1,2,3 ...., N) be the current drift of jth root support cable in Cable Structure.
In this method, represent the knots modification (formula (23) and formula (24) are shown in definition) of the drift of all support cables in Cable Structure with " drift changes vectorial Δ l " (or claiming support cable current slack degree vector).
Δl=[Δl 1Δl 2···Δl j···Δl N] T(23)
Δ l in formula (23) j(j=1,2,3 ...., N) be the knots modification of the drift of jth root support cable in current cable structure, its definition is shown in formula (24), Δ l jbe not 0 rope be slack line, Δ l jnumerical value be the slack of rope, and representing the current slack degree of cable system jth root support cable, is also the long adjustment amount of rope of this rope during adjustment Suo Li.
Δ l j = l j - l oj t - - - ( 24 )
In the method by slack line is carried out with damaged cable the relax level identification that mechanic equivalent carries out slack line, the mechanical condition of equivalence is:
One, the nothing of the rope of two equivalences is lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter and material;
Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical.
When meeting above-mentioned two equivalent conditions, such two support cables mechanics function is in the structure exactly identical, if after namely replacing slack line with the damaged cable of equivalence, any change can not occur Cable Structure, and vice versa.
Obtain cable system current actual damage vector d aafter, d aa jth element d a j(j=1,2,3 ...., N) represent the actual damage value of jth root rope, its definition is shown in formula (15), although by d i jbe called the actual damage value of jth root rope or the actual damage degree of jth root rope, but also may be lax because jth root Suo Keneng is impaired, so d aa jth element d a jthe actual damage value of the jth root rope represented is actually the actual equivalent damage value of jth root rope, when jth root rope is actually impaired, and d a jthe actual damage value of the jth root rope just represented, when jth root rope is actually lax, d a jthe jth root rope just represented with the actual damage value of lax equivalence, for sake of convenience, claim d in the method a jrepresent jth root rope not damaged when being 0, when being 100%, represent that this rope thoroughly loses load-bearing capacity, time between 0 and 100%, represent that the load-bearing capacity of corresponding proportion lost by jth root rope, by cable system current actual damage vector d ajust can identify the support cable that health status goes wrong, but in the support cable that goes wrong of these health status, some is impaired, some relaxes, if a jth support cable is actually relax (its current slack degree Δ l jdefinition), the current slack degree Δ l of a so lax jth support cable j(Δ l jdefinition see formula (23)) with the current actual damage degree d of damaged cable of equivalence a jbetween relation determined by aforementioned two mechanic equivalent conditions.Δ l jsame d a jbetween physical relationship can adopt accomplished in many ways, such as can directly determine (see formula (25)) according to aforementioned equivalent condition, also can adopt after replacing the E in formula (25) to revise based on Ernst equivalent elastic modulus and determine (see formula (26)), other method such as trial and error procedure based on finite element method also can be adopted to determine.
Δ l j t = d j a 1 - d j a F j E j t A j t + F j l oj t - - - ( 25 )
Δ l j t = d j a 1 - d j a F h [ 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 - - - ( 26 )
Formula (25) and the middle E of formula (26) t jthe current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, the elastic modulus of a jth support cable, A t jthe current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, the cross-sectional area of a jth support cable, F jthe current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, the current cable power of a jth support cable, d a jthe current actual damage degree of a jth support cable, ω t jthe current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, the weight of the unit length of a jth support cable, l t jxthe current initial Cable Structure steady temperature data vector T of steady temperature data in Cable Structure t oduring expression, the horizontal range of two supporting end points of a jth support cable.E t jcan obtain according to the characteristic material data looking into or survey a jth support cable, A t jand ω t jcan according to the thermal expansivity of a jth support cable, A oj, ω oj, F j, T oand T t oobtained by Typical physical and Mechanics Calculation.Item in formula (26) in [] is the Ernst equivalent elastic modulus of this support cable, just can determine support cable current slack degree vector Δ l by formula (25) or formula (26).Formula (26) is the correction to formula (25).
The Part III of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (comprise monitored amount monitoring system, horizontal range monitoring system, Cable Structure bearing angular coordinate monitoring system that temperature monitoring system, cable force monitoring system, support cable two support end points), signal picker and computing machine etc.Require that Real-Time Monitoring obtains temperature required measured data, require that the Suo Li of each monitored amount of Real-Time Monitoring, simultaneously each support cable of Real-Time Monitoring simultaneously, simultaneously each support cable two of Real-Time Monitoring support the data of horizontal range, simultaneously each Cable Structure bearing angular coordinate of Real-Time Monitoring of end points.
Software should to complete in this method required, can by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
This method specifically comprises:
A. establish total N root support cable, first determine the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; Determine the point being monitored of specifying, namely point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified points numbering; Determine monitored should the changing direction of point being monitored, and give all monitored strain numberings of specifying; " monitored strain numbering " will be used for generating vector sum matrix in subsequent step; " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity of point being monitored 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 moment;
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 environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of describing hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, especially, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, namely the temperature of the surface point of support cable is only measured, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, 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 the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and 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 ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax,
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum error Δ T that survey calculation obtains above emaxbe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T pmaxΔ T is not more than after deducting 2 degrees Celsius emax, and Cable Structure surface maximum temperature difference Δ T smaxbe not more than Δ T pmax, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining 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 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius, this method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation hthe numerical value obtained, gets Δ T for convenience of describing hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under 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 various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement osynchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time cross-sectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T othe weight of the unit length of all support cables when the cross-sectional area of all support cables and Suo Li are 0 when the length of all support cables, Suo Li are 0 when Suo Li under condition is 0, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule 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 and Initial cable force vector F othe coding rule of element identical; T is obtained in actual measurement owhile, namely at the initial Cable Structure steady temperature data vector T of acquisition othe synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Non-destructive Testing Data of the health status expressing support cable, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing spatial data, initial Cable Structure bearing angular data, initial Cable Structure angle-data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C o; The Non-destructive Testing Data that utilization can express the health status of support cable sets 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 support cable corresponding to this element is intact, no problem, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the load-bearing capacity of corresponding proportion, if there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or when thinking that Cable Structure original state is not damaged state, vectorial d oeach element numerical value get 0; If d othe numerical value of a certain element be not 0, then represent that the support cable corresponding to this element is problematic, this support cable may be impaired in the method also may be lax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the initial equivalent damage degree of the support cable of this its correspondence of element numeric representation; Cable system initial damage vector d othe coding rule of element and Initial cable force vector F othe coding rule of element identical; Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U o;
Temperature variant physical and mechanical properties parameter, the initial Cable Structure bearing angular coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure o, initial Cable Structure steady temperature data vector T owith all Cable Structure data obtained with preceding step, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " o, based on A othe Cable Structure that calculates calculates data must closely 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 angular data be exactly initial Cable Structure bearing angular coordinate vector U o; Corresponding to A osupport cable health status with cable system initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " t o, monitored amount current initial value vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t owith monitored amount current initial value vector C t otime, the current initial mechanical Calculation Basis model A of Cable Structure t ojust equal the initial mechanical Calculation Basis model A of Cable Structure o, monitored amount current initial value vector C t ojust equal monitored amount initial value vector C o; A t ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t otime, T t ojust equal T o; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure t ocable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U t o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t otime, U t ojust equal U o; A t othe initial health of support cable and A othe health status of support cable identical, also use cable system initial damage vector d orepresent, A in cyclic process below t othe initial health of support cable use cable system initial damage vector d all the time orepresent; Work as T o, U oand d oa oparameter time, by A othe initial value of all monitored amount that obtains of Mechanics Calculation result and C othe initial value of all monitored amount represented is identical, therefore alternatively C oby A omechanics Calculation result composition, work as T t o, U t oand d oa t oparameter time, C t oby A t omechanics Calculation result composition; A in the method o, U o, C o, d oand T oconstant;
E. from entering the circulation being walked to o step by e here; In Cable Structure military service process, the current data of " Cable Structure steady temperature data " is constantly obtained according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the difference of the volume coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector othe coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas composition current cable structure actual measurement bearing angular coordinate vector U t;
F. according to current cable structure actual measurement bearing angular coordinate vector U twith current cable structure steady temperature data vector T t, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3 t o, current initial Cable Structure bearing angular coordinate vector U t o, monitored amount current initial value vector C t owith current initial Cable Structure steady temperature data vector T t o;
F1. U is compared respectively twith U t o, T twith T t oif, U tequal U t oand T tequal T t o, then A t o, U t o, C t oand T t oremain unchanged; Otherwise need to follow these steps to A t o, U t oand T t oupgrade;
F2. U is calculated twith U odifference, U twith U odifference be exactly the front holder angular displacement of Cable Structure bearing about initial position, with angular displacement of support vector V represent angular displacement of support, V equals U tdeduct U o, be one-to-one relationship between the element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing; Calculate T twith T odifference, T twith T odifference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T odifference represent with steady temperature change vector S, S equals T tdeduct T o, S represents the change of Cable Structure steady temperature data;
F3. first to A oin Cable Structure bearing apply front holder angular displacement constraint, the numerical value of front holder angular displacement constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A omiddle Cable Structure bearing applies angular displacement of support constraint and to A oin Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A t o, upgrade A t owhile, U t oall elements numerical value also uses U tall elements numerical value correspondence replaces, and namely have updated U t o, T t oall elements numerical value also uses T tall elements numerical value correspondence replace, namely have updated T t o, so just obtain and correctly correspond to A t ot t oand U t o; Upgrade C t omethod be: when renewal A t oafter, obtain A by Mechanics Calculation t oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C t o; A t othe initial health of support cable use cable system initial damage vector d all the time orepresent;
G. at current initial mechanical Calculation Basis model A t obasis on carry out several times Mechanics Calculation according to step g 1 to g4, obtain Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D by calculating u;
G1. Cable Structure unit damage monitored amount transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal t o, current initial Cable Structure bearing angular coordinate vector U t o, monitored amount current initial value vector C t owith current initial Cable Structure steady temperature data vector T t oafterwards, Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D must then be upgraded u;
G2. at the current initial mechanical Calculation Basis model A of Cable Structure t obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity of all ropes, has N root support cable just to have N calculating, each time calculate hypothesis cable system in only have a support cable to have unit damage scalar D uoccur in calculating each time that the rope damaged is different from during other time calculates the rope occurring damaging, calculate the current calculated value of all monitored amounts in Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector oelement number rule identical;
G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C t oobtain a monitored amount change vector; N root support cable is had just to have N number of monitored amount change vector;
G4. the Cable Structure unit damage monitored amount transformation matrices Δ C having N to arrange is made up of successively this N number of monitored amount change vector; Each row of Cable Structure unit damage monitored amount transformation matrices Δ C correspond to a monitored amount change vector;
H. current cable structure steady temperature data vector T is obtained in actual measurement twhile, actual measurement obtains at acquisition current cable structure steady temperature data vector T tthe current measured value of all monitored amount of Cable Structure of synchronization in moment, form monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C t owith monitored amount initial value vector C odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time;
I. cable system current nominal fatigue vector d is defined, the element number of cable system current nominal fatigue vector d equals the quantity of support cable, be one-to-one relationship between the element of cable system current nominal fatigue vector d and support cable, the element numerical value of cable system current nominal fatigue vector d represents the nominal fatigue degree of corresponding support cable or nominal health status; The coding rule of the element of vector d and vectorial d othe coding rule of element identical;
J. according to monitored amount current value vector C with monitored amount current initial value vector C t o, Cable Structure unit damage monitored amount transformation matrices Δ C, unit damage scalar D uand the linear approximate relationship existed between cable system to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, and other amount in formula 1 except d is known, solves formula 1 and just can calculate cable system current nominal fatigue vector d;
C = C o t + 1 D u ΔC · d Formula 1
K. cable system current actual damage vector d is defined a, cable system current actual damage vector d aelement number equal the quantity of support cable, cable system current actual damage vector d aelement and support cable between be one-to-one relationship, cable system current actual damage vector d aelement numerical value represent the actual damage degree of corresponding support cable or actual health status; Vector d athe coding rule of element and vectorial d othe coding rule of element identical;
L. the cable system utilizing formula 2 to express current actual damage vector d aa jth element d a jwith cable system initial damage vector d oa jth element d o jwith a jth element d of cable system current nominal fatigue vector d jbetween relation, calculate cable system current actual damage vector d aall elements;
d j a = 1 - ( 1 - d oj ) ( 1 - d j ) Formula 2
J=1 in formula 2,2,3 ...., N, d a jrepresent when being 0 that jth root support cable is without health problem, d a jnumerical value represents when not being 0 that jth root support cable is the support cable of unsoundness problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, the degree of the lax or damage of its numerical response, cable system current actual damage vector d aelement numerical value be not less than 0, be not more than 100%, cable system current actual damage vector d aelement numerical value represent the degree of injury of corresponding support cable, if cable system current actual damage vector d athe numerical value of a certain element be 0, represent that the support cable corresponding to this element is intact, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that the support cable corresponding to this element is unsoundness problem, the health problem of this support cable may be impaired in the method also may be relax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the support cable of this its correspondence of element numeric representation with the current actual equivalent damage degree of its relax level mechanic equivalent,
M. identify damaged cable in the problematic support cable identified from l step, remaining is exactly slack line;
N. utilize at current cable structure steady temperature data vector T tthe cable system current actual damage vector d obtained in l step under condition athat obtain slack line with current actual equivalent damage degree that is its relax level mechanic equivalent, utilize e step obtain at current cable structure steady temperature data vector T tcurrent cable force vector F under condition and current support cable two support end points horizontal range vector, utilize c step obtain at initial Cable Structure steady temperature data vector T othe initial drift vector of the support cable under condition, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the support cable current steady state temperature data utilizing current cable structure steady temperature data vector Tt to represent, utilize c step obtain at initial Cable Structure steady temperature data vector T othe support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in c step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting 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 the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the non-destructive tests of support cable; During calculating, institute's demand power is provided by current cable force vector F corresponding element;
O. get back to e step, start the circulation next time being walked to o step by e.
Beneficial effect: when the temperature field of Cable Structure is subject to affecting of the factor such as sunshine and environment temperature, the temperature field of Cable Structure is constantly change, the change of temperature field of Cable Structure must affect the monitored amount of Cable Structure, only have and monitored amount is rejected could carry out rational cable structure health monitoring based on monitored amount by temperature profile effect part, and the temperature field measurement of Cable Structure and calculating are very complicated, this method discloses to comprise and is a kind ofly suitable for the simple of cable structure health monitoring, economical, feasible, the cable structure health monitoring method of efficient Cable Structure Calculation Method of Temperature Field, adopt this method when angular displacement appears in Cable Structure bearing, when many ropes of Cable Structure are synchronously impaired and lax, and the temperature of Cable Structure along with time variations time, health status that very monitor assessment identifies cable system (position of all slack lines and damaged cable can be comprised, and relax level or degree of injury), the effective health monitoring of system and method disclosed in this method to cable system is highly profitable.
Embodiment
When having angular displacement of support and temperature variation, for the health monitoring of the cable system of Cable Structure, this method discloses a kind of system and method can monitoring the health status identifying each root rope in cable system in Cable Structure rationally and effectively.The following describes of embodiment of this method is in fact only exemplary, and object is never the application or the use that limit this method.
This method adopts a kind of algorithm, and this algorithm is for monitoring the health status of the cable system in Cable Structure.During concrete enforcement, the following step is the one in the various steps that can take.
The first step: determine " the temperature survey calculating method of the Cable Structure of this method ", the method concrete steps are as follows:
A walks: inquiry or actual measurement (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure) obtain the temperature variant thermal conduction study parameter of environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model (such as finite element model) of Cable Structure.Inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of describing hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual observation record, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ".From the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, if the sea level elevation of such as Cable Structure is between 0m to 200m, so can choose height above sea level 0m, 50m, 100m and height above sea level 200m, crossing with Cable Structure surface with imaginary surface level at the sea level elevation place that each is chosen, obtain intersection, surface level is crossing with Cable Structure obtains cross surface, intersection is the outer edge line of cross surface, 6 points are chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness.In 6 directions of the measurement Cable Structure chosen along the Temperature Distribution of wall thickness, first according to the meteorological data throughout the year in region, Cable Structure position and the physical dimension of Cable Structure, volume coordinate, Cable Structure surrounding environment etc. determines the sunny slope of Cable Structure 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 the sea level elevation place that each is chosen, aforementioned intersection respectively has one section in sunny slope and in the shade, these two sections of intersection respectively have a mid point, cross the outer normal that these two mid points get Cable Structure, these two outer normals are called the sunny slope outer normal of Cable Structure and in the shade outer normal of Cable Structure by this method, these two outer normal directions are called the sunny slope outer normal direction of Cable Structure and in the shade outer normal direction of Cable Structure by this method, the outer normal of obvious sunny slope and the outer normal of in the shade all crossing with aforementioned intersection, also two intersection points are just had, intersection is divided into two line segments by these two intersection points, 2 points are got respectively on two line segments, totally 4 points, each line segment in two of intersection line segments is divided into equal 3 sections of length by taken point, the outer normal on Cable Structure surface is got at these 4 some places, the outer normal on 6 Cable Structure surfaces is just have chosen so altogether at each selected sea level elevation place, the direction of 6 outer normals is exactly " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ".There are two intersection points on the surface of each " measures the direction of Cable Structure along the Temperature Distribution of wall thickness " line and Cable Structure, if Cable Structure is hollow, these, two intersection points are on Cable Structure outside surface, another on an internal surface, if Cable Structure is solid, these two intersection points are all on Cable Structure outside surface, connect these two intersection points and obtain a straight-line segment, straight-line segment is chosen three points again, this straight-line segment is divided into four sections by these three points, measure two end points of Cable Structure at these three points chosen and straight-line segment, amount to the temperature of 5 points, concrete can first hole in Cable Structure, temperature sensor is embedded in this 5 some places, especially, can not hole in support cable, support cable is only measured to the temperature of support cable surface point, in any case, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, especially, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", the number temperature profile data at sea level elevation place " identical sea level elevation Cable Structure is along the temperature profile data of thickness " will chosen at each in this method ".Measuring temperature in Cable Structure location according to meteorology to require to choose a position, meeting obtaining in this position actual observation record the temperature that meteorology measures the Cable Structure place environment of temperature requirement, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable (as long as there was sunrise the same day, this position just should by sunlight), at the flat board (square plate that the wide 3mm of such as 30cm is thick) of this position of sound production one piece of carbon steel material (such as 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 be placed in and meet the top that meteorology temperature measures the wooden thermometer screen required, the one side of this reference plate on the sunny side, be called that sunny slope (such as, time on the Northern Hemisphere, sunny slope faces up towards south, full daytime is all by sunshine, sunny slope should have the suitable gradient to make snow can not accumulate or clear up sunny slope after snow), the sunny slope of reference plate is coarse with dark color (being conducive to accepting solar irradiation), the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, the non-sunny slope of reference plate is covered with insulation material (the thick calcium carbonate insulation material of such as 5mm), Real-Time Monitoring record is obtained the temperature of the sunny slope of reference plate.
B walks, Real-Time Monitoring (can be measured by ordinary temperature measuring method, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) record R the Cable Structure surface temperature measured data obtaining above-mentioned R Cable Structure surface point, Real-Time Monitoring (can be measured by ordinary temperature measuring method simultaneously, thermal resistance is such as used to measure, such as every 10 minutes survey records temperature data) obtain the temperature profile data of previously defined Cable Structure along thickness, Real-Time Monitoring (can be measured by ordinary temperature measuring method simultaneously, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record the temperature record obtaining the Cable Structure place environment meeting the requirement of meteorology measurement temperature, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, such as in the wooden thermometer screen meeting meteorology temperature measurement requirement, lay thermal resistance and measure temperature, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes the same day, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be designated as Δ T emax, (such as first the temperature measured data sequence of Cable Structure place environment is carried out curve fitting by Conventional mathematical calculating by the temperature measured data sequence of Cable Structure place environment, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of Cable Structure place environment, this rate of change is also along with time variations, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, such as use the temperature of the dull and stereotyped sunny slope of thermal resistance witness mark, such as every 10 minutes survey records temperature data) obtain being carved at sunrise the same day measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be designated as Δ T pmax, (can be measured by ordinary temperature measuring method by Real-Time Monitoring, thermal resistance is such as used to measure Cable Structure surface point, such as every 10 minutes survey records temperature data) record obtains being carved at sunrise the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes the same day, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein designated as Δ T smax, (such as first each Cable Structure surface temperature measured data sequence is carried out curve fitting by Conventional mathematical calculating by each Cable Structure surface temperature measured data sequence, then by asking curve to the derivative of time or by asking on curve each point corresponding to survey record data time by numerical method to the rate of change of time) obtain the rate of change of temperature about the time of each Cable Structure surface point, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations.Obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and 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 ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax.
C walks, and survey calculation obtains Cable Structure steady temperature data; First, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology; The a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, Δ T pmaxwith Δ T smaxall be not more than 5 degrees Celsius; The b condition that Section 2 must meet be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the Δ T that survey calculation obtains above emaxbe not more than with reference to temperature difference per day Δ T r, and the Δ T that survey calculation obtains above pmaxdeduct 2 degrees Celsius and be not more than Δ T emax, and the Δ T that survey calculation obtains above smaxbe not more than Δ T pmax; One of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition; Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 4 condition is when obtaining 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 0.1 degree Celsius per hour about the absolute value of the rate of change of time; Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise; Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius.This method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly one in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly moment of the acquisition Cable Structure steady temperature data of this method, then, according to Cable Structure heat transfer characteristic, utilize R Cable Structure surface temperature measured data and " HBE Cable Structure is along the thickness temperature measured data " in the moment obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model (such as finite element model) of Cable Structure, the Temperature Distribution that (such as finite element method) obtains the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ".When the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%; Cable Structure surface comprises support cable surface; Second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point "; " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation hthe numerical value obtained, gets Δ T for convenience of describing hunit be DEG C/m, be m for convenience of describing the unit getting Δ h; " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two; 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure.
Second step: set up initial mechanical Calculation Basis model A o.
If total N root support cable, first determines the coding rule of rope, numbered by ropes all in Cable Structure by this rule, this numbering will be used for generating vector sum matrix in subsequent step.Determine measured point (namely the specified point of all sign Cable Structure strain informations, is provided with K specified point), number to all specified points; Determine that the measured strain of each specified point (establishes the strain of L the assigned direction measuring each specified point, do not require that each specified point has the strain in the designated direction of same number, here the strain of L the assigned direction measuring each specified point is just established in order to describe conveniently), and give all measured strains numbering; Above-mentioned numbering will be used for equally generating vector sum matrix in subsequent step.Each specified point can be exactly a point near the fixed endpoint (being such as the stiff end of drag-line on bridge floor of cable-stayed bridge) of each root rope, and this point should not be generally stress concentration point, to avoid occurring excessive strain measurement value; This numbering will be used for equally generating vector sum matrix in subsequent step.Only can measure the strain in a direction at each specified point, also can measure the strain of multiple directions." the whole monitored strain data of Cable Structure " by K specified point in the Cable Structure determined above, the strain of L assigned direction of crossing each specified point describes, the change that Cable Structure strains is exactly the change of the strain of all assigned directions of all specified points, all appointment straight line.Each total individual strain measurement value of M (M=K × L) or calculated value characterize the strain information of Cable Structure.K and M must not be less than the quantity N of support cable.For simplicity, in the method by " the monitored strain data of Cable Structure " referred to as " monitored amount ".
When Cable Structure is completed, or before setting up health monitoring (damaged cable identification) system, obtain " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation (to measure by ordinary temperature measuring method, thermal resistance is such as used to measure), " Cable Structure steady temperature data " now use vector T orepresent, be called initial Cable Structure steady temperature data vector T o.T is obtained in actual measurement owhile, namely at the synchronization in the moment of the initial Cable Structure steady temperature data vector of acquisition, use the direct survey calculation of conventional method to obtain the initial value of all monitored amount of Cable Structure, form monitored amount initial value vector C o.
Specifically the synchronization in moment of so-and-so Cable Structure steady temperature data vector such as (such as initial or current) can be being obtained according to following method in this method, so-and-so method survey calculation is used to obtain the data of the monitored amount of so-and-so measured amount (all monitored amount of such as Cable Structure): (to comprise the temperature of Cable Structure place environment in survey record temperature, the temperature of the sunny slope of reference plate and Cable Structure surface temperature) while, such as every 10 minutes survey records temperature, so simultaneously equally also every 10 minutes the monitored amount of so-and-so measured amount of survey record (all monitored amount of such as Cable Structure) data.Once determine the moment obtaining Cable Structure steady temperature data, so just be called the synchronization in the moment obtaining Cable Structure steady temperature data, the data of the monitored amount of so-and-so measured amount using so-and-so method survey calculation method to obtain with the data of the monitored amount of so-and-so measured amount (all monitored amount of such as Cable Structure) of the moment synchronization of acquisition Cable Structure steady temperature data.
Conventional method (consult reference materials or survey) is used to obtain temperature variant physical parameter (such as thermal expansivity) and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement osynchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time cross-sectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T othe weight of the unit length of all support cables when the cross-sectional area of all support cables and Suo Li are 0 when the length of all support cables, Suo Li are 0 when Suo Li under condition is 0, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule 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 and Initial cable force vector F othe coding rule of element identical; Initial Cable Structure steady temperature data vector T is obtained at Actual measurement owhile, namely at the synchronization in the moment of acquisition Cable Structure steady temperature data, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Actual measurement data of Cable Structure comprise the measured data such as data, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing spatial data, initial Cable Structure bearing angular data, Cable Structure modal data, Cable Structure strain data, Cable Structure angular coordinate measurement data, Cable Structure volume coordinate measurement data that the Non-destructive Testing Data of support cable etc. can express the health status of rope.Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U o.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 in Cable Structure, and object is the geometric properties according to these coordinate data determination Cable Structure.For cable-stayed bridge, initial geometric data can be the spatial data that the spatial data of the end points of all ropes adds some points on bridge two ends, so-called bridge type data that Here it is.The data utilizing the Non-destructive Testing Data etc. of support cable can express the health status of rope set up cable system initial damage vector d o.If there is no the Non-destructive Testing Data of rope and other are when can express the data of the health status of support cable, or when can think that Cable Structure original state is 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 support cable corresponding to this element is intact, no problem, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the load-bearing capacity of corresponding proportion, if there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or when thinking that Cable Structure original state is not damaged state, vectorial d oeach element numerical value get 0; If d othe numerical value of a certain element be not 0, then represent that the support cable corresponding to this element is problematic, this support cable may be impaired in the method also may be lax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the initial equivalent damage degree of the support cable of this its correspondence of element numeric representation; Cable system initial damage vector d othe coding rule of element and Initial cable force vector F othe coding rule of element identical.The temperature variant physical and mechanical properties parameter of the various materials utilizing the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure to use, initial Cable Structure bearing angular coordinate vector U owith initial Cable Structure steady temperature data vector T o, utilize mechanics method (such as finite element method) to count " Cable Structure steady temperature data " and set up initial mechanical Calculation Basis model A o.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by o, count " Cable Structure steady temperature data " (i.e. initial Cable Structure steady temperature data vector T o), based on A othe Cable Structure that calculates calculates data must closely its measured data, and error generally must not be greater than 5%.Like this can utility A othe Suo Li calculated under the analog case of gained calculates data, strain calculation data, Cable Structure shapometer count certificate and displacement meter counts certificate, Cable Structure angle-data, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A othe health status cable system initial damage vector d of middle support cable orepresent, the initial Cable Structure steady temperature data vector T of Cable Structure steady temperature data orepresent.Due to based on A othe initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A obasis on, carry out Mechanics Calculation obtains, A othe evaluation of each monitored amount form monitored amount initial value vector C o.Corresponding to A o" Cable Structure steady temperature data " be exactly " initial Cable Structure steady temperature data vector T o"; Corresponding to A osupport cable health status with cable system initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent.Corresponding to A ocable Structure support coordinate data initial Cable Structure bearing angular coordinate vector U orepresent; T o, U oand d oa oparameter, C oby A omechanics Calculation result composition.
3rd step: first time sets up current initial mechanical Calculation Basis model A t o, monitored amount current initial value vector C t o" current initial Cable Structure steady temperature data vector T t o", concrete grammar is: at initial time, and namely first time sets up current initial mechanical Calculation Basis model A t owith monitored amount current initial value vector C t otime, A t ojust equal A o, C t ojust equal C o, A t ocorresponding " Cable Structure steady temperature data " are designated as " current initial Cable Structure steady temperature data vector T t o", at initial time, (namely first time sets up A t otime), T t ojust equal T o, vector T t odefinition mode and vector T odefinition mode identical.Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure t ocable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U t o; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t otime, U t ojust equal U o.A t othe health status of support cable and A osupport cable health status (cable system initial damage vector d orepresent) identical, A in cyclic process t othe health status of support cable use cable system initial damage vector d all the time orepresent.T t oand d oa t oparameter, C t oby A t omechanics Calculation result composition.
4th step: in Cable Structure military service process, the current data obtaining " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement (is called " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical).Current cable structure steady temperature data vector T is obtained in actual measurement twhile, namely at acquisition current cable structure steady temperature data vector T tthe synchronization in moment, actual measurement obtains the current measured value of all monitored amount of Cable Structure, composition " monitored amount current value vector C ".Current cable structure steady temperature data vector T is obtained in actual measurement twhile, actual measurement obtains Cable Structure bearing angular coordinate current data, all data composition current cable structure actual measurement bearing angular coordinate vector U t.Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the difference of the volume coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector othe coding rule of element identical.
5th step: according to current cable structure actual measurement bearing angular coordinate vector U twith current cable structure steady temperature data vector T t, upgrade current initial mechanical Calculation Basis model A where necessary t o, current initial Cable Structure bearing angular coordinate vector U t o, monitored amount current initial value vector C t owith current initial Cable Structure steady temperature data vector T t o.Current cable structure actual measurement bearing angular coordinate vector U is obtained in the 4th step actual measurement twith current cable structure steady temperature data vector T tafter, compare U respectively tand U t o, T tand T t oif, U tequal U t oand T tequal T t o, then do not need A t o, U t oand T t oupgrade, otherwise need A t o, U t oand T t oupgrade, update method is carried out to c step by following a step:
A step calculates U twith U odifference, U twith U odifference be exactly the front holder angular displacement of Cable Structure bearing about initial position, with angular displacement of support vector V represent angular displacement of support, V equals U tdeduct U o, be one-to-one relationship between the element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing.
B step calculates T twith T odifference, T twith T odifference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T odifference represent with steady temperature change vector S, S equals T tdeduct T o, S represents the change of Cable Structure steady temperature data.
C step is first to A oin Cable Structure bearing apply front holder angular displacement constraint, the numerical value of front holder angular displacement constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A omiddle Cable Structure bearing applies angular displacement of support constraint and to A oin Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A t o, upgrade A t owhile, U t oall elements numerical value also uses U tall elements numerical value correspondence replaces, and namely have updated U t o, T t oall elements numerical value also uses T tall elements numerical value correspondence replace, namely have updated T t o, so just obtain and correctly correspond to A t ot t oand U t o; Upgrade C t omethod be: when renewal A t oafter, obtain A by Mechanics Calculation t oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C t o.
6th step: at current initial mechanical Calculation Basis model A t obasis on carry out several times Mechanics Calculation, obtain Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D by calculating u.Concrete grammar is: Cable Structure unit damage monitored amount transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal t owith current cable structural bearings angular coordinate vector U t owhile, Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D must be upgraded simultaneously u; At the current initial mechanical Calculation Basis model A of Cable Structure t obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity of all ropes, has N root rope just to have N calculating, each time calculate hypothesis cable system in only have a rope to have unit damage D u(such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), occur in calculating each time that the rope damaged is different from during other time calculates the rope occurring damaging, calculate the current calculated value of all monitored amounts in Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector C; Calculate monitored amount calculation current vector C each time and deduct monitored amount current initial value vector C t oobtain a monitored amount change vector; N root rope is had just to have N number of monitored amount change vector; The unit damage monitored amount transformation matrices Δ C having N to arrange is made up of successively this N number of monitored amount change vector; Each row of unit damage monitored amount transformation matrices correspond to a monitored amount change vector.
7th step: set up linear relationship error vector e and vectorial g.Utilize (the monitored amount current initial value vector C of data above t o, unit damage monitored amount transformation matrices Δ C), while the 6th step calculates each time, namely calculate each time hypothesis cable system in only have a rope to have unit damage D uoccur in calculating each time that the rope damaged is different from during other time calculates the rope occurring damaging, calculate the current value all utilizing mechanics method (such as adopting finite element method) to calculate all monitored amounts in cable system in Cable Structure each time, while calculating the monitored amount calculation current vector C of composition one each time, calculate composition injury vector d each time, this is walked out of existing injury vector d and only uses in this step, only has the numerical value of an element to get D in all elements of this injury vector d u, the numerical value of other element gets 0, and in injury vector d, numerical value is D uelement correspond to this time calculate time unique damaged cable unit damage degree D u; By C, C t o, Δ C, D u, d brings formula (12) into, obtain a linear relationship error vector e, calculate a linear relationship error vector e each time; Having N root rope just to have N calculating, just have N number of linear relationship error vector e, obtain a vector after being added by this N number of linear relationship error vector e, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N.Vector g equals final error vector e.
8th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system is (such as containing strain measurement system, signal conditioner etc.), Cable Structure bearing angular coordinate monitoring system is (containing angle measuring sensor, signal conditioner etc.), Cable Structure temperature monitoring system is (containing temperature sensor, signal conditioner etc.) and Cable Structure ambient temperature measurement system (containing temperature sensor, signal conditioner etc.), cable force monitoring system is (such as containing acceleration transducer, signal conditioner etc.), each support cable two supports the horizontal range monitoring system (such as monitoring with total powerstation) of end points, signal (data) collector, computing machine and communication alert equipment.The bearing angular coordinate of each monitored amount, each Cable Structure, the Suo Li of each root support cable, each root support cable two support the horizontal range of end points, each temperature must arrive by monitored system monitoring, monitoring system by the Signal transmissions that monitors to signal (data) collector; Signal is delivered to computing machine through signal picker; The health monitoring software of the cable system running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring rope and having damage, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.
9th step: by current for monitored amount initial value vector C t o, unit damage monitored amount transformation matrices Δ C, unit damage scalar D uparameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
Tenth step: establishment and on computers installation and operation angular displacement of support temperature variation time based on the slack line recognition methods system software of strain monitoring, the function (in this specific implementation method all work that can complete with computing machine) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method " based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation " required by task and wants
11 step: according to monitored amount current value vector C with monitored amount current initial value vector C t o, unit damage monitored amount transformation matrices Δ C, unit damage scalar D uand cable system current nominal fatigue vector d(be made up of all Suo Dangqian nominal fatigue amounts) between exist linear approximate relationship (formula (8)), calculate the noninferior solution of cable system current nominal fatigue vector d according to multi-objective optimization algorithm, namely can determine the position of damaged cable and the solution of nominal fatigue degree thereof more exactly with reasonable error from all ropes.
The multi-objective optimization algorithm that can adopt has a variety of, such as: the multiple-objection optimization based on genetic algorithm, the multiple-objection optimization based on artificial neural network, the multi-objective optimization algorithm based on population, the multiple-objection optimization based on ant group algorithm, leash law (Constrain Method), weighted method (Weighted SUm Method), Objective Programming (Goal Attainment Method) etc.Because various multi-objective optimization algorithm is all conventional algorithm, can realize easily, this implementation step only provides the process solving current injury vector d for Objective Programming, 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 (8) can transform the multi-objective optimization question shown in an accepted way of doing sth (27) and formula (28), in formula (27), γ is a real number, R is real number field, area of space Ω limits the span (each element of the present embodiment requirements vector d is not less than 0, is not more than 1) of each element of vectorial d.Formula (27) be meant to the minimum real number γ of searching one, formula (28) is met.In formula (28), G (d) is defined by formula (29), and the middle deviation allowed between G (d) and vectorial g of the product representation formula (28) of weighing vector W and γ in formula (28), the definition of g is see formula (13), and its value calculates in the 7th step.During actual computation, vector W can be identical with vectorial g.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of cable system current nominal fatigue vector d.
min imize γ γ ∈ R , d ∈ Ω - - - ( 27 )
G(d)-Wγ≤g (28)
G ( d ) = abs ( 1 D u ΔC · d - C + C o t ) - - - ( 29 )
The element number of cable system current nominal fatigue vector d equals the quantity of rope, be one-to-one relationship between the element of cable system current nominal fatigue vector d and rope, the element numerical value of cable system current nominal fatigue vector d represents the nominal fatigue degree of corresponding rope or nominal health status.Vector d the coding rule of element and vectorial d othe coding rule of element identical.
12 step: definition cable system current actual damage vector d a, cable system current actual damage vector d aelement number equal the quantity of support cable, cable system current actual damage vector d aelement and support cable between be one-to-one relationship, cable system current actual damage vector d aelement numerical value represent the actual damage degree of corresponding support cable or actual health status; Vector d athe coding rule of element and vectorial d othe coding rule of element identical.The cable system utilizing formula (15) to express current actual damage vector d aa jth element d a jwith cable system initial damage vector d oa jth element d ojwith a jth element d of cable system current nominal fatigue vector d jbetween relation, calculate cable system current actual damage vector d aall elements; d a jrepresent when being 0 that jth root support cable is without health problem, d a jnumerical value represents when not being 0 that jth root support cable is the support cable of unsoundness problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, the degree of the lax or damage of its numerical response.Cable system current actual damage vector d aelement numerical value be not less than 0, be not more than 100%, cable system current actual damage vector d aelement numerical value represent the degree of injury of corresponding support cable, if cable system current actual damage vector d athe numerical value of a certain element be 0, represent that the support cable corresponding to this element is intact, without health problem, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that the support cable corresponding to this element is unsoundness problem, the health problem of this support cable may be impaired in the method also may be relax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the support cable of this its correspondence of element numeric representation with the current actual equivalent damage degree of its relax level mechanic equivalent.
13 step: identify damaged cable in the problematic support cable identified from the 12 step, remaining is exactly slack line.Mirror method for distinguishing is varied; can by removing the protective seam of the support cable of unsoundness problem; visual discriminating is carried out to support cable; or carry out visual discriminating by optical imaging apparatus; also can be differentiated whether support cable is impaired by lossless detection method, UT (Ultrasonic Testing) is exactly a kind of now widely used lossless detection method.After differentiating, those do not find that the support cable of the unsoundness problem damaged is exactly there occurs lax rope, need the rope adjusting Suo Li exactly.
14 step: utilize at current cable structure steady temperature data vector T tthe cable system current actual damage vector d obtained in the 12 step under condition athat obtain slack line with current actual equivalent damage degree that is its relax level mechanic equivalent, utilize the 4th step obtain at current cable structure steady temperature data vector T tcurrent cable force vector F under condition and current support cable two support end points horizontal range vector, utilize second step obtain at initial Cable Structure steady temperature data vector T othe initial drift vector of the support cable under 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 trepresent support cable current steady state temperature data, utilize second step obtain at initial Cable Structure steady temperature data vector T othe support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained at second step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical.When meeting 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 the slack line of equivalence, any change can not occur Cable Structure, and vice versa.Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li.Particularly can in the hope of the relax level of these ropes (i.e. the long adjustment amount of rope) according to formula (25) or formula (26).So just achieve lax identification and the non-destructive tests of support cable.During calculating, institute's demand power is provided by current cable force vector F corresponding element.
15 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.Under specified requirements, the computing machine automatic operation communication alert equipment in health monitoring systems is reported to the police to monitor staff, owner and (or) the personnel that specify.
16 step: get back to the 4th step, starts by the circulation of the 4th step to the 16 step.

Claims (1)

1. during angular displacement of support temperature variation based on a slack line recognition methods for strain monitoring, it is characterized in that described method comprises:
A. establish total N root support cable, first determine the coding rule of support cable, support cables all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; Determine the point being monitored of specifying, namely point being monitored characterizes all specified points of Cable Structure strain information, and gives all specified points numbering; Determine monitored should the changing direction of point being monitored, and give all monitored strain numberings of specifying; " monitored strain numbering " will be used for generating vector sum matrix in subsequent step; " the whole monitored strain data of Cable Structure " is made up of above-mentioned all monitored strains; This method by " the monitored strain data of Cable Structure " referred to as " monitored amount "; The quantity of point being monitored 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 moment;
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 environment residing for Cable Structure composition material and Cable Structure, utilize the geometry measured data of the design drawing of Cable Structure, as-constructed drawing and Cable Structure, utilize these data and parameter to set up the Thermodynamic calculation model of Cable Structure, inquiry Cable Structure location is no less than the meteorological data in recent years of 2 years, the statistics cloudy quantity obtained is during this period of time designated as T cloudy day, in the method daytime can not be seen one of the sun and be called the cloudy day all day, statistics obtains 0 highest temperature after sunrise moment next day between 30 minutes and the lowest temperature at each cloudy day in T cloudy day, the sunrise moment refers to the sunrise moment on the meteorology that base area revolutions and revolution rule are determined, do not represent that the same day necessarily can see the sun, data can be inquired about or calculated sunrise moment of each required day by conventional meteorology, 0 highest temperature after sunrise moment next day between 30 minutes at each cloudy day deducts the maximum temperature difference that the lowest temperature is called the daily temperature at this cloudy day, there is T cloudy day, just there is the maximum temperature difference of the daily temperature at T cloudy day, the maximal value of getting in the maximum temperature difference of the daily temperature at T cloudy day is reference temperature difference per day, Δ T is designated as with reference to temperature difference per day r, be no less than between inquiry Cable Structure location and Altitude Region, place temperature that the meteorological data in recent years of 2 years or actual measurement obtain environment residing for Cable Structure in time with delta data and the Changing Pattern of sea level elevation, calculate to be no less than 2 years between Cable Structure location and Altitude Region, place Cable Structure in recent years residing for the temperature of environment about the maximum rate of change Δ T of sea level elevation h, get Δ T for convenience of describing hunit be DEG C/m, the surface of Cable Structure is got " R Cable Structure surface point ", the Specific Principles getting " R Cable Structure surface point " describes in step b3, the temperature of this R Cable Structure surface point will be obtained below by actual measurement, claiming to survey the temperature data obtained is " R Cable Structure surface temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtained the temperature of this R Cable Structure surface point by Calculation of Heat Transfer, just claim the temperature data that calculates for " R Cable Structure land surface pyrometer count certificate ", from the minimum height above sea level residing for Cable Structure to most High aititude, in Cable Structure, uniform choosing is no less than three different sea level elevations, at the sea level elevation place that each is chosen, two points are at least chosen at the intersection place on surface level and Cable Structure surface, from the outer normal of selected point straw line body structure surface, all outer normal directions chosen are called " measuring the direction of Cable Structure along the Temperature Distribution of wall thickness ", measure Cable Structure crossing with " intersection on surface level and Cable Structure surface " along the direction of the Temperature Distribution of wall thickness, in in the shade the outer normal direction of the measurement Cable Structure chosen along the sunny slope outer normal direction and Cable Structure that must comprise Cable Structure in the direction of the Temperature Distribution of wall thickness, three points are no less than along each measurement Cable Structure along direction uniform choosing in Cable Structure of the Temperature Distribution of wall thickness, along each, Cable Structure is measured for support cable and only gets a point along the direction of the Temperature Distribution of wall thickness, only measure the temperature of the surface point of support cable, measure all temperature be selected a little, the temperature recorded 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 the direction of Cable Structure along the Temperature Distribution of wall thickness " to 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 have chosen the individual different sea level elevation of H, at each sea level elevation place, have chosen B and measure the direction of Cable Structure along the Temperature Distribution of wall thickness, in Cable Structure, E point is have chosen along each measurement Cable Structure along the direction of the Temperature Distribution of wall thickness, wherein H and E is not less than 3, B is not less than 2, 1 is equaled for support cable E, what meter Cable Structure " measured the point of Cable Structure along the temperature profile data of thickness " adds up to HBE, the temperature of this HBE " measuring the point of Cable Structure along the temperature profile data of thickness " will be obtained below by actual measurement, claiming to survey the temperature data obtained is " HBE Cable Structure is along thickness temperature measured data ", if utilize the Thermodynamic calculation model of Cable Structure, obtain this HBE by Calculation of Heat Transfer and measure the temperature of Cable Structure along the point of the temperature profile data of thickness, the temperature data calculated just is claimed to be " HBE Cable Structure calculates data along thickness temperature ", measure temperature in Cable Structure location according to meteorology to require to choose a position, obtain meeting the temperature that meteorology measures the Cable Structure place environment of temperature requirement by the actual measurement of this position, in the on-site spaciousness of Cable Structure, unobstructed place chooses a position, this position should each of the whole year day can obtain this ground the most sufficient sunshine of this day getable, at the flat board of this position of sound production one piece of carbon steel material, 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 with dark color, the sunny slope of reference plate should each of the whole year day can obtain one flat plate on this ground the most sufficient sunshine of this day getable, 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 the temperature profile data of previously defined Cable Structure along thickness simultaneously, and Real-Time Monitoring obtains meeting the temperature record that meteorology measures the Cable Structure place environment of temperature requirement simultaneously, the temperature measured data sequence of the Cable Structure place environment after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the temperature measured data sequence of Cable Structure place environment is arranged according to time order and function order by the temperature measured data of the Cable Structure place environment after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the temperature measured data sequence of Cable Structure place environment and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains Cable Structure place environment after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the temperature measured data sequence of Cable Structure place environment, be called environment maximum temperature difference, be designated as Δ T emax, calculated the rate of change of temperature about the time of Cable Structure place environment by Conventional mathematical by the temperature measured data sequence of Cable Structure place environment, this rate of change is also along with time variations, the measured data sequence of the temperature of the sunny slope of the reference plate after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, the measured data sequence of the temperature of the sunny slope of reference plate is arranged according to time order and function order by the measured data of the temperature of the sunny slope of the reference plate after being carved into the sunrise moment next day same day at sunrise between 30 minutes, find the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of the sunny slope of reference plate after sunrise moment next day between 30 minutes is deducted by the maximum temperature in the measured data sequence of the temperature of the sunny slope of reference plate, be called reference plate maximum temperature difference, be designated as Δ T pmax, the Cable Structure surface temperature measured data sequence of all R Cable Structure surface points after sunrise moment next day between 30 minutes is obtained being carved at sunrise the same day by Real-Time Monitoring, R Cable Structure surface point is had just to have R Cable Structure surface temperature measured data sequence, each Cable Structure surface temperature measured data sequence is arranged according to time order and function order by the Cable Structure surface temperature measured data after being carved into the sunrise moment next day same day of a Cable Structure surface point at sunrise between 30 minutes, find the maximum temperature in each Cable Structure surface temperature measured data sequence and minimum temperature, the maximum temperature difference be carved at sunrise on same day that minimum temperature obtains the temperature of each Cable Structure surface point after sunrise moment next day between 30 minutes is deducted by the maximum temperature in each Cable Structure surface temperature measured data sequence, there is R Cable Structure surface point just to have to be carved at sunrise R the same day maximum temperature difference numerical value between 30 minutes after sunrise moment next day, maximal value is wherein called Cable Structure surface maximum temperature difference, be designated as Δ T smax, calculated the rate of change of temperature about the time of each Cable Structure surface point by Conventional mathematical by each Cable Structure surface temperature measured data sequence, the temperature of each Cable Structure surface point about the rate of change of time also along with time variations, obtain being carved at sunrise the same day after sunrise moment next day between 30 minutes by Real-Time Monitoring, at synchronization, after HBE " Cable Structure is along the temperature profile data of thickness ", calculate the difference amounting to maximum temperature in BE " identical sea level elevation Cable Structure is along the temperature profile data of thickness " and 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 ", have chosen H different sea level elevation just to have H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference ", the maximal value in this H " identical sea level elevation place Cable Structure thickness direction maximum temperature difference " is claimed to be " Cable Structure thickness direction maximum temperature difference ", be designated as Δ T tmax,
B3: survey calculation obtains Cable Structure steady temperature data, first, determine the moment obtaining Cable Structure steady temperature data, the condition relevant to the moment determining to obtain Cable Structure steady temperature data has six, Section 1 condition is moment of obtaining Cable Structure steady temperature data between after being carved into sunrise moment next day at sunset between 30 minutes on same day, the sunset moment refers to the sunset moment on base area revolutions and the meteorology determined of revolution rule, can inquire about data or be calculated sunset moment of each required day by conventional meteorology, the a condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day 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 Celsius, the b condition of Section 2 condition be after being carved into sunrise moment next day at sunrise between 30 minutes on same day during this period of time in, the environment maximum temperature difference Δ T that survey calculation obtains above emaxbe not more than with reference to temperature difference per day Δ T r, and reference plate maximum temperature difference Δ T pmaxΔ T is not more than after deducting 2 degrees Celsius emax, and Cable Structure surface maximum temperature difference Δ T smaxbe not more than Δ T pmax, one of only need meet in a condition of Section 2 and b condition is just called and meets Section 2 condition, Section 3 condition is when obtaining Cable Structure steady temperature data, and the temperature of Cable Structure place environment is not more than 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 4 condition is when obtaining 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 0.1 degree Celsius per hour about the absolute value of the rate of change of time, Section 5 condition is when obtaining 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 be carved into the minimal value after sunrise moment next day between 30 minutes the same day at sunrise, Section 6 condition is when obtaining Cable Structure steady temperature data, " Cable Structure thickness direction maximum temperature difference " Δ T tmaxbe not more than 1 degree Celsius, this method utilizes above-mentioned six conditions, any one in following three kinds of moment is called " the mathematics moment obtaining Cable Structure steady temperature data ", the first moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 5 condition, the second moment is moment of the Section 6 condition only met in above-mentioned " condition relevant to determining to obtain moment of Cable Structure steady temperature data ", simultaneously the third moment meets the moment of the Section 1 in above-mentioned " condition relevant to the moment determining to obtain Cable Structure steady temperature data " to Section 6 condition, when the mathematics moment obtaining Cable Structure steady temperature data is exactly a moment in this method in the physical record data moment, the moment obtaining Cable Structure steady temperature data is exactly the mathematics moment obtaining Cable Structure steady temperature data, if the mathematics moment obtaining Cable Structure steady temperature data is not any one moment in this method in the physical record data moment, then getting this method closest to moment of those physical record data in the mathematics moment obtaining Cable Structure steady temperature data is the moment obtaining Cable Structure steady temperature data, the amount being used in the moment survey record obtaining Cable Structure steady temperature data is carried out Cable Structure relevant health monitoring analysis by this method, this method is similar to thinks that the Cable Structure temperature field in the moment obtaining Cable Structure steady temperature data is in stable state, and namely the Cable Structure temperature in this moment does not change in time, and this moment is exactly " obtaining the moment of Cable Structure steady temperature data " of this method, then, according to 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 obtaining Cable Structure steady temperature data, utilize the Thermodynamic calculation model of Cable Structure, the Temperature Distribution of the Cable Structure in the moment obtaining Cable Structure steady temperature data is calculated by conventional heat transfer, now the temperature field of Cable Structure calculates by stable state, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated comprises the accounting temperature of R Cable Structure surface point in Cable Structure, the accounting temperature of R Cable Structure surface point is called that R Cable Structure steady-state surface temperature calculates data, also comprise the accounting temperature of Cable Structure selected HBE " measuring the point of Cable Structure along the temperature profile data of thickness " above, the accounting temperature of HBE " measuring the point of Cable Structure along the temperature profile data of thickness " is called " HBE Cable Structure calculates data along thickness temperature ", when R Cable Structure surface temperature measured data and R Cable Structure steady-state surface temperature calculate data correspondent equal, and when " HBE Cable Structure is along thickness temperature measured data " and " HBE Cable Structure calculates data along thickness temperature " correspondent equal, the temperature profile data of the Cable Structure in the moment in acquisition Cable Structure steady temperature data calculated is called " Cable Structure steady temperature data " in the method, " R Cable Structure surface temperature measured data " is now called " R Cable Structure steady-state surface temperature measured data ", " HBE Cable Structure is along thickness temperature measured data " is called " HBE Cable Structure is along thickness steady temperature measured data ", when the surface of Cable Structure is got " R Cable Structure surface point ", the quantity of " R Cable Structure surface point " must meet three conditions with distribution, first condition is when Cable Structure temperature field is in stable state, when the temperature of Cable Structure any point be on the surface by " R Cable Structure surface point " in obtain with the observed temperature linear interpolation of the Cable Structure point that this arbitrfary point is adjacent on the surface time, the error of the Cable Structure that linear interpolation the obtains temperature of this arbitrfary point and the Cable Structure actual temperature of this arbitrfary point on the surface is on the surface not more than 5%, Cable Structure surface comprises support cable surface, second condition is not less than 4 in the quantity of the point of same sea level elevation in " R Cable Structure surface point ", and uniform along Cable Structure surface at the point of same sea level elevation in " R Cable Structure surface point ", " R Cable Structure surface point " is not more than 0.2 DEG C divided by Δ T along the maximal value Δ h in the absolute value of all differences of the sea level elevation of adjacent Cable Structure surface point between two of sea level elevation hthe numerical value obtained, gets Δ T for convenience of describing hunit be DEG C/m, be m for convenience of describing the unit getting Δ h, " R Cable Structure surface point " along the definition of the adjacent between two Cable Structure surface point of sea level elevation refer to only consider sea level elevation time, in " R Cable Structure surface point ", there is not a Cable Structure surface point, the sea level elevation numerical value of this Cable Structure surface point is between the sea level elevation numerical value of adjacent Cable Structure surface point between two, 3rd condition is inquiry or obtains the rule at sunshine between Cable Structure location and Altitude Region, place by meteorology conventionally calculation, again according to geometric properties and the bearing data of Cable Structure, Cable Structure finds the annual position by sunshine-duration those surface points the most sufficient, in " R Cable Structure surface point ", has at least a Cable Structure surface point to be annual by a point in sunshine-duration the most fully those surface points in Cable Structure,
C. the Cable Structure steady temperature data under original state are obtained according to " the temperature survey calculating method of the Cable Structure of this method " direct survey calculation, Cable Structure steady temperature data under 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 various materials that Cable Structure uses; Initial Cable Structure steady temperature data vector T is obtained in actual measurement osynchronization, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; Obtain the length of all support cables when free state and Suo Li are 0 according to Cable Structure design data, completion data, in free state time cross-sectional area and in free state time the weight of unit length, and the temperature of all support cables when obtaining these three kinds of data, utilize temperature variant physical function parameter and the mechanical property parameters of all support cables on this basis, conveniently physical computing obtains all support cables at initial Cable Structure steady temperature data vector T othe weight of the unit length of all support cables when the cross-sectional area of all support cables and Suo Li are 0 when the length of all support cables, Suo Li are 0 when Suo Li under condition is 0, form the initial drift vector of support cable, the weight vector of the initial free unit length of initial free cross-sectional area vector sum successively, the coding rule 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 and Initial cable force vector F othe coding rule of element identical; T is obtained in actual measurement owhile, namely at the initial Cable Structure steady temperature data vector T of acquisition othe synchronization in moment, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises the Non-destructive Testing Data of the health status expressing support cable, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing spatial data, initial Cable Structure bearing angular data, initial Cable Structure angle-data, initial Cable Structure spatial data; The initial value of all monitored amounts forms monitored amount initial value vector C o; The Non-destructive Testing Data that utilization can express the health status of support cable sets 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 support cable corresponding to this element is intact, no problem, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the load-bearing capacity of corresponding proportion, if there is no the Non-destructive Testing Data of support cable and other are when can express the data of the health status of support cable, or when thinking that Cable Structure original state is not damaged state, vectorial d oeach element numerical value get 0; If d othe numerical value of a certain element be not 0, then represent that the support cable corresponding to this element is problematic, this support cable may be impaired in the method also may be lax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the initial equivalent damage degree of the support cable of this its correspondence of element numeric representation; Cable system initial damage vector d othe coding rule of element and Initial cable force vector F othe coding rule of element identical; Initial Cable Structure bearing angular data forms initial Cable Structure bearing angular coordinate vector U o;
Temperature variant physical and mechanical properties parameter, the initial Cable Structure bearing angular coordinate vector U of the various materials d. used according to the measured data of the design drawing of Cable Structure, as-constructed drawing and initial Cable Structure, the Non-destructive Testing Data of support cable, Cable Structure o, initial Cable Structure steady temperature data vector T owith all Cable Structure data that preceding step obtains, set up the initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " o, based on A othe Cable Structure that calculates calculates data must closely 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 angular data be exactly initial Cable Structure bearing angular coordinate vector U o; Corresponding to A osupport cable health status with cable system initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent; First time sets up the current initial mechanical Calculation Basis model A counting the Cable Structure of " Cable Structure steady temperature data " t o, monitored amount current initial value vector C t o" current initial Cable Structure steady temperature data vector T t o"; Set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t owith monitored amount current initial value vector C t otime, the current initial mechanical Calculation Basis model A of Cable Structure t ojust equal the initial mechanical Calculation Basis model A of Cable Structure o, monitored amount current initial value vector C t ojust equal monitored amount initial value vector C o; A t ocorresponding " Cable Structure steady temperature data " are called " current initial Cable Structure steady temperature data ", are designated as " current initial Cable Structure steady temperature data vector T t o", set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t otime, T t ojust equal T o; Corresponding to the current initial mechanical Calculation Basis model A of Cable Structure t ocable Structure bearing angular data composition current initial Cable Structure bearing angular coordinate vector U t o, set up the current initial mechanical Calculation Basis model A of Cable Structure for the first time t otime, U t ojust equal U o; A t othe initial health of support cable and A othe health status of support cable identical, also use cable system initial damage vector d orepresent, A in cyclic process below t othe initial health of support cable use cable system initial damage vector d all the time orepresent; Work as T o, U oand d oa oparameter time, by A othe initial value of all monitored amount that obtains of Mechanics Calculation result and C othe initial value of all monitored amount represented is identical, therefore alternatively C oby A omechanics Calculation result composition, work as T t o, U t oand d oa t oparameter time, C t oby A t omechanics Calculation result composition; A in the method o, U o, C o, d oand T oconstant;
E. from entering the circulation being walked to o step by e here; In Cable Structure military service process, the current data of " Cable Structure steady temperature data " is constantly obtained according to " the temperature survey calculating method of the Cable Structure of this method " continuous Actual measurement, the current data of " Cable Structure steady temperature data " is called " current cable structure steady temperature data ", is designated as " current cable structure steady temperature data vector T t", vector T tdefinition mode and vector T odefinition mode identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains the rope force data of all support cables in Cable Structure, all these rope force datas composition current cable force vector F, the element of vectorial F and vectorial F othe coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, Actual measurement obtains the volume coordinate of two supporting end points of all support cables, the difference of the volume coordinate component in the horizontal direction of two supporting end points is exactly two supporting end points horizontal ranges, two supporting end points horizontal range data of all support cables form current support cable two and support end points horizontal range vector, and current support cable two supports coding rule and the Initial cable force vector F of the element of end points horizontal range vector othe coding rule of element identical; Current cable structure steady temperature data vector T is obtained in actual measurement tsynchronization, actual measurement obtains Cable Structure bearing angular coordinate current data, all Cable Structure bearing angular coordinate current datas composition current cable structure actual measurement bearing angular coordinate vector U t;
F. according to current cable structure actual measurement bearing angular coordinate vector U twith current cable structure steady temperature data vector T t, upgrade current initial mechanical Calculation Basis model A according to step f1 to f3 t o, current initial Cable Structure bearing angular coordinate vector U t o, monitored amount current initial value vector C t owith current initial Cable Structure steady temperature data vector T t o;
F1. U is compared respectively twith U t o, T twith T t oif, U tequal U t oand T tequal T t o, then A t o, U t o, C t oand T t oremain unchanged; Otherwise need to follow these steps to A t o, U t o, C t oand T t oupgrade;
F2. U is calculated twith U odifference, U twith U odifference be exactly the front holder angular displacement of Cable Structure bearing about initial position, with angular displacement of support vector V represent angular displacement of support, V equals U tdeduct U o, be one-to-one relationship between the element in angular displacement of support vector V and angular displacement of support component, in angular displacement of support vector V, the numerical value of an element corresponds to the angular displacement of an assigned direction of an appointment bearing; Calculate T twith T odifference, T twith T odifference be exactly the changes of current cable structure steady temperature data about initial Cable Structure steady temperature data, T twith T odifference represent with steady temperature change vector S, S equals T tdeduct T o, S represents the change of Cable Structure steady temperature data;
F3. first to A oin Cable Structure bearing apply front holder angular displacement constraint, the numerical value of front holder angular displacement constraint just takes from the numerical value of corresponding element in angular displacement of support vector V, then to A oin Cable Structure apply temperature variation, the numerical value of the temperature variation of applying just takes from steady temperature change vector S, to A omiddle Cable Structure bearing applies angular displacement of support constraint and to A oin Cable Structure apply temperature variation after obtain upgrade current initial mechanical Calculation Basis model A t o, upgrade A t owhile, U t oall elements numerical value also uses U tall elements numerical value correspondence replaces, and namely have updated U t o, T t oall elements numerical value also uses T tall elements numerical value correspondence replace, namely have updated T t o, so just obtain and correctly correspond to A t ot t oand U t o; Upgrade C t omethod be: when renewal A t oafter, obtain A by Mechanics Calculation t oin all monitored amounts, current concrete numerical value, these concrete numerical value composition C t o; A t othe initial health of support cable use cable system initial damage vector d all the time orepresent;
G. at current initial mechanical Calculation Basis model A t obasis on carry out several times Mechanics Calculation according to step g 1 to g4, obtain Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D by calculating u;
G1. Cable Structure unit damage monitored amount transformation matrices Δ C constantly updates, namely at the current initial mechanical Calculation Basis model A of renewal t o, current initial Cable Structure bearing angular coordinate vector U t o, monitored amount current initial value vector C t owith current initial Cable Structure steady temperature data vector T t oafterwards, Cable Structure unit damage monitored amount transformation matrices Δ C and unit damage scalar D must then be upgraded u;
G2. at the current initial mechanical Calculation Basis model A of Cable Structure t obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity of all ropes, has N root support cable just to have N calculating, each time calculate hypothesis cable system in only have a support cable to have unit damage scalar D uoccur in calculating each time that the rope damaged is different from during other time calculates the rope occurring damaging, calculate the current calculated value of all monitored amounts in Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector oelement number rule identical;
G3. the monitored amount calculation current vector calculated each time deducts monitored amount current initial value vector C t oobtain a monitored amount change vector; N root support cable is had just to have N number of monitored amount change vector;
G4. the Cable Structure unit damage monitored amount transformation matrices Δ C having N to arrange is made up of successively this N number of monitored amount change vector; Each row of Cable Structure unit damage monitored amount transformation matrices Δ C correspond to a monitored amount change vector;
H. current cable structure steady temperature data vector T is obtained in actual measurement twhile, actual measurement obtains at acquisition current cable structure steady temperature data vector T tthe current measured value of all monitored amount of Cable Structure of synchronization in moment, form monitored amount current value vector C; Monitored amount current value vector C and monitored amount current initial value vector C t owith monitored amount initial value vector C odefinition mode identical, the same monitored amount of element representation of three vectorial identical numberings is at not concrete numerical value in the same time;
I. cable system current nominal fatigue vector d is defined, the element number of cable system current nominal fatigue vector d equals the quantity of support cable, be one-to-one relationship between the element of cable system current nominal fatigue vector d and support cable, the element numerical value of cable system current nominal fatigue vector d represents the nominal fatigue degree of corresponding support cable or nominal health status; The coding rule of the element of vector d and vectorial d othe coding rule of element identical;
J. according to monitored amount current value vector C with monitored amount current initial value vector C t o, Cable Structure unit damage monitored amount transformation matrices Δ C, unit damage scalar D uand the linear approximate relationship existed between cable system to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, and other amount in formula 1 except d is known, solves formula 1 and just can calculate cable system current nominal fatigue vector d;
C = C o t + 1 D u ΔC · d Formula 1
K. cable system current actual damage vector d is defined a, cable system current actual damage vector d aelement number equal the quantity of support cable, cable system current actual damage vector d aelement and support cable between be one-to-one relationship, cable system current actual damage vector d aelement numerical value represent the actual damage degree of corresponding support cable or actual health status; Vector d athe coding rule of element and vectorial d othe coding rule of element identical;
L. the cable system utilizing formula 2 to express current actual damage vector d aa jth element d a jwith cable system initial damage vector d oa jth element d ojwith a jth element d of cable system current nominal fatigue vector d jbetween relation, calculate cable system current actual damage vector d aall elements;
d j a = 1 - ( 1 - d oj ) ( 1 - d j ) Formula 2
J=1 in formula 2,2,3 ...., N, d a jrepresent when being 0 that jth root support cable is without health problem, d a jnumerical value represents when not being 0 that jth root support cable is the support cable of unsoundness problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, the degree of the lax or damage of its numerical response, cable system current actual damage vector d aelement numerical value be not less than 0, be not more than 100%, cable system current actual damage vector d aelement numerical value represent the degree of injury of corresponding support cable, if cable system current actual damage vector d athe numerical value of a certain element be 0, represent that the support cable corresponding to this element is intact, if its numerical value is 100%, then represent that the support cable corresponding to this element has completely lost load-bearing capacity, if its numerical value is between 0 and 100%, then represent that the support cable corresponding to this element is unsoundness problem, the health problem of this support cable may be impaired in the method also may be relax, when this support cable is impaired, the degree of injury of the support cable of this its correspondence of element numeric representation, if when this support cable is lax, the support cable of this its correspondence of element numeric representation with the current actual equivalent damage degree of its relax level mechanic equivalent,
M. identify damaged cable in the problematic support cable identified from l step, remaining is exactly slack line;
N. utilize at current cable structure steady temperature data vector T tthe cable system current actual damage vector d obtained in l step under condition athat obtain slack line with current actual equivalent damage degree that is its relax level mechanic equivalent, utilize e step obtain at current cable structure steady temperature data vector T tcurrent cable force vector F under condition and current support cable two support end points horizontal range vector, utilize c step obtain at initial Cable Structure steady temperature data vector T othe initial drift vector of the support cable under 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 trepresent support cable current steady state temperature data, utilize c step obtain at initial Cable Structure steady temperature data vector T othe support cable initial steady state temperature data represented, the temperature variant physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in c step to use, count the impact of temperature variation on support cable physics, mechanics and geometric parameter, by by slack line with damaged cable carry out mechanic equivalent calculate slack line, with the relax level of current actual equivalent damage degree equivalence, the mechanical condition of equivalence is: one, two equivalences rope without lax identical with the mechanics parameters of initial drift during not damaged, geometrical property parameter, density and material; Two, after lax or damage, the slack line of two equivalences and the Suo Li of damage rope be out of shape after overall length identical; When meeting 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 the slack line of equivalence, any change can not occur Cable Structure, and vice versa; Try to achieve according to aforementioned mechanic equivalent condition the relax level that those are judged as slack line, relax level is exactly the knots modification of support cable drift, namely determines the long adjustment amount of those ropes that need adjust the support cable of Suo Li; So just achieve lax identification and the non-destructive tests of support cable; During calculating, institute's demand power is provided by current cable force vector F corresponding element;
O. get back to e step, start the circulation next time being walked to o step by e.
CN201210173656.4A 2012-05-30 2012-05-30 Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation Expired - Fee Related CN102735465B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210173656.4A CN102735465B (en) 2012-05-30 2012-05-30 Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210173656.4A CN102735465B (en) 2012-05-30 2012-05-30 Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation

Publications (2)

Publication Number Publication Date
CN102735465A CN102735465A (en) 2012-10-17
CN102735465B true CN102735465B (en) 2015-09-30

Family

ID=46991363

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210173656.4A Expired - Fee Related CN102735465B (en) 2012-05-30 2012-05-30 Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation

Country Status (1)

Country Link
CN (1) CN102735465B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103868724A (en) * 2014-03-10 2014-06-18 东南大学 Space coordinate monitoring damaged cable load angular displacement progressive identification method
CN103852337A (en) * 2014-03-10 2014-06-11 东南大学 Strain monitoring error cable load angular displacement progressive identification method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288431A (en) * 2011-05-13 2011-12-21 东南大学 Health monitoring method of cable system based on strain monitoring in angular displacement of supporting seat
CN102297779A (en) * 2011-05-13 2011-12-28 东南大学 Method for identifying loose cables based on strain monitoring during supporting-base angular displacement

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288431A (en) * 2011-05-13 2011-12-21 东南大学 Health monitoring method of cable system based on strain monitoring in angular displacement of supporting seat
CN102297779A (en) * 2011-05-13 2011-12-28 东南大学 Method for identifying loose cables based on strain monitoring during supporting-base angular displacement

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
混凝土斜拉桥温度场的试验研究;葛耀君等;《中国公路学报》;19960630;第9卷(第2期);第1-3节 *

Also Published As

Publication number Publication date
CN102735465A (en) 2012-10-17

Similar Documents

Publication Publication Date Title
CN102749212B (en) The problem cable of temperature variation hybrid monitoring and generalized displacement of support recognition methods
CN102706594B (en) The problem cable of temperature variation space coordinate monitoring and generalized displacement of support recognition methods
CN102706600B (en) The problem cable of temperature variation angle monitor and generalized displacement of support recognition methods
CN102735465B (en) Based on the slack line recognition methods of strain monitoring during angular displacement of support temperature variation
CN102706610B (en) The problem cable of temperature variation strain monitoring and generalized displacement of support recognition methods
CN102735476B (en) The problem cable of temperature variation strain monitoring and support angular displacement identification method
CN102721551B (en) Problematic cable and support translation identification method based on space coordinate monitoring during temperature change
CN102735471B (en) Generalized displacement of support temperature variation is based on the damaged cable recognition methods of strain monitoring
CN102735469B (en) Based on the damaged cable recognition methods of strain monitoring during angular displacement of support temperature variation
CN102706662B (en) Defective cable and support angular displacement identification method based on angular monitoring of temperature change
CN102706611B (en) The problem cable of temperature variation cable force monitoring and generalized displacement of support recognition methods
CN102706577B (en) Method for identifying problem cable and support translation based on hybrid monitoring during temperature change
CN102706661B (en) Defective cable and support angular displacement identification method based on mixed monitoring of temperature change
CN102735464B (en) Based on the slack line recognition methods of angle monitor during angular displacement of support temperature variation
CN102706621B (en) Loosened cable identification method based on cable tension monitoring at moment of generalized displacement of supporting seat and temperature variation
CN102706617B (en) Loosened cable identification method based on angle monitoring at moment of generalized displacement of supporting seat and temperature variation
CN102735479B (en) Based on the slack line recognition methods of hybrid monitoring during angular displacement of support temperature variation
CN102721552B (en) Slack cable identification method based on hybrid monitoring during temperature change
CN102706587B (en) The problem cable of strain monitoring and support translation identification method during temperature variation
CN102735475B (en) The problem cable of temperature variation space coordinate monitoring and support angular displacement identification method
CN102706578B (en) Method for identifying problem cable and support translation based on angle monitoring during temperature change
CN102735472B (en) Generalized displacement of support temperature variation is based on the damaged cable recognition methods of space coordinate monitoring
CN102706584B (en) The slack line progressive-type recognition method of support settlement and temperature variation space coordinate monitoring
CN102706615B (en) The slack line progressive-type recognition method of generalized displacement of support temperature variation cable force monitoring
CN102735480B (en) Angular displacement of support temperature variation is based on the slack line recognition methods of space coordinate monitoring

Legal Events

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

Granted publication date: 20150930

Termination date: 20180530