книги / Ультразвуковой контроль и регулирование технологических процессов
..pdfнию уровня излучения и отношения информативный сигнал/акустические помехи.
Полученные расчетные формулы применимы и для поперечных 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 |
" |
|
2л |
|
|
|
|
|
— |
= — |
| /7:ехр( - |
ikrx)drxJexp [/ecos(a<#- а л)]х |
|
|
||||
Ps |
4л 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 |
2я |
|
|
|
|
|
|
У0(в) = — |
I exp[zecos(aa |
|
|
|||
и |
|
|
|
|
|
|
|
|
|
|
1 |
2П |
|
|
|
|
|
|
|
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, )]<*;. |
2П |
|
Давление 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. Схема расчета УЗ-давления в параксиальной области поля УЗ-излучателя