Научная статья на тему 'Трансформация и рассеяние поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница акустического контакта'

Трансформация и рассеяние поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница акустического контакта Текст научной статьи по специальности «Нанотехнологии»

CC BY
237
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОВЕРХНОСТНАЯ ВОЛНА / ВОЛНА РЭЛЕЯ / ВОЛНА СТОУНЛИ / ТРАНСФОРМАЦИЯ / ОТРАЖЕНИЕ

Аннотация научной статьи по нанотехнологиям, автор научной работы — Баев А.Р., Пантелеенко Ф.И., Захаренко В.В., Размыслович Г.И., Жаворонков К.Г.

Эффекты трансформации и распространения волн Рэлея и Стоунли представляют значительный интерес для расширения возможностей и совершенствования методов ультразвукового контроля и измерений. Цель данной работы заключалась в установлении влияния геометрических параметров тела акустической нагрузки и его положения на коэффициенты отражения и прохождения волны Стоунли и Рэлея и выявлении возможности использования результатов исследования для практических приложений.На основе анализа акустического тракта и данных эксперимента установлена связь между измеряемыми амплитудными параметрами и коэффициентами прохождения и отражения поверхностных волн, а также отражательной способностью области контакта тела нагрузки в виде призмы через скользящую границу, которая достигает ≈ 32-34 дБ. Впервые определены зависимости указанных коэффициентов от угла наклона одной из боковых граней призмы в диапазоне 0 ± 45°, безразмерной толщины контактной прослойки (0-0,05) и ее ориентации относительно акустической оси.Установлено, что эти коэффициенты преимущественно максимальны, когда призма прямоугольная. Коэффициент же отражательной способности при жестком контакте тел более чем на порядок меньше, а коэффициенты прохождения сравнимы по величине. Показана перспективность использования результатов исследования для контроля качества сцепления материалов при сварке, пайке, склейке, выявления дефектов в труднодоступных местах, а также для определения физико-механических свойств металлов с помощью предложенного способа создания опорного сигнала.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по нанотехнологиям , автор научной работы — Баев А.Р., Пантелеенко Ф.И., Захаренко В.В., Размыслович Г.И., Жаворонков К.Г.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Transformation and scattering of surface waves on the acoustic load to ultrasonic evaluation and measurements. Part 1. The boundary of acoustic contact is sliding

The effects of Rayleigh and Stoneley wave transformation and propagation are of considerable interest for the enhancement and improvement of ultrasonic methods evaluation and measurements. The purpose of this work is to determine the influence of the geometric parameters of the acoustic load body and its position on the coefficients of reflection and propagation of the Stoneley and Rayleigh waves and to identify the possibility of using the results of the study for practical applications.Based on the analysis of the acoustic path and the experimental data, the relationship between the measured amplitude parameters and the coefficients of the propagation and reflection of surface waves, as well as the reflectivity of the contact region of the load body in the form of a prism through the sliding boundary, which reaches up to ≈ 32-34 дБ, is established. For the first time, the dependence of these coefficients on the inclination angle of one of the prism lateral faces in the range of 0 ± 45°, dimensionless thickness of the contact layer (0-0,05) and its orientation relative to the acoustic axis.It is established that these coefficients are mainly maximal when the prism is rectangular. The coefficient of reflectivity in the hard contact of bodies is more than an order of magnitude less, and the coefficients of wave propagation comparable in magnitude. The prospects of using the results of the study to evaluate the quality of adhesion of materials during welding, soldering, gluing, detection of defects in hard-to-reach places, as well as to determine the physical and mechanical properties of metals by the proposed method of creating a reference signal are shown.

Текст научной работы на тему «Трансформация и рассеяние поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница акустического контакта»

УДК 620.179.161. 534.1

Трансформация и рассеяние поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница

акустического контакта

Баев А.Р.1, Пантелеенко Ф.И.2, Захаренко В.В.2, Размыслович Г.И.1, Жаворонков К.Г.1, Гиль Н.Н.1

'Институт прикладной физики Национальной академии наук Беларуси, ул. Академическая, 16, г. Минск 220072, Беларусь 2Белорусский национальный технический университет, пр-т Независимости, 65, г. Минск 220013, Беларусь

Поступила 09.01.2018 Принята к печати 27.02.2018

Эффекты трансформации и распространения волн Рэлея и Стоунли представляют значительный интерес для расширения возможностей и совершенствования методов ультразвукового контроля и измерений. Цель данной работы заключалась в установлении влияния геометрических параметров тела акустической нагрузки и его положения на коэффициенты отражения и прохождения волны Стоунли и Рэлея и выявлении возможности использования результатов исследования для практических приложений.

На основе анализа акустического тракта и данных эксперимента установлена связь между измеряемыми амплитудными параметрами и коэффициентами прохождения и отражения поверхностных волн, а также отражательной способностью области контакта тела нагрузки в виде призмы через скользящую границу, которая достигает ~ 32-34 дБ. Впервые определены зависимости указанных коэффициентов от угла наклона одной из боковых граней призмы в диапазоне 0 ± 45°, безразмерной толщины контактной прослойки (0-0,05) и ее ориентации относительно акустической оси.

Установлено, что эти коэффициенты преимущественно максимальны, когда призма прямоугольная. Коэффициент же отражательной способности при жестком контакте тел более чем на порядок меньше, а коэффициенты прохождения сравнимы по величине. Показана перспективность использования результатов исследования для контроля качества сцепления материалов при сварке, пайке, склейке, выявления дефектов в труднодоступных местах, а также для определения физико-механических свойств металлов с помощью предложенного способа создания опорного сигнала.

Ключевые слова: поверхностная волна, волна Рэлея, волна Стоунли, трансформация, отражение.

DOI: 10.21122/2220-9506-2018-9-1-28-39

Адрес для переписки:

Баев А.Р.

Ин-т прикладной физики Национальной академии наук Беларуси, ул. Академическая, 16, г. Минск 220072, Беларусь e-mail: baev@iaph.bas-net.by

Для цитирования:

Баев А.Р., Пантелеенко Ф.И., Захаренко В.В., Размыслович Г.И., Жаворонков К.Г., Гиль Н.Н.

Трансформация и рассеяния поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница акустического контакта. Приборы и методы измерений. 2018. - Т. 9, № 1. С. 28-39. БОТ: 10.21122/2220-9506-2018-9-1-28-39

