Журнал вычислительной математики и математической физики, 2019, T. 59, № 4, стр. 597-610
Аналитическое решение задачи Лэмба в случае предельного значения коэффициента Пуассона
Х. Х. Ильясов 1, 2, *, А. В. Кравцов 1, 2, **, С. В. Кузнецов 1, 2, ***, С. Я. Секерж-Зенькович 1, 2, ****
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
Аннотация
Рассматривается задача Лэмба для силы, приложенной к границе упругого полупространства в случае, когда коэффициент Пуассона принимает предельное значение 1/2. Решение представляется в виде контурных интегралов, для которых приводятся асимптотические оценки при больших значениях полярного радиуса. Библ. 5.
1. ПОСТАНОВКА ЗАДАЧИ. ИНТЕГРАЛЬНОЕ ПРЕДСТАВЛЕНИЕ РЕШЕНИЯ
Пусть упругая среда занимает полупространство. В случае малых относительных перемещений уравнения Ламэ имеют вид
Согласно теории предельного перехода будем считать, что коэффициент Пуассона $\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$ примут соответственно вид
Представим вектор перемещений в виде $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 $ – заданные постоянные величины, получаем скалярную задачу
Так как нагрузка не зависит от угловой переменной $\theta $, то периодическое по времени решение будем искать в виде
Для комплексной амплитуды $\Phi (r,z)$ получимДалее считаем, что функция $f(r)$ имеет вид
Введем безразмерные переменные и безразмерную комплексную амплитуду по формулам
(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 ,$(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 ,$Покажем, что интеграл (2) удовлетворяет уравнению в (1) и граничному условию при $z = 0$. Для этого убедимся, что данный интеграл можно дважды дифференцировать по параметрам ${{r}_{a}}$ и ${{z}_{a}}$ под знаком интеграла. Представим данный интеграл в виде суммы пяти интегралов: по отрезкам, лучу и полуокружностям. Рассмотрим сначала интеграл по лучу ${{\gamma }_{5}}$ ($\kappa = x \geqslant {{\kappa }_{R}} + {{\rho }_{2}}$). Подынтегральная функция
Покажем, что интегралы от $\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 $ – любое положительное число. Имеем
Таким образом, интеграл вдоль ${{\gamma }_{5}}$ от $h({{r}_{a}},{{z}_{a}},x)$ имеет непрерывные частные производные по ${{r}_{a}}$, ${{z}_{a}}$ первого и второго порядков и они равны интегралам от соответствующих частных производных функции $h({{r}_{a}},{{z}_{a}},x)$.
Рассмотрим функцию трех комплексных переменных
Аналогичным образом доказывается, что интегралы вдоль ${{\gamma }^{{(u)}}}$, ${{\gamma }_{3}}$ также можно дифференцировать два раза по параметрам ${{r}_{a}} > 0$ и ${{z}_{a}} \leqslant 0$ под знаком интеграла. Таким образом,
Покажем, что граничное условие при ${{z}_{a}} = 0$ $({{r}_{a}} > 0)$ также удовлетворяется. Имеем
Для безразмерных комплексных амплитуд перемещений $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$, будем соответственно иметь
Обозначим через ${{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 ]$.
Введем также обозначения
(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.$(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.$(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}}} }},$Покажем, что
(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} $Представим функции Ханкеля в виде
где ${{N}_{m}}(\kappa )$ – функция Неймана индекса $m$. Известно [2], что при $arg\kappa \in ( - \pi ,\pi )$(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)}}},$Интегралы от $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}}].$Подобно тому, как это было сделано выше, можно показать, что
Функции $H_{m}^{{(2)}}({{r}_{a}}\kappa )$ аналитичны в окрестности точки $\kappa = 1$, а следовательно, ограничены там, т.е. для любого $s \in [ - \pi ,\pi ]$, любого достаточно малого ${{\rho }_{5}}$ и фиксированного ${{r}_{a}}$
Перейдем в равенствах (10) к пределу при $R \to + \infty $, ${{\rho }_{4}} \to + 0$, ${{\rho }_{5}} \to + 0$. Учитывая, что
(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$ имеем
На сегменте $[1{\text{/}}\sqrt {{{r}_{a}}} ,\;1]$ при ${{r}_{a}} \gg 1$ значение аргумента $x{{r}_{a}}$ функции Ханкеля $H_{0}^{{(2)}}(x{{r}_{a}})$ также много больше единицы. Поэтому для всех $x$ из указанного сегмента функцию Ханкеля можно представить асимптотической формулой [2]. Для ${{I}_{2}}$ будем иметь
(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 - $(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\},$(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).$Рассуждая аналогичным образом, получаем, что четвертый интеграл в (13) имеет такую же оценку, как и ${{I}_{{21}}}$, а второй и третий, обозначим их ${{I}_{{2k}}}$, $k = 2,3$, оценку
Оценим интегралы ${{I}_{{r(u)}}}$, ${{I}_{{r(l)}}}$. При этом подробных пояснений делать не будем. Представим ${{I}_{{r(u)}}}$ в виде
(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 - $Окончательная оценка для ${{I}_{{r(u)}}}$ имеет вид
В силу полученных ранее равномерных по ${{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)}}}$ имеем
Заметим, что
при ${{(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], получаем
Теорема. Интеграл (2) является решением задачи (1) при ${{r}_{a}} > 0$, ${{z}_{a}} \leqslant 0$.
Для действительных безразмерных перемещений при ${{r}_{a}} \gg 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.
Градштейн И.С., Рыжик И.М. Таблицы интегралов, сумм, рядов и произведений. М.: Физматлит, 1963.
Ильясов Х.Х., Кравцов А.В., Кузнецов С.В., Секерж-Зенькович С.Я. Внешняя пространственная задача Лэмба. Распределенная по поверхности гармоническая нагрузка // Изв. РАН. Механ. твердого тела. 2016. № 1. С. 50–56.
Кравцов А.В., Кузнецов С.В., Секерж-Зенькович С.Я. Конечноэлементные модели в задаче Лэмба // Изв. РАН. Механ. твердого тела. 2011. № 6. С. 160–169.
Свешников А.Г., Тихонов А.Н. Теория функций комплексной переменной. М.: Наука, 1967.
Дополнительные материалы отсутствуют.
Инструменты
Журнал вычислительной математики и математической физики