KR20210006107A - Damage detection method using lamb wave signal energy - Google Patents

Damage detection method using lamb wave signal energy Download PDF

Info

Publication number
KR20210006107A
KR20210006107A KR1020190081963A KR20190081963A KR20210006107A KR 20210006107 A KR20210006107 A KR 20210006107A KR 1020190081963 A KR1020190081963 A KR 1020190081963A KR 20190081963 A KR20190081963 A KR 20190081963A KR 20210006107 A KR20210006107 A KR 20210006107A
Authority
KR
South Korea
Prior art keywords
signal
damage
composite material
equation
energy
Prior art date
Application number
KR1020190081963A
Other languages
Korean (ko)
Other versions
KR102204747B1 (en
Inventor
신철진
박재민
김창식
조진연
박정선
Original Assignee
한국항공우주산업 주식회사
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 한국항공우주산업 주식회사 filed Critical 한국항공우주산업 주식회사
Priority to KR1020190081963A priority Critical patent/KR102204747B1/en
Publication of KR20210006107A publication Critical patent/KR20210006107A/en
Application granted granted Critical
Publication of KR102204747B1 publication Critical patent/KR102204747B1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/14Investigating 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 using acoustic emission techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The present invention relates to a damage detection method using the signal energy of an elastic wave which can accurately predict the position of a damaged portion of a composite material in comparison to a conventional technique. The damage detection method using the signal energy of an elastic wave comprises: (a) a step of generating an input signal which is an elastic wave by an excitation unit installed on a composite material to propagate the input signal along the composite material, and measuring a signal propagated along the composite material by a sensor installed on the composite material; (b) a step of estimating the size of a noise signal in a measurement signal measured by the sensor by using the input signal, and then calculating a correction signal obtained by removing the noise signal from the measurement signal; (c) a step of separating signals included in the correction signal, and determining a first signal arriving at the sensor first from the excitation unit among the signals included in the correction signal; (d) a step of using the distance between the excitation unit and the sensor and the time until the first signal arrives to calculate the speed of the input signal; and (e) a step of using a signal in a case where damage to the composite material is absent and the correction signal to calculate a difference signal before and after damage, and using the time when the difference signal before and after damage arrives and the speed of the input signal calculated in the step (d) to estimate a damage position of the composite material.

Description

탄성파의 신호 에너지를 이용한 손상 탐지 방법{Damage detection method using lamb wave signal energy}Damage detection method using lamb wave signal energy

본 발명은 탄성파의 신호 에너지를 이용한 손상 탐지 방법에 관한 것으로써, 보다 상세히는 항공기의 손상 탐지에 사용될 수 있는 손상 탐지 방법에 관한 것이다.The present invention relates to a damage detection method using signal energy of an acoustic wave, and more particularly, to a damage detection method that can be used to detect damage to an aircraft.

운용중인 항공기는 공력, 항력, 중력 및 다양한 하중을 받으며 운항하고, 이러한 하중으로 인해 항공기를 구성하는 각종 구조물들에 손상이 발생할 수 있다. 이러한 항공기의 손상에 대비하기 위해, 항공기의 주기적인 검사에 대한 각종 세부 규정이 있다. G. Anderson의 논문인 "Providing best value IVHM solution for aging aircraft"에 의하면, 항공기의 유지/보수에 소용되는 비용 중 검사와 관련된 비용은 49%를 차지하고, 항공기의 유지/보수에 필요한 인력 역시 35%를 차지하고 있어, 보다 간소화되고, 항공기의 손상 검출 효율을 향상시킬 수 있으며, 유지/보수비용을 절감할 수 있는 방법이 필요한 실정이다.The aircraft in operation operates under aerodynamic, drag, gravity and various loads, and these loads can cause damage to various structures constituting the aircraft. In order to prepare for such damage to aircraft, there are various detailed regulations for periodic inspection of aircraft. According to G. Anderson's paper "Providing best value IVHM solution for aging aircraft", inspection-related costs account for 49% of the costs used for maintenance/repair of aircraft, and the manpower required for maintenance/repair of aircraft also accounts for 35%. Because it occupies more simplification, it is possible to improve the damage detection efficiency of the aircraft, and there is a need for a method to reduce maintenance / repair costs.

항공기의 주기적인 검사 방법 중 하나로 비파괴검사가 있을 수 있으며, 비파괴검사의 대표적인 예로써 방사선 투과검사(X-ray), 초음파 탐상검사 및 와전류 탐상검사 등이 있다. 최근에는 복합재에 적용할 수 있고, 앞서 일예로 설명한 비파괴 검사 방법들보다 상대적으로 비용과 검사시간이 적고 검사 효율을 높일 수 있는 유도 초음파(탄성파)를 이용한 검사 방법에 대한 연구가 진행되고 있다.One of the periodic inspection methods of aircraft may be non-destructive testing, and representative examples of non-destructive testing include radiographic testing (X-ray), ultrasonic testing, and eddy current testing. In recent years, research on an inspection method using guided ultrasonic waves (elastic waves) that can be applied to composite materials, which has relatively less cost and inspection time than the non-destructive inspection methods described above, and can increase inspection efficiency, is being conducted.

유도 초음파란 보, 파이프, 평판 등과 같이 경계 면으로 싸여 있는 구조물 매체를 따라서 전파되는 탄성파를 말하는 것으로, 평판에 전달되는 유도 초음파는 탄성파(lamb wave)이다. 탄성파를 이용한 비파괴검사 방법은 아래 수식 1과 같이 sine과 cosine의 조합으로 이루어진 정현파 신호를 입력 신호로 사용하며, 탄성파는 버스트 웨이브(burst wave)의 형태를 가지고 있다.Guided ultrasound refers to an acoustic wave propagating along a structure medium surrounded by an interface such as a beam, pipe, or flat plate, and the guided ultrasound transmitted to the flat plate is an acoustic wave (lamb wave). The nondestructive testing method using seismic waves uses a sinusoidal signal composed of a combination of sine and cosine as an input signal as shown in Equation 1 below, and the elastic wave has a form of a burst wave.

[수식 1][Equation 1]

Figure pat00001
Figure pat00001

(여기서 A는 탄성파, w는 중심 주파수, n은 버스트 웨이브의 주기)(Where A is the acoustic wave, w is the center frequency, and n is the period of the burst wave)

도 1은 n이 5일 때, 즉 버스트 웨이브의 주기가 5일 때 수식 1로 결정되는 탄성파이다. n값은 측정하고자 하는 대상(예를 들어 복합재)과 측정하고자 하는 손상 유형에 따라 달리 선택될 수 있으며, n값에 따라 도 1에 도시된 탄성파의 형상이 달라진다. 입력신호로서의 탄성파는 피에조 일렉트릭 액추에이터(Piezo-electric actuator)를 통해 생성되어 복합재상에 전파되고, 복합재상에 위치한 피에조 일렉트릭 센서(Piezo-electric sensor)로 수집된다.1 is an elastic wave determined by Equation 1 when n is 5, that is, when the period of the burst wave is 5. The n value may be selected differently according to the object to be measured (eg, composite material) and the type of damage to be measured, and the shape of the elastic wave shown in FIG. 1 varies according to the n value. The acoustic wave as an input signal is generated through a piezo-electric actuator and propagated on the composite material, and is collected by a piezo-electric sensor located on the composite material.

도 2는 복합재상에 설치된 액추에이터(10)와 센서(11), 그리고 복합재에 손상된 부분(이하 손상위치(12))을 개략적으로 도시한 것이고, 도 3a는 복합재에 손상위치(12)가 없을 때, 센서(11)에서 측정되는 신호의 파형을 도시한 것이고, 도 3b는 복합재에 손상위치(12)가 있을 때, 센서(11)에서 측정되는 신호의 파형을 도시한 것이다.FIG. 2 schematically shows the actuator 10 and the sensor 11 installed on the composite material, and the damaged part (hereinafter, the damage location 12) in the composite material, and FIG. 3A is when there is no damage location 12 in the composite material. , It shows the waveform of the signal measured by the sensor 11, Figure 3b shows the waveform of the signal measured by the sensor 11 when there is a damage position 12 in the composite material.

도 3a에 도시된 바와 같이, 복합재에 손상이 없을 경우, 센서(11)에서는 하나의 신호가 측정되며, 해당 신호인 제1신호(S1)는 도 2의 액추에이터(10)에서 센서(11)로 직접 전달되는 신호이다. 이에 반해, 도 3b에 도시된 바와 같이, 복합재에 손상이 있을 경우, 센서(11)에서는 두 개의 신호가 순차적으로 측정된다. 제1신호(S1)는 상술한 바와 같이 도 2의 액추에이터(10)에서 센서(11)로 직접 전달되는 신호이고, 제1신호(S1)보다 나중에 측정되는 제2신호(S2)는 액추에이터(10)에서 손상위치(12)를 거쳐 센서(11)로 전달되는 신호이다. 제1신호(S1)와 제2신호(S2)가 측정되는 시각이 서로 다른 이유는, 액추에이터(10)에서 입력되는 탄성파의 전달속도는 일정하되, 도 2에 도시된 액추에이터(10)와 센서(11) 사이의 거리와 액추에이터(10), 손상위치(12) 및 센서(11)를 통과하는 거리가 서로 다르기 때문이다. 제2신호(S2)는 손상위치(12)를 통해 센서(11)로 이동하기 때문에, 제1신호(S1)에 비해 진폭이 변화한다.As shown in FIG. 3A, when there is no damage to the composite material, one signal is measured by the sensor 11, and the first signal S1, which is a corresponding signal, is transferred from the actuator 10 of FIG. 2 to the sensor 11. It is a signal transmitted directly. On the other hand, as shown in FIG. 3B, when there is damage to the composite material, the sensor 11 sequentially measures two signals. As described above, the first signal S1 is a signal directly transmitted from the actuator 10 of FIG. 2 to the sensor 11, and the second signal S2 measured later than the first signal S1 is the actuator 10 ) Is a signal transmitted to the sensor 11 through the damage location 12. The reason that the time at which the first signal S1 and the second signal S2 are measured is different is that the transmission speed of the elastic wave input from the actuator 10 is constant, but the actuator 10 and the sensor ( 11) This is because the distance between the actuator 10, the damage location 12, and the distance through the sensor 11 are different. Since the second signal S2 moves to the sensor 11 through the damage location 12, the amplitude changes compared to the first signal S1.