Address for correspondence: Baev A.R.

Institute of Applied Physics of the National Academy of Science of Belarus,

Akademicheskaya str., 16, Minsk 220072, Belarus e-mail: baev@iaph.bas-net.by For citation:

Baev A.R., Panteleenko F.I., Zakharenko V.V., Razmyslovich G.I., Lark K.G., Gill N.N.

[Transformation and scattering of surface waves on the acoustic load

to ultrasonic evaluation and measurements. Part 1. The boundary

of acoustic contact is sliding].

Devices and Methods of Measurements.

2018, vol. 9, no. 1, pp. 28-39 (in Russian).

DOI: 10.21122/2220-9506-2018-9-1-28-39

Transformation and scattering of surface waves on the acoustic load to ultrasonic evaluation and measurements. Part 1. The boundary of acoustic contact is sliding

Baev A.R.1, Panteleenko F.I.2, Zakharenko V.V.2, Razmyslovich G.I.1, Lark K.G.1, Gill N.N.1

1Institute of Applied Physics of the National Academy of Science of Belarus, Akademicheskaya str., 16, Minsk 220072, Belarus 2Belarusian National Technical University, NezavisimostyAve., 65, Minsk 220013, Belarus

Received 09.01.2018

Accepted for publication 27.02.2018

Abstract

The effects of Rayleigh and Stoneley wave transformation and propagation are of considerable interest for the enhancement and improvement of ultrasonic methods evaluation and measurements. The purpose of this work is to determine the influence of the geometric parameters of the acoustic load body and its position on the coefficients of reflection and propagation of the Stoneley and Rayleigh waves and to identify the possibility of using the results of the study for practical applications.

Based on the analysis of the acoustic path and the experimental data, the relationship between the measured amplitude parameters and the coefficients of the propagation and reflection of surface waves, as well as the reflectivity of the contact region of the load body in the form of a prism through the sliding boundary, which reaches up to « 32-34 gE, is established. For the first time, the dependence of these coefficients on the inclination angle of one of the prism lateral faces in the range of 0 ± 45°, dimensionless thickness of the contact layer (0-0,05) and its orientation relative to the acoustic axis.

It is established that these coefficients are mainly maximal when the prism is rectangular. The coefficient of reflectivity in the hard contact of bodies is more than an order of magnitude less, and the coefficients of wave propagation - comparable in magnitude. The prospects of using the results of the study to evaluate the quality of adhesion of materials during welding, soldering, gluing, detection of defects in hard-to-reach places, as well as to determine the physical and mechanical properties of metals by the proposed method of creating a reference signal are shown.

Keywords: surface wave, Rayleigh wave, Stoneley wave, transformation, reflection. DOI: 10.21122/2220-9506-2018-9-1-28-39

Адрес для переписки: Баев А.Р.

Ин-т прикладной физики Национальной академии наук Беларуси, ул. Академическая, 16, г. Минск 220072, Беларусь e-mail: baev@iaph.bas-net.by

Для цитирования:

Баев А.Р., Пантелеенко Ф.И., Захаренко В.В., Размыслович Г.И., Жаворонков К.Г., Гиль Н.Н.

Трансформация и рассеяния поверхностных волн на акустической нагрузке для ультразвукового контроля и измерений. Ч. 1. Скользящая граница акустического контакта. Приборы и методы измерений. 2018. - Т. 9, № 1. С. 28-39. БОТ: 10.21122/2220-9506-2018-9-1-28-39

Address for correspondence: Baev A.R.

Institute of Applied Physics of the National Academy of Science of Belarus,

Akademicheskaya str., 16, Minsk 220072, Belarus e-mail: baev@iaph.bas-net.by For citation:

Baev A.R., Panteleenko F.I., Zakharenko V.V, Razmyslovich G.I., Lark K.G., Gill N.N.

[Transformation and scattering of surface waves on the acoustic load

to ultrasonic evaluation and measurements. Part 1. The boundary

of acoustic contact is sliding].

Devices and Methods of Measurements.

2018, vol. 9, no. 1, pp. 28-39 (in Russian).

DOI: 10.21122/2220-9506-2018-9-1-28-39

Введение

Все более широкое распространение в области неразрушающего контроля и ультразвуковых измерений находят поверхностные волны (ПАВ), включая волны Рэлея (ВР), Стоунли (ВС) и др., относящиеся к классу слабозатухающих волн. Энергия этих волн локализована в поверхностном слое (слоях) исследуемого материала глубиной h, не превышающей их длины волны ^ддз. Амплитудные, фазовые, частотно-фазовые параметры таких ПАВ весьма чувствительны к изменению структуры исследуемого объекта, наличию поверхностных и подповерхностных несплошностей, качеству сцепления контактирующих материалов и др. [1-7]. При определенных условиях изменение граничных условий в области распространения ВР сопровождается появлением краевых волн, обусловленных трансформацией ПАВ в волны Стоунли, а также в объемные моды. В частности, в работе [4] упоминалось, что изменение граничных условий на пути распространения ВР обусловлено наличием на объекте прямоугольного выступа, что позволяет (в качестве физической модели) рассматривать объект как двухслойную среду с жестким акустическим контактом. Причем область в окрестности линии взаимного пересечения контактной поверхности объекта и передней грани выступа с пространственным параметром не более I ~ ^ПДВ служит источником трансформированных и рассеянных в объеме выступа краевых волн.

Отметим, что возбуждение ВР в однородном материале объекта со свободной границей и скоростью СПДВ осуществляется при падении объемной волны, имеющей скорость С на поверхность объекта под углом в = агат^/С^^). Возбуждение же ВС на практике реализуется, как правило, путем трансформации ВР на акустической нагрузке, представляющей собой твердое тело, контактирующее с материалом подложки. Обобщенное граничное условие в области поверхности контакта двух тел г = 0, используемое для моделирования процессов распространения ВС (двумерная постановка), предложено в работах [8, 9] для компонент смещений п(,) и напряжений {о(,) , о(,) , о(,) }:

Х7' 77 7 ХУУ

а(1) - а(2) ; а(1) -а (2)1

lz=0

lz=0

lz=0

lz=0

(1)

где Еt и Еп - коэффициенты, характеризующие эффективность передачи, соответственно, тангенциальной и нормальной составляющей напряжений через границу сред.

Полученное на основе решения 4 волновых уравнений относительно скалярного и вектора потенциала характеристическое уравнение вырождается в уравнение, описывающее распространение ВС при наличии «жесткого» контакта двух тел (Е, Еи^ да) или «скользящего» контакта (Е, ^ 0, ЕпФ 0). Согласно общей методике для определения вещественного корня волнового числа k необходимо выполнение неравенства kt < kti или Ся < С где индекс Т относится к поперечной моде. При реализации соответствующих граничных условий необходимое и достаточное условие существования ВС, заключается в выполнении неравенства:

ДП. 6, Е, , Е„ g, ...)|, > 1, (2)

Ч I М 1 Ь =( &

g - отно-

Ц]

:мл )sJ 0 ; и7\ = uz

x |r=o x ' z=0 z lz=0 z

„(Ц _„(2)

+ CT:

.(2) д

где -

шение коэффициентов Пуассона контактирующих материалов.

Согласно (2) моды ВС существуют и при контакте одинаковых тел через скользящую границу, т.е. при выполнении равенства: 91 = 92 = 93 = g = /ц1 = р2 /р1 =1.

Однако в случае жесткой границы контакта двух одинаковых тел ПАВ трансформируется в объемные моды. При этом в работе [4] примером жесткого контакта тела акустической нагрузки (ТАН) с подложкой служил металлический образец с выступом, плоскость передней грани которого образует с контактной поверхностью угол у = 90°. Там же показано, что превалирующая часть энергии краевых волн, генерируемой непосредственно в области сопряжения указанных поверхностей, является энергией вертикально поляризованной поперечной моды.

Необходимо отметить, что в ряде технологических процессов соединения материалов (склейка, сварка, пайка и др.) граничные условия в области контакта материалов во времени претерпевают изменения - от скользящей до жесткой границы. Причем область в окрестности изменения или разрыва граничных условий можно рассматривать как самостоятельный вторичный источник излучения объемных волн, волн Стоунли и др. [5, 10]. С одной стороны, эти волны несут информацию о состоянии границы раз-

дела сред, а с другой - могут быть использованы для выявления дефектов в труднодоступных для традиционных методик местах. Отраженные от этой области волны могут создавать шумовой фон, препятствующий проведению контроля такого объекта, но в то же время служить для создания опорного сигнала при проведении измерений физико-механических свой твердых тел глубины упрочненных слоев по данным скорости, амплитуды или спектра акустического сигнала ПАВ и др.

Таким образом, установление закономерностей трансформации ПАВ в краевые волны на ТАН, включая ВС, объемные и др., представляет как научный, так и практический интерес. Знание параметров акустического тракта таких систем весьма важны при выполнении методической части работы с целью разработки новых методик ультразвуковой диагностики свойств материалов по данным скорости нескольких мод и их спектральным характеристикам. Попытки аналитического описания указанных процессов трансформации ПАВ на акустической нагрузке наталкиваются на ряд трудностей, связанных с выбором корректной модели механизма преобразования смещений Ъ>. и напряжений ok в области условий проявления сингулярности передаточной функции. В этом случае несомненна роль экспериментального подхода.

В известных же преимущественно экспериментальных исследованиях основное внимание обращено на влияние состояния именно качества контакта сопрягаемых поверхностей при склейке, приложении силы давления. На фундаментальные же зависимости параметров акустического тракта от геометрии тела ТАН не обращалось должного внимания. К ним следует отнести коэффициенты прохождения ПАВ по амплитуде с последующе трансформацией в другую моду согласно схеме ВР^ВС (коэффициент DRS) и ВС^ВР (Dsr), а также коэффициенты отражения волны Стоунли RS и волны Рэлея Rr от соответствующих границ контакта ТАН с подложкой (рисунок 1). Эти вопросы представляют интерес с точки зрения как развития теоретической модели трансформации упругих мод на неоднородной границе, так и разработки усовершенствованных методик и средств ультразвуковых измерений и контроля.

Цель настоящей работы заключалась в разработке методики экспериментального исследования и установлении влияния геометрических параметров ТАН на ее отражательную способ-

ность и коэффициенты отражения, прохождения волны Стоунли и Рэлея, а также выявлении возможности использования результатов исследования для практических приложений, включая создание опорного сигнала при ультразвуковых измерениях и контроле поверхности и структуры материалов, что позволяет существенно повысить надежность и расширить возможности не-разрушающего контроля.

Position I Position II

положение L положение II

Рисунок 1 - Экспериментальная схема исследования отражения и прохождения поверхностных волн по подложке при наличии тела акустической нагрузки в виде призмы: 1 - излучающий пьезоэлектрический преобразователь; 2 - приемный пьезоэлектрический преобразователь; 3 - тело акустической нагрузки;

4 - подложка объекта; 5 - боковые грани призмы; ф - угол наклона боковых граней призмы

Figure 1 - Experimental scheme of investigation of reflection and propagation of surface waves on the basic solid in the presence of the body of the acoustic load in the form of a prism: 1 - emitting piezoelectric probe; 2 - receiving piezoelectric probe; 3 - the body of the acoustic load as a prism; 4 - the basic solid of the object;

5 - the lateral faces of the prism; ф - the angle of the lateral faces of the prism

Краткий анализ акустического тракта применительно к методике исследования

В качестве исследуемого объекта ТАН выбраны располагаемые в области прохождения ПАВ тела с наиболее простой геометрией и представляющие собой стальные призмы, контактирующие со стальной подложкой через контактную смазку, удельное акустическое сопротивление которой Яж = (рС^)ж более чем в 20 раз меньше импеданса контактирующих через нее стальных материалов RT = (рС^ , где индекс L относится к продольной, а Т - поперечной моде (рисунок 1). Эффекты трансформации и рассеяния ПАВ имеют место в окрестности линии пересечения плоскости боковых граней призмы с контактной поверхностью (граница 1 и 2),

т.е. в области разрыва граничных условий. Причем одна из этих граней наклонена под углом к вертикали в варьируемом диапазоне ф = 0 ± 45°, а оппозитная ей грань образует угол ф = 0 с вертикалью. Рассматриваются два положения призмы относительно направления акустической оси, вдоль которой распространяются ПАВ, излучаемые их источником - пьезоэлектрическими преобразователем (ПЭП).

Проанализируем в двумерной постановке акустический тракт применительно к определению указанных выше коэффициентов прохождения и отражения ПАВ, а также введенного коэффициента Ж, характеризующего отражательную способность ТАН. Последний определяется отношением показаний акустических сигналов или амплитуды отраженной ПАВ от 1-й, обращенной к источнику ПАВ границы (А^, и оппозитной ей 2-й границы (А2) (рисунок 1). При этом, как легко показать:

R = A2/A] =fE (Drs)i (Dsr)i (Rs)2 (RrV\

(3)

где индексы 1 и 2 соответствуют указанным границам контакта тела ТАН с подложкой, а/Е - поправочный коэффициент, учитывающий влияние расхождения и диссипацию ПАВ.

Если же прозвучивание объекта производится теневым способом, а Аго и АТ соответствуют амплитудам опорного сигнала в отсутствие ТАН и при наличии последней соответственно, то нормализованный коэффициент прохождения ПАВ через нагрузку примет вид:

NT* AJAm fT (DrJi (Dsr)2

(4)

где /т - корректирующий коэффициент, подобный /

Положим, что измеренные в режиме эхо амплитуды отраженных сигналов от 2-й границы ТАН в отсутствие ее наклона А2Е0 (ф = 0) и с наклоном А2Е (ф Ф 0), а измеренные в теневом режиме - амплитуды Ат0 и Ат соответственно. Тогда из (3) и (4) следует, что нормализованные коэффициенты отражения волны Стоунли а также звукопрозрачности Оя* и Бж*, полученные при развороте призмы на 180°, имеют вид:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

RS* = - (RS)2/RS0 = A2E/A2

(5)

^ = (DrsVDRSO = (AAA; Dr1* = (DSRVDSRo = (A/A™V (6)

Вследствие нарушения однородности граничных условий можно считать, что области, х с ± Ах, 2 с ^ ±Az, где {Аг/ХПдВ, Лх/А^} < 1 соответствуют расположению вторичных источников краевых волн, где х1 соответствует отражающей волну Рэлея границе 1, а х2 - границе, отражающей волну Стоунли и трансформированной из ВР на границе 1. Так что закон сохранения энергии потоков акустических мод в двумерной постановке и пренебрежении диссипативными процессами представляется в виде:

П = (Пт + П + П,) (1 - - (7)

гдеП = С| wdS - поток энергии той или иной

Sr

моды через поверхность граничной области S, V! - плотность этого потока; Пт и П£ - поток ТАН и энергии рассеянных на границах краевых поперечных и продольных волн; Б(0)д и Л(0)д - коэффициенты прохождения потока энергии ПАВ на боковые поверхности и отражения его от зоны разрыва граничных условий соответственно.

Учитывая полученные выше формулы (2-6), устанавливающие связь между упомянутыми выше параметрами акустического тракта и показаниями измеряемых амплитуд акустического сигнала ПАВ, были проведены экспериментальные исследования по определению зависимостей {Ж, Б^*, Бж*} от угла ф, разворота призмы на 180°, а также изменения толщины контактного слоя h.

Методические особенности эксперимента

На рисунке 1 представлена схема экспериментальных исследований. Моделирование граничных условий типа «скользящая граница» осуществляется путем размещения на стальной подложке стального образца в виде призмы через тонкий слой контактной смазки. Проводилась тщательная шлифовка контактирующих поверхностей, так что оцененная толщина контактной смазки h (минерального масла) при установке призмы в рабочее положение под постоянным усилием, создаваемым постоянными магнитами, не превышала 5 мкм. Т.е. выполнялось условие: к = к/А < 0,004, или к/А < 0,01, где А и А - дли-

А ж о ж о

на продольной волны в жидкости и объекте соответственно. Размеры контактной поверхности призм (8 х 19)-10-6 м2. Одна из боковых граней используемых стальных образцов ТАН образует

угол с вертикалью ф = 0, а оппозитная ей грань наклонена к вертикали под углом, варьируемом в диапазоне ф = 0 ± 45°. Для повышения надежности измерений путем устранения шумового фона, вызванного преимущественно краевыми объемными модами, на боковых поверхностях призмы выполнены канавки, играющие роль рас-сеивателей и покрытые поглотителем. Рабочая частота пьезоэлектрического преобразователя (ПЭП) рэлеевских волн 1,8 МГц. При проведении каждого измерения контактная среда, находящаяся вне области контакта ТАН с подложкой, тщательно удалялась.

Источником и приемником недетектиро-ванных импульсных сигналов служат электронные блоки генераторов УД2-12, осциллограф С1-71, измеритель временных интервалов И1-8, используемый для синхронизации и анализа параметров акустического сигнала, а также генератор стандартных сигналов, подаваемых на вторую развертку осциллографа для более точного определения амплитуды сигнала методом сравнения амплитуд.

На первом этапе исследований в режиме эхо измерялась отражающая способность моделируемой области контакта ТАН с подложкой преимущественно в зависимости от угла ф и толщины контактной смазки, а расчет искомых коэффициентов Е, = {Ж, Ож*, О^*, Я*} проведен путем использования формул (2)-(5). Как следует из них, при определении нормализованных коэффициентов представляется возможным исключить влияние корректирующих коэффициенты /Т и/Е на точность измерений.

Высота щелевого зазора регулировалась с помощью специальных подкладок, толщина которых измерялась с точностью 1 мкм. Для определения абсолютного значения коэффициента отражения волны Стоунли ЯS был использован зеркально-теневой метод, позволяющий вначале найти отношение амплитуд ПАВ А12 = А1Т/А2Т, где А соответствует прямой волне, прошедшей через зону контакта ТАН с подложкой на приемный ПЭП, а А2Т - двукратно отраженной от границы 2 и 1. Учтя, как и прежде, ослабление ВС в результате прохождения (при переотражении) ею дополнительного пути в зоне контакта путем ввода некоторого коэффициента у получим, что Я3 = (А12у)0,5. Коэффициент у определялся путем анализа изменения величины отношения амплитуд А в результате варьирования поперечного размера призмы и акустической базы про-

звучивания. Оцененная погрешность измерений не превышала 10-15 %.

Результаты экспериментальных исследований и их обсуждение

Результаты экспериментальных исследований прохождения и трансформации ПАВ на образцах с моделируемой скользящей границей АНГ приведены на рисунках 2-4, где представлены нормализованные зависимости Е,(ф) = {Я^*, Ож*, О^*}, а также Ж. На рисунках 5-7 приведены некоторые предложенные схемы и конструкции устройств для акустических измерений, а также данные, свидетельствующие о возможности использования разработки для измерения глубины упрочненного термообработкой поверхностного слоя.

- Ж.ЗБ 2

/ у 20 -

10 - II.......

-45 -30 -15 0 15 30 45

ф.гроЗ

Рисунок 2 - Зависимость коэффициента отражательной способности тела нагрузки (призмы) R от угла наклона ее боковой грани: 1 - положение призмы I; 2 - положение призмы II (согласно рисунку 1)

Figure 2 - The reflectivity coefficient of the load body (as a prism) ЛЕ depending on the angle of inclination of its side face: 1 - the position of the prism I; 2 - the position II (according to the explanation in Figure1)

Как видно, все представленные впервые зависимости, характеризующие свойства моделируемого объекта как ТАН или отражателя {R, R*}, так и коэффициентов прохождения ПАВ через границы 1 и 2, сопровождающегося трансформацией мод N*={Dsr*, Drs*} от ф, имеют максимум преимущественно в окрестности ф = 0. Как видно, левая и правая ветви этих угловых зависимостей относительно оси ф = 0 асимметричны. Причем для коэффициента, характеризующего отражательную способность тела призмы (формула (2)Х величина R ~ (Rr)2 (Drs)1(dsr)1(rr)1-1

и достигает 32-34 дБ.

Наименьшее же значение Ж ~ 18-20 дБ, имеет место в том случае, когда отражателем служит граница 2, а ф > 0. Как видно из хода функции Ж(ф), превалирующий вклад в ее изменение вносит именно коэффициент отражения ВР, который при ф = 0 равен (К)2 = (К^ = Я8.

i 0 . R*s,N*.aB 2

-4

-8

-12 ......... ......

-45

-30

-15

15

30 45

ф. гроВ.

Рисунок 3 - Влияние угла наклона боковой грани призмы на нормализованный коэффициент отражения волны Стоунли RS* (1) и коэффициент прохождения границ контакта тела акустической нагрузки с подложкой N (усредненный (2)) приустановке призмы в положение I и II, что поясняется рисунком 1

Figure 3 - Influence of the angle of inclination of the side faces of the prism on the normalized reflection coefficient ofthe Stonely wave R* (1) and the transmission coefficient of the contact boundaries of the body of the acoustic load to the basic body N* (averaged (2)) when the prism in position I or II according to the explanation in Figure 1

,^3

o. 31 0, 32 0, 33 0, 14 0, 35 0,

/

/

2

1 —

Рисунок 4 - Влияние безразмерной толщины контактного слоя жидкости h*= h/A,R на амплитуду отраженной поверхностной волны от передней (1), задней (2) границы призмы, а также от ребра подложки (3), выполненной в виде параллелепипеда

Figure 4 - The influence of the dimensionless thickness of the contact layer of the liquid on the amplitude of the reflected surface wave from the front (1), back (2) prism boundary, as well as from the edge of the basic solid (3), made in the form of a parallelepiped

/

О

D

X,

^—*-&-Э

Рисунок 5 - Схемы использования устройства с твердотельным отражателем поверхностных волн для ультразвукового контроля и измерений по данным амплитудных и амплитудно-угловых параметров отраженного сигнала (а) (вид сверху) и по данным скорости поверхностной волны (вид сбоку), определяемой импульсно-фазовым методом с характерной осциллограммой (b): 1 - излучающий преобразователь; 2 - объект контроля; 3 и 4 - грани выступа тела отражателя, формирующие отражающую поверхностные волны границу с вершиной 5; 6 - отражающая волны граница 2, формирующая опорный сигнал; 7 - приемный преобразователь; c: I - зондирующий, II и III - принимаемые преобразователем 7 импульсы

Figure 5 - Schemes of use of the device with a solidstate reflector of surface waves for ultrasonic evaluation and measurements by using amplitude parameters of the reflected signal vs. probe angle (a) (top view) and the velocity of the surface wave (side view) determined by pulse-phase method with characteristic oscillogram (b): 1 - emitting probes; 2 - test object to test; 3 and 4 - lateral sides of the reflector forming a reflecting surface waves boundary with the top 5; 6 - reflecting wave boundary 2, forming the reference signal; 7 - receiving probes; c: I - exciting, II and III - impulses.receiving by probe 7

6

1

a

b

c

Рисунок 6 - Пример схемы использования ультразвукового устройства с твердотельным отражателем поверхностных волн для контроля объекта со сложным рельефом поверхности: 1 - преобразователь; 2 - объект контроля; 3 - поверхностная несплошность; 4 - отражатель поверхностных волн

Figure 6 - Scheme of using an ultrasonic device with a solidstate reflection of surface waves to test an object with a complex surface relief: 1 - probe; 2 - object to be tested; 3 - surface discon-tinuity; 4 - reflector of surface waves

абсолютную величину коэффициента отражения волны Стоунли RS при изменении угла наклона как одной, так и двух граней призмы ТАН в исследованном угловом диапазоне. Это же касается определения абсолютной зависимости коэффициента прохождения ПАВ N(9), определяемой теневым методом.

Что касается зависимостей нормализованных коэффициентов прохождения ПАВ через область с акустической нагрузкой, то N = DRS*(ф) и N* = ^ж*(ф) построены, как указывалось выше, путем разворота призмы ТАН на 180° относительно акустической оси. Установлено, что различие между ними по величине во всем исследованном диапазоне углов наклона одной из граней призмы не превышает половины погрешности измерений. Это позволило отобразить указанные зависимости одной «усредненной кривой» -N*^). Т.е. для углов 45° > ф > 0 зависимость N*^) монотонно убывает не более, чем на 10 %, а для 0 > ф > -45° - на ~ 70-80 %. Таким образом, если угол наклона боковой поверхности призмы 0 > ф > 45°, то можно положить, что:

Drl = (^ )ф=0 = const.

D.

SR

D

SR

Рисунок 7 - Зависимости нормализованной амплитуды поверхностных волн, отраженных отражателем, от угла падения волны в на частоте 1,8 МГц: безразмерная глубина упрочненного слоя h* = h,L = 0 (1); 0,4 (2); 1,5 (3)

Figure 7 - Dependences of the normalized amplitude of reflected surface waves vs. the probe angle в when the depth of work-hardened by the currents of high-purity of steel layer is of h* = h,k = 0 (1); 0,4 (2); 1,5 (3) and when signal frequency is of 1,8 MHz

В результате исследований согласно приведенной выше методике было установлено, что определенная с погрешностью, не превосходящей 10-15 %, величина RS ~ 3,2. Это позволяет по данным, приведенным на рисунке 2, оценить

Что касается абсолютного значения функции прохождения N то, как показывают данные эксперимента, в пределах указанной погрешности измерений (10-15 %) наличие ТАН в виде прямоугольной призмы сопровождается уменьшением N приблизительно в 3,3 раза.

Причина хода указанных выше угловых зависимостей обусловлена сложным характером трансформации ВР и ВС в окрестности границ ТАН, что поясняется уравнением баланса (6). Как следует из экспериментальных данных, представленных на рисунке 2, энергией ВР, отраженных от передней границы ТАН, можно пренебречь. (Однако при проведении прецезион-ных ультразвуковых измерений этот фактор необходимо учитывать.) Как установлено, часть энергии ВР теряется при переходе этих волн на боковую поверхность призмы (граница 1), а часть - трансформируется в объемные моды. Подобный процесс имеет место и при трансформации волны Стоунли в волну Рэлея на границе 2.

Используя подобную указанным методическим особенностям процедуру, были получены данные сравнения коэффициента отражательной способности ТАН и коэффициента отраже-

4

ния ВР от ребра параллелепипеда (ЯЯ)П, которое идентично отражению от бесконечной трещины. Установлено, что величина (ЯЯ)П больше Ж всего приблизительно на 5 дБ. Если же ПАВ проходит через зону контакта ТАН, отражается от ребра параллелепипеда и возвращается обратно, то эта разница достигает приблизительно 18-19 дБ. Увеличение безразмерного значения толщины контактного слоя ^ = Мк (рисунок 4) приводит к росту этой разницы до 47-48 дБ, достигаемой в окрестности минимума исследуемой зависимости при ^ ~ 0,02. В то же время на величину Ж слабо сказывается варьирование ^ в диапазоне от нуля до 0,05. Изменение же скорости ПАВ в контактном зазоре, измеренной импульсно-фа-зовым методом, может составить ~ 4-5 %.

Исходя из анализа акустического тракта и полученных опытных данных о коэффициентах прохождения и отражения ПАВ, можно сделать вывод о возможности проведения на практике контроля дефектности поверхности объектов как непосредственно расположенных в зоне контакта материалов, так и с противоположной стороны технологического объекта.

Необходимо отметить, что подобные указанным выше процессы трансформации волн имеют место и в том случае, когда объектом исследования служит твердое тело с выступом, реализуя модель жесткого сцепления поверхности подложки с ТАН. В этом случае, как показано в работе [4], превалирующая часть потока акустической энергии трансформируется на выступе в поперечную вертикально поляризованную Т-моду. (Более подробный анализ механизма трансформации упругих мод представляется дать во второй части статьи.) Если же выступ имеет форму прямоугольной призмы, то коэффициент прохождения акустического сигнала, включающий трансформацию мод согласно схеме ПАВ^Т-мода^ПАВ, может быть представлен в виде:

N = Н2 = (ОД (ОД/,

где ОЕТ и От - коэффициенты трансформации ВР в Т-моду и обратно; /Т - множитель, обуславливающий ослабление Т-моды или краевой поперечной волны, распространяющейся в пространстве х1 < х < х При работе в режиме эхо аналогичный коэффициент:

ж = (оДОДЯ^/Е [(ЯДГ1.

Сравнительные экспериментальные данные показали, что коэффициент прохождения ПАВ по подложке с ТАН и скользящей контактной границей не более чем на 10 % превышает аналогичную величину при реализации жесткого контакта нагрузки с подложкой. Что касается отражательной способности зоны контакта ТАН с подложкой при жесткой границе между ними, то ее величина составляет всего R ~ 7-8 дБ. При реализации же скользящей границы (рисунок 2) отражательная способность ТАН более чем на порядок (25-26 дБ) больше, а амплитуда принимаемого акустического сигнала, отраженного от 2-й (отражающей) границы, превышает всего на 7-8 дБ.

Таким образом, при соединении одинаковых материалов методом склейки, пайки и др. характер изменения функции R(t) за характерное время t = t может служить высокочувствительным показателем качества технологического процесса соединении одинаковых материалов. При этом отношение величин отражательных способностей ТО

* Wt _0

ТО =- для указанных случаев превышает

TOt=t порядок.

Исходя из анализа механизма отражения и прохождения ПАВ, а также предварительных данных эксперимента, полученных для случая контакта ТАН из дюрали со стальной подложкой (и наоборот), следует предположить эффективность использования временной зависимости функции R(t) для контроля аналогичного технологического процесса и при другом сочетании соединяемых материалов ТАН и подложки.

Создание опорного сигналаприменительно к ультразвуковым измерениям и контролю

При проведении ультразвукового контроля физико-механических свойств твердых тел по данным скорости ПАВ, ее амплитудно-частотным характеристикам, а также при выявлении поверхностных и подповерхностных дефектов требуется высокая точность измерений, существенно зависящая от шумового фона, вызванного указанными выше причинами. Для устранения этого фактора предлагаются некоторые схемы измерений, поясняемые на рисунках 5 и 6. В частности, предлагается использовать ТАН в качестве высокоэффективного отражателя ПАВ при определении, например, структуры чугунов, глубины упрочненного слоя, по данным скорости

ПАВ, коррелирующей с глубиной упрочненного слоя, где СПАВ = СПдВ(й). Последняя может быть определена импульсно-фазовым методом по формуле: СПдВ = L*/At, где L* - акустическая база измерений, At - время распространения ПАВ (рисунок 5), либо по данным угла ß* максимума амплитуды отраженной от отражателя поверхностной волны, где СПАВ = Cj/sinß*.

Дополнительное снижение на десятки децибел шумового фона, попадающего в режиме приема на ПЭП, достигается путем выполнения передней отражающей ПАВ границы ТАН оптимальной формы. Например, эта граница может иметь вид треугольного уступа или выступа, обращенного к излучающему ПЭП (рисунок 6), угол которого 2ф образован двумя линейными участками отражающей ПАВ границы с вершиной в точку S. В этом случае превалирующая часть потока отраженной от границы энергии ПАВ выводится из зоны расположения ПЭП. Используя представление лучевой акустики, представляется возможным в первом приближении оценить соотношение между поперечным размерами призмы нагрузки 2b, преобразователя 2d и его расстоянием s до вершины S. Так, если эта отражающая волну Рэлея граница имеет форму уступа, то указанное соотношение примет вид:

(ä - bctgф)/ d > 2tg2ф, а если граница имеет форму выступа, то: -s / d > tg2ф,

где s* = s - ЬС§ф - расстояние между передней гранью ПЭП и крайними (ближайшими) точками отражающей границы в виде уступа.

Как показывают численные оценки, при одинаковых параметрах {s*/d, b/d) наиболее компактной является конструкция отражателя с выступом, обращенным к излучающему волны ПЭП (рисунок 6). Как видно из приведенных формул, именно в этом случае влияние поперечного размера отражателя на его отражающую ПАВ способность существенно нивелируется. При этом представляется возможным в 1,5-2 раза и более сократить размеры конструкции ультразвукового устройства при одновременном уменьшении на десятки децибел соотношения шум-сигнал. Несомненно, это позволит значительно расширить номенклатуру контролируемых на качество их структуры объектов.

Отметим, что отражающая ПАВ граница 2 может иметь различное направление по отноше-

нию к акустической оси системы, а также разную форму. Так что ТАН может быть использована в качестве «механического девиатора» поверхностных волн, отклоняющего поток энергии ПАВ под некоторым углом к акустической оси. Во втором же случае - осуществлять фокусировку ПАВ перед границей акустической нагрузки 1. При этом отраженные волны от границы 2, выполненной в виде окружности с определенным радиусом r0, концентрируются на расстоянии от этой границы в окрестности своеобразного фокуса F ~ r0/2. Использование таких схем управления направленностью ПАВ представляет интерес для контроля объектов в труднодоступных местах, имеющих изменяемую структуру поверхности и требующих повышенной стабильности и чувствительности измерений.

На рисунке 7 демонстрируется эффективность применения предложенного устройства для оценки глубины упрочненного токами высокой частоты слоя металла по данным корреляционной зависимости скорости ПАВ:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

СПАВ = СПДВ^) ~ (£УPe)0•5, где Ee и Pe рассматриваются как эффективные значения модуля Юнга

и плотности контролируемой среды, локализованной на глубину h ~ Апав. Как видно, измерения могут производиться как импульсно-фазовым, так и амплитудно-угловым методом по формуле СПАВ = С (sin в*)-1, где в* - угол падения волны на металл, обеспечивающий максимум амплитуды возбуждаемой ПАВ. Отметим, что в данном случае нет необходимости усложнять измерительную процедуру, используя два ПЭП, работающих в теневом режиме прозвучивания.

Заключение

На основе анализа акустического тракта получены соотношения, связывающие амплитудные параметры поверхностных волн, отраженных и проходящих через зону акустического контакта (тонкий смазочный слой) создаваемого стальной подложкой с телом акустической нагрузки в виде призмы и коэффициентами прохождения и отражения волн Стоунли через ее границы, а также отражательную способность тела нагрузки. Предложена методика и впервые экспериментально определены зависимости указанных коэффициентов при изменении угла наклона одной из боковых граней призмы в диапазоне ф = 0 ± 45°, безразмерной толщины контактной прослойки (0-0,05) и ее ориентации

вдоль акустической оси измерительной системы и оппозитно ей.

Впервые экспериментально установлено, что представленные зависимости, коэффициенты отражения и прохождения волн Рэлея и Стоун-ли через переднюю (1-я граница) и заднюю (2-я) границы скользящего контакта тела призмы с подложкой (при одновременной трансформации мод), а также коэффициент отражательной способности, определяемый отношением амплитуд принимаемых отраженных акустических сигналов от 2-й и 1-й границы, являются асимметричными функциями от угла наклона боковой поверхности призмы и имеют преимущественно максимум в окрестности ф^-0. Причем это отношение, измеренное в режиме эхо, принимает максимальное значение (~ 32-34 дБ) в случае использования в качестве отражающего тела прямоугольной призмы, что более чем на порядок превышает аналогичный коэффициент, измеренный при наличии жесткого контакта тел, в то время как коэффициенты прохождения акустического сигнала сравнимы по величине.

На основании этих данных предложено использовать коэффициент отражательной способности в качестве высокоинформативного параметра, характеризующего качества соединения материалов пайкой, сваркой, склеиванием и др.

Определенные с погрешностью не более 10-15 % зеркально теневым и теневым способом максимальные значения коэффициенты отражения волны Стоунли и прохождения поверхностной волны через акустическую нагрузку составляют приблизительно 0,31 и 0,33 соответственно. Установлено, что угловые зависимости коэффициента прохождения поверхностной волны в пределах погрешности измерений совпадают при повороте призмы на 1800. Причем для ф > 0 этот коэффициент уменьшается не более чем на 1 дБ, а при ф < 0 - приблизительно на 4-5 дБ. Для коэффициента же отражения волны Стоунли наблюдается обратная зависимость: при ф > 0 его величина падает приблизительно на 12 дБ, а для ф < 0 его изменение не превосходит 1 дБ.

Предложены способы ультразвукового контроля и даны рекомендации по использованию результатов исследования:

- для определения физико-механических свойств поверхности материалов, включая определение глубины упрочненного слоя, по данным скорости поверхностной волны, измеренной

импульсно-фазовым и амплитудно-угловым методом, путем использования отражателей поверхностных волн, что позволяет сократить акустическую базу прозвучивания объекта в 1,5-2 раза и существенно расширить номенклатуру контролируемых объектов;

- для выявления дефектов, расположенных в труднодоступных местах и на объектах со сложным рельефом поверхности и неоднородной структурой, используя возможности фокусировки ПАВ путем выполнения границы 2 с заданным радиусом кривизны либо дефлекции волны механическим способом.

Полученные экспериментально данные могут быть использованы для оценки возможности диагностирования состояния трущихся поверхностей, для определения наличия на них поверхностных дефектов, а также контроля качества сцепления соединяемых материалов. Определенные коэффициенты отражения и прохождения ПАВ представляют интерес при теоретическом моделировании эффектов трансформации поверхностных волн их распространения.

Работа выполнена при частичной поддержке Белорусского Республиканского фонда фундаментальных исследований в рамках договора № Т17-159 от 18.04.2018.

This work was carried out with the partial support of the Belarusian Republican Fund of fundamental investigations within the bounds of contract no. Т17-159 of 18.04.2018.

Список использованных источников

1. Kruger, S.E. Detection of kissing bond in extruded aluminum by laser-ultrasound / A.J. Bakker, S.E. Kruger // Review of Progress in Quantitative Nondestructive Evaluation. - 2007. - Vol. 27A. - P. 279-285.

2. Pecorary, C. Nonlinear interaction of plane ultrasonic waves with an interface between rough surfaces in contact / C. Pecorary // JASA. - 2003. - Vol. 113. -P. 3065-3072.

3. Kim, H.J. Evaluation of thin coating layers using Rayleigh-like waves / S.-J. Song Golden, D.-Y. Kim, S.-D. Kwon, H.-J. Kim // Review of Progress in Quantitative Nondestructive Evaluation. - 2007. -Vol. 27B. - P. 1066-1073.

4. Баев, А.Р. Распространение волны Рэлея в твердых телах с технологическим выступом / А.Р. Баев, О.С. Сергеева,М.В. Асадчая, Г.Е. Коновалов //

Приборы и методы измерений. - 2011. - № 2. -С. 121-129.

5. Jerzak, W. Modeling Rayleigh and Stonely Waves and other Interface and Boundary Effects with the Parabolic Equation / W. Jerzak, W.L. Siegman., H.D. Collinz // JASA. - 2005. - Vol. 117. -P. 3497-3503.

6. Aladeddin, E. Combining NMR and Stonelyn Analysis for Better Estimation and Permeability in Carfonate Reservoirs / E. Aladeddin, H. Tchambaz and N. Al-Adani // Paper E. Proceedings, Formation Evaluation Symposium of Japan Society Petrophysist and Well Log Analysitst, Japan. - 2005. - Chapter, 7 p.

7. Ультразвуковое устройство с отражателем поверхностных волн / А.Р. Баев, О.С. Сергеева, М.В. Асадчая, М.А. Тищенко: пат. 9911 Респ. Беларусь, МПК (1023.09) G 01 N 29/04 / u 20130571; опубл. 28.02.14 // Афщыйны бюл. / Нац. цэнтр штэлектуал. уласнасщ. - 2014. - № 4. - С. 36.

8. Аббакумов, К.Е. Влияние нарушения акустического контакта на распространение волн Стоунли вблизи границы твердых полупространств / К.Е. Абакумов, Р.С. Коновалов // Дефектоскопия. - 2008. - № 3. - С. 52-58.

9. Abbakumov, K.E. Britvin Informative opportunities of wave processes at the boundary of solid mediums in conditions of broken acoustic contact / K.E. Abbakumov, V.A. Konovalov // Abstracts of 10th European conference on Non-destructive testing. -Moscow, 2010. - P. 173-175.

10. Люткевич, А.М. Акустические поля мало-апертурных. Поперечные волны, излучаемые прямоугольным источником нормальной силы / А.М. Люткевич, А.В. Жуков, А.А. Самокрутов, В.Г. Шевалдыкин // Контроль. Диагностика. - 2004. - № 4. - C. 3-8.

References

1. Kruger S.E., Bakker A.J. Detection of kissing bond in extruded aluminum by laser-ultrasound. Review of Progress in Quantitative Nondestructive Evaluation, 2007, vol. 27A, pp. 279-285.

2. Pecorary C. Nonlinear interaction of plane ultrasonic waves with an interface between rough surfaces in contact. JASA, 2003, vol. 113, pp. 3065-3072.

3. Kim H.J., Song Golden S.J., Kim D.Y., Kwon S.D. Evaluation of thin coating layers using Rayleigh-like waves. Review of Progress in Quantitative Nondestructive Evaluation, 2007, vol. 27B, pp. 1066-1073.

4. Baev A.R., Sergeeva O.S., Asadchaya M.V, Konovalov G.E. [Propagation of Rayleigh waves in solids with protrusion]. Devices and methods of measurements, 2011, vol. 2, no. 2, pp. 121-129 (in Russian).

5. Jerzak W., Siegman. W.L., Collinz H.D. Modeling Rayleigh and Stonely Waves and other Interface and Boundary Effects with the Parabolic Equation. JASA, 2005, vol. 117, pp. 3497-3503.

6. Aladeddin E., Tchambaz H., Al-Adani N. Combining NMR and Stonelyn Analysis for Better Estimation and Permeability in Carfonate Reservoirs. Paper E. Proceedings, Formation Evaluation Symposium of Japan Society Petrophysist and Well Log Analysitst, Japan, 2005, Chapter, 7 p.

7. Baev A.R., Sergeeva O.S., Asadchaya M.V., Tishchenko M.A.: Ultraszvukovoye ustroistvo s otraza-telem poverchnostnykh voln [Ultrasonic arrangement with reflector of surface waves]. Patent RB no. u20130571, 2014.

8. Abbakumov K.E., Konovalov V.A. [Influence of the broken acoustical contact on propagation of Stonely waves at the boundary of solid mediums]. Defectoscopiya [Defectoscopy], 2008, no. 3, pp. 52-58 (in Russian).

9. Abbakumov K.E., Konovalov V.A. Britvin Informative opportunities of wave processes at the boundary of solid mediums in conditions of broken acoustic contact. Abstracts of 10th European conference on Nondestructive testing, Moscow, 2010, pp. 173-175.

10. Lutkevich A.M., Zukove A.V., Samokrutov A.A., Shevaldykin V.G. [Acoustical fields of the probes with small aperture. Transverse waves, radiated by rectangular source of the normal force]. Kontrol. Diagnostica [Control. Diagnostics], 2004, no. 4, pp. 3-8 (in Russian).

i Надоели баннеры? Вы всегда можете отключить рекламу.