CN107610112A - A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate - Google Patents
A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate Download PDFInfo
- Publication number
- CN107610112A CN107610112A CN201710819410.2A CN201710819410A CN107610112A CN 107610112 A CN107610112 A CN 107610112A CN 201710819410 A CN201710819410 A CN 201710819410A CN 107610112 A CN107610112 A CN 107610112A
- Authority
- CN
- China
- Prior art keywords
- mrow
- curved
- basic coordinates
- coordinates face
- coordinate
- 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
Landscapes
- Working Measures On Existing Buildindgs (AREA)
Abstract
The invention discloses a kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate, mesh generation is carried out to basic coordinates face first, obtains several pieces curved surface；Then modal equations are established, calculate the camber mode in basic coordinates face, the coordinate surface after measurement structure damage；Finally define damage degree index.The Curvature varying in basic coordinates face differentiates to the degree of injury of curved-surface building before and after the present invention is damaged using building, can be accurately positioned out the position of damage.This method can effecting reaction structure local damage situation, the situation of change of the physical parameters such as its research structure damage stress of near zone, rigidity can be used, and then complete the content that other kinds of monitoring method can not be completed.All key points on building can be monitored simultaneously, improve monitoring efficiency, save cost.
Description
Technical field
The invention belongs to the technical field of building non-destructive tests, and in particular to a kind of song based on digitlization modal coordinate
Face building damnification recognition method.
Background technology
Because the heat dispersion of curved-surface building thing is good, at present, the cooling tower of chemical plant and steam power plant is using curved surface
Profile.Because curved-surface building thing structure can accumulate trickle damage in long service, it is whole that these trickle damages can develop into harm
The macroscopic damage of body structure reliability service and safety.The existing disease using the method that sensor is monitored to curved-surface building thing
Evil degree is difficult to accurate judgement, and needs the bad determination of physical parameter measured, lacks effective monitoring standard.
The content of the invention
In order to solve the above-mentioned technical problem, a kind of curved-surface building thing damage based on digitlization modal coordinate provided by the invention
Hinder recognition methods.
The technical solution adopted in the present invention is：A kind of curved-surface building thing non-destructive tests side based on digitlization modal coordinate
Method, it is characterised in that comprise the following steps：
Step 1：Mesh generation is carried out to basic coordinates face, obtains several pieces curved surface；
Step 2：Establish modal equations；
Step 3：Calculate the camber mode in basic coordinates face；
Step 4：Coordinate surface after measurement structure damage；
Step 5：Damage degree index is defined, quantitatively calculates the degree of injury of building.
It is an advantage of the invention that：
1：The Curvature varying that front and rear basic coordinates face is damaged using building differentiated to the degree of injury of curved-surface building,
The position of damage can be accurately positioned out.This method can effecting reaction structure local damage situation, can use it study knot
Structure damages the situation of change of the physical parameters such as the stress of near zone, rigidity, and then the other kinds of monitoring method of completion can not
The content of completion；
2：All key points on building can be monitored simultaneously, improve monitoring efficiency, save cost.
Brief description of the drawings
Fig. 1 is the installation drawing of the embodiment of the present invention；
Fig. 2 is the division figure in the basic coordinates face of the embodiment of the present invention；
In figure, 1 is image mechanism, and 2 be curved-surface building, and 3 be curved surface where the axis of curved-surface building.
Embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, below in conjunction with the accompanying drawings and embodiment is to this hair
It is bright to be described in further detail, it will be appreciated that implementation example described herein is merely to illustrate and explain the present invention, not
For limiting the present invention.
See Fig. 1 and Fig. 2, a kind of curved-surface building thing non-destructive tests side based on digitlization modal coordinate provided by the invention
Method, comprise the following steps：
Step 1：Mesh generation is carried out to basic coordinates face；
First, using the mode of curved-surface building thing as basic coordinates face, it is the substantially intrinsic state of structure.Then, will
Basic coordinates face is horizontally split into m parts；It is divided into n-layer in short transverse；Now, basic coordinates face is divided into m × n parts；Connect
And 4 summits of every a curved surface are set to key point；Basic coordinates face length and height are respectively l0And h0, per a curved surface
Length and height be respectively l1=l0/ m and h1=h0/n。
Step 2：Establish modal equations；
Model analysis is carried out to every a curved surface in basic coordinates face, oscillatory differential equation is as follows：
Wherein, [M] is mass matrix, and [C] is damping matrix, and [K] is stiffness matrix, and y is displacement,Represent the one of displacement
Order derivative,The second dervative of displacement is represented, F (x) is external force；Solve the above-mentioned differential equation and obtain vibration displacement mode y (x).
Step 3：The camber mode in basic coordinates face is calculated；
Horizontal direction in basic coordinates face, the displacement modes for being obtained using single order central difference method basic coordinates face are led
Number is：
y′1(i, j)=[y (i+1, j)-y (i-1, j)]/(2l1), 2≤i≤m, 1≤j≤n+1；
Then, the second dervative for the displacement modes that recycling Second-Order Central Difference method obtains basic coordinates face is：
y″1(i, j)=[y '1(i+1,j)-y′1(i-1,j)]/(2l1), 3≤i≤m-1,1≤j≤n+1；
Short transverse in basic coordinates face, the displacement modes for being obtained using single order central difference method basic coordinates face are led
Number is：
y′2(i, j)=[y (i, j+1)-y (i, j-1)]/(2h1),2≤j≤n,1≤i≤m+1；
y″2(i, j)=[y '2(i,j+1)-y′2(i,j-1)]/(2h1),3≤j≤n-1,1≤i≤m+1；
To y '1(i, j) and y '2(i, j) carries out vector and tries to achieve y '3(i,j)；To y "1(i, j) and y "2(i, j) is sweared
Y " is tried to achieve in amount summation3(i,j)。
Finally, the curvature matrix that the displacement modes in basic coordinates face are calculated is：
Step 4：Coordinate surface after being damaged using video image measurement structure；
See Fig. 1, target is set on curved-surface building 2 in advance, 3 be curved surface where the axis of curved-surface building 2；Then utilize
Image mechanism 1 carries out shooting imaging to the coordinate surface under different conditions, then changes into video gradation image.With coordinate surface vibrorecord
The multiframe consecutive image intercepted as in is research object, chooses target as tracking target to create contour feature template, passes through
Target is matched to calculate the displacement P (i, j) of key point, 1≤i≤m+1,1≤j≤n+1.Then, the method meter in step 3 is utilized
The curvature matrix for calculating coordinate surface is K2(i,j)。
Step 5：Damage degree index is defined, quantitatively calculates the degree of injury of building；
When the various damage diseases such as crack, steel bar corrosion, carbonization occurs in fabric structure, it will cause the rigidity of structure
Reduction, and then cause the curvature mutation in basic coordinates face somewhere.The coordinate surface after building damage has been calculated in step 4
Curvature matrix be K2(i,j).Defining damage degree index is：
When H (i, j) level off to 1 when, represent degree of injury it is little；When H (i, j) value is deviated considerably from 1 when,
Represent that degree of injury is big.
It should be appreciated that the part that this specification does not elaborate belongs to prior art.
It should be appreciated that the above-mentioned description for preferred embodiment is more detailed, therefore can not be considered to this
The limitation of invention patent protection scope, one of ordinary skill in the art are not departing from power of the present invention under the enlightenment of the present invention
Profit is required under protected ambit, can also be made replacement or deformation, be each fallen within protection scope of the present invention, this hair
It is bright scope is claimed to be determined by the appended claims.
Claims (6)
1. a kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate, it is characterised in that comprise the following steps：
Step 1：Mesh generation is carried out to basic coordinates face, obtains several pieces curved surface；
Step 2：Establish modal equations；
Step 3：Calculate the camber mode in basic coordinates face；
Step 4：Coordinate surface after measurement structure damage；
Step 5：Damage degree index is defined, quantitatively calculates the degree of injury of building.
2. the curved-surface building thing damnification recognition method according to claim 1 based on digitlization modal coordinate, its feature exist
In the specific implementation of step 1 includes following sub-step：
Step 1.1：Using the mode of curved-surface building thing as basic coordinates face, basic coordinates face length and height are respectively l0And h0；
Step 1.2：Basic coordinates face is horizontally split into m parts, is divided into n-layer in short transverse, now, basic coordinates face quilt
It is divided into m × n part curved surfaces；
Step 1.3：4 summits of every a curved surface are set to key point, the length and height per a curved surface are respectively l1=
l0/ m and h1=h0/n。
3. the curved-surface building thing damnification recognition method according to claim 1 based on digitlization modal coordinate, its feature exist
In the specific implementation process of step 2 is：Model analysis, oscillatory differential equation are carried out to every a curved surface in basic coordinates face
For：
<mrow>
<mo>&lsqb;</mo>
<mi>M</mi>
<mo>&rsqb;</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>&rsqb;</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>K</mi>
<mo>&rsqb;</mo>
<mi>y</mi>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, [M] is mass matrix, and [C] is damping matrix, and [K] is stiffness matrix, and y is displacement,Represent that the single order of displacement is led
Number,The second dervative of displacement is represented, F (x) is external force；Solve the above-mentioned differential equation and obtain vibration displacement mode y (x).
4. the curved-surface building thing damnification recognition method according to claim 2 based on digitlization modal coordinate, its feature exist
In the specific implementation of step 3 includes following sub-step：
Step 3.1：Horizontal direction in basic coordinates face, the displacement modes in basic coordinates face are obtained using single order central difference method
Derivative be：
y′1(i, j)=[y (i+1, j)-y (i-1, j)]/(2l1), 2≤i≤m, 1≤j≤n+1；
Then, the second dervative for the displacement modes that recycling Second-Order Central Difference method obtains basic coordinates face is：
y″1(i, j)=[y '1(i+1,j)-y′1(i-1,j)]/(2l1), 3≤i≤m-1,1≤j≤n+1；
Step 3.2：Short transverse in basic coordinates face, the displacement modes in basic coordinates face are obtained using single order central difference method
Derivative be：
y′2(i, j)=[y (i, j+1)-y (i, j-1)]/(2h1),2≤j≤n,1≤i≤m+1；
y″2(i, j)=[y '2(i,j+1)-y′2(i,j-1)]/(2h1),3≤j≤n-1,1≤i≤m+1；
To y '1(i, j) and y '2(i, j) carries out vector and tries to achieve y '3(i,j)；To y "1(i, j) and y "2(i, j) carries out vector and asked
With try to achieve y "3(i,j)；
Step 3.3：The curvature matrix that the displacement modes in basic coordinates face are calculated is：
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msubsup>
<mi>y</mi>
<mn>3</mn>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msubsup>
<mi>y</mi>
<mn>3</mn>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mfrac>
<mo>,</mo>
<mn>3</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
<mo>&le;</mo>
<mi>j</mi>
<mo>&le;</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1.</mn>
</mrow>
5. the curved-surface building thing damnification recognition method according to claim 4 based on digitlization modal coordinate, its feature exist
In the specific implementation of step 4 includes following sub-step：
Step 4.1：Target is set on curved-surface building thing, the coordinate surface under different conditions shot into using image mechanism
Picture, then change into video gradation image；
Step 4.2：Using the multiframe consecutive image intercepted in coordinate surface vibrational image as research object, target is chosen as tracking mesh
Mark to create contour feature template, the displacement P (i, j) of key point, 1≤i≤m+1,1≤j≤n+ are calculated by matching target
1；
Step 4.3：The curvature matrix K in coordinates computed face2(i,j)。
6. the curved-surface building thing damnification recognition method according to claim 5 based on digitlization modal coordinate, its feature exist
In：In step 5, defining damage degree index is：
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
When H (i, j) level off to 1 when, represent degree of injury it is little；When H (i, j) value is deviated considerably from 1 when, represent
Degree of injury is big.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710819410.2A CN107610112A (en) | 2017-09-12 | 2017-09-12 | A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710819410.2A CN107610112A (en) | 2017-09-12 | 2017-09-12 | A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107610112A true CN107610112A (en) | 2018-01-19 |
Family
ID=61063039
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710819410.2A Pending CN107610112A (en) | 2017-09-12 | 2017-09-12 | A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107610112A (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928247A (en) * | 2012-10-26 | 2013-02-13 | 重庆大学 | Equivalent modal strain energy damage identification method |
| US20150035950A1 (en) * | 2012-04-18 | 2015-02-05 | Drexel University | Integration of Digital Image Correlation with Acoustic Emission |
| CN106897543A (en) * | 2017-04-25 | 2017-06-27 | 湘潭大学 | The girder construction damnification recognition method of On Modal Flexibility Curvature matrix norm |
| CN106896156A (en) * | 2017-04-25 | 2017-06-27 | 湘潭大学 | By cross uniform load face curvature difference girder construction damnification recognition method |
| CN107543581A (en) * | 2017-08-09 | 2018-01-05 | 国网山东省电力公司经济技术研究院 | Multi-functional substation framework health monitoring and damnification recognition method |
-
2017
- 2017-09-12 CN CN201710819410.2A patent/CN107610112A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150035950A1 (en) * | 2012-04-18 | 2015-02-05 | Drexel University | Integration of Digital Image Correlation with Acoustic Emission |
| CN102928247A (en) * | 2012-10-26 | 2013-02-13 | 重庆大学 | Equivalent modal strain energy damage identification method |
| CN106897543A (en) * | 2017-04-25 | 2017-06-27 | 湘潭大学 | The girder construction damnification recognition method of On Modal Flexibility Curvature matrix norm |
| CN106896156A (en) * | 2017-04-25 | 2017-06-27 | 湘潭大学 | By cross uniform load face curvature difference girder construction damnification recognition method |
| CN107543581A (en) * | 2017-08-09 | 2018-01-05 | 国网山东省电力公司经济技术研究院 | Multi-functional substation framework health monitoring and damnification recognition method |
Non-Patent Citations (2)
| Title |
|---|
| 白羽: "梁、板、网架结构损伤诊断研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 黄文龙: "模型拱桥模态识别及损伤诊断研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104880248B (en) | A kind of contribution amount quantitative identification method in electric machine structure noise excitation source | |
| CN106897543B (en) | Beam structure damage identification method of modal compliance curvature matrix norm | |
| CN104457681B (en) | Girder structure dynamic deflection monitoring method based on strain mode | |
| CN107844622A (en) | A kind of simply supported beam damage recognition methods based on faulted condition uniform load face curvature | |
| CN101833596B (en) | Two-damage quantitative identification method of rectangular and round section beam structure | |
| CN106248335B (en) | The monitoring system and real-time diagnosis method of non-navigation bridge pier anticollision | |
| CN104239736A (en) | Structure damage diagnosis method based on power spectrum and intelligent algorithms | |
| CN102759573B (en) | Based on the construction damage positioning of frequency change and the appraisal procedure of degree of injury | |
| CN112084583B (en) | Rotor blade life prediction method and system based on digital twin | |
| CN111488678A (en) | Wind turbine generator tower accumulated fatigue damage assessment system and method | |
| CN106052743B (en) | A kind of sensor mass of assessing influences frequency response function the method for size | |
| CN104634785A (en) | Novel concrete plane crack prediction, evaluation and graphical output method | |
| CN117212077B (en) | Wind wheel fault monitoring method, device and equipment of wind turbine and storage medium | |
| JP3837099B2 (en) | Structure damage estimation system and program | |
| CN107609270A (en) | A kind of ellipsoid building damnification recognition method based on digitlization modal coordinate | |
| CN107454925A (en) | For the method for the remaining life for determining wind energy plant | |
| US20230038640A1 (en) | Method for calculating service life of material under action of thermal shock load | |
| CN110780128A (en) | Sensitive equipment voltage sag fault probability evaluation method | |
| CN107345875A (en) | The rigidity appraisal procedure of Reinforced Concrete Members with Flexure based on sound emission Ring-down count | |
| CN107064315A (en) | Blade cracks position and depth recognition method based on wavelet analysis | |
| CN107633510A (en) | A kind of hyperboloid building damnification recognition method based on digitlization modal coordinate | |
| CN111563342A (en) | Strain gauge measuring point layout optimization method for rotor blade | |
| CN102601166A (en) | Straightening method based on support vector machine for workpiece deformation | |
| CN105241589A (en) | Robot arm strain test data processing method | |
| CN107610112A (en) | A kind of curved-surface building thing damnification recognition method based on digitlization modal coordinate |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180119 |