CN104990735A - Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring - Google Patents

Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring Download PDF

Info

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

Links

Landscapes

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

Abstract

Provided is a simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring. Based on space coordinate monitoring, a mechanical calculation reference model of a cable structure is established, and calculation is performed based on the model to obtain a unit damage monitored variable value change matrix. The non-inferior solution of an evaluated object current nominal damage vector is worked out based on the approximate linear relationship of a monitored variable current value vector, a monitored variable current initial value vector, the unit damage monitored variable value change matrix and the to-be-solved evaluated object current nominal damage vector, and accordingly, the health status of a core evaluated object can be identified.

Description

Simplify cable force monitoring loading problem rope generalized displacement progressive-type recognition method
Technical field
Cable-stayed bridge, suspension bridge, the structures such as 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 the rod member only bearing tensile load are support unit, for simplicity, such structure representation is " Cable Structure " by this method, and by all ropeway carrying-ropes of Cable Structure, carrying cable, and all rod members (being also called two power rod members) only bearing axial tension or axial compression load, be collectively referred to as simplicity " cable system ", ropeway carrying-rope is censured with " support cable " this noun in this method, carrying cable and only bear the rod member of axial tension or axial compression load, sometimes referred to as " rope ", so when using " rope " this word below, two power rod members are just referred to truss-frame structure reality.Impaired and the lax pair Cable Structure of support cable is safely a significant threat, and damaged cable and slack line are referred to as the support cable of unsoundness problem, referred to as problem cable by this method.In structure military service process, the correct identification of the health status of support cable or cable system is related to the safety of whole Cable Structure.When there is health problem in support cable, in Cable Structure military service process, may generalized displacement be there is in Cable Structure bearing, the load that Cable Structure is born also may change simultaneously, even if in fact the health status of support cable does not change, the load that Cable Structure is born also may change separately, at this complex condition, this method identifies generalized displacement of support and problem cable based on cable force monitoring (monitored Suo Li is called " monitored amount " by this method), belongs to engineering structure health monitoring field.
Background technology
Reject load change, Cable Structure generalized displacement of support on the impact of Cable Structure health status recognition result, thus identifies the change of the health status of structure exactly, is problem in the urgent need to address at present; Rejecting load change, the impact of Cable Structure health status change on Cable Structure generalized displacement of support recognition result, thus identify Cable Structure generalized displacement of support exactly, is also problem in the urgent need to address at present; This method discloses a kind of effective, cheap method solving this two problems.
Summary of the invention
Technical matters: this method discloses a kind of method, under the condition that cost is lower, achieve two kinds of functions, be respectively, one, rejecting generalized displacement of support and load change are on the impact of Cable Structure health status recognition result, thus identify the health status of support cable exactly; Two, this method can also reject load change and the impact of Cable Structure health status change on Cable Structure generalized displacement of support recognition result, thus identifies Cable Structure generalized displacement of support exactly.
Suo Changdu under support cable free state (now rope tensility also claims Suo Li to be 0) (is called drift, this method specially refers to that support cable two supports the drift of that section of rope between end points) can change, one of object of this method will identify the support cable that drift there occurs change exactly, and identify the knots modification of their drift, this knots modification is that the cable force adjustment of this rope provides direct basis, conveniently, the support cable that drift changes by this method is referred to as slack line.
Technical scheme: in the method, the coordinate of bearing about the X, Y, Z axis of Descartes's rectangular coordinate system is censured with " bearing volume coordinate ", alternatively becoming is the volume coordinate of bearing about X, Y, Z axis, bearing is called the volume coordinate component of bearing about this axle about the concrete numerical value of the volume coordinate of some axles, and a volume coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the volume coordinate of some axles; The angular coordinate of bearing about X, Y, Z axis is censured with " bearing angular coordinate ", bearing is called the angular coordinate component of bearing about this axle about the concrete numerical value of the angular coordinate of some axles, and an angular coordinate component also with bearing in this method expresses the concrete numerical value of bearing about the angular coordinate of some axles; Censure bearing angular coordinate and bearing volume coordinate entirety with " bearing generalized coordinate ", a generalized coordinate component also with bearing in this method expresses the volume coordinate of bearing about an axle or the concrete numerical value of angular coordinate; Bearing is called support wire displacement about the change of the coordinate of X, Y, Z axis, and alternatively the change of bearing volume coordinate is called support wire displacement, and a translational component also with bearing in this method expresses the concrete numerical value of bearing about the displacement of the lines of some axles; Bearing is called angular displacement of support about the change of the angular coordinate of X, Y, Z axis, and an angular displacement component also with bearing in this method expresses the concrete numerical value of bearing about the angular displacement of some axles; Generalized displacement of support censures support wire displacement and angular displacement of support is all, and a generalized displacement component also with bearing in this method expresses the displacement of the lines of bearing about some axles or the concrete numerical value of angular displacement; Support wire displacement also can be described as translational displacement, and support settlement is support wire displacement or the translational displacement component at gravity direction.
The external force that object, structure are born can be described as load, and load comprises face load and volume load.Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds.Volume load be continuous distribution in the load of interior of articles each point, as deadweight and the inertial force of object.
Centre-point load is divided into concentrated force and concentrated couple two kinds, in a coordinate system, such as in Descartes's rectangular coordinate system, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component are called a load, the now change of load is embodied as the change of a concentrated force component or a concentrated couple component.
Distributed load is divided into line distributed load and EDS maps load, the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, the size distribution intensity of distributed load is expressed, distribution intensity distribution characteristics is (such as uniform, sine function equal distribution feature) and amplitude is expressed, and (such as two distributed loads are all uniform, but its amplitude is different, can well-distributed pressure be example so that the concept of amplitude to be described: same structure bears two different well-distributed pressures, two distributed loads are all uniformly distributed loads, but the amplitude of a distributed load is 10MPa, the amplitude of another distributed load is 50MPa).If load is actually distributed load, when this method talks about the change of load, in fact refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of distributed load and distribution intensity is constant.In a coordinate system, a distributed load can resolve into several components, if the amplitude of the respective distribution intensity of several components of this distributed load changes, and the ratio of change is all not identical, so in the method the component of these several distributed loads is regarded as the independently distributed load of same quantity, now load just represents the component of a distributed load, also component identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized a distributed load or be called a load.
Volume load is that continuous distribution is in the load of interior of articles each point, as deadweight and the inertial force of object, the description of volume load at least comprises the zone of action of volume load and the size of volume load, the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics is (such as uniform, linear function equal distribution feature) and amplitude is expressed, and (such as two individual stow lotuses are all uniform, but its amplitude is different, can conduct oneself with dignity for example is to illustrate the concept of amplitude: the material of two parts of same structure is different, therefore density is different, so although this volume load suffered by two parts is all uniform, but the amplitude of the volume load suffered by a part may be 10kN/m 3, the amplitude of the volume load suffered by another part is 50kN/m 3).If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes.In a coordinate system, one individual stow lotus can resolve into several components (such as in Descartes's rectangular coordinate system, volume load can resolve into the component of three axles about coordinate system, that is, in Descartes's rectangular coordinate system, volume load can resolve into three components), if the amplitude of the respective distribution intensity of several components of this volume load changes, and the ratio of change is all not identical, so in the method the component of this several body stow lotus is regarded as the independently load of same quantity, also the volume sharing part of the load identical for the amplitude changing ratio of the intensity that wherein distributes can be synthesized an individual stow lotus or be called a load.
When load is embodied as centre-point load, in the method, " load unit change " in fact refers to " unit change of centre-point load ", similar, " load change " specifically refers to " change of the size of centre-point load ", " load change amount " specifically refers to " variable quantity of the size of centre-point load ", " load change degree " specifically refers to " intensity of variation of the size of centre-point load ", " the actual change amount of load " refers to " the actual change amount of the size of centre-point load ", " load changed " refers to " centre-point load that size changes ", briefly, now " so-and-so load so-and-so change " refers to " size of so-and-so centre-point load so-and-so change ".
When load is embodied as distributed load, in the method, " load unit change " in fact refers to " unit change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, similar, " load change " specifically refers to " change of the amplitude of the distribution intensity of distributed load ", and the distribution characteristics of distributed load is constant, " load change amount " specifically refers to " variable quantity of the amplitude of the distribution intensity of distributed load ", " load change degree " specifically refers to " intensity of variation of the amplitude of the distribution intensity of distributed load ", " the actual change amount of load " specifically refers to " the actual change amount of the amplitude of the distribution intensity of distributed load ", " load changed " refers to " distributed load that changes of amplitude of distribution intensity ", briefly, now " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so distributed load so-and-so change ", and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant.
When load is embodied as volume load, in the method, " load unit change " in fact refers to " unit change of the amplitude of the distribution intensity of volume load ", similar, " load change " refers to " change of the amplitude of the distribution intensity of volume load ", " load change amount " refers to " variable quantity of the amplitude of the distribution intensity of volume load ", " load change degree " refers to " intensity of variation of the amplitude of the distribution intensity of volume load ", " the actual change amount of load " refers to " the actual change amount of the amplitude of the distribution intensity of volume load ", " load changed " refers to " the volume load that changes of amplitude of distribution intensity ", briefly, " so-and-so load so-and-so change " refers to " amplitude of the distribution intensity of so-and-so volume load so-and-so change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.
This method specifically comprises:
Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, namely the quantity of evaluation object is N; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; This method title " core evaluation object " specially refers to evaluated support cable in " evaluation object " and generalized displacement of support component, if the quantity sum of evaluated support cable and generalized displacement of support component is P, namely the quantity of core evaluation object is P, and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; If total M in cable system 1root support cable, Cable Structure rope force data comprises this M 1the Suo Li of root support cable, this method is at the whole M of monitoring 1on the basis of root supporting cable force, Cable Structure artificially increases M 2root rope, is called sensing rope, by this M newly increased of monitoring in cable structure health monitoring process 2the Suo Li of root sensing rope; Comprehensive above-mentioned monitored amount, M the Suo Li that whole Cable Structure has M root rope is monitored, and namely have M monitored amount, wherein M is M 1with M 2sum; M should be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; The M newly increased 2the rigidity of root sensing rope is compared with the rigidity of any support cable of Cable Structure, should be much smaller; The M newly increased 2the Suo Li of each sensing rope of root sensing rope should be more much smaller than the Suo Li of any support cable of Cable Structure, even if can ensure this M newly increased like this 2root sensing rope has occurred damaging or lax, very little on the impact of the stress of other components of Cable Structure, strain, distortion; The M newly increased 2on the xsect of root sensing rope, normal stress should be less than its fatigue limit, and these requirements can ensure the M newly increased 2can not fatigue damage be there is in root sensing rope; The M newly increased 2the two ends of root sensing rope should fully anchoring, ensures to there will not be lax; The M newly increased 2root sensing rope should obtain sufficient anti-corrosion protection, ensures the M newly increased 2damage can not be there is and relax in root sensing rope; For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount "; To M monitored amount serial number, this method represents this numbering, j=1,2,3 with variable j ..., M, this numbering will be used for generating vector sum matrix in subsequent step; The M newly increased in the method 2root sensing rope is as a part for Cable Structure, and when hereinafter mentioning Cable Structure again, Cable Structure comprises increases M 2cable Structure before root sensing rope and the M newly increased 2root sensing rope, refers to comprise the M newly increased when that is hereinafter mentioning Cable Structure 2the Cable Structure of root sensing rope; Measure the M obtaining newly increasing 2the method of the Suo Li of root sensing rope is same as the M of Cable Structure 1the measuring method of the Suo Li of root support cable, hands over later no longer one by one; When carrying out any measurement to the support cable of Cable Structure, simultaneously to the M newly increased 2root sensing rope carries out same measurement, hands over no longer one by one later; The M newly increased 2root sensing rope except there is not damage and be lax, to the M newly increased 2the requirement of the quantity of information of root sensing rope is identical with preparation method with the requirement of the quantity of information of the support cable of Cable Structure with preparation method, hands over no longer one by one later; When setting up the various mechanical model of Cable Structure later, by the M newly increased 2the support cable that Cable Structure treated as by root sensing rope is treated; Below, except the damage of mentioning support cable and lax occasion, the support cable that said support cable comprises Cable Structure when mentioning support cable and the M newly increased 2root sensing rope; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centre-point load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated force component or a concentrated couple component; Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the one-component of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;
B. survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses; Actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; According to comprising Cable Structure design data, the data of completion data obtain the length of all support cables when free state and Suo Li are 0, in free state time cross-sectional area and in free state time the weight of unit length, 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;
C. actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Non-destructive Testing Data comprising support cable, and the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C o, monitored amount initial value vector C othe coding rule of coding rule and M monitored amount identical, the initial generalized displacement measurement data of support cable initial health data, Cable Structure bearing and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d o, vectorial d orepresent with initial mechanical Calculation Basis model A othe initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d oelement number equal N, d oelement and evaluation object be one-to-one relationship, vectorial d othe coding rule of element identical with the coding rule of evaluation object, if d oevaluation object corresponding to some elements be support cable, so a d in cable system othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses 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 d oevaluation object corresponding to some elements be some generalized displacement components of some bearings, so d othe numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d oevaluation object corresponding to some elements be some load, get d in this method othis element numerical value be 0, the initial value representing the change of this load is 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d oin each element numerical value relevant to Cable Structure generalized displacement of support get 0, 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 can think structure original state be not damaged without relaxed state time, vectorial d oin each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data refer to the bearing generalized coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A otime, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point,
All Cable Structure data that the physical and mechanical properties parameter 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, support cable initial health data, Cable Structure bearing initial generalized displacement measurement data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure and preceding step obtain, set up the initial mechanical Calculation Basis model A of Cable Structure 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 oevaluation object health status with evaluation object initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent; d oa oparameter, 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, A in the method o, C oand d oconstant;
E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting i o; The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d i o, d i ocable Structure A when representing that this circulation starts i othe health status of evaluation object, d i odefinition mode and d odefinition mode identical, d i oelement and d oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C i orepresent, vectorial C i odefinition mode and vectorial C odefinition mode identical, C i oelement and C oelement one_to_one corresponding, monitored amount current initial value vector C i orepresent and correspond to A i othe concrete numerical value of all monitored amount; d i oa i ocharacterisitic parameter, C i oby A i omechanics Calculation result composition; When first time, circulation started, A i obe designated as A 1 o, set up A 1 omethod for making A 1 oequal A o; When first time, circulation started, d i obe designated as d 1 o, set up d 1 omethod for making d 1 oequal d o; When first time, circulation started, C i obe designated as C 1 o, set up C 1 omethod for making C 1 oequal C o;
F. from entering the circulation being walked to p step by f here;
G. in structure military service process, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C i, vectorial C idefinition mode and vectorial C odefinition mode identical, C ielement and C oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time, monitored amount current value vector C is obtained in actual measurement iwhile, to the M newly increased 2root sensing rope carries out Non-Destructive Testing, therefrom identify and occur damage or lax sensing rope, according to monitored amount coding rule, the element that the sensing rope that to damage with the appearance identified or relax is corresponding is removed from each vector according to monitored amount coding rule numbering occurred before this method, also the element damaged with the appearance identified or lax sensing rope is corresponding is no longer there is in each vector sum matrix occurred after this method, no longer comprise when mentioning sensing rope after this method being identified here and occur damage or lax sensing rope, the Suo Li being identified and occurring damage or lax sensing rope is no longer comprised here when mentioning monitored amount after this method, identify several from Cable Structure and occur damage or lax sensing rope, just by M 2same quantity is reduced with M, monitored amount current value vector C is obtained in actual measurement isynchronization, actual measurement obtain all M in Cable Structure 1the rope force data of root support cable, all these rope force datas composition current cable force vector F i, vectorial F ielement and vectorial F othe coding rule of element identical, monitored amount current value vector C is obtained in actual measurement isynchronization, Actual measurement obtains all M 1the volume coordinate of two supporting end points of root support cable, 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,
H. at current initial mechanical Calculation Basis model A i obasis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C iwith evaluation object unit change vector D i u;
H1., when i-th circulation starts, directly Δ C is obtained by method listed by step h2 to step h4 iand D i u; In other moment, when in step g to A i oafter upgrading, Δ C must be regained by method listed by step h2 to step h4 iand D i uif, not to A in step g i oupgrade, then directly proceed to step I herein and carry out follow-up work;
H2. at current initial mechanical Calculation Basis model A i obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating, according to the coding rule of evaluation object, calculate successively, calculating hypothesis each time only has an evaluation object on the basis of original damage or generalized displacement or load, increase unit damage or unit generalized displacement or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is the generalized displacement component in a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, if this evaluation object is a load, just suppose that this load increases load unit change again, use D i ukrecord the unit damage of this increase or unit generalized displacement or load unit change, wherein k represents the numbering of the evaluation object increasing unit damage or unit generalized displacement or load unit change, D i ukevaluation object unit change vector D i uan element, evaluation object unit change vector D i uthe coding rule of element and vectorial d othe coding rule of element identical, the evaluation object increasing unit damage or unit generalized displacement or load unit change in calculating each time is again different from during other time calculates the evaluation object increasing unit damage or unit generalized displacement or load unit change again, calculate the current calculated value all utilizing mechanics 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, when suppose a kth evaluation object increase again unit damage or unit generalized displacement or load unit change time, use C i tkrepresent corresponding " monitored amount calculation current vector ", when giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target, C i tkdefinition mode and vectorial C odefinition mode identical, C i tkelement and C oelement one_to_one corresponding,
H3. the vectorial C calculated each time i tkdeduct vectorial C i oobtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated divided by this unit damage or unit generalized displacement or load unit change numerical value supposed i k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ";
H4. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively i"; Unit damage monitored numerical quantity transformation matrices Δ C ieach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ievery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of row and vectorial d othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of coding rule and M monitored amount of row identical;
I. current nominal fatigue vector d is defined i cwith current actual damage vector d i, d i cand d ielement number equal the quantity of evaluation object, d i cand d ielement and evaluation object between be one-to-one relationship, d i celement numerical value represent the nominal fatigue degree of corresponding evaluation object or nominal generalized displacement or nominal load variable quantity, d i cand d iwith evaluation object initial damage vector d oelement number rule identical, d i celement, d ielement and d oelement be one-to-one relationship;
J. according to monitored amount current value vector C iwith " monitored amount current initial value vector C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " current nominal fatigue vector d i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1 i cother outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d i c;
formula 1
K. the current actual damage vector d utilizing formula 2 to express ia kth element d i kwith evaluation object current initial damage vector d i oa kth element d i okwith current nominal fatigue vector d i ca kth element d i ckbetween relation, calculate current actual damage vector d iall elements;
K=1 in formula 2,2,3 ..., N; d i krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system i krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d i kthe degree of the lax or damage of this support cable of numerical response; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, evaluation object current actual damage vector d iin correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent; If this evaluation object is load, so a d i krepresent the actual change amount of this load; If this evaluation object is generalized displacement component, so a d of a bearing i krepresent its current actual generalized displacement numerical value; By current for evaluation object actual damage vector d iin with M 1the M that root support cable is relevant 1individual element takes out, composition support cable current actual damage vector d ci, support cable current actual damage vector d cithe coding rule of element and Initial cable force vector F othe coding rule of element identical; Support cable current actual damage vector d cih element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M 1; Support cable current actual damage vector d cimiddle numerical value be not 0 element correspond to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line; The support cable corresponding with damaged cable current actual damage vector d ciin element numerical expression be the current actual damage of this damaged cable, represent when element numerical value is 100% that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; After identifying slack line, utilize support cable current actual damage vector d cithese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize current cable force vector F isupport end points horizontal range vector with current support cable two, utilize initial drift vector, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the vectorial F of Initial cable force of the support cable obtained in b step othe physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in b step to use, 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, mechanic equivalent condition 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 mechanic 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 by current cable force vector F icorresponding element provides; Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far this method achieve reject generalized displacement of support, load change impact, the problem cable identification of Cable Structure, achieve simultaneously reject generalized displacement of support and support cable health status variable effect, the identification of load change amount, also achieve reject load change and support cable health status variable effect, the identification of generalized displacement of support; So far this method achieves the accurate identification of the health status of core evaluation object with a kind of effective, cheap method; Exact value is departed to the recognition result of the health status of secondary evaluation object more, therefore will not accept and believe, only require the correct health status identifying core evaluation object in the method;
L. current nominal fatigue vector d is tried to achieve i cafter, set up mark vector B according to formula 3 i, formula 4 gives mark vector B ithe definition of a kth element;
formula 3
Element B in formula 4 i kmark vector B ia kth element, D i ukevaluation object unit change vector D i ua kth element, d i ckevaluation object current nominal fatigue vector d i ca kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;
If m. mark vector B ielement be 0 entirely, then get back to step f continue this circulation; If mark vector B ielement be not 0 entirely, then enter next step, i.e. step n;
N. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the i-th+1 time circulation i+1 oeach element;
D in formula 5 i+1 okthe evaluation object current initial damage vector d next time, namely needed for the i-th+1 time circulation i+1 oa kth element, d i okthis, i.e. the evaluation object of i-th circulation current initial damage vector d i oa kth element, D i ukthe evaluation object unit change vector D of i-th circulation i ua kth element, B i kthe mark vector B of i-th circulation ia kth element, k=1 in formula 5,2,3 ..., N;
O. at initial mechanical Calculation Basis model A obasis on, make the health status of rope be d i+1 oafter obtain be exactly next time, namely the i-th+1 time circulation needed for Mechanics Calculation benchmark model A i+1; Obtain A i+1after, obtain A by Mechanics Calculation i+1in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C i+1 o;
P. get back to step f, start to circulate next time.
Beneficial effect: structural healthy monitoring system is first by using sensor to carry out long-term on-line monitoring to structural response, after obtaining Monitoring Data, (or off-line) analysis is online carried out to it and obtain structural health conditions data, due to the complicacy of structure, structural healthy monitoring system needs to use a large amount of sensor equipment to carry out monitoring structural health conditions, therefore its cost is usually quite high, and therefore cost problem is a subject matter of limit structural health monitoring technique application.On the other hand, the correct identification of the health status of core evaluation object (such as suspension cable) is the indispensable ingredient of the correct identification of structural health conditions, or even they are whole, and the impact of correct identification on the correct identification of the health status of Cable Structure of the change (such as by the change of the quality and quantity of the automobile of cable-stayed bridge) of secondary evaluation object (load that such as structure is born) is very little, or even unwanted.But the quantity of the quantity of secondary evaluation object and core evaluation object is normally suitable, the quantity of secondary evaluation object is also usually greater than the quantity of core evaluation object, and the quantity of such evaluation object is usually many times of the quantity of core evaluation object.When secondary evaluation object (load) changes, in order to accurately identify core evaluation object, conventional method requires that the quantity of monitored amount (using sensor device measuring to obtain) must be more than or equal to the quantity of evaluation object, when the number ratio of the secondary evaluation object changed is larger (in fact often so), the quantity of the sensor equipment required for structural healthy monitoring system is very huge, therefore the cost of structural healthy monitoring system will become very high, unacceptablely even high.Inventor studies discovery, in the secondary evaluation object (normal load that such as structure is born, the normal load of structure refers to that the load that structure is being born is no more than the structure allowable load limited according to structural design book or structure completion book) change less time (be exactly that structure only bears normal load for load, whether the load that structure is born is normal load, can be observed by methods such as naked eyes and determine, if find that the load that structure is born is not normal load, so artificially remove, after removing improper load, structure just only bears normal load), the amplitude of variation (this instructions is called " secondary response ") of the structural response caused by them much smaller than core evaluation object change (such as support cable is impaired) caused by the amplitude of variation (this instructions is called " core response ") of structural response, secondary response and core respond total change (this instructions is called " global response ") that sum is structural response, obvious core response dominate in global response, based on this, find to choose when determining monitored amount quantity to be a bit larger tham core evaluation object quantity even if inventor studies, but much smaller than the numerical value (this method is exactly do like this) of evaluation object quantity, even if that is adopt the relatively few a lot of sensor equipment of quantity, still the state of health data of core evaluation object can accurately be obtained, meet the core demand of structural health conditions monitoring, therefore this method cost of structural healthy monitoring system of advising is more much lower than the cost of the structural healthy monitoring system required by conventional method apparently, that is this method can realize to the health status of the core evaluation object of Cable Structure with the much lower condition of cost assessment, can this benefit be used structural health monitoring technology is very important.
Embodiment
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.
The first step: the quantity first confirming the load that may change that Cable Structure is born.According to the feature of the load that Cable Structure is born, confirm wherein " load likely changed ", or all load is considered as " load likely changed ", if total JZW the load that may change, i.e. total JZW secondary evaluation object.
If the quantity sum of the quantity of the quantity of generalized displacement of support component of Cable Structure, the support cable of Cable Structure and JZW " load likely changed " is N, i.e. total N number of evaluation object.To evaluation object serial number, this numbering will be used for generating vector sum matrix in subsequent step.
If the quantity sum of evaluated support cable and generalized displacement of support component is P, namely the quantity of core evaluation object is P, if the quantity of evaluated generalized displacement of support component is Z, if the quantity of evaluated support cable is M 1, structure rope force data comprises this M 1the Suo Li of root support cable, this method is at the whole M of monitoring 1on the basis of root supporting cable force, increase M 2other monitored amounts individual.
The M increased 2other monitored amounts individual remain Suo Li, are described below:
Before structural healthy monitoring system is started working, first in Cable Structure, artificially increase M 2(M 2be not less than 4) root rope, be called sensing rope, the M newly increased 2the rigidity of root sensing rope is compared with the rigidity of any support cable of Cable Structure, should be much little, such as little 20 times, the M newly increased 2the Suo Li of root sensing rope should be less, and such as its xsect normal stress should be less than its fatigue limit, and these requirements can ensure the M newly increased 2can not there is fatigue damage in root sensing rope, the M newly increased 2the two ends of root sensing rope should fully anchoring, ensures to there will not be lax, the M newly increased 2root sensing rope should obtain sufficient anti-corrosion protection, ensures the M newly increased 2can not be there is damage and relax in root sensing rope, will this M newly increased of monitoring in monitoring structural health conditions process 2the Suo Li of root sensing rope.
The mode of many increase sensing ropes can also be adopted to ensure the reliability of health monitoring, such as, make M 2be not less than 8, in monitoring structural health conditions process, only selecting the rope force data of intact sensing rope wherein, (be called actual operable monitored amount together with the rope force data of support cable, recording its quantity is that K, K must not be less than 4+M 1) and corresponding Cable Structure monitored amount unit change matrix Δ C carry out health state evaluation.By this M newly increased of monitoring in monitoring structural health conditions process 2the Suo Li of root sensing rope.The M newly increased 2the position that root sensing rope should be installed structurally, personnel are easy to arrival, is convenient to personnel and carries out Non-Destructive Testing to it.
The M newly increased in the method 2root sensing rope is as a part for Cable Structure, and when hereinafter mentioning Cable Structure again, Cable Structure comprises increases M 2cable Structure before root sensing rope and the M newly increased 2root sensing rope, refers to comprise the M newly increased when that is hereinafter mentioning Cable Structure 2the Cable Structure of root sensing rope.Therefore hereinafter mention when obtaining " Cable Structure steady temperature data " according to " the temperature survey calculating method of the Cable Structure of this method " survey calculation, Cable Structure wherein comprises the M newly increased 2root sensing rope, " the Cable Structure steady temperature data " that obtain comprise the M newly increased 2the steady temperature data of root sensing rope, obtain the M newly increased 2the method of the steady temperature data of root sensing rope is same as the M of Cable Structure 1the preparation method of the steady temperature data of root support cable, hands over later no longer one by one; Measure the M obtaining newly increasing 2the method of the Suo Li of root sensing rope is same as the M of Cable Structure 1the measuring method of the Suo Li of root support cable, hands over later no longer one by one; When carrying out any measurement to the support cable of Cable Structure, simultaneously to the M newly increased 2root sensing rope carries out same measurement, hands over no longer one by one later; The M newly increased 2root sensing rope except there is not damage and be lax, the M newly increased 2the quantity of information of root rope is identical with the quantity of information of the support cable of Cable Structure, hands over no longer one by one later.When setting up the various mechanical model of Cable Structure later, by the M newly increased 2the M of Cable Structure treated as by root sensing rope 1root support cable is treated, and except mentioning the damage of support cable and lax occasion, comprises the M newly increased when other occasions mention support cable 2root rope.
Comprehensive above-mentioned monitored amount, whole Cable Structure has M (M=M 1+ M 2) M the monitored amount of root rope, the quantity that M must not be less than core evaluation object adds the quantity that 4, M is less than evaluation object.
For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount ".To M monitored amount serial number, this numbering will be used for generating vector sum matrix in subsequent step.This method represents this numbering, j=1,2,3 with variable j ..., M.
Second step: set up initial mechanical Calculation Basis model A o.
When Cable Structure is completed, or before setting up health monitoring systems, 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.
Obtaining monitored amount initial value vector C owhile, use conventional method (consult reference materials or survey) to obtain physical parameter and the mechanical property parameters (such as elastic modulus, Poisson ratio) of the various materials that Cable Structure uses.
Obtaining monitored amount initial value vector C owhile, utilize the elastic modulus data composition support cable initial elastic modulus vector E of support cable wherein o; Actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; According to comprising Cable Structure design data, the data of completion data obtain the length of all support cables when free state and Suo Li are 0, in free state time cross-sectional area and in free state time the weight of unit length, form the initial drift vector l of support cable successively o, initial free cross-sectional area vector A owith the weight vector ω of initial free unit length o, support cable initial elastic modulus vector E o, support cable initial drift vector l o, initial free cross-sectional area vector A owith the weight vector ω of initial free unit length othe coding rule of element and Initial cable force vector F othe coding rule of element identical.
Obtaining monitored amount initial value vector C owhile, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Non-destructive Testing Data etc. that the Actual measurement data of Cable Structure comprise support cable can express the data of the health status of rope, the initial geometric data of Cable Structure, rope force data, draw-bar pull data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure bearing generalized coordinate data, the initial generalized displacement measurement data of Cable Structure bearing, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure modal data, structural strain data, structural point measurement data, the measured datas such as structure space measurement of coordinates data.Initial Cable Structure bearing generalized coordinate data refer to the bearing generalized coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A otime, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point 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.Utilize the Non-destructive Testing Data etc. of support cable can express the data of the health status of support cable, the initial generalized displacement measurement data of Cable Structure bearing and Cable Structure load measurement data and set up evaluation object initial damage vector d o, use d orepresent that Cable Structure is (with initial mechanical Calculation Basis model A orepresent) the initial health of evaluation object.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 can think structure original state be not damaged without relaxed state time, vectorial d oin each element numerical value relevant to support cable get 0; If when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d oin each element numerical value relevant to Cable Structure generalized displacement of support get 0; If d oevaluation object corresponding to some elements be some load, get d in this method othis element numerical value be 0, the initial value representing the change of this load is 0.The 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, utilizes mechanics method (such as finite element method) to 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, 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 evaluation object initial damage vector d of middle support cable 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 oevaluation object health status with evaluation object initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent.D oa oparameter, C oby A omechanics Calculation result composition.
3rd step: in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting i o; The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d i o, d i ocable Structure A when representing that this circulation starts i othe health status of evaluation object, d i odefinition mode and d odefinition mode identical, d i oelement and d oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C i orepresent, vectorial C i odefinition mode and vectorial C odefinition mode identical, C i oelement and C oelement one_to_one corresponding, monitored amount current initial value vector C i orepresent and correspond to A i othe concrete numerical value of all monitored amount; d i oa i ocharacterisitic parameter; C i oby A i omechanics Calculation result composition; When first time, circulation started, A i obe designated as A 1 o, set up A 1 omethod for making A 1 oequal A o; When first time, circulation started, d i obe designated as d 1 o, set up d 1 omethod for making d 1 oequal d o; When first time, circulation started, C i obe designated as C 1 o, set up C 1 omethod for making C 1 oequal C o.
4th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (such as containing cable force measurement system, signal conditioner etc.), signal (data) collector, computing machine and communication alert equipment.Each monitored amount 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 evaluation object running Cable Structure is then responsible for by computing machine, comprises the signal that the transmission of tracer signal collector comes; When monitoring evaluation object health status and changing, computer control communication panalarm is reported to the police to monitor staff, owner and (or) the personnel that specify.
5th step: establishment the system software of installation and operation this method on computers, the function (i.e. all work that can complete with computing machine in this specific implementation method) such as monitoring, record, control, storage, calculating, notice, warning that this software will complete this method required by task and wants.
6th step: step starts circulation running thus, in structure military service process, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C i, vectorial C idefinition mode and vectorial C odefinition mode identical, C ielement and C oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time.
Monitored amount current value vector C is obtained in actual measurement isynchronization, actual measurement obtain all M in Cable Structure 1the rope force data of root support cable, all these rope force datas composition current cable force vector F i, vectorial F ielement and vectorial F othe coding rule of element identical; Monitored amount current value vector C is obtained in actual measurement isynchronization, Actual measurement obtains all M 1the volume coordinate of two of root support cable supporting end points, the difference of the volume coordinate of two supporting end points component is in the horizontal direction exactly two supporting end points horizontal ranges, all M 1two supporting end points horizontal range data of root support cable form current support cable two and support end points horizontal range vector l i x, current support cable two supports end points horizontal range vector l i xthe coding rule of element and Initial cable force vector F othe coding rule of element identical.
Monitored amount current value vector C is obtained in actual measurement iwhile, to the M newly increased 2root sensing rope carries out Non-Destructive Testing, such as UT (Ultrasonic Testing), visual examination, infrared imaging checks, therefrom identify and occur damage or lax sensing rope, according to monitored amount coding rule, the element that the sensing rope that to damage with the appearance identified or relax is corresponding is removed from each vector according to monitored amount coding rule numbering occurred before this method, also the element damaged with the appearance identified or lax sensing rope is corresponding is no longer there is in each vector sum matrix occurred after this method, no longer comprise when mentioning sensing rope after this method being identified here and occur damage or lax sensing rope, the Suo Li being identified and occurring damage or lax sensing rope is no longer comprised here when mentioning monitored amount after this method, identify several from Cable Structure and occur damage or lax sensing rope, just by M 2same quantity is reduced with M.
7th step: at current initial mechanical Calculation Basis model A i obasis on, carry out several times Mechanics Calculation according to step a to steps d, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C iwith evaluation object unit change vector D i u.
A., when i-th circulation starts, Δ C is obtained by method listed by step b to steps d iand D i u; In other moment, directly proceed to the 8th step and carry out follow-up work.
B. at current initial mechanical Calculation Basis model A i obasis on carry out several times Mechanics Calculation, vectorial d i orepresent A i othe health status of evaluation object, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating; Calculate hypothesis each time and only have an evaluation object at vectorial d i othe basis of the health status of the evaluation object represented there is unit damage or unit generalized displacement or load unit to change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable is at vectorial d i othe basis that this support cable represented has a damage there is again unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is the generalized displacement component in a direction of a bearing, just suppose this bearing at this sense of displacement at vectorial d i othere is unit generalized displacement again (if such as this evaluation object is the translational component in the x direction of a bearing in the basis that this bearing represented has a generalized displacement, just suppose that this bearing has unit line displacement in x direction, such as get 1mm, if this evaluation object is the angular displacement component around x-axis of a bearing, just suppose that this bearing is around the angular displacement of xZhou You unit, such as get 100,000/radian), if this evaluation object is a load, just suppose that this load is at vectorial d i othe basis that this load represented has a variable quantity increases again load unit change (if this load is distributed load, and this distributed load is line distributed load, load unit change can get 1kN/m, 2kN/m, 3kN/m or 1kNm/m, 2kNm/m, 3kNm/m etc. for unit change; If this load is distributed load, and this distributed load is EDS maps load, and load unit change can get 1MPa, 2MPa, 3MPa or 1kNm/m 2, 2kNm/m 2, 3kNm/m 2deng be unit change; If this load is centre-point load, and this centre-point load is couple, and load unit change can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centre-point load, and this centre-point load is concentrated force, and load unit change can get 1kN, 2kN, 3kN etc. for unit change; If this load is volume load, load unit change can get 1kN/m 3, 2kN/m 3, 3kN/m 3deng be unit change), use D i ukrecord this unit damage or unit generalized displacement or load unit to change, wherein k represents the numbering of the evaluation object that unit damage or unit generalized displacement or load unit change occur, D i ukevaluation object unit change vector D i uan element, evaluation object unit change vector D i uthe coding rule of element and vectorial d othe coding rule of element identical; Occur in calculating each time that the evaluation object of unit damage or unit generalized displacement or load unit change is different from during other time calculates the evaluation object occurring unit damage or unit generalized displacement or load unit change, calculate the current calculated value all utilizing mechanics 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; When supposing that a kth evaluation object has unit damage or unit generalized displacement or load unit to change, available C i tkrepresent corresponding " monitored amount calculation current vector "; When giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target; C i tkdefinition mode and vectorial C odefinition mode identical, C i tkelement and C oelement one_to_one corresponding.
C. the vectorial C calculated each time i tkdeduct vectorial C i oobtain a vector, then the unit damage of supposition during each element of this vector is calculated divided by this or unit generalized displacement or load unit change numerical value D i ukafter obtain " the numerical value change of a monitored amount vector δ C i k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ".
D. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively i"; Unit damage monitored numerical quantity transformation matrices Δ C ieach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ievery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of row and vectorial d othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of coding rule and M monitored amount of row identical.
8th step: set up linear relationship error vector e iwith vectorial g i.Utilize (" the monitored amount current initial value vector C of data above i o", " unit damage monitored numerical quantity transformation matrices Δ C i"); while the 8th step calculates each time; namely while calculating each time in hypothesis evaluation object and only having the increase unit damage of an evaluation object or unit generalized displacement or load unit change; when hypothesis kth (k=1,2,3; ..., N), when individual evaluation object increases unit damage or unit generalized displacement or load unit change, calculate composition injury vector each time, use d i tkrepresent this injury vector, corresponding monitored amount calculation current vector is C i tk(see the 8th step), injury vector d i tkelement number equal the quantity of evaluation object, vectorial d i tkall elements in only have the numerical value of an element to get to calculate each time in hypothesis increase the unit damage of the evaluation object of unit damage or unit generalized displacement or load unit change or unit generalized displacement or load unit changing value, d i tkthe numerical value of other element get 0, that be not numbering and the supposition of the element of 0 increase the evaluation object that unit damage or unit generalized displacement or load unit change corresponding relation, be identical with the element of the same numbering of other vectors with the corresponding relation of this evaluation object; d i tkwith evaluation object initial damage vector d oelement number rule identical, d i tkelement and d oelement be one-to-one relationship.By C i tk, C i o, Δ C i, d i tkbring formula (1) into, obtain a linear relationship error vector e i k, calculate a linear relationship error vector e each time i k; e i ksubscript k represent kth (k=1,2,3 ..., N) and individual evaluation object increases unit damage or unit generalized displacement or load unit change.There is N number of evaluation object just to have N calculating, just have N number of linear relationship error vector e i k, by this N number of linear relationship error vector e i kobtaining a vector after addition, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N i.Vector g iequal final error vector e i.By vectorial g ibe kept on the hard disc of computer of operation health monitoring systems software, for health monitoring systems software application.
9th step: define current nominal fatigue vector d i cwith current actual damage vector d i, d i cand d ielement number equal the quantity of evaluation object, d i cand d ielement and evaluation object between be one-to-one relationship, d i cand d ielement numerical value represent the degree of injury of corresponding evaluation object or unit generalized displacement or load change degree, d i cand d iwith evaluation object initial damage vector d oelement number rule identical, d i celement, d ielement and d oelement be one-to-one relationship.
Tenth step: according to monitored amount current value vector C iwith " monitored amount current initial value vector C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " current nominal fatigue vector d i c" between exist linear approximate relationship, this linear approximate relationship can be expressed as formula (2), according to multi-objective optimization algorithm calculate current nominal fatigue vector d i cnoninferior solution, namely can reflect the solution of the change of the health status of evaluation object more exactly with reasonable error.
The Objective Programming in multi-objective optimization algorithm (Goal Attainment Method) can be adopted to solve current nominal fatigue vector d i c, according to Objective Programming, formula (2) can transform the multi-objective optimization question shown in an accepted way of doing sth (3) and formula (4), and in formula (3), γ is a real number, and R is real number field, and area of space Ω limits vectorial d i cspan (the present embodiment requirements vector d of each element i ceach element be not less than 0, be not more than 1).Formula (3) be meant to the minimum real number γ of searching one, formula (4) is met.G (d in formula (4) i c) defined by formula (5), the middle G (d of the product representation formula (4) of weighing vector W and γ in formula (35) i c) and vectorial g ibetween allow deviation.During actual computation vector W can with vectorial g iidentical.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of current nominal fatigue vector d i c.
minimize γ
(3)
11 step: according to evaluation object current actual damage vector d idefinition and the definition of its element calculate current actual damage vector d ieach element, thus can by d idetermine the health status of evaluation object.Current actual damage vector d ia kth element d i krepresent the current actual health status of a kth evaluation object in i-th circulation.
D i krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system i krepresent its current actual damage, d i krepresent when being 0 that the support cable of its correspondence is without health problem, d i knumerical value represents when not being 0 that the support cable of its correspondence 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.
D i krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is generalized displacement component, so a d of a bearing i krepresent its current actual generalized displacement numerical value.
D i krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is load, so a d i krepresent its current real load change numerical value.
By current for evaluation object actual damage vector d iin the M relevant to support cable 1individual element takes out, composition support cable current actual damage vector d ci, support cable current actual damage vector d cithe coding rule of element and Initial cable force vector F othe coding rule of element identical.Support cable current actual damage vector d cih element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M 1; Support cable current actual damage vector d cimiddle numerical value be not 0 element correspond to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, 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 support cable that damage, unsoundness problem is exactly there occurs lax rope, and needing the rope adjusting Suo Li exactly, is exactly slack line, and these ropes that need adjust Suo Li are at support cable current actual damage vector d ciin corresponding element numerical value (such as one of them element can use d ci hrepresent) represent the degree of injury with the relax level mechanic equivalent of these support cables, just determine slack line thus, the computing method of concrete slack are described below.The support cable corresponding with damaged cable current actual damage vector d ciin element numerical expression be the current actual damage of this damaged cable, represent when element numerical value is 100% that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; According to support cable current actual damage vector d ci, identify slack line from the support cable of unsoundness problem after, remaining is exactly damaged cable, and these damaged cables are at support cable current actual damage vector d cithe numerical value of the element of middle correspondence just represents its degree of injury, the numerical value of corresponding element represents when being 100% that this support cable thoroughly loses load-bearing capacity, represent time between 0 and 100% that this support cable loses the load-bearing capacity of corresponding proportion, so far just have identified damaged cable and degree of injury thereof.
Utilize support cable current actual damage vector d cithe slack line of expressing with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize and obtain monitored amount current value vector C in the 6th step, in actual measurement isynchronization obtain current cable force vector F iend points horizontal range vector l is supported with current support cable two i x, to utilize in second step, actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while (also can be setting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization) obtain support cable initial elastic modulus vector E o, support cable initial drift vector l o, initial free cross-sectional area vector A owith the weight vector ω of initial free unit length othe physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained at second step to use, 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, particularly can in the hope of the relax level of these ropes (i.e. the long adjustment amount of rope) according to formula (6) or formula (7).So just achieve the lax identification of support cable.So far just damaged cable and slack line is all identified.
Formula (6) and the middle E of formula (7) ohsetting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization, the elastic modulus of h root support cable, A ohsetting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization, the cross-sectional area of h root support cable, ω ohsetting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization, the weight of the unit length of h root support cable, l ohsetting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization, the drift of h root support cable in Cable Structure, E ohvectorial E oh element, be setting up the initial mechanical Calculation Basis model A of Cable Structure osynchronization, the elastic modulus of h root support cable in Cable Structure, F i hobtain monitored amount current value vector C in actual measurement isynchronization, the current cable power of h root support cable, d ci hthe current actual damage degree of h root support cable, l i xhobtain monitored amount current value vector C in actual measurement isynchronization, the horizontal range of two of h root support cable supporting end points, l i xhthat current support cable two supports end points horizontal range vector l i xan element, current support cable two support end points horizontal range vector l i xthe coding rule of element and initial drift vector l othe coding rule of element identical.Item in formula (7) in [] is the Ernst equivalent elastic modulus of this support cable, just can determine support cable current slack degree vector Δ l by formula (6) or formula (7) i.Formula (7) is the correction to formula (6).
So far this method achieves the accurate identification of the health status of core evaluation object with a kind of effective, cheap method.May exact value be departed to the recognition result of the health status of secondary evaluation object more, only require the correct health status identifying core evaluation object in the method.
12 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.
13 step: 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.
14 step: according to method establishment mark vector B in technical scheme iif, mark vector B ielement be 0 entirely, then get back to the 6th step and proceed health monitoring to cable system and calculating; If mark vector B ielement be not 0, then after completing subsequent step entirely, enter and circulate next time.
15 step: according to method in technical scheme calculate next time (namely the i-th+1 time, i=1,2,3,4 ...) needed for circulation initial damage vector d i+1 oeach element d i+1 ok(k=1,2,3 ..., N); The second, at initial mechanical Calculation Basis model A obasis on, make the health status of rope be d i+1 oafter obtain be exactly next time, namely the i-th+1 time (i=1,2,3,4 ...) Mechanics Calculation benchmark model A needed for circulation i+1.Obtain A i+1, d i+1 oafter, obtain A by Mechanics Calculation i+1in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C i+1 o.
16 step: get back to the 6th step, starts by the circulation of the 6th step to the 16 step.

