Решение сфероидических треугольников

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

Федеральное государственное автономное образовательное учреждение

высшего профессионального образования

«Дальневосточный федеральный университет»

 

Инженерная школа

Кафедра геодезии, землеустройства и кадастра

 

 

РЕШЕНИЕ СФЕРОИДИЧЕСКИХ ТРЕУГОЛЬНИКОВ

 

 

КУРСОВАЯ РАБОТА

 

 

Регистрационный №  ________   ___________  _______________       подпись                      И.О.Фамилия « _____» ___________________ 2014 г.

Студент _____________________                   Руководитель  к.т.н., профессор   _______________ Ф.И.О.     Оценка  ______________________   ___________  _________________           подпись                            И.О.Фамилия «_____» _______________ 2014 г.

 

Владивосток

2014 г.

 

 

Содержание работы:

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

2. Описать последовательность решения сферических треугольников с применением теоремы Лежандра и по способу аддитаментов.

Решить треугольники своей сети по способу аддитаментов, а затем, используя вычисленные стороны, решить эти же треугольники как линейные с применением теоремы Лежандра.

 

Контрольные вопросы

1. Что такое сфероидический треугольник?
2. При каких размерах сторон сфероидические треугольники можно решить как сферические, если требуется определить элементы треугольника с точностью 1·10-6?
3. В чем отличие решения сфероидических и сферических  треугольников?
4. Что такое аддитамент стороны, и как он вычисляется?
5. Сформулировать теорему Лежандра и привести формулу перехода от угла сфероидического треугольника к плоскому при больших сторонах.
6. Как вычисляется сферический избыток ε при сторонах меньших и больших 90 км?

Каковы возможные теоретические пределы изменения ε?

 

 

 

 

 

Решение сфероидических треугольников

 

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

 

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

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

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

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

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

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

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

Решение этой задачи выполним с использованием отображения части поверхности эллипсоида на шар, радиус которого примем равным среднему радиусу кривизны эллипсоида

 

 в некоторой точке О (рис. 1),

Рис. 1

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

Приняв точку О за полюс системы полярных координат So и А, отобразим часть поверхности эллипсоида на шар таким образом, чтобы полярные координаты точки Q1' на шаре не изменялись.

Тогда, при таком способе изображения линейные искажения в точке Q1' в направлении Q1`o`(дуги большого круга) будут отсутствовать, а в перпендикулярном направлении Q1` Q2`(дуги малого круга) будут наибольшими.

Обозначая длины элементарных дуг Q1Q2 иQ1`Q2,. как это показано на рис 1, можно найти наибольшие относительные линейные искажения ΔS:S, как:

 (1)

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

Рис. 2

На шаре (рис. 2) приведенной длине дуги большого круга ( с полюсом в точке О') будет соответствовать радиус кривизны геодезической окружности (малого круга ). Поэтому, для шара, непосредственно из чертежа (рис. 2), можно написать

(2)

Для поверхности эллипсоида приведенная длина геодезической линииmэ не имеет такой простой геометрической интерпретации как для сферы, поэтому, полагая, что mэ есть функция длины геодезической линии So, можно написать:

Очевидно f(o) = m0 есть приведенная длина геодезической линии, вычисленная для точки 0 (рис. 1) и, тогда:

 

    (3)

Для получения производных приведенной длины геодезической линии по длине So можно воспользоваться формулой (2), в которой для поверхности эллипсоида следует радиус считать величиной переменной.

Дифференцируя выражение (2) по So последовательно, находим:

 

и т.д.

В этих формулах через "к" обозначена полная кривизна поверхности эллипсоида.

(4)

Приведенная длина геодезической линии и ее производные в формуле (3) должны вычисляться по аргументам точки 0, для которой So=0. Но при So=0 , m0 как функция расстояния Sо, очевидно, также должно быть равно нулю, а производные примут следующие значения:

Подставляя производные в формулу приведенной длины (3), находим

 (5)

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

Так, например, для плоскости К0=0 и, поэтому, приведенная длина для плоскости равна самой длине линии.

Для сферы Ко=1/Ro2, а производные полной кривизны будут равны нулю, отсюда для сферы имеем:

(6)

Если в формуле приведенной длины дуги большого круга (2) синус заменить рядом, то, с точностью до членов пятого порядка малости, получим формулу (6).

Для получения формулы приведенной длины геодезической линии поверхности эллипсоида вначале найдем производную полной кривизны:

Продифференцировав формулу полной кривизны (4) по широте, а затем умножив полученное равенство на выражение производной dB/dS, находим

Подставив производную К0', а также полную кривизну поверхности эллипсоида (4) в выражение (5), получаем окончательно формулу вычисления приведенной длины геодезической линии на поверхности эллипсоида

 

(7)

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

Подставляя в числитель формулы (1) выражения (6) и(7) , а в знаменателе с достаточной точностью можно ограничиться mэ~So, находим

Из этой формулы видно, что наибольшие линейные искажения будут при Во = 45° и Ао = 0°. Следовательно,

 

(8)

 

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

Если, ориентируясь на точность первоклассных геодезических построений, принять (ΔS/S) max< 1*10-8 , то по формуле (8) находим радиус геодезической окружности, равный 133км. А так как вписать в окружность радиуса 133км можно треугольник со сторонами порядка 250-270км то, следовательно, сфероидические треугольники со сторонами, не превышающими 270 км, можно решать как сферические, при этом относительные искажения их элементов не будут превышать 1*10-8. Радиус сферы, при решении таких треугольников, следует принимать равным среднему радиусу кривизны для центра тяжести сфероидического треугольника.

 

Решение сферических треугольников

 

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

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

Рис. 3

Выражая стороны сферического треугольника (рис.3) в частях радиуса сферы:

при заданных углах  А, В,D и стороне а, находим:

 

                       (9)

 

 

или

 (10)

Если в треугольнике известны все стороны, то на основании теоремы косинуса стороны, будем иметь:

(11)

или

 (12)

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

На практике решение треугольников непосредственно по формулам сферической тригонометрии удобно и оправдано в том случае, если это решение выполняется на ЭВМ. Если же оно ведется в ручную - не по программе на ЭВМ, а с использованием настольных средств вычислительной техники, то решение, непосредственно, по формулам сферической тригонометрии становится практически громоздким. Действительно, в этом случае приходится с большой степенью точности вычислять ряд вспомогательных величин (R, a/R, sin (a/R), sin (b/R)), которые в конечном итоге не нужны.

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

 

Способ аддитаментов

 

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

Рассмотрим теоретические основы этого способа.

Полагая, что стороны сферического треугольника - малые величины (S< 200 км),по сравнению с радиусом сферы, разложим синусы сторон в выражении (9) в ряд, ограничиваясь членами пятого порядка малости: