Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике В данном классе задач разложения в рассматриваемой области имеют особенности, которые не содержатся в точном решении. Более того, в членах высших порядков особенности не только сохраняются, но и становятся более выраженными. 2.4.1. Сдвиг особенности В качестве первого примера указанного класса рассмотрим задачу, изучавшуюся Лайтхиллом [1949a]: Это уравнение имеет особенность на прямой $x=-\varepsilon y$. Граничное условие, однако, обеспечивает положительность точного решения $y(x)$ для $x \geqslant 0$; следовательно, $y(x)$ не имеет особенностей при $0 \leqslant x<\infty$. Подстановка (2.4.2) в (2.4.1), разложение по степеням $\varepsilon$ и приравнивание коэффициентов при $\varepsilon^{0}$ и $\varepsilon$ даст нам Решение задачи нулевого порядка имеет вид Подставив $y_{0}$ в (2.4.4) и решив полученное уравнение, будем иметь При $x \rightarrow 0$ имеем $y_{0}=O\left(x^{-2}\right)$, в то время как $y_{1}=O\left(x^{-5}\right)$. Таким образом, хотя и точное решение не имеет особенности при $x=0$, решение нулевого порядка имеет особенность в точке $x=0$, и эта особенность далее становится сильней. Далее мы рассмотрим движение космического корабля массы $m$ в гравитационном поле двух фиксированных притягивающих центров. Масса $M_{e}$ Земли много больше массы $M_{m}$ Луны. В де- картовой прямоугольной системе координат, показанной на рис. 2.7 , безразмерные уравнения движения имеют вид ${ }^{1}$ ) где Для введения безразмерных расстояния и времени использованы соответственно расстояние $d$ между притягивающими центрами и величина где $G$-постоянная всемирного тяготения. Лагерстром и Кеворкян [1966б] изучали эту задачу при начальных условиях где $h$-полная энергия корабля. Подстановка (2.4.12) и (2.4.13) в (2.4.7)-(2.4.9) и приравнивание коэффициентов при одинаковых степенях $\mu$ дают Решения этих уравнений при начальных условиях (2.4.10) и (2.4.11) имеют вид Таким образом, приведенное выше разложение при $x \rightarrow 1$ нарушается, поскольку $t_{1}$ имеет логарифмическую особенность. Можно показать, что в высших приближениях особенности в окрестности $x=1$ усиливаются. В самом деле, 2.4.3. Термоупругие поверхностные волны Рассмотрим влияние теплопроводности на распространение волн по поверхности изотропного упругого полупространства. Вместо закона теплопроводности Фурье будем использовать модифицированный закон Максвелла, чтобы учесть то малое время, которое необходимо для установления стационарной теплопроводности после внезапного возникновения градиента температуры в твердом теле. Будем предполагать, таким образом, что поток тепла $\mathbf{h}$ определяется соотношением Здесь $\theta$ означает изменение исходной абсолютной температуры $\theta_{0}$, $k$-коэффициент теплопроводности, $\tau_{0}$-время термической релаксации. В уравнении (2.4.21) предполагается, что тепловые сигналы распространяются не мгновенно, а имеют конечную скорость распространения. В наших рассмотрениях здесь мы будем следовать Найфэ и Неммат-Нассеру [1971]. Поскольку материал предполагается изотропным, мы будем рассматривать двумерное движение в плоскости $(x, y)$, обозначая смещения соответственно через $u$ и $v$. Ось $x$ лежит на свободной поверхности, а ось $y$ перпендикулярна к ней и направлена внутрь тела. Пусть $\beta^{2}=(\lambda+2 \mu) / \mu$, где $\lambda$ и $\mu$-коэффициенты упругости Ламе, пусть $b=[2+(3 \lambda / \mu)] \alpha \theta_{0}$, где $\alpha-$ коэффициент линейного теплового расширения, и пусть $g=\alpha(3 \lambda+2 \mu) / \rho c_{v}$, где $\rho$-плотность материала, $c_{v}$-удельная теплоемкость при постоянном объеме. Тогда совместные уравнения теории упругости и теплопроводности запишутся в виде В приведенных уравнениях для введения безразмерных времени и длины использованы соответственно величины $1 / \omega^{*}$ и $v_{p} / \omega^{*}$, где На свободной поверхности обращаются в нуль нормальные и касательные напряжения и градиент температуры, т. е. имеем при $y=0$ Решение задачи (2.4.22)-(2.4.28) будем искать в виде где $\alpha$ имеет положительную вещественную часть. Тогда волновая скорость будет задаваться равенством $c=\omega / \operatorname{Re} q$, а коэффициент затухания – равенством $s=\operatorname{Im} q$. Подстановка (2.4.9) в $(2.4 .22)-(2.4 .25)$ и приравнивание нулю определителя линейной системы уравнений относительно $a_{k}$ дают следующие три решения для $\alpha$, выраженные через $q$ и $\omega$ : Здесь принято обозначение $\varepsilon=b g / \beta^{2}$. Для каждого $\alpha$ существует собственный вектор вида и вида Подставим собственные функции, соответствующие собственным векторам (2.4.33) и (2.4.34), в уравнения (2.4.26) – (2.4.28) и приравняем нулю определитель системы уравнений относительно $a_{1 k}$. Устремляя $\omega \rightarrow \infty$, получаем где Положив в (2.4.35) и (2.4.36) $\varepsilon=0$, можно получить классическую скорость волн Рэлея. При $\beta^{2}=3$ (соответствует коэф- фициенту Пуассона $1 / 4$ ) имеем $c^{2}=0,2817$. Для малого $\varepsilon$ можно попытаться построить такое разложение для $c^{2}$, первый член которого совпадает с решением Рэлея. Итак, положим Подставив (2.4.37) в (2.4.35) и (2.4.36), разложив по степеням $\varepsilon$ и приравняв коэффициенты при одинаковых степенях в обеих частях, получим где Из равенства (2.4.38) видно, что при значениях $\tau$, близких к $c_{R}^{-2}$, коэффициент $c_{1}$ становится неограниченным, и разложение (2.4.37) нарушается. Можно показать, что при $\tau=c_{R}^{2}$ высшие приближения более сингулярны, чем второй член. При $\tau>c_{R}^{2}$ разложение (2.4.37) следует модифицировать; в противном случае величина $A^{2}$ окажется отрицательной, а в силу (2.4.32) станет отрицательной действительная часть одного из слагаемых, составляющих $\alpha$. Модифицированное разложение имеет вид Функции $c_{1}$ и $c_{2}$ можно определить, подставив (2.4.40) в (2.4.35) и приравняв коэффициенты при $\varepsilon$ и $\varepsilon^{2}$ в обеих частях. Проделав это, получим где $M$ является корнем уравнения Разложение (2.4.40) является особенным при $\tau=c_{R}^{-2}$, поскольку в этой точке коэффициент при $M^{2}$ обращается в нуль. Можно показать, что в высшем приближении эта особенность усложняется. Таким образом, первое из вышеприведенных разложений справедливо при $\tau<c_{R}^{-2}$, второе-при $\tau>c_{R}^{2}$; оба разложения нарушаются вблизи $\tau=c_{R}^{-2}$. Разложение, пригодное в окрестности этой особенности, построено в п. 4.1.6 с помощью метода сращивания асимптотических разложений. 2.4.4. Задача с точкой возврата В качестве последнего примера рассмотрим асимптотические разложения решений уравнения при больших $\lambda$. В области $|x|<1$ решения этого уравнения имеют колебательный характер, а при $|x|>1$-экспоненциальный. Это обстоятельство наводит на мысль о разложении вида где Подставив это разложение в (2.4.43) и приравняв коэффициенты при одинаковых степенях $\lambda$, получим Решениями этих уравнений являются функции Следовательно, и где $a_{i}$ и $b_{i}$-постоянные. врата. Неравномерность в этом примере возникла вследствие того, что решения представлены в элементарных функциях, а именно в экспоненциальных и тригонометрических функциях. Из разложений видно, что при переходе через значения $|x|=1$ характер решения меняется с колебательного на экспоненциальный. Следовательно, для представления решений нужны функции, имеющие аналогичное качественное поведение. Подходящими для этой цели являются функции Эйри (п. 7.3.1).
|
1 |
Оглавление
|