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

Автор: Шайдуров Владимир Викторович, Непомнящий Олег Владимирович, Вейсов Евгений Алексеевич

Журнал: Сибирский аэрокосмический журнал @vestnik-sibsau

Рубрика: Математика, механика, информатика

Статья в выпуске: 1-2 (22), 2009 года.

Бесплатный доступ

Рассмотрены методы расчета ультразвуковых преобразователей. Предложена математическая модель ультразвукового преобразователя для сканирования высокотемпературных изделий. Даны практические рекомендации по реализации микропроцессорной системы ультразвукового сканирования этих изделий.

Ультразвук, модель, микропроцессор, датчик

Короткий адрес: https://sciup.org/148175844

IDR: 148175844

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

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

Особое внимание при выборе преобразователя следует уделить условиям контроля, в первую очередь свойствам материала и условиям контакта преобразователя с объектом контроля [1].

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

Особую сложность вызывает ультразвуковое сканирование высокотемпературных изделий, поскольку с повышением температуры в исследуемом объекте происходит искажение распространения звуковых волн. Разрушение преобразователя имеет место и при контакте с поверхностью, нагретой более 50 ° С [2].

В качестве примера построения модели излучателя рассмотрим ультразвуковое исследование нагретого до 300 ° С стального образца (сталь 4340), имеющего толщину до 200 мм.

В общем случае под ультразвуком понимают излучение на частоте выше 20 кГц, однако при ультразвуковом исследовании используют диапазон от 100 кГц до 50 мГц. Частоту излучения определяют как

f = р CD где f – частота; Т – период незатухающей волны.

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

X = с, X = cT, (2)

где X - длина волны; с - скорость звука в материале; f – частота; Т – период.

Для определения частоты излучателя воспользуемся табличными данными (табл. 1). Однако изначально сле- дует определить характер излучаемой волны (продольная или поперечная). При ультразвуковом исследовании используются и продольные, и поперечные волны, а также в специальных случаях – поверхностные (релеевские) и плоские волны (волны Лэмба).

Чем меньше длина волны, тем больше на нее влияет среда распространения, и поэтому при анализе материалов используется максимально допустимая для данного материала частота.

Известно [3], что при применении для ультразвукового исследования поперечной волны чувствительность к мелким дефектам повышается, а затухание с расстоянием увеличивается, следовательно в случае измерения толщины требуется применение излучателя продольной волны, который менее чувствителен к мелким дефектам и имеет меньший коэффициент затухания.

Таким образом, используя соотношения (1) и (2) и рассматривая единичную волну (приняв X = 1), для исследования стали 4340 получаем частоту излучателя продольной волны порядка 5…6 мГц.

При проведении ультразвукового исследования измеряются две величины: время пролета (время прохождения ультразвуковой волны через материал) и амплитуда принятого сигнала. Эти параметры определены следующим соотношением:

Т = f •

где Т – толщина материала; с – скорость звука в материале; t – время пролета. Используя это соотношение, рассчитаем время пролета для стального бруска толщиной 200 мм. Оно составит 0,001 с.

Для измерения плотности материала по величине затухания или для выявления дефектов используют данные об относительном изменении амплитуды сигнала, которое выражают в децибелах и определяют логарифмичес- ки исходя из соотношения двух сигналов:

dB = 20log10

f A ) —L

l A 2 J

где dB – децибелы; А 1 – амплитуда первого сигнала;

А 2 – амплитуда второго сигнала.

Примеры коэффициента ослабления сигнала для раз-A личных соотношений 1 приведены в табл. 2.

A 2

Основываясь на вышеизложенном, с целью уточнения результатов сделаем дополнительные расчеты частоты используемого излучателя.

При проведении ультразвукового сканирования получена следующая форма волны (рис. 1) [4].

Форма волны датчика и спектральный анализ выполняются в соответствии с определениями и условиями проведения теста ASTM E1065. Обычно для частотного анализа в качестве единиц измерения используются мегагерцы, для анализа формы волны – микросекунды и для определения пика амплитуды – децибелы [1].

Рассмотрим пик частоты прохождения ультразвука в исследуемом материале. С увеличением частоты амплитуда выходного сигнала растет. Здесь следует обратить внимание на верхний и нижний пределы частоты. При достижении нижнего предела коэффициент ослабления становится равным 2, т. е. мы имеем 6-децибеловый сигнал (рис. 2). При достижении верхнего предела коэффициент начинает монотонно падать. Таким образом, нас интересует только этот 6-децибеловый диапазон.

Известно соотношение между мегагерцевой полосой пропускания и длительностью прохождения волны (рис. 3) [1]. В этом соотношении разброс на –40дБ больше, так как 1%-й задний конец формы волны содержит очень небольшое количество звуковой энергии и поэтому оказывает очень малое влияние на анализ полосы пропускания. Из-за рассеивания это наиболее подходящий вариант для точного определения формы волны за промежуток времени, мкс, и спектров в частотной области.

Таблица 1

Акустические свойства материалов [1]

Материал

Скорость продольной волны, м/сек

Скорость поперечной волны, м/сек

Акустический импеданс, кг/(м2 с еск)∙1088)

Алюминий

6 320

3 130

17,06

Бериллий

12 900

8 880

23,5

Латунь морская

4 430

2 120

37,30

Кадмий

2 780

1 500

24.02

Ниобий

4 920

2 100

42,16

Медь

4 660

2 260

41,61

Золото

3 240

1 200

62,60

Инконел (жаропрочный и жаростойкий сплав на никелевой основе)

5 820

3 020

49,47

Железо

5 900

3 230