도 3b에 도시된 제1신호(S1)와 제2신호(S2)를 통해 액추에이터(10) 및 센서(11)에서 손상위치(12)를 계산할 수 있다. 보다 구체적으로, 액추에이터(10)와 센서(11) 사이의 거리는 알고 있으므로, 액추에이터(10)에서 탄성파를 입력한 후, 센서(11)에서 제1신호(S1)를 측정할 때까지의 시간을 통해 탄성파의 속도를 구할 수 있고, 탄성파의 속도와 액추에이터(10)에서 탄성파를 입력한 후, 센서(11)에서 제2신호(S2)를 측정할 때까지의 시간을 곱해 액추에이터(10)-손상위치(12)-센서(11) 사이의 거리 d1을 구할 수 있다. 따라서 손상위치(12)가 있을 수 있는 후보위치는 액추에이터(10)와 센서(11)를 정점으로 거리의 합이 d1인 지점인 타원으로 나타낼 수 있으며, 도 2에 도시된 센서(11)와 다른 위치 적어도 두 개 이상의 센서에서 상술한 과정을 반복하여, 적어도 세 개 이상의 타원을 구해 해당 타원들의 교점을 손상위치(12)로 계산할 수 있다.The damage location 12 in the actuator 10 and the sensor 11 may be calculated through the first signal S1 and the second signal S2 shown in FIG. 3B. More specifically, since the distance between the actuator 10 and the sensor 11 is known, the time until the sensor 11 measures the first signal S1 after inputting the acoustic wave from the actuator 10 The speed of the elastic wave can be obtained, and the time until the second signal S2 is measured by the sensor 11 after the speed of the elastic wave and the actuator 10 is input, and the actuator 10-the damage location (12)-The distance d1 between the sensors 11 can be obtained. Therefore, the candidate location where the damage location 12 may be may be represented by an ellipse, which is a point where the sum of the distances d1 to the actuator 10 and the sensor 11 as a vertex is different from that of the sensor 11 shown in FIG. Position By repeating the above-described process in at least two or more sensors, at least three or more ellipses may be obtained and the intersection of the ellipses may be calculated as the damage position 12.

한편, 손상 탐지를 위해 사용되는 탄성파는 중심 주파수에 따라 전파 속도가 변하는 멀티 모드 특성 및 복합재상에 전파해 가는 동안 형상이 변화하는 이산 특성을 가진다. 또한, 탄성파가 평판형상의 복합재상에 전파될 경우 탄성파는 평판의 중심을 기준으로 상하의 움직임이 대칭을 이루는 Symmetric Mode(이하 S모드)와 비대칭을 이루는 Anti-Symmetric Mode(이하 A모드)로 이루어진다. 탄성파는 가진되는 중심 주파수에 따라, S모드와 A모드가 여러 형태로 중첩되어 나타날 수 있으며, 이러한 탄성파의 특성으로 인하여 평판형의 복합재에서 획득된 신호는 이산특성으로 인해 분화된 신호의 조합과 반사된 신호의 조합 그리고 S모드와 A모드가 조합된 다양한 신호의 합으로 측정되며, 이에 더해 외부 잡음(White-noise) 성분 또한 함께 측정될 수 있다. 도 4a는 액추에이터(10)에서 입력되는 탄성파를, 도 4b는 상술한 바와 같은 다양한 신호가 혼합되어 센서(11)에서 측정되는 신호의 예를 도시한 것으로, 도 4b에 도시된 바와 같이 센서(11)에서는 다양한 신호가 혼합되어 측정되기 때문에, 측정되는 신호 중 손상위치(12)에 따른 손상신호 또는 손상정보를 추출하는 방법이 중요하다.Meanwhile, the acoustic wave used for damage detection has a multi-mode characteristic in which the propagation speed changes according to the center frequency, and a discrete characteristic in which the shape changes while propagating on the composite material. In addition, when an elastic wave propagates on a flat composite material, the elastic wave is composed of a Symmetric Mode (hereinafter referred to as S mode) in which the vertical movement is symmetrical with respect to the center of the plate, and an Anti-Symmetric Mode (hereinafter referred to as A mode) that is asymmetric. Depending on the center frequency to be excitation, the S-mode and A-mode may be overlapped in various forms. Due to the characteristics of these elastic waves, the signals obtained from the flat-type composite material are combined and reflected by the differentiated signals due to discrete characteristics. It is measured by the combination of the signal and the sum of the various signals in which the S mode and A mode are combined. In addition, the white-noise component can also be measured. FIG. 4A shows an elastic wave input from the actuator 10, and FIG. 4B shows an example of a signal measured by the sensor 11 by mixing various signals as described above. As shown in FIG. 4B, the sensor 11 ), since various signals are mixed and measured, a method of extracting the damage signal or damage information according to the damage location 12 among the measured signals is important.

종래 센서(11)에서 측정되는 신호 중 잡음 성분을 줄이기 위해 1979년에 Spectral Subtraction 방법이 음성 신호 처리에서 제안되었으며, 이 방법은 센서(11)에서 측정된 신호 x(t)에서 잡음 성분의 신호 v(t)를 추정하고, 주파수 영역에서 추정된 잡음 신호의 에너지를 제거함으로써 센서(11)에서 측정되는 신호 중, 입력된 탄성파에 의한 신호 s(t)를 추정하는 방법이다. 이와 같은 방법은 잡음의 에너지를 줄여, 잡음의 크기를 감소시키는 방법이며, 이러한 방법은 도 5와 같이 도시할 수 있다.In order to reduce the noise component of the signal measured by the conventional sensor 11, a Spectral Subtraction method was proposed in speech signal processing in 1979, and this method is the signal v of the noise component in the signal x(t) measured by the sensor 11 This is a method of estimating (t) and estimating the signal s(t) due to the input acoustic wave among the signals measured by the sensor 11 by removing the energy of the noise signal estimated in the frequency domain. This method is a method of reducing the energy of the noise and reducing the size of the noise, and this method may be illustrated as shown in FIG. 5.

도 5에 도시된 종래의 방법을 사용하기 위해서는, 잡음 신호 v(t)의 에너지 크기를 추정해야 한다. 그러나 손상 탐지를 위해 측정되는 신호 중 잡음 신호의 크기는 대상물의 형태와 센서의 부착 상태에 따라 달라질 수 있으며, 잡음 신호의 에너지 크기를 추정하는 방법은 현실적으로 한계가 있는 실정이다.In order to use the conventional method shown in Fig. 5, it is necessary to estimate the energy level of the noise signal v(t). However, the size of the noise signal among signals measured for damage detection may vary depending on the shape of the object and the attachment state of the sensor, and the method of estimating the energy level of the noise signal is practically limited.

