Методика определения конечно-элементной модели гравитационного поля Земли тема диссертации и автореферата по ВАК РФ 25.00.32, кандидат технических наук Елагин, Александр Викторович

  • Елагин, Александр Викторович
  • кандидат технических науккандидат технических наук
  • 2007, Новосибирск
  • Специальность ВАК РФ25.00.32
  • Количество страниц 252
Елагин, Александр Викторович. Методика определения конечно-элементной модели гравитационного поля Земли: дис. кандидат технических наук: 25.00.32 - Геодезия. Новосибирск. 2007. 252 с.

Оглавление диссертации кандидат технических наук Елагин, Александр Викторович

Введение.

1 Краевые задачи теории потенциала и методы их решения.

1.1 Потенциалы тяготения и силы тяжести.

1.2 Типы краевых задач.

1.3 Методы решения краевых задач.

1.3.1 Общие сведения о приближенных методах решения.

1.3.2 Метод взвешенных невязок.

1.3.3 Метод коллокаций.

1.3.4 Метод конечных разностей.

1.3.5 Метод коллокаций с подобластями.

1.3.6. Метод Бубнова - Галеркина.

1.3.7 Метод наименьших квадратов.

1.3.8 Слабая формулировка краевой задачи.

1.3.9 Метод конечных элементов.

1.3.10 Метод граничных интегральных уравнений.

1.3.11 Решение уравнения Лапласа методом разделения переменных

1.3.12 Метод точечных масс.

1.4 Обоснование выбора метода конечных элементов.

2 Методика определения конечно-элементной модели гравитационного поля Земли по измерениям.

2.1 Постановка задачи и доказательство единственности её решения.

2.2 Дискретизация области определения геопотенциала.

2.3 Алгоритм поиска адреса конечного элемента.

2.4 Преобразование координат и производных от возмущающего потенциала из общеземной в локальную и субпараметрическую системы координат.

2.5 Аппроксимация возмущающего потенциала внутри субпараметрического конечного элемента алгебраическими полиномами.

2.6 Определение базисных функций субпараметрического конечного конечного элемента.

2.7 Вычисление коэффициентов и свободных членов линейных уравнений.

2.8 Определение граничных условий в спутниковом шаровом слое и на земной поверхности.

2.8.1 Построение траекторий движения ИСЗ.

2.8.2 Вариационные уравнения для производных от векторов положения и скорости ИСЗ по узловым параметрам конечных элементов.

2.8.3 Вариационные уравнения для изохронных производных.

2.8.4 Аналитические формулы Суханова А. А. для вычисления изохронных производных.

2.8.5 Параметрические уравнения поправок для дальностей и скорости изменения дальностей.

2.8.6 Параметрические уравнения поправок для градиентометричес-ких измерений.

2.8.7 Определение граничных условий на земной поверхности.

3 Методика определения конечно-элементной модели гравитационного поля Земли по заданному разложению геопотенциала в ряд по шаровым функциям.

3.1 Определение геопотенциала тяготения и его производных с помощью моделей разложения в ряд по шаровым функциям.

3.2 Определение нормального потенциала силы тяжести и его производных.

3.3 Вычисление возмущающего потенциала и его производных.

3.4 Построение конечно-элементной модели.

4 Описание компьютерных программ, алгоритмов и вычислительных экспериментов.

4.1 Программа построения спутниковых траекторий GENORX.

4.2 Программа вычисления потенциала тяготения и его первых и вторых производных GRAVP.

4.3 Программа вычисления высот квазигеоида GEOID.

4.4 Система высот на основе модели EGM96.

4.5 Программа вычисления уклонения отвесной линии UOL.

4.6 Программа определения возмущающего потенциала и высот квазигеоида внутри нивелирного полигона методом конечных элементов МКЕ.

4.7 Вычислительные эксперименты по определению возмущающего потенциала внутри нивелирного полигона методом конечных элементов.

4.8 Вычислительный эксперимент по построению отдельного конечного элемента с использованием модели EGM96.

Рекомендованный список диссертаций по специальности «Геодезия», 25.00.32 шифр ВАК

Введение диссертации (часть автореферата) на тему «Методика определения конечно-элементной модели гравитационного поля Земли»

Актуальность проблемы. Проблема построения глобальной высокоточной модели гравитационного поля Земли (ГПЗ) остается до сих пор нерешенной. Это связано с недостатком гравиметрической информации на большей части континентов, а также в приполярных областях. Достаточно плотной и качественной гравиметрической съёмкой покрыты лишь территории бывшего Советского Союза, США, Европы и Австралии. Кроме этого, серьёзной проблемой является согласование наземных гравиметрических сетей разных стран. Например, при сравнении гравитационных данных Советского Союза и Китая обнаружилось, что, в местах перекрытия, рассогласование этих данных огромно и достигает 70 мГал [1].

Альтиметрические измерения также не решают проблемы, так как ими не охвачены суша и приполярные зоны и по ним определяются не высоты геоида, а высоты морской поверхности (отличие достигает 1,5 м).

По траекторным спутниковым измерениям с наземных пунктов уверенно находятся лишь низкочастотные составляющие разложения ГПЗ в ряд по шаровым функциям. Следовательно, построить высокоточную модель гравитационного поля Земли только по наземным гравиметрическим, альтиметрическим и траекторным измерениям невозможно.

В настоящее время с появлением новых спутниковых градиентометриче-ских и спутник-спутниковых средств измерений, которые чувствительны к высокочастотной составляющей разложения ГПЗ в ряд по шаровым функциям, открылась возможность осуществить уточнение модели ГПЗ. Измерения уже выполняются, и будут в дальнейшем выполняться, глобально и плотно над всей поверхностью Земли на высотах от 240 до 500 км. Предлагаются и частично реализованы за рубежом три научных проекта, основанных на этих новых видах измерений [1-6]:

1. CHAMP (Challeging Minisatellite Payload) - совместный проект немецких, французских и американских ученых. Спутник запущен с космодрома Плесецк 15 июля 2000 г. Орбита круговая, высотой 450 км и с наклонением 87,27°. Для определения ГПЗ используются измерения трёхосного бортового акселерометра и спутник-спутниковые измерения «высокий - низкий» (GPS-измерения). По шестнадцатисуточным данным уже построена модель 0SU02A до пятидесятой степени и порядка. По точности эта модель, для волн длиной 800 км, на 10 % выше ранее полученных моделей, особенно в приполярных областях [7]. По трехмесячным акселерометрическим и спутник-спутниковым измерениям программы CHAMP построена также германо-французская модель длинноволнового гравитационного поля EIGEN-1S, которая по точности в два раза выше предыдущих моделей. Геоид, построенный по этой модели, характеризуется погрешностью не более 20 см. [8].

2. GRACE (Gravity Recovered and Climate Experiment) - проект, разработанный американскими учеными в Национальном управлении по аэронавтике и исследованию космического пространства (NASA). Пара искусственных спутников Земли (ИСЗ) проекта была запущена 17 марта 2002 г. с космодрома Плесецк. В апреле этого же года началось тестирование и калибровка научной аппаратуры ИСЗ. Полученные данные доказали работоспособность научной концепции проекта и подтвердили отличные характеристики измерительных устройств [9]. Спутники запущены на одну и туже близполярную орбиту высотой 500 км и на удалении друг от друга от 170 до 270 км. Расстояние между двумя ИСЗ измеряются с субмикронной точностью. На ИСЗ производятся GPS-измерения и акселерометрические измерения [10]. Погрешность измерения относительных ускорений 1 мкм/с, а погрешность GPS-измерений 1 см [11]. Предполагается, что, после обработки измерительной информации, модель статического гравитационного поля будет на несколько порядков точней существующих моделей [12].

3. GOCE (Gravity field and steady - state Ocean Circulation Explorer) - проект Европейского космического агентства [13-18]. Запуск планируется осуществить в ближайшем будущем с космодрома Плесецк. Вес спутника 1200 кг. Орбита ИСЗ близкруговая, солнечно-синхронная, с наклонением 96,5° и высотой от 240 до 270 км. Для точного определения орбиты на ИСЗ устанавливается высокоточный GPS-приемник. Основным инструментом для определения ГПЗ будет гравитационный градиентометр [19, 20], разработанный французской фирмой ONERF. Градиентометр представляет собой систему из шести трехосных электростатических акселерометров. Основой каждого акселерометра является пробная масса в 320 г, поддерживаемая в подвешенном состоянии внутри блока, вблизи его центра, с помощью электростатического поля, создаваемого электродами, размещенными на стенках блока. Перемещения этой массы ежесекундно измеряются. Результаты измерений фильтруются и преобразуются в гравитационные градиенты по восемнадцати различным направлениям. Точность измерения градиентов на уровне 1-10~3 Этвеш/^[Гц [21]. Предполагается создать новую модель ГПЗ, по которой можно будет вычислять аномалии силы тяжести с погрешностью 1 мГал, а высоты геоида с погрешностью 1 см и с пространственным разрешением 70 км [18,22].

Для будущих космических проектов по изучению гравитационного поля Земли зарубежные ученые разрабатывают лазерный дальномер со стабилизированными длинами волн [23]. Он предназначен для измерения расстояний с погрешностью 2,7 ■ Ю~9 м/-^Гц между двумя близкими ИСЗ, движущимися по одной и той же орбите. По разности расстояний за определённые промежутки времени определяются вторые производные от геопотенциала по линии спутник-спутник.

В нашей стране также ведутся научно-исследовательские работы по использованию градиентометрических и спутник-спутниковых измерений для определения гравитационного поля Земли [13, 24]. Научно-производственное объединение прикладной механики предлагает включить в Федеральную космическую программу до 2010 года проект «Геопот» [25]. Он предназначен для определения параметров ГПЗ и их изменений во времени. Масса ИСЗ 800 кг, высота орбиты 400 км, наклонение 90°, срок службы - не менее семи лет.

Однако на высотах от 240 до 500 км, где выполняются спутник-спутниковые и спутниковые градиентометрические измерения, гравитационное поле имеет более гладкую структуру, чем у поверхности Земли. Поэтому, определяемые по спутниковым измерениям коэффициенты разложения ГПЗ в ряд по шаровым функциям будут отличаться от коэффициентов, которые были бы получены, если такого же типа измерения выполнить глобально на всей поверхности Земли. Отсюда возникает пока не решенная проблема продолжения вниз.

Частично эта проблема может быть решена, если поставить и решить краевую задачу с двумя граничными поверхностями. Одна граничная поверхность должна принадлежать спутниковому шаровому слою со спутниковыми градиентометрическими и/или спутник-спутниковыми измерениями. Другая граничная поверхность - земная поверхность, с максимально возможным набором наземных типов измерений, выполненных с достаточной точностью. Наземных измерений будет явно не хватать, поэтому в некоторых районах земной поверхности придется выполнять аппроксимацию и интерполирование функционалов от потенциала силы тяжести. Неизбежным является принятие каких-либо гипотез и допущений. Например, на море, где не выполнены гравиметрические съемки, считать морскую поверхность, вычисленную при помощи аль-тиметрических измерений, поверхностью геоида, с заданным значением потенциала силы тяжести Wq. На суше необходимо использовать нивелирные, гравиметрические, градиентометрические и GPS-измерения для вычисления потенциала силы тяжести и первых производных от него.

При построении моделей гравитационного поля Земли обычно решаются различные внешние краевые задачи теории потенциала в области, ограниченной земной поверхностью и бесконечной внешней границей. До недавнего времени целью исследований было совместное определение внутренней границы области, т. е. физической поверхности Земли и внешнего ГПЗ. В настоящее время GPS-измерения и спутниковая альтиметрия позволяют, принципиально, новыми геометрическими методами решать первую часть задачи - определять геометрию земной поверхности. Актуальной проблемой остается вторая часть задачи - уточнение внешнего гравитационного поля Земли из совместной обработки наземных и новых спутник-спутниковых и спутниковых градиенто-метрических измерениий.

Гравитационное поле Земли можно описать глобально при помощи разложения в ряд по шаровым функциям, набором точечных масс или мультипо-лей, аномалиями силы тяжести и т. д. В рассмотренных выше космических проектах уже определяются и планируется дальнейшее уточнение коэффициентов Спт ,Snm разложения гравитационного поля в ряд по шаровым функциям. Недостаток такого разложения и других перечисленных способов представления гравитационного поля состоит в том, что для вычисления какой-либо характеристики поля в заданной точке на Земле или вне Земли с большой точностью требуется выполнить значительный объём вычислительной работы (просуммировать все члены ряда) и хранить большой объем информации. Кроме этого, при использовании математического аппарата гармонического анализа приходится относить выполненные на земной поверхности или в космосе измерения к сфере или решать проблему редуцирования измерений на сферу.

Избежать редуцирования и существенно сократить затраты на вычисления возможно, если разбить область определения ГПЗ на малые подобласти -конечные элементы и представлять в них гравитационный потенциал и функции от него локальными алгебраическими или тригонометрическими полиномами низкой степени. В этом случае придется сузить область определения геопотенциала и определять его в слое между земной поверхностью и спутниковым шаровым слоем. Сужение области - с одной стороны недостаток, но с другой стороны оно может повлиять на точность определения геопотенциала, потому что геопотенциальная функция, оказывается «зажатой» между двумя близкими поверхностями с известными значениями потенциала и его производных.

Преимущество таких локальных представлений гравитационного поля Земли по сравнению с глобальными, в задаче определения орбит спутников, показано в работах [27 - 29].

Если решать краевую задачу методом конечных элементов, то проще аппроксимировать сложную форму граничной поверхности Земли. Построенная путем решения краевой задачи конечно-элементная модель ГПЗ значительно удобней для практического применения, по сравнению с другими моделями. Из всего набора конечных элементов выбираются для использования только те ячейки, которые непосредственно попадают в район практической деятельности. Разработке такой конечно-элементной модели ГПЗ и посвящено данное исследование.

Какими оптимальными по форме и размерам должны быть конечные элементы? Какого вида, и какой степени должны быть полиномы? Как по координатам заданной точки определить в каком конечном элементе она находится? Каким образом подключить к решению краевой задачи разнородные наземные измерения и новые спутник-спутниковые и спутниковые градиентомет-рические измерения и определить значения потенциала в узлах конечных элементов? Как выполнить интерполирование и найти потенциал и его первые, и вторые производные внутри конечного элемента? Эти и другие актуальные вопросы являются предметом исследования диссертации.

Степень разработанности проблемы. Решению проблемы определения ГПЗ посвятили свои труды многие отечественные и зарубежные ученые: Молоденский М. С., Юркина М. И., Еремеев В. Ф., Бровар В. В., Шимбирев Б. П., Макаров Н. П., Пеллинен JI. П., Жонголович И. Д., Грушинский Н. П., Машимов М. М., Мещеряков Г. А., Марченко А. Н., Петровская М. С., Нейман Ю. М., Бывшев В. А., Остач О. М., Непоклонов В. Б., Плешаков Д. И., Бузук В. В., Вовк И. Г., Канушин В. Ф., Суздалев А. С., Панаев Г. А., Сурнин Ю. В., Кужелев С. В., Рапп Р., Чернинг С., Краруп Т., Шварц К., Руммель Р., Хейсканен В., Мориц Г., Кау-ла У. М., Бьерхаммар А., Бурша М. и др. Однако, решаемая проблема уточнения моделей ГПЗ по новым спутник-спутниковым и спутниковым градиенто-метрическим измерениям является ещё недостаточно разработанной. Реализация проектов CHAMP и GRACE ещё не завершена, получены только первые предварительные результаты. Запуск ИСЗ программы GOCE ещё только планируется произвести и выполняются лишь модельные исследования. Трехмерные конечно-элементные модели ГПЗ, в области между земной поверхностью и спутниковым шаровым слоем, ранее не определялись.

Целью исследования являлась разработка методики определения конечно-элементной модели ГПЗ по разного вида спутниковым и наземным измерениям и по уже существующим разложениям геопотенциала в ряд по шаровым функциям.

Задачи исследования. Для достижения поставленной цели требовалось решить следующие задачи:

- постановка краевой задачи;

- доказательство единственности ее решения;

- определение граничных значений производных геопотенциала в спутниковом шаровом слое;

- определение граничных значений геопотенциала и его производных на земной поверхности;

- дискретизация области определения геопотенциала;

- замена интегральных уравнений системой линейных уравнений;

- решение системы линейных уравнений для определения параметров модели ГПЗ;

- вычисление геопотенциала и его производных по известным моделям разложения в ряд по шаровым функциям;

- определение нормального потенциала силы тяжести и его производных;

- вычисление возмущающего потенциала в узлах конечных элементов.

Объект исследования - внешнее гравитационное поле Земли в области ограниченной земной поверхностью и спутниковым шаровым слоем.

Предмет исследования - методика построения конечно-элементной модели гравитационного поля Земли.

Теоретическая, методологическая и информационная база исследований. Рассмотрены различные численные методы решения задач математической физики. Из сопоставления методов, в качестве основы для решения проблемы уточнения ГПЗ, выбран метод конечных элементов. Экспериментальные исследования выполнены методом математического моделирования. Реальное поле заменено моделью гравитационного поля EGM96 до 360-й степени и порядка. Эта модель позволяет вычислить и считать «эталонными» значения потенциала и его производных в любой точке на и вне поверхности Земли, а также моделировать измерения, которые требуется обработать по предлагаемой методике. Достоинство метода математического моделирования состоит в том, что заранее известны параметры гравитационного поля, которые должны быть получены с помощью построенной конечно-элементной модели. Этот метод применяется для проверки правильности работы новых алгоритмов и программ, предназначенных для обработки реальных измерений. В диссертационной работе основные аспекты предлагаемой методики апробированы на модельных измерениях. Обработка реальных измерений возможна в будущем после того, как появится доступ к спутниковой и наземной измерительной информации.

Научная новизна работы:

- внешний гравитационный геопотенциал определяется в области, ограниченной земной поверхностью и спутниковым шаровым слоем;

- область определения геопотенциала дискретизируется шаровыми слоями и усечёнными треугольными пирамидами (конечными элементами) внутри слоёв;

- разработан алгоритм поиска номера конечного элемента, в котором расположена точка с заданными координатами;

- получены формулы прямых и обратных преобразований координат точек и производных от возмущающего геопотенциала из общеземной системы координат в локальную систему координат субпараметрического конечного элемента;

- для определения граничных условий в спутниковом шаровом слое получены формулы связи спутник-спутниковых и спутниковых градиентометри-ческих измерений с возмущающим потенциалом в узлах конечных элементов;

- для определения граничных условий на земной поверхности получены уравнения Лапласа в криволинейной геодезической системе координат и в системе плоских прямоугольных координат Гаусса - Крюгера и геодезической высоты;

- разработана методика определения возмущающего потенциала внутри нивелирного полигона по спутниковым, гравиметрическим и градиентометри-ческим измерениям.

- получены формулы для вычисления первых и вторых производных от нормального потенциала силы тяжести по прямоугольным пространственным координатам.

Теоретическая и практическая значимость работы. Предложенная в работе методика дает возможность определить удобную для практического использования трехмерную конечно-элементную модель ГПЗ. Разработанные для решения краевой задачи алгоритмы и компьютерные программы могут быть также использованы для вычисления орбит ИСЗ, высот квазигеоида и уклонений отвесных линий, для определения возмущающего потенциала внутри нивелирного полигона по гравиметрическим и градиентометрическим измерениям.

В первом разделе диссертации приведены общие сведения о потенциалах тяготения и силы тяжести, перечислены типы краевых задач, даны общие сведения о приближенных методах решения краевых задач. Рассмотрены: метод взвешенных невязок, метод коллокаций, метод конечных разностей, метод коллокаций с подобластями, метод Бубнова - Галеркина, метод наименьших квадратов, метод конечных элементов, метод граничных интегральных уравнений, решение уравнения Лапласа методом разделения переменных, метод точечных масс. Показаны достоинства и недостатки перечисленных методов.

В результате анализа различных методов в качестве основы для исследований был выбран метод конечных элементов.

Во втором разделе рассмотрена методика определения конечно-элементной модели ГПЗ по измерениям, выполненным на земной поверхности и в спутниковом шаровом слое. Дана постановка задачи и доказана единственность её решения. Предложено дискретизировать область определения возмущающего потенциала шаровыми слоями и усеченными треугольными пирамидами (конечными элементами), построенными на основе многогранников. Получены формулы для вычисления базисных функций конечных элементов, а также коэффициентов и свободных членов линейных уравнений, из решения которых находятся значения возмущающего потенциала в узлах конечных элементов. Рассмотрены вопросы об определении граничных условий на внутренней поверхности спутникового шарового слоя и на земной поверхности.

В третьем разделе представлена методика определения конечно-элементной модели ГПЗ по заданному разложению геопотенциала в ряд по шаровым функциям. В основе методики лежат модифицированные автором формулы Пайнза С. для вычисления потенциала тяготения и его первых и вторых производных, а также формулы Молоденского М. С. для вычисления нормального потенциала. Они позволяют определять возмущающий потенциал в узлах конечных элементов и создавать конечно-элементные модели ГПЗ. Получены формулы для вычисления первых и вторых производных от нормального потенциала по направлению осей общеземной прямоугольной пространственной системы координат.

В четвертом разделе диссертации дано описание компьютерных программ, алгоритмов и вычислительных экспериментов.

Основные положения, выносимые на защиту следующие:

1) математическое описание конечно-элементной модели гравитационного поля Земли;

2) методика определения параметров конечно-элементной модели гравитационного поля Земли по измерениям из решения краевой задачи;

3) методика определения параметров конечно-элементной модели гравитационного поля Земли по заданному разложению геопотенциала в ряд по шаровым функциям.

Похожие диссертационные работы по специальности «Геодезия», 25.00.32 шифр ВАК

Заключение диссертации по теме «Геодезия», Елагин, Александр Викторович

В результате проведенных исследований в диссертации разработана ме тодика определения конечно-элементной модели ГПЗ. Рассмотрено два случая. Первый, наиболее важный - определение модели из решения краевой задачи

по измерениям. Второй - определение модели по априорно заданному разло жению геопотенциала в ряд по шаровым функциям. Второй случай реализуется

проще, однако все погрешности первоначального представления складываются

с погрешностями конечно-элементной аппроксимации. В первом случае, по сравнению с классической постановкой решения

внешней краевой задачи, область определения геопотенциала сужена. Она ог раниченна земной поверхностью и спутниковым шаровым слоем. Доказано,

что, если измерены потенциал и ускорения силы тяжести на поверхность Зем ли, а в спутниковом шаровом слое измерены ускорения силы тяготения, то

краевая задача имеет единственное решение. Показано, что решение уравнения Лапласа, с заданными граничными ус ловиями, эквивалентно решению задачи минимизации функционала (59). Ре шать задачу минимизации предпочтительнее, так как в этом случае требования

к гладкости аппроксимирующей функции понижаются. Для работы с глобальной конечно-элементной моделью необходимо рас полагать алгоритмом поиска адреса конечного элемента по заданным коорди натам точки. Такой алгоритм в диссертации разработан на основе многогран ников. Показано, что если проекциями узлов конечных элементов на единич ной сфере будут узлы икосододекаэдра и икосаэдра, то единичная сфера делит ся на два полушария (северное и южное), а каждое полушарие на однотипные

Для дальнейшей дискретизации предложено отображать «большой» сфериче 218 ский треугольник на плоскость в гномонической проекции и дискретизировать

его отображение одинаковыми «малыми» треугольниками заданного размера. В результате такого преобразования упрощается алгоритм поиска номера ко нечного элемента по координатам точки. Разработан также алгоритм определе ния номера слоя, в котором располагается точка. В качестве конечного элемента выбрана усеченная треугольная пирами да. Её размеры и форма могут варьироваться, что позволяет гибко аппроксими ровать форму земной поверхности. Если использовать такие конечные элемен ты напрямую, то для каждого элемента пришлось бы определять свои базисные

функции, зависящие от координат узлов и хранить их. Поэтому введены в рас смотрение субпараметрические конечные элементы, в которых базисные функ ции для всех элементов одинаковые, хотя возникают неудобства при пересче тах координат и производных от возмущающего потенциала из общеземной

системы координат в субпараметрическую систему и обратно. Формулы преобразования координат и производных внутри конечного

элемента осуществляют отображение усеченной треугольной пирамиды в стан дартную треугольную призму. Однако основной причиной введения стандарт ных по форме субпараметрических конечных элементов явилось то, что они

позволили получить формулы для вычисления объемных и поверхностных ин тегралов численным методом Гаусса. Для упрощения решения краевой задачи первые производные от возму щающего потенциала на внутренней границе спутникового шарового слоя

предложено определять отдельно. Спутниковый слой, также как и вся область

определения, делится на усеченные треугольные пирамиды (конечные элемен ты). В диссертации получены уравнения связи спутник-спутниковых и спутни ковых градиентометрических измерений с определяемым возмущающим по тенциалом в узлах конечных элементов и сделан вывод параметрических урав нений поправок. Для определения граничных условий на земной поверхности были полу чены уравнения Лапласа в криволинейной геодезической системе координат и в системе плоских прямоугольных координат Гаусса - Крюгера и геодезиче ской высоты. Представлена методика определения возмущающего потенциала

внутри нивелирного полигона методом конечных разностей и методом конеч ных элементов по спутниковым, гравиметрическим и градиентометрическим

измерениям. Для определения конечно-элементной модели ГПЗ по заданному разло жению геопотенциала тяготения в ряд по шаровым функциям, в качестве ис ходных взяты формулы Пайнза, не имеющие особенностей на полюсах. Они

были модифицированы в части рекуррентных формул для вычисления аналога

нормированных присоединенных функций Лежандра. Для вычисления нор мального потенциала силы тяжести была использована формула

Молоденского М.С. Путем дифференцирования этой формулы по прямоуголь ным координатам найдены выражения для вычисления первых и вторых произ водных. В результате получены формулы, которые позволяют по заданным

моделям ГПЗ вычислять возмущающий потенциал в узлах конечных элементов

и создавать конечно-элементные модели гравитационного поля Земли. Цель исследований достигнута - разработана методика определения ко нечно-элементной модели ГПЗ по разнородным спутниковым и наземным из мерениям, а также по уже существующим разложениям геопотенциала в ряд по

шаровым функциям.

Список литературы диссертационного исследования кандидат технических наук Елагин, Александр Викторович, 2007 год

1. Кабан, М. К. Новые возможности гравитационного моделирования с использованием данных спутников CHAMP и GRACE Текст. / М. К. Кабан, К. Райгбер//Физика Земли.-2005. -№ 11. -С. 101-109.

2. Forsberg, R. Arctic gravity project Text. / R. Forsberg, S. Kenyon // IUGG99: XXII General Assembly of the International Union of Geodesy and Geophysics, Birmingham, 19-24 July, 1999: Abstr. Week A. Birmingham. 1999. P. A.47.

3. Sunkel, H. G + +. G + +Text. / H. Sunkel // VGI: Osterr. Z. Vermess. Und Geoinf. -2000. 88 ,№3, P. 169-172.

4. Balmino, G. New Space missions for mapping the Earth's gravity field Text. / G. Balmino // C. r. Acad. sci. Ser. 4. 2001. 2, № 9, P. 1353-1359.

5. Han, S.C. Efficient gravity field recovery using in situ disturbing potential observables from CHAMP Text. / S.C. Han, J. Christopher, C.K. Shum // Geophysical Research Letters. 2002. - Vol. 29. - No. 16. - P. 36/1-36/4.

6. Лисов, И. GRACE: работа началась Текст. / И. Лисов // Новости космонавтики. 2002. - 12, № 11. - С. 49.

7. Heb, D. Gradiometric mit GRACE. Teil I. Fehleranalyse kunstlicher Gradiometerdaten Text. / D. Heb, W. Keller // Z. Vermessungsw. 1999.124, № 5. -РЛ 37—144.

8. Tapley, B.D. The gravity recovered and climate experiment (GRACE): Abstr. 23rd Gen. Assem. Eur. Geophys. Soc., Nice, 1998 Text. / B.D. Tapley // Ann. Geophys. 1998. 16. Suppl.-P. 367.

9. Нейман, Ю.М. Модель цветного шума градиентометра миссии GOCE Текст. / Ю.М. Нейман, Л.С. Сугаипова // Изв. вузов. Геодезия и аэрофотосъёмка. 2005. -№ 2. - С. 31-42.

10. Bruce, В. The Four Candidate Earth Explorer Core Missions: Reports for Mission Selection. Rept 1. Gravity Field and Steady state Ocean Circulation Mission Text. / B. Bruce. Noordwijk: ESTEC. 1999. 217. ESA SP. № 1233.

11. Van den IJssel, J. Simulation of GOCE gravity field recovery: Abstr. 24 Gen. Assemb. Eur. Geophys. Soc. includ. Symp. "Solid Earth Geophys. and Geod.", 1999 Text. / J. Van den IJssel, R. Koop, P. Visser // Geophys. Res. Abstr. 1999. 1. № l.P. 197.

12. Visser, P.N.A.M., The gravity field mission GOCE: a revolution for the geosciences Text. / P.N.A.M. Visser, R. Rummel, G. Balmino, H. Sunkel, J. Johannessen, M. Aguirre, P.L. Woodworth, C. Le Provost, C.C. Tscherning, R. Sabadini //

13. G International Symposium on Gravity, Geoid and Geodynamics 2000, Banff, July 31-Aug. 4,2000; GGG2000 Technical Program and Abstracts. Banff. 2000. 155.

14. Rummel, R. Gravitacios gradiometria: Eotvos Lorandtol a modern urkor-szakig Text. / R. Rummel // Magy. Geofiz. 2002. 43, № 4. P. 145-150.

15. Pail, R. GOCE gravity field processing strategy Text. / R. Pail, G. Plank // Stud. Geophys. Et geod. 2004. 48, № 2. P. 289-309.

16. Петровская, М.С. Методика использования спутниковой градиенто-метрии для изучения гравитационного поля Земли Текст. / М.С. Петровская // Космические исследования. 1999. -Т.З7, № 5. - С. 488-497.

17. Место малых космических аппаратов в решении задач Федеральной космической программы России Текст. / В.В. Алавердов, Б.В. Бодин, А.В. Головко, Г.Д. Голубев, А.Н. Мальченко // Космонавтика и ракетостроение. 2003. -2(31).

18. Медведев, П.П. Изучение топографической поверхности мирового океана Текст. / П.П. Медведев // Итоги науки и техники. Сер. Геодезия и аэросъёмка. М: ВИНИТИ, 1988. - Т.26. Проблемы морской геодезии. -С. 76-129.

19. Junkins, J. L. Investigation of finite-element representations of the geopo-tential Text. / J. L. Junkins // «А1АА Journal». 1976, 14, № 6. - P. 803-808.

20. Junkins, J. L. The finite element approach in gravity modeling Text. / J. L. Junkins, R. C. Engels // «Manuscr. Geod.». 1979, 4, № 2. - P. 185-206.

21. Парамзин, A.B. Разработка и исследование методов представления и уточнения параметров геопотенциала Текст.: Автореф. дис. канд. техн. наук / А.В. Парамзин. М.: МИИГАиК, 1984.

22. Грушинский, Н. П. Теория фигуры Земли Текст. / Н. П. Грушинский.- 2-е изд.,перераб. и доп. М.: Наука, 1976. - 512 с.

23. Молоденский, М.С. Основные вопросы геодезической гравиметрии Текст. / М.С. Молоденский // Тр. ЦНИИГАиК. 1945. - Вып.42. - С. 3-107.

24. Heiskanen, W. A. Physical Geodesy Text. / W. A. Heiskanen, H. Moritz.- W. H. Freeman and company, San Francisco and London, 1967.

25. Бровар, В. В. Теория фигуры Земли Текст. / В. В. Бровар, В. А. Магницкий, Б. П. Шимбирев; под ред. В. А. Магницкого. М.: Изд-во геодезической литературы, 1961. — 256 с.

26. Шимбирев, Б. П. Теория фигуры Земли Текст. / Б. П. Шимбирев. -М.: Недра, 1975.-432 с.

27. Антонов, В.А. Введение в теорию ньютоновского потенциала Текст. / В.А. Антонов, Е.И. Тимошкова, К.В. Холщевников. -М.: Наука, 1988.-272 с.

28. Власова, Е.А. Приближенные методы математической физики: учеб. для вузов Текст. / Е.А. Власова, B.C. Зарубин, Г.Н. Кувыркин. М.: Изд-во МГТУ им. Н.Э. Баумана, 2001. - 700 с. (Сер. Математика в техн. ун-те; Вып. XIII).

29. Зенкевич, О. Конечные элементы и аппроксимация Текст. / О. Зенкевич, К. Морган; пер. с англ. М.: Мир, 1986. - 318 с.

30. Бреббия, К. Методы граничных элементовТекст. / К. Бреббия, Ж. Теллес, JI. Вроубел; пер. с англ. М.: Мир, 1987. - 524 с.

31. Бахвалов, Н. С. Численные методы Текст. / Н. С. Бахвалов, Н. П. Жидков, Г. М. Кобельков. 2-е изд. - М.: Лаборатория базовых знаний, 2001.- 632 с.

32. Годунов, С. К. Разностные схемы (введение в теорию) Текст. / С. К. Годунов, В. С. Рябенький. -М.: Наука, 1977. 439 с.

33. Годунов, С. К. Уравнения математической физики Текст. / С. К. Годунов. -М.: Наука, 1971. 416 с.

34. Канторович, Л. В. Приближенные методы высшего анализа Текст. / Л. В. Канторович, В. И. Крылов. 5-е изд.,испр. - Л.: Физматгиз, 1962. - 708 с.

35. Марчук, Г. И. Методы вычислительной математики Текст. / Г. И. Марчук. 3-е изд., перераб. и доп. - М.: Наука, 1989. - 608 с.

36. Михлин, С. Г. Приближенные методы решения дифференциальных и интегральных уравнений Текст. / С. Г. Михлин, X. Л. Смолицкий. М.: Наука, 1965.-384 с.

37. Флетчер, К. Численные методы на основе метода Галёркина Текст. / К. Флетчер; пер. с англ. М.: Мир, 1988. - 352 с.

38. Михлин, С. Г. Вариационные методы в математической физике Текст. / С. Г. Михлин. 2-е изд.,перераб. и доп. - М.: Наука, 1970. - 512 с.

39. Зенкевич, О. Метод конечных элементов в теории сооружений и в механике сплошных сред Текст. / О. Зенкевич, И. Чанг; пер. с англ. М.: Недра, 1974.-240 с.

40. Зенкевич, О. Метод конечных элементов в технике Текст. / О. Зенкевич; пер. с англ. М.: Мир, 1975. - 541 с.

41. Оден, Дж. Конечные элементы в нелинейной механике сплошных сред Текст. / Дж. Оден; пер. с англ. М.: Мир, 1976. - 464 с.

42. Деклу, Ж. Метода конечных элементов Текст. / Ж. Деклу; пер. с франц. М: Мир, 1976. - 96 с.

43. Стренг, Г. Теория метода конечных элементов Текст. / Г. Стренг, Дж. Фикс; пер. с англ. М.: Мир, 1977. - 349 с.

44. Сегерлинд, Л. Применение метода конечных элементов Текст. / JI. Сегерлинд; пер. с англ. М: Мир, 1979. - 392 с.

45. Норри, Д. Введение в метод конечных элементов Текст. / Д. Норри, Ж. де Фриз; пер. с англ. М: Мир, 1981. - 304 с.

46. Митчелл, Э. Метод конечных элементов для уравнений с частными производными Текст. / Э. Митчелл, Р. Уэйт; пер. с англ. М.: Мир, 1981. -216с.

47. Марчук, Г.И. Введение в проекционно-сеточные методы Текст. / Г.И. Марчук, В.И. Агошков. -М.: Наука, 1981. -416 с.

48. Бате, К. Численные методы анализа и метод конечных элементов Текст. / К. Бате, Е. Вилсон; пер. с англ. М: Стройиздат, 1982.-448 с.

49. Сабоннадьер, Ж.-К. Метод конечных элементов и САПР Текст. / Ж.-К. Сабоннадьер, Ж.-Л. Кулон; пер. с франц. М.: Мир, 1989. - 190 с.

50. Шайдуров, В.В. Многосеточные методы конечных элементов Текст. / В. В. Шайдуров. М.: Наука, Гл. ред. Физ.-мат. лит., 1989. - 288 с.

51. Fredholm, I. Sur une classe d'equations fonctionelles Text. /1. Fredholm // Acta Math. 1903. No. 27. - P. 365-390.

52. Пеллинен, Л. П. Физическая геодезия Текст. / Л. П. Пеллинен, Ю. М. Нейман // Итоги науки и техники. Сер. Геодезия и аэросъемка. М.: ВИНИТИ, 1980.-Т.18.-132 с.

53. Bjerhammar, A. On the discrete boundary value problem in physical geodesy Text. / A. Bjerhammar // «Proc. Symp. Earth's Gravit. Field and Secul. Variat. Posit. Syndey, 1973». Sudney. 1974. - P.475^88.

54. Bjerhammar, A. Discrete approaches to the solution of the boundary value problem in physical geodesy Text. / A. Bjerhammar // «Boll. Geod. E sci. affini». -1975, №2.-P. 185-240.

55. Bjerhammar, A. A review of discrete methods in physical geodesy Text. / A. Bjerhammar // «Approximation Methods in Geodesy», ed. H. Moritz, H. Sunkel. Herbert Wichmann Verlag, Karlsruhe. -1978.

56. Мориц, Г. Современная физическая геодезия Текст. / Г. Мориц; пер. с англ. М.: Недра, 1983. - 392 с.

57. Мещеряков, Г. А. Задачи теории потенциала и обобщенная Земля Текст. / Г. А. Мещеряков. М.: Наука, 1991. - 216 с.

58. Елагин, А. В. Методика определения трёхмерной конечно-элементной модели гравитационного поля Земли Текст. / А. В. Елагин // Геодезия и картография. 2007 . - № 7. - С. 37-41.

59. Вахрамеева, JI. А. Картография Текст.: учебник для вузов / Л. А. Вахрамеева. М.: Недра, 1981. - 224 с.

60. Бронштейн, И. Н. Справочник по математике для инженеров и учащихся втузов Текст. / И. Н. Бронштейн, К. А. Семендяев. М.: Наука, 1986. -544 с.

61. Будак, Б. М. Кратные интегралы и ряды Текст. / Б. М. Будак, С. В. Фомин. 2-е изд. - М.: Наука, 1967. - 608 с.

62. Крылов, В. И. Приближенное вычисление интегралов Текст. / В. И. Крылов. 2-е изд., доп. - М.: Наука, 1967. - 500 с.

63. Абрамович, М. Справочник по специальным функциям с формулами, графиками и математическими таблицами Текст. / М. Абрамович, И. Стиган; пер. с англ. М.: Наука, 1979. - 832 с.

64. Хейгеман, JT. Прикладные итерационные методы Текст. / J1. Хейгеман, Д. Янг; пер. с англ. М.: Мир, 1986. - 448 с.

65. Кужелев, С. В. Алгоритм вычисления частных производных первого и второго порядков от геопотенциала по координатам ИСЗ Текст. / С.В. Кужелев // Изв.* вузов. Геодезия и аэрофотосъемка. 1982. - №5. - С. 79-86.

66. Кужелев, С. В. Разработка и исследование алгоритмов построения и оценки точности спутниковых траекторий в орбитальном методе космической геодезии Текст. : Автореф. дис. . канд. техн. наук / С.В. Кужелев. Новосибирск: НИИГАиК, 1983. - 24 с.

67. Программа прогнозирования движения геодезических искусственных спутников Земли Текст. / Ю. В. Сурнин, С. В. Кужелев, В. А. Ащеулов, Ю. В. Дементьев // Наблюдения искусственных спутников Земли. София, 1978.-№16.-С. 154-174.

68. Сурнин, 10. В. Модели движения ИСЗ и точность численного прогнозирования орбит Текст. / Ю.В. Сурнин, С.В. Кужелев // Геодезия и картография.-1982.-№10.-С. 8-13.

69. McCarthy, D.D. IERS Standards Text. / D.D. McCarthy (ed.) // IERS Technical Note 13, Observatoire de Paris, July, 1992.

70. Справочное руководство по небесной механике и астродинамике Текст. / В. К. Абалакин, Е. П. Аксенов, Е. А. Гребенников и др. 2-е изд., доп. и перераб. - М.: Наука, 1976. - 864 с.

71. Урмаев, М. С. Орбитальные методы космической геодезии Текст. / М. С. Урмаев. М.: Недра, 1981. - 256 с.

72. Everhart, Е. Implicit sincle sequcuce methods for integretingorbits Text. / E. Everhart // Celestial Mechanics, 1974,10. - P. 35-56.

73. Пайнз, С. Представление потенциала поля тяготения и его производных, свободное от особых точек Текст. / С. Пайнз // Ракетная техника и космонавтика. 1973. - Т.2, №Ц. С.65-69.

74. Эльясберг, П. Е. Определение движения по результатам измерений Текст. / П. Е. Эльясберг. -М.: Наука, 1976. 416 с.

75. Жданюк, Б. Ф. Основы статистической обработки измерений Текст. / Б. Ф. Жданюк. М.: Сов. радио, 1978. - 384 с.

76. Чарный, В. И. Об изохронных производных Текст. / В. И. Чарный // Искусственные спутники Земли. М.: Изд-во АН СССР, 1963. - Вып. 16. -С. 226-237.

77. Суханов, А. А. Об изохронных производных невозмущённого движения Текст. / А. А. Суханов. М.: ИКИ АН СССР, Препринт, 1989.

78. Телеганов, Н. А. Высшая геодезия и основы координатно-временных систем: учеб. пособие Текст. / Н. А. Телеганов, А. В. Елагин. Новосибирск: СГГА, 2004.-238 с.

79. Медведев, П. П. Глобальные космические навигационные системы (геодезическое использование) Текст. / П. П. Медведев, И. С. Баранов // Итоги науки и техники. Сер. Геодезия и аэросъёмка. М: ВИНИТИ, 1988. - Т. 29. -160 с.

80. Нейман, Ю. М. Вариационный метод физической геодезии Текст. / Ю. М. Нейман. М.: Недра, 1979. - 200 с.

81. Торге, В. Гравиметрия Текст. / В. Торге; пер. с англ. М., Мир, 1999.-429 с.

82. Гиенко, Е. Г. Определение уклонения отвесной линии и астрономических координат по наземным и GPS-измерениям Текст. / Е. Г. Гиенко, А. В. Елагин // Вестник Сибирской государственной геодезической академии. -Новосибирск, 2000.-Вып. 5.-С. 16-19.

83. Елагин, А. В. Применение метода конечных элементов для решения краевой задачи теории потенциала внутри нивелирного полигона Текст. / А.В. Елагин // Вестник Сибирской государственной геодезической академии. -Новосибирск, 2005.-Вып. 10.-С. 25-30.

84. Елагин, А. В. Методика определения конечно-элементной модели гравитационного поля Земли по спутниковым и наземным измерениям Текст. / А. В. Елагин // Вестник Сибирской государственной геодезической академии. -Новосибирск, 2006. Вып. 11. - С.46-52.

85. Zhong, D. Robust estimation and optimal selection of polynomial parameters for the interpolation of GPS geoid heights Text. / D. Zhong // Journal of Geodesy / Springer-Verlag. 1997, №71. - P.552-561.

Обратите внимание, представленные выше научные тексты размещены для ознакомления и получены посредством распознавания оригинальных текстов диссертаций (OCR). В связи с чем, в них могут содержаться ошибки, связанные с несовершенством алгоритмов распознавания. В PDF файлах диссертаций и авторефератов, которые мы доставляем, подобных ошибок нет.