45.43

Чугун (с малой скоростью)

3 500

2 200

25,00

Чугун (с высокой скоростью)

5 600

3 220

40,00

Свинец

2 160

700

24,49

Молибден

6 250

3 350

63,75

Никель чистый

5 630

2 960

49.99

Платина

3 960

1 670

84,74

Серебро

3 600

1 590

37,76

Сталь 1020

5 890

3 240

45,63

Сталь 4340

5 850

3 240

45,63

Сталь   марки   302   аустенитная

нержавеющая

5 660

3 120

45,45

Сталь 347 аустенитная

5 740

3 090

45,40

Жесть

3 320

1 670

24,20

Титан марки Ti 150A

6 100

3 120

27,69

Вольфрам

5 180

2 870

99,72

Цинк

4 170

2 410

29,61

Таблица 2

Коэффициенты ослабления сигнала

А A 1,1,%%//А А 2, %

Коэффициент

dB

100

70,71

1,414 2

3

100

50

2

6

100

25

4

12

100

10

10

20

100

1

100

40

Основываясь на данных рис. 3, можно определить частоту датчика.

Рис. 1. Волна на уровне 6 дБ (50%-й пик амплитуды)

Рис. 2. Пик частоты прохождения для исследуемого материала

Рис. 3. Длина волны и полоса пропускания

Если, например, требуется – 14-децибеловая волна за 1 мкс, то какой следует выбрать датчик? На графике (см. рис. 3) полоса пропускания от 1 до 1,2 МГц соотносится примерно с 1 мкс – 14-децибеловой протяженности длины волны. Принимая номинальный 50%-й фракционный широкополосный датчик, получаем среднюю номинальную частоту от 2 до 2,4 МГц. Поэтому может быть использован датчик частотой 2,25 МГц или 3,5 МГц.

Не менее важны исследования акустического импеданса материала, который является показателем противодействия смещению его частиц звуком:

Z = P c,

где Z - акустический импеданс; P - плотность материа- ла; с – скорость звука в материале.

Граница между двумя материалами с разным акустическим импедансом называется акустической границей . При нормальном падении звука в акустическую границу часть звуковой энергии отражается, а другая часть продолжает свое перемещение в материале. Потеря энергии, дБ, при прохождении сигнала из среды 1 в среду 2 рассчитывается следующим образом:

( dBloss = 1010g1ol 7-1-2-2 I ,

V ( Z 1 + Z 2 ) )

где Z 1 – акустический импеданс первого материала; Z 2 – акустический импеданс второго материала.

Необходимо учесть потерю энергии децибелового сигнала в среде 1 при отражении от акустической границы со средой 2:

(            2)

dB loss = 10log io | (_1    _2) I .                (7)

v ( Z 1 + Z 2 ) )

Например, потеря сигнала при его распространении из воды ( Z = 1,48) в сталь 1020 ( Z = 45,41) составит –9,13 дБ. Потеря децибелового отраженного эхо-сигнала в стали 1020 составит –0,57 дБ, что также соответствует потере эхо-сигнала при распространении в обратном направлении. При Z 2 Z 1 происходит переворот формы волны.

В материале происходит увеличение затухания ультразвука. При отсутствии значительных отражателей вы- деляют три варианта затухания:

– преломление;

– рассеивание;

– поглощение.

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

Звуковое поле датчика (рис. 4) описывается следующими зависимостями:

N = Df,                (8)

4 C

N = D 2 ,                   (9)

4X где N – протяженность ближнего поля; D – диаметр элемента; f – частота; с – скорость звука в материале; X - длина волны.

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

1, 02 Fc

= fD ’

BD = 0,256 8 DS F ,              (11)

где BD – диаметр луча; F – фокусное расстояние; с – скорость звука в материале; f – частота, D – диаметр элемента; SF – нормированное фокусное расстояние.

Для определения фокусной зоны датчика (рис. 5) рассмотрим начальную и конечную точки, которые расположены в точке падения амплитуды сигнала «импульс – эхо» до –6-децибеловой амплитуды в точке фокуса (рис. 6). Длина фокусной зоны вычисляется по формуле

F = NS^ f ( 2 + 0,5 S f ) , (12)

где F – фокусная зона; N – ближнее поле; SF – нормированное фокусное расстояние.

Рис. 4. Звуковое поле датчика

Рис. 5. Фокусная зона датчика: ZВ – начало фокусной зоны;

FZ – фокусная зона; ZE – конец фокусной зоны; D – диаметр элемента

Рис. 6. Конечная и начальная точки фокусной зоны –6 дБ при измерении амплитуды сигнала

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

Угол распространения луча плоского датчика –6-де-цибелового импульса рассчитывается следующим образом:

sin |a| = 0,514 —,              (13)

12 J fD где a / 2 - половина угла рассеивания луча между точка ми –6 дБ.

Рис. 7. Рассеивание луча и половинный угол

Таким образом, угол рассеивания луча можно варьировать, изменяя частоту и диаметр датчика, на чем основана работа датчиков ультразвукового контроля.

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

Таким образом, математическая модель преобразователя ультразвукового сигнала адекватно описывается выражениями (1)…(12).

При ультразвуковом сканировании высокотемпературных изделий, в частности бруска стали, нагретого до 200.. .600 ° С, согласно приведенным выше расчетам, целесообразно применять ультразвуковой датчик с рабочей частотой порядка 5 MГц и диаметром 25 мм. Кроме того, необходима дальнейшая разработка методики проведения исследований и реализация экспериментального образца микропроцессорной системы на основе высокопроизводительного сигнального процессора DSP.

Статья научная