한국 공개특허공보 제10-2009-0012818호("비파괴 검사장치", 공개일 2009.02.04.)Korean Laid-Open Patent Publication No. 10-2009-0012818 ("Non-destructive inspection device", published date 2009.02.04.)

본 발명은 상기한 바와 같은 문제점을 해결하기 위해 안출된 것으로써, 본 발명에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법의 목적은 종래의 방식보다 정확하게 복합재상의 손상부분의 위치를 예측할 수 있는 탄성파의 신호 에너지를 이용한 손상 탐지 방법을 제공함에 있다.The present invention was conceived to solve the above problems, and an object of the method for detecting damage using signal energy of an elastic wave according to the present invention is to predict the location of a damaged part on a composite material more accurately than the conventional method. It is to provide a damage detection method using signal energy.

상기한 바와 같은 문제점을 해결하기 위한 본 발명에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법은 a) 복합재에 설치된 가진부에서 탄성파인 입력신호를 발생시켜 상기 복합재를 따라 전파되게 하고, 복합재에 설치된 센서에서 상기 복합재를 따라 전파되는 신호를 측정하는 단계; b) 상기 센서에서 측정된 측정신호에서 상기 입력신호를 이용해 잡음신호의 크기를 추정한 후, 상기 측정신호에서 상기 잡음신호가 제거된 보정신호를 계산하는 단계; c) 상기 보정신호에 포함된 신호들을 분리하고, 상기 보정신호에 포함된 신호들 중, 상기 가진부에서 상기 센서로 가장 먼저 도착한 제1신호를 판별하는 단계; d) 상기 가진부와 상기 센서의 간격과 상기 제1신호가 도달하기까지의 시간을 이용해 상기 입력신호의 속도를 계산하는 단계; 및 e) 상기 보정신호와 상기 복합재에 손상이 없을 경우의 신호를 이용해 손상 전후 차이 신호를 계산하고, 상기 손상 전후 차이 신호가 도달한 시간과 상기 d) 단계에서 계산된 상기 입력신호의 속도를 이용해 상기 복합재의 손상위치를 추정하는 단계;를 포함하는 것을 특징으로 한다.In order to solve the above problems, the method for detecting damage using signal energy of an elastic wave according to the present invention is a) generating an input signal, which is an elastic wave, from an excitation part installed in a composite material to propagate along the composite material, and a sensor installed in the composite material. Measuring a signal propagating along the composite material at; b) estimating a magnitude of a noise signal using the input signal from the measurement signal measured by the sensor, and then calculating a correction signal from which the noise signal is removed from the measurement signal; c) separating signals included in the correction signal, and determining a first signal that arrives first from the excitation unit to the sensor from among signals included in the correction signal; d) calculating the speed of the input signal using the distance between the excitation unit and the sensor and the time until the first signal arrives; And e) calculating a difference signal before and after damage using the correction signal and a signal when there is no damage to the composite material, and using the time at which the difference signal before and after damage arrives and the speed of the input signal calculated in step d). It characterized in that it comprises; estimating the damage location of the composite material.

또한, 상기 b) 단계는, b-1) 아래 수식 5를 이용해 상기 잡음신호의 크기를 추정하는 단계; 를 포함하는 것을 특징으로 한다.In addition, step b) may include: b-1) estimating the magnitude of the noise signal using Equation 5 below; It characterized in that it comprises a.

[수식 5][Equation 5]

Figure pat00002
Figure pat00002

(

Figure pat00003
는 추정된 상기 잡음신호의 크기,
Figure pat00004
는 상기 측정신호의 에너지 크기,
Figure pat00005
는 상기 입력신호의 에너지 크기)(
Figure pat00003
Is the estimated size of the noise signal,
Figure pat00004
Is the energy level of the measurement signal,
Figure pat00005
Is the energy level of the input signal)

또한, 상기 b) 단계는, b-2) 상기 b-1) 단계 이후에 수행되며, 아래 수식 6을 이용해 시간영역에서의 상기 보정신호를 계산하는 단계; 를 더 포함하는 것을 특징으로 한다.In addition, step b) is performed after step b-2) and step b-1), and calculating the correction signal in the time domain using Equation 6 below; It characterized in that it further comprises.

[수식 6][Equation 6]

Figure pat00006
Figure pat00006

(상기

Figure pat00007
는 상기 측정신호를 STFT(Short Time Fourier Transform)한 신호, 상기
Figure pat00008
는 주파수 영역에서의 상기 보정신호,
Figure pat00009
는 상기
Figure pat00010
을 inverse STFT한 시간영역에서의 보정신호)(remind
Figure pat00007
Is a signal obtained by short-time fourier transform (STFT) of the measurement signal,
Figure pat00008
Is the correction signal in the frequency domain,
Figure pat00009
Is reminded
Figure pat00010
Correction signal in time domain inverse STFT)

또한, 상기 b) 단계는, b-3) 상기 b-2) 단계 이후에 수행되고, 상기 수식 5의

Figure pat00011
을 상기 수식 6을 통해 구해진
Figure pat00012
로 대체하여 상기 보정신호에 포함된 잡음신호의 크기를 추정하고, 상기 보정신호의 크기 대비 상기 보정신호에 포함된 잡음신호의 크기가 기준치 이상일 경우, 상기 수식 6의
Figure pat00013
을 상기
Figure pat00014
로 대체하여 잡음신호의 크기가 줄어든 보정신호를 생성하는 단계; 를 더 포함하는 것을 특징으로 한다.In addition, step b) is performed after step b-3) b-2), and
Figure pat00011
Obtained through Equation 6 above
Figure pat00012
In the case where the size of the noise signal included in the correction signal is estimated and the size of the noise signal included in the correction signal is greater than or equal to the reference value,
Figure pat00013
Remind
Figure pat00014
Generating a correction signal having a reduced noise signal by replacing it with It characterized in that it further comprises.

또한, 상기 b-3) 단계는 생성된 보정신호 내에 포함된 잡음신호의 크기가 기준치 이하가 될 때 까지 반복 수행되는 것을 특징으로 한다.In addition, the step b-3) is characterized in that it is repeatedly performed until the size of the noise signal included in the generated correction signal becomes less than or equal to a reference value.

또한, 상기 c) 단계는 c-1) 아래 수식 7을 이용해 상기 보정신호에 포함된 신호 중 에너지가 가장 큰 신호를 분리해, 분리된 신호를 상기 제1신호로 판별하는 단계; 를 포함하는 것을 특징으로 한다.In addition, step c) may include: c-1) separating a signal having the highest energy among signals included in the correction signal using Equation 7 below, and determining the separated signal as the first signal; It characterized in that it comprises a.

[수식 7][Equation 7]

Figure pat00015
Figure pat00015

(r(t)은 보정신호에 포함된 신호 중 에너지가 가장 큰 신호,

Figure pat00016
는 시간영역에서의 보정신호를 STFT한 신호,
Figure pat00017
는 입력신호)(r(t) is the signal with the largest energy among the signals included in the correction signal,
Figure pat00016
Is the signal obtained by STFT of the correction signal in the time domain,
Figure pat00017
Is the input signal)

또한, 상기 e) 단계는, 아래 수식 7을 이용해 상기 손상 전후 차이 신호에 포함된 신호 중 하나를 분리하고, 분리된 신호가 도달한 시간과 상기 d) 단계에서 계산된 상기 입력신호의 속도를 이용해 상기 복합재의 손상위치를 추정하는 것을 특징으로 한다.In addition, in step e), one of the signals included in the difference signal before and after damage is separated using Equation 7 below, and the time when the separated signal arrives and the speed of the input signal calculated in step d) are used. It characterized in that the damage location of the composite material is estimated.

[수식 7][Equation 7]

Figure pat00018
Figure pat00018

(r(t)은 상기 손상 전후 차이 신호에 포함된 신호 중 에너지가 가장 큰 신호,

Figure pat00019
는 시간영역에서의 상기 손상 전후 차이 신호를 STFT한 신호,
Figure pat00020
는 입력신호)(r(t) is the signal with the largest energy among the signals included in the difference signal before and after damage,
Figure pat00019
Is a signal obtained by STFT of the difference signal before and after damage in the time domain,
Figure pat00020
Is the input signal)