Claims (1)

1. simplify cable force monitoring loading problem rope generalized displacement progressive-type recognition method, it is characterized in that described method comprises:
Though the load of a. bearing when Cable Structure changes, when the load that Cable Structure is being born does not exceed Cable Structure initial allowable load, this method is suitable for; The initial allowable load of Cable Structure refers to the allowable load of Cable Structure when being completed, and can be obtained by conventional Mechanics Calculation; This method unitedly calls evaluated generalized displacement of support component, support cable and load to be evaluation object, if the quantity sum of the evaluated quantity of generalized displacement of support component, the quantity of support cable and load is N, namely the quantity of evaluation object is N; Determine the coding rule of evaluation object, evaluation objects all in Cable Structure numbered by this rule, this numbering will be used for generating vector sum matrix in subsequent step; This method variable k represents this numbering, k=1,2,3 ..., N; This method title " core evaluation object " specially refers to evaluated support cable in " evaluation object " and generalized displacement of support component, if the quantity sum of evaluated support cable and generalized displacement of support component is P, namely the quantity of core evaluation object is P, and this method title " secondary evaluation object " specially refers to the evaluated load in " evaluation object "; If total M in cable system 1root support cable, Cable Structure rope force data comprises this M 1the Suo Li of root support cable, this method is at the whole M of monitoring 1on the basis of root supporting cable force, Cable Structure artificially increases M 2root rope, is called sensing rope, by this M newly increased of monitoring in cable structure health monitoring process 2the Suo Li of root sensing rope; Comprehensive above-mentioned monitored amount, M the Suo Li that whole Cable Structure has M root rope is monitored, and namely have M monitored amount, wherein M is M 1with M 2sum; M should be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; The M newly increased 2the rigidity of root sensing rope is compared with the rigidity of any support cable of Cable Structure, should be much smaller; The M newly increased 2the Suo Li of each sensing rope of root sensing rope should be more much smaller than the Suo Li of any support cable of Cable Structure, even if can ensure this M newly increased like this 2root sensing rope has occurred damaging or lax, very little on the impact of the stress of other components of Cable Structure, strain, distortion; The M newly increased 2on the xsect of root sensing rope, normal stress should be less than its fatigue limit, and these requirements can ensure the M newly increased 2can not fatigue damage be there is in root sensing rope; The M newly increased 2the two ends of root sensing rope should fully anchoring, ensures to there will not be lax; The M newly increased 2root sensing rope should obtain sufficient anti-corrosion protection, ensures the M newly increased 2damage can not be there is and relax in root sensing rope; For simplicity, in the method by " monitored all parameters of Cable Structure " referred to as " monitored amount "; To M monitored amount serial number, this method represents this numbering, j=1,2,3 with variable j ..., M, this numbering will be used for generating vector sum matrix in subsequent step; The M newly increased in the method 2root sensing rope is as a part for Cable Structure, and when hereinafter mentioning Cable Structure again, Cable Structure comprises increases M 2cable Structure before root sensing rope and the M newly increased 2root sensing rope, refers to comprise the M newly increased when that is hereinafter mentioning Cable Structure 2the Cable Structure of root sensing rope; Measure the M obtaining newly increasing 2the method of the Suo Li of root sensing rope is same as the M of Cable Structure 1the measuring method of the Suo Li of root support cable, hands over later no longer one by one; When carrying out any measurement to the support cable of Cable Structure, simultaneously to the M newly increased 2root sensing rope carries out same measurement, hands over no longer one by one later; The M newly increased 2root sensing rope except there is not damage and be lax, to the M newly increased 2the requirement of the quantity of information of root sensing rope is identical with preparation method with the requirement of the quantity of information of the support cable of Cable Structure with preparation method, hands over no longer one by one later; When setting up the various mechanical model of Cable Structure later, by the M newly increased 2the support cable that Cable Structure treated as by root sensing rope is treated; Below, except the damage of mentioning support cable and lax occasion, the support cable that said support cable comprises Cable Structure when mentioning support cable and the M newly increased 2root sensing rope; The external force that object, structure are born can be described as load, and load comprises face load and volume load; Face load, also known as surface load, is the load acting on body surface, comprises centre-point load and distributed load two kinds; Volume load be continuous distribution in the load of interior of articles each point, comprise deadweight and the inertial force of object; Centre-point load is divided into concentrated force and concentrated couple two kinds, tie up in interior coordinate system comprising Descartes's rectangular coordinate, a concentrated force can resolve into three components, same, a concentrated couple also can resolve into three components, if load is actually centre-point load, in the method a concentrated force component or a concentrated couple component being counted or added up is a load, and the now change of load is embodied as the change of a concentrated force component or a concentrated couple component; Distributed load is divided into line distributed load and EDS maps load, and the description of distributed load at least comprises the zone of action of distributed load and the size of distributed load, and the size distribution intensity of distributed load is expressed, and distribution intensity distribution characteristics and amplitude are expressed; If load is actually distributed load, when this method talks about the change of load, in fact refer to the change of the amplitude of distributed load distribution intensity, and the distribution characteristics of the zone of action of all distributed loads and distribution intensity is constant; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, a distributed load can resolve into three components, if the amplitude of the respective distribution intensity of three of this distributed load components changes, and the ratio of change is all not identical, so in the method three of this distributed load components being counted or added up is three distributed loads, and now load just represents the one-component of distributed load; Volume load be continuous distribution in the load of interior of articles each point, the description of volume load at least comprises the zone of action of volume load and the size of volume load, and the size distribution intensity of volume load is expressed, distribution intensity distribution characteristics and amplitude express; If load is actually volume load, actual treatment is the change of the amplitude of volume load diatibution intensity in the method, and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant, in fact the change of the amplitude of the distribution intensity of volume load is referred to when now mentioning the change of load in the method, now, the load changed refers to the volume load that the amplitude of those distribution intensities changes; Tie up in interior coordinate system comprising Descartes's rectangular coordinate, one individual stow lotus can resolve into three components, if the amplitude of the respective distribution intensity of three of this volume load components changes, and the ratio of change is all not identical, so in the method three of this volume load components being counted or added up is three distributed loads;
B. survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses; Actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the Initial cable force of all support cables, composition Initial cable force vector F o; According to comprising Cable Structure design data, the data of completion data obtain the length of all support cables when free state and Suo Li are 0, in free state time cross-sectional area and in free state time the weight of unit length, 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;
C. actual measurement or consult reference materials obtain the various materials that Cable Structure uses physical and mechanical properties parameter while, direct survey calculation obtains the measured data of initial Cable Structure, and the measured data of initial Cable Structure comprises Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, the initial generalized displacement measurement data of Cable Structure bearing, the initial value of all monitored amounts, the Initial cable force data of all support cables, initial Cable Structure modal data, initial Cable Structure strain data, initial Cable Structure geometric data, initial Cable Structure bearing generalized coordinate data, initial Cable Structure angle-data, initial Cable Structure spatial data is in interior measured data, while the measured data obtaining initial Cable Structure, survey calculation obtains the data can expressing the health status of support cable of the Non-destructive Testing Data comprising support cable, and the data can expressing the health status of support cable are now called support cable initial health data, the initial value of all monitored amounts forms monitored amount initial value vector C o, monitored amount initial value vector C othe coding rule of coding rule and M monitored amount identical, the initial generalized displacement measurement data of support cable initial health data, Cable Structure bearing and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d o, vectorial d orepresent with initial mechanical Calculation Basis model A othe initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d oelement number equal N, d oelement and evaluation object be one-to-one relationship, vectorial d othe coding rule of element identical with the coding rule of evaluation object, if d oevaluation object corresponding to some elements be support cable, so a d in cable system othe numerical value of this element represent the initial damage degree of corresponding support cable, if the numerical value of this element is 0, represent that the support cable corresponding to this element is intact, do not damage, if its numerical value is 100%, then represent that the support cable corresponding to this element completely loses 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 d oevaluation object corresponding to some elements be some generalized displacement components of some bearings, so d othe numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d oevaluation object corresponding to some elements be some load, get d in this method othis element numerical value be 0, the initial value representing the change of this load is 0, if when there is no Cable Structure bearing initial generalized displacement measurement data or can think that the initial generalized displacement of Cable Structure bearing is 0, vectorial d oin each element numerical value relevant to Cable Structure generalized displacement of support get 0, 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 can think structure original state be not damaged without relaxed state time, vectorial d oin each element numerical value relevant to support cable get 0, initial Cable Structure bearing generalized coordinate data refer to the bearing generalized coordinate data under Cable Structure design point, and Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A otime, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point,
All Cable Structure data that the physical and mechanical properties parameter 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, support cable initial health data, Cable Structure bearing initial generalized displacement measurement data, Cable Structure centre-point load measurement data, Cable Structure distributed load measurement data, Cable Structure volume load measurement data, Cable Structure and preceding step obtain, set up the initial mechanical Calculation Basis model A of Cable Structure 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 oevaluation object health status with evaluation object initial damage vector d orepresent; Corresponding to A othe initial value monitored amount initial value vector C of all monitored amount orepresent; d oa oparameter, 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, A in the method o, C oand d oconstant;
E. in the method, alphabetical i is except representing the place of number of steps significantly, and alphabetical i only represents cycle index, i.e. i-th circulation; I-th circulation needs the current initial mechanical Calculation Basis model of Cable Structure that is that set up or that set up to be designated as current initial mechanical Calculation Basis model A when starting i o; The current initial damage vector of evaluation object that i-th circulation needs when starting is designated as d i o, d i ocable Structure A when representing that this circulation starts i othe health status of evaluation object, d i odefinition mode and d odefinition mode identical, d i oelement and d oelement one_to_one corresponding; When i-th circulation starts, the initial value of all monitored amounts, with monitored amount current initial value vector C i orepresent, vectorial C i odefinition mode and vectorial C odefinition mode identical, C i oelement and C oelement one_to_one corresponding, monitored amount current initial value vector C i orepresent and correspond to A i othe concrete numerical value of all monitored amount; d i oa i ocharacterisitic parameter, C i oby A i omechanics Calculation result composition; When first time, circulation started, A i obe designated as A 1 o, set up A 1 omethod for making A 1 oequal A o; When first time, circulation started, d i obe designated as d 1 o, set up d 1 omethod for making d 1 oequal d o; When first time, circulation started, C i obe designated as C 1 o, set up C 1 omethod for making C 1 oequal C o;
F. from entering the circulation being walked to p step by f here;
G. in structure military service process, actual measurement obtains the currency of all monitored amounts in Cable Structure, and all these numerical value form monitored amount current value vector C i, vectorial C idefinition mode and vectorial C odefinition mode identical, C ielement and C oelement one_to_one corresponding, represent that identical monitored amount is at not numerical value in the same time, monitored amount current value vector C is obtained in actual measurement iwhile, to the M newly increased 2root sensing rope carries out Non-Destructive Testing, therefrom identify and occur damage or lax sensing rope, according to monitored amount coding rule, the element that the sensing rope that to damage with the appearance identified or relax is corresponding is removed from each vector according to monitored amount coding rule numbering occurred before this method, also the element damaged with the appearance identified or lax sensing rope is corresponding is no longer there is in each vector sum matrix occurred after this method, no longer comprise when mentioning sensing rope after this method being identified here and occur damage or lax sensing rope, the Suo Li being identified and occurring damage or lax sensing rope is no longer comprised here when mentioning monitored amount after this method, identify several from Cable Structure and occur damage or lax sensing rope, just by M 2same quantity is reduced with M, monitored amount current value vector C is obtained in actual measurement isynchronization, actual measurement obtain all M in Cable Structure 1the rope force data of root support cable, all these rope force datas composition current cable force vector F i, vectorial F ielement and vectorial F othe coding rule of element identical, monitored amount current value vector C is obtained in actual measurement isynchronization, Actual measurement obtains all M 1the volume coordinate of two supporting end points of root support cable, 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,
H. at current initial mechanical Calculation Basis model A i obasis on, carry out several times Mechanics Calculation according to step h1 to step h4, by calculate set up unit damage monitored numerical quantity transformation matrices Δ C iwith evaluation object unit change vector D i u;
H1., when i-th circulation starts, directly Δ C is obtained by method listed by step h2 to step h4 iand D i u; In other moment, when in step g to A i oafter upgrading, Δ C must be regained by method listed by step h2 to step h4 iand D i uif, not to A in step g i oupgrade, then directly proceed to step I herein and carry out follow-up work;
H2. at current initial mechanical Calculation Basis model A i obasis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity N of all evaluation objects, has N number of evaluation object just to have N calculating, according to the coding rule of evaluation object, calculate successively, calculating hypothesis each time only has an evaluation object on the basis of original damage or generalized displacement or load, increase unit damage or unit generalized displacement or load unit change again, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable increases unit damage again, if this evaluation object is the generalized displacement component in a direction of a bearing, just suppose that this bearing increases unit generalized displacement again at this sense of displacement, if this evaluation object is a load, just suppose that this load increases load unit change again, use D i ukrecord the unit damage of this increase or unit generalized displacement or load unit change, wherein k represents the numbering of the evaluation object increasing unit damage or unit generalized displacement or load unit change, D i ukevaluation object unit change vector D i uan element, evaluation object unit change vector D i uthe coding rule of element and vectorial d othe coding rule of element identical, the evaluation object increasing unit damage or unit generalized displacement or load unit change in calculating each time is again different from during other time calculates the evaluation object increasing unit damage or unit generalized displacement or load unit change again, calculate the current calculated value all utilizing mechanics 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, when suppose a kth evaluation object increase again unit damage or unit generalized displacement or load unit change time, use C i tkrepresent corresponding " monitored amount calculation current vector ", when giving the element number of each vector in this step, same coding rule should be used with other vector in this method, to ensure any one element in this step in each vector, with in other vector, number identical element, have expressed the relevant information of same monitored amount or same target, C i tkdefinition mode and vectorial C odefinition mode identical, C i tkelement and C oelement one_to_one corresponding,
H3. the vectorial C calculated each time i tkdeduct vectorial C i oobtain a vector, then obtain " numerical value change vector δ a C for monitored amount after each element of this vector is calculated divided by this unit damage or unit generalized displacement or load unit change numerical value supposed i k"; N number of evaluation object is had just to have N number of " the numerical value change vector of monitored amount ";
H4. by this is N number of " the numerical value change vector of monitored amount " according to the coding rule of N number of evaluation object, " the unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively i"; Unit damage monitored numerical quantity transformation matrices Δ C ieach row correspond to a monitored amount unit change vector; Unit damage monitored numerical quantity transformation matrices Δ C ievery a line correspond to the different unit change amplitude of same monitored amount when different evaluation object increases unit damage or unit generalized displacement or load unit change; Unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of row and vectorial d othe coding rule of element identical, unit damage monitored numerical quantity transformation matrices Δ C ithe coding rule of coding rule and M monitored amount of row identical;
I. current nominal fatigue vector d is defined i cwith current actual damage vector d i, d i cand d ielement number equal the quantity of evaluation object, d i cand d ielement and evaluation object between be one-to-one relationship, d i celement numerical value represent the nominal fatigue degree of corresponding evaluation object or nominal generalized displacement or nominal load variable quantity, d i cand d iwith evaluation object initial damage vector d oelement number rule identical, d i celement, d ielement and d oelement be one-to-one relationship;
J. according to monitored amount current value vector C iwith " monitored amount current initial value vector C i o", " unit damage monitored numerical quantity transformation matrices Δ C i" and " current nominal fatigue vector d i c" between the linear approximate relationship that exists, this linear approximate relationship can be expressed as formula 1, except d in formula 1 i cother outer amount is known, solves formula 1 and just can calculate current nominal fatigue vector d i c;
C i = C o i + ΔC i · d c i Formula 1
K. the current actual damage vector d utilizing formula 2 to express ia kth element d i kwith evaluation object current initial damage vector d i oa kth element d i okwith current nominal fatigue vector d i ca kth element d i ckbetween relation, calculate current actual damage vector d iall elements;
formula 2
K=1 in formula 2,2,3 ..., N; d i krepresent the current actual health status of a kth evaluation object in i-th circulation, if this evaluation object is support cable, so a d in cable system i krepresent the order of severity of its current health problem, the support cable of unsoundness problem may be slack line, also may be damaged cable, d i kthe degree of the lax or damage of this support cable of numerical response; From the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line, evaluation object current actual damage vector d iin correspond to slack line element numerical expression be the current actual equivalent damage degree with slack line relax level mechanic equivalent; If this evaluation object is load, so a d i krepresent the actual change amount of this load; If this evaluation object is generalized displacement component, so a d of a bearing i krepresent its current actual generalized displacement numerical value; By current for evaluation object actual damage vector d iin with M 1the M that root support cable is relevant 1individual element takes out, composition support cable current actual damage vector d ci, support cable current actual damage vector d cithe coding rule of element and Initial cable force vector F othe coding rule of element identical; Support cable current actual damage vector d cih element representation Cable Structure in the current actual damage amount of h root support cable, h=1,2,3 ...., M 1; Support cable current actual damage vector d cimiddle numerical value be not 0 element correspond to the support cable of unsoundness problem, from the support cable of these unsoundness problems, identify damaged cable, remaining is exactly slack line; The support cable corresponding with damaged cable current actual damage vector d ciin element numerical expression be the current actual damage of this damaged cable, represent when element numerical value is 100% that this support cable thoroughly loses load-bearing capacity, time between 0 and 100%, represent the load-bearing capacity losing corresponding proportion; After identifying slack line, utilize support cable current actual damage vector d cithese slack lines of expressing, with the current actual equivalent damage degree of its relax level mechanic equivalent, utilize current cable force vector F isupport end points horizontal range vector with current support cable two, utilize initial drift vector, the weight vector of the initial free unit length of initial free cross-sectional area vector sum, the vectorial F of Initial cable force of the support cable obtained in b step othe physical and mechanical properties parameter of the various materials utilizing the Cable Structure obtained in b step to use, 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, mechanic equivalent condition 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 mechanic 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 by current cable force vector F icorresponding element provides; Damaged cable and slack line are referred to as the support cable of unsoundness problem by this method, referred to as problem cable, so far this method achieve reject generalized displacement of support, load change impact, the problem cable identification of Cable Structure, achieve simultaneously reject generalized displacement of support and support cable health status variable effect, the identification of load change amount, also achieve reject load change and support cable health status variable effect, the identification of generalized displacement of support; So far this method achieves the accurate identification of the health status of core evaluation object with a kind of effective, cheap method; Exact value is departed to the recognition result of the health status of secondary evaluation object more, therefore will not accept and believe, only require the correct health status identifying core evaluation object in the method;
L. current nominal fatigue vector d is tried to achieve i cafter, set up mark vector B according to formula 3 i, formula 4 gives mark vector B ithe definition of a kth element;
B i = B 1 i B 2 i · · · B k i · · · B N i T Formula 3
Element B in formula 4 i kmark vector B ia kth element, D i ukevaluation object unit change vector D i ua kth element, d i ckevaluation object current nominal fatigue vector d i ca kth element, they all represent the relevant information of a kth evaluation object, k=1 in formula 4,2,3 ..., N;
If m. mark vector B ielement be 0 entirely, then get back to step f continue this circulation; If mark vector B ielement be not 0 entirely, then enter next step, i.e. step n;
N. calculate next time according to formula 5, evaluation object current initial damage vector d namely needed for the i-th+1 time circulation i+1 oeach element;
formula 5
D in formula 5 i+1 okthe evaluation object current initial damage vector d next time, namely needed for the i-th+1 time circulation i+1 oa kth element, d i okthis, i.e. the evaluation object of i-th circulation current initial damage vector d i oa kth element, D i ukthe evaluation object unit change vector D of i-th circulation i ua kth element, B i kthe mark vector B of i-th circulation ia kth element, k=1 in formula 5,2,3 ..., N;
O. at initial mechanical Calculation Basis model A obasis on, make the health status of rope be d i+1 oafter obtain be exactly next time, namely the i-th+1 time circulation needed for Mechanics Calculation benchmark model A i+1; Obtain A i+1after, obtain A by Mechanics Calculation i+1in all monitored amounts, current concrete numerical value, these the monitored amounts of concrete numerical value composition next time, namely needed for the i-th+1 time circulation current initial value vector C i+1 o;
P. get back to step f, start to circulate next time.
CN201510437824.XA 2015-07-23 2015-07-23 Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring Pending CN104990735A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510437824.XA CN104990735A (en) 2015-07-23 2015-07-23 Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510437824.XA CN104990735A (en) 2015-07-23 2015-07-23 Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring

Publications (1)

Publication Number Publication Date
CN104990735A true CN104990735A (en) 2015-10-21

Family

ID=54302577

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510437824.XA Pending CN104990735A (en) 2015-07-23 2015-07-23 Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring

Country Status (1)

Country Link
CN (1) CN104990735A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05306971A (en) * 1992-04-30 1993-11-19 Fujikura Ltd Cable behavior testing device
CN102221478A (en) * 2011-05-31 2011-10-19 东南大学 Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support
CN102323084A (en) * 2011-05-31 2012-01-18 东南大学 Progressive method for identifying damaged cable and support generalized displacement on the basis of mixed monitoring
CN103868744A (en) * 2014-03-10 2014-06-18 东南大学 Progressing identification method for problem cable load generalized displacement through cable force monitoring
CN103913328A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement hybrid monitoring damaged cable load progressive identification method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05306971A (en) * 1992-04-30 1993-11-19 Fujikura Ltd Cable behavior testing device
CN102221478A (en) * 2011-05-31 2011-10-19 东南大学 Progressive health monitoring method for cable system based on hybrid monitoring during generalized displacement of support
CN102323084A (en) * 2011-05-31 2012-01-18 东南大学 Progressive method for identifying damaged cable and support generalized displacement on the basis of mixed monitoring
CN103868744A (en) * 2014-03-10 2014-06-18 东南大学 Progressing identification method for problem cable load generalized displacement through cable force monitoring
CN103913328A (en) * 2014-03-10 2014-07-09 东南大学 Generalized displacement hybrid monitoring damaged cable load progressive identification method

