Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Ультразвуковой контроль и регулирование технологических процессов

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
10.93 Mб
Скачать

нию уровня излучения и отношения информативный сигнал/акустические помехи.

Полученные расчетные формулы применимы и для поперечных t колебаний при вводе t в индексы обозначений параметров. Они получены для непрерыв­ ных колебаний и пакетных длительностью, превышающей 12 полуциклов ко­ лебаний, в широком частотном диапазоне, включая круговые частоты: соо на ПВР и со. на ЭМР (последний рассмотрен в гл. 6). Для часто применяемого возбуждения ПП посредством ЭИ в виде двух полуциклов синусоидальных колебаний или прямоугольного импульса полуцикловой длительности давле­ ние излучения уменьшается почти вдвое в сравнении с непрерывным излуче­ нием. При повышенных V возбуждения нарушается линейность зависимости p/V из-за роста противо-ЭДС пьезореакции ПП.

1.2. Общие параметры поля излучателя

Принимаемый информативный сигнал зависит также от структуры УЗ-поля в средах и описывается системой, состоящей из известного уравнения Г. Гель­ мгольца для давления р и уравнения волновой функции \|/:

\|f =[exp(-ikr)\/ г ,

(1.13)

где г — расстояние от точки излучения А (рис. 1.3) круглого (радиусом а) из­ лучающего ПП до точки В в поле, где определяют искомое давление р; к = 2ти/ X— волновое число; X — длина УЗ-волны в среде с импедансом Zh.

Упрощенное решение системы в зависимости от координаты z для rm= (z + + а2)0,5 дается формулой Г. Штенцеля [41,42, 62]:

p(z) = р0[ехр( - ikz) - ехр( - ikrm)]exp(ia>0.

(1.14)

Рис. 1.3. Поле акустического излучателя

Для коротковолновых (ка » 1) гармонических колебаний и длительных им­ пульсов полное решение системы волновых уравнений можно выполнить на основе теоремы Г. Грина:

ДОКуЭри^ - (pdy/dn)Si]<й,= 0,

(1.14 а)

где \|/ — потенциал точечных источников; п — внутренняя нормаль к элементу dSxзамкнутой поверхности интегрирования Sv В точке поля при г = 0 подин­ тегральная функция (1.13) терпит разрыв (\|/ = «>), формально указывая на на­ личие в ней источника колебаний с давлением р. Поэтому Sxразделяем на сфе­ ру с радиусом-вектором г'—>0 (поверхность S/, нормаль ri) и охватывающую ее замкнутую поверхность S с радиусом-вектором г и нормалью и к ее элемен­ ту dS.

Вследствие этого из формулы Г. Грина получим уравнение:

Д О М р и), -(рЭ\|//Эи*)]</5+ Jjlim[ у ' (др/дп) йВ"] -

s

S'

(1.15)

 

S'

где il/7 = (/) 4 ехр(-/Лг') — потенциал на элементах dS' сферы, dS' =(rfdCl9Cl— телесный угол (от 0 до 4л), ps — давление на элементе поверхности сферы. Так как \\f'dS' = /dClexp(-ik/)9то 2-й интеграл в выражении (1.15) при г = 0 обращается в нуль. Определение градиента (Э\|//Эл)5,потенциала \|/ на сфере S' по внутренней нормали ri (встречно радиус-вектору г вышеуказанного источ­ ника р) приводит к следующему результату:

(Э\р/Эя)5, dS ' = —(Э\|//Эг ') (г f)2dCl = (ihr'+1)dClexp ( - ikr')

и, с учетомps,~p при /

= 0, соответственно к выражению для 3-го интеграла в

(1.15):

 

 

jpd£l =4np.

S

0

Подстановка в (1.15) этого выражения, градиентов давления psи потенциала

\|f, составляющих для коротковолновых излучателей: bpjbn = ikps и ду/дп = = (дг/дп)(ду/дг) = z r l(dy/dr)9с учетом На|/ = -r~ldy(n\f)/dr, дает волновой интег­

рал Г. Гельмгольца - Н. Бражникова [46, 63]:

Преобразуем представленный в выражении (1.16) интеграл с учетом, согласно

(1.13), zr~'d\\i / d r - -ik^z / г1 + z / (ikr3)^exp(-ikr) и равномерности распреде­

ления на поверхности S ее колебательной скорости, как и давленияps, к следу­ ющему виду

р = — ■JJ/fc(r"* + zr 2 + z(Jkr3 )"1)exp( - ikr)dS s

Представим это выражение в цилиндрических координатах: (0, ра, аи) для точки излученияЛ на элементе акустического преобразователя, ось которо­ го (рис. 1.3) совмещена с координатной осью Z, и , Ь, аь) для точки приема В. Отсчет а а, ah ведем относительно параллельных прямых X и Хг Обозначив dS = padpadaa и пределы интегрирования по раот 0 до а и по а в от 0 до 2л, получим:

P = Y~ )ikp„dp„ j(r" ' +zr’2 + z(ikr3y')exp(-ikr)d<xa,

(1 .1 7 )

471 О

О

 

где г = [г,2 - 26pf/cos(ocrt -

а,,)]°’5, г,2 = z2 + b1 +р2 ; b — расстояние от точки В

определения давления до оси ПЭ.

Ограничение 2-мя членами разложения г в ряд Маклорена, подстановка ве­

личины , получаемой из дифференцирования

в (1.17), и замена пределов ин­

тегрирования по с

на пределы интегрирования от г. = (z2 + 62)0,5до г = . 2 +

2 0 5

^

 

 

 

 

 

®

o n

+ а ) ’ по переменной г, для преобразуют (1.17) к виду:

 

р

1

"

 

 

 

 

 

 

= —

| /7:ехр( -

ikrx)drxJexp [/ecos(a<#- а л)]х

 

 

Ps

r

 

 

о

 

 

 

 

(1.18)

 

 

L

 

 

 

 

 

 

 

 

 

П

e

, 2 Г

з

'l

>da.

 

 

х<П + —

1 + —

+ —

1 + - 2 + —

cos(aa - a t )

 

 

 

rx V

i k r )

kr, .

Ъ К

ikr*J

 

Для фиксированного af)при daa= d(ocfj - a /f), функции Бесселя имеют вид:

 

 

 

J

 

 

 

 

 

 

У0(в) = —

I exp[zecos(aa

 

 

и

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

J ,(e )= —

[exp[iecos(aa -<Xj)]cos(afl- a b)daa .

 

 

 

2iJt ;

 

 

 

 

 

При этом из (1.18) следует:

 

 

 

 

 

 

 

 

l + i |

2 + - 3

Jx(e)exp( - ih\ )drx

(1.19)

 

 

-PsM - Îh -r

 

ib\

 

 

 

где 4 ~

i U

1 + i

exp( -

ikrx)

 

 

 

 

2 ■'Эг,

 

 

 

 

 

 

 

rb

1

 

 

 

 

 

 

Интегрирование Л, по частям, с учетом

— [J 0 (e)] =

e k2b2

J,(e ), дает

 

 

 

 

 

 

dr

 

 

выражение:

 

 

 

 

 

 

 

A = 2r,

J о (e)exP( -

) ,

Д

 

(e)exp(—/AT,)dr, +

 

r , j 2 e i

 

 

 

 

 

 

-

f— 1 + - Jj (e)exp( -

/fo;)c/r,.

 

 

 

 

 

 

M

ny

 

 

Третье слагаемое этого выражения сокращает второе слагаемое в (1.19) до малой величины, которой можно пренебречь. Поэтому:

1

Z

P = Ps ^ h ~ Jo(е)ехр( - Her)

— [exp(-i*r, )]<*;.

 

Давление ps на колеблющейся поверхности ПП имеет амплитуду р0и зави­ симость от частоты со и времени V.ps =pQехр(/со/). (1.20)

В первом приближении для Еа= кЪа!гаполученная формулар имеет вид:

р = 0,5р0[(1+z / гь)ехр( - ikrb) - (1+ z / ra)J0(га)] ехр(гсм),

(1.21)

Распределение осевого давления в звукопроводе

Z/ZB

0,1

0,22

0,42

0,96

1,57

3,14

(ро/тс

4

3

2

1

0,63

0,32

№ У р о \

0,38

1,74

0,13

1,96

1,65

0,95

которая на оси ЯЭ, когда b = sa= 0, J0(ea) = 1 и rb= z, r= rm= (z2 + я2)0,5, для осевого давления принимает вид формулы Н. Бражникова [46, 63]:

P(z) = A [ l _ 0.5(1 + z / r(J1)ехр (-кр0 )] ехр [/(со/ -к£)\ ,

(1.22)

где ф0= k(rm-z). Распределениеp(z) в ЗПвышеуказанного расходомера по (1.22) от приведенного расстояния z/zBи фазы <р0 приведено в табл. 1.2.

В ближней зоне УЗ-поля протяженностью, согласно (1.22), давления по мере удаления от ПЭ в общем случае могут изменяться в пределах:

1,5р0 < \р\тах <2р0 и 0<|p|min < 0,5р0.Чередующиеся экстремумы давления в сре­

де излучения имеют место при фазе давления <р0, кратной 71 (минимумы — четное тс, максимумы — нечетное).

УЗ-давление в поле ПЭ есть результатом интерференции как бы 2-х волн, излучаемых центром {га= 0, гх= гь) и периметром (ро = я, г{= га) ПП. Вслед­ ствие различия давлений 2-х волн по фазе, амплитуде и знакопеременности 70(е) в (1.21) при b > 0 также имеют место экстремумы давления. При этом каждый минимум давления на оси оказывается окруженным кольцевым мак­ симумом вне оси. Аналогично, каждый осевой максимум (кроме последнего при z =zB), окружен кольцевым минимумом.

Полученные выражения компонент эквивалентного контура, электрическо­ го импеданса, давления излучения и его распределения в поле ПП использу­ ются при проектировании акустических устройств для систем УЗКР ТП.

Глава 2. СТРУКТУРА ПОЛЯ В СИСТЕМЕ УЗРК ТП

Рассмотрена структура поля (используемого в УЗ-контроле жидкостных по­ токов [11-20], газообразных сред [64] и поверхностной плотности [46-49, 65]) коротковолнового УЗ-излучения непрерывных гармонических колебаний и длительных импульсов. Уточнена формула давления излучения в параксиаль­ ной области поля и на поверхности акустического цилиндра излучения. При­ водятся дифракционные поправки на УЗ-давление и скорость волны.

2.1. Концепции распределения давления в ультразвуковом поле

Поле ПП описывается системой волновых уравнений. Их общее решение для упругого потенциала ср поля по данным Дж. Рэлея [66] первым выполнил Г. Гельмгольц с помощью интегральной теоремы Г. Грина в функции потенци­ ала у точечного источника.

Такое решение системы уравнений для потенциала ф (и соответствующего ему давления р) представляет собой совокупность воздействий, распределен­ ных по колеблющейся поверхности S излучателя: точечных источников (мо­ нополей) с потенциалом \|/ и дипольных источников с потенциалом ф колеба­ тельной скорости поверхности S .

На практике применяют три концепции распределения давления с представ­ лением его в функции: 1) одних монополей, 2) одних диполей, 3) совокупнос­ ти монополей и диполей.

1. Концепция одних монополей. Согласно этой концепции [62,63] давление определяется интегралом Дж. Рэлея:

(2. 1)

Вывод (2.1), как упрощенного решения системы волновых уравнений на базе теоремы Г. Грина, был выполнен Е. Скучиком в 1954 г. и С. Ржевкиным в 1959 г. [67]. Из интегрирования (2.1) следует используемая до настоящего времени формула Г. Штенцеля [41,42,68] для давления p(z) на оси ПП в виде (1.14), где

р0exp(ffltf) = рсх>п

(2.2)

представляет собой начальное давление излучения коротковолновым ПП не­ прерывных колебаний или длительных (в сравнении с периодом колебаний) импульсов колебаний в функции колебательной скорости мп ПП.

Выражение (1.14) может быть записано в следующем виде:

где фаза давления cp0 определяется выражением

% = K r m- z ) .

(2.4)

2. Концепция одних диполей. Эта концепция базируется на представлении потенциала <р поля (или давления) в функции только одних диполей, выведена Е.Скучиком из общего решения системы волновых уравнений (также на базе теоремы Г. Грина) в виде интегрального выражения:

Ф=

Цф^г-'йВЭф/Э/-

(2.5)

 

2к s

 

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

(2.6)

Для вывода выражения давления на оси поршневого ПП начало координат (ра, 2, а) располагаем в центре ПП. При этом z = z2 + рв2. Из дифференцирова­ ния этого выражения следует:

rdr = padpa и dS = padpada = rdrda.

(2.7)

Пределы изменения независимых переменных составляют: для а от 0 до 2я, a для г от z до гти тогда выражение (2.6) записывается в виде:

Р(г) = ~

Tps ^ Э г

|Э а

(2 .8)

2 я

Эг

0J

 

или, с учетом значения psиз (1.20):

1 - — е х р ( - г ф 0) expi((ùt-fo)

(2.9)

. 2»

 

3. Концепция совокупности монополей и диполей. Решение системы вол­ новых уравнений на основе указанной выше теоремы Г. Грина в виде зависи­ мости давления от совокупности монополей и диполей для коротковолновых

ПП определяется [46, 63] интегралом Г. Гельмгольца - Н. Бражникова (1.16), которое может быть представлено в виде:

+

(2.10)

При указанном выше использовании координат из интегрирования (2.10) с учетом у из (1.13) и ps из (1.20) следует формула Н. Бражникова [46, 61] для давления на оси ПП:

Р(г) = />0[ l - ° ,5 (1 + z / rm)ехр( - г<р0)] ехр/(со/ - kz).

(2.11)

Осевое давление в ближней зоне поля ПП в зависимости от относительного расстояния z/a испытывает ряд амплитудных экстремумов при кратной я фазе <р0. Расчетный параметр z/ruсвязан с z/a соотношением

z / r m= [l + ( z / a ) 2] 5z / a

(2.12)

При вычислении экстремумов (фаза ср0 тогда кратна к) расстояния z/a, на которых они находятся от поверхности ПП с известным волновым парамет­ ром ка, определяются из выражения (2.4), записываемого в виде:

%/ка = (rm- z)/a = (1- z2/a 2)05 - z/a,

откуда

2z/ а = te(p0 1- (te )"1ф0.

Т а б л и ц а 2.1 Расчетные зависимости модуля отношения давлений р/р0по формулам концепций 1, 2

и 3 для ПП с ка = 16 на оси = 0) и в параксиальной области 0,1а) поля для расстояний zJa

 

Приb == 0 поконцепции(формуле)

Приb = 0,1 а поконцепции(формуле)

da

1

2

з

1

3

 

Г. Штенцеля

Е. Скучика

Н. Бражникова(2.11)

А. Шоха

 

 

(2.3)

(2.9)

(2.13)

(2.18)

 

 

0,018

2

1,018

1,509

 

 

0,244

0

0,764

0,382

0,519

0,702

0,554

2

1,484

1,742

1,567

1,421

1,077

0

0,268

0,134

0,275

0,372

2,448

2

1,926

1,963

1,910

1,877

4

1,666

1,742

1,654

1,635

1,623

6

1,230

1,220

1,225

1,219

1,215

8

0,955

0,953

0,954

0,950

0,949

Результаты расчета экстремумов модулей осевого давления \p(z)\/pQв функ­ ции da по формулам концепций 1, 2 и 3 в пределах ближней зоны поля ПП с ка = 16 приведены в левой части табл. 2.1.

Максимумы осевого давления возникают при <р0, равной нечетному числу к, причем наиболее удаленному от ПП данного типа максимуму (при ка = 16) соответствует окончание ближней зоны. Ближний к ПП максимум (при ср = я) находится на расстоянии 0,02zla от него. При фазе ср0, равной четному числу к, возникают минимумы давления. Наиболее удаленный от ПП минимум (при ср0 = 2я) находится от ПП данного типа на расстоянии 1,08а, а ближайший к ПП — на расстоянии 0,24а.

По концепции 1 все минимумы осевого давления равны нулю, а максимумы равны 2р0. Это не соответствует экспериментальным данным. По концепции 3 на оси ПП минимумы давления по мере роста da не принимают значения 0: для ПП с ка = 16 они спадают от 0,38/?0 при ф0 = 4л до 0,13/?0 при ф0 = 2л. Максимумы с ростом увеличиваются от 1,51р0при ф0 = 5л вблизи ПП до 1,96р0 при ф0 = л в конце ближней зоны. Здесь расчеты соответствуют эксперимен­ тальным данным [69] на частоте 350 кГц (ка = 59) и нашим исследованиям жидкостного звукопровода фазового измерителя скорости потока на частоте 500 кГц (ка= 16) при их сопоставлении по обобщенной координате .

По концепции 2 минимумы давления завышают минимумы концепции 3, а максимумы — занижают максимумы концепции 3.

Выражения концепций 1 и 2 являются упрощенными решениями системы волновых уравнений. Такое определение могло бы быть поставлено под со­ мнение, если бы полученные из одной системы волновых уравнений выраже­ ния осевых давлений по (2.9) концепции 2 и по (2.3) концепции 1 были бы равны друг другу для любого da. В действительности же существует их раз­ ность, определяемая выражением:

p(z)2 - p(z\ =P o ( l - Z/гт)expi((ùt-krm).

Разность модулей этих давлений имеет величину, по амплитуде убывающую (для ка = 16) от 0,98р0вблизи ПП до 0,07р0 в конце ближней зоны.

Следовательно, только интеграл и формулу концепции 3 можно считать полными решениями системы волновых уравнений соответственно давления рьв точке В поля и p(z) на оси коротковолнового ПП. Использование концеп­ ций 1 и 2 приемлемо лишь в дальней зоне, а концепции 3 — во всех зонах поля

ПП.

2.2. Параксиальная область поля

Из интегрирования (2Л) следует известная [62] формула А. Шоха для давле­ ния в точке, отстоящей на расстоянии b от оси ПП и расстоянии z от колеблю­ щейся поверхности ПП:

P = Ps [ехР( - ikz)- J o (гь)ехр( - ikrm)],

(2.13)

где

 

гь = Ька/гт.

(2.14)

Для определения давления в параксиальной области (0 < b < 0,2а) поля нача­ ло системы координат (z, р^, а) разместим (рис. 2.1) в проекции Вхточки при­ ема (отстоящей на расстоянии b от оси ПП) на плоскость ПП.

Параметр раизменяется от 0 до рА, зависящего от b и определяемого рассто­ янием от центра В{системы координат до точки К в центре элемента dSk, сколь-

Рис. 2.1. Схема расчета УЗ-давления в параксиальной области поля УЗ-излучателя