또한, 상기 e) 단계는 반복 수행되어 상기 손상 전후 차이 신호에 포함된 소정 개수의 신호를 분리하고, 상기 e) 단계가 반복 수행될 때, 상기 수식 7에서 상기

Figure pat00021
신호는 이전에 수행된 e) 단계의
Figure pat00022
신호에서 이전에 수행된 e) 단계에서 분리된 신호를 뺀 신호인 것을 특징으로 한다.In addition, when step e) is repeatedly performed to separate a predetermined number of signals included in the difference signal before and after damage, and when step e) is repeatedly performed, in Equation 7
Figure pat00021
The signal was previously performed in step e).
Figure pat00022
It is characterized in that it is a signal obtained by subtracting the signal separated in step e) previously performed from the signal.

또한, 상기 e) 단계는 반복 수행되어 상기 손상 전후 차이 신호에 포함된 소정 개수의 신호를 분리하여, 손상위치의 영역을 추정하는 것을 특징으로 한다.In addition, the step e) is repeatedly performed to separate a predetermined number of signals included in the difference signal before and after damage to estimate an area of the damage location.

또한, 상기 복합재에 설치된 센서는 복수개인 것을 특징으로 한다.In addition, it is characterized in that a plurality of sensors installed in the composite material.

상기한 바와 같은 본 발명에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법에 의하면, b) 단계에서 미리 알 수 있는 입력신호의 에너지 크기를 이용하여 잡음신호의 크기를 추정하여 측정신호에서 보정신호를 제거하기 때문에, 종래 부정확한 측정신호의 크기를 추정하여 측정신호에서 잡음신호를 제거하는 방식보다 정확한 신호를 얻을 수 있어, 복합재의 손상부분의 위치를 보다 정확하게 예측할 수 있는 효과가 있다.According to the damage detection method using the signal energy of the acoustic wave according to the present invention as described above, the correction signal is removed from the measurement signal by estimating the magnitude of the noise signal using the energy magnitude of the input signal known in advance in step b). Therefore, it is possible to obtain a more accurate signal than the conventional method of removing the noise signal from the measurement signal by estimating the size of the incorrect measurement signal, there is an effect of more accurately predicting the location of the damaged portion of the composite.

또한, 본 발명에 의하면 수식 7을 이용하여 측정신호에서 제1신호를 분리하므로, 보다 정확하게 제1신호를 분리할 수 있어 복합재의 손상부분의 위치를 보다 정확하게 예측할 수 있는 효과가 있다.In addition, according to the present invention, since the first signal is separated from the measurement signal by using Equation 7, it is possible to more accurately separate the first signal, so that the location of the damaged part of the composite material can be more accurately predicted.

도 1은 sine과 cosine의 조합으로 이루어지는 탄성파의 파형.
도 2는 복합재상에 설치된 액추에이터, 센서, 손상위치의 개략도.
도 3은 복합재에 손상위치의 유무에 따른 탄성파 신호의 변화의 파형.
도 4는 입력신호 및 탄성파의 특성에 따라 다수개의 신호가 조합되어 측정되는 측정신호의 일예.
도 5는 종래 측정신호에서 잡음신호의 에너지 크기를 추정하여 노이즈를 제거하는 방식의 개략도.
도 6은 복합재에 설치된 가진부와 센서들의 개략도.
도 7은 본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법의 b) 단계에서 잡음신호의 크기를 추정하는 방법의 개략도.
도 8은 본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법의 c) 단계에서 제1신호를 분리한 개략도.
도 9는 본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법에서 사용되는 손상 전후 차이 신호의 파형.
도 10은 손상 전후 차이 신호의 확대 파형 및 분리신호의 파형.
1 is a waveform of an acoustic wave composed of a combination of sine and cosine.
Figure 2 is a schematic diagram of an actuator, a sensor, and a damage location installed on a composite material.
3 is a waveform of a change in an elastic wave signal according to the presence or absence of a damage location in the composite material.
4 is an example of a measurement signal that is measured by combining a plurality of signals according to characteristics of an input signal and an acoustic wave.
5 is a schematic diagram of a method of removing noise by estimating the energy level of a noise signal from a conventional measurement signal.
6 is a schematic diagram of the excitation unit and sensors installed in the composite material.
7 is a schematic diagram of a method of estimating the magnitude of a noise signal in step b) of a method for detecting damage using signal energy of an acoustic wave according to an embodiment of the present invention.
8 is a schematic diagram of separating a first signal in step c) of a method for detecting damage using signal energy of an acoustic wave according to an embodiment of the present invention.
9 is a waveform of a difference signal before and after damage used in a damage detection method using signal energy of an acoustic wave according to an embodiment of the present invention.
10 is an enlarged waveform of a difference signal before and after damage and a waveform of a split signal.

이하 첨부된 도면을 참고하여 본 발명에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법의 바람직한 실시예에 관하여 상세히 설명한다.Hereinafter, a preferred embodiment of a damage detection method using signal energy of an acoustic wave according to the present invention will be described in detail with reference to the accompanying drawings.

본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법은 순차적으로 수행되는 a) 단계, b) 단계, c) 단계, d) 단계 및 e) 단계를 포함할 수 있다.A method of detecting damage using signal energy of an elastic wave according to an embodiment of the present invention may include steps a), b), c), d) and e) sequentially performed.

먼저, 배경기술에서도 설명했듯 본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법은 복합재의 일부분이 손상되었을 때, 손상된 부분의 위치를 파악하기 위해 안출된 것으로, 복합재에 설치된 가진부와 센서를 이용할 수 있다. 복합재에 설치된 가진부는 특정 중심주파수를 가지는 탄성파를 생성하여 복합재를 따라 전파되도록 할 수 있으며, 가진부에 의해 생성된 탄성파는 배경기술에서 설명한 것과 같이 주기가 n인 버스트 웨이브(burst wave) 형태이다.First, as described in the background art, the damage detection method using the signal energy of an elastic wave according to an embodiment of the present invention was devised to determine the location of the damaged part when a part of the composite material is damaged. And sensors are available. The excitation unit installed in the composite material can generate an elastic wave having a specific center frequency to propagate along the composite material, and the elastic wave generated by the excitation unit is in the form of a burst wave having a period of n as described in the background art.

도 6은 복합재(100)에 설치된 가진부(200)와 센서들을 개략적으로 도시한 것이다.6 schematically shows the excitation part 200 and sensors installed in the composite material 100.

도 6에 도시된 바와 같이, 본 발명의 일실시예에 의한 탄성파의 신호 에너지를 이용한 손상 탐지 방법을 적용하기 위해서는 하나의 복합재(100)에 가진부(200)와 복수개의 센서들이 설치되어 있다. 보다 상세히, 도 6에 도시된 바와 같이 복합재(100)의 표면상에는 제1센서(310), 제2센서(320) 및 제3센서(330)가 설치되어 있는데, 이는 적어도 세 개 이상의 센서에서 신호를 측정해야 손상부분의 정확한 위치를 파악할 수 있기 때문이며, 이에 관해서는 후술한다.As shown in FIG. 6, in order to apply a damage detection method using signal energy of an elastic wave according to an embodiment of the present invention, an excitation unit 200 and a plurality of sensors are installed in one composite material 100. In more detail, as shown in FIG. 6, a first sensor 310, a second sensor 320, and a third sensor 330 are installed on the surface of the composite material 100, which are signals from at least three sensors. This is because it is possible to grasp the exact position of the damaged part only by measuring the value, which will be described later.

a) 단계는 복합재(100)에 설치된 가진부(200)에서 탄성파인 입력신호를 발생시켜, 상기 입력신호가 복합재(100)를 따라 전파되게 한다. 복합재(100)에서 전파된 입력신호는 도 6에 도시된 바와 같이 각각의 센서들에 집적 전달되거나, 복합재(100)의 손상부분(110)을 거쳐 각각의 센서들에 전달될 수 있으며, 센서는 복합재를 통해 전달되는 신호를 측정할 수 있다.Step a) generates an input signal that is an elastic wave in the excitation unit 200 installed in the composite material 100, so that the input signal propagates along the composite material 100. The input signal propagated from the composite material 100 may be integratedly transmitted to each of the sensors as shown in FIG. 6 or may be transmitted to each of the sensors through the damaged portion 110 of the composite material 100, and the sensor is You can measure the signal transmitted through the composite.

