KR20100106120A - Laser ultrasonic measuring apparatus and laser ultrasonic measuring method - Google Patents
Laser ultrasonic measuring apparatus and laser ultrasonic measuring methodInfo
- Publication number
- KR20100106120A KR20100106120A KR1020090024585A KR20090024585A KR20100106120A KR 20100106120 A KR20100106120 A KR 20100106120A KR 1020090024585 A KR1020090024585 A KR 1020090024585A KR 20090024585 A KR20090024585 A KR 20090024585A KR 20100106120 A KR20100106120 A KR 20100106120A
- Authority
- KR
- South Korea
- Prior art keywords
- laser
- wave
- ultrasonic
- frequency
- phase transformation
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/008—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means by using ultrasonic waves
Abstract
Description
본 발명은 레이저를 피검사체에 조사하여 발생하는 초음파를 검출하는 레이저 초음파 측정장치 및 측정방법에 관한 것이며, 특히 재료의 상변태율을 계측하는 레이저 초음파 측정장치 및 측정방법에 관한 것이다.BACKGROUND OF THE
레이저 초음파법은 피검사체에 대해 펄스레이저광을 조사하여 발생하는 열적 응력 또는 표면 근방의 피검사체 자체의 기화에 의해 발생하는 기화 응력을 이용하여 초음파를 송신하고 연속 발진하는 다른 레이저광을 수신점에 조사하여, 그 직진성이나 가간섭성을 이용하여 초음파에 의해 유기되는 변위를 수신하여 피검사체 내부를 전파한 초음파를 검출하는 기술이다.The laser ultrasonography transmits an ultrasonic wave by using thermal stress generated by irradiating pulsed laser beam to an inspected object or vaporization stress generated by vaporization of the inspector itself near the surface, and transmits other laser beams to the receiving point continuously. It is a technique of detecting the ultrasonic wave propagated inside the subject under test by receiving the displacement induced by the ultrasonic wave by using the straightness or coherence.
레이저 초음파법은 초음파를 사용하여 재료의 균열이나 내재하는 결함의 검출 또는 재료특성의 평가를 비접촉으로 행하는 것이 가능하며, 여러 재료평가분야로의 응용이 기대되고 있다.The laser ultrasonic method is capable of performing non-contact detection of material cracks and inherent defects or evaluation of material properties using ultrasonic waves, and application to various material evaluation fields is expected.
도 13을 참조하여 레이저 초음파법의 원리에 대해 설명한다. With reference to FIG. 13, the principle of the laser ultrasonic method is demonstrated.
높은 에너지의 펄스레이저인 초음파 발생용 레이저를 예를 들어, 피검사체인 강재 표면에 조사하면 그 충격으로 금속 표면에 생기는 열팽창 및 수축에 의해 왜곡이 발생한다. 그리고, 발생한 왜곡이 초음파로 강재 내부를 전파한다. 다음으로, 초음파 검출용 단일 주파수의 연속 레이저광을 금속 표면에 조사하면 이 반사광은 전파한 초음파에 의한 표면 진동에 따른 주파수 변화(도플러 시프트)를 한다. 이하에 도플러 시프트양(Δf)을 나타내는 식을 나타낸다.When the laser for generating an ultrasonic wave, which is a high-energy pulse laser, is irradiated to the surface of a steel, for example, an object to be inspected, distortion occurs due to thermal expansion and contraction generated on the metal surface by the impact. The generated distortion propagates inside the steel by ultrasonic waves. Next, when a continuous laser beam of a single frequency for ultrasonic detection is irradiated onto the metal surface, the reflected light causes a frequency change (Doppler shift) according to the surface vibration by the propagated ultrasonic waves. The formula showing the Doppler shift amount Δf is shown below.
Δf = 2V/λΔf = 2V / λ
여기서, V = 표면변위속도, λ = 레이저파장Where V = surface displacement velocity and λ = laser wavelength
레이저 초음파법을 이용한 계측기는 파브리 페롯 간섭계 등의 레이저 간섭계를 구비하고 있다. 파브리페롯 간섭계는 특정주파수만을 진동시켜 투과시키는 필터로 동작한다. 예를 들면, 강재 내부에 결함부가 있는 경우 표면 진동이 통상의 강재와 다르기 때문에 도플러 시프트양은 통상적인 강재와 다른 값을 나타낸다. 이 때문에 파브리페롯 간섭계를 투과하는 투과 광량이 변화하여 검사재의 균열이나 결함의 검사 또는 재료평가를 행할 수 있다.The measuring instrument using the laser ultrasonic method is equipped with a laser interferometer, such as a Fabry-Perot interferometer. The Fabriferot interferometer acts as a filter that vibrates and transmits only a specific frequency. For example, when there are defects inside the steel, the Doppler shift amount shows a different value from the conventional steel because the surface vibration is different from that of the conventional steel. For this reason, the amount of transmitted light passing through the Fabry-Perot interferometer changes, and inspection or material evaluation of cracks and defects in the inspection material can be performed.
레이저 초음파법으로는 레이저광이 닿는 범위라면 검사대상에 아무것도 접촉시키지 않고 검사대상의 상태를 관측할 수 있다. 레이저광은 광섬유나 밀러를 이용함으로써 접촉이나 접근이 곤란한 피검사체에 대해서도 비교적 용이하게 초음파 검사를 행할 수 있다.In the laser ultrasound method, the state of the inspection object can be observed without bringing anything into contact with the inspection object as long as it is within the reach of the laser light. By using an optical fiber or a mirror, a laser beam can perform ultrasonic inspection relatively easily also to the to-be-tested object which is difficult to contact or access.
도 14를 참조하여 레이저광 조사에 의한 초음파의 발생 원리에 대해 설명한다. 도 14(a)는 레이저광을 피검사체에 조사하여 피검사체의 일부를 증발시키는 현상인 어블레이션에 의해 초음파를 발생시키는 경우이며, 도 14(b)는 레이저광을 피검사체에 조사하여 피검사체에 열탄성 변화를 발생시켜 초음파를 발생시키는 경우이다.The generation principle of the ultrasonic wave by laser light irradiation is demonstrated with reference to FIG. FIG. 14 (a) is a case where ultrasonic waves are generated by ablation, which is a phenomenon in which a laser beam is irradiated to a test object to evaporate a part of the test object, and FIG. 14 (b) is a test object irradiated with laser light to the test object. This is the case where ultrasonic waves are generated by generating a thermoelastic change.
어블레이션이 발생하는 경우에는 증발응력에 의해 탄성파를 발생시키기 때문에 재료에 조사흔이 발생한다. 한편, 열탄성 경우에는 레이저에 의한 급속가열로 인한 열팽창 및 수축에 따른 탄성파를 발생시키나 재료에 조사흔은 발생하지 않는다. 열탄성의 경우에는 재료에 손상을 주지 않지만, 어블레이션 조건에 비해 열탄성 경우는 플루언스 레벨이 2자리수 가깝게 낮기 때문에 발생하는 초음파 음압은 이 이상의 비율로 감소한다.In the case of ablation, irrigation stresses cause seismic waves, causing irradiation marks on the material. On the other hand, in the case of thermoelasticity, elastic waves are generated due to thermal expansion and contraction due to rapid heating by a laser, but no irradiation traces occur on the material. In the case of thermoelasticity, the material is not damaged, but in the case of thermoelasticity compared with the ablation conditions, the ultrasonic sound pressure generated by the low fluence level is nearly two orders of magnitude lower.
종래의 레이저 초음파법은 대상재에 손상을 주어 초음파를 발생시키는 어블레이션 영역에서 사용하는 방법이다. 이 방법은 발생초음파의 강도가 강하고 초음파의 종파나 횡파의 음속을 검출한다(하기, 비특허문헌 1~5). 또한, 검출한 종파나 횡파의 음속을 활용하여 상변태율을 계측에 활용하는 방법(하기, 특허문헌 1 및 2)이 제안되어 있다.Conventional laser ultrasonication is a method used in an ablation area that damages an object and generates ultrasonic waves. In this method, the intensity of the generated ultrasonic waves is strong and the sound velocity of the longitudinal waves or the transverse waves of the ultrasonic waves is detected (
[비특허문헌 1] M. Dubois, A. Moreau, M. Militzer, and J. F. Bussiere, “Laser-Ultrasonic Monitoring of Phase Transformations in Steels”, Scripta Materialia, Vol. 39, No. 6, p.735-741, 1998.[Non-Patent Document 1] M. Dubois, A. Moreau, M. Militzer, and J. F. Bussiere, “Laser-Ultrasonic Monitoring of Phase Transformations in Steels”, Scripta Materialia, Vol. 39, No. 6, p. 735-741, 1998.
[비특허문헌 2] M. Ericsson, E.Lindh-Ulmgren, D. Artymowicz and B. Hutchison, “’Laser-ultrasonics (LUS) for microstructure characterization”. Research report IM-2003-113, Swedish institute for metals research, 2003.[Non-Patent Document 2] M. Ericsson, E. Lindh-Ulmgren, D. Artymowicz and B. Hutchison, “'Laser-ultrasonics (LUS) for microstructure characterization”. Research report IM-2003-113, Swedish institute for metals research, 2003.
[비특허문헌 3] 「J.D.Achenbach, “Wave propagation in elastic solids”, Amsterdam, North-Holland, 1973」[Non-Patent Document 3] J.D. Achenbach, “Wave propagation in elastic solids”, Amsterdam, North-Holland, 1973
[비특허문헌 4] C.B. Scruby and L.E. Drain, “Laser Ultrasonics -Techniques and Applications” , ISBN0-7503-0050-7, Adam Hilger, p.289, 1990.[Non-Patent Document 4] C.B. Scruby and L.E. Drain, “Laser Ultrasonics-Technologies and Applications”, ISBN 0-7503-0050-7, Adam Hilger, p. 289, 1990.
[비특허문헌 5] J. Huang, Y. Nagata, S. Krishnaswamy, and J. D. Achenbach,“Laser-based ultrasonics for flaw detection”, Proc. of IEEE Ultrasonics Symp., p1205-1209, 1994.[Non-Patent Document 5] J. Huang, Y. Nagata, S. Krishnaswamy, and J. D. Achenbach, “Laser-based ultrasonics for flaw detection”, Proc. of IEEE Ultrasonics Symp., p 1205-1209, 1994.
[특허문헌 1] 특개평11-166918호 공보[Patent Document 1] Japanese Patent Laid-Open No. 11-166918
[특허문헌 2] 특개평11-166919호 공보[Patent Document 2] Japanese Patent Application Laid-Open No. 11-166919
대상재의 표면에 손상을 주지 않고 초음파를 발생시키는 열탄성영역에서는 초음파의 발생기구가 다르기 때문에 발생한 전파방향의 지향성이 크게 다르고 결과적으로 종파 및 횡파의 음속을 계측하는 것이 곤란했다. 이 때문에 상변태율을 구하기 위해서는 어블레이션 영역에서 검출한 종파나 횡파의 음속을 이용할 필요가 있었다.In the thermoelastic region in which ultrasonic waves are generated without damaging the surface of the target material, since the ultrasonic generating mechanisms are different, the directivity of the generated propagation direction is greatly different, and as a result, it is difficult to measure the sound velocity of longitudinal and transverse waves. For this reason, it was necessary to use the sound velocity of the longitudinal wave or the transverse wave detected in the ablation region in order to obtain the phase transformation rate.
상술한 문제점을 감안하여, 종파나 횡파의 음속이 아니라 대상재의 표면에 주는 손상을 경감할 수 있는 열탄성 영역에서 군속도 제로의 판파초음파 주파수를 계측함으로써 상변태율을 계측하는 재료의 상변태율의 계측장치 및 계측방법을 제공하는 것을 목적으로 한다.In view of the above-described problems, a device for measuring the phase transformation rate of a material for measuring the phase transformation rate by measuring the plate ultrasonic wave frequency of zero group velocity in the thermoelastic region where the damage to the surface of the target material can be reduced, not the sound velocity of longitudinal or transverse waves. And an object of the measurement method.
상기 과제를 해결하기 위해 제공되는 본 발명의 레이저 초음파 측정장치 및 레이저 초음파 측정방법은 하기 (1)~(14)에 기재된 바와 같다.The laser ultrasonic measuring apparatus and laser ultrasonic measuring method of the present invention provided to solve the above problems are as described in the following (1) to (14).
(1) 피검사체에 초음파 발생용 레이저광원에 의해 펄스레이저광을 조사하여 초음파를 발생시키고 초음파 검출용 레이저광원에 의해 연속파 레이저광을 이 피검사체에 조사하여 그 레이저 반사광을 간섭계에서 간섭시켜서 상기 초음파에 의한 도플러 시프트에 의한 광강도 변화를 검지해서 상기 피검사체 내부를 전파한 초음파를 검출하는 레이저 초음파 측정장치에 있어서,(1) Irradiating pulsed laser light to an inspected object with a laser light source for generating an ultrasonic wave to generate an ultrasonic wave, and irradiating a continuous wave laser light to the inspected object with an ultrasonic detection laser light source and interfering the laser reflected light with an interferometer. In the laser ultrasonic measuring apparatus for detecting the ultrasonic wave propagated inside the test object by detecting a change in light intensity due to the Doppler shift
상기 펄스레이저광을 상기 피검사체 상에 점상 스폿으로 조사하는 조사 광학계와,An irradiation optical system for irradiating the pulsed laser light on a test subject with a point spot;
상기 점상 스폿에 의해 발생하여 피검사체의 두께방향으로 전파한 군속도 제로의 판파초음파를 검출하도록 상기 연속파 레이저광을 피검사체 상에 조사하여 그 레이저 반사광을 상기 간섭계에 입사시켜 간섭시키는 검출 광학계와, A detection optical system for irradiating the continuous wave laser light onto the inspected object so as to detect the plate wave ultrasonic wave generated by the point spot and propagating in the thickness direction of the inspected object, and causing the laser reflected light to enter the interferometer and interfering therewith;
상기 간섭계를 투과한 레이저광을 수광하여 광강도 변화를 전기신호로 출력하는 광검출부와, A light detector for receiving a laser beam transmitted through the interferometer and outputting a change in light intensity as an electric signal;
상기 광검출부로부터 출력된 전기신호가 입력되어 피검사체의 상기 판파초음파의 주파수를 도출하는 신호처리부를 구비하고,An electric signal outputted from the photodetector is input to provide a signal processor for deriving a frequency of the plate wave ultrasonic wave of the inspected object;
상기 신호처리부는, The signal processing unit,
상기 조사 광학계에 의해 발생한 판파초음파에 의한 광강도 변화의 전기신호에 기초하여 주파수 해석을 하여 대칭판파 주파수를 산출하는 주파수 산출부와, A frequency calculator for calculating a symmetric wave frequency by performing frequency analysis on the basis of an electrical signal of light intensity change caused by the wave wave ultrasonic wave generated by the irradiation optical system;
상기 대칭판파 주파수와 미리 작성된 대칭판파 주파수와 피검사체 두께와의 곱과 상변태율의 관계를 나타내는 교정선으로부터 상기 피검사체의 상변태율을 산출하는 상변태율 산출부를 구비하는 것을 특징으로 하는 레이저 초음파 측정장치를 제공한다.And a phase transformation rate calculating unit for calculating a phase transformation rate of the inspected object from a calibration line indicating a relationship between the product of the symmetric wave frequency, the pre-made symmetric wave frequency, the thickness of the inspected object, and the phase transformation rate. To provide.
(2) 상기 신호처리부는 별도로 측정된 상기 피검사체의 상변태율이 입력되고상기 대칭판파 주파수와 피검사체 두께와의 곱과 상기 상변태율의 관계를 나타내는 교정선을 작성하는 교정선 작성부를 추가로 구비하는 것을 특징으로 하는 (1)에 기재된 레이저 초음파 측정장치를 제공한다.(2) The signal processing unit further includes a calibration line preparation unit for inputting a phase transformation rate of the test subject separately measured and preparing a calibration line indicating a product of the symmetric wave frequency and the thickness of the test subject and the phase transformation rate. The laser ultrasonic measuring apparatus described in (1) is provided.
(3) 상기 연속파 레이저광을 상기 펄스 레이저광의 조사 스폿영역 내에 조사하는 것을 특징으로 하는 (1) 또는 (2)에 기재된 레이저 초음파 측정장치를 제공한다.(3) The laser ultrasonic measuring apparatus according to (1) or (2), wherein the continuous wave laser light is irradiated into the irradiation spot region of the pulse laser light.
(4) 상기 주파수 해석에서는 복수 주기의 판파초음파를 검출하는 수록시간이 설정되는 것을 특징으로 하는 (1)~(3) 중 어느 한 항에 기재된 레이저 초음파 측정장치를 제공한다.(4) The frequency analysis provides the laser ultrasonic measuring apparatus according to any one of (1) to (3), wherein the recording time for detecting the plural cycle wave ultrasonic waves is set.
(5) 상기 대칭판파 주파수는 군속도 제로의 대칭 판파모드 주파수인 것을 특징으로 하는 (1)~(4) 중 어느 한 항에 기재된 레이저 초음파 측정장치를 제공한다.(5) Said symmetric plate wave frequency is the symmetric plate wave mode frequency of group speed zero, The laser ultrasonic measuring apparatus in any one of (1)-(4) characterized by the above-mentioned.
(6) 상기 교정선의 작성에 사용되는 상변태율은 현미경조직관찰 또는 열팽창률측정에 의해 측정되는 것을 특징으로 하는 (1)~(5) 중 어느 한 항에 기재된 레이저 초음파 측정장치를 제공한다.(6) The laser ultrasonic measuring apparatus according to any one of (1) to (5), wherein the phase transformation rate used for preparing the correction line is measured by microscopic tissue observation or thermal expansion coefficient measurement.
(7) 상기 교정선은 상기 상변태율이 0%일 때의 주파수와 상변태율이 100%일 때의 주파수로부터 작성되는 선형 교정선인 것을 특징으로 하는 (1)~(6) 중 어느 한 항에 기재된 레이저 초음파 측정장치를 제공한다.(7) The correction line is a linear correction line created from a frequency when the phase transformation rate is 0% and a frequency when the phase transformation rate is 100%. Provided is a laser ultrasonic measuring device.
(8) 피검사체에 초음파 발생용 레이저광원에 의해 펄스레이저광을 조사하여 초음파를 발생시키고 초음파 검출용 레이저광원에 의해 연속파 레이저광을 이 피검사체에 조사하여 그 레이저 반사광을 간섭계에서 간섭시켜서 상기 초음파에 의한 도플러 시프트에 의한 광강도 변화를 검지하여 상기 피검사체 내부를 전파한 초음파를 검출하는 레이저 초음파 측정장치를 이용한 레이저 초음파 측정방법에 있어서,(8) Irradiating pulsed laser light to an inspected object with a laser light source for generating an ultrasonic wave to generate ultrasonic waves, and irradiating the continuous wave laser light to the inspected object with an ultrasonic detection laser light source and interfering the laser reflected light with an interferometer. In the laser ultrasonic measuring method using a laser ultrasonic measuring device for detecting the change in the light intensity due to the Doppler shift by the ultrasonic wave to detect the ultrasonic wave propagated inside the object,
상기 펄스레이저광을 상기 피검사체 상에 점상 스폿으로 조사하는 레이저광조사공정과,A laser light irradiation step of irradiating the pulsed laser light on a test subject with a point spot;
상기 점상 스폿에 의해 발생하여 피검사체의 두께방향으로 전파한 군속도 제로의 판파초음파를 검출하도록 상기 연속파 레이저광을 피검사체 상에 조사하여 이 레이저 반사광을 상기 간섭계에 입사시켜서 간섭시키는 초음파 검출공정과,An ultrasonic detection step of irradiating the continuous wave laser light onto the inspected object so as to cause interference with the interferometer by detecting the continuous wave laser light generated by the point spot and propagating in the thickness direction of the inspected object in the plate wave ultrasonic wave;
상기 간섭계를 투과한 레이저광을 수광하여 광강도 변화를 전기신호로 출력하는 광검출공정과,A photodetection step of receiving a laser beam transmitted through the interferometer and outputting a change in light intensity as an electric signal;
상기 전기신호에 기초하여 주파수 해석을 하여 대칭판파 주파수를 산출하는 주파수 산출공정과,A frequency calculating step of calculating a symmetric wave frequency by performing frequency analysis on the basis of the electrical signal;
상기 대칭판파 주파수와 미리 작성된 대칭판파 주파수와 피검사체 두께와의 곱과 상변태율의 관계를 나타내는 교정선으로부터 상기 피검사체의 상변태율을 산출하는 상변태율 산출부공정을 가지는 것을 특징으로 하는 레이저 초음파 측정방법을 제공한다.And a phase transformation rate calculating part step of calculating a phase transformation rate of the subject under test from a calibration line indicating a relationship between the product of the symmetric wave frequency, the pre-created symmetric wave frequency and the thickness of the subject, and the phase transformation rate. Provide a method.
(9) 별도로 측정된 상기 피검사체의 상변태율이 입력되고 상기 대칭판파 주파수와 피검사체 두께와의 곱과 상기 상변태율의 관계를 나타내는 교정선을 작성하는 교정선 작성공정을 추가로 가지는 것을 특징으로 하는 (8)에 기재된 레이저 초음파 측정방법을 제공한다.(9) further comprising a calibration line preparation process for inputting a phase transformation rate of the test object separately measured and preparing a calibration line indicating a relationship between the product of the symmetric wave frequency and the thickness of the test object and the phase transformation rate; The laser ultrasonic measuring method described in (8) is provided.
(10) 초음파 검출공정을 상기 펄스 레이저광의 조사 스폿영역 내에 조사하는 것을 특징으로 하는 (8) 또는 (9)에 기재된 레이저 초음파 측정방법을 제공한다.(10) The laser ultrasonic measuring method according to (8) or (9), wherein the ultrasonic detecting step is irradiated into the irradiation spot region of the pulsed laser light.
(11) 상기 주파수 해석에서는 복수 주기의 판파초음파를 검출하는 수록시간 이 설정되는 것을 특징으로 하는 (8)~(10) 중 어느 한 항에 기재된 레이저 초음파 측정방법을 제공한다.(11) The frequency analysis provides the laser ultrasonic measuring method according to any one of (8) to (10), wherein a recording time for detecting a plurality of cycle wave ultrasonic waves is set.
(12) 상기 대칭판파 주파수는 군속도 제로의 대칭 판파모드 주파수인 것을 특징으로 하는 (8)~(11) 중 어느 한 항에 기재된 레이저 초음파 측정방법을 제공한다.(12) The symmetric plate wave frequency is a symmetric plate wave mode frequency of group velocity zero, the laser ultrasonic measurement method according to any one of items (8) to (11).
(13) 상기 교정선 작성에 사용되는 상변태율은 현미경조직관찰 또는 열팽창률측정에 의해 측정되는 것을 특징으로 하는 (8)~(12) 중 어느 한 항에 기재된 레이저 초음파 측정방법을 제공한다.(13) The method for measuring the ultrasonic wave according to any one of (8) to (12), wherein the phase transformation rate used for preparing the correction line is measured by microscopic tissue observation or thermal expansion coefficient measurement.
(14) 상기 교정선은 상기 상변태율이 0%일 때의 주파수와 상변태율이 100%일 때의 주파수로부터 작성되는 선형 교정선인 것을 특징으로 하는 (8)~(13) 중 어느 한 항에 기재된 레이저 초음파 측정방법을 제공한다.(14) The correction line described in any one of (8) to (13), wherein the correction line is a linear correction line created from a frequency when the phase transformation rate is 0% and a frequency when the phase transformation rate is 100%. It provides a laser ultrasonic measuring method.
상기 레이저 초음파 측정장치 및 레이저 초음파 측정방법은 피검사체의 표면으로의 손상이 없는 펄스에너지밀도가 낮은 레이저광에 의해 군속도 제로의 판파초음파의 주파수측정을 높은 정밀도로 행할 수 있다. 그 결과, 보다 정밀도가 높은 검사재의 균열이나 결함의 검사 또는 재료평가를 비파괴, 비접촉으로 행할 수 있다.In the laser ultrasonic measuring apparatus and the laser ultrasonic measuring method, frequency measurement of a plate wave ultrasonic wave with a group velocity of zero can be performed with high precision by a laser beam having a low pulse energy density without damage to the surface of an inspected object. As a result, inspection or material evaluation of cracks and defects of inspection materials with higher precision can be performed non-destructively and non-contact.
이하, 첨부된 도면을 참조하여 본 발명의 실시형태를 설명한다. 그러나, 본 발명의 실시형태는 여러 가지의 다른 형태로 변형될 수 있으며, 본 발명의 범위가 이하 설명하는 실시형태로만 한정되는 것은 아니다. 도면에서의 요소들의 형상 및 크기 등은 보다 명확한 설명을 위해 과장될 수 있으며, 도면상의 동일한 부호로 표시되는 요소는 동일한 요소이다.Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. However, embodiments of the present invention may be modified in various other forms, and the scope of the present invention is not limited to the embodiments described below. The shape and the size of the elements in the drawings may be exaggerated for clarity and the same elements are denoted by the same reference numerals in the drawings.
상기 비특허문헌 3에 기재된 판파의 이론에서는, 판파에는 각각이 복수의 모드를 갖는 대칭판파와 비대칭판파가 있다. 여기서, 판파는 초음파 탐상 등의 초음파 이용분야에서 주지인 초음파 전파 현상이며, 종파와 횡파가 일체가 되어 전달됨으로써 판 전체를 진동시키는 것처럼 보이는 파이다. 비대칭판파는 판 두께의 중심부터 비대칭한 변위를 발생시키면서 전달하는 판파이고, 대칭판파는 판 두께의 중심부터 대칭한 변위를 발생시키면서 전달하는 판파이다.In the theory of the plate wave described in the said
도 1(a)는 피검사체(5)의 측면에서 본 레이저 초음파 측정장치(10)의 계략 구성도이다. 도 1(b)는 피검사체(5)를 레이저 초음파 측정장치(10)에서 본 도면이다.FIG. 1A is a schematic configuration diagram of the laser
도 1(a)에 나타낸 레이저 초음파 측정장치(10)는 초음파 발생용 레이저광원(11), 초음파 검출용 레이저광원(12), 레이저 간섭계(13), 신호처리부(60)를 갖는다. 도 1(b)에는 발생용 레이저광을 조사하는 초음파 발생용 레이저 스폿영역(6), 검출용 레이저광을 조사하는 초음파 검출용 레이저 스폿영역(7)이 나타나 있다.The laser
초음파 발생용 레이저광원(11)은 광학계를 통해 펄스레이저광을 초음파 발생 용 레이저 스폿영역(6)에 조사한다. 초음파 검출용 레이저광원(12)은 연속파인 레이저광을 광학계를 통해 초음파 발생용 레이저 스폿영역(6)보다 작은 초음파 검출용 레이저 스폿영역(7)에 조사하면 초음파 발생용 레이저 스폿영역(6)에 생긴 판파의 진동에 의해 검출용 레이저광의 파장에 도플러 시프트가 생긴다. 그리고, 도플러 시프트가 생긴 검출용 레이저의 레이저 반사광은 광학계(미도시)에 의해 페브리페롯 간섭계 등의 레이저 간섭계(13)에 안내되어, 해당 레이저 간섭계(13) 안에서 간섭된다. 이때, 레이저 간섭계(13)의 공진조건은 검출용 레이저의 원래 주파수에 대응한 공진조건에서 약간 벗어나게 한다. 일시적으로 도플러 시프트된 검출용 레이저 반사광을 레이저 간섭계(13)를 투과시키면 투과한 레이저광에서는 도플러 시프트에 의한 파장변화가 광강도의 변화로 변환된다.The ultrasonic wave generation
그리고, 투과된 레이저광을 해당 레이저광의 파장에 양호한 감도를 갖는 공지된 광검출기로 구성하는 광검출부(미도시)에서 검지하여 전기신호로 변환함으로써, 판파에 의한 진동파형(이하, 판파의 파형이라고도 함)을 전기신호(전압파형)으로 검출할 수 있다. 이러한 레이저 초음파법에서의 초음파 검출에 대해서는 예를 들면, 상기 비특허문헌 5에 기재되어 있다. 신호처리부(60)는 이 전기신호를 이하에서 설명 하는 바와 같이 신호처리를 함으로써 원하는 측정정보를 도출할 수 있다.Then, the transmitted laser light is detected by a photodetector (not shown) constituted by a known photodetector having a good sensitivity to the wavelength of the laser light and converted into an electrical signal. Can be detected as an electrical signal (voltage waveform). Ultrasonic detection in such a laser ultrasonic method is described in the
신호처리부(60)는 검출한 판파의 파형의 전기신호를 A/D 변환기로 디지털신호로 변환한 후, 예를 들면, 주파수 해석으로 FFT(Fast Fourier Transform)를 행하여 각 주파수 성분의 크기의 분포를 도출한다. 해당 분포에서 피크 주파수를 구해 피검사체에 발생한 군속도 제로의 판파 주파수를 검출한다.The
도 2 및 도 3을 참조하면, 레이저 초음파 측정장치(10)가 군속도 제로의 판파를 검출이 가능한 것을 나타낸다.2 and 3, it is shown that the laser
도 2는 레이저 초음파 측정장치(10)에 의해 펄스레이저조사에 의해 판파를 발생시키고 이를 검출했을 때의 전기신호의 파형을 나타낸 도면이다. 도 2(a)는 펄스레이저광을 조사했을 때의 검출한 시간영역의 검출파형을 나타낸다. 도 2(b)는 도 2(a)에 나타낸 시간파형을 푸리에 변환한 주파수 특성을 나타낸다.FIG. 2 is a diagram showing waveforms of electric signals when a wave wave is generated and detected by pulsed laser irradiation by the laser
도 2(a)의 시간(T1)은 초음파 발생용 펄스레이저를 조사한 시간이다. 도 2(b)에 나타낸 주파수 파형은 도 2(a)의 펄스레이저 조사 후부터 90㎲ 사이에 검출된 시간파형에 대해 FFT를 실행하여 얻은 것이다. 도 2(b)에 나타낸 주파수 스펙트럼 상의 피크(F1)는 FFT에 의해 얻어지는 가장 큰 주파수 성분이다.The time T1 of FIG. 2 (a) is a time when the ultrasonic wave generation laser beam is irradiated. The frequency waveform shown in FIG. 2 (b) is obtained by performing FFT on the time waveform detected between 90 kHz after the pulse laser irradiation of FIG. The peak F1 on the frequency spectrum shown in Fig. 2B is the largest frequency component obtained by the FFT.
해당 측정의 측정조건은 다음과 같다.The measurement conditions of the measurement are as follows.
피검사체 : 두께(d) 2㎜, SUS304Subject under test: Thickness (d) 2mm, SUS304
측정시 온도 : 실온Temperature during measurement: room temperature
FFT를 실행한 시간영역 : 90㎲Time domain with FFT: 90 ms
초음파 발생용 레이저광원 : YAG 레이저(최대 펄스에너지 450mJ/pulse(ND필터에서 감쇠), 펄스간격 10nsUltrasonic generation laser light source: YAG laser (maximum pulse energy 450mJ / pulse (damped by ND filter), pulse interval 10ns
레이저펄스 에너지 밀도 : 1.8mJ/㎟(열탄성영역)Laser pulse energy density: 1.8 mJ / mm2 (thermoelastic region)
초음파 검출용 레이저광원 : 파워 500㎽, 파장 532㎚Ultrasonic Detection Laser Light Source:
초음파 발생용 레이저빔 스폿 : 직경 5㎜Laser Beam Spot for Ultrasonic Generation:
검출 레이저빔 스폿 : 직경 1㎜Detection laser beam spot:
상기 비특허문헌 3에는 판파의 위상속도를 계산하기 위한 계산식이 기재되어있다. 대칭모드의 판파(Sn)(n은 정수)는 하기에 나타내는 수학식 1에 의해 「속도」와 「진동수×두께」의 관계로 나타난다. 또한, 군속도(Cg)는 수학식 2에 의해 구해진다.
도 3은 논리계산에 의한 군속도 제로의 S1 모드의 주파수를 나타낸 도면이다. 도 3에는 수학식 1에 나타낸 위상속도(c)와 주파수(f = ω/2π)의 관계 및 수학식 2에 의해 구해지는 군속도(Cg)와 주파수(f = ω/2π)와의 관계가 나타난다.3 is a diagram showing the frequency of the S1 mode of the group speed zero by logical calculation. 3 shows the relationship between the phase speed c and the frequency f = ω / 2π shown in
도 3에 나타낸 논리계산에 의한 S1 모드의 군속도 제로에서의 주파수는 검출된 주파수 1.3MHz(F1)와 일치한다. 이에 의해, 검출된 주파수(F1)는 군속도 제로의 S1 모드 판파 주파수인 것을 알 수 있다. 따라서, 도 2 및 도 3의 결과로 실온에서 열탄성영역(표면손상이 없는 조사조건)에서의 군속도 제로의 S1 모드 판파 주파수가 계측이 가능하다는 것을 알 수 있다.The frequency at group speed zero in S1 mode by the logic calculation shown in FIG. 3 coincides with the detected frequency 1.3 MHz (F1). Thus, it can be seen that the detected frequency F1 is the S1 mode fan wave frequency of group speed zero. Thus, it can be seen from the results of Figs. 2 and 3 that the S1 mode plate wave frequency at zero group velocity in the thermoelastic region (irradiation condition without surface damage) can be measured at room temperature.
또한, 피검사체(5)인 SUS304에서는 수학식 1에 사용되는 종파음속(VL) 및 횡파음속(VS)은 기지이며, 판파의 종파음속은 5,900[m/s]이고, 판파의 횡파음속은 3,200[m/s]이며, 도 3의 논리계산에서는 이 값을 이용하고 있다. 또한, 피검사체(5)인 SUS304는 두께가 2mm이고, 도 3의 논리계산에서는 이값을 이용하고 있다.In addition, in SUS304, which is the
군속도 제로의 판파는 발생개소에 멈추어 판 두께 방향에서 진동을 계속하기 때문에, 복수 주기를 검출할 수 있도록 예를 들면, 도 2(b)와 같이 90㎲의 장시간의 파형에 대해 FFT를 행함으로써 그 주파수를 계측할 수 있는 것을 알 수 있다.Since the plate wave with the group velocity of zero stops at the generation point and continues to oscillate in the plate thickness direction, the FFT is performed for a long time waveform of 90 Hz as shown in FIG. It can be seen that the frequency can be measured.
도 4는 레이저 초음파 측정장치(10)에 의해 측정 온도조건을 바꿔 S1 판파 모드의 판파를 발생시키고 이를 검출했을 때의 전기신호 파형의 FFT결과를 나타낸 도이다. 해당 측정 예에서의 측정조건은 이하에 나타내는 조건을 제외하고는 도 2를 참조하여 설명한 측정조건과 동일하다.FIG. 4 is a diagram showing an FFT result of an electric signal waveform when generating a wave wave in an S1 wave wave mode by changing the measurement temperature condition by the laser
측정 온도조건 : 실온에서 900℃까지 승온(1℃/1s)하여 5분 동안 유지하고, 500℃까지 강온(-20℃/분)한다. 측정은 900℃에서 샘플 공랭 중에 계측하고 온도는 샘플표면에 붙인 열전대의 계측값이다.Measurement temperature condition: It heated up to 900 degreeC from room temperature (1 degreeC / 1s), hold | maintains for 5 minutes, and it temperature-falls to 500 degreeC (-20 degreeC / min). Measurements are taken during sample air cooling at 900 ° C, and temperature is the measured value of the thermocouple attached to the sample surface.
피검사체 : 고탄소강(C = 0.85%), 두께 2mmSubject: High carbon steel (C = 0.85%), thickness 2mm
FFT를 실행한 시간영역 : 초음파 발생용 펄스레이저의 조사 타이밍 이후의 200㎲Time domain with FFT: 200 ms after irradiation timing of pulse generator for ultrasonic generation
F2, F3 및 F4에서 나타내는 피크는, 온도 529℃, 609℃, 711℃일 때 각각 검출한 피크 주파수이다. 이들 값은 고온시의 군속도 제로의 S1 모드 판파 주파수이다. F2~F4에서 나타내는 바와 같이 측정온도와 함께 군속도 제로의 S1 모드 판파 주파수가 변화하는 것을 알 수 있다.The peaks shown by F2, F3, and F4 are the peak frequencies detected when the temperatures are 529 ° C, 609 ° C, and 711 ° C, respectively. These values are the S1 mode wave frequency at group speed zero at high temperatures. As shown in F2 to F4, it can be seen that the S1 mode plate wave frequency of the group speed zero changes with the measured temperature.
도 5는 피검사체의 온도와 계측한 군속도 제로의 S1 모드 판파 주파수와의 관계를 나타낸 도면이다. △로 나타낸 검출 값은 실온에서 900℃까지 승온했을 때 측정한 군속도 제로의 S1 모드 판파 주파수이다. ○로 나타낸 검출 값은 900℃에서 500℃까지 강온했을 때 측정한 군속도 제로의 S1 모드 판파 주파수이다.It is a figure which shows the relationship between the temperature of a test subject and the S1 mode plate | board frequency of the measured group speed zero. The detection value shown by (triangle | delta) is S1 mode pan wave frequency of group speed zero measured when heated up to 900 degreeC from room temperature. The detection value shown by (circle) is S1 mode fan wave frequency of group speed zero measured when it cooled from 900 degreeC to 500 degreeC.
예를 들면, γ철에서 α철로 변태시킬 때의 비율을 강재의 상변태율이라 한다. 또한, 상변태율은 전자현미경으로 촬영한 화상에 대해 화상처리를 행함으로써 상변태율을 측정하는 방법, 또는 열팽창을 이용하여 상변태율을 산출하는 방법이 있다.For example, the ratio at which iron is transformed from γ iron to α iron is called the phase transformation rate of steel. As the image transformation rate, there is a method of measuring the image transformation rate by performing image processing on an image photographed with an electron microscope, or a method of calculating the image transformation rate using thermal expansion.
전자현미경을 이용하여 화상에서 상변태를 관측하는 것은, 인간이 화상을 보고 변태율을 판단해도 되고, 미세한 구조조직을 가지는 특징이 있으므로 화상인식 에 의해 상변태율을 구해도 된다. 예를 들면, SUS인 경우, Ti2O3 등의 탄화물을 이용하여 γ→α 변태시에 페라이트 조직이 미세화되기 때문에 이 페라이트의 미세화를 화상인식하여 변태율을 구할 수 있다.Observing the image transformation in the image using an electron microscope may be performed by the human being to determine the transformation rate by looking at the image, or may have a characteristic structure with fine structure, so the image transformation rate may be determined by image recognition. For example, in the case of SUS, the ferrite structure becomes fine at the time of?-? Transformation using carbides such as Ti 2 O 3 , so that the transformation rate can be obtained by image recognition of the miniaturization of the ferrite.
열팽창을 이용하여 상변태율을 구하는 방법은, 화상처리 측정이나 레이저 측정 등의 방법이 있다. 예를 들면, 화상처리 측정으로는 피검사체의 화상에 대해 화상처리로 측정개소의 변위를 측정함으로써 열팽창률을 취득한다.As a method of obtaining the phase transformation rate using thermal expansion, there are methods such as image processing measurement and laser measurement. For example, in the image processing measurement, the thermal expansion coefficient is obtained by measuring the displacement of the measurement location with the image processing on the image of the inspected object.
도 5에 나타낸 S1f100은 상변태율이 100%일 때의 측정 주파수이고, 도 5에 나타낸 S1f0은 상변태율이 0%일 때의 측정 주파수이다.S1f100 shown in FIG. 5 is a measurement frequency when a phase transformation rate is 100%, and S1f0 shown in FIG. 5 is a measurement frequency when a phase transformation rate is 0%.
온도를 변화시킴으로써 피검사체인 강재의 상변태율이 변화한다. 또한, 이 상변태율의 변화는 승온한 경우과 강온한 경우의 경향이 다른 것을 알 수 있다.By changing the temperature, the phase transformation rate of the steel to be inspected changes. In addition, it turns out that the change of this phase transformation rate differs in the tendency at the time of temperature rising from the case of temperature rising.
또한, 군속도 제로의 S1 모드 판파 주파수는 비특허문헌 6에 의해 수학식 3과 같이 되는 것이 나타나 있다. 한편, 종파음속(VL)과 상변태율은 특허문헌 1, 2 등에 나타나 있는 바와 같이 1 : 1 대응이 있다. 따라서, 군속도 제로의 S1 모드 판파 주파수와 피검사체 두께의 곱의 값이 상변태율의 정보를 갖는 것을 알 수 있다.In addition, the
[비특허문헌 6] Dominique Clorennec, etc, ‘Local and non-contact measurements of bulk acoustic wave velocities in thin isotropic plates and shells using zero group velocity Lamb modes,’ Journal of Applied Physics, 101, 034908, 2007.[Non-Patent Document 6] Dominique Clorennec, etc, 'Local and non-contact measurements of bulk acoustic wave velocities in thin isotropic plates and shells using zero group velocity Lamb modes,' Journal of Applied Physics, 101, 034908, 2007.
여기서, S1f는 군속도 제로의 S1 모드 판파 주파수이고, 은 정수이다.Where S1f is the S1 mode wave frequency at group velocity zero, Is an integer.
또한, S1f100×d0과 S1f0×d0 사이의 상변태율에서 상기와 같은 상변태율의 관측 또는 산출에 의해 샘플 데이터를 복수 개 취함으로써, 도 6(a) 또는 도 6(b)에 나타낸 바와 같이 교정선을 구할 수 있다. 다만, dO는 교정선을 구하기 위해 사용한 피검사체의 두께이다. 도 6(a)는 상변태율과 군속도 제로의 S1 모드 판파 주파수와의 관계를 선형의 교정선으로 나타낸 도면이며, 도 6(b)는 상변태율과 군속도 제로의 S1 모드 판파 주파수와의 관계를 비선형의 교정선으로 나타낸 도면이다. 또한, S1f100×d0과 S1f0×d0 사이의 상변태율에서 샘플 데이터를 복수 개 취하여 교정선을 얻은 결과, 상변태율과 군속도 제로의 S1 모드 판파 주파수와의 관계가 계측 정밀도의 관점에서 선형으로 간주해도 되는 경우는 도 6(a)의 선형의 교정선을 이용하고, 선형으로 간주할 수 없는 경우는 도 6(b)의 비선형의 교정선을 이용한다.Further, as shown in Fig. 6 (a) or 6 (b), a plurality of pieces of sample data are taken by observation or calculation of the above-described phase transformation rate at the phase transformation rate between S1f100xd0 and S1f0xd0. Can be obtained. However, dO is the thickness of the subject used to obtain the calibration line. Fig. 6 (a) is a diagram showing the relationship between the phase transformation rate and the S1 mode wave frequency at group speed zero with a linear calibration line, and Fig. 6 (b) is a nonlinear relationship between the phase transformation rate and S1 mode wave frequency at group speed zero. It is a figure shown with the correction line of. In addition, as a result of obtaining a calibration line by taking a plurality of sample data at a phase transformation rate between S1f100xd0 and S1f0xd0, the relationship between the phase transformation rate and the S1 mode wave frequency of group speed zero may be regarded as linear from the viewpoint of measurement accuracy. In this case, the linear calibration line of FIG. 6 (a) is used, and when it cannot be regarded as linear, the nonlinear calibration line of FIG. 6 (b) is used.
예를 들면, 도 6(a)에 나타낸 바와 같이 선형 교정선을 이용하는 경우, 임의의 피검사체(두께(d))의 측정 주파수 S1f의 상변태율(P)은 S1f100×d0 및 S1f0× d0를 이용하여 하기 수학식 4에 의해 구할 수 있다.For example, as shown in Fig. 6A, when the linear calibration line is used, the phase transformation rate P of the measurement frequency S1f of an arbitrary subject (thickness d) is used by S1f100xd0 and S1f0xd0. It can be obtained by the following equation (4).
또한, 도 6(b)에 나타낸 비선형 교정선을 이용하는 경우에는 측정된 주파수와 피검사체 두께와의 곱과 상변태율과의 관계를 다항식 등의 함수로부터 핏팅하거나 선형조사함으로써 하기 수학식 5에 나타낸 바와 같은 관계를 도출한다.In addition, in the case of using the nonlinear calibration line shown in Fig. 6 (b), the relationship between the product of the measured frequency, the thickness of the subject, and the phase transformation rate is obtained by fitting or linearly irradiating from a function such as a polynomial equation. Draw the same relationship.
여기서, F(S1f×d)는 S1f×d의 함수인 것을 나타낸다. 상변태율의 측정을 행하는 경우에는 군속도 제로의 S1 모드 판파의 측정 주파수 S1f를 측정하여 피검사체 두께와 함께 이 함수(F)에 대입함으로써 상변태율을 구할 수 있다.Here, F (S1fxd) shows that it is a function of S1fxd. In the case of measuring the phase transformation rate, the phase transformation rate can be obtained by measuring the measurement frequency S1f of the S1 mode plate wave at the group speed zero and substituting this function (F) together with the thickness of the subject.
이와 같이 상변태율과 군속도 제로의 S1 모드 판파 주파수×두께의 교정선을 구하여 상변태율을 수학식에 의해 산출할 수 있도록 함으로써, 군속도 제로의 S1 모드 판파 주파수(S1f)를 계측하여 S1f의 주파수와 두께로부터 상변태율을 산출할 수 있다.Thus, by calculating the phase transformation rate by the equation by calculating the correction line of the phase transformation rate and the group speed zero S1 mode wave frequency × thickness, the S1 mode wave frequency (S1f) of the group speed zero is measured to measure the frequency and thickness of S1f. The phase transformation rate can be calculated from
도 7을 참조하여, 피검사체를 이동하면서 레이저 초음파 측정방법을 행하는 예에 대해 설명한다. 피검사체(5)는 콘베이어 등의 반송설비에 반송된다. 측정조건은 하기에 나타낸 측정조건을 제외하고는 도 2를 참조하여 설명한 측정조건과 동일한다.With reference to FIG. 7, the example which performs a laser ultrasonic measuring method while moving a test subject is demonstrated. The inspected
반송속도 : V(㎜/s)Carrying Speed: V (mm / s)
S1f 계측을 위한 FFT 실행용 파형수록시간 : 200㎲Waveform recording time for FFT execution for S1f measurement: 200 ms
또한, FFT 실행용 파형수록시간은 짧으면 주파수 분해능이 저하하고 길면 노이즈성분이 지배적이 되어 SN이 저하한다. 피검사체가 이동함으로써 초음파 발생용 레이저 스폿영역(7a)을 초음파 검출용 레이저 스폿영역(6a)이 횡단하듯이 움직이기 때문에, FFT 실행용 파형수록시간은 도 2를 참조하여 설명한 측정조건에서의 파형수록시간보다 길게 100㎲~200㎲ 정도가 바람직하다.In addition, if the waveform recording time for FFT execution is short, the frequency resolution decreases, and if it is long, the noise component becomes dominant, and the SN decreases. Since the object under test moves as the ultrasonic spot laser spot region 7a for ultrasonic generation crosses the ultrasonic spot laser region 6a, the waveform recording time for the FFT execution is the waveform under the measurement conditions described with reference to FIG. It is preferable that the length is about 100 ms to about 200 ms longer than the recording time.
피검사체(5)는 반송설비에 반송되어 있기 때문에 초음파 검출용 레이저 스폿영역(6a)은 도시된 바와 같이 직사각형이 된다. 이때, 초음파 발생용 레이저 스폿영역(7a)은 이 스폿영역 내에 초음파 검출용 레이저 스폿영역(6a)을 넣어 두기 위해 파형수록시간 안에 초음파 검출용 레이저 스폿영역이 이동하는 길이(X)보다 길이 방향으로 긴 것이 바람직하다.Since the inspected
X = V×T = V(㎜/s)×FFT 실행을 위한 파형수록시간(s)X = V × T = V (mm / s) × Waveform recording time for executing FFT (s)
초음파 발생용 레이저의 조사 스폿영역의 길이방향의 길이는 X 이상으로 하는 것이 바람직하다. 또한, 초음파 발생용 레이저 조사시의 초음파 검출용 레이저스폿영역(6a)의 위치는 초음파 발생용 레이저 스폿영역(7a)에서의 피검사체이동의 하류측에 설정할 필요가 있다.It is preferable that the length of the longitudinal direction of the irradiation spot area | region of an ultrasonic wave generation laser shall be more than X. In addition, it is necessary to set the position of the ultrasonic detection laser spot area | region 6a at the time of the ultrasonic generation laser irradiation downstream of the test subject movement in the ultrasonic generation laser spot area | region 7a.
도 8 내지 도 10을 참조하면, 피검사체(5)의 이동속도와 피검사체(5)로부터 검출한 주파수 파형과의 관계를 나타낸다. 도 8은 피검사체(5)를 3[m/s]로 이동한 경우의 검출한 주파수 파형이며, 도 9는 피검사체(5)를 5[m/s]로 이동한 경우의 검출한 주파수 파형이고, 도 10은 피검사체(5)를 10[m/s]로 이동한 경우의 검출한 주파수 파형이다.8-10, the relationship between the moving speed of the
측정조건은 상기 반송속도가 상이한 것 외에는 도 4에서 설명한 측정조건이다. 군속도 제로의 S1 모드 판파 주파수인 약 1.3MHz의 주파수 스펙트럼 F5, F6, F7은 각각 피검사체의 이동속도가 3[m/s](도 8), 5[m/s](도 9), 10[m/s](도 10)이어도 검출된 것을 알 수 있다.Measurement conditions are the measurement conditions described with reference to FIG. 4 except that the conveyance speeds are different. In the frequency spectrums F5, F6, and F7 of about 1.3 MHz, which are S1 mode wave frequencies of group speed zero, the moving speeds of the subject were 3 [m / s] (Fig. 8), 5 [m / s] (Fig. 9), and 10, respectively. It can be seen that even m / s (FIG. 10) was detected.
이와 같이 도 8 내지 도 10에 나타낸 바와 같이, 피검사체(5)를 이동시킨 경우라도, 군속도 제로의 S1 모드 판파 주파수가 측정 가능한 것을 알 수 있다. 이 결과로, 실제로 강재가 콘베이어 등에 의해 이동한 경우라도 군속도 제로의 S1 모드 판파 주파수를 측정하는 것이 가능하며, 이와 같이 하여 이동한 강재로부터 검출한 군속도 제로의 S1 모드 판파 주파수를 이용하여 도 5 및 도 6에 나타낸 바와 같이 상변태율의 측정을 행할 수 있다.As shown in Fig. 8 to Fig. 10, even when the
상기와 같이 본 실시형태에 관한 레이저 초음파 측정장치 및 그 방법은 피검사체의 표면으로의 손상이 없는 펄스에너지 밀도가 낮은 레이저광에 의해 상변태율의 측정을 높은 정밀도로 행할 수 있다.As described above, the laser ultrasonic measuring apparatus and the method according to the present embodiment can measure the phase transformation rate with high precision by laser light having a low pulse energy density without damage to the surface of the inspected object.
또한, 본 발명은 이와 같이 약한 펄스레이저의 사용에만 한정되는 것이 아니며 어블레이션과 열탄성변화영역의 공존영역, 어블레이션 영역에서도 이용할 수 있다.In addition, the present invention is not limited to the use of such a weak pulse laser, but can also be used in the coexistence region and ablation region of the ablation and thermoelastic change regions.
이상 설명한 바와 같이, 본 실시형태에 관한 레이저 초음파 측정장치 및 그 방법은 대상재의 표면에 손상을 주지 않고 초음파를 발생시키는 열탄성영역에서, 검출한 초음파파형에 대해 FFT를 실행하여 주파수 피크를 산출하고 그 주파수값으로부터 변태율을 계측할 수 있다.As described above, the laser ultrasonic measuring apparatus and the method according to the present embodiment calculate the frequency peak by performing FFT on the detected ultrasonic waveform in the thermoelastic region that generates ultrasonic waves without damaging the surface of the target material. The transformation rate can be measured from the frequency value.
상변태율을 계측하는 것은 금속의 제조프로세스에서 재질을 만드는데 중요한 정보가 된다. 예를 들면, 고강도 강판의 제조에서는 γ철에서 α철로 변태시키는 조직제어를 행하여 고강도 만드는 제조방법이 있으나, 이 제조에서는 γ철에서 α철로 변태시킬 때의 비율, 즉, 변태율을 잘 조절함으로써 강도를 동일한 정도로 강도의 편차가 적은 강판을 제조할 수 있게 된다. 이와 같은 관점에서 상변태율을 용이하게 계측할 수 있는 것, 또한 실제의 금속 제조 프로세스와 같이 반송 중의 강재에 대해 상변태율을 구할 수 있는 것은 매우 유용하다.Measuring phase transformation is important information to make materials in the metal manufacturing process. For example, in the manufacture of high strength steel sheet, there is a manufacturing method for making high strength by performing the structure control of transforming from iron to α iron, but in this manufacturing, the strength by controlling the ratio, i.e. It is possible to produce a steel sheet with less variation in strength to the same extent. From this point of view, it is very useful to be able to easily measure the phase transformation rate, and to be able to obtain the phase transformation rate with respect to the steel material being conveyed like the actual metal manufacturing process.
도 11은 상변태율의 산출처리를 나타내는 흐름도이다. 도 11(a)는 주파수 스펙트럼 상의 주파수 피크와 피검사체 두께와의 곱과 상변태율의 관계를 나타내는 교정선을 작성하는 흐름이다. 도11(b)는 검출한 주파수 스펙트럼 상의 주파수 피크와 피검사체 두께와의 곱과 교정선으로부터 상변태율을 산출하는 흐름이다.11 is a flowchart showing the calculation process of the phase transformation rate. FIG. 11 (a) is a flow which produces the correction line which shows the relationship between the product of the frequency peak in a frequency spectrum, the thickness of a test subject, and a phase transformation rate. Fig. 11B is a flow for calculating the phase transformation rate from the product of the frequency peak on the detected frequency spectrum and the thickness of the subject and the calibration line.
도 11(a)에 나타낸 흐름에서, 레이저 초음파 측정장치(10)는 도2 및 도4에서 설명한 바와 같이, 피검사체로부터 발생한 군속도 제로의 S1 모드 판파 주파수를 산출한다(S101). 다음으로, 레이저 초음파 측정장치(10)는 피검사체의 상변태율을 얻는다(S102). 여기서는 상술한 바와 같이 미리 전자현미경으로 측정 또는 열팽창률로 산출한 상변태율을 얻는다.In the flow shown in FIG. 11A, the laser
스텝 S102에서 소정의 상변태율이 산출되면, 도 6에 나타낸 바와 같이 레이저 초음파 측정장치(10)는 피검사체의 상변태율과 주파수 스펙트럼 상의 주파수 피크 및 피검사체 두께로 교정선을 작성한다(S103). 이와 같이 하여 교정선이 작성된다.When the predetermined phase transformation rate is calculated in step S102, as shown in FIG. 6, the laser
도 11(b)에 나타낸 흐름에서, 레이저 초음파 측정장치(10)는 도 2 및 도 4에서 설명한 바와 같이, 피검사체에서 발생한 군속도 제로의 S1 모드 판파 주파수를 산출한다(S151). 다음으로, 레이저 초음파 측정장치(10)는 스텝 S103에서 작성된 교정선 및 피검사체 두께를 이용하여 검출된 주파수 스펙트럼 상의 주파수 피크에 상당하는 상변태율을 산출한다(S152).In the flow shown in FIG. 11B, the laser
도 12는 일 실시형태에 관한 레이저 초음파 측정장치(10a)의 전체도이다. 도시한 바와 같이 레이저 초음파 측정장치(10a)는 초음파 발생용 레이저광원(11), 초음파 검출용 레이저광원(12), 레이저 간섭계(13) 및 신호처리부(60)를 갖는다.12 is an overall view of a laser
초음파 발생용 레이저광원(11)은 피검사체(5)에 초음파를 발생시키기 위해 고출력 펄스 레이저광(EL)을 밀러(31a, 31b)를 통해 피검사체(5)에 조사한다. 펄스 레이저광(EL)은 렌즈(33)에 의해 집속되어 점상 스폿으로 초음파 검출용 레이저 스폿영역(7)에 조사된다. 피검사체(5)의 펄스 레이저광(EL)이 조사되는 부분은 열팽창하고, 그 후 수축함으로써 왜곡이 발생하여 피검사체(5)에 초음파가 전달된다.The ultrasonic
초음파 검출용 레이저광원(12)으로부터는 연속파 레이저광(DL)이 밀러(41a, 41b)를 통해 점상 스폿으로 초음파 검출용 레이저 스폿영역(7)에 조사된다. 연속파 레이저광(DL)은 피검사체(5) 상의 계측점에 조사된다. 피검사체(5)로부터의 연속파 레이저광(DL)의 레이저 반사광(RL)의 주파수는 표면진동에 의해 도플러 시프트된다. 레이저 반사광(RL)은 밀러(41c, 41d)를 통해 렌즈(43)에 의해 집속되어 레이저 간섭계(13)에 입사한다.The continuous wave laser beam DL is irradiated from the ultrasonic wave detection
레이저 간섭계(13)는 입사한 레이저광을 밀러 13a와 13b 사이에서 왕복시켜 간섭(공진)시키고, 밀러(13b)로 일부 광량을 투과시킨다. 레이저 간섭계(13)를 투과하는 레이저광의 스펙트럼은 매우 좁은 파장영역의 급격한 피크와 양측의 슬로프로 구성된다. 밀러 13a와 13b의 간격을 레이저광의 파장의 공진조건에서 약간 벗어나게 하고 레이저광의 파장이 해당 슬로프의 파장영역에 위치하도록 설정하여 둔다. 이때, 레이저 반사광의 도플러 시프트는 레이저 간섭계(13)를 투과할 때에 광강도 변화로 변환된다. 그리고, 레이저 간섭계(13)를 투과한 광은 애벌란시 포토다이오드 등으로 이루어진 고속 응답이 가능한 광검출기(20)에 입사한다. 광검출기(20)에서는 투과광의 광강도는 전기신호(E1)로 변환되어 신호처리부(60)에 입력 한다.The
신호처리부(60)에서 주파수 산출부(61)는 도 11의 스텝 S101 및 S151에서 설명한 처리를 실행할 수 있다. 상변태율 격납부(62)는 미리 측정 또는 산출한 상변태율을 격납한다. 두께정보 격납부(63)는 미리 측정한 피검사체 두께의 값을 격납한다. 또한, 이값은 방사선계측기에 의한 계측 및 메저(measure)류를 사용한 계측에 의해 가능하다. 교정선 작성부(64)는, 도 11의 스텝 S102에서 설명한 바와 같이, 상변태율 격납부(62)에 격납된 상변태율을 얻어 도 11의 스텝 S103에서 설명한 바와 같이 교정선을 작성할 수 있다. 상변태율 산출부(65)는 도 11의 스텝 S152에서 설명한 처리를 실행할 수 있다.In the
신호처리부(60)는 프로세서로 구성되며 상기 주파수 산출부(61), 교정선 작성부(64), 상변태율 산출부(65)는 프로세서가 도시되어 있지 않는 메모리나 외부기억장치에 격납된 프로그램을 실행함으로써 실장해도 된다. 상변태율 격납부(62)는 메모리나 외부기억장치로 실장할 수 있다. 상변태율 산출부(65)의 산출결과는 표시장치(19)에 출력할 수 있다. 또한, 신호처리부(60)에는 도시하지 않은 입출력장치, 통신장치가 접속될 수 있다. 측정한 상변태율은 입력장치 또는 통신장치를 통해 상변태율 격납부(62)에 격납할 수 있다.The
또한, 레이저 초음파 측정장치(10a)는 금속 제조 프로세스에 배치되는 경우, 교정선의 작성기능을 가지지 않고, 미리 작성된 교정선을 이용하여 상변태율을 구해도 된다. 이 경우, 교정선 작성부(64)의 기능을 레이저 초음파 측정장치(10a)와 데이터 접속된 별도의 장치(프로세서, 메모리, 통신수단, 외부기억장치를 갖는다) 에 실행해도 된다. 그리고, 신호처리부(60)는 상변태율 격납부(62) 및 교정선 작성부(64)를 가지지 않고, 상기 별도의 장치에서 네트워크 또는 기록매체를 통해 수신한 교정선의 데이터를 이용하고 상변태율 산출부(65)는 도 11의 스텝 S152에서 설명한 처리를 실행할 수 있다.When the laser
본 발명은 상술한 실시형태 및 첨부된 도면에 의해 한정되지 아니한다. 첨부된 청구범위에 의해 권리범위를 한정하고자 하며, 청구범위에 기재된 본 발명의 기술적 사상을 벗어나지 않는 범위 내에서 다양한 형태의 치환, 변형 및 변경이 가능하다는 것은 당 기술분야의 통상의 지식을 가진 자에게 자명할 것이다.The present invention is not limited by the above-described embodiment and the accompanying drawings. It will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims, .
도 1(a)는 피검사체(5)의 측면에서 본 레이저 초음파 측정장치(10)의 계략 구성도이며, 도 1(b)는 피검사체(5)를 레이저 초음파 측정장치(10)에서 본 도면이다.FIG. 1 (a) is a schematic configuration diagram of the laser
도 2(a)는 펄스레이저광을 조사했을 때의 검출파형을 나타내고, 도2(b)는 도 2(a)에 나타낸 시간파형을 푸리에 변환한 주파수특성을 나타낸 도면이다. Fig. 2 (a) shows the detection waveform when the pulsed laser light is irradiated, and Fig. 2 (b) shows the frequency characteristic obtained by Fourier transforming the time waveform shown in Fig. 2 (a).
도 3은 논리계산에 의한 S1 모드 판파 주파수와 위상속도 및 군속도의 관계를 나타낸 도면이다. 3 is a diagram showing the relationship between the S1 mode wave frequency, phase velocity, and group velocity by logic calculation.
도 4는 레이저 초음파 측정장치(10)에 의해 측정온도조건을 바꿔서검출한 파형의 FFT결과를 나타낸 도면이다. 4 is a view showing the FFT result of the waveform detected by changing the measurement temperature conditions by the laser
도 5는 피검사체의 온도와 계측한 군속도 제로의 S1 모드 판파 주파수와의 관계를 나타낸 도면이다. It is a figure which shows the relationship between the temperature of a test subject and the S1 mode plate | board frequency of the measured group speed zero.
도 6(a)는 상변태율과 군속도 제로의 S1 모드 판파 주파수와 피검사체 두께와의 곱의 관계를 선형의 교정선으로 나타낸 도면이며, 도 6(b)는 상변태율과 S군속도 제로의 S1 모드 판파 주파수와 피검사체 두께와의 곱의 관계를 비선형 교정선으로 나타낸 도면이다. Fig. 6 (a) is a diagram showing a linear calibration line showing the relationship between the phase transformation rate and the S1 mode plate frequency of group velocity zero and the thickness of the subject, and Fig. 6 (b) shows the S1 mode of phase transformation rate and S group velocity zero. It is a figure which shows the relationship of the product of a plate wave frequency and the thickness of a subject with a nonlinear correction line.
도 7은 피검사체를 이동하면서 레이저 초음파측정법을 적용하는 경우의 초음파 발생용 레이저 스폿과 초음파 검출용 레이저 스폿의 위치관계를 나타낸 도면이다. FIG. 7 is a view showing the positional relationship between the laser spot for ultrasonic generation and the laser spot for ultrasonic detection when the laser ultrasonic measuring method is applied while moving the inspected object. FIG.
도 8은 피검사체의 이동속도가 3m/s인 경우의 검출파형의 주파수 특성을 나 타낸 도면이다. 8 is a diagram showing the frequency characteristics of the detected waveform when the moving speed of the subject is 3 m / s.
도 9는 피검사체의 이동속도가 5m/s인 경우의 검출파형의 주파수 특성을 나타낸 도면이다. 9 is a diagram showing the frequency characteristics of the detection waveform when the moving speed of the subject is 5 m / s.
도 10은 피검사체의 이동속도가 10m/s인 경우의 검출파형의 주파수 특성을 나타낸 도면이다. 10 is a diagram showing the frequency characteristics of the detected waveform when the moving speed of the object under test is 10 m / s.
도 11(a)는 주파수 스펙트럼 상의 주파수 피크와 피검사체 두께와의 곱과 상변태율의 관계를 나타내는 교정선을 작성하는 흐름도이며, 도 11(b)는 검출한 주파수 스펙트럼 상의 주파수 피크, 피검사체 두께와 교정선으로부터 상변태율을 산출하는 흐름도이다.Fig. 11A is a flowchart for creating a calibration line showing the relationship between the product of the frequency peak on the frequency spectrum and the thickness of the test object and the phase transformation rate, and Fig. 11B is the frequency peak and test thickness on the detected frequency spectrum. Is a flowchart for calculating the phase transformation rate from and the calibration line.
도 12는 레이저 초음파 측정장치의 전체 구성을 설명하는 도면이다.It is a figure explaining the whole structure of a laser ultrasonic measuring apparatus.
도 13는 레이저 초음파법의 측정원리에 대해 설명하는 도면이다.It is a figure explaining the measuring principle of a laser ultrasonic method.
도 14는 레이저광조사에 의한 초음파의 2종류의 발생 원리를 설명하기 위한 도면이다.It is a figure for demonstrating the two types of generation | occurrence principle of the ultrasonic wave by laser beam irradiation.
<도면의 주요부분에 대한 부호의 설명> <Description of the symbols for the main parts of the drawings>
5 : 피검사체 5:: test subject
6, 6a : 초음파 발생용 레이저 스폿영역6, 6a: Laser spot area for ultrasonic generation
7, 7a : 초음파 검출용 레이저 스폿영역 7, 7a: Laser spot area for ultrasonic detection
10, 10a : 레이저 초음파 측정장치10, 10a: laser ultrasonic measuring device
11 : 초음파 발생용 레이저광원 11:: laser light source for ultrasonic generation
12 : 초음파 검출용 레이저광원12 ms: laser light source for ultrasonic detection
13 : 레이저 간섭계 13: Laser interferometer
13a, 13b, 31a, 41a, 41b, 41c : 밀러13a, 13b, 31a, 41a, 41b, 41c
19 : 표시장치 19: Display device
20 : 광검출기 20: Photodetector
33, 43 : 렌즈33, 43 ': Lens
60 : 신호처리부 60Hz: Signal processing part
61 : 주파수 산출부 61 Hz: Frequency output section
62 : 상변태율 격납부62: phase transformation rate storing part
63 : 두께정보 격납부63: thickness information storage
64 : 교정선 작성부 64: Calibration line preparation department
65 : 상변태율 산출부65 °: Phase transformation rate calculation unit
Claims (14)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090024585A KR101053415B1 (en) | 2009-03-23 | 2009-03-23 | Laser Ultrasonic Measuring Device and Measuring Method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090024585A KR101053415B1 (en) | 2009-03-23 | 2009-03-23 | Laser Ultrasonic Measuring Device and Measuring Method |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20100106120A true KR20100106120A (en) | 2010-10-01 |
KR101053415B1 KR101053415B1 (en) | 2011-08-01 |
Family
ID=43128542
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020090024585A KR101053415B1 (en) | 2009-03-23 | 2009-03-23 | Laser Ultrasonic Measuring Device and Measuring Method |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101053415B1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101242873B1 (en) * | 2011-06-02 | 2013-03-12 | 니뽄스틸코포레이션 | Method and apparatus for measuring stress of strip sheet |
KR101253913B1 (en) * | 2011-08-01 | 2013-04-16 | 주식회사 포스코 | Laser ultrasonic apparatus having multiple beam head and method of processing signal using thereof |
KR101253927B1 (en) * | 2011-08-01 | 2013-04-16 | 주식회사 포스코 | Method and apparatus for measuring quality of steel plate using laser ultrasonic |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102267862B1 (en) * | 2015-03-31 | 2021-06-23 | 삼성전자주식회사 | Examining apparatus and examinig method thereof |
JP6776181B2 (en) * | 2017-05-31 | 2020-10-28 | 株式会社神戸製鋼所 | Stress measurement method |
JP6815933B2 (en) * | 2017-05-31 | 2021-01-20 | 株式会社神戸製鋼所 | Stress measurement method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4386709B2 (en) | 2003-11-14 | 2009-12-16 | 関西電力株式会社 | Material nondestructive inspection method and apparatus by laser ultrasonic wave |
KR100733539B1 (en) * | 2004-12-21 | 2007-06-28 | 재단법인 포항산업과학연구원 | Apparatus and method of laser-ultrasonic measurement for hot object |
JP4649280B2 (en) | 2005-07-07 | 2011-03-09 | 株式会社東芝 | Laser ultrasonic receiving apparatus and laser ultrasonic receiving method |
-
2009
- 2009-03-23 KR KR1020090024585A patent/KR101053415B1/en active IP Right Grant
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101242873B1 (en) * | 2011-06-02 | 2013-03-12 | 니뽄스틸코포레이션 | Method and apparatus for measuring stress of strip sheet |
KR101253913B1 (en) * | 2011-08-01 | 2013-04-16 | 주식회사 포스코 | Laser ultrasonic apparatus having multiple beam head and method of processing signal using thereof |
KR101253927B1 (en) * | 2011-08-01 | 2013-04-16 | 주식회사 포스코 | Method and apparatus for measuring quality of steel plate using laser ultrasonic |
Also Published As
Publication number | Publication date |
---|---|
KR101053415B1 (en) | 2011-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6057927A (en) | Laser-ultrasound spectroscopy apparatus and method with detection of shear resonances for measuring anisotropy, thickness, and other properties | |
JP6034468B1 (en) | Real-time non-contact non-destructive thickness measurement system using terahertz waves | |
KR101053415B1 (en) | Laser Ultrasonic Measuring Device and Measuring Method | |
JP4783263B2 (en) | Ultrasonic multi-echo measurement device | |
JP5058109B2 (en) | Method and apparatus for measuring longitudinal and transverse sound velocities in materials by laser ultrasonic method | |
JP5528385B2 (en) | Poisson's ratio measuring method and measuring device | |
JP5580426B2 (en) | Metal structure and material measuring device and measuring method | |
JP4106400B2 (en) | Thickness measuring device and thickness measuring method | |
Li et al. | Non destructive testing thickness measurement by laser ultrasound under high temperature | |
KR101180151B1 (en) | Measuring method of poisson's ratio and measuring apparatus thereof | |
KR100496826B1 (en) | Apparatus and method of noncontact measurement of crystal grain size | |
KR100993989B1 (en) | Laser ultrasonic measuring device and laser ultrasonic measuring method | |
JP5072789B2 (en) | Method and apparatus for measuring longitudinal and transverse sound velocities in materials by laser ultrasonic method | |
JP5058196B2 (en) | Apparatus and method for measuring phase transformation rate of material | |
KR101257203B1 (en) | Apparatus and method for evaluating micro damage of materials using nonlinear laser-generated surface wave | |
JP5528395B2 (en) | Thin plate stress measuring method and measuring apparatus | |
JP3704843B2 (en) | Non-contact non-destructive material evaluation method and apparatus, elastic wave excitation method and elastic wave excitation apparatus | |
KR100994037B1 (en) | Laser ultrasonic measuring device and laser ultrasonic measuring method | |
González et al. | Generation of sub-microsecond quasi-unipolar pressure pulses | |
JP4471714B2 (en) | Crystal grain size distribution measuring method and apparatus | |
JP5419677B2 (en) | Poisson's ratio measuring method and measuring device | |
KR100733539B1 (en) | Apparatus and method of laser-ultrasonic measurement for hot object | |
Oe et al. | Nondestructive internal defect detection using photoacoustic and self‐coupling effect | |
KR101538908B1 (en) | Safety monitoring system of nuclear equipment with intergrating system of laser ultrasonic and shearography | |
Imano | Detection of Unbonded Defect under Surface of Material Using Phase Information of Rayleigh and A0 Mode Lamb Waves. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20140728 Year of fee payment: 4 |
|
FPAY | Annual fee payment |
Payment date: 20150727 Year of fee payment: 5 |
|
FPAY | Annual fee payment |
Payment date: 20160726 Year of fee payment: 6 |
|
FPAY | Annual fee payment |
Payment date: 20170726 Year of fee payment: 7 |
|
FPAY | Annual fee payment |
Payment date: 20180726 Year of fee payment: 8 |
|
FPAY | Annual fee payment |
Payment date: 20190724 Year of fee payment: 9 |