Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике Рассмотрим плоское движение материальной точки в гравитационном поле Земли. Если ее движение происходит вне атмосферы, то основной силой, дей- ствующей на точку, будет сила притяжения Земли. Поэтому движение рассматриваемой точки будет описываться уравнением где $\boldsymbol{r}$-радиус-вектор точки, $\mu$-гравитационная постоянная, $\boldsymbol{\varepsilon} \boldsymbol{f}$ – вектор возмущающих сил, $\varepsilon$-малый параметр (обозначения даны на рис. 32). Первое слагаемое, стоящее в правой части уравнения (8.1), – это напряженность гравитационного поля одного притягивающего центра. В действительности, силы притяжения Земли несколько отличаются от этой величины, поскольку Земля не представляет собой точечную массу. Разность вектора напряженности реального поля Земли и величины $\left(\mu / r^{3}\right) \boldsymbol{r}$ является одной из составляющих вектора $\varepsilon f$. Кроме того, к числу составляющих этого вектора мы отнесем притяжение Луны и планет, сопротивление атмосферы и т. д. Может оказаться также, что аппарат (который мы схематизируем как материальную точку) находится под дей. Рис. 32. ствием тяги двигателя «малой тяги». В этом случае ускорение, которое создается этим двигателем, мы также отнесем к числу возмущающих сил. Перепишем уравнение (8.1) в проекциях на оси полярной системы координат Здесь $v_{r}$ и $v_{\varphi}$ – проекции вектора скорости на направление радиуса-вектора и трансверсаль Величина $\boldsymbol{h}$, определенная равенством носит название секториальной скорости. Используя равенства (8.3) и (8.4), легко найдем, что Вводя обозначение можно выражение угловсй скорости $d \varphi / d t$ записать так: Рассмотрим первое из уравнений (8.2). Заметим прежде всего, что, используя (8.5), его можно представить в виде Составим теперь уравнение, которому удовлетворяет функция $u(\varphi)$. Для этого проделаем следующие вычисления: Используя эти вспомогательные выражения и уравнение (8.6), получим Составим далее уравнение, которому удовлетворяет величина $h$; так как то второе из уравнений системы (8.2) мы можем переписать в виде Переходя к переменному $\varphi$, найдем Из этих двух выражений получим искомое уравнение Используя в свою очередь уравнение (8.8), мы можем исключить из уравнения (8.7) производную $d h / d \varphi$ Уравнения (8.8) и (8.9) образуют замкнутую систему уравнений относительно $u$ и $h$. В результате ее интегрирования мы получаем элементы орбиты как функции полярного угла ч. Для того чтобы перейти к аргументу $t$, необходимо еще вычислить квадратуру При $\varepsilon=0$ система уравнений (8.8), (8.9) имеет точное решение где $x$ и $h_{0}$ – некоторые постоянные. где $x$ и $y$ – неизвестные функции угла $\varphi$. Из (8.12) находим Дифференцируя (8.13), в силу уравнений (8.8) и (8.9) получим $\frac{d x}{d \varphi}=\frac{1}{x}\left[-x \cos y \frac{d}{d \varphi}\left(\frac{\mu}{h^{2}}\right)-\frac{d u}{d \varphi} \varepsilon\left(\frac{f_{r}}{u^{2} h^{2}}+\frac{d u}{d \varphi} \frac{f_{\varphi}}{u^{3} h^{2}}\right)\right]=$ и окончательно Продифференцируем теперь равенство (8.14) откуда после замены (9.12) получим или окончательно Система уравнений (8.8), (8.15) и (8.16) может быть представлена в виде где Если функции $f_{\varphi}, f_{r}$ являются периодическими функциями полярного угла $\varphi$, го система уравнений (8.17) относится к типу уравнений, изучаемому в этой главе. Она содержит две медленные переменные $x$ и $h$ и две быстрые $y$ и $\varphi$, поэтому для ее исследования могут быть использованы изложенные выше методы анализа. При изучении этой системы может встретиться вся гамма случаев, изученных в этой главе. В системе (8.17) могут возникать резонансные ситуации. В том случае, когда возмущающие силы не зависят от полярного угла $\varphi$, система (8.17) – это система с одной вращающейся фазой. Укороченные уравнения Ван-дер-Поля можно записать для нее в следующем виде: Переход от системы (8.17) к системе (8.18) связан с операцией усреднений по быстрой переменной $y$. Заметим, что $\bar{F}_{i}$ в общем случае, когда $f_{i}$ зависит от $x, h$ и $y$, – это некоторые довольно сложные квадратуры. Таким образом, переход от системы (8.17) к (8.18) приводит к отделению медленных переменных, что позволяет при той же точности проводить интегрирование с большим шагом и, следовательно, уменьшить число шагов численного интегрирования. В то же время на каждый шаг численного интегрирования в системе (8.18) затрачивается значительно большее число машинных операций, поскольку он требует на каждом шаге вычисления квадратур. Отсюда следует, что вопрос о целесообразности перехода от одной системы к другой в каждом отдельном случае должен быть предметом специального исследования. Примечание. Метод численного интегрирования, по идее близкий к изложенному, но не использующий непосредственно приемов теории усреднений движения, был разработан еще в 1960 г. Г. П. Таратыновой *). Метод осреднений в его стандартном виде в траекторных задачах теории искусственных спутников Земли, по-видимому, впервые был применен Лассом и Лореллом **). Рассмотрим снова систему (8.17) и предположим, что требования точности определяют величину шага по фазовым переменным. Пусть шаг $\Delta x \sim \delta, \Delta y \sim \delta, \Delta h \sim \delta$. Тогда шаг по времени $\Delta t$, если принять простейшую разностную схему, определится из равенства которое справедливо для достаточно малых значений $\varepsilon$. Таким образом, число шагов $N_{1}$ определится так: где $T$ – промежуток интегрирования. Рассмотрим теперь систему (8.18) и поставим себе ту же задачу: вычислить $x(T)$ и $h(T)$. где Обозначая через $k_{2}$ число машинных операций, необходимых для реализации одного шага численного интегрирования системы (8.18), мы получаем такую оценку $R_{2}$ – числа операций: Вопрос о целесообразности замены в процессе численного интегрирования системы (8.17) системой (8.18) может быть в дан ном случае решен сравнением чисел $R_{1}$ и $R_{2}$. Примечания. 1. Для изучения задачи, рассмотренной в данном параграфе, мы составили уравнения для величин $x$ и $h$, которые в случае $\varepsilon=0$ являются интегралами изучаемой системы уравнений. Такие величины в небесной механике носят название оскулирующих элементов. Выбранные оскулирующие элементы были очень удобны для демонстрации процедуры применения методов осреднения. Однако уравнения для этих элементов оказались довольно сложными: при усреднении возникала проблема вычисления сложных вспомогательных квадратур. Поэтому в прикладных расчетах используют обычно другие оскулирующие элементы.
|
1 |
Оглавление
|