b) 단계는 센서들 각각에서 측정된 측정신호에서 잡음신호가 제거된 보정신호를 계산한다. 보다 구체적으로, 가진부(200)에서 발생한 입력신호는 복합재(100)를 따라 진행하면서 진폭 또는 주파수가 변화할 수 있으며, 외란에 의한 잡음신호(white-noise)가 포함될 수 있다. 입력신호에 잡음신호가 포함될 경우, 정확한 본래 입력신호에 의한 측정신호의 분산(variance)을 증가시키며, 이는 정확한 손상부분(110)의 위치를 파악하지 못하는 요인이 될 수 있으며, 측정신호에서 잡음신호를 제거하는 것이 중요하다.Step b) calculates a correction signal from which the noise signal is removed from the measurement signal measured by each of the sensors. More specifically, the input signal generated by the excitation unit 200 may change in amplitude or frequency while traveling along the composite material 100, and a white-noise signal due to disturbance may be included. When a noise signal is included in the input signal, the variance of the measurement signal due to the exact original input signal is increased, which may be a factor in which the location of the damaged part 110 cannot be accurately identified. It is important to eliminate.

본 실시예의 b) 단계는 상기한 바와 같은 기술적 과제를 해결하기 위해, 입력신호를 이용해 잡음신호의 크기를 추정한 후, 잡음신호를 측정신호에서 제거하며, b-1) 단계, b-2) 단계 및 b-3) 단계를 포함할 수 있다.In step b) of this embodiment, in order to solve the technical problem as described above, after estimating the magnitude of the noise signal using the input signal, the noise signal is removed from the measurement signal, and step b-1), b-2) It may include steps and b-3).

센서에서 측정된 측정신호인

Figure pat00023
은 잡음이 없는 탄성파의 성분
Figure pat00024
와 잡음성분
Figure pat00025
의 구성으로 표현할 수 있으며, 이는 아래와 수식 2와 같다.The measurement signal measured by the sensor
Figure pat00023
Is the component of the acoustic wave without noise
Figure pat00024
And noise component
Figure pat00025
It can be expressed as the composition of, which is shown in Equation 2 below.

[수식 2][Equation 2]

Figure pat00026
Figure pat00026

수식 2는 하기 수식 3과 같이 STFT(Short Time Fourier Transform)로 표현 가능하다.Equation 2 can be expressed as a Short Time Fourier Transform (STFT) as shown in Equation 3 below.

[수식 3][Equation 3]

Figure pat00027
Figure pat00027

또한, 상기한 수식은 이산화된 형태의 다음 수식 4로 표현 가능하다.In addition, the above equation can be expressed by the following equation 4 in a discrete form.

[수식 4][Equation 4]

Figure pat00028
Figure pat00028

본 실시예에서 신호를 STFT를 할 때 창함수(window function)는 Hann window function을 사용하였다.In this embodiment, the Hann window function was used as the window function when STFT the signal.

b-1) 단계는 아래 수식 5를 이용해, 센서에서 측정된 측정신호에 포함된 잡음신호의 크기를 추정한다. Step b-1) uses Equation 5 below to estimate the size of the noise signal included in the measured signal measured by the sensor.

[수식 5][Equation 5]

Figure pat00029
Figure pat00029

(

Figure pat00030
는 추정된 상기 잡음신호의 크기,
Figure pat00031
는 상기 측정신호의 에너지 크기,
Figure pat00032
는 상기 입력신호의 에너지 크기)(
Figure pat00030
Is the estimated size of the noise signal,
Figure pat00031
Is the energy level of the measurement signal,
Figure pat00032
Is the energy level of the input signal)

상기 수식을 포함하여, 이하 다른 수식에서 대문자는 주파수영역에서의 신호를 말하며, 소문자는 시간영역에서의 신호를 말한다. 즉, 본 실시예는 수식 5 사용하기 위해 센서에서 측정되는 시간영역에서의 측정신호인

Figure pat00033
을 수식 4를 이용해 STFT해서 주파수 영역에서의 신호인
Figure pat00034
을 구하고, 가진부(100)에서 버스트 웨이브 형태의 탄성파로 생성되는 입력신호인
Figure pat00035
을 마찬가지로 STFT하여
Figure pat00036
을 구한다. 수식 5에서 사용되는 p값과
Figure pat00037
값은 복합재의 재질에 따라 달리 결정되는 값이다.Including the above formula, in other formulas hereinafter, uppercase letters refer to signals in the frequency domain, and lowercase letters refer to signals in the time domain. That is, this embodiment is a measurement signal in the time domain measured by the sensor to use Equation 5
Figure pat00033
STFT using Equation 4 to get the signal in the frequency domain
Figure pat00034
Is obtained, and the input signal generated as a burst wave type elastic wave in the excitation unit 100
Figure pat00035
Likewise by STFT
Figure pat00036
Find The p-value used in Equation 5 and
Figure pat00037
The value is determined differently depending on the material of the composite material.

이 과정에서 입력신호인

Figure pat00038
의 길이와 측정신호인
Figure pat00039
의 길이는 서로 다르므로, 수식 4의 STFT와 수식 5를 이용하여 입력신호와 측정신호 각각의 에너지를 계산하였으며, STFT의 창함수의 크기는 입력신호의 크기와 동일하게 설정하였다. 여기서 STFT의 창함수의 크기를 입력신호의 크기와 동일하게 설정한다는 뜻은, 시간영역에서의 신호의 길이를 동일하게 한다는 의미이다. 상기한 수식 4 및 수식 5를 이용하여 잡음신호의 에너지 크기를 추정하는 과정은 도 7에 개략적으로 도시되어 있다.In this process, the input signal
Figure pat00038
Length and measurement signal
Figure pat00039
Since the lengths of are different, energy of each input signal and measurement signal was calculated using the STFT of Equation 4 and Equation 5, and the size of the window function of the STFT was set equal to the size of the input signal. Here, setting the size of the window function of the STFT equal to the size of the input signal means that the length of the signal in the time domain is the same. A process of estimating the energy level of a noise signal using Equations 4 and 5 is schematically illustrated in FIG. 7.

b-2) 단계는 상기 b-1) 단계 이후에 수행되며, 아래 수식 6을 이용해 시간영역에서의 상기 보정신호를 계산한다.Step b-2) is performed after step b-1), and the correction signal in the time domain is calculated using Equation 6 below.

[수식 6][Equation 6]

Figure pat00040
Figure pat00040

(

Figure pat00041
는 주파수 영역에서의 상기 보정신호,
Figure pat00042
는 상기
Figure pat00043
을 inverse STFT한 시간영역에서의 보정신호)(
Figure pat00041
Is the correction signal in the frequency domain,
Figure pat00042
Is reminded
Figure pat00043
Correction signal in time domain inverse STFT)

b-3) 단계는 상기한 b-1) 단계 및 b-2) 단계를 통해 구해진

Figure pat00044
에 포함된 보정신호의 비율이 기준치 이상일 경우, b-1) 단계와 b-2) 단계를 반복 수행하는 단계이다. 즉, b-1) 단계와 b-2) 단계가 한 번 수행되었을 때 계산되는
Figure pat00045
로 수식 5의
Figure pat00046
을 대체한 후,
Figure pat00047
에 포함되는 잡음신호의 크기가 기준치 이상일 경우, b-1) 단계 및 b-2) 단계를 반복 수행할 수 있으며, b-2) 단계가 반복 수행될 때, 수식 6에서
Figure pat00048
Figure pat00049
로 대체될 수 있다.Step b-3) is obtained through the above steps b-1) and b-2).
Figure pat00044
If the ratio of the correction signal included in is greater than or equal to the reference value, steps b-1) and b-2) are repeatedly performed. That is, when steps b-1) and b-2) are performed once,
Figure pat00045
As of Formula 5
Figure pat00046
After replacing
Figure pat00047
If the size of the noise signal included in is greater than or equal to the reference value, steps b-1) and b-2) may be repeatedly performed, and when step b-2) is repeatedly performed, in Equation 6
Figure pat00048
Is
Figure pat00049
Can be replaced with

