Журнал вычислительной математики и математической физики, 2019, T. 59, № 4, стр. 597-610

Аналитическое решение задачи Лэмба в случае предельного значения коэффициента Пуассона

Х. Х. Ильясов 12*, А. В. Кравцов 12**, С. В. Кузнецов 12***, С. Я. Секерж-Зенькович 12****

1 Ин-т проблем механ.
117526 Москва, пр-т Вернадского, 101, Россия

2 МГУ, физ. ф-т
119992 Москва, Ленинские горы, Россия

* E-mail: ilyasov@ipmnet.ru
** E-mail: avkravtsow@rambler.ru
*** E-mail: kuzn-sergey@yandex.ru
**** E-mail: sekerzh@gmail.com

Поступила в редакцию 03.10.2018
После доработки 14.11.2018
Принята к публикации 14.11.2018

Полный текст (PDF)

Аннотация

Рассматривается задача Лэмба для силы, приложенной к границе упругого полупространства в случае, когда коэффициент Пуассона принимает предельное значение 1/2. Решение представляется в виде контурных интегралов, для которых приводятся асимптотические оценки при больших значениях полярного радиуса. Библ. 5.

Ключевые слова: упругая среда, уравнения Ламэ, коэффициент Пуассона, полюс Рэлея, асимптотические оценки контурных интегралов.

1. ПОСТАНОВКА ЗАДАЧИ. ИНТЕГРАЛЬНОЕ ПРЕДСТАВЛЕНИЕ РЕШЕНИЯ

Пусть упругая среда занимает полупространство. В случае малых относительных перемещений уравнения Ламэ имеют вид

$(\lambda + \mu )\operatorname{grad} \operatorname{div} u + \mu \Delta u = \rho \frac{{{{\partial }^{2}}u}}{{\partial {{t}^{2}}}}.$
Здесь $u$ – вектор перемещения, $\lambda $, $\mu $, $\rho $ – соответственно параметры Ламэ и плотность упругой среды. Пусть к свободной поверхности $S$ упругой среды приложена внешняя нагрузка $np$, где $n$ – вектор внешней нормали к $S$. Следуя [1], граничные условия на свободной поверхности зададим в виде
$2\mu \frac{{\partial u}}{{\partial n}} + \lambda n\operatorname{div} u + \rho gn(n,u) + \mu [n,\operatorname{rot} u] = np,$
где $g$ – ускорение силы тяжести (вектор силы тяжести противоположен вектору $n$). Считаем, что $p$ зависит от точки свободной поверхности и времени.

Согласно теории предельного перехода будем считать, что коэффициент Пуассона $\nu = \frac{\lambda }{{2(\lambda + \mu )}} \to \frac{1}{2} - 0$, а модуль Юнга $E = (3\lambda + 2\mu )(1 - 2\nu ) \to + 0$ так, что отношение $\frac{E}{{1 - 2\nu }}$ остается конечным, не равным нулю. Тем самым $\mu = \frac{E}{{2(1 + \nu )}} \to + 0$, что означает отсутствие в упругой среде волн сдвига. Тогда уравнения для перемещений и граничные условия на $S$ примут соответственно вид

$\lambda \operatorname{grad} \operatorname{div} u = \rho \frac{{{{\partial }^{2}}u}}{{\partial {{t}^{2}}}},\quad \lambda \operatorname{div} u + \rho g{{u}_{z}} = p.$

Представим вектор перемещений в виде $u = grad\varphi $ и введем цилиндрическую систему координат $(r,\theta ,z)$, для которой поверхность $S$ совпадает с плоскостью $z = 0$ и орт оси $z$ сонаправлен с вектором $n$. Считая, что внешняя нагрузка имеет вид $p(r,t) = {{p}_{0}}f(r)cos\omega t$, где $f(r)$ –заданная функция, допускающая разложение в интеграл Фурье–Бесселя, а ${{p}_{0}}$, $\omega $ – заданные постоянные величины, получаем скалярную задачу

$\frac{{{{\partial }^{2}}\varphi }}{{\partial {{t}^{2}}}} = {{a}^{2}}\Delta \varphi ,\quad r > 0,\quad z \leqslant 0,\quad t > 0,\quad a = \sqrt {\frac{\lambda }{\rho }} ,$
${{\left. {\rho \frac{{{{\partial }^{2}}\varphi }}{{\partial {{t}^{2}}}}} \right|}_{{z = 0}}} + \rho g{{\left. {\frac{{\partial \varphi }}{{\partial z}}} \right|}_{{z = 0}}} = {{p}_{0}}f(r)cos\omega t,\quad r > 0,\quad t > 0,$
которая должна быть дополнена граничными условиями на бесконечности. А именно, потребуем, чтобы функция $\varphi $ описывала уходящую в бесконечность волну и $\left| \varphi \right| \to + 0$ при ${{({{r}^{2}} + {{z}^{2}})}^{{1/2}}} \to + \infty $.

Так как нагрузка не зависит от угловой переменной $\theta $, то периодическое по времени решение будем искать в виде

$\varphi (r,z,t) = {\text{Re}}\{ \Phi (r,z){{{\text{e}}}^{{i\omega t}}}\} .$
Для комплексной амплитуды $\Phi (r,z)$ получим
$\Delta \Phi + \frac{{{{\omega }^{2}}}}{{{{a}^{2}}}}\Phi = 0,\quad r > 0,\quad z \leqslant 0,$
$\rho g{{\left. {\frac{{\partial \Phi }}{{\partial z}}} \right|}_{{z = 0}}} - {{\left. {\rho {{\omega }^{2}}\Phi } \right|}_{{z = 0}}} = {{p}_{0}}f(r),\quad r > 0,\quad \left| \Phi \right| \to + 0,\quad {{({{r}^{2}} + {{z}^{2}})}^{{1/2}}} \to + \infty ,\quad z \leqslant 0,$
$\frac{{\partial \Phi }}{{\partial r}} + i{{k}_{R}}\Phi = o\left( {\frac{1}{{\sqrt r }}} \right),\quad r \to + \infty ,\quad z \leqslant 0,\quad {{k}_{R}} = \frac{\omega }{a}\sqrt {1 + \frac{{{{\omega }^{2}}{{a}^{2}}}}{{{{g}^{2}}}}} .$
Величина ${{k}_{R}}$ носит название размерного полюса Рэлея.

Далее считаем, что функция $f(r)$ имеет вид

$f(r) = \frac{{{{l}^{3}}}}{{{{{({{r}^{2}} + {{l}^{2}})}}^{{3/2}}}}},\quad l > 0.$

Введем безразмерные переменные и безразмерную комплексную амплитуду по формулам

${{r}_{a}} = \frac{\omega }{a}r,\quad {{z}_{a}} = \frac{\omega }{a}z,\quad {{\Phi }^{{(a)}}} = \frac{{{{\omega }^{2}}}}{{{{a}^{2}}}}\Phi .$
Обозначим
${{l}_{a}} = \frac{\omega }{a}l,\quad \beta = \frac{g}{{\omega a}},\quad {{\kappa }_{R}} = \sqrt {1 + \frac{1}{{{{\beta }^{2}}}}} .$
Тогда для ${{\Phi }^{{(a)}}}$ получим задачу
$\Delta {{\Phi }^{{(a)}}} + {{\Phi }^{{(a)}}} = 0,\quad {{r}_{a}} > 0,\quad {{z}_{a}} \leqslant 0,$
(1)
$\beta {{\left. {\frac{{\partial {{\Phi }^{{(a)}}}}}{{\partial {{z}_{a}}}}} \right|}_{{{{z}_{a}} = 0}}} - {{\left. {{{\Phi }^{{(a)}}}} \right|}_{{{{z}_{a}} = 0}}} = \frac{{{{p}_{0}}}}{{\rho {{a}^{2}}}}f({{r}_{a}}),\quad {{r}_{a}} > 0,\quad \left| {{{\Phi }^{{(a)}}}} \right| \to + 0,\quad {{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty ,$
$\frac{{\partial {{\Phi }^{{(a)}}}}}{{\partial {{r}_{a}}}} + i{{\kappa }_{R}}{{\Phi }^{{(a)}}} = o\left( {\frac{1}{{\sqrt {{{r}_{a}}} }}} \right),\quad {{r}_{a}} \to + \infty ,\quad z \leqslant 0,$
где
$f({{r}_{a}}) = \frac{{l_{a}^{3}}}{{{{{(r_{a}^{2} + l_{a}^{2})}}^{{3/2}}}}}.$
Заметим, что для образа Фурье–Бесселя данной функции имеет место формула [2]
$\int\limits_0^{ + \infty } {{J}_{0}}(x{{r}_{a}})f({{r}_{a}}){{r}_{a}}d{{r}_{a}} = l_{a}^{2}{{{\text{e}}}^{{ - {{l}_{a}}x}}},\quad x > 0.$
Поэтому будем искать решение задачи (1) в виде комплексного интеграла
(2)
${{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}}) = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L \frac{{{{J}_{0}}(\kappa {{r}_{a}})\kappa }}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}}d\kappa ,$
где $L$ – кривая, состоящая из отрезка ${{\gamma }^{{(u)}}}$: $[0,1 - {{\rho }_{1}}]$ на верхнем берегу разреза $[ - 1,1]$, отрезка ${{\gamma }_{3}}$: $[1 + {{\rho }_{1}},{{\kappa }_{R}} - {{\rho }_{2}}]$, луча ${{\gamma }_{5}}$: $[{{\kappa }_{R}} + {{\rho }_{2}}, + \infty )$ и полуокружностей ${{C}_{1}}$: $\kappa = 1 + {{\rho }_{1}}{{{\text{e}}}^{{is}}}$, ${{C}_{2}}$: $\kappa = {{\kappa }_{R}} + {{\rho }_{2}}{{{\text{e}}}^{{is}}}$, $s \in [0,\pi ]$, пробегаемых по часовой стрелке, ${{J}_{0}}(\kappa {{r}_{a}})$ – функция Бесселя нулевого индекса аргумента $\kappa {{r}_{a}}$, а под $\sqrt {{{\kappa }^{2}} - 1} : = y(\kappa )$ понимается однозначная аналитическая ветвь двузначной функции $\sqrt {{{\kappa }^{2}} - 1} $, определенная на всей комплексной плоскости с разрезом $[ - 1,1]$ и положительная при действительном $\kappa > 1$. На ${{\gamma }^{{(u)}}}$ функция $y(\kappa )$ доопределена своими предельными значениями.

Покажем, что интеграл (2) удовлетворяет уравнению в (1) и граничному условию при $z = 0$. Для этого убедимся, что данный интеграл можно дважды дифференцировать по параметрам ${{r}_{a}}$ и ${{z}_{a}}$ под знаком интеграла. Представим данный интеграл в виде суммы пяти интегралов: по отрезкам, лучу и полуокружностям. Рассмотрим сначала интеграл по лучу ${{\gamma }_{5}}$ ($\kappa = x \geqslant {{\kappa }_{R}} + {{\rho }_{2}}$). Подынтегральная функция

$h({{r}_{a}},{{z}_{a}},x) = \frac{{{{J}_{0}}(x{{r}_{a}})x}}{{\beta \sqrt {{{x}^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{x}^{2}} - 1} - {{l}_{a}}x}}}$
непрерывна по совокупности переменных $({{r}_{a}},x)$ при фиксированном ${{z}_{a}}$ и по $({{z}_{a}},x)$ при фиксированном ${{r}_{a}}$ вместе с частными производными $\frac{{\partial h}}{{\partial {{r}_{a}}}}$, $\frac{{\partial h}}{{\partial {{z}_{a}}}}$, $\frac{{{{\partial }^{2}}h}}{{\partial r_{a}^{2}}}$, $\frac{{{{\partial }^{2}}h}}{{\partial z_{a}^{2}}}$. Интеграл вдоль ${{\gamma }_{5}}$ от $h({{r}_{a}},{{z}_{a}},x)$ сходится при любом фиксированном ${{r}_{a}} > 0$ по признаку сравнения
$\mathop {lim}\limits_{x \to + \infty } \left| {h({{r}_{a}},{{z}_{a}},x)} \right|{{x}^{2}} = 0.$
Интегралы вдоль ${{\gamma }_{5}}$ от $\frac{{\partial h}}{{\partial {{r}_{a}}}}$, $\frac{{\partial h}}{{\partial {{z}_{a}}}}$ сходятся равномерно соответственно по ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$ по признаку Вейерштрасса. Это следует из оценок

$\left| {\frac{{\partial h}}{{\partial {{r}_{a}}}}} \right| = \frac{{\left| {{{J}_{1}}(x{{r}_{a}})} \right|{{x}^{2}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{x}^{2}} - 1} - {{l}_{a}}x}}} < \frac{{{{x}^{2}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}},$
$\left| {\frac{{\partial h}}{{\partial {{z}_{a}}}}} \right| = \frac{{\left| {{{J}_{0}}(x{{r}_{a}})} \right|x\sqrt {{{x}^{2}} - 1} }}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{x}^{2}} - 1} - {{l}_{a}}x}}} < \frac{{x\sqrt {{{x}^{2}} - 1} }}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}.$

Покажем, что интегралы от $\frac{{{{\partial }^{2}}h}}{{\partial r_{a}^{2}}}$, $\frac{{{{\partial }^{2}}h}}{{\partial z_{a}^{2}}}$ вдоль ${{\gamma }_{5}}$ сходятся равномерно соответственно по ${{r}_{a}} \geqslant \delta $, ${{z}_{a}} \leqslant 0$, где $\delta $ – любое положительное число. Имеем

$\begin{gathered} \left| {\frac{{{{\partial }^{2}}h}}{{\partial r_{a}^{2}}}} \right| \leqslant \frac{{\left| {J_{1}^{'}(x{{r}_{a}})} \right|{{x}^{3}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}{{{\text{e}}}^{{ - {{l}_{a}}x}}} \leqslant \left\{ {\left| {{{J}_{0}}(x{{r}_{a}})} \right| + \frac{{\left| {{{J}_{1}}(x{{r}_{a}})} \right|}}{{x{{r}_{a}}}}} \right\}\frac{{{{x}^{3}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}} < \\ < \;\left( {1 + \frac{1}{{x{{r}_{a}}}}} \right)\frac{{{{x}^{3}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}} \leqslant \left( {1 + \frac{1}{{x\delta }}} \right)\frac{{{{x}^{3}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}, \\ \end{gathered} $
$\left| {\frac{{{{\partial }^{2}}h}}{{\partial z_{a}^{2}}}} \right| \leqslant \frac{{\left| {{{J}_{0}}(x{{r}_{a}})} \right|x({{x}^{2}} - 1)}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}{{{\text{e}}}^{{ - {{l}_{a}}x}}} < \frac{{{{x}^{3}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{\left| {\beta \sqrt {{{x}^{2}} - 1} - 1} \right|}}.$
Мы воспользовались соотношениями [2] $J_{0}^{'}(x) = - {{J}_{1}}(x)$, $J_{1}^{'}(x) = {{J}_{0}}(x) - \frac{1}{x}{{J}_{1}}(x)$.

Таким образом, интеграл вдоль ${{\gamma }_{5}}$ от $h({{r}_{a}},{{z}_{a}},x)$ имеет непрерывные частные производные по ${{r}_{a}}$, ${{z}_{a}}$ первого и второго порядков и они равны интегралам от соответствующих частных производных функции $h({{r}_{a}},{{z}_{a}},x)$.

Рассмотрим функцию трех комплексных переменных

$h(\tilde {r},\tilde {z},\kappa ) = \frac{{{{J}_{0}}(\kappa \tilde {r})\kappa }}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{\tilde {z}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}},$
где $\tilde {r}$, $\tilde {z}$ изменяются на всей комплексной плоскости $С $, а $\kappa \in C_{k}^{ - }$ $(k = 1,2)$ (индекс $a$ для краткости опускаем). Она аналитична по переменной $\tilde {r}$ на всей комплексной плоскости при фиксированных $\tilde {z} \in \mathbb{C}$, $\kappa \in C_{k}^{ - }$ и по переменной $\tilde {z}$ на всей комплексной плоскости при фиксированных $\tilde {r} \in \mathbb{C}$, $\kappa \in C_{k}^{ - }$. То же самое относится к частным производным $\frac{{\partial h}}{{\partial{ \tilde {r}}}}$, $\frac{{\partial h}}{{\partial{ \tilde {z}}}}$, $\frac{{{{\partial }^{2}}h}}{{\partial {{{\tilde {r}}}^{2}}}}$, $\frac{{{{\partial }^{2}}h}}{{\partial {{{\tilde {z}}}^{2}}}}$. Кроме того, и сама функция $h(\tilde {r},\tilde {z},\kappa )$ и ее указанные частные производные непрерывны по совокупности двух переменных $(\tilde {r},\kappa )$, $(\tilde {z},\kappa )$ при фиксированной третьей. Поэтому интегралы вдоль $C_{k}^{ - }$ от $h(\tilde {r},\tilde {z},\kappa )$ аналитичны на всей комплексной плоскости по переменным $\tilde {r}$ и $\tilde {z}$, а частные производные первого и второго порядков от указанных интегралов равны интегралам от частных производных функции $h(\tilde {r},\tilde {z},\kappa )$ всюду на $\mathbb{C}$, а стало быть и при $\tilde {r} = {{r}_{a}} > 0$, $\tilde {z} = {{z}_{a}} \leqslant 0$.

Аналогичным образом доказывается, что интегралы вдоль ${{\gamma }^{{(u)}}}$, ${{\gamma }_{3}}$ также можно дифференцировать два раза по параметрам ${{r}_{a}} > 0$ и ${{z}_{a}} \leqslant 0$ под знаком интеграла. Таким образом,

$\begin{gathered} \Delta {{\Phi }^{{(a)}}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L \left\{ { - \frac{\kappa }{{{{r}_{a}}}}{{J}_{1}}(\kappa {{r}_{a}}) - {{\kappa }^{2}}{{J}_{0}}(\kappa {{r}_{a}})} \right. + \left. {\frac{\kappa }{{{{r}_{a}}}}{{J}_{1}}(\kappa {{r}_{a}}) + ({{\kappa }^{2}} - 1){{J}_{0}}(\kappa {{r}_{a}})} \right\}\frac{{\kappa {{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}}}}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}d\kappa = \\ = \; - {\kern 1pt} \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L {{J}_{0}}(\kappa {{r}_{a}})\frac{{\kappa {{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}}}}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}d\kappa = - {{\Phi }^{{(a)}}}, \\ \end{gathered} $
т.е. $\Delta {{\Phi }^{{(a)}}} + {{\Phi }^{{(a)}}} = 0$.

Покажем, что граничное условие при ${{z}_{a}} = 0$ $({{r}_{a}} > 0)$ также удовлетворяется. Имеем

$\begin{gathered} \beta {{\left. {\frac{{\partial {{\Phi }^{{(a)}}}}}{{\partial {{z}_{a}}}}} \right|}_{{{{z}_{a}} = 0}}} - {{\left. {{{\Phi }^{{(a)}}}} \right|}_{{{{z}_{a}} = 0}}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L {{J}_{0}}(\kappa {{r}_{a}})\kappa {{{\text{e}}}^{{ - {{l}_{a}}\kappa }}}d\kappa = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\left\{ {\int\limits_{{{\gamma }^{{(u)}}}} \ldots + \int\limits_{C_{1}^{ - }} \ldots + \int\limits_{{{\gamma }_{3}}} \ldots + \int\limits_{C_{2}^{ - }} \ldots + \int\limits_{{{\gamma }_{5}}} \ldots } \right\} = \\ = \;\frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_0^{ + \infty } {{J}_{0}}(x{{r}_{a}})x{{{\text{e}}}^{{ - {{l}_{a}}x}}}dx = \frac{{{{p}_{0}}l_{a}^{3}}}{{\rho {{a}^{2}}}}\frac{1}{{{{{(r_{a}^{2} + l_{a}^{2})}}^{{3/2}}}}} = \frac{{{{p}_{0}}}}{{\rho {{a}^{2}}}}f({{r}_{a}}). \\ \end{gathered} $
Мы заменили интегралы по $C_{1}^{ - }$, $C_{2}^{ - }$ интегралами по отрезкам $[1 - {{\rho }_{1}},\;1 + {{\rho }_{1}}]$, $[{{\kappa }_{R}} - {{\rho }_{2}},\;{{\kappa }_{R}} + {{\rho }_{2}}]$ соответственно в силу аналитичности функции ${{J}_{0}}(\kappa {{r}_{a}})\kappa {{{\text{e}}}^{{ - {{l}_{a}}\kappa }}}$ по переменной $\kappa $, а затем воспользовались соотношением [2]

$\int\limits_0^{ + \infty } {{J}_{0}}(x{{r}_{a}})x{{{\text{e}}}^{{ - {{l}_{a}}x}}}dx = \frac{{{{l}_{a}}}}{{{{{(r_{a}^{2} + l_{a}^{2})}}^{{3/2}}}}}.$

Для безразмерных комплексных амплитуд перемещений $U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}})$ и $U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}})$ получаем выражения

(3)
$\begin{gathered} U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}}) = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L \frac{{{{J}_{0}}(\kappa {{r}_{a}})\kappa \sqrt {{{\kappa }^{2}} - 1} }}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}}d\kappa , \\ U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\int\limits_L \frac{{{{J}_{1}}(\kappa {{r}_{a}}){{\kappa }^{2}}}}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}}d\kappa . \\ \end{gathered} $

Заметим, что формальное интегральное представление решения задачи Лэмба в случае распределенной нагрузки для $0 < \nu < 1{\text{/}}2$ получено в [3]. Там же проведено сравнение аналитического и численного решений. В работе [4] начально-краевая задача Лэмба для полупространства решалась методом конечных элементов.

2. АСИМПТОТИЧЕСКИЕ ОЦЕНКИ ДЛЯ КОМПЛЕКСНЫХ АМПЛИТУД

Функция $y(\kappa )$ положительна при $\kappa = x > 1$. Тогда для ее значений при $\kappa = \pm i\sigma $, $\sigma > 0$, и на берегах разреза $[ - 1,1]$, т.е. при $\kappa = x \pm i0$, $ - 1 < x < 1$, будем соответственно иметь

$y( \pm i\sigma ) = \pm i\sqrt {{{\sigma }^{2}} + 1} ,\quad y(x \pm i0) = \pm i\sqrt {1 - {{x}^{2}}} .$
Проведем на комплексной плоскости $\kappa $ дополнительный разрез по лучу $( - \infty , - 1]$. Функции Бесселя в подынтегральных выражениях в (2) и (3) представим в виде полусуммы функций Ханкеля I и II рода
${{J}_{m}}(\kappa {{r}_{a}}) = \frac{1}{2}\{ H_{m}^{{(1)}}(\kappa {{r}_{a}}) + H_{m}^{{(2)}}(\kappa {{r}_{a}})\} ,\quad m = 0,\;1,$
под которыми понимаем такие однозначные аналитические в области $ - \pi < arg\kappa < \pi $ функции, что при $\kappa = x > 0$ их полусумма действительна. Итак, подынтегральные функции комплексной переменной $\kappa $ мы будем рассматривать на всей комплексной плоскости с разрезом вдоль луча $( - \infty ,1]$. Это подобласть области определения функции $y(\kappa )$. Заметим, что значения функции $H_{m}^{{(n)}}(\kappa {{r}_{a}})$, $m = 0,\;1$, $n = 1,\;2$, при $\kappa = x \pm i0$, $0 < x < 1$ совпадают в силу ее аналитичности в области $ - \pi < arg\kappa < \pi $.

Обозначим через ${{C}_{R}}$ и $C_{R}^{'}$ дуги окружностей $\kappa = R{{{\text{e}}}^{{is}}}$, $s \in [0,\pi {\text{/}}2]$ и $\kappa = R{{{\text{e}}}^{{is}}}$, $s \in [ - \pi {\text{/}}2,0]$, через ${{l}_{R}}$, $l_{R}^{'}$, ${{\gamma }_{R}}$ отрезки $[i{{\rho }_{3}},iR]$, $[ - iR, - i{{\rho }_{4}}]$, $[{{\kappa }_{R}} + {{\rho }_{2}},R]$ соответственно через $\gamma _{1}^{{(u)}}$, $\gamma _{2}^{{(u)}}$ отрезки $[{{\rho }_{3}},\;1 - {{\rho }_{1}}]$, $[1 - {{\rho }_{1}},\;1 - {{\rho }_{5}}]$ (${{\rho }_{5}} < {{\rho }_{1}}$) на верхнем берегу разреза $( - \infty ,1]$ соответственно через $\gamma _{4}^{{(l)}}$ отрезок $[{{\rho }_{4}},\;1 - {{\rho }_{5}}]$ на нижнем берегу разреза $( - \infty ,1]$, проходимый справа налево, через ${{C}_{3}}$ и ${{C}_{4}}$ дуги окружностей $\kappa = {{\rho }_{3}}{{{\text{e}}}^{{is}}}$, $s \in [0,\pi {\text{/}}2]$ и $\kappa = {{\rho }_{4}}{{{\text{e}}}^{{is}}}$, $s \in [ - \pi {\text{/}}2,0]$ и через ${{C}_{5}}$ окружность $\kappa = 1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}}$, $s \in [ - \pi ,\pi ]$.

Введем также обозначения

$h_{0}^{{(n)}}(\kappa ) = \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\frac{{H_{0}^{{(n)}}(\kappa {{r}_{a}})\kappa \sqrt {{{\kappa }^{2}} - 1} }}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}},$
$h_{1}^{{(n)}}(\kappa ) = - \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\frac{{H_{1}^{{(n)}}(\kappa {{r}_{a}}){{\kappa }^{2}}}}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}},$
${{h}^{{(n)}}}(\kappa ) = \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\frac{{H_{0}^{{(n)}}(\kappa {{r}_{a}})\kappa }}{{\beta \sqrt {{{\kappa }^{2}} - 1} - 1}}{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{\kappa }^{2}} - 1} - {{l}_{a}}\kappa }}},\quad n = 1,\;2.$
Интегралы от $h_{m}^{{(n)}}$ и ${{h}^{{(n)}}}(\kappa )$ вдоль $L$ будем понимать как предел при $R \to + \infty $, ${{\rho }_{3}} \to + 0$ интегралов
$\int\limits_{{{L}_{1}}} \,h_{m}^{{(n)}}(\kappa )d\kappa ,\quad \int\limits_{{{L}_{1}}} \,{{h}^{{(n)}}}(\kappa )d\kappa ,\quad {{L}_{1}} = C_{3}^{ - } \cup \gamma _{1}^{{(u)}} \cup C_{1}^{ - } \cup {{\gamma }_{3}} \cup C_{2}^{ - } \cup {{\gamma }_{R}}.$
Пусть ${{\Gamma }_{1}}$ – замкнутая кусочно гладкая кривая вида ${{\Gamma }_{1}} = {{L}_{1}} \cup {{C}_{R}} \cup l_{R}^{ - }$, проходимая против часовой стрелки. По теореме Коши имеем
(4)
$\int\limits_{{{\Gamma }_{1}}} \,h_{m}^{{(1)}}(\kappa )d\kappa = 0,\quad \int\limits_{{{\Gamma }_{1}}} \,{{h}^{{(1)}}}(\kappa )d\kappa = 0,\quad m = 0,1.$
Покажем, что
(5)
$\mathop {lim}\limits_{R \to + \infty } \int\limits_{{{C}_{R}}} \,h_{m}^{{(1)}}(\kappa )d\kappa = 0,\quad m = 0,1.$
Пусть $\kappa = R{{{\text{e}}}^{{is}}}$, $s \in [0,\pi {\text{/}}2]$. Тогда
$\left| {{{I}_{{Rm}}}} \right|: = \left| {\int\limits_{{{C}_{R}}} h_{m}^{{(1)}}(\kappa )d\kappa } \right| \leqslant \int\limits_0^{\pi /2} \left| {h_{m}^{{(1)}}(R{{{\text{e}}}^{{is}}})} \right|Rds.$
Для достаточно большого $R$ и любого $s \in [0,\pi {\text{/}}2]$ справедливы неравенства
$\left| {\sqrt {{{R}^{2}}{{{\text{e}}}^{{2is}}} - 1} } \right| < 2R,\quad \left| {{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{R}^{2}}{{{\text{e}}}^{{2is}}} - 1} }}}} \right| \leqslant 1.$
Мы учли, что $\arg (\sqrt {{{R}^{2}}{{{\text{e}}}^{{2is}}} - 1} )$ изменяется на сегменте $[0,\pi ]$. Поэтому
(6)
$\left| {{{I}_{{Rm}}}} \right| < \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}}}\frac{{{{R}^{3}}}}{{\beta \sqrt {{{R}^{2}} - 1} - 1}}\int\limits_0^{\pi /2} \left| {H_{m}^{{(1)}}({{r}_{a}}R{{{\text{e}}}^{{is}}})} \right|{{{\text{e}}}^{{ - {{l}_{a}}Rcoss}}}ds.$
Считая, что $R \gg 1$, воспользуемся асимптотикой функции Ханкеля [2]. Вместо (6) будем иметь
(7)
$\left| {{{I}_{{Rm}}}} \right| < \frac{{{{R}^{{5/2}}}}}{{\beta \sqrt {{{R}^{2}} - 1} - 1}}\left( {{{A}_{1}} + \frac{{{{A}_{2}}}}{R}} \right)\int\limits_0^{\pi /2} {{{\text{e}}}^{{ - {{r}_{a}}Rsins - {{l}_{a}}Rcoss}}}ds < \frac{{A{{R}^{{5/2}}}}}{{\beta \sqrt {{{R}^{2}} - 1} - 1}}\int\limits_0^{\pi /2} {{{\text{e}}}^{{ - {{r}_{a}}Rsins - {{l}_{a}}Rcoss}}}ds,\quad {{A}_{1}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{\rho {{a}^{2}}\sqrt {\pi {{r}_{a}}} }},$
${{A}_{2}}$ – некоторое положительное число, не зависящее от $R$, $A = {{A}_{1}} + {{A}_{2}}$. Интеграл в правой части (7) оценивается также, как при доказательстве леммы Жордана [5]
$\int\limits_0^{\pi /2} {{{\text{e}}}^{{ - {{r}_{a}}Rsins - {{l}_{a}}Rcoss}}}ds \leqslant {{{\text{e}}}^{{ - R{{l}_{a}}}}}\int\limits_0^{\pi /2} {{{\text{e}}}^{{2R({{l}_{a}} - {{r}_{a}})s/\pi }}}ds = \frac{\pi }{{2R}}({{{\text{e}}}^{{ - R{{r}_{a}}}}} - {{{\text{e}}}^{{ - R{{l}_{a}}}}}).$
Отсюда
$0 \leqslant \left| {{{I}_{{Rm}}}} \right| < \frac{\pi }{2}\frac{{A{{R}^{{3/2}}}}}{{\beta \sqrt {{{R}^{2}} - 1} - 1}}({{{\text{e}}}^{{ - R{{r}_{a}}}}} - {{{\text{e}}}^{{ - R{{l}_{a}}}}}) \to + 0,\quad R \to + \infty ,$
т.е. ${{I}_{{Rm}}} \to + 0$ при $R \to + \infty $.

Покажем, что

$\mathop {lim}\limits_{{{\rho }_{3}} \to + 0} \int\limits_{C_{3}^{ - }} h_{m}^{{(1)}}(\kappa )d\kappa = 0,\quad m = 0,1.$
Пусть $\kappa = {{\rho }_{3}}{{{\text{e}}}^{{is}}}$, $s \in [0,\pi {\text{/}}2]$. Тогда
$\left| {{{I}_{{3m}}}} \right|: = \left| {\int\limits_{C_{3}^{ - }} \,h_{m}^{{(1)}}(\kappa )d\kappa } \right| \leqslant \int\limits_0^{\pi /2} \left| {h_{m}^{{(1)}}({{\rho }_{3}}{{{\text{e}}}^{{is}}})} \right|{{\rho }_{3}}ds.$
Для ${{\rho }_{3}} \in (0,1)$ и любого $s \in [0,\pi {\text{/}}2]$ справедливы неравенства
$\left| {{{{\text{e}}}^{{{{z}_{a}}\sqrt {\rho _{3}^{2}{{{\text{e}}}^{{2is}}} - 1} }}}} \right| \leqslant 1,\quad \left| {{{{\text{e}}}^{{ - {{l}_{a}}{{\rho }_{3}}{{{\text{e}}}^{{is}}}}}}} \right| \leqslant 1,\quad \left| {\sqrt {\rho _{3}^{2}{{{\text{e}}}^{{2is}}} - 1} } \right| \leqslant \sqrt {1 + \rho _{3}^{2}} < 2.$
Оценим $\beta \sqrt {\rho _{3}^{2}{{{\text{e}}}^{{2is}}} - 1} - 1$ при достаточно малом ${{\rho }_{3}}$ и любом $s \in [0,\pi {\text{/}}2]$. Для абсолютной величины будем иметь
$\left| {\beta \sqrt {\rho _{3}^{2}{{{\text{e}}}^{{2is}}} - 1} - 1} \right| = \left| {\beta d(s){{{\text{e}}}^{{i(\pi /2 - \psi (s))}}} - 1} \right| = {{\{ {{\beta }^{2}}{{d}^{2}}(s) - 2\beta d(s)sin\psi (s) + 1\} }^{{1/2}}},$
где
$d(s) = {{(\rho _{3}^{4} - 2\rho _{3}^{2}cos2s + 1)}^{{1/4}}},\quad \psi (s) = \frac{1}{2}arctg\left( {\frac{{\rho _{3}^{2}sin2s}}{{1 - \rho _{3}^{2}cos2s}}} \right).$
Следовательно,
(8)
$\begin{gathered} \left| {\beta \sqrt {\rho _{3}^{2}{{{\text{e}}}^{{2is}}} - 1} - 1} \right| > {{\{ {{\beta }^{2}}{{d}^{2}}(s) - \beta d(s) + 1\} }^{{1/2}}} \geqslant {{\{ {{\beta }^{2}}{{(\rho _{3}^{4} - 2\rho _{3}^{2} + 1)}^{{1/2}}} - \beta (\rho _{3}^{4} + 2\rho _{2}^{2} + 1) + 1\} }^{{1/2}}} = \\ = \{ {{\beta }^{2}}{{(1 - \rho _{3}^{2})}^{2}} - \beta {{(\rho _{3}^{2} + 1)}^{2}} + {{1\} }^{{1/2}}}. \\ \end{gathered} $
Отметим, что правая часть (8) стремится к ${{({{\beta }^{2}} - \beta + 1)}^{{1/2}}} > 0$ при ${{\rho }_{3}} \to + 0$.

Представим функции Ханкеля в виде

$H_{m}^{{(1)}}(\kappa ) = {{J}_{m}}(\kappa ) + i{{N}_{m}}(\kappa ),\quad m = 0,1,$
где ${{N}_{m}}(\kappa )$ – функция Неймана индекса $m$. Известно [2], что при $arg\kappa \in ( - \pi ,\pi )$
${{N}_{0}}(\kappa ) = \frac{2}{\pi }{{J}_{0}}(\kappa )\left( {ln\frac{\kappa }{2} + C} \right) + f(\kappa ),\quad {{N}_{1}}(\kappa ) = \frac{2}{\pi }{{J}_{1}}(\kappa )\left( {ln\frac{\kappa }{2} + C} \right) - \frac{2}{{\pi \kappa }} + g(\kappa ),$
где $f(\kappa )$, $g(\kappa )$ – целые функции, $ln\kappa = ln\left| \kappa \right| + iarg\kappa $, а $C$ – постоянная Эйлера. Тогда для достаточно малого ${{\rho }_{3}}$ и всех $s \in [0,\pi {\text{/}}2]$
$\left| {H_{0}^{{(1)}}({{r}_{a}}{{\rho }_{3}}{{{\text{e}}}^{{is}}})} \right| \leqslant {{C}_{0}}({{r}_{a}})\left\{ {1 + \frac{2}{\pi }\left( {\left| {ln{{\rho }_{3}}} \right| + \left| {ln\frac{{{{r}_{a}}}}{2}} \right| + \frac{\pi }{2} + C} \right)} \right\} + {{C}_{2}}({{r}_{a}}): = {{A}_{0}}({{r}_{a}},{{\rho }_{3}}),$
$\left| {H_{1}^{{(1)}}({{r}_{a}}{{\rho }_{3}}{{{\text{e}}}^{{is}}})} \right| \leqslant {{C}_{1}}({{r}_{a}})\left\{ {1 + \frac{2}{\pi }\left( {\left| {ln{{\rho }_{3}}} \right| + \left| {ln\frac{{{{r}_{a}}}}{2}} \right| + \frac{\pi }{2} + C} \right)} \right\} + \frac{2}{{\pi {{\rho }_{3}}}} + {{C}_{3}}({{r}_{a}}): = {{A}_{1}}({{r}_{a}},{{\rho }_{3}}),$
где ${{C}_{1}}({{r}_{a}})$, ${{C}_{2}}({{r}_{a}})$ не зависят от ${{\rho }_{3}}$. Мы воспользовались тем, что функции ${{J}_{m}}(\kappa {{r}_{a}})$, $f(\kappa {{r}_{a}})$, $g(\kappa {{r}_{a}})$ ограничены в окрестности точки $\kappa = 0$ при любом фиксированном ${{r}_{a}}$. Следовательно,
$0 \leqslant \left| {{{I}_{{3m}}}} \right| < \frac{{{{p}_{0}}l_{a}^{2}\rho _{3}^{{2 + m}}{{A}_{m}}({{r}_{a}},{{\rho }_{3}})}}{{\rho {{a}^{2}}{{{\left\{ {{{\beta }^{2}}{{{(1 - \rho _{3}^{2})}}^{2}} - \beta {{{(1 + \rho _{3}^{2})}}^{2}} + 1} \right\}}}^{{1/2}}}}} \to + 0,\quad {{\rho }_{3}} \to + 0,$
т.е. ${{I}_{{3m}}} \to + 0$ при ${{\rho }_{3}} \to + 0$. Аналогичным образом доказывается, что интегралы по ${{C}_{R}}$ $C_{3}^{ - }$ от ${{h}^{{(1)}}}(\kappa )$ также стремятся к нулю. Перейдем в (4) к пределу при $R \to + \infty $, ${{\rho }_{3}} \to + 0$. Тогда
(9)
$\int\limits_L \,h_{0}^{{(1)}}(\kappa )d\kappa = - I_{z}^{{(1)}},\quad \int\limits_L \,h_{1}^{{(1)}}(\kappa )d\kappa = - I_{r}^{{(1)}},\quad \int\limits_L \,{{h}^{{(1)}}}(\kappa )d\kappa = - {{I}^{{(1)}}},$
где $I_{z}^{{(1)}}$, $I_{r}^{{(1)}}$, ${{I}^{{(1)}}}$ есть интегралы по неотрицательной части мнимой оси. Они имеют вид
$I_{z}^{{(1)}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{0}}(\sigma {{r}_{a}})\sigma \sqrt {{{\sigma }^{2}} + 1} }}{{i\beta \sqrt {{{\sigma }^{2}} + 1} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} - i{{l}_{a}}\sigma }}}d\sigma ,$
$I_{r}^{{(1)}} = i\frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{1}}(\sigma {{r}_{a}}){{\sigma }^{2}}}}{{i\beta \sqrt {{{\sigma }^{2}} + 1} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} - i{{l}_{a}}\sigma }}}d\sigma ,$
${{I}^{{(1)}}} = - i\frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{0}}(\sigma {{r}_{a}})\sigma }}{{i\beta \sqrt {{{\sigma }^{2}} + 1} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} - i{{l}_{a}}\sigma }}}d\sigma ,$
${{K}_{0}}({{r}_{a}}\sigma )$, ${{K}_{1}}({{r}_{a}}\sigma )$ – функции Макдональда нулевого и первого индексов. Мы воспользовались соотношениями [2]

