CN108828076A - A kind of bridge quality analysis system and method based on image analysis - Google Patents
A kind of bridge quality analysis system and method based on image analysis Download PDFInfo
- Publication number
- CN108828076A CN108828076A CN201810597256.3A CN201810597256A CN108828076A CN 108828076 A CN108828076 A CN 108828076A CN 201810597256 A CN201810597256 A CN 201810597256A CN 108828076 A CN108828076 A CN 108828076A
- Authority
- CN
- China
- Prior art keywords
- equation
- point
- value
- bridge
- section
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4472—Mathematical theories or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
Abstract
The invention belongs to bridge approach road technical fields, disclose a kind of bridge quality analysis system and method based on image analysis, bridge quality analysis system based on image analysis is provided with alarm and acoustic detector, the alarm is provided with shell, LED light, jingle bell, and the detector is provided with operator control panel, handle, stock, detecting head.The jingle bell is embedding in the shell, the LED light around jingle bell and it is embedding the upper end of the stock connects handle by nut on the shell, one end of the handle is fixedly connected operator control panel by nut, and the bottom of the stock is welded with detecting head.The configuration of the present invention is simple, easy to operate, acoustic detector can accurately detect the situation inside bridge, when there are quality problems, alarm can timely issue prompting, give operator's information, the safety for ensuring bridge greatly reduces the generation of accident.
Description
Technical field
The invention belongs to bridge approach road technical field more particularly to a kind of bridge quality analysis systems based on image analysis
And method.
Background technique
Currently, the prior art commonly used in the trade is such：
Currently, highway communication is national economy main artery, and road and bridge engineering structure is then the throat thoroughfare of highway communication,
Ensure that highway plays a crucial role in unobstructed.But after road and bridge are open to the traffic, over time, Zhong Zhongyin
Element can be such that road and bridge degree of safety is declined, and such as former design is not up to requirement, and construction meet design requirements, and road and bridge have disease
Evil, material aging, corrosion and conserve not in time, vehicular load increase or the volume of traffic increase severely, Highway-bridge Expansion Joints damage or bridge floor not
It is smooth to road and bridge structure bring adverse effect etc..At present, only increase from vehicular load and the volume of traffic increases severely this factor
From the point of view of, many road and bridge there is huge security risk.
In conclusion problem of the existing technology is：
Traditional detection method is artificial vision's recognition methods, and the searching in usual crack and form are passed through aobvious by investigator
Micro mirror measurement record, close traffic, detection work operations amount is big, detection time is long, testing cost there is needing for this method
The shortcomings that high, testing result human factor is affected, while testing staff need to go deep into crack institute by detection vehicle or scaffold
In region detection, this is also that detection brings many unpredictable security risks.
The Dispersion Equation of wave is generally a binary transcendental equation about wave number and frequency, when solve complex wave number field in more
When the solution for the relationship of dissipating, equation becomes more complicated ternary transcendental equation, and the coefficient of Dispersion Equation may contain plural number,
Therefore the solution of such issues that is highly difficult, generally can only solve disperse relationship to the case where extremely special very simple, this is right
In the analysis of the sonic sensor of various different structures be far from being enough.
The safety of bridge is evaluated in the prior art, depends on human factor evaluation, deviation is big.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of bridge quality analysis system based on image analysis
System.
The invention is realized in this way a kind of bridge mass analysis method based on image analysis, including：
Using the surface for the detecting head ultrasonic listening bridge that acoustic detector integrates, in detection, detecting head is internally integrated
Analysis module to obtain the deformation data on bridge surface, crack data modulus value restrains and solves deformation data, crack data
Transcendental equation；For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 etc.
Valence；In modulus, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, at section [d, e]
Inside there is zero point s；With modulus value restrain solve f (x)=0 equation when, first find out | f (x) | minimum point；Either segment in x-axis
10 parts of equal part of section [m, n] is had 11 Along ents, respectively by minizone [m, n]Compare on this 11 nodes | f (x) | value, find out
Along ent x locating for minimum value0；If not having in [m, n] | f (x) | minimum point, then x0=m or x0=n；In [if m, n]
Have | f (x) | minimum point, then x0It takesIt is therein some
Value；Work as x0≠ m and x0≠ n has in [m, n] | f (x) | minimum point；With with x0Adjacent Along ent is endpoint, takes sectionBy 10 parts of its equal part, Along ent x locating for minimum value is equally found out1, and | f (x1)|≤|
f(x0)|；Above procedure is constantly repeated, a series of point x is obtained0,x1,x2,...；These points are minimum point in the section [m, n]
Numerical solution under different accuracy；
The data detected are sent to the computer being attached thereto, computer carries out servicing bridges Failure Factors and comments
It is fixed：Failure Assessment Curves equation is：The equation of vertical line is：Value depend on Reinforced Concrete Materials characteristic：It is calculated as followsValue：Computer also carries out the fatigue assessment of planar disfigurement, is first depending on fatigue crack
Spreading rate da/dN and crack tip stress intensity factor amplitude of variation Δ K relational expressionDetermination is providing
Cycle period in fatigue crack propagation and final size；Then according to the criterion provided, judge the planar disfigurement
Whether can deformation occurs and fatigue fracture；
Computer in the form of images shows data that, when dangerous parameter occurs, alarm device is alarmed.
Further, the analysis module that detecting head is internally integrated is to deformation data, the crack data mould for obtaining bridge surface
Value convergence solves the transcendental equation of deformation data, crack data, specifically includes：
1) it is restrained with modulus value and solves transcendental equation
For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 etc.
Valence；
In modulus, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, in section
There is zero point s in [d, e]；
With modulus value restrain solve the equation when, first find out | f (x) | minimum point；Either segment minizone in x-axis
10 parts of equal part of section [m, n] is had 11 Along ents, respectively by [m, n]Compare this
On 11 nodes | f (x) | value, find out Along ent x locating for minimum value0；If not having in [m, n] | f (x) | minimum point,
So x0=m or x0=n；If having in [m, n] | f (x) | minimum point, then x0It takes
Some value therein；Work as x0≠ m and x0≠ n has in [m, n] | f (x) | minimum point；With with x0Adjacent Along ent is end
Point, takes sectionBy 10 parts of its equal part, Along ent x locating for minimum value is equally found out1, and |
f(x1)|≤|f(x0)|；Above procedure is constantly repeated, a series of point x is obtained0,x1,x2,...；These points are in the section [m, n]
Numerical solution of the minimum point under different accuracy；
2) after finding out minimum point, zero point and non-zero points then are distinguished from these minimum points；
When minimum point is non-zero points, section [a, b] equal part is obtained into a series of point x0,x1,x2..., and|f(x0)|≥|f(x1)|≥|f(x2)|≥...≥|f(c)|,| f (c) | it is big
In zero constant,For a limited big constantWherein
When minimum point is zero point, section [d, e] equal part is obtained into a series of point x0,x1,x2..., and|f(x0)|≥|f(x1)|≥|f(x2)|≥...≥0,
3) analysis for comparing zero point and non-zero points, is obtaining minimum point x0,x1,x2... after, according toValue
Judge whether minimum is zero point, chooses a particular value M, after restraining n step, ifThen extreme point is
Zero point；M and convergent step number are depending on different problems；In the case where 10 parts of equal part,It is zero point；
Modulus value | f (xn) | convergence rate and initial section [m, n] equal part quantity of size and [m, n] for being taken it is related, original area
Between it is smaller, equal part is more, restrain it is faster；
Further, the analysis module that detecting head is internally integrated is to deformation data, the crack data mould for obtaining bridge surface
Value convergence solves the transcendental equation of deformation data, crack data, further comprises：
The smeared cracking solved in real wave number and empty wave number plane is restrained with modulus value：
The Dispersion Equation that wave is propagated in different structure has different concrete forms, and form is about frequencies omega and wave number ξ
Transcendental equation：
G (ω, ξ)=0
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is real number and pure imaginary number；Since Dispersion Equation is in frequency
Solution in the plane of rate ω and wave number ξ is curve, with line element to entire flat scanning；
Any of first fixed frequency ω and wave number ξ, if fixed wave number is ξ0, equation g (ω, ξ)=0 becomes at this time：
F (ω)=g (ω, ξ0)=0
With the line infinitesimal of a length of 3t to straight line ξ=ξ0It is scanned, scan start point ω0；In section [ω0,ω0+3t]
It is interior, comparison node ω0,ω0+t,ω0+2t,ω0+ 3t locates modulus value | f (ω0)|,|f(ω0+t)|,|f(ω0+2t)|,|f(ω0+
3t) | size, find out minimum value；In ω0+ t or ω0+ 2t obtains minimum value, then further equal part, and judges whether it is zero point；
If minimum value takes on endpoint, enters next section and scan；To prevent minimum point to be by chance in endpoint ω0The place+3t, under
One section is taken as [ω0+2t,ω0+ 5t], this section includes ω0+3t；Adjacent section is separated by 2t；If dividing n parts for section etc.,
Then adjacent interval is separated by n-1 parts；When to straight line ξ=ξ0The end of scan, into lower straight line ξ=ξ0+ Δ ξ scanning, multiple scanning
Repeatedly, the non trivial solution in entire ω, ξ plane is obtained；In order to enable the solution curve arrived is more complete, scans and tie in fixed wave number
Shu Hou fixed frequency can be scanned again.
The smeared cracking solved in complex wave number field is restrained with modulus value：
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is plural number；ξ=a+bi being enabled, wherein a, b is real number,
Then equation g (ω, ξ)=0 is turned to：
H (ω, a, b)=g (ω, a+bi)
Since solution of the Dispersion Equation in the space that real frequency ω and complex wave number ξ is formed is curve, with face element to entire sky
Between scan；Any of first fixed frequency ω, wave number real part a, wave number imaginary part b, if fixed frequency is ω0, equation h at this time
(ω, a, b)=g (ω, a+bi) becomes：
Q (a, b)=h (ω0,a,b)
In plane ω=ω0It is interior, plane is scanned with face infinitesimal [a, a+3t] × [b, b+3s]；Scan start point is
(a0,b0)：
In face infinitesimal [a0,a0+3t]×[b0,b0+ 3s] in, compare 16 node (a0,b0), (a0+t,b0) ..., (a0+
3t,b0+ 3s) at | q (a, b) | size, find out minimum value；If at the internal node of face element, i.e. (a0+t,b0+ s), (a0+
2t,b0+ s), (a0+t,b0+ 2s), (a0+2t,b0+ 2s), minimum value is obtained, then further equal part, and judge whether it is zero point；
If minimum value takes on the boundary node of face element, enters next section and scan；To prevent minimum point by chance
It is in face element boundary, next section is taken as [a0+2t,a0+5t]×[b0,b0+ 3s], when along the scanning of a axis, adjacent face element
It is separated by 2t；If face element divides n parts along a axis etc., adjacent face element is separated by n-1 parts；It scans, obtains after a axis end of scan, then along b axis
To entire plane ω=ω0Interior equation q (a, b)=h (ω0, a, b) solution；
To plane ω=ω0The end of scan, into next plane ω=ω0+ Δ ω scanning, multiple scanning are obtained whole
A ω, a, equation q (a, b)=h (ω in the space b0, a, b) solution.
Further, the analysis module that detecting head is internally integrated is to deformation data, the crack data mould for obtaining bridge surface
Value convergence solves the transcendental equation of deformation data, crack data, further comprises：
The solution of any first transcendental equation is solved with modulus value convergence：
Assuming that n member transcendental equation is：
f(x1,x2,...,xn)=0
To have scanning infinitesimal and solution curve and only one intersection point, if solution curve is that m ties up (m≤n), infinitesimal choosing is scanned
For n-m dimension, fixed x1,x2,...,xnIn m amount；If fixed x1,x2,...,xm, with [xm+1,xm+1+Δxm+1]×[xm+2,
xm+2+Δxm+2]×...×[xn,xn+Δxn] as scanning infinitesimal scanning；When solution curve is that m ties up (m≤n), if scanning infinitesimal dimension
Degree, which is greater than n-m, to be tieed up, then the intersection point of each infinitesimal and solution curve is not unique, and many solutions are missed；If scanning infinitesimal dimension is less than n-m
Dimension, then each infinitesimal hardly intersects with solution curve；In three dimensions, solution curve is one-dimensional curve, then uses plane infinitesimal
Scanning；If solution curve is Two-dimensional Surfaces, equation q (a, b)=h (ω is fixed0, a, b) in two amount, with line infinitesimal scan；
If solution curve is discrete point, scanned with body infinitesimal；If solution curve is three-dimensional space, each point in direct clarifying space
Equation modulus value whether converge to zero；When practical solution n member transcendental equation, if the value of uncertain solution curve dimension m, will scan micro-
The dimension of member is tested from high to low, and obtained solution reduces and fine and close suddenly disappearing up to solution with scanning infinitesimal dimension, is scanned
Dimension is n-m-1 dimension, determines the value of m.
Further, in computer progress servicing bridges evaluation, further comprise：
Establish set of factors：
The various parameters compositing factor set for influencing servicing bridges, takes the second level factor U={ u1,u2,u3}={ defect, material
Matter, internal environment }, u1={ u11,u12The fracture failure of }={ planar disfigurement, the fatigue failure of planar disfigurement }, u2={ u21,
u22,u23}={ material processing quality, the mechanical property of materials, design rationality }, u3={ u31,u32,u33}={ pressure, temperature are rotten
Corrosion }；
Establish evaluate collection：
In order to quantitative analysis is carried out to each evaluation index it needs to be determined that each index evaluate collection, evaluated using 5 grades of hundred-mark system
Evaluate collection V is divided 5 opinion ratings, i.e. V={ v1, v2, v3, v4, v5}={ is minimum, very little, small, larger, big }, wherein v1For
Servicing bridges multifactor failure risk is minimum, and scoring section is 90~100, intermediate value 95；v2For risk very little, area of scoring
Between be 80~89, intermediate value 84.5；The rest may be inferred；Select the intermediate value in each section as the parameter of grade, then corresponding to 5 grades
Parameter be { 95,84.5,74.5,64.5,49.5 }, parameter column vector be ν={ 95,84.5,74.5,64.5,49.5 }T；
Establish weight sets：
(1) recursive hierarchy structure is established：
According to servicing bridges factor of evaluation collection, that is, assessment indicator system of foundation, each factor for being included by problem is divided into four
A level：First layer is the general objective layer G of evaluation, i.e. servicing bridges comprehensive safety；The second layer is rule layer C, i.e. defect, important
Component, material and internal environment；Finally using a specific targets as third layer, i.e. indicator layer P；
(2) multilevel iudge matrix two-by-two is constructed：
In terms of inviting bridge security, importance degree tax is successively carried out to each element according to 1~9 scaling law between any two
Value, Judgement Matricies U=(uij)n×n, wherein uijExpression factor uiAnd ujRelative to the importance value of rule layer, having property of matrix U
Matter：uii=1, uij=1/uji, i, j=1,2 ..., n obtain judgment matrix：By matrix X1~X5By row normalization, i.e.,：
Calculating matrix Y is：
(3) under single criterion element relative weighting calculating：
Y matrix by rows is added, by formulaIt obtains：
W1=(2.652 0.686 0.253 0.409)T
W2=(1 1)T
W3=(1.273 0.371 0.221 2.135)T
W4=(1.9 0.319 0.781)T
W5=(2.121 0.604 0.275)T
Obtain and vector is normalized, by formulaWeight vector can be obtained：
(4) consistency check of judgment matrix：
Calculate the Maximum characteristic root λ of judgment matrixmax, by formulaIt is calculated:
According to formulaConsistency check is carried out, is obtained：
CI1=0.019
CI2=0
CI3=0.031
CI4=0.020
CI5=0.048
By formula?：
CR1=0.022
CR2=0
CR3=0.035
CR4=0.038
CR5=0.092
CR < 0.1 is all satisfied coherence request, therefore the relative weighting of each factor
Degree of membership calculates：
Multidigit uses frequency number analysis, the danger to the indices being evaluated by evaluate collection to servicing bridges indices
Degree is graded, and the degree of membership of set of factors is obtained：
It determines and judges Subject Matrix：
By obtaining the relative defects matrix of k-th of set of factors：
Wherein：
In formula：RkThe relative defects matrix of-k-th set of factors；
rkijI-th of factor of-k-th set of factors belongs to the degree of membership of the j in evaluate collection；
pkij- group membership is rated the frequency of j to i-th of factor index of k-th of set of factors；
Construct fuzzy matrix for assessment：
Fuzzy matrix for assessment B can be constructed by the weight vector W and matrix R of each index,
Calculate Comprehensive Evaluation result：
By the parameter column vector of fuzzy matrix for assessment B and evaluate collection, Comprehensive Evaluation result Z can be acquired；
Z=BV
Arrive as available from the above equation fuzzy overall evaluation further according to opinion rating as a result, provides, it is more can to evaluate servicing bridges
Factor failure risk size.
Another object of the present invention is to provide the bridge mass analysis methods described in a kind of realize based on image analysis
Computer program.
Another object of the present invention is to provide the bridge mass analysis methods described in a kind of realize based on image analysis
Information data processing terminal.
Another object of the present invention is to provide a kind of computer readable storage mediums, including instruction, when it is in computer
When upper operation, so that computer executes the bridge mass analysis method based on image analysis.
Another object of the present invention is to provide a kind of bridge quality analysis system based on image analysis, is provided with：Report
Alert device and acoustic detector；
Alarm is provided with shell, LED light, jingle bell；Jingle bell is embedding in the shell, and LED light is around jingle bell and is embedded in shell
On；
Detector is provided with operator control panel, handle, stock, detecting head；
The upper end of stock connects handle by nut；One end of handle is fixedly connected operator control panel by nut, stock
Bottom is welded with detecting head.
Further, the external computer of acoustic detector, for the data that detected to be shown in the form of images；
Alarm connects power supply and is built-in with reserve battery, for working on when accidental power failure.
Advantages of the present invention and good effect are：
The configuration of the present invention is simple, easy to operate, acoustic detector can accurately detect the situation inside bridge, when out
When existing quality problems, alarm can timely issue prompting, give operator's information, it is ensured that the safety of bridge, significantly
Reduce the generation of accident.
Detecting head of the invention, after obtaining some node of equation modulus value minimalization in scanning infinitesimal, with this node
Centered on, adjacent node is boundary node, forms new infinitesimal, and suitable step-length is taken to divide this infinitesimal, calculates new infinitesimal node
On equation modulus value, compare the node for obtaining and being minimized.It repeats the above process, a series of available modulus value are successively decreased minimum
Value node, if the modulus value of initial minimum value node is intended to infinite, this minimum section than the modulus value of upper newest minimum value node
Point is zero point, this shows in this detecting head acoustic detection head, and wave can be propagated according to the wave number of the point with frequency.If
It is intended to a limited big constant, then this minimum value node is not zero a little, this shows that in this acoustic detection head, wave can not
It can be propagated according to the wave number of the point with frequency.Can use can when convergent step number controls this acoustic detection head medium wave propagation
The wave number of energy and the solving precision of frequency.The Dispersion Equation of wave is generally a binary transcendental equation about wave number and frequency,
When solving the solution of disperse relationship in complex wave number field, equation becomes more complicated ternary transcendental equation, and Dispersion Equation is
Number may containing plural number, therefore such issues that solution it is highly difficult, generally can only be to extremely special very simple the case where
Disperse relationship is solved, this is far from being enough for the analysis of the acoustic detection head of various different structures；And utilize the present invention
The method of offer efficiently, can be solved widely in the structures such as various table sound waves or bulk acoustic wave resonator, filter and sensor
The dispersion equation and frequency characteristic of wave propagation problem.After solution obtains disperse relationship, it can be easy to solve bridge corresponding
Displacement field or other physical change values.
Bridge quantifying defects comprehensive safety evaluation method provided by the invention is come using the calculation formula with degree of precision
The quantization for realizing defect, to keep the result detected more acurrate, true, reliable；And it calculates and is damaged caused by servicing bridges；
Obtain maximum pressure, crack length, crack depth, residual compression, the remanent fatigue life that can be born that servicing bridges can be born
With the technical indicators such as corresponding safe score, the difficulty for being unable to dynamic detection disaster trend is overcome, it can be more preferable, more accurately in time
It was found that disaster, accomplishes to prevent in advance；Using comprehensive evaluation system, quantification is combined with qualitative analysis, and incorporation engineering is real
Evaluate collection is established on border, establishes overall merit judgment matrix, according to each Failure Factors to the damage ratio of bridge and its to bridge peace
The total damage ratio of the weight calculation that umbra is rung is abandoned by the way of the evaluation of single angle, excessively dependence or field data, comprehensive
Consider to influence all principal elements of servicing bridges safety, and clearly each influence connects each other, makes on this basis comprehensive
The safety evaluation of conjunction property；Can not only correctly obtain whether can trouble free service conclusion, moreover it is possible to solve the problems, such as safe coefficient；Simplification is commented
Valence process, the subjectivity for eliminating evaluation is random, and it is practical to be applied to engineering convenient for common engineers and technicians.Of the invention is reliable
Property it is high, operability is good, reflect engineering reality with enabling assessment result more objective reality.
Detailed description of the invention
Fig. 1 is that the structure of alarm device in the bridge quality analysis system provided in an embodiment of the present invention based on image analysis is shown
It is intended to；
Fig. 2 is that the structure of detector in the bridge quality analysis system provided in an embodiment of the present invention based on image analysis is shown
It is intended to；
In figure：1, shell；2, LED light；3, jingle bell；4, operator control panel；5, handle；6, stock；7, detecting head.
Fig. 3 is that the phenogram of planar disfigurement provided in an embodiment of the present invention illustrates intention；
Fig. 4 is stress range distribution straight line schematic diagram provided in an embodiment of the present invention.
Specific embodiment
In order to further understand the content, features and effects of the present invention, the following examples are hereby given, and cooperate attached drawing
Detailed description are as follows.
As depicted in figs. 1 and 2, the bridge quality analysis system provided in an embodiment of the present invention based on image analysis is provided with
Alarm and acoustic detector, the alarm are provided with shell 1, LED light 2, jingle bell 3, and the detector is provided with operator control panel
4, handle 5, stock 6, detecting head 7.For the jingle bell 3 in shell 1, the LED light 2 is around jingle bell 3 and is embedded in shell 1
On, the upper end of the stock 6 connects handle 5 by nut, and one end of the handle 5 is fixedly connected operator control panel 4 by nut,
The bottom of the stock 6 is welded with detecting head 7.
The external computer of acoustic detector, can show in the form of images the data that detected,
The alarm connects power supply and is built-in with reserve battery, can work on when accidental power failure.
In the use of the present invention,
Alarm device external power supply, when using detector, by adjusting suitable function to the operation in operator control panel 4,
Then hold handle 5, with the surface of 7 proximity detection bridge of detecting head, detector the data detected can be sent to
Connection computer, computer can show in the form of images data, when dangerous parameter occurs, alarm device
It can react, LED light 2 flashes, and jingle bell 3 sounds an alarm, and reminds operator, and bridge quality is problematic herein, needs in time
It repairs, avoids that accident occurs, and reserve battery can be in the case where having a power failure, so that alarm device still works.
Below with reference to concrete analysis, the invention will be further described.
Bridge mass analysis method provided in an embodiment of the present invention based on image analysis, including：
Using the surface for the detecting head ultrasonic listening bridge that acoustic detector integrates, in detection, detecting head is internally integrated
Analysis module to obtain the deformation data on bridge surface, crack data modulus value restrains and solves deformation data, crack data
Transcendental equation；For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 etc.
Valence；In modulus, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, at section [d, e]
Inside there is zero point s；With modulus value restrain solve f (x)=0 equation when, first find out | f (x) | minimum point；Either segment in x-axis
10 parts of equal part of section [m, n] is had 11 Along ents, respectively by minizone [m, n]Than
Compared on this 11 nodes | f (x) | value, find out Along ent x locating for minimum value0；If not having in [m, n] | f (x) | minimum
Point, then x0=m or x0=n；If having in [m, n] | f (x) | minimum point, then x0It takesSome value therein；Work as x0≠ m and x0≠ n has in [m, n] | f
(x) | minimum point；With with x0Adjacent Along ent is endpoint, takes sectionBy its equal part
10 parts, equally find out Along ent x locating for minimum value1, and | f (x1)|≤|f(x0)|；Above procedure is constantly repeated, a system is obtained
The point x of column0,x1,x2,...；These points are numerical solution of the minimum point under different accuracy in the section [m, n]；
The data detected are sent to the computer being attached thereto, computer carries out servicing bridges Failure Factors and comments
It is fixed：Failure Assessment Curves equation is：The equation of vertical line is： Value depend on Reinforced Concrete Materials characteristic：It is calculated as followsValue：Computer also carries out the fatigue assessment of planar disfigurement, is first depending on fatigue crack
Spreading rate da/dN and crack tip stress intensity factor amplitude of variation Δ K relational expressionDetermination is providing
Cycle period in fatigue crack propagation and final size；Then according to the criterion provided, judge the planar disfigurement
Whether can deformation occurs and fatigue fracture；
Computer in the form of images shows data that, when dangerous parameter occurs, alarm device is alarmed.
The analysis module that detecting head is internally integrated asks the deformation data on acquisition bridge surface, the convergence of crack data modulus value
The transcendental equation for solving deformation data, crack data, specifically includes：
1) it is restrained with modulus value and solves transcendental equation
For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 etc.
Valence；
In modulus, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, in section
There is zero point s in [d, e]；
With modulus value restrain solve the equation when, first find out | f (x) | minimum point；Either segment minizone in x-axis
10 parts of equal part of section [m, n] is had 11 Along ents, respectively by [m, n]
Compare on this 11 nodes | f (x) | value, find out Along ent x locating for minimum value0；If not having in [m, n] | f (x) | it is minimum
It is worth point, then x0=m or x0=n；If having in [m, n] | f (x) | minimum point, then x0It takesSome value therein；Work as x0≠ m and x0≠ n has in [m, n] | f
(x) | minimum point；With with x0Adjacent Along ent is endpoint, takes sectionBy its equal part
10 parts, equally find out Along ent x locating for minimum value1, and | f (x1)|≤|f(x0)|；Above procedure is constantly repeated, a system is obtained
The point x of column0,x1,x2,...；These points are numerical solution of the minimum point under different accuracy in the section [m, n]；
2) after finding out minimum point, zero point and non-zero points then are distinguished from these minimum points；
When minimum point is non-zero points, section [a, b] equal part is obtained into a series of point x0,x1,x2..., and|f(x0)|≥|f(x1)|≥|f(x2)|≥...≥|f(c)|,| f (c) | it is big
In zero constant,For a limited big constantWherein
When minimum point is zero point, section [d, e] equal part is obtained into a series of point x0,x1,x2..., and|f(x0)|≥|f(x1)|≥|f(x2)|≥...≥0,
3) analysis for comparing zero point and non-zero points, is obtaining minimum point x0,x1,x2... after, according toValue
Judge whether minimum is zero point, chooses a particular value M, after restraining n step, ifThen extreme point is
Zero point；M and convergent step number are depending on different problems；In the case where 10 parts of equal part,It is zero point；
Modulus value | f (xn) | convergence rate and initial section [m, n] equal part quantity of size and [m, n] for being taken it is related, original area
Between it is smaller, equal part is more, restrain it is faster；
The analysis module that detecting head is internally integrated asks the deformation data on acquisition bridge surface, the convergence of crack data modulus value
Solve deformation data, crack data transcendental equation, further comprise：
The smeared cracking solved in real wave number and empty wave number plane is restrained with modulus value：
The Dispersion Equation that wave is propagated in different structure has different concrete forms, and form is about frequencies omega and wave number ξ
Transcendental equation：
G (ω, ξ)=0
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is real number and pure imaginary number；Since Dispersion Equation is in frequency
Solution in the plane of rate ω and wave number ξ is curve, with line element to entire flat scanning；
Any of first fixed frequency ω and wave number ξ, if fixed wave number is ξ0, equation g (ω, ξ)=0 becomes at this time：
F (ω)=g (ω, ξ0)=0
With the line infinitesimal of a length of 3t to straight line ξ=ξ0It is scanned, scan start point ω0；In section [ω0,ω0+3t]
It is interior, comparison node ω0,ω0+t,ω0+2t,ω0+ 3t locates modulus value | f (ω0)|,|f(ω0+t)|,|f(ω0+2t)|,|f(ω0+
3t) | size, find out minimum value；In ω0+ t or ω0+ 2t obtains minimum value, then further equal part, and judges whether it is zero point；
If minimum value takes on endpoint, enters next section and scan；To prevent minimum point to be by chance in endpoint ω0The place+3t, under
One section is taken as [ω0+2t,ω0+ 5t], this section includes ω0+3t；Adjacent section is separated by 2t；If dividing n parts for section etc.,
Then adjacent interval is separated by n-1 parts；When to straight line ξ=ξ0The end of scan, into lower straight line ξ=ξ0+ Δ ξ scanning, multiple scanning
Repeatedly, the non trivial solution in entire ω, ξ plane is obtained；In order to enable the solution curve arrived is more complete, scans and tie in fixed wave number
Shu Hou fixed frequency can be scanned again.
The smeared cracking solved in complex wave number field is restrained with modulus value：
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is plural number；ξ=a+bi being enabled, wherein a, b is real number,
Then equation g (ω, ξ)=0 is turned to：
H (ω, a, b)=g (ω, a+bi)
Since solution of the Dispersion Equation in the space that real frequency ω and complex wave number ξ is formed is curve, with face element to entire sky
Between scan；Any of first fixed frequency ω, wave number real part a, wave number imaginary part b, if fixed frequency is ω0, equation h at this time
(ω, a, b)=g (ω, a+bi) becomes：
Q (a, b)=h (ω0,a,b)
In plane ω=ω0It is interior, plane is scanned with face infinitesimal [a, a+3t] × [b, b+3s]；Scan start point is
(a0,b0)：
In face infinitesimal [a0,a0+3t]×[b0,b0+ 3s] in, compare 16 node (a0,b0), (a0+t,b0) ..., (a0+
3t,b0+ 3s) at | q (a, b) | size, find out minimum value；If at the internal node of face element, i.e. (a0+t,b0+ s), (a0+
2t,b0+ s), (a0+t,b0+ 2s), (a0+2t,b0+ 2s), minimum value is obtained, then further equal part, and judge whether it is zero point；
If minimum value takes on the boundary node of face element, enters next section and scan；To prevent minimum point by chance
It is in face element boundary, next section is taken as [a0+2t,a0+5t]×[b0,b0+ 3s], when along the scanning of a axis, adjacent face element
It is separated by 2t；If face element divides n parts along a axis etc., adjacent face element is separated by n-1 parts；It scans, obtains after a axis end of scan, then along b axis
To entire plane ω=ω0Interior equation q (a, b)=h (ω0, a, b) solution；
To plane ω=ω0The end of scan, into next plane ω=ω0+ Δ ω scanning, multiple scanning are obtained whole
A ω, a, equation q (a, b)=h (ω in the space b0, a, b) solution.
The analysis module that detecting head is internally integrated asks the deformation data on acquisition bridge surface, the convergence of crack data modulus value
Solve deformation data, crack data transcendental equation, further comprise：
The solution of any first transcendental equation is solved with modulus value convergence：
Assuming that n member transcendental equation is：
f(x1,x2,...,xn)=0
To have scanning infinitesimal and solution curve and only one intersection point, if solution curve is that m ties up (m≤n), infinitesimal choosing is scanned
For n-m dimension, fixed x1,x2,...,xnIn m amount；If fixed x1,x2,...,xm, with [xm+1,xm+1+Δxm+1]×[xm+2,
xm+2+Δxm+2]×...×[xn,xn+Δxn] as scanning infinitesimal scanning；When solution curve is that m ties up (m≤n), if scanning infinitesimal dimension
Degree, which is greater than n-m, to be tieed up, then the intersection point of each infinitesimal and solution curve is not unique, and many solutions are missed；If scanning infinitesimal dimension is less than n-m
Dimension, then each infinitesimal hardly intersects with solution curve；In three dimensions, solution curve is one-dimensional curve, then uses plane infinitesimal
Scanning；If solution curve is Two-dimensional Surfaces, equation q (a, b)=h (ω is fixed0, a, b) in two amount, with line infinitesimal scan；
If solution curve is discrete point, scanned with body infinitesimal；If solution curve is three-dimensional space, each point in direct clarifying space
Equation modulus value whether converge to zero；When practical solution n member transcendental equation, if the value of uncertain solution curve dimension m, will scan micro-
The dimension of member is tested from high to low, and obtained solution reduces and fine and close suddenly disappearing up to solution with scanning infinitesimal dimension, is scanned
Dimension is n-m-1 dimension, determines the value of m.
Computer carries out：
Establish set of factors：
The various parameters compositing factor set for influencing servicing bridges, takes the second level factor U={ u1,u2,u3}={ defect, material
Matter, internal environment }, u1={ u11,u12The fracture failure of }={ planar disfigurement, the fatigue failure of planar disfigurement }, u2={ u21,
u22,u23}={ material processing quality, the mechanical property of materials, design rationality }, u3={ u31,u32,u33}={ pressure, temperature are rotten
Corrosion }；
Establish evaluate collection：
In order to quantitative analysis is carried out to each evaluation index it needs to be determined that each index evaluate collection, evaluated using 5 grades of hundred-mark system
Evaluate collection V is divided 5 opinion ratings, i.e. V={ v1, v2, v3, v4, v5}={ is minimum, very little, small, larger, big }, wherein v1For
Servicing bridges multifactor failure risk is minimum, and scoring section is 90~100, intermediate value 95；v2For risk very little, area of scoring
Between be 80~89, intermediate value 84.5；The rest may be inferred；Select the intermediate value in each section as the parameter of grade, then corresponding to 5 grades
Parameter be { 95,84.5,74.5,64.5,49.5 }, parameter column vector be ν={ 95,84.5,74.5,64.5,49.5 }T；
Establish weight sets：
(1) recursive hierarchy structure is established：
According to servicing bridges factor of evaluation collection, that is, assessment indicator system of foundation, each factor for being included by problem is divided into four
A level：First layer is the general objective layer G of evaluation, i.e. servicing bridges comprehensive safety；The second layer is rule layer C, i.e. defect, important
Component, material and internal environment；Finally using a specific targets as third layer, i.e. indicator layer P；
(2) multilevel iudge matrix two-by-two is constructed：
In terms of inviting bridge security, importance degree tax is successively carried out to each element according to 1~9 scaling law between any two
Value, Judgement Matricies U=(uij)n×n, wherein uijExpression factor uiAnd ujRelative to the importance value of rule layer, having property of matrix U
Matter：uii=1, uij=1/uji, i, j=1,2 ..., n obtain judgment matrix：By matrix X1~X5By row normalization, i.e.,：
Calculating matrix Y is：
(3) under single criterion element relative weighting calculating：
Y matrix by rows is added, by formulaIt obtains：
W1=(2.652 0.686 0.253 0.409)T
W2=(1 1)T
W3=(1.273 0.371 0.221 2.135)T
W4=(1.9 0.319 0.781)T
W5=(2.121 0.604 0.275)T
Obtain and vector is normalized, by formulaWeight vector can be obtained：
(4) consistency check of judgment matrix：
Calculate the Maximum characteristic root λ of judgment matrixmax, by formulaIt is calculated:
According to formulaConsistency check is carried out, is obtained：
CI1=0.019
CI2=0
CI3=0.031
CI4=0.020
CI5=0.048
By formula?：
CR1=0.022
CR2=0
CR3=0.035
CR4=0.038
CR5=0.092
CR < 0.1 is all satisfied coherence request, therefore the relative weighting of each factor
Degree of membership calculates：
Multidigit uses frequency number analysis, the danger to the indices being evaluated by evaluate collection to servicing bridges indices
Degree is graded, and the degree of membership of set of factors is obtained：
It determines and judges Subject Matrix：
By obtaining the relative defects matrix of k-th of set of factors：
Wherein：
In formula：RkThe relative defects matrix of-k-th set of factors；
rkijI-th of factor of-k-th set of factors belongs to the degree of membership of the j in evaluate collection；
pkij- group membership is rated the frequency of j to i-th of factor index of k-th of set of factors；
Construct fuzzy matrix for assessment：
By the weight vector of each indexFuzzy matrix for assessment B can be constructed with matrix R,
Calculate Comprehensive Evaluation result：
By the parameter column vector of fuzzy matrix for assessment B and evaluate collection, Comprehensive Evaluation result Z can be acquired；
Z=BV
Arrive as available from the above equation fuzzy overall evaluation further according to opinion rating as a result, provides, it is more can to evaluate servicing bridges
Factor failure risk size.
Fig. 3 is that the phenogram of planar disfigurement provided in an embodiment of the present invention illustrates intention；
Fig. 4 is stress range distribution straight line schematic diagram provided in an embodiment of the present invention.
In the above-described embodiments, can come wholly or partly by software, hardware, firmware or any combination thereof real
It is existing.When using entirely or partly realizing in the form of a computer program product, the computer program product include one or
Multiple computer instructions.When loading on computers or executing the computer program instructions, entirely or partly generate according to
Process described in the embodiment of the present invention or function.The computer can be general purpose computer, special purpose computer, computer network
Network or other programmable devices.The computer instruction may be stored in a computer readable storage medium, or from one
Computer readable storage medium is transmitted to another computer readable storage medium, for example, the computer instruction can be from one
A web-site, computer, server or data center pass through wired (such as coaxial cable, optical fiber, Digital Subscriber Line (DSL)
Or wireless (such as infrared, wireless, microwave etc.) mode is carried out to another web-site, computer, server or data center
Transmission).The computer-readable storage medium can be any usable medium or include one that computer can access
The data storage devices such as a or multiple usable mediums integrated server, data center.The usable medium can be magnetic Jie
Matter, (for example, floppy disk, hard disk, tape), optical medium (for example, DVD) or semiconductor medium (such as solid state hard disk Solid
State Disk (SSD)) etc..
The above is only the preferred embodiments of the present invention, and is not intended to limit the present invention in any form,
Any simple modification made to the above embodiment according to the technical essence of the invention, equivalent variations and modification, belong to
In the range of technical solution of the present invention.
Claims (10)
1. a kind of bridge mass analysis method based on image analysis, which is characterized in that the bridge matter based on image analysis
Analysis method includes：
Using the surface for the detecting head ultrasonic listening bridge that acoustic detector integrates, in detection, point that detecting head is internally integrated
To obtaining, the deformation data on bridge surface, the convergence of crack data modulus value solves deformation data to analysis module, crack data surmounts
Equation；For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 is of equal value；It takes
In mould, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, has zero in section [d, e]
Point s；With modulus value restrain solve f (x)=0 equation when, first find out | f (x) | minimum point；Either segment minizone in x-axis
10 parts of equal part of section [m, n] is had 11 Along ents, respectively by [m, n]
Compare on this 11 nodes | f (x) | value, find out Along ent x locating for minimum value0；If not having in [m, n] | f (x) | it is minimum
It is worth point, then x0=m or x0=n；If having in [m, n] | f (x) | minimum point, then x0It takesSome value therein；Work as x0≠ m and x0≠ n has in [m, n] | f
(x) | minimum point；With with x0Adjacent Along ent is endpoint, takes sectionBy its equal part
10 parts, equally find out Along ent x locating for minimum value1, and | f (x1)|≤|f(x0)|；Above procedure is constantly repeated, a system is obtained
The point x of column0,x1,x2,...；These points are numerical solution of the minimum point under different accuracy in the section [m, n]；
The data detected are sent to the computer being attached thereto, computer carries out the evaluation of servicing bridges Failure Factors：It loses
Imitating Assessment Curves equation is：The equation of vertical line is： Value depend on Reinforced Concrete Materials characteristic：It is calculated as followsValue：
Computer also carries out the fatigue assessment of planar disfigurement, is first depending on fatigue crack growth rate da/dN and crack tip stress is strong
Spend factor variations amplitude, ao K relational expressionDetermine the propagation of the fatigue crack in defined cycle period
And final size；Then according to the criterion that provides, judge whether the planar disfigurement can deformation occurs and fatigue fracture；
Computer in the form of images shows data that, when dangerous parameter occurs, alarm device is alarmed.
2. the bridge mass analysis method based on image analysis as described in claim 1, which is characterized in that detecting head is internally integrated
Analysis module to obtain the deformation data on bridge surface, crack data modulus value restrains and solves deformation data, crack data
Transcendental equation specifically includes：
1) it is restrained with modulus value and solves transcendental equation
For unitary transcendental equation f (x)=0, the mould of f (x) | f (x) |, Equation f (x)=0 and equation | f (x) |=0 is of equal value；
In modulus, | f (x) | perseverance is more than or equal to zero, in section [a, b], | f (x) | there is non-zero minimum point c, at section [d, e]
Inside there is zero point s；
With modulus value restrain solve the equation when, first find out | f (x) | minimum point；Either segment minizone [m, n] in x-axis,
By 10 parts of equal part of section [m, n], there are 11 Along ents, respectivelyThan
Compared on this 11 nodes | f (x) | value, find out Along ent x locating for minimum value0；If not having in [m, n] | f (x) | minimum
Point, then x0=m or x0=n；If having in [m, n] | f (x) | minimum point, then x0It takesSome value therein；Work as x0≠ m and x0≠ n has in [m, n] | f
(x) | minimum point；With with x0Adjacent Along ent is endpoint, takes sectionBy its equal part
10 parts, equally find out Along ent x locating for minimum value1, and | f (x1)|≤|f(x0)|；Above procedure is constantly repeated, a system is obtained
The point x of column0,x1,x2,...；These points are numerical solution of the minimum point under different accuracy in the section [m, n]；
2) after finding out minimum point, zero point and non-zero points then are distinguished from these minimum points；
When minimum point is non-zero points, section [a, b] equal part is obtained into a series of point x0,x1,x2..., and
|f(x0)|≥|f(x1)|≥|f(x2)|≥...≥|f(c)|,| f (c) | it is normal greater than zero
Number,For a limited big constantWherein
When minimum point is zero point, section [d, e] equal part is obtained into a series of point x0,x1,x2..., and|f
(x0)|≥|f(x1)|≥|f(x2)|≥...≥0,
3) analysis for comparing zero point and non-zero points, is obtaining minimum point x0,x1,x2... after, according toValue sentence
Whether disconnected minimum is zero point, chooses a particular value M, after restraining n step, ifThen extreme point is zero point；
M and convergent step number are depending on different problems；In the case where 10 parts of equal part,It is zero point；Modulus value |
f(xn) | convergence rate and initial section [m, n] equal part quantity of size and [m, n] for being taken it is related, initial section is got over
It is small, equal part is more, restrain it is faster.
3. the bridge mass analysis method based on image analysis as claimed in claim 2, which is characterized in that detecting head is internally integrated
Analysis module to obtain the deformation data on bridge surface, crack data modulus value restrains and solves deformation data, crack data
Transcendental equation further comprises：
The smeared cracking solved in real wave number and empty wave number plane is restrained with modulus value：
The Dispersion Equation that wave is propagated in different structure has a different concrete forms, and form is about the super of frequencies omega and wave number ξ
More equation：
G (ω, ξ)=0
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is real number and pure imaginary number；Since Dispersion Equation is in frequencies omega
It is curve with the solution in the plane of wave number ξ, with line element to entire flat scanning；
Any of first fixed frequency ω and wave number ξ, if fixed wave number is ξ0, equation g (ω, ξ)=0 becomes at this time：
F (ω)=g (ω, ξ0)=0
With the line infinitesimal of a length of 3t to straight line ξ=ξ0It is scanned, scan start point ω0；In section [ω0,ω0+ 3t] in, than
Compared with node ω0,ω0+t,ω0+2t,ω0+ 3t locates modulus value | f (ω0)|,|f(ω0+t)|,|f(ω0+2t)|,|f(ω0+ 3t) |
Size finds out minimum value；In ω0+ t or ω0+ 2t obtains minimum value, then further equal part, and judges whether it is zero point；If minimum
Value takes on endpoint, then enters next section and scan；To prevent minimum point to be by chance in endpoint ω0The place+3t, next area
Between be taken as [ω0+2t,ω0+ 5t], this section includes ω0+3t；Adjacent section is separated by 2t；It is adjacent if dividing n parts for section etc.
Section is separated by n-1 parts；When to straight line ξ=ξ0The end of scan, into lower straight line ξ=ξ0+ Δ ξ scanning, multiple scanning is multiple,
Obtain the non trivial solution in entire ω, ξ plane；
The smeared cracking solved in complex wave number field is restrained with modulus value：
The numerical solution equation g (ω, ξ)=0 in the case where wave number ξ is plural number；ξ=a+bi is enabled, wherein a, b is real number, then side
Journey g (ω, ξ)=0 is turned to：
H (ω, a, b)=g (ω, a+bi)
Since solution of the Dispersion Equation in the space that real frequency ω and complex wave number ξ is formed is curve, entire space is swept with face element
It retouches；Any of first fixed frequency ω, wave number real part a, wave number imaginary part b, if fixed frequency is ω0, at this time equation h (ω, a,
B)=g (ω, a+bi) becomes：
Q (a, b)=h (ω0,a,b)
In plane ω=ω0It is interior, plane is scanned with face infinitesimal [a, a+3t] × [b, b+3s]；Scan start point is (a0,
b0)：
In face infinitesimal [a0,a0+3t]×[b0,b0+ 3s] in, compare 16 node (a0,b0), (a0+t,b0) ..., (a0+3t,b0+
3s) place | q (a, b) | size, find out minimum value；If at the internal node of face element, i.e. (a0+t,b0+ s), (a0+2t,b0+
S), (a0+t,b0+ 2s), (a0+2t,b0+ 2s), minimum value is obtained, then further equal part, and judge whether it is zero point；
If minimum value takes on the boundary node of face element, enters next section and scan；To prevent minimum point to be by chance in
Face element boundary, next section are taken as [a0+2t,a0+5t]×[b0,b0+ 3s], when along the scanning of a axis, adjacent face element is separated by
2t；If face element divides n parts along a axis etc., adjacent face element is separated by n-1 parts；It scans, obtains whole after a axis end of scan, then along b axis
A plane ω=ω0Interior equation q (a, b)=h (ω0, a, b) solution；
To plane ω=ω0The end of scan, into next plane ω=ω0+ Δ ω scanning, multiple scanning are obtained entire
Equation q (a, b)=h (ω in the space ω, a, b0, a, b) solution.
4. the bridge mass analysis method based on image analysis as claimed in claim 3, which is characterized in that detecting head is internally integrated
Analysis module to obtain the deformation data on bridge surface, crack data modulus value restrains and solves deformation data, crack data
Transcendental equation further comprises：
The solution of any first transcendental equation is solved with modulus value convergence：
Assuming that n member transcendental equation is：
f(x1,x2,...,xn)=0
To have scanning infinitesimal and solution curve and only one intersection point, if solution curve is that m ties up (m≤n), scans infinitesimal and be selected as n-
M dimension, fixed x1,x2,...,xnIn m amount；If fixed x1,x2,...,xm, with [xm+1,xm+1+Δxm+1]×[xm+2,xm+2+Δ
xm+2]×...×[xn,xn+Δxn] as scanning infinitesimal scanning；When solution curve is that m ties up (m≤n), if scanning infinitesimal dimension is greater than
N-m dimension, then the intersection point of each infinitesimal and solution curve is not unique, and many solutions are missed；If scanning infinitesimal dimension to tie up less than n-m,
Each infinitesimal hardly intersects with solution curve；In three dimensions, solution curve is one-dimensional curve, then is scanned with plane infinitesimal；
If solution curve is Two-dimensional Surfaces, equation q (a, b)=h (ω is fixed0, a, b) in two amount, with line infinitesimal scan；If solution is bent
Line is discrete point, then is scanned with body infinitesimal；If solution curve is three-dimensional space, the equation of each point in direct clarifying space
Whether modulus value converges to zero；When practical solution n member transcendental equation, if the value of uncertain solution curve dimension m, will scan the dimension of infinitesimal
Degree is tested from high to low, and obtained solution reduces and fine and close suddenly disappearing up to solution with scanning infinitesimal dimension, and scanning dimension is
N-m-1 dimension, determines the value of m.
5. the bridge mass analysis method based on image analysis as described in claim 1, which is characterized in that computer carries out in-service
In bridge evaluation, further comprise：
Establish set of factors：
The various parameters compositing factor set for influencing servicing bridges, takes the second level factor U={ u1,u2,u3}={ defect, material are interior
Portion's environment }, u1={ u11,u12The fracture failure of }={ planar disfigurement, the fatigue failure of planar disfigurement }, u2={ u21,u22,u23}
={ material processing quality, the mechanical property of materials, design rationality }, u3={ u31,u32,u33}={ pressure, temperature, corrosivity }；
Establish evaluate collection：
In order to quantitative analysis is carried out to each evaluation index it needs to be determined that each index evaluate collection, using the evaluation of 5 grades of hundred-mark system commenting
Valence collection V divides 5 opinion ratings, i.e. V={ v1, v2, v3, v4, v5}={ is minimum, very little, small, larger, big }, wherein v1It is in-service
Bridge multifactor failure risk is minimum, and scoring section is 90~100, intermediate value 95；v 2For risk very little, section of scoring is
80~89, intermediate value 84.5；The rest may be inferred；Select the intermediate value in each section as the parameter of grade, then ginseng corresponding to 5 grades
Number is { 95,84.5,74.5,64.5,49.5 }, and parameter column vector is ν={ 95,84.5,74.5,64.5,49.5 }T；
Establish weight sets：Construct multilevel iudge matrix two-by-two：
In terms of bridge security, importance degree assignment is successively carried out to each element according to 1~9 scaling law between any two, constructed
Judgment matrix U=(uij)n×n, wherein uijExpression factor uiAnd ujRelative to the importance value of rule layer, matrix U has property：uii
=1, uij=1/uji, i, j=1,2 ..., n obtain judgment matrix：Matrix X1~X5 is pressed into row normalization, i.e.,：
Degree of membership calculates：Using frequency number analysis, to the indices being evaluated by evaluate collection to servicing bridges indices
Degree of danger is graded, and the degree of membership of set of factors is obtained；
Construct fuzzy matrix for assessment：
By the weight vector of each indexFuzzy matrix for assessment B can be constructed with matrix R,
Calculate Comprehensive Evaluation result：
By the parameter column vector of fuzzy matrix for assessment B and evaluate collection, Comprehensive Evaluation result Z can be acquired；
Z=BV
Arrive as available from the above equation fuzzy overall evaluation further according to opinion rating as a result, provides, it is multifactor can to evaluate servicing bridges
Fail risk size.
6. a kind of computer for realizing the bridge mass analysis method described in Claims 1 to 5 any one based on image analysis
Program.
7. a kind of Information Number for realizing the bridge mass analysis method described in Claims 1 to 5 any one based on image analysis
According to processing terminal.
8. a kind of computer readable storage medium, including instruction, when run on a computer, so that computer is executed as weighed
Benefit requires the bridge mass analysis method described in 1-5 any one based on image analysis.
9. a kind of bridge based on image analysis for realizing the bridge mass analysis method described in claim 1 based on image analysis
Quality analysis system, which is characterized in that the bridge quality analysis system based on image analysis is provided with：Alarm and sound wave
Detector；
Alarm is provided with shell, LED light, jingle bell；Jingle bell is embedding in the shell, LED light around jingle bell and it is embedding on the shell；
Detector is provided with operator control panel, handle, stock, detecting head；
The upper end of stock connects handle by nut；One end of handle is fixedly connected operator control panel, the bottom of stock by nut
It is welded with detecting head.
10. the bridge quality analysis system based on image analysis as claimed in claim 9, which is characterized in that acoustic detector
External computer, for the data that detected to be shown in the form of images；
Alarm connects power supply and is built-in with reserve battery, for working on when accidental power failure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810597256.3A CN108828076A (en) | 2018-06-12 | 2018-06-12 | A kind of bridge quality analysis system and method based on image analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810597256.3A CN108828076A (en) | 2018-06-12 | 2018-06-12 | A kind of bridge quality analysis system and method based on image analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108828076A true CN108828076A (en) | 2018-11-16 |
Family
ID=64145062
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810597256.3A Pending CN108828076A (en) | 2018-06-12 | 2018-06-12 | A kind of bridge quality analysis system and method based on image analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108828076A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112069447A (en) * | 2020-08-25 | 2020-12-11 | 南京航空航天大学 | Numerical method for rapidly drawing dispersion curve in complex wave number domain in waveguide structure |
| CN112069447B (en) * | 2020-08-25 | 2024-04-26 | 南京航空航天大学 | Numerical method for rapidly drawing dispersion curve in complex wave number domain in waveguide structure |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103868985A (en) * | 2014-03-20 | 2014-06-18 | 西南石油大学 | Defect quantified comprehensive safety judgment method for on-service pressure container |
| CN104504769A (en) * | 2014-12-03 | 2015-04-08 | 上海市政工程设计研究总院（集团）有限公司 | Electronic disease collection system for rail transit bridge and tunnel structure and collection and analysis method |
| CN105021699A (en) * | 2015-07-16 | 2015-11-04 | 无锡市崇安区科技创业服务中心 | Bridge pier crack nondestructive detection device |
| CN105067239A (en) * | 2015-07-31 | 2015-11-18 | 电子科技大学 | Beam crack fault detection apparatus and apparatus based on frequency sweep frequency sweep excitation vibration |
| CN105590025A (en) * | 2015-12-16 | 2016-05-18 | 南京航空航天大学 | Method for solving numerical value of dispersion curve in sonic sensor |
| CN106442720A (en) * | 2016-08-29 | 2017-02-22 | 重庆交通大学 | Vibroacoustic railway bridge healthiness monitoring apparatus, system and method |
| CN206497078U (en) * | 2017-02-13 | 2017-09-15 | 南阳博湘源工贸有限公司 | A kind of ultrasonic flaw detecting device |
| JP6236413B2 (en) * | 2015-06-22 | 2017-11-22 | 東電設計株式会社 | Deformation monitoring method for road bridge deck |
-
2018
- 2018-06-12 CN CN201810597256.3A patent/CN108828076A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103868985A (en) * | 2014-03-20 | 2014-06-18 | 西南石油大学 | Defect quantified comprehensive safety judgment method for on-service pressure container |
| CN104504769A (en) * | 2014-12-03 | 2015-04-08 | 上海市政工程设计研究总院（集团）有限公司 | Electronic disease collection system for rail transit bridge and tunnel structure and collection and analysis method |
| JP6236413B2 (en) * | 2015-06-22 | 2017-11-22 | 東電設計株式会社 | Deformation monitoring method for road bridge deck |
| CN105021699A (en) * | 2015-07-16 | 2015-11-04 | 无锡市崇安区科技创业服务中心 | Bridge pier crack nondestructive detection device |
| CN105067239A (en) * | 2015-07-31 | 2015-11-18 | 电子科技大学 | Beam crack fault detection apparatus and apparatus based on frequency sweep frequency sweep excitation vibration |
| CN105590025A (en) * | 2015-12-16 | 2016-05-18 | 南京航空航天大学 | Method for solving numerical value of dispersion curve in sonic sensor |
| CN106442720A (en) * | 2016-08-29 | 2017-02-22 | 重庆交通大学 | Vibroacoustic railway bridge healthiness monitoring apparatus, system and method |
| CN206497078U (en) * | 2017-02-13 | 2017-09-15 | 南阳博湘源工贸有限公司 | A kind of ultrasonic flaw detecting device |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112069447A (en) * | 2020-08-25 | 2020-12-11 | 南京航空航天大学 | Numerical method for rapidly drawing dispersion curve in complex wave number domain in waveguide structure |
| CN112069447B (en) * | 2020-08-25 | 2024-04-26 | 南京航空航天大学 | Numerical method for rapidly drawing dispersion curve in complex wave number domain in waveguide structure |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Omar et al. | Condition assessment of reinforced concrete bridges: Current practice and research challenges | |
| CN107036581B (en) | bridge deformation online monitoring system and method based on MEMS gyroscope | |
| Mariniello et al. | Structural damage detection and localization using decision tree ensemble and vibration data | |
| Nair et al. | Acoustic emission pattern recognition in CFRP retrofitted RC beams for failure mode identification | |
| CN113900381B (en) | Steel structure remote health monitoring platform based on Internet of things and application method | |
| Yang et al. | Damage identification of frame structure based on approximate Metropolis–Hastings algorithm and probability density evolution method | |
| KR102098888B1 (en) | System and method for disaster prediction and analysis of structures, and a recording medium having computer readable program for executing the method | |
| CN109490072A (en) | A kind of civil engineering work detection system and its detection method | |
| CN108956638A (en) | A kind of evaluation detection system for civil engineering structure visual defects | |
| Gomes et al. | Inverse structural damage identification problem in CFRP laminated plates using SFO algorithm based on strain fields | |
| CN105426646A (en) | Bridge intelligent degree comprehensive evaluation index selection and quantification method | |
| CN109799155A (en) | A kind of existing bituminous pavement damage detecting method | |
| CN106257541A (en) | A kind of modification method of bridge finite element model | |
| CN113379238A (en) | Risk assessment method and device and electronic equipment | |
| Schoefs et al. | Quantitative evaluation of contactless impact echo for non-destructive assessment of void detection within tendon ducts | |
| CN110188399B (en) | Dam safety monitoring single-measuring-point evaluation method based on multiple correlation sequences | |
| CN108168923A (en) | A kind of assay method of concrete gravity danger dam dam break risk | |
| Suryanita et al. | Application of neural networks in bridge health prediction based on acceleration and displacement data domain | |
| Wu et al. | Dam health diagnosis and evaluation | |
| KR101227776B1 (en) | Designing method for displacement estimation on construction using ann and estimating method for displacement of construction | |
| CN108458752A (en) | A kind of control system in damping bridge | |
| CN108828076A (en) | A kind of bridge quality analysis system and method based on image analysis | |
| CN112883480B (en) | Bridge load test evaluation method using inclination angle and displacement combined test | |
| CN111339710B (en) | Concrete solid structure early strength integral judgment method | |
| CN115183742A (en) | Building stability monitoring method and device based on acceleration tilt angle sensor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181116 |
|
| RJ01 | Rejection of invention patent application after publication |