c) 단계는 상기 보정신호에 포함된 신호들을 분리하고, 상기 보정신호에 포함된 신호들 중, 상기 가진부에서 상기 센서로 가장 먼저 도착한 제1신호를 판별하고, d) 단계는 가진부(100)와 센서들 사이의 간격, 그리고 가진부(100)에서 입력신호가 생성된 후, 보정신호의 제1신호가 도달하기까지의 시간을 이용해 입력신호의 속도를 계산하며, e) 단계는 보정신호와 상기 복합재에 손상이 없을 경우의 신호를 이용해 손상 전후 차이 신호를 계산하고, 상기 손상 전후 차이 신호가 도달한 시간과 상기 d) 단계에서 계산된 상기 입력신호의 속도를 이용해 상기 복합재의 손상위치를 추정한다. 상술한 c) 단계와 e) 단계의 개략적인 동작은 종래 배경기술에서 설명한 것과 동일하며, 이하 c) 단계 및 d) 단계에서 신호를 분리하는 과정에 관하여 상세히 설명한다.Step c) separates the signals included in the correction signal, and among the signals included in the correction signal, the first signal that arrives first from the excitation unit to the sensor is determined, and in step d), the excitation unit 100 ) And the distance between the sensors, and the time until the first signal of the correction signal arrives after the input signal is generated by the excitation unit 100, and calculates the speed of the input signal, step e) is the correction signal The difference signal before and after damage is calculated using the signal when there is no damage to the composite material, and the damage location of the composite material is determined using the time the difference signal before and after damage arrives and the speed of the input signal calculated in step d). Estimate. The schematic operations of steps c) and e) described above are the same as those described in the background art, and a detailed description will be given of the process of separating signals in steps c) and d).

c) 단계는 상술한 바와 같은 동작을 위해, c-1) 단계를 포함할 수 있다.Step c) may include step c-1) for the above-described operation.

c-1) 단계는 아래 수식 7을 이용해 상기 보정신호에 포함된 신호 중 신호의 에너지가 입력신호와 가장 유사한 신호를 분리하고, 분리된 신호를 제1신호로 판별한다.In step c-1), among the signals included in the correction signal, a signal whose energy is most similar to an input signal is separated using Equation 7 below, and the separated signal is determined as a first signal.

[수식 7][Equation 7]

Figure pat00050
Figure pat00050

(r(t)은 보정신호에 포함된 신호 중 입력신호와 에너지 크기가 가장 유사한 신호,

Figure pat00051
는 시간영역에서의 보정신호를 STFT한 신호)(r(t) is the signal that has the most similar energy to the input signal among the signals included in the correction signal,
Figure pat00051
Is the signal obtained by STFT of the correction signal in the time domain)

상술한 수식 7에서 입력신호와 에너지 크기가 가장 유사한 제1신호를 판별하는 이유는, 가진부(100)에서 센서로 가장 먼저 도달한 제1신호의 경우, 다른 경로(가진부-센서-센서 또는 가진부-손상부분-센서)로 도달한 신호들보다 진폭의 감쇠가 가장 적게 일어나기 때문이다. r(t)은 R을 인버스 STFT한 시간영역에서의 신호로, 수식 7에서는 Inverse FFT라고 기재되어 있지만, 이는 인버스 STFT와 같은 의미이며, 수식 7에서 사용되는 p값과

Figure pat00052
값은 복합재의 재질에 따라 달리 결정되는 값이다.The reason for determining the first signal having the most similar energy level with the input signal in Equation 7 above is that in the case of the first signal that first reaches the sensor from the excitation unit 100, a different path (excitation unit-sensor-sensor or This is because the attenuation of the amplitude occurs less than the signals reached by the excitation part-damaged part-sensor). r(t) is a signal in the time domain in which R is inverse STFT, and it is described as Inverse FFT in Equation 7, but it has the same meaning as the inverse STFT, and the p value used in Equation 7 and
Figure pat00052
The value is determined differently depending on the material of the composite material.

도 8은 가진부(100)에서 생성한 입력신호와, 센서에서 측정한 측정신호를 개략적으로 도시하고, 측정신호에서 제1신호를 분리한 것을 도시한 것이다.FIG. 8 schematically shows the input signal generated by the excitation unit 100 and the measurement signal measured by the sensor, and shows the separation of the first signal from the measurement signal.

도 8에 도시된 바와 같이, 센서에서 측정된 측정신호에는 가진부(100)에서 입력신호를 생성한 시점에 입력신호와 동일한 파형의 신호가 포함되어 있다. 이는 전자기 간섭 현상(EMI, Electro Magnetic Interference)로 발생한 잡음성분으로, c-1) 단계는 상기한 수식 7을 적용하기 이전에, 보정신호에서 입력신호를 뺀 신호를

Figure pat00053
로 하여 수식 7을 적용할 수 있다.As shown in FIG. 8, the measurement signal measured by the sensor includes a signal having the same waveform as the input signal at the time when the excitation unit 100 generates the input signal. This is a noise component generated by an electromagnetic interference (EMI, Electro Magnetic Interference). In step c-1), before applying Equation 7 above, the signal minus the input signal from the correction signal is
Figure pat00053
Equation 7 can be applied by

도 8에 도시된 바와 같이, 가진부(100)에서 입력신호가 가진된 이후 제1신호(S1)는 소정 시간이 지난 후, 가장 먼저 센서에서 측정되며, d) 단계는 제1신호(S1)가 측정되기까지의 시간과, 미리 알고 있는 가진부(100)와 특정 센서간의 거리를 통해, 입력신호가 전달되는 속도를 계산한다.As shown in FIG. 8, after the excitation unit 100 excites the input signal, the first signal S1 is first measured by the sensor after a predetermined time has passed, and the step d) is the first signal S1. The speed at which the input signal is transmitted is calculated through the time until is measured and the distance between the excitation unit 100 and a specific sensor known in advance.

도 9는 손상부분의 유무만 다르고, 다른 조건이 동일한 센서에서 측정된 각각의 신호와 신호의 차이를 도시한 것이다. 보다 구체적으로, 도 9a는 복합재에 손상부분이 없을 때, 센서에서 측정된 측정신호이고, 도 9b는 복합재에 손상부분이 있을 때 동일한 센서에서 측정된 신호이며, 도 9c는 도 9a와 도 9b에 도시된 신호의 차, 즉 손상 전후 차이 신호를 도시한 것이다.9 shows the difference between each signal and the signal measured by the sensor having different conditions different from only the presence or absence of a damaged part. More specifically, FIG. 9A is a measurement signal measured by the sensor when there is no damaged part in the composite material, FIG. 9B is a signal measured by the same sensor when there is a damaged part in the composite material, and FIG. 9C is shown in FIGS. 9A and 9B. It shows the difference of the illustrated signal, that is, the difference signal before and after damage.

도 9c에 도시된 신호는 복합재에 손상된 부분이 있음으로써 변형된 신호이므로, 입력신호가 가진부(100)와 손상부분을 통해 센서로 이동한 신호임을 알 수 있으므로, 도 9c에 도시된 신호가 도달한 시간을 이용해 가진부(100)-손상부분-센서의 거리를 알 수 있다.Since the signal shown in FIG. 9c is a signal that has been transformed by the presence of a damaged part in the composite material, it can be seen that the input signal is a signal that has moved to the sensor through the excitation part 100 and the damaged part, so that the signal shown in FIG. Using one hour, the distance between the excitation unit 100-the damaged part-the sensor can be known.

도 10은 도 9c에 도시된 손상 전후 차이 신호를 확대 도시한 후, 손상 전후 차이 신호를 다수개의 신호로 분리한 것이다.10 is an enlarged view of the difference signal before and after damage shown in FIG. 9C, and then the difference signal before and after damage is divided into a plurality of signals.

도 10에 도시된 바와 같이, 손상 전후 차이 신호는 다수개의 신호가 조합된 형태의 신호인데, 이는 손상부분(110)이 점이 아니라, 일정 영역에 걸쳐서 형성되기 때문이다. 따라서, 도 10에서 분리된 다수개의 신호 각각이 도달한 시간과 d) 단계에서 구한 입력신호의 속도를 이용하면, 도 10에 도시된 바와 같이 가진부(100)-손상부분(110)의 주변 중 일부위치-센서간의 거리를 구할 수 있고, 가진부(100)와 센서를 중심으로 하는 타원을 그려 손상부분(110)의 영역을 추정할 수 있다.As shown in FIG. 10, the difference signal before and after damage is a signal in the form of a combination of a plurality of signals, because the damaged part 110 is formed over a certain area, not a point. Therefore, using the time when each of the plurality of signals separated in FIG. 10 arrives and the speed of the input signal obtained in step d), as shown in FIG. 10, the area around the excitation unit 100-damaged part 110 It is possible to obtain the distance between the partial position-sensor, and to estimate the area of the damaged part 110 by drawing an ellipse centered on the excitation part 100 and the sensor.

