CN105115756A  Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change  Google Patents
Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change Download PDFInfo
 Publication number
 CN105115756A CN105115756A CN201510438096.4A CN201510438096A CN105115756A CN 105115756 A CN105115756 A CN 105115756A CN 201510438096 A CN201510438096 A CN 201510438096A CN 105115756 A CN105115756 A CN 105115756A
 Authority
 CN
 China
 Prior art keywords
 load
 cable
 evaluation
 vector
 initial
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Pending
Links
 238000006073 displacement reaction Methods 0.000 title claims abstract description 167
 238000004364 calculation method Methods 0.000 claims abstract description 76
 239000011159 matrix material Substances 0.000 claims abstract description 12
 238000011156 evaluation Methods 0.000 claims description 271
 238000009826 distribution Methods 0.000 claims description 73
 238000005259 measurement Methods 0.000 claims description 63
 230000003862 health status Effects 0.000 claims description 54
 230000000875 corresponding Effects 0.000 claims description 43
 230000001131 transforming Effects 0.000 claims description 27
 239000000463 material Substances 0.000 claims description 13
 238000009659 nondestructive testing Methods 0.000 claims description 12
 206010022114 Injury Diseases 0.000 claims description 9
 239000000203 mixture Substances 0.000 claims description 7
 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 6
 240000002853 Nelumbo nucifera Species 0.000 claims description 6
 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 6
 230000000694 effects Effects 0.000 claims description 5
 230000001771 impaired Effects 0.000 claims description 5
 238000000547 structure data Methods 0.000 claims description 2
 239000002965 rope Substances 0.000 description 20
 238000007796 conventional method Methods 0.000 description 9
 238000000034 method Methods 0.000 description 9
 238000005457 optimization Methods 0.000 description 9
 238000004891 communication Methods 0.000 description 3
 239000000725 suspension Substances 0.000 description 3
 230000005540 biological transmission Effects 0.000 description 2
 238000007906 compression Methods 0.000 description 2
 238000005516 engineering process Methods 0.000 description 2
 238000003860 storage Methods 0.000 description 2
 241001269238 Data Species 0.000 description 1
 206010043431 Thinking abnormal Diseases 0.000 description 1
 230000004075 alteration Effects 0.000 description 1
 238000004458 analytical method Methods 0.000 description 1
 230000015572 biosynthetic process Effects 0.000 description 1
 238000005755 formation reaction Methods 0.000 description 1
 230000002068 genetic Effects 0.000 description 1
 230000005484 gravity Effects 0.000 description 1
 239000004615 ingredient Substances 0.000 description 1
 238000009434 installation Methods 0.000 description 1
 230000001537 neural Effects 0.000 description 1
 239000000700 tracer Substances 0.000 description 1
 230000017105 transposition Effects 0.000 description 1
 238000005303 weighing Methods 0.000 description 1