$H_{0}^{{(1)}}( - i\sigma ) = - i\frac{2}{\pi }{{K}_{0}}(\sigma ),\quad H_{1}^{{(1)}}( - i\sigma ) = - \frac{2}{\pi }{{K}_{1}}(\sigma ).$

Интегралы от $h_{m}^{{(2)}}(\kappa )$ вдоль ${{L}_{1}}$ представим в виде суммы интегралов вдоль ${{L}_{2}} = C_{3}^{ - } \cup \gamma _{1}^{{(u)}}$ и ${{L}_{3}} = C_{1}^{ - } \cup {{\gamma }_{3}} \cup C_{2}^{ - } \cup {{\gamma }_{R}}$. Пусть ${{\Gamma }_{2}}$ – замкнутая кусочно гладкая кривая вида ${{\Gamma }_{2}} = \gamma _{2}^{{(u)}} \cup C_{5}^{ - } \cup \gamma _{4}^{{(l)}} \cup $ $ \cup \;C_{4}^{ - } \cup l_{R}^{'} \cup C_{R}^{'} \cup L_{3}^{ - }$, проходимая против часовой стрелки. По основной теореме теории вычетов

(10)
$\int\limits_{{{\Gamma }_{2}}} \,h_{m}^{{(2)}}(\kappa )d\kappa = 2\pi iRes[h_{m}^{{(2)}}(\kappa ),{{\kappa }_{R}}],\quad \int\limits_{{{\Gamma }_{2}}} \,{{h}^{{(2)}}}(\kappa )d\kappa = 2\pi iRes[{{h}^{{(2)}}}(\kappa ),{{\kappa }_{R}}].$

Подобно тому, как это было сделано выше, можно показать, что

$\mathop {lim}\limits_{R \to + \infty } \int\limits_{C_{R}^{'}} h_{m}^{{(2)}}(\kappa )d\kappa = \mathop {lim}\limits_{R \to + \infty } \int\limits_{C_{R}^{'}} {{h}^{{(2)}}}(\kappa )d\kappa = 0,$
$\mathop {lim}\limits_{{{\rho }_{4}} \to + 0} \int\limits_{C_{4}^{ - }} h_{m}^{{(2)}}(\kappa )d\kappa = \mathop {lim}\limits_{{{\rho }_{4}} \to + 0} \int\limits_{C_{4}^{ - }} {{h}^{{(2)}}}(\kappa )d\kappa = 0,\quad m = 0,1.$
Покажем, что
$\mathop {lim}\limits_{{{\rho }_{5}} \to + 0} \int\limits_{C_{5}^{ - }} h_{m}^{{(2)}}(\kappa )d\kappa = 0,\quad m = 0,1.$
Пусть $\kappa = 1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}}$, $s \in [ - \pi ,\pi ]$. Тогда имеем
$\left| {{{I}_{{5m}}}} \right|: = \left| {\int\limits_{{{C}_{5}}} h_{m}^{{(2)}}(\kappa )d\kappa } \right| \leqslant \int\limits_{ - \pi }^\pi \left| {h_{m}^{{(2)}}(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})} \right|{{\rho }_{5}}ds.$
Для достаточно малого ${{\rho }_{5}}$ и любого $s \in [ - \pi ,\pi ]$ справедливы неравенства

$\left| {{{{(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})}}^{2}}} \right| < {{(1 + {{\rho }_{5}})}^{2}} < 4,\quad \left| {\sqrt {(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}}^{2}) - 1} } \right| \leqslant \sqrt {2{{\rho }_{5}} + \rho _{5}^{2}} < 2,$
$\left| {\beta \sqrt {{{{(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})}}^{2}} - 1} - 1} \right| \geqslant 1 - \beta \sqrt {2{{\rho }_{5}} + \rho _{5}^{2}} ,\quad \left| {(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})\sqrt {{{{(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})}}^{2}} - 1} } \right| < 2(1 + {{\rho }_{5}}) < 4,$
$\left| {{{{\text{e}}}^{{ - {{l}_{a}}(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})}}}} \right| \leqslant {{{\text{e}}}^{{ - {{l}_{a}}(1 - {{\rho }_{5}})}}} < 1,\quad \left| {{{{\text{e}}}^{{{{z}_{a}}\sqrt {{{{(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}})}}^{2}} - 1} }}}} \right| \leqslant {{{\text{e}}}^{{ - {{z}_{a}}\sqrt {2{{\rho }_{5}} + \rho _{5}^{2}} }}} < {{{\text{e}}}^{{ - 2{{z}_{a}}}}}.$

Функции $H_{m}^{{(2)}}({{r}_{a}}\kappa )$ аналитичны в окрестности точки $\kappa = 1$, а следовательно, ограничены там, т.е. для любого $s \in [ - \pi ,\pi ]$, любого достаточно малого ${{\rho }_{5}}$ и фиксированного ${{r}_{a}}$