e) 단계에서 가진부와 다수개의 센서를 이용하여 손상부분의 위치를 추정하는 방식은 가진부-제1센서, 가진부-제2센서, 가진부-제3센서 각각에서 추정된 가진부-손상부분-센서의 거리를 이용해, 가진부와 센서를 중심으로하는 타원을 그리고, 세 개의 타원이 겹쳐지는 부분을 손상부분으로 추정하는 방식으로, 배경기술에서 설명한 종래의 방식과 동일하다.The method of estimating the location of the damaged part using the excitation part and multiple sensors in step e) is the excitation part-damage estimated by each of the excitation part-first sensor, the excitation part-second sensor, and the excitation part-third sensor. By using the distance of the part-sensor, an ellipse centered on the excitation part and the sensor is drawn, and a part where the three ellipses overlap is estimated as a damaged part, which is the same as the conventional method described in the background art.

본 발명의 상기한 실시 예에 한정하여 기술적 사상을 해석해서는 안된다. 적용범위가 다양함은 물론이고, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당업자의 수준에서 다양한 변형 실시가 가능하다. 따라서 이러한 개량 및 변경은 당업자에게 자명한 것인 한 본 발명의 보호범위에 속하게 된다.The technical idea should not be interpreted as limited to the above-described embodiments of the present invention. As well as a variety of application ranges, various modifications may be made at the level of those skilled in the art without departing from the gist of the present invention claimed in the claims. Therefore, these improvements and changes will fall within the scope of protection of the present invention as long as it is apparent to those skilled in the art.

S1 : 제1신호
S2 : 제2신호
10 : 액추에이터
11 : 센서
12 : 손상위치
100 : 복합재
110 : 손상부분
200 : 가진부
310 : 제1센서
320 : 제2센서
330 : 제3센서
S1: first signal
S2: second signal
10: actuator
11: sensor
12: Damage location
100: composite
110: damaged part
200: excitation part
310: first sensor
320: second sensor
330: third sensor

Claims (10)

a) 복합재에 설치된 가진부에서 탄성파인 입력신호를 발생시켜 상기 복합재를 따라 전파되게 하고, 복합재에 설치된 센서에서 상기 복합재를 따라 전파되는 신호를 측정하는 단계;
b) 상기 센서에서 측정된 측정신호에서 상기 입력신호를 이용해 잡음신호의 크기를 추정한 후, 상기 측정신호에서 상기 잡음신호가 제거된 보정신호를 계산하는 단계;
c) 상기 보정신호에 포함된 신호들을 분리하고, 상기 보정신호에 포함된 신호들 중, 상기 가진부에서 상기 센서로 가장 먼저 도착한 제1신호를 판별하는 단계;
d) 상기 가진부와 상기 센서의 간격과 상기 제1신호가 도달하기까지의 시간을 이용해 상기 입력신호의 속도를 계산하는 단계; 및
e) 상기 보정신호와 상기 복합재에 손상이 없을 경우의 신호를 이용해 손상 전후 차이 신호를 계산하고, 상기 손상 전후 차이 신호가 도달한 시간과 상기 d) 단계에서 계산된 상기 입력신호의 속도를 이용해 상기 복합재의 손상위치를 추정하는 단계;
를 포함하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
a) generating an input signal that is an elastic wave from an excitation part installed in the composite material to propagate along the composite material, and measuring a signal propagating along the composite material by a sensor installed in the composite material;
b) estimating a magnitude of a noise signal using the input signal from the measurement signal measured by the sensor, and then calculating a correction signal from which the noise signal is removed from the measurement signal;
c) separating signals included in the correction signal, and determining a first signal that arrives first from the excitation unit to the sensor from among signals included in the correction signal;
d) calculating the speed of the input signal using the distance between the excitation unit and the sensor and the time until the first signal arrives; And
e) The difference signal before and after damage is calculated using the correction signal and the signal when the composite material is not damaged, and the difference signal before and after damage arrives and the speed of the input signal calculated in step d) Estimating the damage location of the composite material;
Damage detection method using the signal energy of the acoustic wave comprising a.
제1항에 있어서,
상기 b) 단계는,
b-1) 아래 수식 5를 이용해 상기 잡음신호의 크기를 추정하는 단계;
를 포함하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
[수식 5]
Figure pat00054

(
Figure pat00055
는 추정된 상기 잡음신호의 크기,
Figure pat00056
는 상기 측정신호의 에너지 크기,
Figure pat00057
는 상기 입력신호의 에너지 크기)
The method of claim 1,
Step b),
b-1) estimating the magnitude of the noise signal using Equation 5 below;
Damage detection method using the signal energy of the acoustic wave comprising a.
[Equation 5]
Figure pat00054

(
Figure pat00055
Is the estimated size of the noise signal,
Figure pat00056
Is the energy level of the measurement signal,
Figure pat00057
Is the energy level of the input signal)
제2항에 있어서,
상기 b) 단계는,
b-2) 상기 b-1) 단계 이후에 수행되며, 아래 수식 6을 이용해 시간영역에서의 상기 보정신호를 계산하는 단계;
를 더 포함하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
[수식 6]
Figure pat00058

(상기
Figure pat00059
는 상기 측정신호를 STFT(Short Time Fourier Transform)한 신호, 상기
Figure pat00060
는 주파수 영역에서의 상기 보정신호,
Figure pat00061
는 상기
Figure pat00062
을 inverse STFT한 시간영역에서의 보정신호)
The method of claim 2,
Step b),
b-2) performing the step b-1) and calculating the correction signal in the time domain using Equation 6 below;
Damage detection method using the signal energy of the acoustic wave, characterized in that it further comprises.
[Equation 6]
Figure pat00058

(remind
Figure pat00059
Is a signal obtained by short-time fourier transform (STFT) of the measurement signal,
Figure pat00060
Is the correction signal in the frequency domain,
Figure pat00061
Is reminded
Figure pat00062
Correction signal in time domain inverse STFT)
제3항에 있어서,
상기 b) 단계는,
b-3) 상기 b-2) 단계 이후에 수행되고, 상기 수식 5의
Figure pat00063
을 상기 수식 6을 통해 구해진
Figure pat00064
로 대체하여 상기 보정신호에 포함된 잡음신호의 크기를 추정하고, 상기 보정신호의 크기 대비 상기 보정신호에 포함된 잡음신호의 크기가 기준치 이상일 경우, 상기 수식 6의
Figure pat00065
을 상기
Figure pat00066
로 대체하여 잡음신호의 크기가 줄어든 보정신호를 생성하는 단계;
를 더 포함하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
The method of claim 3,
Step b),
b-3) is performed after step b-2), and in Equation 5
Figure pat00063
Obtained through Equation 6 above
Figure pat00064
In the case where the size of the noise signal included in the correction signal is estimated and the size of the noise signal included in the correction signal is greater than or equal to the reference value,
Figure pat00065
Remind
Figure pat00066
Generating a correction signal having a reduced noise signal by replacing it with
Damage detection method using the signal energy of the acoustic wave, characterized in that it further comprises.
제4항에 있어서,
상기 b-3) 단계는 생성된 보정신호 내에 포함된 잡음신호의 크기가 기준치 이하가 될 때 까지 반복 수행되는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
The method of claim 4,
Step b-3) is a damage detection method using signal energy of an acoustic wave, characterized in that it is repeatedly performed until the magnitude of the noise signal included in the generated correction signal becomes less than a reference value.
제1항에 있어서,
상기 c) 단계는
c-1) 아래 수식 7을 이용해 상기 보정신호에 포함된 신호 중 에너지가 가장 큰 신호를 분리해, 분리된 신호를 상기 제1신호로 판별하는 단계;
를 포함하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
[수식 7]
Figure pat00067

(r(t)은 보정신호에 포함된 신호 중 에너지가 가장 큰 신호,
Figure pat00068
는 시간영역에서의 보정신호를 STFT한 신호,
Figure pat00069
는 입력신호)
The method of claim 1,
Step c)
c-1) separating a signal having the largest energy among signals included in the correction signal using Equation 7 below, and determining the separated signal as the first signal;
Damage detection method using the signal energy of the acoustic wave comprising a.
[Equation 7]
Figure pat00067

(r(t) is the signal with the largest energy among the signals included in the correction signal,
Figure pat00068
Is the signal obtained by STFT of the correction signal in the time domain,
Figure pat00069
Is the input signal)
제1항에 있어서,
상기 e) 단계는, 아래 수식 7을 이용해 상기 손상 전후 차이 신호에 포함된 신호 중 하나를 분리하고, 분리된 신호가 도달한 시간과 상기 d) 단계에서 계산된 상기 입력신호의 속도를 이용해 상기 복합재의 손상위치를 추정하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
[수식 7]
Figure pat00070

(r(t)은 상기 손상 전후 차이 신호에 포함된 신호 중 에너지가 가장 큰 신호,
Figure pat00071
는 시간영역에서의 상기 손상 전후 차이 신호를 STFT한 신호,
Figure pat00072
는 입력신호)
The method of claim 1,
In the step e), one of the signals included in the difference signal before and after the damage is separated using Equation 7 below, and the composite material is performed using the time the separated signal arrives and the speed of the input signal calculated in step d). Damage detection method using the signal energy of the acoustic wave, characterized in that estimating the damage location of.
[Equation 7]
Figure pat00070

(r(t) is the signal with the largest energy among the signals included in the difference signal before and after damage,
Figure pat00071
Is a signal obtained by STFT of the difference signal before and after damage in the time domain,
Figure pat00072
Is the input signal)
제7항에 있어서,
상기 e) 단계는 반복 수행되어 상기 손상 전후 차이 신호에 포함된 소정 개수의 신호를 분리하고,
상기 e) 단계가 반복 수행될 때, 상기 수식 7에서 상기
Figure pat00073
신호는 이전에 수행된 e) 단계의
Figure pat00074
신호에서 이전에 수행된 e) 단계에서 분리된 신호를 뺀 신호인 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
The method of claim 7,
Step e) is repeatedly performed to separate a predetermined number of signals included in the difference signal before and after damage,
When step e) is repeatedly performed, in Equation 7
Figure pat00073
The signal was previously performed in step e).
Figure pat00074
A damage detection method using signal energy of an acoustic wave, characterized in that the signal is a signal obtained by subtracting the signal separated in step e) previously performed from the signal.
제8항에 있어서,
상기 e) 단계는 반복 수행되어 상기 손상 전후 차이 신호에 포함된 소정 개수의 신호를 분리하여, 손상위치의 영역을 추정하는 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
The method of claim 8,
The step e) is repeatedly performed to separate a predetermined number of signals included in the difference signal before and after the damage, and estimate the area of the damage location.
제1항에 있어서,
상기 복합재에 설치된 센서는 복수개인 것을 특징으로 하는 탄성파의 신호 에너지를 이용한 손상 탐지 방법.
The method of claim 1,
Damage detection method using signal energy of an acoustic wave, characterized in that a plurality of sensors installed in the composite material.
KR1020190081963A 2019-07-08 2019-07-08 Damage detection method using lamb wave signal energy KR102204747B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190081963A KR102204747B1 (en) 2019-07-08 2019-07-08 Damage detection method using lamb wave signal energy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190081963A KR102204747B1 (en) 2019-07-08 2019-07-08 Damage detection method using lamb wave signal energy

Publications (2)

Publication Number Publication Date
KR20210006107A true KR20210006107A (en) 2021-01-18
KR102204747B1 KR102204747B1 (en) 2021-01-19

Family

ID=74236729

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190081963A KR102204747B1 (en) 2019-07-08 2019-07-08 Damage detection method using lamb wave signal energy

Country Status (1)

Country Link
KR (1) KR102204747B1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113533516A (en) * 2021-07-12 2021-10-22 北京交通大学 Damage detection method in steel rail based on Lamb wave time field method
KR102382152B1 (en) * 2021-09-03 2022-04-04 주식회사 솔리드아이티 Image Based Damage Localization Estimation Method With Selection of signal having specific frequency and Definition of Damage from Signal for a Panel Structure in Structure Health Monitoring Based on Guided Elastic Wave
KR102391019B1 (en) * 2021-11-30 2022-04-27 주식회사 솔리드아이티 Image Based Damage Localization Estimation Method With Definition of Damage from Signal for a Curved Structure in Structure Health Monitoring Based on Guided Elastic Wave
CN115876883A (en) * 2022-12-29 2023-03-31 南京航空航天大学 Detection method and detection system for layered damage position of composite laminated plate

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20090012818A (en) 2007-07-31 2009-02-04 한국전력공사 Non-destruction apparatus
KR101157811B1 (en) * 2010-08-17 2012-06-22 한국표준과학연구원 Method for locating source of damage in large composite material structure
KR20160110722A (en) * 2015-03-11 2016-09-22 주식회사 하이드로넷 Method and System for detecting water leak location based on elastic wave velocity measured by section in pipe

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20090012818A (en) 2007-07-31 2009-02-04 한국전력공사 Non-destruction apparatus
KR101157811B1 (en) * 2010-08-17 2012-06-22 한국표준과학연구원 Method for locating source of damage in large composite material structure
KR20160110722A (en) * 2015-03-11 2016-09-22 주식회사 하이드로넷 Method and System for detecting water leak location based on elastic wave velocity measured by section in pipe

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113533516A (en) * 2021-07-12 2021-10-22 北京交通大学 Damage detection method in steel rail based on Lamb wave time field method
KR102382152B1 (en) * 2021-09-03 2022-04-04 주식회사 솔리드아이티 Image Based Damage Localization Estimation Method With Selection of signal having specific frequency and Definition of Damage from Signal for a Panel Structure in Structure Health Monitoring Based on Guided Elastic Wave
KR20230035202A (en) * 2021-09-03 2023-03-13 주식회사 솔리드아이티 Damage Localization Estimation Method and Damage Localization Estimation Device
KR102391019B1 (en) * 2021-11-30 2022-04-27 주식회사 솔리드아이티 Image Based Damage Localization Estimation Method With Definition of Damage from Signal for a Curved Structure in Structure Health Monitoring Based on Guided Elastic Wave
CN115876883A (en) * 2022-12-29 2023-03-31 南京航空航天大学 Detection method and detection system for layered damage position of composite laminated plate
CN115876883B (en) * 2022-12-29 2024-03-29 南京航空航天大学 Method and system for detecting layered damage position of composite material laminated plate

Also Published As

Publication number Publication date
KR102204747B1 (en) 2021-01-19

Similar Documents

Publication Publication Date Title
KR102204747B1 (en) Damage detection method using lamb wave signal energy
Prosser et al. Time-frequency analysis of the dispersion of Lamb modes
Nakamura et al. Mode conversion behavior of SH guided wave in a tapered plate
Edwards et al. Depth gauging of defects using low frequency wideband Rayleigh waves
US6360609B1 (en) Method and system for interpreting and utilizing multimode dispersive acoustic guided waves
JP2008185426A (en) Integrity evaluation support device, method, and program of concrete pile
CA2616900C (en) Method for error-free checking of tubes for surface faults
Ohtsu et al. Principles of the acoustic emission (AE) method and signal processing
Ma et al. The reflection of guided waves from simple dents in pipes
KR20090017769A (en) Method for non-destructive testing of concretestructure
JP2013130572A (en) Ultrasonic thickness measurement method and instrument
US11193912B2 (en) Ultrasonic measurement apparatus and ultrasonic measurement method
CN108426948A (en) A kind of electromagnet ultrasonic changer and its working method of the single mode Lamb wave of excitation
RU2480742C1 (en) Method for acoustic emission control
JP5562118B2 (en) Ultrasonic nondestructive measuring method, ultrasonic nondestructive measuring device, and program
JP4764921B2 (en) Ultrasonic exploration method using resonance phenomenon
JP5507279B2 (en) Ultrasonic inspection method and apparatus
KR20010038725A (en) Method for non-destructive testing of concrete structure
Draudvilienė et al. Validation of dispersion curve reconstruction techniques for the A0 and S0 modes of Lamb waves
JP4519852B2 (en) Ultrasonic exploration method and apparatus using resonance phenomenon
JP6606342B2 (en) Ultrasonic measuring device
US9518959B2 (en) Structural health monitoring system and method
JP2003149214A (en) Nondestructive inspecting method and its apparatus using ultrasonic sensor
CN208383812U (en) A kind of electromagnet ultrasonic changer exciting single mode Lamb wave
JP5268686B2 (en) Measuring apparatus and measuring method by electromagnetic ultrasonic method

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant