CN103592369A - In-situ crack testing method for disk parts of aero-engine - Google Patents
In-situ crack testing method for disk parts of aero-engine Download PDFInfo
- Publication number
- CN103592369A CN103592369A CN201310574570.7A CN201310574570A CN103592369A CN 103592369 A CN103592369 A CN 103592369A CN 201310574570 A CN201310574570 A CN 201310574570A CN 103592369 A CN103592369 A CN 103592369A
- Authority
- CN
- China
- Prior art keywords
- peak
- measuring point
- broad sense
- crackle
- local entropy
- 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.)
- Granted
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
The invention discloses an in-situ crack testing method for disk parts of an aero-engine. The in-situ crack testing method comprises the following steps: (1), arranging a plurality of measuring points on the disk parts of the aero-engine, and carrying out the standard modal shape testing to obtain data sequences of modal shapes of first four orders and form corresponding curves of modal shapes of different orders; (2), calculating general local entropies of modal shape data points di on the curves of the modal shapes of different orders to obtain general local entropy crack index curves of the modal shapes of different orders; (3), judging whether cracks exist at the measuring points corresponding to the maximum peak or not according to the relations among the maximum peak, the minimum peak and the secondary maximum peak of the general local entropy crack index curves of the modal shapes of different orders. The in-situ crack testing method has high in-situ performance, and is convenient for site instruction of in-situ dynamic nondestructive testing of the disk parts such as a turbine disk of the aero-engine.
Description
Technical field
The present invention relates to Structural Dynamics Dynamic Non-Destruction Measurement, be specifically related to a kind of aeromotor disk-like accessory crack detecting method in place.
Background technology
The characteristic of aeromotor is the deciding factor of performance, reliability and the cost of aircraft, and its lifetime expense accounts for the 20%-40% of whole aircraft.Moreover, because aeromotor cost is high, any technical small improvement all can cause great economic benefit, so the status monitoring of aircraft engine and Fault Quick Diagnosis are important and difficult research directions.
Nearly ten years, Europe is under the promotion and subsidy of European Union's the 6th scientific and technological framework planning and the 7th scientific and technological framework planning, successively investment surpasses 100,000,000,000 Euros and carries out science and technology innovation, in order to improve the reliability of aircraft and a plurality of research centres that the security ,Yi EASA of EASA is core, carry out one after another the investigation of Liao Dui EASA member state military and civilian aircraft culprit and treating method.Compare western developed country, aircraft industry is the short slab of China, particularly aeromotor life problems.The J-79 of the General Corporation engine using with meter Ge-21 U.S. of the same period F-4 fighter plane, its life-cycle can reach about 4000 hours, and some advanced aircrafts, as F-110 and the F-136 engine that F-15, F-22 were used, its life-cycle has reached respectively 6000-8000 and more than 13000 hours especially.By contrast, China's military secret engine life-cycle generally about 1000-2000 hour, and relies on foreign import to a certain extent.
For these problems, China is to improve operational reliability at the national economy and social development 11th Five-Year Plan outline > > of the < < People's Republic of China (PRC) by great product and great installation forecasting technique in life span, security, maintainability is classified gordian technique research as, in < < ' 12 ' outline > >, also clearly emphasized " the detection of Important Project health status, monitoring and diagnosis and dispose " importance of research.The disk-like accessories such as the turbine disk, as the core component of aeromotor, have material impact to engine life, therefore carry out aeromotor disk-like accessory monitoring structural health conditions in place and are significant.
Summary of the invention
The object of the present invention is to provide a kind of aeromotor disk-like accessory crack detecting method in place, the crackle of the disk-like accessories such as aero-engine turbine disk is carried out to dynamics nondestructive examination.
For realizing above goal of the invention, technical scheme of the present invention is as follows:
A disk-like accessory crack detecting method in place, comprises the following steps:
(1) a plurality of measuring points are set on aeromotor disk-like accessory, carry out normal modal test, obtain front quadravalence Mode Shape data sequence, and form each corresponding rank Mode Shape curve;
(2) calculate each Mode Shape data point d on each rank Mode Shape curve
ibroad sense local entropy (GLE), obtain the broad sense local entropy crackle indicative curve of each rank Mode Shape;
(3), according to the relation of peak-peak, minimum peak and inferior peak value on the broad sense local entropy crackle indicative curve of each rank Mode Shape, judge whether measuring point place corresponding to this peak-peak occurs crackle.
Technique effect of the present invention is:
1) take information entropy theory as basis, formed a kind of without the dynamics lossless detection method as a reference of data under structural health conditions;
2) utilize broad sense local entropy peak value can clearly judge damage position;
3) computing of the present invention property in place is good, simple, is convenient to the dynamics Non-Destructive Testing in place of the disk-like accessories such as on-the-spot guidance aero-engine turbine disk.
Accompanying drawing explanation
Fig. 1 is experiment test scheme schematic diagram;
Fig. 2 is for to arrange schematic diagram at bit test measuring point;
Fig. 3 is measuring point 3 test datas;
Fig. 4 is full measuring point 200-2500Hz frequency range frequency response function;
Fig. 5 is front quadravalence Mode Shape GLE crackle indicative curve.
Embodiment
Below in conjunction with accompanying drawing, content of the present invention is described in further detail:
The aeromotor disk-like accessory crack detecting method in place that the present invention proposes, comprises the following steps:
(1) a plurality of measuring points are set on aeromotor disk-like accessory, carry out normal modal test, obtain front quadravalence Mode Shape data sequence, and form each corresponding rank Mode Shape curve.
Wherein, normal modal test can employing power be hammered tap test, resonant method survey mode etc. into shape.In this specific embodiment, the hammering mode test macro that test macro is standard, as shown in Figure 1.Experimental configuration comprises one, universal portable computer, domestic Econ AVANT MI-7008 data collecting instrument is a set of, and domestic Dong Huali hammer and each of PCB-300A12 type single-axis acceleration sensors, Validity Test frequency field 10-10000Hz, for obtaining as far as possible many information, sample frequency is set to 10240Hz, so analysis frequency is about below 4000Hz, and selected sensor decay can meet test analysis requirement.
With reference to Fig. 2, at bit test measuring point, schematic diagram is set.It may be noted that this specific embodiment take turbine disk body as example, but it will be appreciated by those skilled in the art that other disk-like accessories for aeromotor, the present invention is applicable too.At turbine disk body, each tongue-and-groove place respectively arranges a measuring point, in this specific embodiment, take 47 tongue-and-grooves as indexed configuration measuring point, make one by one firmly hammer carry out mode test with acceleration transducer: with power hammer, to knock downside disk body, sensor is arranged in upside disk body, now power hammer hammer force will be passed to sensor by the multiple position of crackle, with this, obtain the quadravalence vibration shape before turbine disk body.After this measuring point has been tested, move to next measuring point and test.Each measuring point knocks firmly hammer 10 times, does linear averaging to reduce noise effect.The signal collecting is power hammer pumping signal time series and as response signal seasonal effect in time series acceleration signal.Fig. 3 (a) and Fig. 3 (b) are respectively pumping signal time series and the response signal time series of measuring point 3.By Fig. 3 (a) is done to ratio with the inputoutput data shown in Fig. 3 (b), can obtain the frequency response function of structure, as shown in Fig. 3 (c).
Can see, in analyzed frequency range, co-exist in quadravalence Mode Shape, respectively 1-47 measuring point be extracted to its 1-4 rank Mode Shape peak value, can obtain its Mode Shape curve.Wherein, whole frequency response function corresponding peaks of measuring point 1-47 are connected, obtain the quadravalence Mode Shape curve shown in Fig. 4.
(2) calculate each Mode Shape data point d on each rank Mode Shape curve
ibroad sense local entropy (GLE), obtain the broad sense local entropy crackle indicative curve of each rank Mode Shape.
Wherein, the definition of described broad sense local entropy and to obtain the detailed process of broad sense local entropy crackle indicative curve as follows:
A) in the concept of entropy, the probability of happening p of definition subevent i
iby this vibration shape amplitude and consecutive point amplitude with assess, be formulated as:
D in formula
irepresent measuring point i vibration shape amplitude, utilize P=(p
0, p
1..., p
n) represent the limited discrete probability distribution system that single-degree-of-freedom number is n, according to property of probability, obtain the probability of happening p of each subevent
i>=0 (i=0,1 ..., n), and have
i.e. all subevents probability and be 1.In the concept of entropy, probable value p
ithe entropy that is called as P distributes, entropy H (P) and p
ibetween pass be:
To after above formula normalization, obtain:
In the definition of local entropy, entropy function H is limited in the subsystem L of system P, and to event i, broad sense local entropy is expressed as:
Wherein N represents the data length of subsystem L.Thereby can obtain the broad sense local entropy expression formula for structural crack location:
D in formula
ithe vibration shape amplitude that represents measuring point i, d
jthe vibration shape amplitude that represents measuring point j adjacent with measuring point i in standard window, N is standard window length.
Consider two measuring point impacts that measuring point left and right is adjacent, choose N=3 long as standard window, if need, consider four measuring point impacts that measuring point left and right is adjacent, choose N=5 long as standard window, by that analogy.N value can effectively suppress to measure noise greatly, but larger N value chooses and can cause in mode test experiments measuring point number to become many, should choose as the case may be, and in the present embodiment, it is long as standard window that N=3 is chosen in suggestion, to reduce the measuring point data of needs.
B) Mode Shape data point d normal modal test being obtained
idifference substitution broad sense local entropy calculating formula (5), horizontal ordinate is measuring point numbering, ordinate is the local entropy of broad sense, obtains broad sense local entropy crackle indicative curve.Figure 5 shows that the broad sense local entropy crackle indicative curve of front quadravalence Mode Shape.
(3), according to the relation of peak-peak, minimum peak and inferior peak value on the broad sense local entropy crackle indicative curve of each rank Mode Shape, judge whether measuring point place corresponding to this peak-peak occurs crackle.
Particularly, this decision criteria can be expressed as: if the peak-peak of broad sense local entropy crackle indicative curve and its minimum peak difference are time peak value and the more than 1.2 times of minimum peak difference, think that measuring point place corresponding to this peak-peak place is candidate's crack position.Wherein, because low order mode in mode test is more accurate, in the broad sense local entropy crackle indicative curve of front quadravalence Mode Shape, if single order or second order meet above-mentioned decision criteria, judge that this measuring point place exists crackle, iff three rank or quadravalence meet, only using this measuring point as object of suspicion.
In Fig. 5, peak-peak and minimum peak be poor has all surpassed inferior peak and poor 1.2 times of minimum peak, therefore think that measuring point 4 and measuring point 27 places may exist crackle, but considering in quadravalence Mode Shape only has the 3rd rank Mode Shape result to provide crackle indication to measuring point 27, and one, two and quadravalence Mode Shape all measuring point 4 is provided to crackle indication, therefore assert that measuring point 4 exists crackle, and measuring point 27 is listed in object of suspicion.By Non-Destructive Testing, find really to have crackle at measuring point 4 places, conform to testing result of the present invention.
The foregoing is only preferred embodiment of the present invention, not in order to limit the present invention, all any modifications of doing within the spirit and principles in the present invention, be equal to and replace and improvement etc., within all should being included in protection scope of the present invention.
Claims (4)
1. an aeromotor disk-like accessory crack detecting method in place, comprises the following steps:
(1) a plurality of measuring points are set on aeromotor disk-like accessory, carry out normal modal test, obtain front quadravalence Mode Shape data sequence, and form each corresponding rank Mode Shape curve;
(2) calculate each Mode Shape data point d on each rank Mode Shape curve
ibroad sense local entropy (GLE), obtain the broad sense local entropy crackle indicative curve of each rank Mode Shape;
(3), according to the relation of peak-peak, minimum peak and inferior peak value on the broad sense local entropy crackle indicative curve of each rank Mode Shape, judge whether measuring point place corresponding to this peak-peak occurs crackle.
2. method according to claim 1, wherein, described normal modal test comprises that employing power hammer tap test or resonant method survey mode etc.
3. method according to claim 1, wherein, described step (2) is specially:
(a) the broad sense local entropy calculating formula for structural crack location is:
D in formula
ithe vibration shape amplitude that represents measuring point i, d
jthe vibration shape amplitude that represents measuring point j adjacent with measuring point i in standard window, N is standard window length;
(b) front quadravalence Mode Shape data point d normal modal test being obtained
isubstitution broad sense local entropy calculating formula (5), is numbered horizontal ordinate with measuring point respectively, and the local entropy of the broad sense of take is ordinate, obtains the broad sense local entropy crackle indicative curve of each rank Mode Shape.
4. method according to claim 1, wherein, in described step (3), decision criteria is: if the peak-peak of broad sense local entropy crackle indicative curve and its minimum peak difference are time peak value and the more than 1.2 times of minimum peak difference, think that measuring point place corresponding to this peak-peak place is candidate's crack position, wherein, for the broad sense local entropy crackle indicative curve of front quadravalence Mode Shape, if single order or second order meet above-mentioned decision criteria, judge that this measuring point place exists crackle; Iff three rank or quadravalence meet, only using this measuring point as object of suspicion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310574570.7A CN103592369B (en) | 2013-11-14 | 2013-11-14 | A kind of disc part of aero-engine crack detecting method in place |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310574570.7A CN103592369B (en) | 2013-11-14 | 2013-11-14 | A kind of disc part of aero-engine crack detecting method in place |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103592369A true CN103592369A (en) | 2014-02-19 |
CN103592369B CN103592369B (en) | 2016-01-06 |
Family
ID=50082593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310574570.7A Active CN103592369B (en) | 2013-11-14 | 2013-11-14 | A kind of disc part of aero-engine crack detecting method in place |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103592369B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105784936A (en) * | 2016-03-09 | 2016-07-20 | 西安交通大学 | Method and system for quickly detecting damage to composite material plate |
CN106370419A (en) * | 2016-08-17 | 2017-02-01 | 湖南科技大学 | Vibration response non-linearity based transmission shaft crack positioning and detecting method |
CN106546396A (en) * | 2016-11-24 | 2017-03-29 | 中国航空综合技术研究所 | A kind of reconstructing method for cracks can spread size in ferrimagnet |
CN107016156A (en) * | 2016-01-28 | 2017-08-04 | 波音公司 | The method that test data for failure mode in composite layer simplifies |
CN107976482A (en) * | 2017-11-23 | 2018-05-01 | 西安交通大学 | The crackle in-situ monitor system and method for disk-like accessory in life assessment test |
CN108857302A (en) * | 2018-07-28 | 2018-11-23 | 国营芜湖机械厂 | A kind of deployment method of metal parts crackle |
CN110031203A (en) * | 2019-05-15 | 2019-07-19 | 哈尔滨电气股份有限公司 | A kind of experimental verification device of wheel disc crack propagation life |
CN113252794A (en) * | 2021-06-03 | 2021-08-13 | 沈阳工业大学 | Acoustic emission crack monitoring method and system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701938A (en) * | 2009-11-23 | 2010-05-05 | 中国计量学院 | Audio detection method of defect type of workpiece adopting vibrating-mode frequency combination |
US20100262390A1 (en) * | 2009-04-10 | 2010-10-14 | University Of South Carolina | System and method for modal identification using smart mobile sensors |
JP4631515B2 (en) * | 2005-04-07 | 2011-02-16 | 横浜ゴム株式会社 | Non-destructive inspection method for FRP structures |
CN102759573A (en) * | 2012-07-25 | 2012-10-31 | 中国海洋石油总公司 | Frequency change-based structure damage positioning and damage degree evaluating method |
CN103076394A (en) * | 2013-01-06 | 2013-05-01 | 中国海洋石油总公司 | Safety evaluation method for ocean platform based on integration of vibration identification frequencies and vibration mode |
-
2013
- 2013-11-14 CN CN201310574570.7A patent/CN103592369B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4631515B2 (en) * | 2005-04-07 | 2011-02-16 | 横浜ゴム株式会社 | Non-destructive inspection method for FRP structures |
US20100262390A1 (en) * | 2009-04-10 | 2010-10-14 | University Of South Carolina | System and method for modal identification using smart mobile sensors |
CN101701938A (en) * | 2009-11-23 | 2010-05-05 | 中国计量学院 | Audio detection method of defect type of workpiece adopting vibrating-mode frequency combination |
CN102759573A (en) * | 2012-07-25 | 2012-10-31 | 中国海洋石油总公司 | Frequency change-based structure damage positioning and damage degree evaluating method |
CN103076394A (en) * | 2013-01-06 | 2013-05-01 | 中国海洋石油总公司 | Safety evaluation method for ocean platform based on integration of vibration identification frequencies and vibration mode |
Non-Patent Citations (3)
Title |
---|
WEI-XIN REN ET AL: "Output-only modal parameter identification of civil engineering structures", 《STRUCTURAL ENGINEERING AND MECHANICS》 * |
曾贤斌: "基于独立分量分析的桥梁动力参数识别研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
程林: "基于结构振动模态分析的损伤识别数值与试验研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107016156A (en) * | 2016-01-28 | 2017-08-04 | 波音公司 | The method that test data for failure mode in composite layer simplifies |
CN105784936A (en) * | 2016-03-09 | 2016-07-20 | 西安交通大学 | Method and system for quickly detecting damage to composite material plate |
CN106370419A (en) * | 2016-08-17 | 2017-02-01 | 湖南科技大学 | Vibration response non-linearity based transmission shaft crack positioning and detecting method |
CN106370419B (en) * | 2016-08-17 | 2018-11-23 | 湖南科技大学 | Transmission shaft crackle position finding and detection method based on vibratory response nonlinearity |
CN106546396A (en) * | 2016-11-24 | 2017-03-29 | 中国航空综合技术研究所 | A kind of reconstructing method for cracks can spread size in ferrimagnet |
CN107976482A (en) * | 2017-11-23 | 2018-05-01 | 西安交通大学 | The crackle in-situ monitor system and method for disk-like accessory in life assessment test |
CN107976482B (en) * | 2017-11-23 | 2020-07-10 | 西安交通大学 | System and method for monitoring cracks of disc parts in situ in service life assessment test |
CN108857302A (en) * | 2018-07-28 | 2018-11-23 | 国营芜湖机械厂 | A kind of deployment method of metal parts crackle |
CN110031203A (en) * | 2019-05-15 | 2019-07-19 | 哈尔滨电气股份有限公司 | A kind of experimental verification device of wheel disc crack propagation life |
CN110031203B (en) * | 2019-05-15 | 2021-07-27 | 哈尔滨电气股份有限公司 | Experimental verification device for service life of wheel disc crack propagation |
CN113252794A (en) * | 2021-06-03 | 2021-08-13 | 沈阳工业大学 | Acoustic emission crack monitoring method and system |
Also Published As
Publication number | Publication date |
---|---|
CN103592369B (en) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103592369B (en) | A kind of disc part of aero-engine crack detecting method in place | |
CN103076394B (en) | Safety evaluation method for ocean platform based on integration of vibration identification frequencies and vibration mode | |
CN103063528B (en) | Method for fast evaluating high temperate member residual service life on spot | |
CN109839440B (en) | Bridge damage positioning method based on static vehicle test | |
Jang et al. | Impact source localization for composite structures under external dynamic loading condition | |
KR101692595B1 (en) | An apparatus for monitoring metal pipes using pulsed eddy current detection and a method thereof | |
CN103063808A (en) | Assessment method and system for building damage after fire | |
Mohan et al. | Studies on damage detection using frequency change correlation approach for health assessment | |
Sekine et al. | Structural health monitoring of cracked aircraft panels repaired with bonded patches using fiber Bragg grating sensors | |
US20180306677A1 (en) | Structure abnormality detection device, structure abnormality detection method, storage medium, and structure abnormality detection system | |
CN105259398A (en) | Method for analyzing random error characteristic of optical fiber current transformer based on total variance | |
Zhou et al. | Improved reliability analysis method based on the failure assessment diagram | |
CN102778632A (en) | Double normalization recognition method for directly forecasting and recognizing transformer winding fault type | |
CN109670143B (en) | Method for detecting statistical law of vibration frequency domain response signals of civil engineering structure under environmental excitation | |
CN105784936A (en) | Method and system for quickly detecting damage to composite material plate | |
CN105651537A (en) | High-damage-sensitivity truss structure damage real-time monitoring system | |
CN105424554A (en) | Method for determining fatigue damage degree of metal material based on specific gravity measurement | |
Elwalwal et al. | Crack inspection using guided waves (GWs)/structural health monitoring (SHM) | |
CN104181230A (en) | Composite material plate structure damage monitoring method | |
Tao et al. | Nondestructive online detection of welding defects in track crane boom using acoustic emission technique | |
Vovna et al. | Improving efficiency of information measurement system of coal mine air gas protection | |
CN107144388B (en) | Global peak searching method for flexible rope vibration frequency | |
CN101982795B (en) | Method and system for detecting accuracy of gamma ray spectrometer | |
CN104236834A (en) | High-precision method for measuring frequency response functions of fragile parts online | |
CN103604615B (en) | Performance verification method of contact copying protection device of satellite large-aperture engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |