книги / Ультразвуковой контроль и регулирование технологических процессов
..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. Схема расчета УЗ-давления в параксиальной области поля УЗ-излучателя
