CN104267097A - Determination method of fan blade crack position - Google Patents
Determination method of fan blade crack position Download PDFInfo
- Publication number
- CN104267097A CN104267097A CN201410469897.2A CN201410469897A CN104267097A CN 104267097 A CN104267097 A CN 104267097A CN 201410469897 A CN201410469897 A CN 201410469897A CN 104267097 A CN104267097 A CN 104267097A
- Authority
- CN
- China
- Prior art keywords
- blade
- crackle
- radial position
- natural frequency
- crack
- 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
Abstract
The invention discloses a determination method of a fan blade crack position. The determination method comprises the following steps: S1. simulating crack faults at different positions of a blade and calculating front three-order inherent frequency; S2. drawing an inherent frequency resolving surface diagram by taking the inherent frequency of each order, the relative depth of a crack and the relative radial position of the crack as coordinate axes, and constructing a fault sample library; S3. calculating the front three-order inherent frequency of the blade when the blade is subjected to fault diagnosis, and reversely inputting into the fault sample library; and S4. solving by a contouring method to obtain the depth and the radial position of the crack. According to the method, the position of the crack of the blade can be solved by the contouring method, the fan blade crack can be timely diagnosed, the position of the fan blade crack can be quantitatively analyzed, the diagnosis speed is high and the result is reliable, and therefore, a powerful guarantee is provided for equipment maintenance.
Description
Technical field
The present invention relates to Diagnosis Technique field, more specifically relate to a kind of fan blade crack position defining method.
Background technology
Blower fan, because of features such as its blast are stable, air quantity is large, efficiency is high, noise is little, is widely used in factory, mine, tunnel, cooling tower, the ventilation of boats and ships and buildings, dust discharge and cooling, for safety in production provides basic guarantee.Draught fan impeller system is the core component of blower fan, and which is responsible electric energy conversion is the vital task of mechanical energy.Impeller subjects the combined action of centrifugal force, hydrodynamic force, vibration, the temperature difference, actuating medium etc. in the course of the work, stress state more complicated.Special in real work, the working environment of blade is more severe, not only will experience startup, shut down, also run under off-design behaviour of being everlasting.These stable and the air-flow exciting forces of the complexity that astable factor causes and the centrifugal force of change make blade produce destruction of resonating because the natural frequency of vibration is consistent with driving source frequency, in addition because of the combined influence of the factors such as manufacture, alignment error and machining stress, reduce the anti-fatigue ability of blade, cause the fatigue breakdown of blade, make it to crack, and serious threat is brought to the safe operation of whole unit, even cause the generation of major accident.
Often have very strong ground unrest in the working environment of blower fan, the status signal of the blower fan utilizing sensor to get is often along with the interference of neighbourhood noise.Existing blade cracks diagnostic method, no matter be carry out spectrum analysis based on acoustic emission or dependence to the vibration signal got, blade can only be diagnosed out whether to there is crackle, and can not be concrete the position that crackle occurs detected, the economic loss that can not effectively avoid blade cracks to bring, also increases the difficulty of plant maintenance simultaneously.
Summary of the invention
(1) technical matters that will solve
The technical problem to be solved in the present invention how to determine whether blade crackle occurs, and the position that crackle occurs.
(2) technical scheme
In order to solve the problems of the technologies described above, the invention provides a kind of fan blade crack position defining method, comprising the following steps:
S1, simulation blade diverse location generation crack fault, and calculate its first three rank natural frequency;
S2, with natural frequency, the relative depth of crackle, the fractional radial position of crackle described in every rank for coordinate axis, draw natural frequency solution surface chart, structure fault sample storehouse;
When S3, blade fault diagnosing, by calculating first three the rank natural frequency obtaining blade, it is oppositely input to described fault sample storehouse;
S4, utilize contouring method to solve, obtain radial position and the degree of depth at crackle place.
Preferably, in described step S2, the fractional radial position α of crackle is the radial position l of crackle and the ratio α=l/L of length of blade L, wherein, described radial position l calculates from the root of blade near impeller, and span is 0<l<L.
Preferably, in described step S2, the relative depth β of crackle is the degree of depth b of crackle and the ratio β=b/h of vane thickness h, and wherein said degree of depth b calculates from the suction surface of blade, and span is 0<b<h.
Preferably, in described step S2, with described crackle fractional radial position α for x-axis, with described crackle relative depth β for y-axis, with natural frequency ω described in the every rank of blade for z-axis.
Preferably, in described step S2, natural frequency solution surface is for by natural frequency ω described in the i-th rank
iwith the relation of described fractional radial position α and described relative depth β, and with described fractional radial position α and described relative depth β for variable, the curved surface fitted to.
Preferably, in described step S4, the process that contouring method solves is, natural frequency described in every rank is updated in described fault sample storehouse, draw out corresponding frequency contour map, then drawn at grade by frequency contour map described in three width, the fractional radial position of the corresponding crackle of the intersecting point coordinate obtained and relative depth, finally calculate radial position and the degree of depth of crackle.
(3) beneficial effect
The invention provides a kind of fan blade crack position defining method, compared with prior art, the present invention utilizes contouring method to solve the position of blade cracks generation, not only achieve the timely diagnosis to fan blade crackle, also achieve the quantitative test to fan blade crack position, fast and the reliable results of diagnosis speed, for plant maintenance provides powerful guarantee, for early warning and the quantitatively detection of fan blade crack fault, reduce the unplanned property stopping production because fault causes, thus reduce economic loss.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the process flow diagram of a kind of fan blade crack position defining method of the present invention;
Fig. 2 a is the structural representation of fan blade;
Fig. 2 b is A-A sectional view in Fig. 2 a;
Fig. 3 is the natural frequency solution surface schematic diagram of a kind of fan blade crack position defining method of the present invention;
Fig. 4 is the level line schematic three dimensional views of a kind of fan blade crack position defining method of the present invention;
Fig. 5 is the level line floor map of a kind of fan blade crack position defining method of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.Following examples for illustration of the present invention, but can not be used for limiting the scope of the invention.
Fig. 1 is the process flow diagram of a kind of fan blade crack position defining method of the present invention; The method that this law buys that comprises the following steps:
S1, simulation blade diverse location generation crack fault, and calculate its first three rank natural frequency;
S2, with natural frequency, the relative depth of crackle, the fractional radial position of crackle described in every rank for coordinate axis, draw natural frequency solution surface chart, structure fault sample storehouse;
When S3, blade fault diagnosing, by calculating first three the rank natural frequency obtaining blade, it is oppositely input to described fault sample storehouse;
S4, utilize contouring method to solve, obtain the degree of depth and the radial position at crackle place.
The present invention utilizes contouring method to solve the position of blade cracks generation, not only achieve the timely diagnosis to fan blade crackle, also achieve the quantitative test to fan blade crack position, diagnose speed fast and reliable results, for plant maintenance provides powerful guarantee.
In described step S2, the fractional radial position α of crackle is the radial position l of crackle and the ratio α=l/L of length of blade L, and wherein, radial position l calculates from the root of blade 1 near impeller, span is 0<l<L, as shown in Figure 2 a; The relative depth β of crackle is the degree of depth b of crackle and the ratio β=b/h of vane thickness h, and wherein said degree of depth b calculates from the suction surface 2 of blade, and span is 0<b<h, and as shown in Figure 2 b, 3 is pressure face; With described crackle fractional radial position α for x-axis, with described crackle relative depth β for y-axis, with natural frequency ω described in the every rank of blade for z-axis, utilize spline surface fitting technique matching natural frequency solution surface chart; Natural frequency solution surface is for by natural frequency ω described in the i-th rank
iwith the relation of described fractional radial position α and described relative depth β, and with described fractional radial position α and described relative depth β for variable, the curved surface fitted to, as shown in Figure 3.
When calculating fault Natural Frequency of Blade, blade can be reduced to rectangular cross section beam, then the Theory Solution of its first three rank natural frequency is:
Wherein, E is the elastic modulus of blade, and L is the length of blade, and I is the cross sectional moment of inertia of blade,
for the linear mass of blade.First three the rank natural frequency calculating fault blade can also adopt other computing method, is not limited to said method.
Preferably, the process that in described step S4, contouring method solves is, natural frequency described in first three rank is updated in described fault sample storehouse, draw out corresponding frequency contour map, as shown in Figure 4, then frequency contour map described in three width is drawn at grade, as shown in Figure 5, now x-axis represents the fractional radial position α of crackle, y-axis represents the relative depth β of crackle, ω axle is natural frequency, in plan view shown in Fig. 5, an intersection point characterizing this Crack Parameters must be had, the fractional radial position of the corresponding crackle of described intersecting point coordinate and relative depth, finally calculate radial position and the degree of depth of crackle, if do not obtain intersection point, representative does not have crackle.
Method provided by the invention achieves the quantitative test to crack position, fast and the reliable results of diagnosis speed, for plant maintenance provides powerful guarantee, method of the present invention can be used for the early warning of fan blade crack fault and quantitatively detects, reduce the unplanned property stopping production because fault causes, thus reduce economic loss.
Above embodiment is only for illustration of the present invention, but not limitation of the present invention.Although with reference to embodiment to invention has been detailed description, those of ordinary skill in the art is to be understood that, various combination, amendment or equivalent replacement are carried out to technical scheme of the present invention, do not depart from the spirit and scope of technical solution of the present invention, all should be encompassed in the middle of right of the present invention.
Claims (6)
1. a fan blade crack position defining method, is characterized in that, comprises the following steps:
S1, simulation blade diverse location generation crack fault, and calculate its first three rank natural frequency;
S2, with natural frequency, the relative depth of crackle, the fractional radial position of crackle described in every rank for coordinate axis, draw natural frequency solution surface chart, structure fault sample storehouse;
When S3, blade fault diagnosing, by calculating first three the rank natural frequency obtaining blade, it is oppositely input to described fault sample storehouse;
S4, utilize contouring method to solve, obtain radial position and the degree of depth at crackle place.
2. method according to claim 1, it is characterized in that, in described step S2, the fractional radial position α of crackle is the radial position l of crackle and the ratio α=l/L of length of blade L, wherein, described radial position l calculates from the root of blade near impeller, and span is 0<l<L.
3. method according to claim 1, it is characterized in that, in described step S2, the relative depth β of crackle is the degree of depth b of crackle and the ratio β=b/h of vane thickness h, wherein said degree of depth b calculates from the suction surface of blade, and span is 0<b<h.
4. the method according to claim 1 or 2 or 3, is characterized in that, in described step S2, with described crackle fractional radial position α for x-axis, with described crackle relative depth β for y-axis, with natural frequency ω described in the every rank of blade for z-axis.
5. method according to claim 5, is characterized in that, in described step S2, natural frequency solution surface is for by natural frequency ω described in the i-th rank
iwith the relation of described fractional radial position α and described relative depth β, and with described fractional radial position α and described relative depth β for variable, the curved surface fitted to.
6. method according to claim 1, it is characterized in that, in described step S4, the process that contouring method solves is, natural frequency described in every rank is updated in described fault sample storehouse, draws out corresponding frequency contour map, then frequency contour map described in three width is drawn at grade, the fractional radial position of the corresponding crackle of the intersecting point coordinate obtained and relative depth, finally calculate radial position and the degree of depth of crackle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410469897.2A CN104267097A (en) | 2014-09-15 | 2014-09-15 | Determination method of fan blade crack position |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410469897.2A CN104267097A (en) | 2014-09-15 | 2014-09-15 | Determination method of fan blade crack position |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104267097A true CN104267097A (en) | 2015-01-07 |
Family
ID=52158638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410469897.2A Pending CN104267097A (en) | 2014-09-15 | 2014-09-15 | Determination method of fan blade crack position |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104267097A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792516A (en) * | 2015-05-13 | 2015-07-22 | 中国科学院工程热物理研究所 | Device and method for testing structural fatigue of H-shaped vertical shaft wind turbine blade |
| CN105783799A (en) * | 2016-03-03 | 2016-07-20 | 四川升拓检测技术股份有限公司 | Ballastless track plate seam depth non-destructive detection method and equipment based on vibration |
| CN106501369A (en) * | 2016-10-11 | 2017-03-15 | 北京工业大学 | A kind of fan blade crack position determines method |
| CN107064315A (en) * | 2017-03-24 | 2017-08-18 | 北京工业大学 | Blade cracks position and depth recognition method based on wavelet analysis |
| CN108956075A (en) * | 2018-08-31 | 2018-12-07 | 天津大学 | Movable vane piece crackle inline diagnosis method |
| CN113504302A (en) * | 2021-06-30 | 2021-10-15 | 上海电气风电集团股份有限公司 | Method and system for monitoring fan blade state, electronic equipment and storage medium |
| CN114167026A (en) * | 2021-11-23 | 2022-03-11 | 哈尔滨工程大学 | Experimental device for turbine blade crack quantity online identification |
| CN114184763A (en) * | 2021-11-23 | 2022-03-15 | 哈尔滨工程大学 | Experimental device and method for online identification of crack positions of turbine blade |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1595146A (en) * | 2004-07-08 | 2005-03-16 | 西安交通大学 | Rotor transversal crack prediction method based on wavelet finite element model |
| CN1837780A (en) * | 2006-04-10 | 2006-09-27 | 西安交通大学 | Quantitative diagnosis method for rotor transverse crack by B-spline wavelet on the interval |
| CN102914594A (en) * | 2012-10-30 | 2013-02-06 | 西安建筑科技大学 | Method for quantitatively diagnosing transverse cracks of pipeline on basis of frequency contour lines |
| CN103984813A (en) * | 2014-05-09 | 2014-08-13 | 西安交通大学 | Vibration modeling and analyzing method of crack impeller structure of centrifugal compressor |
-
2014
- 2014-09-15 CN CN201410469897.2A patent/CN104267097A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1595146A (en) * | 2004-07-08 | 2005-03-16 | 西安交通大学 | Rotor transversal crack prediction method based on wavelet finite element model |
| CN1837780A (en) * | 2006-04-10 | 2006-09-27 | 西安交通大学 | Quantitative diagnosis method for rotor transverse crack by B-spline wavelet on the interval |
| CN102914594A (en) * | 2012-10-30 | 2013-02-06 | 西安建筑科技大学 | Method for quantitatively diagnosing transverse cracks of pipeline on basis of frequency contour lines |
| CN103984813A (en) * | 2014-05-09 | 2014-08-13 | 西安交通大学 | Vibration modeling and analyzing method of crack impeller structure of centrifugal compressor |
Non-Patent Citations (7)
| Title |
|---|
| 刘明军等: "基于多表达式编程与频率等高线相结合的结构裂纹检测", 《山东大学学报(理学版)》 * |
| 卢艳辉等: "轴裂纹识别的频率法", 《长春工业大学学报(自然科学版)》 * |
| 李兵等: "基于小波有限元法的悬臂梁裂纹识别的试验研究", 《机械工程学报》 * |
| 王可等: "基于小波有限元技术的单螺杆压缩机螺杆轴的裂纹识别", 《中国工程机械学报》 * |
| 陈雪峰等: "基于区间B样条小波有限元的转子裂纹定量识别", 《机械工程学报》 * |
| 陈雪峰等: "梁类结构多裂纹微弱损伤的小波有限元定量检测方法", 《机械工程学报》 * |
| 高丽宏等: "基于有限元模型的梁结构损伤识别技术研究", 《茂名学院学报》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104792516A (en) * | 2015-05-13 | 2015-07-22 | 中国科学院工程热物理研究所 | Device and method for testing structural fatigue of H-shaped vertical shaft wind turbine blade |
| CN105783799A (en) * | 2016-03-03 | 2016-07-20 | 四川升拓检测技术股份有限公司 | Ballastless track plate seam depth non-destructive detection method and equipment based on vibration |
| CN105783799B (en) * | 2016-03-03 | 2018-09-11 | 四川升拓检测技术股份有限公司 | A kind of fragment-free track slab gap depth lossless detection method and equipment based on vibration |
| CN106501369A (en) * | 2016-10-11 | 2017-03-15 | 北京工业大学 | A kind of fan blade crack position determines method |
| CN107064315A (en) * | 2017-03-24 | 2017-08-18 | 北京工业大学 | Blade cracks position and depth recognition method based on wavelet analysis |
| CN107064315B (en) * | 2017-03-24 | 2019-09-27 | 北京工业大学 | Blade cracks position and depth recognition method based on wavelet analysis |
| CN108956075A (en) * | 2018-08-31 | 2018-12-07 | 天津大学 | Movable vane piece crackle inline diagnosis method |
| CN113504302A (en) * | 2021-06-30 | 2021-10-15 | 上海电气风电集团股份有限公司 | Method and system for monitoring fan blade state, electronic equipment and storage medium |
| CN114167026A (en) * | 2021-11-23 | 2022-03-11 | 哈尔滨工程大学 | Experimental device for turbine blade crack quantity online identification |
| CN114184763A (en) * | 2021-11-23 | 2022-03-15 | 哈尔滨工程大学 | Experimental device and method for online identification of crack positions of turbine blade |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104267097A (en) | Determination method of fan blade crack position | |
| CN107064315B (en) | Blade cracks position and depth recognition method based on wavelet analysis | |
| CN101571120B (en) | Hierarchical cluster aviation pump multiple fault diagnostic method based on frequency multiplication relative energy sum | |
| CN104006908B (en) | fan energy consumption monitoring method and system | |
| CN106777783B (en) | Method for predicting blade cracks of aircraft engine | |
| CN104198138A (en) | Early warning method and system for abnormal vibration of wind driven generator | |
| Purarjomandlangrudi et al. | Acoustic emission condition monitoring: an application for wind turbine fault detection | |
| CN106501369A (en) | A kind of fan blade crack position determines method | |
| Konig et al. | Aeroacoustic Phenomena In High Pressure Centrifugal Compressors-A Possible Root Cause For Impeller Failures. | |
| D’Souza et al. | A new experimental facility for characterizing bladed disk dynamics at design speed | |
| Kumar et al. | A review: multiplicative faults and model-based condition monitoring strategies for fault diagnosis in rotary machines | |
| Biancolini et al. | An efficient approach to simulating ice accretion on 2D and 3D airfoils | |
| Ning et al. | Blade forced response prediction for industrial gas turbines: Part 2—verification and application | |
| Krivosheev et al. | Complex diagnostic index for technical condition assessment for GTE | |
| Schrade et al. | Experimental and numerical investigation of erosive change of shape for high-pressure compressors | |
| Huang et al. | System eigenvalue identification of mistuned bladed disks using least-squares complex frequency-domain method | |
| CN1288440C (en) | Rotor transversal crack prediction method based on wavelet finite element model | |
| Ramakrishnan et al. | Multi-blade row interactions in a low pressure ratio centrifugal compressor stage with a vaned diffuser | |
| CN113779728A (en) | Fatigue life acceleration test method suitable for electric air compressor | |
| Cicciotti et al. | Online performance monitoring of industrial compressors using meanline modelling | |
| Richards et al. | Unsteady acoustic forcing on an impeller due to coupled blade row interactions | |
| Petrie-Repar et al. | Analytical maps of aerodynamic damping as a function of operating condition for a compressor profile | |
| Jang et al. | Optimal design of splitters attached to turbo blower impeller by RSM | |
| CN204043936U (en) | A kind of hovermarine propulsion system test platform | |
| CN113420473A (en) | Method for predicting turbine wheel life |
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 |
Application publication date: 20150107 |
|
| RJ01 | Rejection of invention patent application after publication |