Abstract
An identification method for a damaged cable and generalized displacement based on angle monitoring through simplified load change establishes a mechanical calculation reference model of a cable structure based on the angle monitoring, and calculates and obtains a unit damage monitored quantity valuevariable matrix based on the model. Noninferior solution of a current nominal damage vector of an evaluated object is calculated according to an approximate linear relation existing in a current value vector of monitored quantity, an initial value vector of the monitored quantity, the unit damage monitored quantity valuevariable matrix, and the current nominal damage vector of the estimated object to be calculated. Hereby a health state of a core estimated object can be identified.
Description
Technical field
Cablestayed bridge, suspension bridge, the structures such as trussframe structure have a common ground, be exactly that they have many parts bearing tensile load, as suspension cable, main pushtowing 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 carryingropes 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 carryingrope 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 trussframe structure reality.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 support cable is impaired, 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 Cable Structure 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 damaged cable based on angle monitor (monitored angle 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.
Technical scheme: this method is made up of three parts.Method, knowledge based storehouse (containing parameter) and actual measurement the structural health conditions appraisal procedure of monitored amount, the software and hardware part of health monitoring systems of setting up knowledge base needed for structural healthy monitoring system and parameter respectively.
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 centrepoint 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.
Centrepoint 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 centrepoint 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 welldistributed pressure be example so that the concept of amplitude to be described: same structure bears two different welldistributed 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 centrepoint load, in the method, " load unit change " in fact refers to " unit change of centrepoint load ", similar, " load change " specifically refers to " change of the size of centrepoint load ", " load change amount " specifically refers to " variable quantity of the size of centrepoint load ", " load change degree " specifically refers to " intensity of variation of the size of centrepoint load ", " the actual change amount of load " refers to " the actual change amount of the size of centrepoint load ", " load changed " refers to " centrepoint load that size changes ", briefly, now " soandso load soandso change " refers to " size of soandso centrepoint load soandso 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 " soandso load soandso change " refers to " amplitude of the distribution intensity of soandso distributed load soandso 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, " soandso load soandso change " refers to " amplitude of the distribution intensity of soandso volume load soandso change ", and the distribution characteristics of the zone of action of all volume load and distribution intensity is constant.
First the quantity of the load that may change that Cable Structure is born is confirmed.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}.
" the whole monitored angledata of structure " by the specified point of K in structure, cross L of each specified point appointment straight line, H angle coordinate component of each appointment straight line describe, the change of structural point is exactly the change of all angle coordinate components of specifying of all specified points, all appointment straight line.Each total M (M=K × L × H) individual angle coordinate component measurement value or calculated value carry out the angle information of characterisation of structures.
Comprehensive abovementioned monitored amount, whole Cable Structure has M monitored amount, and M should be greater than the quantity of core evaluation object, and M is less than the quantity of 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.
The Part I of this method: the method setting up knowledge base needed for structural healthy monitoring system and parameter.Specific as follows:
1. set up the initial mechanical Calculation Basis model A of Cable Structure
_{o}the method of (such as finite element benchmark model), sets up and A
_{o}corresponding monitored amount initial value vector C
_{o}method.Monitored amount initial value vector C
_{o}the coding rule of coding rule and M monitored amount identical.Set up A
_{o}and C
_{o}method as follows.
When Cable Structure is completed, or before setting up structural healthy monitoring system, the direct survey calculation of conventional method is used to obtain the initial value of all monitored amount of Cable Structure, 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 simultaneously, conventional method Actual measurement also must be used to obtain the Actual measurement data of Cable Structure.First the Actual measurement data of Cable Structure are the data can expressing the health status of rope of the Nondestructive Testing Data comprising support cable, and the Actual measurement data of Cable Structure still comprise the measured data of Cable Structure bearing initial generalized displacement measurement data, the initial geometric data of Cable Structure, rope force data, drawbar pull data, initial Cable Structure bearing generalized coordinate data, Cable Structure modal data, structural strain data, structural point measurement data, structure space measurement of coordinates data and load data.The initial geometric data of Cable Structure can be the spatial data that the spatial data of the end points of all ropes adds a series of point in structure, and object is the geometric properties according to these coordinate data determination Cable Structure.Cable Structure bearing initial generalized displacement measurement data refers to setting up initial mechanical Calculation Basis model A
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point.For cablestayed 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, socalled bridge type data that Here it is.The variable quantity of " load likely changed " is setting up initial mechanical Calculation Basis model A
_{o}time be all 0, the variable quantity of " load likely changed " that is identifies below is relative to setting up initial mechanical Calculation Basis model A
_{o}time the structure variable quantity of corresponding load that bears.Utilize the Nondestructive Testing Data etc. of support cable can express the data of the health status of support cable, the variable quantity data of the initial generalized displacement measurement data of Cable Structure bearing and " load likely changed " set up evaluation object initial damage vector d
_{o}(such as formula (1) Suo Shi), uses d
_{o}represent that Cable Structure is (with initial mechanical Calculation Basis model A
_{o}represent) the initial health of evaluation object.If there is no the Nondestructive 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
_{o}in each element numerical value relevant to support cable get 0.Vector d
_{o}in each element numerical value relevant to the variable quantity of load 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
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0.The physical and mechanical properties parameter of the various materials utilizing the Nondestructive Testing Data of the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, Cable Structure bearing initial generalized displacement measurement data, support cable and Cable Structure to use, utilizes mechanics method (such as finite element method) to set up initial mechanical Calculation Basis model A
_{o}.
d
_{o}＝[d
_{o1}d
_{o2}···d
_{ok}···d
_{oN}]
^{T}(1)
D in formula (1)
_{ok}(k=1,2,3 ...., N) represent initial mechanical Calculation Basis model A
_{o}in the original state of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d
_{ok}represent its initial damage, d
_{ok}represent not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; If this evaluation object is generalized displacement component, so a d of a bearing
_{ok}represent its initial displacement numerical value; If this evaluation object is one " load that may change ", so d
_{ok}represent its initial value, d
_{ok}be 0, the variable quantity of " load likely changed " that is identifies below is relative to setting up initial mechanical Calculation Basis model A
_{o}time the structure variable quantity of corresponding load that bears.Subscript T represents the transposition (same afterwards) of vector.
At acquisition A
_{o}while, the initial value of all monitored amount of the Cable Structure using the direct survey calculation of conventional method to obtain, forms monitored amount initial value vector C
_{o}(see formula (2)).Monitored amount initial value vector C
_{o}represent and correspond to A
_{o}the concrete numerical value of " monitored amount ".Because of subject to the foregoing, the Calculation Basis model based on Cable Structure calculates the monitored amount of gained reliably close to the measured data of initial monitored amount, in describing below, will represent this calculated value and measured value with prosign.
C
_{o}＝[C
_{o1}C
_{o2}···C
_{oj}···C
_{oM}]
^{T}(2)
C in formula (2)
_{oj}(j=1,2,3 ...., M) be the original bulk of jth monitored amount in Cable Structure, this component corresponds to a specific jth monitored amount according to coding rule.
No matter which kind of method to obtain initial mechanical Calculation Basis model A by
_{o}, based on A
_{o}the 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
_{o}the 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 angledata, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A
_{o}the health status evaluation object initial damage vector d of middle evaluation object
_{o}represent.Due to based on A
_{o}the initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A
_{o}basis on, carry out Mechanics Calculation obtains, A
_{o}the evaluation of each monitored amount form monitored amount initial value vector C
_{o}.D
_{o}a
_{o}parameter, alternatively C
_{o}by A
_{o}mechanics Calculation result composition.
In Cable Structure, the currency of all monitored amounts forms monitored amount current value vector C (formula (3) is shown in definition).
C＝[C
_{1}C
_{2}···C
_{j}···C
_{M}]
^{T}(3)
C in formula (3)
_{j}(j=1,2,3 ...., M) be the currency of jth monitored amount in Cable Structure, this component C
_{j}according to coding rule and C
_{oj}corresponding to same " monitored amount ".In Cable Structure military service process, constantly actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C.
2. set up the method for Cable Structure unit damage monitored numerical quantity transformation matrices Δ C.
The concrete grammar setting up Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is as follows:
At the initial mechanical Calculation Basis model A of Cable Structure
_{o}basis on carry out several times calculating, calculation times numerically equals the quantity of all evaluation objects.Calculating hypothesis each time only has an evaluation object (to use vectorial d in original damage or generalized displacement or load
_{o}corresponding element represent) basis on increase again unit generalized displacement, unit damage or load unit change, concrete, if this evaluation object is a support cable in cable system, so just suppose that this support cable has unit damage (such as getting 5%, 10%, 20% or 30% equivalent damage is unit damage), if this evaluation object is a load, just suppose that this load is at vectorial d
_{o}the 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
^{2}deng be unit change; If this load is centrepoint load, and this centrepoint load is couple, and load unit change can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centrepoint load, and this centrepoint 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
^{3}deng be unit change), 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
_{o}there 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 xaxis of a bearing, just suppose that this bearing is around the angular displacement of xZhou You unit, such as get 100,000/radian), use D
_{uk}record this unit generalized displacement, unit damage or load unit change, wherein k represents the numbering that the evaluation object that unit generalized displacement, unit damage or load unit change occurs.With " evaluation object unit change vector D
_{u}" (such as formula (4) Suo Shi) record all unit generalized displacements, unit damage or load unit change.Unit generalized displacement is there is in calculating each time, the evaluation object of unit damage or load unit change is different from during other time calculates and occurs unit generalized displacement, the evaluation object of unit damage or load unit change, calculate the current calculated value all utilizing mechanics method (such as finite element method) to calculate all monitored amount of Cable Structure each time, the current calculated value of all monitored amount calculated each time forms a monitored amount calculation current vector (when a hypothesis kth evaluation object has unit generalized displacement, when unit damage or load unit change, available formula (5) represents monitored amount calculation current vector C
_{t} ^{k}), calculate monitored amount calculation current vector C each time
_{t} ^{k}deduct monitored amount current initial value vector C
^{t} _{o}after again divided by this time calculate suppose unit generalized displacement, unit damage or load unit change numerical value D
_{uk}gained vector be exactly under this condition (with have unit generalized displacement, unit damage or load unit change evaluation object be numbered mark) monitored amount unit change vector (when a kth evaluation object has unit generalized displacement, unit damage or load unit to change, use δ C
_{k}represent monitored amount unit change vector, formula (6) is shown in definition), the Unit alteration amount of the monitored amount corresponding to this element that each element representation of monitored amount unit change vector causes owing to suppose there is unit generalized displacement, the unit generalized displacement of that evaluation object of unit damage or load unit change, unit damage or load unit change when calculating, N number of evaluation object is had just to have N number of monitored amount unit change vector, owing to there being M monitored amount, so each monitored amount unit change vector has M element, the definition of monitored amount unit change matrix Δ C, the Δ C that have M × N number of element is made up of successively such as formula shown in (6) this N number of monitored amount unit change vector.
D
_{u}＝[D
_{u1}D
_{u2}···D
_{uk}···D
_{uN}]
^{T}(4)
Evaluation object unit change vector D in formula (4)
_{u}element D
_{uk}(k=1,2,3 ...., N) represent the unit generalized displacement of a kth evaluation object, unit damage or load unit change numerical value.
Elements C in formula (5)
_{tj} ^{k}(k=1,2,3 ...., N; J=1,2,3 ...., M) represent due to a kth evaluation object have unit generalized displacement, unit damage or load unit change time, according to the current calculated amount of the individual monitored amount of the jth corresponding to coding rule.
Δ C in formula (7)
_{j,k}(k=1,2,3 ...., N; J=1,2,3,., M) represent only due to a kth evaluation object have unit generalized displacement, unit damage or load unit change to cause, according to the unit change (algebraic value) of the calculating current value of the monitored amount of the jth corresponding to coding rule, the vectorial δ C of monitored amount unit change
_{k}be actually the row in matrix Δ C.
3. monitored amount current value vector C (calculating or actual measurement) is with monitored amount initial value vector C
_{o}, linear approximate relationship between unit damage monitored numerical quantity transformation matrices Δ C and evaluation object current nominal fatigue vector d, shown in (8) or formula (9).The definition of evaluation object current nominal fatigue vector d is see formula (10).
C＝C
_{o}+ΔC·d(8)
CC
_{o}＝ΔC·d(9)
d＝[d
_{1}d
_{2}···d
_{k}···d
_{N}]
^{T}(10)
D in formula (10)
_{k}(k=1,2,3 ...., N) be the current health state of a kth evaluation object in Cable Structure, if this evaluation object is a rope (or pull bar) in cable system, so d
_{k}represent its current damage, d
_{k}represent not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, represent the loadbearing capacity losing corresponding proportion time between 0 and 100%, if this evaluation object is generalized displacement component, so a d of a bearing
_{k}represent its current shift value, if this evaluation object is load, so a d
_{k}represent its variable quantity.
The error of the linear relationship error vector e expression (8) that available formula (11) defines or the shown linear relationship of formula (9).
e＝abs(ΔC·dC+C
_{o})(11)
In formula (11), abs () is the function that takes absolute value, and takes absolute value to each element of the vector of trying to achieve in bracket.
The Part II of this method: the structural health conditions appraisal procedure of knowledge based storehouse (containing parameter) and the monitored amount of actual measurement.
There is certain error in the linear relationship represented by formula (8) or formula (9), therefore simply can not carry out direct solution according to formula (8) or formula (9) and actual measurement monitored amount current value vector C and obtain evaluation object current nominal fatigue vector d.If this has been doned, the element in the evaluation object obtained current nominal fatigue vector d even there will be larger negative value, namely negative damage, and this is obviously irrational.Therefore the acceptable solution of evaluation object current nominal fatigue vector d is obtained (namely with reasonable error, but position and the degree of injury thereof of damaged cable can be determined more exactly from cable system, also can determine load change numerical value more exactly, also can determine generalized displacement of support numerical value more exactly) become a rational solution, available formula (12) expresses this method.
abs(ΔC·dC+C
_{o})≤g(12)
In formula (12), abs () is the function that takes absolute value, and vectorial g describes the legitimate skew departing from ideal linearity relation (formula (8) or formula (9)), is defined by formula (13).
g＝[g
_{1}g
_{2}···g
_{j}···g
_{M}]
^{T}(13)
G in formula (13)
_{j}(j=1,2,3 ...., M) describe the maximum allowable offset departing from formula (8) or the ideal linearity relation shown in formula (9).The error vector e tentative calculation that vector g can define according to formula (11) is selected.
At monitored amount current initial value vector C
^{t} _{o}, unit damage monitored numerical quantity transformation matrices Δ C, survey monitored amount current value vector C known time, suitable algorithm (such as multiobjective optimization algorithm) can be utilized to solve formula (12), obtain the acceptable solution of evaluation object current nominal fatigue vector d.
Definition evaluation object current actual damage vector d
^{a}(see formula (14)), can by d
^{a}determine the health status of evaluation object.
D in formula (14)
^{a} _{k}(k=1,2,3,., N) represent eliminate load change on after the impact of health status recognition result, the current actual health status of a kth evaluation object, if this evaluation object is a rope (or pull bar) in cable system, so d
^{a} _{k}represent its current actual damage, its definition is shown in formula (15), d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this rope thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; If this evaluation object is a generalized displacement component of a bearing, its definition is shown in formula (15), so d
^{a} _{k}represent its current actual generalized displacement numerical value; If this evaluation object is a load, its definition is shown in formula (15), so d
^{a} _{k}represent that it is relative to setting up initial mechanical Calculation Basis model A
_{o}time the structure variable quantity of corresponding load that bears.Vector d
^{a}the coding rule of element and formula (1) in vectorial d
_{o}the coding rule of element identical.
D in formula (15)
_{ok}(k=1,2,3 ...., N) be vectorial d
_{o}a kth element, d
_{k}it is a kth element of vectorial d.
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.
The Part III of this method: the software and hardware part of health monitoring systems.
Hardware components comprises monitoring system (comprising monitored amount monitoring system), signal picker and computing machine etc.Require each monitored amount of RealTime Monitoring.
Software should to complete in this method required, can by functions such as computer implemented monitoring, record, control, storage, calculating, notice, warnings.
This method specifically comprises:
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
_{1}root support cable; Determine the measured point of specifying, number to all specified points; Determined the measured straight line of each measurement point, to all measured straight line numberings of specifying; Determine the measured angle coordinate component of each measured straight line, to all measured angle coordinate component numberings; Abovementioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored angledata of Cable Structure " is made up of abovementioned all measured angle coordinate components; For simplicity, in the method by " the monitored angledata of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M should be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; 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 centrepoint 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; Centrepoint 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 centrepoint 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 onecomponent 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 calculation obtains the measured data of initial Cable Structure, the measured data of initial Cable Structure comprises Cable Structure centrepoint 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 bearing generalized coordinate data, initial Cable Structure angledata and initial Cable Structure spatial data are 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 Nondestructive Testing Data comprising support cable, 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
_{o}the coding rule of coding rule and M monitored amount identical, the initial generalized displacement measurement data of Cable Structure bearing, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the 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 loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some generalized displacement components of some bearings, so d
_{o}the numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this 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
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0, if there is no the Nondestructive 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
_{o}in 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
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point, obtain the measured data of initial Cable Structure at survey calculation while, survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses,
All Cable Structure data that the physical and mechanical properties parameter of the various materials c. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, support cable initial health data, Cable Structure bearing initial generalized displacement measurement data, Cable Structure centrepoint 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
_{o}the Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition;
D. from entering the circulation being walked kth step by d here;
E. at initial mechanical Calculation Basis model A
_{o}basis on carry out several times Mechanics Calculation according to step e1 to e3, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u};
E1. at the initial mechanical Calculation Basis model A of Cable Structure
_{o}basis 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 is at vectorial d
_{o}the basis that this support cable represented has a damage 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 is at vectorial d
_{o}the basis that this load represented has a variable quantity increases load unit change again, use D
_{uk}record 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
_{uk}evaluation object unit change vector D
_{u}an element, evaluation object unit change vector D
_{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or unit generalized displacement or load unit change in calculating each time is different from during other time calculates the evaluation object increasing 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, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical;
E2. the monitored amount calculation current vector calculated each time deducts monitored amount initial value vector C
_{o}obtain a vector, again each element of this vector is calculated the unit damage or unit generalized displacement or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;
E3. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change degree of same monitored amount when different evaluation object increases unit damage or unit generalized displacement or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
F. actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C; Monitored amount current value vector C and monitored amount initial value vector C
_{o}definition mode identical, the same monitored amount of element representation of two vectorial identical numberings is at not concrete numerical value in the same time;
G. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents the nominal fatigue degree of corresponding evaluation object or nominal generalized displacement or nominal load variable quantity; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical;
H. according to monitored amount current value vector C with monitored amount initial value vector C
_{o}, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;
C=C
_{o}+ Δ Cd formula 1
I. evaluation object current actual damage vector d is defined
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent the actual damage degree of corresponding evaluation object or actual generalized displacement or real load variable quantity; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical;
J. the evaluation object utilizing formula 2 to express current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements;
K=1 in formula 2,2,3 ...., N, d
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; If this evaluation object is generalized displacement component, so a d of a bearing
^{a} _{k}represent its current actual generalized displacement numerical value; If this evaluation object is load, so a d
^{a} _{k}represent the actual change amount of this load; So according to evaluation object current actual damage vector d
^{a}impaired and the degree of injury of which support cable can be defined, define which bearing and there occurs generalized displacement and numerical value thereof, define the numerical value which load there occurs change and change thereof; So far this method achieves damaged cable identification that reject the impact of generalized displacement of support and load change, Cable Structure, achieve reject load change and support cable health status variable effect, the identification of generalized displacement of support, achieve reject generalized displacement of support and support cable health status variable effect, the identification of load change amount; 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;
K. get back to d step, start the circulation next time being walked kth step by d.
Beneficial effect: structural healthy monitoring system is first by using sensor to carry out longterm online monitoring to structural response, after obtaining Monitoring Data, (or offline) 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 cablestayed 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}.
" the whole monitored angledata of structure " by the specified point of K in structure, cross L of each specified point appointment straight line, H angle coordinate component of each appointment straight line describe, the change of structural point is exactly the change of all angle coordinate components of specifying of all specified points, all appointment straight line.Each total M (M=K × L × H) individual angle coordinate component measurement value or calculated value carry out the angle information of characterisation of structures.
Comprehensive abovementioned monitored amount, whole Cable Structure has M monitored amount, and the quantity that M must not be less than core evaluation object adds the quantity N 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
_{o}while, 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
_{o}while, use conventional method Actual measurement to obtain the Actual measurement data of Cable Structure.The Nondestructive 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, drawbar 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 centrepoint 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
_{o}time, 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 cablestayed 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, socalled bridge type data that Here it is.Utilize the Nondestructive 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}(such as formula (1) Suo Shi), uses d
_{o}represent that Cable Structure is (with initial mechanical Calculation Basis model A
_{o}represent) the initial health of evaluation object.If there is no the Nondestructive 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
_{o}in 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
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0; If d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this 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 design drawing of Cable Structure, asconstructed drawing and the initial measured data of Cable Structure, the Nondestructive Testing Data of support cable and 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
_{o}the 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
_{o}the 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 angledata, Cable Structure spatial data etc., measured data when reliably truly occurring close to institute's analog case.Model A
_{o}the health status evaluation object initial damage vector d of middle support cable
_{o}represent.Due to based on A
_{o}the initial value (actual measurement obtains) of the evaluation calculating all monitored amounts closely all monitored amounts, so also can be used in A
_{o}basis on, carry out Mechanics Calculation obtains, A
_{o}the evaluation of each monitored amount form monitored amount initial value vector C
_{o}.Corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent.D
_{o}a
_{o}parameter, C
_{o}by A
_{o}mechanics Calculation result composition.
3rd step: in Cable Structure military service process, actual measurement obtains the current measured value of all monitored amount of Cable Structure, composition " monitored amount current value vector C ".
4th step: at initial mechanical Calculation Basis model A
_{o}basis on carry out several times Mechanics Calculation, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u}.Concrete grammar is: at the initial mechanical Calculation Basis model A of Cable Structure
_{o}basis on carry out several times Mechanics Calculation, calculation times numerically equals the quantity of all evaluation objects, N number of evaluation object is had just to have N calculating, calculating hypothesis each time only has an evaluation object on the basis of original damage or original generalized displacement component 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 is at vectorial d
_{o}the basis that this support cable represented has a damage increases 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
_{o}there 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 xaxis 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
_{o}the 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
^{2}deng be unit change; If this load is centrepoint load, and this centrepoint load is couple, and load unit change can get 1kNm, 2kNm, 3kNm etc. for unit change; If this load is centrepoint load, and this centrepoint 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
^{3}deng be unit change), use D
_{uk}record this unit damage or load unit change, wherein k represents the numbering of the evaluation object that unit damage or unit generalized displacement or the change of generation load unit occur; 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 C, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical; The monitored amount calculation current vector C calculated each time deducts monitored amount initial value vector C
_{o}after calculate the unit damage supposed or unit generalized displacement or load unit change numerical value divided by this time again, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector; The unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively by this N number of monitored amount unit change vector; Each row of unit damage monitored numerical quantity transformation matrices correspond to a monitored amount unit change vector, and every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change amplitude of same monitored amount when different evaluation object generation unit damage or unit generalized displacement or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount.
5th step: set up linear relationship error vector e and vectorial g.Utilize (the monitored amount initial value vector C of data above
_{o}, unit damage monitored numerical quantity transformation matrices Δ C), while the 4th step calculates each time, namely calculating the increase unit damage or unit generalized displacement load unit change D of only having an evaluation object in hypothesis evaluation object each time
_{uk}, the evaluation object increasing unit damage or unit generalized displacement or load unit change in calculating each time is different from during other time calculates the evaluation object increasing unit damage or unit generalized displacement or load unit change, calculate the current value all utilizing mechanics method (such as adopting finite element method) to calculate all monitored amounts in Cable Structure each time, while calculating the monitored amount calculation current vector C of composition one each time, calculate composition injury vector d each time, originally walk out of existing injury vector d only to use in this step, the numerical value of an element is only had to get D in all elements of injury vector d
_{uk}, the numerical value of other element gets 0, the coding rule of the element of injury vector d and vectorial d
_{o}the coding rule of element identical, by C, C
^{t} _{o}, Δ C, D
_{u}, d brings formula (12) into, obtain a linear relationship error vector e, calculate a linear relationship error vector e each time, N number of evaluation object is had just to have N calculating, just there is N number of linear relationship error vector e, obtaining a vector after being added by this N number of linear relationship error vector e, is exactly final linear relationship error vector e by each element of this vector divided by the new vector obtained after N.Vector g equals final error vector e.
6th step: the hardware components of pass line structural healthy monitoring system.Hardware components at least comprises: monitored amount monitoring system (such as containing angle 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.
7th step: by monitored amount initial value vector C
_{o}, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D
_{u}parameter is kept on the hard disc of computer of operation health monitoring systems software in the mode of data file.
8th step: establishment installation and operation this method system software 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
9th step: according to monitored amount current value vector C with monitored amount initial value vector C
_{o}, unit damage monitored numerical quantity transformation matrices Δ C, evaluation object unit change vector D
_{u}and the linear approximate relationship (formula (8)) existed between evaluation object current nominal fatigue vector d (being made up of all Suo Dangqian nominal fatigue amounts), calculate the noninferior solution of evaluation object current nominal fatigue vector d according to multiobjective optimization algorithm, namely can reflect the solution of the change of the health status of evaluation object more exactly with reasonable error.
The multiobjective optimization algorithm that can adopt has a variety of, such as: the multipleobjection optimization based on genetic algorithm, the multipleobjection optimization based on artificial neural network, the multiobjective optimization algorithm based on population, the multipleobjection optimization based on ant group algorithm, leash law (ConstrainMethod), weighted method (WeightedSumMethod), Objective Programming (GoalAttainmentMethod) etc.Because various multiobjective optimization algorithm is all conventional algorithm, can realize easily, this implementation step only provides the process solving current injury vector d for Objective Programming, the specific implementation process of other algorithm can realize in a similar fashion according to the requirement of its specific algorithm.
According to Objective Programming, formula (8) can transform the multiobjective optimization question shown in an accepted way of doing sth (16) and formula (17), in formula (16), γ is a real number, R is real number field, area of space Ω limits the span (each element of the present embodiment requirements vector d is not less than 0, is not more than 1) of each element of vectorial d.Formula (16) be meant to the minimum real number γ of searching one, formula (17) is met.In formula (17), G (d) is defined by formula (18), the middle deviation allowed between G (d) and vectorial g of the product representation formula (17) of weighing vector W and γ in formula (17), the definition of g is see formula (13), and its value calculates in the 5th step.During actual computation, vector W can be identical with vectorial g.The concrete programming realization of Objective Programming has had universal program directly to adopt.Use Objective Programming just can in the hope of evaluation object current nominal fatigue vector d.
minimizeγ(16)
γ∈R,d∈Ω
G(d)Wγ≤g(17)
G(d)＝abs(ΔC·dC+C
_{o})(18)
The element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents the nominal generalized displacement of corresponding evaluation object, nominal fatigue degree or nominal load intensity of variation; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical.
Tenth step: definition evaluation object current actual damage vector d
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent the actual damage degree of corresponding evaluation object or actual generalized displacement or real load intensity of variation; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical.The evaluation object utilizing formula (15) to express current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements.
D
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion.
D
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is generalized displacement component, so a d of a bearing
^{a} _{k}represent its current actual generalized displacement numerical value.
D
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is a load, its definition is shown in formula (15), so d
^{a} _{k}represent that it is relative to setting up initial mechanical Calculation Basis model A
_{o}time the structure variable quantity of corresponding load that bears; So according to evaluation object current actual damage vector d
^{a}impaired and the degree of injury of which support cable can be defined, which load can be defined simultaneously and there occurs change and numerical value thereof.
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.
11 step: the computing machine in health monitoring systems regularly generates cable system health condition form automatically or by human users's health monitoring systems.
12 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.
13 step: get back to the 3rd step, starts by the circulation of the 3rd step to the 13 step.
Claims (1)
1. simplify angle monitor load damaged cable generalized displacement recognition methods, 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
_{1}root support cable; Determine the measured point of specifying, number to all specified points; Determined the measured straight line of each measurement point, to all measured straight line numberings of specifying; Determine the measured angle coordinate component of each measured straight line, to all measured angle coordinate component numberings; Abovementioned numbering will be used for generating vector sum matrix in subsequent step; " the whole monitored angledata of Cable Structure " is made up of abovementioned all measured angle coordinate components; For simplicity, in the method by " the monitored angledata of Cable Structure " referred to as " monitored amount "; The quantity sum of all monitored amounts is designated as M, and M should be greater than the quantity of core evaluation object, and M is less than the quantity of evaluation object; 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 centrepoint 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; Centrepoint 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 centrepoint 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 onecomponent 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 calculation obtains the measured data of initial Cable Structure, the measured data of initial Cable Structure comprises Cable Structure centrepoint 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 bearing generalized coordinate data, initial Cable Structure angledata and initial Cable Structure spatial data are 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 Nondestructive Testing Data comprising support cable, 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
_{o}the coding rule of coding rule and M monitored amount identical, the initial generalized displacement measurement data of Cable Structure bearing, support cable initial health data and Cable Structure load measurement data are utilized to set up evaluation object initial damage vector d
_{o}, vectorial d
_{o}represent with initial mechanical Calculation Basis model A
_{o}the initial health of the evaluation object of the Cable Structure represented, evaluation object initial damage vector d
_{o}element number equal N, d
_{o}element and evaluation object be onetoone relationship, vectorial d
_{o}the coding rule of element identical with the coding rule of evaluation object, if d
_{o}evaluation object corresponding to some elements be support cable, so a d in cable system
_{o}the 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 loadbearing capacity, if its numerical value is between 0 and 100%, then represent that this support cable loses the loadbearing capacity of corresponding proportion, if d
_{o}evaluation object corresponding to some elements be some generalized displacement components of some bearings, so d
_{o}the numerical value of this element represent the initial value of this generalized displacement component of this bearing, if d
_{o}evaluation object corresponding to some elements be some load, get d in this method
_{o}this 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
_{o}in each element numerical value relevant to Cable Structure generalized displacement of support get 0, if there is no the Nondestructive 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
_{o}in 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
_{o}time, the generalized displacement that Cable Structure bearing occurs relative to the bearing under Cable Structure design point, obtain the measured data of initial Cable Structure at survey calculation while, survey or consult reference materials and obtain the physical and mechanical properties parameter of the various materials that Cable Structure uses,
All Cable Structure data that the physical and mechanical properties parameter of the various materials c. used according to the measured data of the design drawing of Cable Structure, asconstructed drawing and initial Cable Structure, support cable initial health data, Cable Structure bearing initial generalized displacement measurement data, Cable Structure centrepoint 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
_{o}the Cable Structure that calculates calculates data must closely its measured data, and difference therebetween must not be greater than 5%; Corresponding to A
_{o}evaluation object health status with evaluation object initial damage vector d
_{o}represent; Corresponding to A
_{o}the initial value monitored amount initial value vector C of all monitored amount
_{o}represent; T
_{o}and d
_{o}a
_{o}parameter, by A
_{o}the initial value of all monitored amount that obtains of Mechanics Calculation result and C
_{o}the initial value of all monitored amount represented is identical, therefore alternatively C
_{o}by A
_{o}mechanics Calculation result composition;
D. from entering the circulation being walked kth step by d here;
E. at initial mechanical Calculation Basis model A
_{o}basis on carry out several times Mechanics Calculation according to step e1 to e3, by calculate obtain Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object unit change vector D
_{u};
E1. at the initial mechanical Calculation Basis model A of Cable Structure
_{o}basis 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 is at vectorial d
_{o}the basis that this support cable represented has a damage 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 is at vectorial d
_{o}the basis that this load represented has a variable quantity increases load unit change again, use D
_{uk}record 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
_{uk}evaluation object unit change vector D
_{u}an element, evaluation object unit change vector D
_{u}the coding rule of element and vectorial d
_{o}the coding rule of element identical; The evaluation object increasing unit damage or unit generalized displacement or load unit change in calculating each time is different from during other time calculates the evaluation object increasing 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, element number rule and the monitored amount initial value vector C of monitored amount calculation current vector
_{o}element number rule identical;
E2. the monitored amount calculation current vector calculated each time deducts monitored amount initial value vector C
_{o}obtain a vector, again each element of this vector is calculated the unit damage or unit generalized displacement or load unit change numerical value supposed divided by this time, obtain a monitored amount unit change vector, have N number of evaluation object just to have N number of monitored amount unit change vector;
E3. by the vectorial coding rule according to N number of evaluation object of this N number of monitored amount unit change, the Cable Structure unit damage monitored numerical quantity transformation matrices Δ C having N to arrange is formed successively; Each row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C correspond to a monitored amount unit change vector; Every a line of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C corresponds to the different unit change degree of same monitored amount when different evaluation object increases unit damage or unit generalized displacement or load unit change; The coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and vectorial d
_{o}the coding rule of element identical, the coding rule of the row of Cable Structure unit damage monitored numerical quantity transformation matrices Δ C is identical with the coding rule of M monitored amount;
F. actual measurement obtains the current measured value of all monitored amount of Cable Structure, forms monitored amount current value vector C; Monitored amount current value vector C and monitored amount initial value vector C
_{o}definition mode identical, the same monitored amount of element representation of two vectorial identical numberings is at not concrete numerical value in the same time;
G. evaluation object current nominal fatigue vector d is defined, the element number of evaluation object current nominal fatigue vector d equals the quantity of evaluation object, be onetoone relationship between the element of evaluation object current nominal fatigue vector d and evaluation object, the element numerical value of evaluation object current nominal fatigue vector d represents the nominal fatigue degree of corresponding evaluation object or nominal generalized displacement or nominal load variable quantity; The coding rule of the element of vector d and vectorial d
_{o}the coding rule of element identical;
H. according to monitored amount current value vector C with monitored amount initial value vector C
_{o}, the linear approximate relationship that exists between Cable Structure unit damage monitored numerical quantity transformation matrices Δ C and evaluation object to be asked current nominal fatigue vector d, this linear approximate relationship can be expressed as formula 1, other amount in formula 1 except d is known, solves formula 1 and just can calculate evaluation object current nominal fatigue vector d;
C=C
_{o}+ Δ Cd formula 1
I. evaluation object current actual damage vector d is defined
^{a}, evaluation object current actual damage vector d
^{a}element number equal the quantity of evaluation object, evaluation object current actual damage vector d
^{a}element and evaluation object between be onetoone relationship, evaluation object current actual damage vector d
^{a}element numerical value represent the actual damage degree of corresponding evaluation object or actual generalized displacement or real load variable quantity; Vector d
^{a}the coding rule of element and vectorial d
_{o}the coding rule of element identical;
J. the evaluation object utilizing formula 2 to express current actual damage vector d
^{a}a kth element d
^{a} _{k}with evaluation object initial damage vector d
_{o}a kth element d
_{ok}with a kth element d of evaluation object current nominal fatigue vector d
_{k}between relation, calculate evaluation object current actual damage vector d
^{a}all elements;
K=1 in formula 2,2,3 ...., N, d
^{a} _{k}represent the current actual health status of a kth evaluation object, if this evaluation object is support cable, so a d in cable system
^{a} _{k}represent its current actual damage, d
^{a} _{k}represent not damaged when being 0, when being 100%, represent that this support cable thoroughly loses loadbearing capacity, time between 0 and 100%, represent the loadbearing capacity losing corresponding proportion; If this evaluation object is generalized displacement component, so a d of a bearing
^{a} _{k}represent its current actual generalized displacement numerical value; If this evaluation object is load, so a d
^{a} _{k}represent the actual change amount of this load; So according to evaluation object current actual damage vector d
^{a}impaired and the degree of injury of which support cable can be defined, define which bearing and there occurs generalized displacement and numerical value thereof, define the numerical value which load there occurs change and change thereof; So far this method achieves damaged cable identification that reject the impact of generalized displacement of support and load change, Cable Structure, achieve reject load change and support cable health status variable effect, the identification of generalized displacement of support, achieve reject generalized displacement of support and support cable health status variable effect, the identification of load change amount; 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;
K. get back to d step, start the circulation next time being walked kth step by d.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201510438096.4A CN105115756A (en)  20150723  20150723  Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201510438096.4A CN105115756A (en)  20150723  20150723  Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change 
Publications (1)
Publication Number  Publication Date 

CN105115756A true CN105115756A (en)  20151202 
Family
ID=54663797
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201510438096.4A Pending CN105115756A (en)  20150723  20150723  Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change 
Country Status (1)
Country  Link 

CN (1)  CN105115756A (en) 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

JP2006337144A (en) *  20050601  20061214  Kawasaki Heavy Ind Ltd  Fatigue life diagnostic method and diagnostic support device of bridge 
CN103616240A (en) *  20131209  20140305  东南大学  Identifying method for damaged cable concentrated loads and generalized displacement based on angle monitoring 
CN103852334A (en) *  20140310  20140611  东南大学  Anglemonitoring load generalized displacement identification method for damaged cable 
CN103884525A (en) *  20140310  20140625  东南大学  Method for recognizing loads of damaged rod based on generalized displacement angle monitoring 

2015
 20150723 CN CN201510438096.4A patent/CN105115756A/en active Pending
Patent Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

JP2006337144A (en) *  20050601  20061214  Kawasaki Heavy Ind Ltd  Fatigue life diagnostic method and diagnostic support device of bridge 
CN103616240A (en) *  20131209  20140305  东南大学  Identifying method for damaged cable concentrated loads and generalized displacement based on angle monitoring 
CN103852334A (en) *  20140310  20140611  东南大学  Anglemonitoring load generalized displacement identification method for damaged cable 
CN103884525A (en) *  20140310  20140625  东南大学  Method for recognizing loads of damaged rod based on generalized displacement angle monitoring 
Similar Documents
Publication  Publication Date  Title 

CN101813570B (en)  Health monitoring method for recognizing damaged cable and support displacement based on mixed monitoring  
CN101806666B (en)  Health monitoring method for identifying damaged cable and support displacement based on space coordinate monitoring  
CN101813569B (en)  Health monitoring method for identifying damaged cable and support displacement based on strain monitoring  
CN101799373B (en)  Health monitor method for identifying damaged cable and support seat displacement based on angle monitor  
CN105115756A (en)  Identification method for damaged cable and generalized displacement based on angle monitoring through simplified load change  
CN105067301A (en)  Simplified strainmonitoring load damaged cable generalized displacement identifying method  
CN105067310A (en)  Method for recognizing damaged cable generalized displacement of load based on streamlined space coordinate monitoring process  
CN104990753A (en)  Simplified method of recognizing load damage cable and generalized displacement through mixed monitoring  
CN101793628B (en)  Cable structure health monitoring method based on hybrid monitoring  
CN105067309A (en)  Simplified cableforcemonitoring load damaged cable generalized displacement identifying method  
CN105115748A (en)  Identification method for damaged cable based on angle monitoring through simplified generalized displacement and load change  
CN104990723A (en)  Simplified load damaged cable identification method through hybrid monitoring under generalized displacement  
CN105115752A (en)  Method for identifying damaged cable based on space coordinate monitoring through simplified generalized displacement and load change  
CN105067349A (en)  Method for recognizing damaged cable of load based on cable force monitoring process of streamlined and generalized displacement  
CN104990733A (en)  Simplified load problematical cable generalized displacement identification method through space coordinate monitoring  
CN105115771A (en)  Progressive identification method for load, damaged cable and generalized displacement through simplified space coordinate monitoring  
CN104990577A (en)  Simplified load damaged cable generalized displacement progressiveidentification method through hybrid monitoring  
CN105136501A (en)  Simplified angle monitoring load faulty cable generalized displacement identification method  
CN105067317A (en)  Method for recognizing damaged cable generalized displacement of load based on streamlined hybrid monitoring process  
CN104990724A (en)  Simplified load damaged cable generalized displacement progressiveidentification method through strain monitoring  
CN104990587A (en)  Simplified method of recognizing problematic load cable generalized displacement through strain monitoring  
CN105004557A (en)  Method for recognizing damaged cable of load based on streamlined angle monitoring process of generalized displacement  
CN105067303A (en)  Simplified generalized displacement anglemonitoring load damaged cable progressive identifying method  
CN105004563A (en)  Method for recognizing damaged cable of load based on streamlined hybrid monitoring process of generalized displacement  
CN105115767A (en)  Progressive identification method for load, damaged cable and generalized displacement through simplified cable force monitoring 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
RJ01  Rejection of invention patent application after publication  
RJ01  Rejection of invention patent application after publication 
Application publication date: 20151202 