CC BY
24
Ф1З1КА
УДК 534.8:535.5
Г. В. Кулак, Т. В. Николаенко, В. Г. Гуделев
ОПТИКО-АКУСТИЧЕСКОЕ ВОЗБУЖДЕНИЕ ГИПЕРЗВУКА В ТВЕРДЫХ ТЕЛАХ В РЕЖИМЕ ЛАЗЕРНОГО ИСПАРЕНИЯ
Введение
Среди источников ультразвука в твердых телах оптико-акустические имеют ряд преимуществ: отсутствие контакта со средой, возможность легкого изменения геометрических параметров оптико-акустической антенны и диапазона возбуждаемых частот звука, создание источников, движущихся с определенной скоростью [1]. В работе [2] получено и решено общее уравнение фотоакустического эффекта в твердых телах. Экспериментальному исследованию лазерной генерации звука в нелинейном режиме испарения металла посвящена работа [3]. В [4] изучена оптико-акустическая диагностика изотропных твердых тел в условиях термоупругого
эффекта. При интенсивностях оптической волны 10 > 107 -108 Вт / см2 проявляются нелинейные эффекты при лазерном возбуждении ультразвука [1, 5].
В настоящей работе с использованием формализма тензорной функции Грина для векторного уравнения Ламе [6] исследованы особенности генерации продольных и сдвиговых УЗ волн в твердых телах, включая металлы, в условиях лазерного испарения.
Основная часть
Сравнительно просто удается выполнить приближенные расчеты генерации звука при очень больших интенсивностях света, когда происходит оптический пробой продуктов испарения и образования плазменного факела [3]. В этом случае в основе явления лежит процесс поглощения света в плазме и образования бегущей навстречу лазерному лучу детонационной (ударной) волны, за фронтом которой образуется область повышенного давления, действующего на поверхность металла.
Давление, оказываемое на поверхность в условиях плазменного режима, дается соотношением [3]
р f = р V2 » ЕЯ - » Е (Е / р') ~3/5 г-6/5, (1)
где рf - давление позади фронта ударной волны, р' - плотность окружающей среды, Е - объемная плотность энергии УЗ волны, Vf - скорость ударной волны, Яf - радиус ударной волны, возникающей в некоторой точке на расстоянии ё от поверхности мишени. Характерно, что величина импульса силы практически не зависит от рода материала (использовались алюминий, сталь, вольфрам). При этом давление достигает р{ ~ 100 МПа [1, 3].
Импульс УЗ смещений й (г, г) имеет ширину полосы ДО, которая определяется
длительностью падающего светового импульса т ~ 1 / ДО.. Предположим, что центральная частота УЗ импульса равна О. Тогда уравнение Ламе, описывающее импульсное возбуждение ультразвука за счет распределенной объемной силы f (г') на поверхности металла, вызванной упругими напряжениями, имеет вид [6, 7]
Ей = -}5(г - г'), (3)
где 5 (х) - дельта-функция. Предполагается в дальнейшем, что в рамках линейной теории вектор УЗ смещений и ~ ехр(-Юг).
Оператор дается соотношением:
2
= р. Д + (1 + ¡л)graddiv + р О , (4)
где 1, л - постоянные Ламе.
Решение стационарного уравнения (3) имеет вид [6, 7]
u, (r) = \ G. (r - r'|) f, (r') dr',
(5)
где О.. - компоненты тензорной функции Грина, г' - произвольные внешние параметры.
У
Компоненты тензорной функции Грина даются соотношениями:
G.. =
г,
4ppW
iktR
8 ..k. г1 t
R dx. dx.
г,
ik¡R iktR
e e
R
R
(6)
Здесь введены обозначения: R = | r - r k{ = W / u{, kt = W /ut, где u{ = д/ (l + 2m) / p
ut = д/ m / p - фазовые скорости продольной и сдвиговой УЗ волн (p - плотность материала
звукопровода); произвольные внешние параметры r' связаны с областью лазерного воздействия на поверхности металла.
В режиме лазерного испарения на свободной поверхности металла образуется область радиусом а с резкими границами. При этом импульс давления возбуждает УЗ волну с цилиндрическим фронтом. Данный режим возбуждения ультразвука аналогичен возбуждению УЗ волны поршневым пьезоэлектрическим преобразователем [8]. Смещение в продольной УЗ волне находится с использованием соотношений (5), (6). При этом в соответствии с результатами работы [8] следует положить R = r - r'sinв • cos(y - y'), где r, y -цилиндрические координаты в точке наблюдения УЗ поля, r', y' - цилиндрические координаты в плоскости лазерного воздействия X' Y' (рис. 1).
2
9
1
e
Рис. 1. Схема возбуждения ультразвука для области воздействия круглой формы (а) и схема расположения цилиндрических координат (б)
Для волнового вектора ультразвуковой волны, лежащего в плоскости Х2, смещение и 3 дается соотношением
22
и 3( г, в) = и у к{а) 008 в
1
+ — и у ^ (к(а)
2
J1 (ка бшв) (ка бшв)
J1 (кга бш в) (к а бшв)
шгг п \2 е - (к, а)
22
е ^ + и у ^(к,а) бш в
J1 (к а бшв) (к а бшв)
J1 (к а бшв) (к а бшв)
е ^ +
(7)
где и у = ру /2рО г, J1 (х) - функция Бесселя первого порядка.
Смещение и1 для сдвиговой составляющей УЗ волны находим из соотношения
22
и,(г ,в) = и у ка) бш в'
+ и у { (к,а )2 соб 2 в'
J1 (к(а бш в')
(кга бш в') J1 (к, а бш в')
е г > +
(к ¡а бш в')
е ' > +
+ — игу ^ (к(а)
J1 (кга бш в') (к а бшв')
е г - (к а)
J1 (ка бшв') (к, а бш в')
(8)
где иу » и у = ру / 2рО г, в' = в - вг.
На рис. 2 представлены зависимости амплитуды УЗ смещения |и3| продольной УЗ волны, рассчитанные по формуле (7), от полярного угла в при различных радиусах области возбуждения поверхности круглой формы.
¡кг г
2
1
г
е
1 г - у»
_V , у- _
Г» 0.5 0.75 1.0 0. град
Рис. 2. Зависимость амплитуды смещений продольной УЗ волны |и3| от полярного угла в для области возбуждения круглой формы при различных значениях параметра радиуса а: 1 (1) 1,5 (2), 3 (3) мм (т = 10-9с, р{ = 100 МПа).
Из рис. 2 следует, что амплитуда смещения и3 увеличивается с увеличением размеров области возбуждения ультразвука. Величина смещения и3 ~ 1-4,5 нм.
Предположим, что световой импульс падает на свободную поверхность твердого тела в форме прямоугольника со сторонами, ориентированными вдоль осей X' и У' соответственно: -И < у' < И, - 1< х' <1 (рис. 3).
Рис. 3. Схема возбуждения ультразвука полоской прямоугольной формы
Для дифракции в дальней зоне следует положить
R = r - (х 'cosa + y 'cos b + z 'cos/).
(9)
где r' = (х',y', z'), cosa = х/r, cos b = y /r, cos у = z /r, z' << х', y', r = | r | - расстояние от центра лазерного воздействия до точки наблюдения УЗ поля. Предполагается, что распределение давления однородно в пределах области воздействия. Подставив выражения (6), (7) в (5) и выполнив интегрирование по области воздействия, получим выражение для проекции вектора смещений продольной УЗ волны u3, лежащей в плоскости XZ:
2 2
u3(r,в) = uf(2lh/p)ík cos в
sin(kh sine)
(kh sine)
iktr
+ Uf (2lh / p)í kl sin в
sin(klh sin в) (k¡h sin в)
+ u lf (lh / р^ш(2в^ kt
sin(kth sin в) (kh sin в)
iktr i 2 в - kl
sin(klh sin в) (klh sin в)
(10)
Проекцию вектора смещений и1 на плоскость XX для сдвиговой УЗ волны находим из соотношения
2
u1(r,в) = uf(2lh/p)íkt cos в'
sin(kth sin в') (kh sin в)
e'ktr > +
2
+ uf (2lh/p)í k1 sin в'
+ uf (lh / p)sin(2в')í kt
sin(klh sin в') (klh sinв')
sin(kth sin в') (kh sin в')
e ' ^ +
iktr , 2 e - kl
sin(klh sin в') (k,h sin в')
(11)
На рис. 4 представлена зависимость амплитуды смещения |и3| продольной УЗ волны, рассчитанная на основании выражения (10), от угла рассеяния в по отношению к нормали плоской поверхности металла (железо).
+
iklr
+
e
iklr
e
2
ik r
iklr
e
u .„lim "lti
il 1
ч Фаз
Рис. 4. Зависимость амплитуды смещений продольной УЗ волны |и3| от полярного угла в для полоски прямоугольной формы при различных значениях параметра д: 100 (1), 50 (2), 25 (3) (т = 10-9с, р = 100 МПа).
Из рис. 4 следует, что максимальное смещение происходит в направлении, близком к в = 0 . С увеличением ширины полоски к = дЛ (Л - длина УЗ волны, д - параметр) амплитуда смещений увеличивается и положение первого минимума (и3 = 0) смещается в сторону больших углов в . Узость диаграммы направленности (до первого минимума) объясняется высокими частотами возбуждаемого ультразвука /0 »1ГГц.
Величина смещения рассеянной УЗ волны дается очевидным соотношением [8]
m
'sinô'-Vn '2 - cos2 в
m 'sine '+
cos2 в '
4x '
cos(e ')
(12)
где и - амплитуда смещений УЗ волны.
Для оптико-акустического излучателя круглой формы амплитуду рассеянной УЗ волны |их | находим под углом в' = 8,50 к оси 2, для излучателя прямоугольной формы в' = 1,660. С использованием выражений (7) и (12) для г » x, » 1см и плотного рассеивателя [см. 4] получаем соответственно: ~0,4 нм (круг), ~0,85 нм (прямоугольник).
Заключение
Полученные результаты показывают, что рассеянные ультразвуковые волны могут быть обнаружены методом оптического гетеродинирования и методом РИой-ЕМР [9, 10]. Применение системы сфазированных оптико-акустических источников гиперзвука позволит значительно увеличить амплитуду возбужденных волн и, тем самым, увеличить глубину расположения исследуемых дефектов. Система датчиков, расположенных на поверхности материала, позволит определить положение дефекта в пространстве.
Литература
1. Лямшев, Л. М. Оптико-акустические источники звука / Л. М. Лямшев // Успехи физ. наук. - 1981. -Т. 135, вып. 4. - С. 636-669.
2. Lin, Gu. Theory of the photoacoustical effect in condenced medium / Lin Gu // Appl. Opt. - 1982. -V. 21, № 5. - P. 955-960.
3. Pirri, A. N. Momentum transfer and plasma formation above a surface with and highe-power CO2 laser / A. N. Pirri, R. Schlier, D. Northam // Appl. Phys. Lett. - 1972. - V. 21, № 3. - P. 79-81.
4. Кулак, Г. В. Особенности оптико-акустической диагностики твердых тел в условиях термоупругости / Г. В. Кулак, Т. В. Николаенко // Весшк Мазырскага дзяржаунага педагапчнага ушверспэта. - 2006. - № 1. -С. 25-29.
5. Experimental study of laser-generated shear waves using interferometry / S. Y. Zang, M. Paul, S. Fassbendtr, U. Schleichert and W. Arnold // Res. Nondestr. Eval. - 1990. - V. 2. - P. 143-155.
6. Ландау, Л. Д. Теория упругости / Л. Д. Ландау, Е. М. Лифшиц. - М. : Наука, 1987. - 244 с.
7. Муратиков, К. Л. О формировании изображений неоднородных объектов лазерным фотоакустическим методом / К. Л. Муратиков // Письма в ЖТФ. - 2004. - Т. 30, вып. 22. - С. 58-64.
8. Кайно, Г. Акустические волны. Устройства, визуализация и аналоговая обработка сигналов / Г. Кайно. - М. : Мир, 1990. - 652 с.
9. Paul, M. Interferometric detection of ultrasound at rough surfaces using optical phase conjugation / M. Paul, B. Betz, W. Arnold // Appl. Phys. Lett. - 1987. - V. 22, № 1. - P. 1569-1571.
10. Petrov, M. High-frequensy branch of space-charge waves in photorefractive crystals / M. Petrov, V. Bryksin, A. Emgrunt [et. all.] // J. Opt. Soc. Am. - 2005. - V. 22, № 7. - P. 1529-1537.
Summary
Optical - acoustical ultrasound sources of circle form and rectangular form, disposed on the surface of solid state body is considered. Angular distribution of the longitudinal and shear ultrasonic waves in the laser ablation regime from polar angle relatively to the normal of plane boundary of the solid state body is investigated. The possibility of the effective ultrasound excitation with narrow angle distribution for short laser pulses with a 1 ns duration are established. It is shown that the scattered ultrasonic waves, excited from the metal surface may be detected by the optical methods.
Поступила в редакцию 10.11.06.
УДК 535.42
В. В. Шепелевич, А. Е. Загорский, Р. Коваршик, А. Кислинг, В. Матусевич
ИССЛЕДОВАНИЕ РАСПРОСТРАНЕНИЯ ДВУМЕРНЫХ СВЕТОВЫХ ПУЧКОВ В КРИСТАЛЛЕ В1128Ю20 СРЕЗА (ТТ0)
Введение
Уже более десятилетия внимание исследователей привлекают процессы самофокусировки двумерных световых пучков в фоторефрактивных кристаллах [1-2]. Получаемые в таких кристаллах двумерные солитоны [2-4] могут быть использованы для адресной локализации световых пучков, а также для их бездифракционного распространения в динамически формируемых волноводах. Известно, что в одноосных кристаллах типа 8БМ [1-4] можно учитывать лишь одну компоненту электрооптического тензора, пренебрегая остальными компонентами из-за их малости. В кубическом кристалле Б^28Ю20 электрооптический тензор имеет несколько одинаковых компонент [5]. Поэтому уравнения, описывающие распространение световых пучков в кубических кристаллах, должны быть несколько изменены. В связи со сложностью вычислений, возникающих при использовании подхода [1, 2] к описанию квазисолитоного распространения двумерных световых пучков, в некоторых работах [6, 7] применялись более простые, приближенные способы описания этого явления.
В данной статье используется модифицированная система скалярных дифференциальных уравнений, описывающая распространение двумерного гауссового светового пучка в фоторефрактивном, оптически активном кристалле Б^28Ю20 кристаллографического среза (1 10). При этом в системе дифференциальных уравнений учитывается распределение потенциала ф как
вдоль приложенного к кристаллу внешнего электрического поля E0, так и перпендикулярно ему. Потенциал ф вычисляется с использованием соотношений [1, 2, 6].
Основная часть
Для описания распространения двумерного светового пучка в кубических фоторефрактивных оптически активных кристаллах будем использовать следующую систему уравнений, полученную в параксиальном приближении на базе уравнений Максвелла и основных уравнений фоторефрактивного эффекта [8]:
дЛ 1 i— +-
Ад2 Л д2 ЛЛ
2 2 dz 2k0n0 удх ду
Kn,
3
(ЛгЕ) + ip[ez,Л] = 0, (1)