$\left| {H_{m}^{{(2)}}({{r}_{a}}(1 + {{\rho }_{5}}{{{\text{e}}}^{{is}}}))} \right| \leqslant {{C}_{{2m}}}({{r}_{a}}).$
Поэтому
$0 \leqslant \left| {{{I}_{{5m}}}} \right| < C{\text{'}}({{r}_{a}})\frac{{4{{p}_{0}}l_{a}^{2}\pi }}{{\rho {{a}^{2}}}}\frac{{{{\rho }_{5}}{{{\text{e}}}^{{ - 2{{z}_{a}}}}}}}{{1 - \beta \sqrt {2{{\rho }_{5}} + \rho _{5}^{2}} }} \to + 0,\quad {{\rho }_{5}} \to + 0,$
где $C{\text{'}}({{r}_{a}}) = max\left\{ {{{C}_{{20}}}({{r}_{a}}),{{C}_{{21}}}({{r}_{a}})} \right\}$, т.е. ${{I}_{{5m}}} \to + 0$ при ${{\rho }_{5}} \to + 0$. Точно также доказывается, что интеграл по ${{C}_{5}}$ от ${{h}^{{(2)}}}(\kappa )$ стремится к нулю.

Перейдем в равенствах (10) к пределу при $R \to + \infty $, ${{\rho }_{4}} \to + 0$, ${{\rho }_{5}} \to + 0$. Учитывая, что

$\mathop {lim}\limits_{{{\rho }_{3}} \to + 0} \int\limits_{{{L}_{2}}} ... + \mathop {lim}\limits_{R \to + \infty } \int\limits_{{{L}_{3}}} ... = \int\limits_L ...,$
окончательно получаем
$\int\limits_L \,h_{0}^{{(2)}}(\kappa )d\kappa = I_{z}^{{(2)}} + {{I}_{{z(u)}}} - {{I}_{{z(l)}}} - 2\pi iRes[h_{0}^{{(2)}}(\kappa ),{{\kappa }_{R}}],$
$\int\limits_L \,h_{1}^{{(2)}}(\kappa )d\kappa = I_{r}^{{(2)}} + {{I}_{{r(u)}}} - {{I}_{{r(l)}}} - 2\pi iRes[h_{1}^{{(2)}}(\kappa ),{{\kappa }_{R}}],$
$\int\limits_L \,{{h}^{{(2)}}}(\kappa )d\kappa = {{I}^{{(2)}}} + {{I}_{{(u)}}} + {{I}_{{(l)}}} - 2\pi iRes[{{h}^{{(2)}}}(\kappa ),{{\kappa }_{R}}],$
где $I_{z}^{{(2)}}$, $I_{r}^{{(2)}}$, ${{I}^{{(2)}}}$ – интегралы по неположительной мнимой оси, ${{I}_{{z(u)}}}$, ${{I}_{{z(l)}}}$, ${{I}_{{r(u)}}}$, ${{I}_{{r(l)}}}$, ${{I}_{{(u)}}}$, ${{I}_{{(l)}}}$ – интегралы по отрезкам на берегах разреза. Указанные интегралы имеют вид
$I_{z}^{{(2)}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{0}}(\sigma {{r}_{a}})\sigma \sqrt {{{\sigma }^{2}} + 1} }}{{i\beta \sqrt {{{\sigma }^{2}} + 1} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} + i{{l}_{a}}\sigma }}}d\sigma ,$
$I_{r}^{{(2)}} = i\frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{1}}(\sigma {{r}_{a}}){{\sigma }^{2}}}}{{i\beta \sqrt {{{\sigma }^{2}} + 1} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} + i{{l}_{a}}\sigma }}}d\sigma ,$
${{I}^{{(2)}}} = - i\frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}\int\limits_0^{ + \infty } \frac{{{{K}_{0}}(\sigma {{r}_{a}})\sigma }}{{i\beta \sqrt {{{\sigma }^{2}} + 1} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {{{\sigma }^{2}} + 1} + i{{l}_{a}}\sigma }}}d\sigma ,$
${{I}_{{z(u)}}} = i\frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{0}^{{(2)}}(x{{r}_{a}})x\sqrt {1 - {{x}^{2}}} }}{{i\beta \sqrt {1 - {{x}^{2}}} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{I}_{{z(l)}}} = i\frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{0}^{{(2)}}(x{{r}_{a}})x\sqrt {1 - {{x}^{2}}} }}{{i\beta \sqrt {1 - {{x}^{2}}} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{I}_{{r(u)}}} = - \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{1}^{{(2)}}(x{{r}_{a}}){{x}^{2}}}}{{i\beta \sqrt {1 - {{x}^{2}}} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{I}_{{r(l)}}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{1}^{{(2)}}(x{{r}_{a}}){{x}^{2}}}}{{i\beta \sqrt {1 - {{x}^{2}}} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{I}_{{(u)}}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{0}^{{(2)}}(x{{r}_{a}})x}}{{i\beta \sqrt {1 - {{x}^{2}}} - 1}}{{{\text{e}}}^{{i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{I}_{{(l)}}} = \frac{{{{p}_{0}}l_{a}^{2}}}{{2\rho {{a}^{2}}}}\int\limits_0^1 \frac{{H_{0}^{{(2)}}(x{{r}_{a}})x}}{{i\beta \sqrt {1 - {{x}^{2}}} + 1}}{{{\text{e}}}^{{ - i{{z}_{a}}\sqrt {1 - {{x}^{2}}} - {{l}_{a}}x}}}dx,$
${{K}_{0}}({{r}_{a}}\sigma )$, ${{K}_{1}}({{r}_{a}}\sigma )$ – функции Макдональда нулевого и первого индексов. Мы воспользовались соотношениями [2] в виде
$H_{0}^{{(2)}}( - i\sigma ) = i\frac{2}{\pi }{{K}_{0}}(\sigma ),\quad H_{1}^{{(2)}}( - i\sigma ) = - i\frac{2}{\pi }{{K}_{1}}(\sigma ).$
Отметим, что интегралы по отрезкам $[0,\;1 - {{\rho }_{1}}]$ и $[1 - {{\rho }_{1}},\;1]$ на верхнем берегу разреза мы объединили в один по отрезку $[0,1]$. Для комплексных амплитуд получаем формулы

(11)
${{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}}) = - {{I}^{{(1)}}} + {{I}^{{(2)}}} + {{I}_{{(u)}}} - {{I}_{{(l)}}} - 2\pi iRes[{{h}^{{(2)}}}(\kappa ),{{\kappa }_{R}}],$
(12)
$\begin{gathered} U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - I_{z}^{{(1)}} + I_{z}^{{(2)}} + {{I}_{{z(u)}}} - {{I}_{{z(l)}}} - 2\pi iRes[h_{0}^{{(2)}}(\kappa ),{{\kappa }_{R}}], \\ U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - I_{r}^{{(1)}} + I_{r}^{{(2)}} + {{I}_{{r(u)}}} - {{I}_{{r(l)}}} - 2\pi iRes[h_{1}^{{(2)}}(\kappa ),{{\kappa }_{R}}]. \\ \end{gathered} $

Покажем, что при ${{r}_{a}} \gg 1$ интегралы в правых частях (11), (12) асимптотически малы по сравнению с вычетами в полюсе Рэлея ${{\kappa }_{R}}$.

Оценим сначала интегралы $I_{z}^{{(m)}}$, $I_{r}^{{(m)}}$, ${{I}^{{(m)}}}$ ($m = 1,2$). При любых ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$ имеем

$\begin{gathered} \left| {I_{z}^{{(m)}}} \right| \leqslant B\int\limits_0^{ + \infty } {{K}_{0}}({{r}_{a}}\sigma )\sigma \sqrt {{{\sigma }^{2}} + 1} d\sigma \leqslant \sqrt 2 B\int\limits_0^1 {{K}_{0}}({{r}_{a}}\sigma )\sigma d\sigma + \sqrt 2 B\int\limits_1^{ + \infty } {{K}_{0}}({{r}_{a}}\sigma ){{\sigma }^{2}}d\sigma < \\ < \;\sqrt 2 B\int\limits_0^{ + \infty } {{K}_{0}}({{r}_{a}}\sigma )\sigma d\sigma + \sqrt 2 B\int\limits_0^{ + \infty } {{K}_{0}}({{r}_{a}}\sigma ){{\sigma }^{2}}d\sigma = \frac{{\sqrt 2 B}}{{r_{a}^{2}}} + \frac{{\pi B}}{{\sqrt 2 r_{a}^{3}}} < \frac{{{{C}_{1}}}}{{r_{a}^{2}}}, \\ \end{gathered} $
$\left| {I_{r}^{{(m)}}} \right| \leqslant B\int\limits_0^{ + \infty } {{K}_{1}}({{r}_{a}}\sigma ){{\sigma }^{2}}d\sigma = \frac{{2B}}{{r_{a}^{3}}},\quad \left| {{{I}^{{(m)}}}} \right| \leqslant B\int\limits_0^{ + \infty } {{K}_{0}}({{r}_{a}}\sigma )\sigma d\sigma = \frac{B}{{r_{a}^{2}}},\quad B = \frac{{{{p}_{0}}l_{a}^{2}}}{{\pi \rho {{a}^{2}}}}.$
Мы воспользовались формулой [2] в виде
$\int\limits_0^{ + \infty } {{K}_{\nu }}({{r}_{a}}\sigma ){{\sigma }^{\mu }}d\sigma {{ = 2}^{{\mu - 1}}}r_{a}^{{ - \mu - 1}}\Gamma \left( {\frac{{\mu + 1 + \nu }}{2}} \right)\Gamma \left( {\frac{{\mu + 1 - \nu }}{2}} \right)$
при $\nu = 0$, $\mu = 1,2$ и при $\nu = 1$, $\mu = 2$. Отсюда получаем, что
$I_{z}^{{(m)}} = O\left( {\frac{1}{{r_{a}^{2}}}} \right),\quad I_{r}^{{(m)}} = O\left( {\frac{1}{{r_{a}^{3}}}} \right),\quad {{I}^{{(m)}}} = O\left( {\frac{1}{{r_{a}^{2}}}} \right)$
при любых ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$, а следовательно, при ${{r}_{a}} \gg 1$ и любом ${{z}_{a}} \leqslant 0$. Оценим интегралы по берегам разреза. Рассмотрим, например, интеграл ${{I}_{{z(u)}}}$. Представим его в виде
${{I}_{{z(u)}}} = iB\frac{\pi }{2}\left\{ {\int\limits_0^{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}} ... + \int\limits_{1/\sqrt {{{r}_{a}}} }^1 ...} \right\}: = iB\frac{\pi }{2}\left( {{{I}_{1}} + {{I}_{2}}} \right),\quad {{r}_{a}} \gg 1.$
Тогда
$\left| {{{I}_{{z(u)}}}} \right| \leqslant B\frac{\pi }{2}\left( {\left| {{{I}_{1}}} \right| + \left| {{{I}_{2}}} \right|} \right),\quad \left| {{{I}_{1}}} \right| \leqslant \int\limits_0^{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}} \left| {H_{0}^{{(2)}}(x{{r}_{a}})} \right|xdx = \frac{1}{{r_{a}^{2}}}\int\limits_0^{\sqrt {{{r}_{a}}} } \left| {H_{0}^{{(2)}}(s)} \right|sds < \frac{{{{C}_{2}}}}{{r_{a}^{{5/4}}}}.$
Последнее неравенство справа имеет место, поскольку для непрерывной на сегменте $[0,\sqrt {{{r}_{a}}} ]$ функции $\left| {H_{0}^{{(2)}}(s)} \right|s$

$\mathop {lim}\limits_{{{r}_{a}} \to + \infty } \frac{{\int\limits_0^{\sqrt {{{r}_{a}}} } \left| {H_{0}^{{(2)}}(s)} \right|sds}}{{r_{a}^{{3/4}}}} = \mathop {lim}\limits_{{{r}_{a}} \to + \infty } \frac{{2\left| {H_{0}^{{(2)}}\left( {\sqrt {{{r}_{a}}} } \right)} \right|}}{{3r_{a}^{{ - 1/4}}}} = \frac{{{{2}^{{3/2}}}}}{{3{{\pi }^{{1/2}}}}}.$

На сегменте $[1{\text{/}}\sqrt {{{r}_{a}}} ,\;1]$ при ${{r}_{a}} \gg 1$ значение аргумента $x{{r}_{a}}$ функции Ханкеля $H_{0}^{{(2)}}(x{{r}_{a}})$ также много больше единицы. Поэтому для всех $x$ из указанного сегмента функцию Ханкеля можно представить асимптотической формулой [2]. Для ${{I}_{2}}$ будем иметь

