- •Оглавление
- •Введение
- •Методические указания к выполнению расчетно-графических работ
- •О вычислениях с приближенными числами
- •1 Теоретическое введение. Движение искусственных спутников земли
- •1.1 Системы отсчета в космической съемке
- •1.2 Невозмущенное движение искусственных спутников Земли
- •1.2.1 Законы Кеплера и элементы орбиты
- •1.2.2 Дифференциальные уравнения невозмущенного движения искусственных спутников Земли
- •1.3 Возмущенное движение исз
- •1.3.1 Возмущения, обусловленные отличием гравитационного поля Земли от центрального
- •1.3.2 Возмущения, обусловленные сопротивлением атмосферы
- •1.4 Трасса космического летательного аппарата
- •1.5 Выбор орбит для дистанционного зондирования
- •2 Расчетная часть ргр Содержание работы
- •Примеры выполнения заданий
- •2.1 Вычисление координат и составляющих скорости исз в заданный момент съемки
- •2.2 Построение орбитального эллипса вокруг сферической Земли
- •2.3Пересчет инерциальных декартовых координат в гринвичские и в криволинейные геодезические (широту, долготу, высоту)
- •2.4 Построение трассы спутника в координатах долгота-широта
- •Контрольные вопросы
- •Варианты исходных данных к заданиям
- •Литература
1.3.2 Возмущения, обусловленные сопротивлением атмосферы
При движении по орбите, близкой к круговой, спутник вследствие сопротивления атмосферы теряет свою энергию и непрерывно спускается по спирали. Период обращения уменьшается примерно на 0,1 с за виток, что мало сказывается на результатах космической съемки в течение достаточно длительного периода. Однако высота полета постепенно уменьшается и достигает критического значения — 150 км. Ниже этой высоты КЛА без двигательной установки может существовать на орбите очень малое время. Среднее время существования картографических и ресурсных спутников на низких орбитах — 23 года.
Для расчета влияния атмосферы на движение ИСЗ существенное значение имеют форма корпуса и ориентация ИСЗ, а также плотность атмосферы.
Наиболее сложной проблемой является создание адекватной модели атмосферы, поскольку плотность атмосферы является функцией высоты полета, расположения Солнца относительно перигея орбиты, географической широты ИСЗ, солнечной активности и вращения Солнца вокруг своей оси. На плотность атмосферы влияют ее разогревание солнечным излучением, геомагнитные бури и другие факторы.
Сила сопротивления атмосферы определяется по формуле
, (1.
где СD — аэрогидродинамический коэффициент сопротивления;S — сечение миделя — наибольшее поперечное сечение ИСЗ;— плотность атмосферы;v — скорость ИСЗ. Для сферыСD = 2,102,15.
Сила сопротивления не является потенциальной и ее учет в уравнениях движения значительно усложняет их решение. Действие сопротивления атмосферы приводит к вековым возмущениям таких элементов орбит как большая полуось и эксцентриситет.
Прочие возмущающие факторы — притяжение Луны, Солнца, солнечное излучение и другие — учитываются в особо точных расчетах и измерениях, например, по изучению параметров и аномалий гравитационного поля Земли.
1.4 Трасса космического летательного аппарата
Трасса полета КЛА (проекция орбиты на поверхность вращающейся Земли) может быть рассмотрена следующим образом.
Рассчитывается широта Ви долготаL(географические) для ряда точек по формулам:
; (1.
, (1.
где J – наклонение орбиты; - начальное значение долготы восходящего узла;S(0) – звездное время в момент наблюдения спутника, выражены в градусах, то естьS(0) = 15,0411S(h); = 0,2506845 град·мин–1 – угловая скорость вращения Земли;t – текущее время полета;T – период обращения спутника, мин.
— прецессия узла орбиты за один оборот, вычисляется по формуле:
, (1.
где A = 2,378·107градкм2;p – фокальный параметр орбиты, км;
T – период обращения спутника.
Для расчета трассы задаются шагом по времени от 0,1 до 10 минут и по приведенным формулам находят координаты точек проекции. Полученные точки наносят на карту и соединяют плавной линией, которая и является трассой полета.
Обычно рассчитывают трассу только одного витка, например, первого. Трассы последующих и предыдущих витков могут быть получены путем смещения на долготу начала n-го витка
, (1.
где L1 – долгота начала первого витка;n – номер витка;ae – большая полуось орбиты.
Расстояние Lмежду зонами обзора с двух соседних витков определяется из выражения
, (1.
где vэ — линейная скорость произвольной точки, лежащей на экваторе;
B – широта места съемки.