RU2540942C1 - Method to monitor dynamics of variation of controlled object wall thickness - Google Patents
Method to monitor dynamics of variation of controlled object wall thickness Download PDFInfo
- Publication number
- RU2540942C1 RU2540942C1 RU2013149631/28A RU2013149631A RU2540942C1 RU 2540942 C1 RU2540942 C1 RU 2540942C1 RU 2013149631/28 A RU2013149631/28 A RU 2013149631/28A RU 2013149631 A RU2013149631 A RU 2013149631A RU 2540942 C1 RU2540942 C1 RU 2540942C1
- Authority
- RU
- Russia
- Prior art keywords
- wall thickness
- frequency
- value
- lamb wave
- thickness
- Prior art date
Links
Images
Landscapes
- Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
Заявленное изобретение относится к области технической диагностики и неразрушающего контроля промышленных объектов и используется для контроля за динамикой изменения минимального значения толщины стенки тонкостенных и листовых изделий, а также других изделий, в которых могут распространяться волны Лэмба, например трубопроводов, резервуаров, сосудов, цистерн.The claimed invention relates to the field of technical diagnostics and non-destructive testing of industrial facilities and is used to control the dynamics of changes in the minimum wall thickness of thin-walled and sheet products, as well as other products in which Lamb waves can propagate, for example pipelines, tanks, vessels, tanks.
Уровень техники.The level of technology.
Известен способ измерения толщины изделия (патент РФ №2231753, опубл. 27.06.2004) с помощью ультразвуковых импульсов, заключающийся в том, что с помощью ультразвукового преобразователя излучают в материал изделия ультразвуковой импульс, принимают из него эхо-импульсы, отраженные от противоположных поверхностей материала, измеряют время распространения эхо-импульсов от одной поверхности до другой и обратно и вычисляют толщину как произведение половины этого времени на известную скорость распространения ультразвуковых импульсов в материале.A known method of measuring the thickness of the product (RF patent No. 2231753, publ. June 27, 2004) using ultrasonic pulses, which consists in the fact that using an ultrasonic transducer emit an ultrasonic pulse into the product material, receive echo pulses from it reflected from opposite surfaces of the material measure the propagation time of the echo pulses from one surface to another and vice versa and calculate the thickness as the product of half this time by the known velocity of propagation of ultrasonic pulses in the material.
Необходимость обеспечения физического доступа к поверхности всех контролируемых участков изделия и высокая трудоемкость контроля больших по площади изделий снижают эффективность известного способа.The need to provide physical access to the surface of all controlled areas of the product and the high complexity of control of large area products reduce the effectiveness of the known method.
Известен способ (патент WO 2007139389, опубл. 06.12.2007) измерения величины потери средней толщины материала твердой структуры, способной поддерживать распространение акустических волн и во время работы контактирующей с одним или несколькими веществами, приводящими к изменению толщины структуры, заключающийся в использовании множества датчиков любого типа, установленных в контакте с поверхностью указанной твердой структуры, и отличающийся тем, что включен один или несколько из следующих шагов:A known method (patent WO 2007139389, published 06.12.2007) for measuring the loss of the average thickness of a material of a solid structure capable of supporting the propagation of acoustic waves and in contact with one or more substances leading to a change in the thickness of the structure, which consists in using a variety of sensors of any type installed in contact with the surface of the specified solid structure, and characterized in that it includes one or more of the following steps:
а) выполнение один или большее число раз этапа передачи акустического сигнала по крайней мере от одного из множества акустических датчиков в стенку указанной структуры и прием переданного акустического сигнала множеством указанных датчиков после прохождения дистанции в твердой структуре;a) performing one or more times the step of transmitting an acoustic signal from at least one of the plurality of acoustic sensors to a wall of said structure and receiving a transmitted acoustic signal by a plurality of said sensors after covering a distance in a solid structure;
б) сбор сигналов с указанного множества датчиков в блок обработки;b) collecting signals from the specified set of sensors in the processing unit;
в) анализ и сравнение принятых сигналов, сгенерированных на шаге (а), и/или результатов, соответствующих более ранним измерениям;c) analysis and comparison of the received signals generated in step (a) and / or the results corresponding to earlier measurements;
г) вычисление степени потери толщины материала указанной твердой структуры по результатам анализа, выполненного на шаге (в),g) calculating the degree of loss of thickness of the material of the specified solid structure according to the results of the analysis performed in step (c),
причем блок обработки данных способен вычислять среднее снижение истинной толщины стенки в пределах области твердой структуры путем анализа сигналов, принятых преобразователями.moreover, the data processing unit is able to calculate the average decrease in the true wall thickness within the region of the solid structure by analyzing the signals received by the converters.
Недостатком известного способа является необходимость измерения средней толщины при установке датчиков на объект для последующего получения текущих значений средней толщины, что значительно снижает эффективность известного способа.The disadvantage of this method is the need to measure the average thickness when installing sensors on the object for the subsequent receipt of the current values of the average thickness, which significantly reduces the effectiveness of the known method.
Известен способ (патент US 5965818, опубл. 12.10.1999) измерения минимального значения толщины стенки трубы, лежащей на опоре, при помощи ультразвуковых волн Лэмба, включающийA known method (patent US 5965818, publ. 12.10.1999) measuring the minimum thickness of the pipe wall lying on a support using ultrasonic Lamb waves, including
а) установку передающего преобразователя и приемного преобразователя на указанной трубе на противоположных сторонах от опоры указанной трубы;a) installing a transmitting converter and a receiving converter on the specified pipe on opposite sides of the support of the specified pipe;
б) излучение в указанную трубу указанным передающим преобразователем высокодисперсной и монотонной части волны Лэмба S0, у которой значения произведения частоты на толщину равны или меньше значения, соответствующего минимуму групповой скорости;b) radiation into the specified tube by the specified transmitting transducer of the finely dispersed and monotonous part of the Lamb wave S 0 , in which the values of the product of the frequency by the thickness are equal to or less than the value corresponding to the minimum group velocity;
в) прием указанных волн Лэмба указанным приемным преобразователем;c) receiving said Lamb waves by said receiving transducer;
г) измерение времени, требуемого для того, чтобы указанная излученная волна Лэмба была принята указанным приемным преобразователем, указанное время далее называется время пролета;d) measuring the time required for the specified radiated Lamb wave to be received by the specified receiving transducer, the specified time is hereinafter referred to as the time of flight;
д) повтор шагов от а) до г) для получения реперного значения времени пролета без изменения настроек прибора, но с перемещением указанных преобразователей в точки, расположенные рядом с указанной опорой трубы;e) repeating steps a) to d) to obtain a reference value of the time of flight without changing the settings of the device, but with the movement of these transducers to points located next to the specified pipe support;
е) определение изменения величины времени пролета, вызванного локальными коррозионными потерями в стенке, путем вычитания времени пролета, измеренного на опоре трубы, из реперного значения времени пролета иf) determining the change in the value of the transit time caused by local corrosion losses in the wall by subtracting the transit time measured on the pipe support from the reference value of the transit time and
ж) определение толщины трубы на указанной опоре трубы с использованием указанного изменения величины времени пролета.g) determining the thickness of the pipe on the specified pipe support using the specified change in the value of the time of flight.
Недостатком известного способа является невозможность определения минимального значения толщины стенки за пределами области соприкосновения трубы и опоры, что значительно снижает эффективность способа.The disadvantage of this method is the impossibility of determining the minimum value of the wall thickness outside the contact area of the pipe and support, which significantly reduces the effectiveness of the method.
Наиболее близким по технической сущности является способ (патент US 6360609, опубл. 26.03.2002) измерения толщины в многомодовых дисперсионных условиях, заключающийся в том, что излучают в первой точке в среду многомодовую обладающую дисперсией акустическую нормальную волну, принимают во второй точке многомодовую обладающую дисперсией акустическую нормальную волну, с использованием частотно-временного анализа проводят частотное разложение многомодовой обладающей дисперсией акустической нормальной волны, выбирают частоту, которая определяет групповую задержку tg нормальной волны на каждой имеющейся на этой частоте моде, выбирают моду и определяют групповые задержки по сравнению с первой точкой, определяют расстояние L между первой и второй точкой, выбирают моду и определяют групповую задержку tg по сравнению с первой точкой, вычисляют групповую скорость cg=L/tg, определяют толщину по дисперсионному соотношению.The closest in technical essence is the method (patent US 6360609, publ. March 26, 2002) for measuring the thickness under multimode dispersion conditions, namely, that a multimode dispersive acoustic normal wave is emitted at a first point in a medium, and a multimode dispersion is received at a second point acoustic normal wave, using time-frequency analysis, carry out the frequency decomposition of a multimode dispersive acoustic normal wave, select the frequency that determines the group delay t g of the normal wave at each mode available at this frequency, select the mode and determine group delays compared to the first point, determine the distance L between the first and second points, select the mode and determine the group delay t g compared to the first point, calculate group speed c g = L / t g , determine the thickness by the dispersion relation.
Недостатком известного способа является условие постоянства на всем пути между первой и второй точкой величины групповой скорости выбранной моды на выбранной частоте, которое не выполняется при наличии вариаций значений толщины между первой и второй точкой, к примеру при наличии на объекте очаговой коррозии.The disadvantage of this method is the condition of constancy on the whole path between the first and second point of the group velocity of the selected mode at the selected frequency, which is not satisfied if there are variations in the thickness values between the first and second point, for example, if there is focal corrosion on the object.
Раскрытие изобретения.Disclosure of the invention.
Технический результат, достигаемый от реализации заявленного решения, заключается в снижении трудоемкости, упрощении и ускорении контроля толщины стенки больших по площади объектов, устранении необходимости получения физического доступа ко всей поверхности объекта, подлежащей контролю толщины, обеспечении возможности контроля толщины при наличии вариаций значений толщины, к примеру при наличии на объекте очаговой коррозии, обеспечении возможности определения среднеарифметического и минимального значений толщины на участке между двумя акустическими преобразователями, установленными на объекте в произвольных точках, без проведения предварительного измерения толщины в тех же точках.The technical result achieved from the implementation of the claimed solution is to reduce the complexity, simplify and accelerate the control of the wall thickness of large objects, eliminating the need for physical access to the entire surface of the object to be controlled by thickness, providing the possibility of controlling the thickness in the presence of variations in the thickness values, to example, if there is focal corrosion at the facility, providing the ability to determine the arithmetic mean and minimum thickness values in the area between vumya acoustic transducers mounted at the object at arbitrary points, without a preliminary thickness measurements at the same points.
Сущность способа заключается в том, что размещают на его поверхности на известном расстоянии друг от друга, по крайней мере, один акустический преобразователь для излучения волн Лэмба и, по крайней мере, один преобразователь для их приема, излучают в заданный момент времени импульсный сигнал, производят расчет зависимости спектральной плотности мощности принятого сигнала от времени, выбирают волны Лэмба и частоты, определяют разность между временем приема выбранной частотной составляющей выбранной волны Лэмба и временем излучения сигнала, определяют значения групповой скорости выбранной частотной составляющей выбранной волны Лэмба по известному значению расстояния между преобразователями и значению разности между временем приема выбранной частотной составляющей выбранной волны Лэмба и временем излучения сигнала, определяют значения толщины стенки по полученному значению групповой скорости, выбранному значению частоты и эталонной зависимости групповой скорости выбранной волны Лэмба от произведения толщины стенки и частоты, при этом устанавливают минимальную величину толщины стенки по полученным среднеарифметической величине толщины стенки и дисперсии значений толщины стенки объекта, причем среднеарифметическое значение толщины стенки определяют с выбором симметричной волны Лэмба нулевого порядка и, по крайней мере, одной частоты, на которой эталонная зависимость обратной величины групповой скорости симметричной волны Лэмба нулевого порядка от произведения толщины стенки и частоты близка к линейной, а дисперсию значений толщины стенки определяют с выбором антисимметричной волны Лэмба нулевого порядка и частоты, на которой эталонная зависимость обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка от произведения толщины стенки и частоты существенно нелинейна.The essence of the method lies in the fact that at least one acoustic transducer for emitting Lamb waves and at least one transducer for receiving them is placed on its surface at a known distance from each other, they emit a pulse signal at a given time, produce Calculation of the dependence of the spectral power density of the received signal on time, Lamb waves and frequencies are selected, the difference between the reception time of the selected frequency component of the selected Lamb wave and the signal emission time is determined ala, determine the group velocity values of the selected frequency component of the selected Lamb wave from the known distance between the transducers and the difference between the reception time of the selected frequency component of the selected Lamb wave and the signal emission time, determine the wall thickness values from the obtained group velocity value, the selected frequency value and the reference the dependence of the group velocity of the selected Lamb wave on the product of the wall thickness and frequency, while setting the minimum the wall thickness value from the obtained arithmetic mean wall thickness and variance of the object wall thickness values, the arithmetic mean wall thickness being determined with the choice of a symmetric Lamb wave of zero order and at least one frequency at which the reference dependence of the reciprocal of the group velocity of a symmetric Lamb wave of zero of the order of the product of the wall thickness and frequency is close to linear, and the variance of the wall thickness values is determined with the choice of the antisymmetric Lamb wave Nya Ullevi order and the frequency at which the reference dependence of the reciprocal of the group velocity of the antisymmetric Lamb wave of zero order of the product of wall thickness and frequency is strongly nonlinear.
Краткое описание чертежей.A brief description of the drawings.
На фиг.1 приведены графики зависимостей групповых скоростей симметричной волны Лэмба нулевого порядка S0 и антисимметричной волны Лэмба нулевого порядка A0 от произведения частоты f и толщины h стенки объекта. На фиг.2 приведен пример принятого акустического сигнала. На фиг.3 приведен пример зависимости спектральной плотности мощности принятого сигнала от времени. На фиг.4 приведены пример выявления прихода симметричной волны Лэмба нулевого порядка S0 и антисимметричной волны Лэмба нулевого порядка A0 на графике зависимости спектральной плотности мощности сигнала от времени, пример двух частотных диапазонов, в пределах которых зависимость обратной величины групповой скорости симметричной волны Лэмба нулевого порядка S0 от частоты близка по форме к линейной зависимости, и пример частотного диапазона, в пределах которого зависимость обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка A0 от частоты существенно нелинейна. Также на фиг.4 приведены примеры определения времен приема выбранных волн Лэмба на частотах, принадлежащих указанным участкам. На фиг.5 приведен пример графика значений производной второго порядка обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка по произведению частоты f и толщины h стенки. На фиг.6 приведен пример нахождения значения произведения толщины стенки и частоты, при котором эталонное значение групповой скорости симметричной волны Лэмба нулевого порядка S0 равняется полученному путем анализа принятого сигнала значению групповой скорости.Figure 1 shows graphs of the dependences of the group velocities of a symmetric zero-order Lamb wave S 0 and an antisymmetric zero-order Lamb wave A 0 from the product of the frequency f and the thickness h of the object wall. Figure 2 shows an example of a received acoustic signal. Figure 3 shows an example of the dependence of the spectral power density of the received signal on time. Figure 4 shows an example of detecting the arrival of a zero-order symmetric Lamb wave S 0 and a zero-order antisymmetric Lamb wave A 0 on a graph of the signal power spectral density versus time, an example of two frequency ranges within which the reciprocal of the group velocity of a symmetric Lamb wave of zero S 0 order of the frequency is close in form to the linear relationship, and an example of the frequency range within which the dependence of the reciprocal of the group velocity of the antisymmetric wave Le ba A zero order 0 of frequency substantially linear. Figure 4 also shows examples of determining the reception times of the selected Lamb waves at frequencies belonging to these sections. Figure 5 shows an example of a graph of the values of the second-order derivative of the reciprocal of the group velocity of the antisymmetric Lamb wave of zero order with respect to the product of the frequency f and the wall thickness h. Figure 6 shows an example of finding the product of the wall thickness and the frequency at which the reference value of the group velocity of the symmetric Lamb wave of zero order S 0 is equal to the group speed value obtained by analyzing the received signal.
Осуществление изобретения.The implementation of the invention.
Способ осуществляют следующим образом.The method is as follows.
Размещают на поверхности стенки объекта, например трубопровода, на известном расстоянии L друг от друга, по крайней мере, один акустический преобразователь для излучения волн Лэмба и, по крайней мере, один преобразователь для их приема.At least one acoustic transducer for emitting Lamb waves and at least one transducer for receiving them is placed on a wall surface of an object, such as a pipeline, at a known distance L from each other.
Излучатель испускает в заданный момент времени импульсный акустический сигнал.The emitter emits a pulsed acoustic signal at a given point in time.
Принимают сигнал приемником волн Лэмба.Receive the signal by the Lamb wave receiver.
Выполняют расчет зависимости спектральной плотности мощности принятого сигнала от времени.The calculation of the dependence of the spectral power density of the received signal on time.
Определяют минимальное значение толщины стенки объекта hmin на участке между преобразователями. Для этого определяют среднеарифметическое значение M[h] толщины стенки объекта на участке между преобразователями и дисперсию значений D[h] толщины стенки объекта на участке между преобразователями. Для этого определяют среднеарифметическое значение M[h] толщины стенки объекта на участке между преобразователями. Для этого выбирают первое значение частоты f1 такое, что в его окрестности эталонная зависимость обратной величины групповой скорости симметричной волны Лэмба нулевого порядка S0 от произведения толщины стенки и частоты близка к линейной зависимости.The minimum value of the object wall thickness h min is determined in the area between the transducers. To do this, determine the arithmetic mean value M [h] of the wall thickness of the object in the area between the converters and the variance of the values D [h] of the wall thickness of the object in the area between the converters. To do this, determine the arithmetic mean value M [h] of the wall thickness of the object in the area between the transducers. To do this, choose the first value of the frequency f 1 such that in its vicinity the reference dependence of the reciprocal of the group velocity of the symmetric Lamb wave of zero order S 0 on the product of the wall thickness and frequency is close to a linear dependence.
Определяют первую разницу времен t1 между временем приема частотной составляющей симметричной волны Лэмба нулевого порядка S0 на первой выбранной частоте f1 и временем излучения по зависимости спектральной плотности мощности принятого сигнала от времени.The first time difference t 1 is determined between the time of receiving the frequency component of the symmetric Lamb wave of zero order S 0 at the first selected frequency f 1 and the radiation time according to the dependence of the spectral power density of the received signal on time.
По значению первой разницы времен t1 и по значению расстояния L между преобразователями определяют значение групповой скорости СГР выбранной частотной составляющей симметричной волны Лэмба нулевого порядка S0 по формулеThe value of the first time difference t 1 and the value of the distance L between the transducers determine the value of the group velocity С ГР of the selected frequency component of the symmetric Lamb wave of zero order S 0 according to the formula
По известным формулам рассчитывают эталонные зависимости от произведения частоты f и толщины h значений групповой скорости
Подбирают значение толщины M[h], при котором значение эталонной групповой скорости симметричной волны Лэмба нулевого порядка S0 на первой выбранной частоте f1 равняется полученному путем анализа принятого сигнала значению групповой скоростиThe thickness value M [h] is selected at which the value of the reference group velocity of the symmetric Lamb wave of zero order S 0 at the first selected frequency f 1 is equal to the group speed value obtained by analyzing the received signal
Полученное значение толщины M[h] считают среднеарифметическим значением толщины стенки на участке между преобразователями.The obtained value of the thickness M [h] is considered the arithmetic mean value of the wall thickness in the area between the transducers.
Затем определяют дисперсию значений D[h] толщины стенки объекта на участке между преобразователями. Для этого выбирают второе значение частоты f2 такое, что в его окрестности эталонная зависимость обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка A0 от произведения толщины стенки и частоты существенно нелинейна.Then determine the variance of the values D [h] of the wall thickness of the object in the area between the transducers. To do this, choose the second value of the frequency f 2 such that in its vicinity the reference dependence of the reciprocal of the group velocity of the zero-order antisymmetric Lamb wave A 0 on the product of the wall thickness and frequency is substantially non-linear.
Определяют вторую разницу времен t2 между временем приема частотной составляющей антисимметричной волны Лэмба нулевого порядка A0 на второй выбранной частоте f2 и временем излучения по зависимости спектральной плотности мощности принятого сигнала от времени.The second time difference t 2 is determined between the time of receiving the frequency component of the zero-order antisymmetric Lamb wave A 0 at the second selected frequency f 2 and the radiation time according to the dependence of the power spectral density of the received signal on time.
Определяют дисперсию D[h] значений толщины стенки объекта по значениям второй выбранной частоты f2, второй разницы времен t2, расстояния L между преобразователями, среднеарифметическому значению толщины M[h], эталонным значениям групповой скорости антисимметричной волны Лэмба нулевого порядка A0 и второй производной ее обратной величины по произведению частоты и толщины, к примеру, по формулеDetermine the dispersion D [h] values of the second selected frequency values of the object wall thickness at f 2, the second difference time t 2, the distance L between the transducers, the arithmetic mean value of thickness M [h], the reference values of group antisymmetric wave velocity Lamb zero order A 0 and a second its inverse derivative of the product of frequency and thickness, for example, by the formula
Определяют минимальное значение толщины стенки объекта hmin по среднеарифметическому значению толщины M[h] и дисперсии значений толщины D[h], к примеру, по формулеThe minimum value of the object wall thickness h min is determined by the arithmetic mean value of the thickness M [h] and the variance of the thickness values D [h], for example, by the formula
Пример конкретного выполнения способа.An example of a specific implementation of the method.
Предлагаемый способ испытан на стальной пластине длиной L=1 м. На 50% длины пластины ее толщина равнялась 10 мм, на других 50% длины пластины ее толщина в результате воздействия коррозии была уменьшена до 8 мм. Таким образом, минимальное значение толщины составляло 8 мм, среднеарифметическое значение толщины равнялось 9 мм.The proposed method was tested on a steel plate with a length of L = 1 m. At 50% of the plate length, its thickness was 10 mm, on the other 50% of the plate length, its thickness was reduced to 8 mm due to corrosion. Thus, the minimum thickness value was 8 mm, the arithmetic mean value of the thickness was 9 mm.
В качестве акустических преобразователей для излучения и приема волн Лэмба использовались пьезоакустические преобразователи акустической эмиссии. Их установили на дальних краях пластины. Для излучения сигнала подали в заданный момент времени на электрический вход излучающего преобразователя импульсный электрический сигнал. Излучатель испустил импульсный акустический сигнал. Сигнал распространился по пластине в виде комбинации симметричной и антисимметричной волн Лэмба нулевого порядка S0 и A0, обладающих зависимостью групповой скорости от частоты и толщины стенки (фиг.1). Приемник принял сигнал U(t) (фиг.2). Из осциллограммы принятого сигнала U(t) путем преобразования Чои-Вильямса была получена зависимость спектральной плотности мощности принятого сигнала от времени W(f, t) (фиг.3).Piezoacoustic transducers of acoustic emission were used as acoustic transducers for the emission and reception of Lamb waves. They were installed at the far edges of the plate. To emit a signal, a pulsed electrical signal was applied to the electrical input of the emitting transducer at a given point in time. The emitter emitted a pulsed acoustic signal. The signal propagated through the plate in the form of a combination of symmetric and antisymmetric zero-order Lamb waves S 0 and A 0 , having a dependence of the group velocity on the frequency and wall thickness (Fig. 1). The receiver received a signal U (t) (figure 2). From the oscillogram of the received signal U (t) by the Choi-Williams transform, the dependence of the spectral power density of the received signal on the time W (f, t) was obtained (Fig. 3).
Определили первой значение частоты f1, в окрестности которой эталонная зависимость обратной величины групповой скорости симметричной волны Лэмба нулевого порядка S0 от произведения частоты и толщины близка к линейной зависимости. Для этого воспользовались тем фактом, что у объектов из стали подходящими являются значения первой частоты f1 из диапазона от 1.7 МГц·мм/М[h] до 2.4 МГц·мм/М[h] или из диапазона от 2.7 МГц·мм/М[h] до 3.4 МГц·мм/М[h], где M[h] - среднеарифметическое значение толщины (фиг.4). Исходя из начального значения толщины (10 мм) стенки, выбрали значение первой частоты f1, равное f1=230 кГц (фиг.4).The first value of the frequency f 1 was determined, in the vicinity of which the reference dependence of the reciprocal of the group velocity of the symmetric Lamb wave of zero order S 0 on the product of frequency and thickness is close to a linear dependence. For this, we took advantage of the fact that for steel objects, the values of the first frequency f 1 from the range from 1.7 MHz · mm / M [h] to 2.4 MHz · mm / M [h] or from the range from 2.7 MHz · mm / M are suitable [h] up to 3.4 MHz · mm / M [h], where M [h] is the arithmetic mean of the thickness (Fig. 4). Based on the initial value of the thickness (10 mm) of the wall, we chose the value of the first frequency f 1 equal to f 1 = 230 kHz (figure 4).
По зависимости спектральной плотности мощности принятого сигнала от времени W(f, t) определили первую разницу времен t1 между временем приема симметричной волны Лэмба нулевого порядка S0 на первой выбранной частоте f1 и временем излучения.Using the dependence of the spectral power density of the received signal on time W (f, t), we determined the first time difference t 1 between the time of reception of a symmetric Lamb wave of zero order S 0 at the first selected frequency f 1 and the radiation time.
Для этого провели определение времени приема следующим способом: выбрали уровень порога, равный 1% от максимального абсолютного значения зависимости спектральной плотности мощности принятого сигнала от времени, выявили на зависимости спектральной плотности мощности сигнала от времени W(f, t) точки, абсолютная величина спектральной плотности в которых превышает порог. Множество выявленных точек образовало на зависимости спектральной плотности 2 кривые характерной формы, соответствующие временам приема различных частотных составляющих симметричной S0 и антисимметричной A0 волн Лэмба нулевого порядка. Выявили точки, соответствующие симметричной волне Лэмба нулевого порядка S0 и первой выбранной частоте f1, выбрали из них точку, имеющую наибольшее соответствующее ей абсолютное значение спектральной плотности мощности принятого сигнала от времени. Соответствующее данной точке значение времени использовали в качестве времени приема симметричной волны Лэмба нулевого порядка S0 на первой выбранной частоте f1. В качестве времени излучения использовали время подачи управляющего электрического импульса на излучающий пьезоакустический преобразователь акустической эмиссии.To do this, we determined the reception time in the following way: we chose a threshold level equal to 1% of the maximum absolute value of the dependence of the power spectral density of the received signal on time, revealed the points, the absolute value of the spectral density on the dependence of the signal power spectral density on time W (f, t) in which exceeds the threshold. A number of identified points formed characteristic curves corresponding to the times of reception of various frequency components of the symmetric S 0 and antisymmetric A 0 zero-order Lamb waves in the dependence of spectral density 2. The points corresponding to the symmetric Lamb wave of zero order S 0 and the first selected frequency f 1 were identified, the point having the highest absolute value of the spectral power density of the received signal from time corresponding to it was selected. The time value corresponding to this point was used as the time of reception of a symmetric Lamb wave of zero order S 0 at the first selected frequency f 1 . As the radiation time, the time of supply of the control electric pulse to the radiating piezoelectric acoustic emission transducer was used.
Определенная указанным способом первая разница времен t1 оказалась равной t1=404.9 мкс.The first time difference t 1 determined by this method turned out to be t 1 = 404.9 μs.
По значению первой разницы времен t1 и по значению расстояния L между преобразователями определили значение групповой скорости СГР выбранной частотной составляющей симметричной волны Лэмба нулевого порядка S0 по формулеUsing the value of the first time difference t 1 and the value of the distance L between the transducers, we determined the group velocity С ГР of the selected frequency component of the symmetric Lamb wave of zero order S 0 according to the formula
По известным формулам рассчитали эталонные зависимости от произведения частоты f и толщины h значений групповой скорости
Подобрали значение толщины M[h], при котором эталонное значение групповой скорости
Для этого сначала по эталонной зависимости групповой скорости симметричной волны Лэмба нулевого порядка S0 от произведения частоты и толщины определили значение произведения толщины на частоту f1·M[h], при котором выполняется искомое равенство (фиг.6)For this, first, based on the reference dependence of the group velocity of a symmetric Lamb wave of zero order S 0 on the product of frequency and thickness, we determined the value of the product of the thickness by the frequency f 1 · M [h], at which the desired equality is satisfied (Fig.6)
f1·М[h]=2124.5 кГц·мм,f 1 · M [h] = 2124.5 kHz · mm,
затем определили само значение толщины M[h]then determined the very value of the thickness M [h]
M[h]=2124.5 кГц·м/f1=9.24 мм.M [h] = 2124.5 kHz · m / f 1 = 9.24 mm.
Полученное значение M[h] стали считать среднеарифметическим значением толщины стенки участка объекта между преобразователями.The obtained value of M [h] was considered the arithmetic mean value of the wall thickness of the site of the object between the transducers.
Определили значение второй частоты f2, в окрестности которой эталонная зависимость обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка A0 от произведения частоты и толщины существенно нелинейна. Для этого воспользовались тем фактом, что у объектов из стали подходящими являются значения второй частоты f2 из диапазона от 0 МГц·мм/М[h] до 1.0 МГц·мм/М[h] (фиг.4). Исходя из начального значения толщины (10 мм) стенки, выбрали значение второй частоты, равное f2=25 кГц (фиг.4).The value of the second frequency f 2 was determined, in the vicinity of which the reference dependence of the reciprocal of the group velocity of the zero-order antisymmetric Lamb wave A 0 from the product of frequency and thickness is substantially non-linear. For this, we took advantage of the fact that for steel objects, the values of the second frequency f 2 from the range from 0 MHz · mm / M [h] to 1.0 MHz · mm / M [h] are suitable (Fig. 4). Based on the initial value of the thickness (10 mm) of the wall, we chose the value of the second frequency equal to f 2 = 25 kHz (figure 4).
По зависимости спектральной плотности мощности принятого сигнала от времени W(f, t) определили вторую разницу времен t2 между временем приема антисимметричной волны Лэмба нулевого порядка A0 на второй выбранной частоте f2 и временем излучения t2=424.4 мкс (фиг.4).From the dependence of the spectral power density of the received signal on time W (f, t), we determined the second time difference t 2 between the time of reception of the antisymmetric zero-order Lamb wave A 0 at the second selected frequency f 2 and the radiation time t 2 = 424.4 μs (Fig. 4) .
Взяли значение второй производной
Определили дисперсию D[h] значений толщины стенки объекта по значениям второй выбранной частоты f2, второй разницы времен t2, расстояния L между преобразователями, среднеарифметическому значению толщины M[h], эталонным значениям групповой скорости антисимметричной волны Лэмба нулевого порядка A0 и второй производной ее обратной величины по произведению частоты и толщины по формулеIdentified dispersion D [h] values of the second selected frequency values of the object wall thickness at f 2, the second difference time t 2, the distance L between the transducers, the arithmetic mean value of thickness M [h], the reference values of group antisymmetric wave velocity Lamb zero order A 0 and a second the derivative of its reciprocal of the product of frequency and thickness according to the formula
Определили минимальное значение толщины стенки hmin объекта на участке между преобразователями по среднеарифметическому значению толщины M[h] на участке между преобразователями и дисперсии значений толщины D[h] на участке между преобразователями по формулеThe minimum value of the wall thickness h min of the object in the area between the transducers was determined by the arithmetic mean of the thickness M [h] in the area between the transducers and the variance of the thicknesses D [h] in the area between the transducers by the formula
Таким образом, предлагаемый способ является высокоэффективным и может при широком внедрении принести значительный экономический эффект за счетThus, the proposed method is highly effective and can, with widespread implementation, bring significant economic benefits due to
снижения трудоемкости, упрощения и ускорения контроля толщины больших по площади объектов,reduce the complexity, simplify and speed up the control of the thickness of large objects,
устранения необходимости получения физического доступа ко всей поверхности объекта, подлежащей контролю толщины,eliminating the need for physical access to the entire surface of the object subject to thickness control,
обеспечения возможности контроля толщины при наличии вариаций значений толщины, к примеру при наличии на объекте очаговой коррозии,providing the ability to control the thickness in the presence of variations in the thickness values, for example, if there is focal corrosion at the facility,
обеспечения возможности определения среднеарифметического и минимального значений толщины на участке между двумя акустическими преобразователями, установленными на объекте в произвольных точках, без проведения предварительного измерения толщины в тех же точках.providing the ability to determine the arithmetic mean and minimum thickness values in the area between two acoustic transducers installed on the object at arbitrary points, without preliminary thickness measurement at the same points.
Математический расчет.Mathematical calculation.
Теоретические исследования и анализ литературных источников показали, что в стенке объекта акустический сигнал на частотах, на которых длины объемных волн сравнимы или превосходят толщину стенки, распространяется в виде волн Лэмба. Значение групповой скорости СГР волны Лэмба зависит (фиг.1) от произведения частоты f и толщины стенки объекта h, а также от значений скоростей объемной продольной CL и объемной поперечной CT волн в материале стенки, типа G волны Лэмба (симметричная S или антисимметричная A) и порядкового номера i волны ЛэмбаTheoretical studies and analysis of literature showed that in the wall of an object an acoustic signal at frequencies at which body wavelengths are comparable or greater than the wall thickness propagates in the form of Lamb waves. The value of the group velocity C of the GR of the Lamb wave depends (Fig. 1) on the product of the frequency f and the wall thickness of the object h, as well as on the values of the velocities of the longitudinal longitudinal C L and volumetric transverse C T waves in the wall material, type G Lamb waves (symmetric S or antisymmetric A) and serial number i of the Lamb wave
Значение групповой скорости
Известно, что излучение акустического импульса эквивалентно синхронному излучению множества частотных составляющих. В этом случае различные частотные составляющие различных волн Лэмба приходят (фиг.2) на приемный преобразователь в различные моменты времени
где f - частота, x - координата вдоль прямой линии, соединяющей преобразователи, начало координат соответствует точке, в которой установлен излучающий преобразователь, h(x) - зависимость толщины стенки объекта от координаты x, L - расстояние между преобразователями, время
Расчет зависимости спектральной плотности мощности принятого сигнала от времени (фиг.3) позволяет определить время приема
Разложение функции
Следовательно, время приема равноTherefore, the reception time is
Поскольку , тоInsofar as then
Следовательно,Hence,
где величинаwhere is the value
является дисперсией значений толщины стенки объекта на участке между преобразователями, величина M[h2] равна среднеарифметическому значению квадрата толщины стенки
Как правило, на контролируемых объектах дисперсия значений толщины стенки мала по сравнению с квадратом толщины стенкиAs a rule, on controlled objects, the variance of the wall thickness is small compared to the square of the wall thickness
Если при некотором значении частоты f1, значении толщины, равном среднеарифметическому значению толщины стенки M[h], типе волны Лэмба G и порядковом номере волны Лэмба i справедливо, что зависимость значений обратной величины эталонной групповой скорости волны Лэмба от произведения толщины и частоты близка к линейной или, другими словами, значение производной второго порядка обратной величины эталонной групповой скорости волны Лэмба по произведению толщины и частоты
В этом случае выражение (1) упрощается доIn this case, expression (1) is simplified to
Определение по известным значениям
Если при некотором значении частоты f2, значении толщины, равном среднеарифметическому значению толщины стенки M[h], типе волны Лэмба G′ и порядковом номере волны Лэмба j справедливо, что зависимость значений обратной величины эталонной групповой скорости волны Лэмба от произведения толщины и частоты существенно нелинейна или, другими словами, значение производной второго порядка обратной величины эталонной групповой скорости волны Лэмба по произведению толщины и частоты
Из (4) следует, что имея среднеарифметическое значение толщины стенки M[h], значения второй частоты f2, времени
Наличие среднеарифметического значения толщины M[h] и дисперсии значений толщины D[h] позволяет установить минимальное значение толщины стенки объекта hmin, к примеру, по формулеThe presence of the arithmetic mean of the thickness M [h] and the variance of the thickness D [h] allows you to set the minimum value of the wall thickness of the object h min , for example, by the formula
Claims (1)
размещение на его поверхности на известном расстоянии друг от друга, по крайней мере, одного акустического преобразователя для излучения волн Лэмба и, по крайней мере, одного преобразователя для их приема,
излучение в заданный момент времени импульсного сигнала,
расчет зависимости спектральной плотности мощности принятого сигнала от времени,
выбор волны Лэмба и частоты,
определение разности между временем приема выбранной частотной составляющей выбранной волны Лэмба и временем излучения сигнала,
определение значения групповой скорости выбранной частотной составляющей выбранной волны Лэмба по известному значению расстояния между преобразователями и значению разности между временем приема выбранной частотной составляющей выбранной волны Лэмба и временем излучения сигнала,
определение значения толщины стенки по полученному значению групповой скорости, выбранному значению частоты и эталонной зависимости групповой скорости выбранной волны Лэмба от произведения толщины стенки и частоты,
отличающийся тем, что
устанавливают минимальную величину толщины стенки
по полученным среднеарифметической величине толщины стенки и дисперсии значений толщины стенки объекта, причем
среднеарифметическое значение толщины стенки определяют с выбором симметричной волны Лэмба нулевого порядка и, по крайней мере, одной частоты, на которой эталонная зависимость обратной величины групповой скорости симметричной волны Лэмба нулевого порядка от произведения толщины стенки и частоты близка к линейной, а
дисперсию значений толщины стенки определяют с выбором антисимметричной волны Лэмба нулевого порядка и частоты, на которой эталонная зависимость обратной величины групповой скорости антисимметричной волны Лэмба нулевого порядка от произведения толщины стенки и частоты существенно нелинейна. A method of monitoring the dynamics of changes in the wall thickness of a controlled object, including
placing on its surface at a known distance from each other at least one acoustic transducer for emitting Lamb waves and at least one transducer for receiving them,
radiation at a given point in time of a pulse signal,
calculation of the dependence of the spectral power density of the received signal on time,
Lamb wave and frequency selection,
determination of the difference between the reception time of the selected frequency component of the selected Lamb wave and the signal emission time,
determining the group velocity value of the selected frequency component of the selected Lamb wave from the known distance between the transducers and the difference between the reception time of the selected frequency component of the selected Lamb wave and the signal emission time,
determining the wall thickness value from the obtained group velocity value, the selected frequency value and the reference dependence of the group velocity of the selected Lamb wave on the product of the wall thickness and frequency,
characterized in that
set the minimum wall thickness
according to the arithmetic mean wall thickness and variance of the object wall thickness values, and
the arithmetic mean value of the wall thickness is determined with the choice of a symmetric Lamb wave of zero order and at least one frequency at which the reference dependence of the inverse of the group velocity of a symmetric Lamb wave of zero order on the product of the wall thickness and frequency is close to linear, and
the variance of the wall thickness values is determined with the choice of the zero-order antisymmetric Lamb wave and frequency, at which the reference dependence of the inverse group velocity of the zero-order antisymmetric Lamb wave on the product of the wall thickness and frequency is substantially non-linear.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013149631/28A RU2540942C1 (en) | 2013-11-07 | 2013-11-07 | Method to monitor dynamics of variation of controlled object wall thickness |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2013149631/28A RU2540942C1 (en) | 2013-11-07 | 2013-11-07 | Method to monitor dynamics of variation of controlled object wall thickness |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2540942C1 true RU2540942C1 (en) | 2015-02-10 |
Family
ID=53287020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2013149631/28A RU2540942C1 (en) | 2013-11-07 | 2013-11-07 | Method to monitor dynamics of variation of controlled object wall thickness |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2540942C1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007308A (en) * | 2017-12-04 | 2019-07-12 | 新日本无线株式会社 | Ultrasonic sensor and vehicle control system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4372163A (en) * | 1981-02-03 | 1983-02-08 | Rockwell International Corporation | Acoustic measurement of near surface property gradients |
DE3611967A1 (en) * | 1985-04-17 | 1986-10-30 | Rolls-Royce Ltd., London | DETERMINATION OF THE SPECTRAL CONTENT OF HIKING WAVE STRESSES |
RU2153163C1 (en) * | 1999-11-29 | 2000-07-20 | Долгих Владимир Иванович | Method of intratube ultrasonic diagnostics of condition of pipe-line |
RU2177612C2 (en) * | 1999-09-01 | 2001-12-27 | Пермикин Владимир Сергеевич | Method of ultrasonic inspection of state of metal working under conditions of creepage forecasting its residual resource and acoustic unit for its implementation |
US6360609B1 (en) * | 2000-02-23 | 2002-03-26 | Massachusetts Institute Of Technology | Method and system for interpreting and utilizing multimode dispersive acoustic guided waves |
RU2480742C1 (en) * | 2011-10-10 | 2013-04-27 | Федеральное Государственное Унитарное Предприятие "Научно-Производственное Объединение "Техномаш" | Method for acoustic emission control |
-
2013
- 2013-11-07 RU RU2013149631/28A patent/RU2540942C1/en not_active IP Right Cessation
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4372163A (en) * | 1981-02-03 | 1983-02-08 | Rockwell International Corporation | Acoustic measurement of near surface property gradients |
DE3611967A1 (en) * | 1985-04-17 | 1986-10-30 | Rolls-Royce Ltd., London | DETERMINATION OF THE SPECTRAL CONTENT OF HIKING WAVE STRESSES |
RU2177612C2 (en) * | 1999-09-01 | 2001-12-27 | Пермикин Владимир Сергеевич | Method of ultrasonic inspection of state of metal working under conditions of creepage forecasting its residual resource and acoustic unit for its implementation |
RU2153163C1 (en) * | 1999-11-29 | 2000-07-20 | Долгих Владимир Иванович | Method of intratube ultrasonic diagnostics of condition of pipe-line |
US6360609B1 (en) * | 2000-02-23 | 2002-03-26 | Massachusetts Institute Of Technology | Method and system for interpreting and utilizing multimode dispersive acoustic guided waves |
RU2480742C1 (en) * | 2011-10-10 | 2013-04-27 | Федеральное Государственное Унитарное Предприятие "Научно-Производственное Объединение "Техномаш" | Method for acoustic emission control |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110007308A (en) * | 2017-12-04 | 2019-07-12 | 新日本无线株式会社 | Ultrasonic sensor and vehicle control system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2479119C (en) | Self calibrating apparatus and method for ultrasonic determination of fluid properties | |
JP6438957B2 (en) | System and method for monitoring defects | |
US6634233B2 (en) | Method for determining the wall thickness and the speed of sound in a tube from reflected and transmitted ultrasound pulses | |
EP2799820B1 (en) | Liquid surface level measurement device, method, and program | |
US10458871B2 (en) | Apparatus and method for measuring the pressure inside a pipe or container | |
Rahiman et al. | Design and modelling of ultrasonic tomography for two-component high-acoustic impedance mixture | |
US20100131246A1 (en) | Ultrasonic surface monitoring | |
US9383238B2 (en) | Apparatus, system and process for characterizing multiphase fluids in a fluid flow stream | |
US20210293947A1 (en) | Continuous wave ultrasound or acoustic non-destructive testing | |
RU2540942C1 (en) | Method to monitor dynamics of variation of controlled object wall thickness | |
RU2661455C1 (en) | Method for determining the viscoelastic properties of liquid and solid media and the device for its implementation | |
JP5059344B2 (en) | Plate thickness measuring apparatus and measuring method | |
JP4904099B2 (en) | Pulse signal propagation time measurement device and ultrasonic flow measurement device | |
RU2687086C1 (en) | Method of ultrasonic monitoring of pipeline wall thickness | |
RU2545065C2 (en) | Method to measure acoustic speed in water | |
Pal | Fourier transform ultrasound spectroscopy for the determination of wave propagation parameters | |
JP4704536B2 (en) | Liquid level measuring device in pipe and liquid level measuring method | |
RU2437066C1 (en) | Method for ultrasonic measurement of level of liquid in reservoirs and apparatus for ultrasonic measurement of level of liquid in reservoirs | |
Liu et al. | Acoustic method for obtaining the pressure reflection coefficient using a half-wave layer | |
JP4674007B2 (en) | Liquid level measuring device in pipe and liquid level measuring method | |
RU2570097C1 (en) | Method of ultrasonic echo-pulse thickness gauging | |
RU2425362C2 (en) | Method of determining location of acoustic emission sources using one receiver | |
RU2313068C2 (en) | Mode of measuring gas consumption in main pipelines and an arrangement for its execution | |
RU2589751C2 (en) | Method of determining the average grain diameter of metal products and device for its implementation | |
RU2586087C1 (en) | Method of locating defects during acoustic-emission control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20181108 |