$\begin{gathered} {{I}_{2}} = \frac{{\sqrt 2 }}{{\sqrt {\pi {{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 \frac{{\sqrt x \sqrt {1 - {{x}^{2}}} }}{{i\beta \sqrt {1 - {{x}^{2}}} - 1}}{{{\text{e}}}^{{i({{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \pi /4)}}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}dx + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right) = \\ = - {\kern 1pt} \frac{1}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{1}}(x)cos\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx + \\ \end{gathered} $
(13)
$ + \frac{1}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{2}}(x)sin\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx - $
$\begin{gathered} - \frac{i}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{2}}(x)cos\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx - \\ - \;\frac{i}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{1}}(x)sin\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right), \\ \end{gathered} $
где
${{f}_{1}}(x) = \frac{{\sqrt 2 }}{{\sqrt \pi }}\frac{{\sqrt x \sqrt {1 - {{x}^{2}}} }}{{1 + {{\beta }^{2}}(1 - {{x}^{2}})}}{{{\text{e}}}^{{ - {{l}_{a}}x}}},\quad {{f}_{2}}(x) = \beta \frac{{\sqrt 2 }}{{\sqrt \pi }}\frac{{\sqrt x (1 - {{x}^{2}})}}{{1 + {{\beta }^{2}}(1 - {{x}^{2}})}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}.$
Интегралы в (13) оцениваются однотипно. Рассмотрим, например, первый, обозначив его через ${{I}_{{21}}}$. Имеем
(14)
${{I}_{{21}}} = - \frac{1}{{\sqrt {{{r}_{a}}} }}\left\{ {\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{g}_{1}}(x)cos\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)dx + \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{g}_{2}}(x)sin\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)dx} \right\},$
где
${{g}_{1}}(x) = {{f}_{1}}(x)cos({{z}_{a}}\sqrt {1 - {{x}^{2}}} ),\quad {{g}_{2}}(x) = {{f}_{1}}(x)sin({{z}_{a}}\sqrt {1 - {{x}^{2}}} ).$
Отметим, что функции ${{g}_{1}}(x)$, ${{g}_{2}}(x)$ непрерывны на $[0,1]$ и непрерывно дифференцируемы соответственно на $(0,1)$ и $(0,1]$. Кроме того, производные $g_{k}^{'}(x)$, $k = 1,2$ абсолютно интегрируемы на $[0,1]$, поскольку при $x \in (0,1)$ имеют место оценки
(15)
$\left| {g_{k}^{'}(x)} \right| < \frac{{\sqrt 2 }}{{\sqrt \pi }}\left( {\frac{1}{{2\sqrt x }} + \frac{1}{{\sqrt {1 - {{x}^{2}}} }} + {{l}_{a}} + \left| {{{z}_{a}}} \right| + 2{{\beta }^{2}}} \right).$
Поэтому для ${{I}_{{21}}}$ справедлива формула интегрирования по частям. Тогда вместо (14) будем иметь
$\begin{gathered} {{I}_{{21}}} = \frac{1}{{r_{a}^{{3/2}}}}\left\{ { - \left. {{{g}_{1}}(x)sin\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)} \right|_{{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}}^{1} + \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 g_{1}^{'}(x)sin\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)dx} \right\} + \\ + \;\frac{1}{{r_{a}^{{3/2}}}}\left\{ {{{g}_{2}}(x)\left. {cos\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)} \right|_{{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}}^{1} - \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 g_{2}^{'}(x)cos\left( {x{{r}_{a}} - \frac{\pi }{4}} \right)dx} \right\}. \\ \end{gathered} $
Так как функции ${{g}_{k}}(x)$, $k = 1,2,$ непрерывны на сегменте $[0,1]$, то ограничены на нем, а следовательно, равномерно ограничены на сегменте $[1{\text{/}}\sqrt {{{r}_{a}}} ,\;1]$, т.е. $\left| {{{g}_{k}}(x)} \right| \leqslant {{C}_{3}}$ для любого $x \in [1{\text{/}}\sqrt {{{r}_{a}}} ,\;1]$ и любого ${{r}_{a}} > 1$. Отсюда с учетом (15) имеем
$\begin{gathered} \left| {{{I}_{{21}}}} \right| \leqslant \frac{1}{{r_{a}^{{3/2}}}}\left\{ {2{{C}_{3}} + \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 \left| {g_{1}^{'}(x)} \right|dx + \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 \left| {g_{2}^{'}(x)} \right|dx} \right\} < \frac{1}{{r_{a}^{{3/2}}}}\left\{ {2{{C}_{3}} + \int\limits_0^1 \left| {g_{1}^{'}(x)} \right|dx + \int\limits_0^1 \left| {g_{2}^{'}(x)} \right|dx} \right\} < \\ < \;\frac{2}{{r_{a}^{{3/2}}}}\left\{ {{{C}_{3}} + \frac{{\sqrt 2 }}{{\sqrt \pi }}\left( {1 + \frac{\pi }{2} + {{l}_{a}} + \left| {{{z}_{a}}} \right| + 2{{\beta }^{2}}} \right)} \right\} = \frac{{{{C}_{4}}({{z}_{a}})}}{{r_{a}^{{3/2}}}}. \\ \end{gathered} $
Мы воспользовались оценками (15).

Рассуждая аналогичным образом, получаем, что четвертый интеграл в (13) имеет такую же оценку, как и ${{I}_{{21}}}$, а второй и третий, обозначим их ${{I}_{{2k}}}$, $k = 2,3$, оценку

$\left| {{{I}_{{2k}}}} \right| < \frac{{{{C}_{6}}}}{{r_{a}^{{3/2}}}},\quad {{C}_{6}} = 2\left\{ {{{C}_{5}} + \beta \frac{{\sqrt 2 }}{{\sqrt \pi }}(3 + {{l}_{a}} + \left| {{{z}_{a}}} \right| + 2{{\beta }^{2}})} \right\},$
где ${{C}_{5}}$ – постоянная, которая равномерно ограничивает функции
${{f}_{2}}(x)cos({{z}_{a}}\sqrt {1 - {{x}^{2}}} ),\quad {{f}_{2}}(x)sin({{z}_{a}}\sqrt {1 - {{x}^{2}}} ).$
Тем самым имеем
$\left| {{{I}_{2}}} \right| < \frac{{{{C}_{7}}({{z}_{a}})}}{{r_{a}^{{3/2}}}},\quad {{C}_{7}} = 2\left\{ {{{C}_{4}}({{z}_{a}}) + {{C}_{6}}({{z}_{a}})} \right\},$
$\left| {{{I}_{{z(u)}}}} \right| < \frac{{C{\text{'}}}}{{r_{a}^{{5/4}}}},\quad C{\text{'}} = B\frac{\pi }{2}\left\{ {{{C}_{2}} + {{C}_{7}}({{z}_{a}})} \right\},$
т.е. ${{I}_{{z(u)}}} = O(1{\text{/}}r_{a}^{{5/4}})$ при ${{r}_{a}} \gg 1$ и любом фиксированном ${{z}_{a}} \leqslant 0$. Точно также можно показать, что ${{I}_{{z(l)}}} = O(1{\text{/}}r_{a}^{{5/4}})$, ${{I}_{{(u)}}} = O(1{\text{/}}r_{a}^{{5/4}})$, ${{I}_{{(l)}}} = O(1{\text{/}}r_{a}^{{5/4}})$ при тех же ${{r}_{a}}$ и ${{z}_{a}}$.

Оценим интегралы ${{I}_{{r(u)}}}$, ${{I}_{{r(l)}}}$. При этом подробных пояснений делать не будем. Представим ${{I}_{{r(u)}}}$ в виде

${{I}_{{r(u)}}} = B\frac{\pi }{2}\left\{ {\int\limits_0^{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}} \ldots + \int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 \ldots } \right\}: = B\frac{\pi }{2}\left( {{{I}_{3}} + {{I}_{4}}} \right),\quad {{r}_{a}} \gg 1.$
Откуда имеем
$\left| {{{I}_{{r(u)}}}} \right| \leqslant B\frac{\pi }{2}\left( {\left| {{{I}_{3}}} \right| + \left| {{{I}_{4}}} \right|} \right),\quad \left| {{{I}_{3}}} \right| \leqslant \int\limits_0^{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}} \left| {H_{1}^{{(2)}}(x{{r}_{a}})} \right|{{x}^{2}}dx = \frac{1}{{r_{a}^{3}}}\int\limits_0^{\sqrt {{{r}_{a}}} } \left| {H_{1}^{{(2)}}(s)} \right|{{s}^{2}}ds < \frac{{C_{2}^{'}}}{{r_{a}^{{7/4}}}},$
так как верно
$\mathop {lim}\limits_{{{r}_{a}} \to + \infty } \frac{{\int\limits_0^{\sqrt {{{r}_{a}}} } \left| {H_{1}^{{(2)}}(s)} \right|{{s}^{2}}ds}}{{r_{a}^{{5/4}}}} = \mathop {lim}\limits_{{{r}_{a}} \to + \infty } \frac{{2\left| {H_{1}^{{(2)}}(\sqrt {{{r}_{a}}} )} \right|}}{{5r_{a}^{{ - 1/4}}}} = \frac{{{{2}^{{3/2}}}}}{{5{{\pi }^{{1/2}}}}}.$
Далее имеем
$\begin{gathered} {{I}_{4}} = i\frac{{\sqrt 2 }}{{\sqrt {\pi {{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 \frac{{{{x}^{{3/2}}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{i\beta \sqrt {1 - {{x}^{2}}} - 1}}{{{\text{e}}}^{{i({{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \pi /4)}}}dx + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right) = \\ = \; - {\kern 1pt} \frac{i}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{3}}(x)cos\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx - \\ \end{gathered} $
(16)
$ - \;\frac{i}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{4}}(x)sin\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx - $
$\begin{gathered} - \;\frac{1}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{4}}(x)cos\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx + \\ + \;\frac{1}{{\sqrt {{{r}_{a}}} }}\int\limits_{{1 \mathord{\left/ {\vphantom {1 {\sqrt {{{r}_{a}}} }}} \right. \kern-0em} {\sqrt {{{r}_{a}}} }}}^1 {{f}_{3}}(x)sin\left( {{{z}_{a}}\sqrt {1 - {{x}^{2}}} - x{{r}_{a}} + \frac{\pi }{4}} \right)dx + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right), \\ \end{gathered} $
где
${{f}_{3}}(x) = \frac{{\sqrt 2 }}{{\sqrt \pi }}\frac{{{{x}^{{3/2}}}{{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{1 + {{\beta }^{2}}(1 - {{x}^{2}})}},\quad {{f}_{4}}(x) = \beta \frac{{\sqrt 2 }}{{\sqrt \pi }}\frac{{{{x}^{{3/2}}}\sqrt {1 - {{x}^{2}}} {{{\text{e}}}^{{ - {{l}_{a}}x}}}}}{{1 + {{\beta }^{2}}(1 - {{x}^{2}})}}.$
Интегралы в (16) оцениваются также, как в (13). Приведем конечный результат:
$\left| {{{I}_{4}}} \right| < \frac{{C_{7}^{'}}}{{r_{a}^{{3/2}}}},\quad C_{7}^{'} = 2[C_{4}^{'}({{z}_{a}}) + C_{6}^{'}({{z}_{a}})],$
$C_{4}^{'}({{z}_{a}}) = 2\left\{ {C_{3}^{'} + \frac{{\sqrt 2 }}{{\sqrt \pi }}\left( {\frac{3}{2} + {{l}_{a}} + \frac{\pi }{2}\left| {{{z}_{a}}} \right| + 2{{\beta }^{2}}} \right)} \right\},$
$C_{6}^{'}({{z}_{a}}) = 2\left\{ {C_{5}^{'} + \beta \frac{{\sqrt 2 }}{{\sqrt \pi }}\left( {\frac{3}{2} + \frac{\pi }{2} + {{l}_{a}} + \left| {{{z}_{a}}} \right| + 2{{\beta }^{2}}} \right)} \right\},$
а $C_{3}^{'}$, $C_{5}^{'}$ – постоянные, которыми при $x \in [0,1]$ равномерно ограничены соответственно функции

${{f}_{3}}(x)cos({{z}_{a}}\sqrt {1 - {{x}^{2}}} ),\quad {{f}_{3}}(x)sin({{z}_{a}}\sqrt {1 - {{x}^{2}}} ),$
${{f}_{4}}(x)cos({{z}_{a}}\sqrt {1 - {{x}^{2}}} ),\quad {{f}_{4}}(x)sin({{z}_{a}}\sqrt {1 - {{x}^{2}}} ).$

Окончательная оценка для ${{I}_{{r(u)}}}$ имеет вид

$\left| {{{I}_{{r(u)}}}} \right| < \frac{{C{\text{''}}({{z}_{a}})}}{{r_{a}^{{3/2}}}},\quad C{\text{''}}({{z}_{a}}) = B\frac{\pi }{2}\{ C_{2}^{'} + C_{7}^{'}({{z}_{a}})\} ,$
т.е. ${{I}_{{r(u)}}} = O(1{\text{/}}r_{a}^{{3/2}})$ при ${{r}_{a}} \gg 1$ и любом фиксированном ${{z}_{a}} \leqslant 0$. Аналогичная оценка имеет место для ${{I}_{{r(l)}}}$ при тех же ${{r}_{a}}$ и ${{z}_{a}}$. Рассуждая так же, как и выше, можно получить, что ${{I}_{{(u)}}} = O(1{\text{/}}r_{a}^{{5/4}})$, ${{I}_{{(u)}}} = O(1{\text{/}}r_{a}^{{5/4}})$, где ${{r}_{a}} \gg 1$ и ${{z}_{a}} \leqslant 0$ фиксировано.

В силу полученных ранее равномерных по ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$ оценок интегралы $I_{z}^{{(m)}}$, $I_{r}^{{(m)}}$, ${{I}^{{(m)}}}$ стремятся к нулю при ${{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty $. Покажем, что все интегралы по берегам разреза также стремятся к нулю при ${{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty $. Для них полученные ранее оценки уже неприменимы, так как они неравномерны относительно ${{z}_{a}}$. Получим равномерные. Например, для ${{I}_{{z(u)}}}$ имеем

$\left| {{{I}_{{z(u)}}}} \right| \leqslant \frac{{\pi B}}{2}\int\limits_0^1 \left| {H_{0}^{{(2)}}(x{{r}_{a}})} \right|xdx = \frac{{\pi B}}{{2r_{a}^{2}}}\int\limits_0^{{{r}_{a}}} \left| {H_{0}^{{(2)}}(s)} \right|sds \leqslant \frac{{C{\text{'}}}}{{\sqrt {{{r}_{a}}} }} \to + 0,\quad {{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty .$
Аналогичным образом доказывается, что и остальные интегралы по берегам разреза равномерно стремятся к нулю. Вычисляя вычеты в (11), (12), для ${{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}})$, $U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}})$ и $U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}})$ будем иметь

${{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}}) = - i\frac{{{{\pi }^{2}}B}}{{{{\beta }^{2}}}}H_{0}^{{(2)}}\left( {\frac{{{{r}_{a}}}}{\beta }\sqrt {1 + {{\beta }^{2}}} } \right){{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{5/4}}}}} \right),$
$U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - i\frac{{{{\pi }^{2}}B}}{{{{\beta }^{3}}}}H_{0}^{{(2)}}\left( {\frac{{{{r}_{a}}}}{\beta }\sqrt {1 + {{\beta }^{2}}} } \right){{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{5/4}}}}} \right),$
$U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}}) = i{{\pi }^{2}}B\frac{{\sqrt {1 + {{\beta }^{2}}} }}{{{{\beta }^{3}}}}H_{1}^{{(2)}}\left( {\frac{{{{r}_{a}}}}{\beta }\sqrt {1 + {{\beta }^{2}}} } \right){{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right).$

Заметим, что

при ${{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty $. Отсюда и из доказанного выше следует, что $\left| {{{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}})} \right| \to + 0$ при ${{(r_{a}^{2} + z_{a}^{2})}^{{1/2}}} \to + \infty $.

Используя асимптотику функций Ханкеля на бесконечности [2], получаем

${{\Phi }^{{(a)}}}({{r}_{a}},{{z}_{a}}) = - i{{\pi }^{{3/2}}}B\frac{{\sqrt 2 }}{{\sqrt {{{r}_{a}}} }}\frac{{{{{\text{e}}}^{{i(\pi /4 - {{r}_{a}}\sqrt {1 + {{\beta }^{2}}} {{\beta }^{{ - 1}}})}}}}}{{{{\beta }^{{3/2}}}{{{(1 + {{\beta }^{2}})}}^{{1/4}}}}}{{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{5/4}}}}} \right),$
$U_{z}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - i{{\pi }^{{3/2}}}B\frac{{\sqrt 2 }}{{\sqrt {{{r}_{a}}} }}\frac{{{{{\text{e}}}^{{i(\pi /4 - {{r}_{a}}\sqrt {1 + {{\beta }^{2}}} {{\beta }^{{ - 1}}})}}}}}{{{{\beta }^{{5/2}}}{{{(1 + {{\beta }^{2}})}}^{{1/4}}}}}{{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{5/4}}}}} \right),$
$U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}}) = - {{\pi }^{{3/2}}}B\frac{{\sqrt 2 }}{{\sqrt {{{r}_{a}}} }}\frac{{{{{(1 + {{\beta }^{2}})}}^{{1/4}}}}}{{{{\beta }^{{5/2}}}}}{{{\text{e}}}^{{i(\pi /4 - {{r}_{a}}\sqrt {1 + {{\beta }^{2}}} {{\beta }^{{ - 1}}})}}}{{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){{\beta }^{{ - 1}}}}}} + O\left( {\frac{1}{{r_{a}^{{3/2}}}}} \right).$
Поэтому третье условие (1) также выполняется. Таким образом, доказана следующая