Similar Documents

Publication Publication Date Title
CN104990735A (en) Simplified load problematical cable generalized displacement progressive-identification method through cable force monitoring
CN104990736A (en) Simplified load problematical cable generalized displacement progressive-identification method through strain monitoring
CN104990722A (en) Simplified load problematical cable generalized displacement progressive-identification method through space coordinate monitoring
CN105115746A (en) Progressive identification method for problematic cable and generalized displacement based on angle monitoring through simplified load change
CN105067296A (en) Simplified generalized displacement cable-force-monitoring load defective cable progressive identifying method
CN105115770A (en) Progressive identification method for load, problematic cable and generalized displacement through simplified hybrid monitoring
CN104990740A (en) Simplified load problematical cable progressive-identification method through strain monitoring under generalized displacement
CN105004546A (en) Method for progressively recognizing damaged cable of load based on streamlined generalized displacement and angle monitoring process
CN105067351A (en) Method for progressively recognizing damaged cable of load based on space coordinate monitoring process of streamlined and generalized displacement
CN105067345A (en) Method for recognizing damaged cable generalized displacement of load based on streamlined cable force monitoring process
CN104990587A (en) Simplified method of recognizing problematic load cable generalized displacement through strain monitoring
CN105067317A (en) Method for recognizing damaged cable generalized displacement of load based on streamlined hybrid monitoring process
CN104990733A (en) Simplified load problematical cable generalized displacement identification method through space coordinate monitoring
CN105115767A (en) Progressive identification method for load, damaged cable and generalized displacement through simplified cable force monitoring
CN105136501A (en) Simplified angle monitoring load faulty cable generalized displacement identification method
CN105067340A (en) Method for recognizing damaged cable of load based on space coordinate monitoring process of streamlined and generalized displacement
CN104990724A (en) Simplified load damaged cable generalized displacement progressive-identification method through strain monitoring
CN105067299A (en) Simplified generalized displacement cable-force-monitoring load damaged cable progressive identifying method
CN105067309A (en) Simplified cable-force-monitoring load damaged cable generalized displacement identifying method
CN104990577A (en) Simplified load damaged cable generalized displacement progressive-identification method through hybrid monitoring
CN105115771A (en) Progressive identification method for load, damaged cable and generalized displacement through simplified space coordinate monitoring
CN105067331A (en) Method for recognizing damaged cable of load based on cable force monitoring process of streamlined and generalized displacement
CN105136504A (en) Simplified generalized displacement strain monitoring load faulty cable progressive identification method
CN105067301A (en) Simplified strain-monitoring load damaged cable generalized displacement identifying method
CN105115756A (en) Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20151021