Теорема. Интеграл (2) является решением задачи (1) при ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$.

Для действительных безразмерных перемещений при ${{r}_{a}} \gg 1$ будем иметь следующие выражения:

$\begin{gathered} u_{r}^{{(a)}}({{r}_{a}},{{z}_{a}},t) = {\text{Re}}\{ U_{r}^{{(a)}}({{r}_{a}},{{z}_{a}}){{{\text{e}}}^{{i\omega t}}}\} = - {{\pi }^{{3/2}}}B\frac{{\sqrt 2 }}{{\sqrt {{{r}_{a}}} }}\frac{{{{{(1 + {{\beta }^{2}})}}^{{1/4}}}}}{{{{\beta }^{{5/2}}}}}{{{\text{e}}}^{{({{z}_{a}} - {{l}_{a}}\sqrt {1 + {{\beta }^{2}}} ){\kern 1pt} {{\beta }^{{ - 1}}}}}} \times \\ \times \;cos\left( {\omega t - \frac{{{{r}_{a}}}}{\beta }\sqrt {1 + {{\beta }^{2}}} + \frac{\pi }{4}} \right) + {{f}_{2}}({{r}_{a}},{{z}_{a}},t), \\ \end{gathered} $
где $\left| {{{f}_{1}}({{r}_{a}},{{z}_{a}},t)} \right| \leqslant {{C}_{1}}{\text{/}}r_{a}^{{5/4}}$, $\left| {{{f}_{2}}({{r}_{a}},{{z}_{a}},t)} \right| \leqslant {{C}_{2}}{\text{/}}r_{a}^{{3/2}}$ для любого $t > 0$ и любого фиксированного ${{z}_{a}} \leqslant 0$.

Возвращаясь к размерным переменным, получаем

${{u}_{z}}(r,z,t) = \sqrt {\frac{{2\pi }}{{{{k}_{R}}r}}} \frac{{{{p}_{0}}}}{\rho }{{\left( {\frac{{l\omega }}{g}} \right)}^{2}}\frac{{{{\omega }^{2}}}}{g}{{{\text{e}}}^{{{{\omega }^{2}}z/g - l{{k}_{R}}}}}sin\left( {\omega t - {{k}_{R}}r + \frac{\pi }{4}} \right) + {{g}_{1}}(r,z,t),$
${{u}_{r}}(r,z,t) = \sqrt {\frac{{2\pi {{k}_{R}}}}{r}} \frac{{{{p}_{0}}}}{\rho }{{\left( {\frac{{l\omega }}{g}} \right)}^{2}}{{{\text{e}}}^{{{{\omega }^{2}}z/g - l{{k}_{R}}}}}cos\left( {\omega t - {{k}_{R}}r + \frac{\pi }{4}} \right) + {{g}_{2}}(r,z,t),$
где ${{g}_{k}}(r,z,t)$ имеют ту же оценку по переменной $r$, что и ${{f}_{k}}({{r}_{a}},{{z}_{a}},t)$ по переменной ${{r}_{a}}$ $(k = 1,2)$.

Список литературы

  1. Bromwich T.J.I`A. On the influence of gravity on elastic waves and, in particular, on the vibrations of an elastic Globe // Proc. London Math. Soc. 1898. V. 30. P. 98–120.

  2. Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Физматлит, 1963.

  3. Ильясов Х.Х., Кравцов А.В., Кузнецов С.В., Секерж-Зенькович С.Я. Внешняя пространственная задача Лэмба. Распределенная по поверхности гармоническая нагрузка // Изв. РАН. Механ. твердого тела. 2016. № 1. С. 50–56.

  4. Кравцов А.В., Кузнецов С.В., Секерж-Зенькович С.Я. Конечноэлементные модели в задаче Лэмба // Изв. РАН. Механ. твердого тела. 2011. № 6. С. 160–169.

  5. Свешников А.Г., Тихонов А.Н. Теория функций комплексной переменной. М.: Наука, 1967.

Дополнительные материалы отсутствуют.