Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
geokniga-геодезия-попов-вн-чекалин-ви-2007.pdf
Скачиваний:
3640
Добавлен:
16.03.2015
Размер:
39.37 Mб
Скачать

Если условные уравнения изначально существенно нелинейны и при разложении в ряд Тейлора, вообще говоря, недостаточно ограничиваться первыми членами разложения, то условия (16.9) могут не выполниться. В этом случае производят второе приближение уравнивания, считая уравненные из первого приближения значения xi' измеренными, а свободными членами Wj – остаточные невязки в уравнениях (16.9).

В§ 153 рассмотрены примеры уравнивания различных геодезических построений коррелатным способом.

§153. Примеры коррелатного способа уравнивания

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

153.1.Уравнивание углов в полигоне

Вполигоне, состоящем из четырех вершин (рис. 16.7), неравноточно

измерены горизонтальные углы: А = β1 , В = β2 , С = β3 , D = β4 (табл. 16.6). Выполнить уравнивание углов без учета измерения длин сторон.

Предварительно найдем веса pi и обратные веса qi, приняв sе = [s4i ] ≈ 645 м

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

pi = se2

 

s2 s2

 

=

1

 

 

 

1 2

 

 

 

,

s

2

+ s

2

q

i

 

 

2

 

 

 

1

 

 

 

 

где s1 и s2 – стороны, образующие данный угол.

Шаг 1. Общее число измеренных величин n = измерений k = 3, число избыточных измерений r = 1.

(16.93)

4, число необходимых

415

Шаг 2. Составим условное уравнение (условие сумм углов полигона).

β1' + β2' + β3' + β4' – 3600 = 0

Всего одно уравнение, поскольку r = 1.

Шаг 3. Приводим условное уравнение к линейному виду, для чего продифференцируем его и найдем частные производные функции по аргументам βi . Очевидно, что

а11 = +1, а21 = +1, а31 = +1, а41 = +1.

Составим матрицу коэффициентов aij со строкой обратных весов qi (таблица 16.7).

 

 

Рис. 16.7. Уравнивание углов в полигоне

 

 

 

 

 

 

Таблица 16.6

Обозначени

Значение угла

 

Вес pi

Обратный вес qi

е

 

 

 

 

 

β1

800

16' 44,3"

 

0,221

4,520

β2

910

45' 00,7"

 

0,459

2,181

β3

690

25' 56,8"

 

0,473

2,113

β4

1180 32' 25,2"

 

0,225

4,452

 

 

 

 

 

Таблица 16.7

i→

 

1

2

3

4

j↓

 

 

 

 

 

1

+ 1

+ 1

+ 1

+ 1

рi

0,221

0,459

0,473

0,225

qi

4,520

2,181

2,113

4,452

Свободный член уравнения

W1 = [βi] – 3600 = 3600 00' 07,0" – 3600 = + 7".

Шаг 4. Найдем коэффициенты bjj нормальных уравнений (в данном случае – уравнений коррелат):

т.е.

[qa1a1 ]k1 + W1 = 0

(16.94)

2 + q2 a21

2 + q3a31

2 + q4 a41

2 )k1 + W1 = 0

(16.95)

(q1a11

416

Для приведенного примера, с учетом значений aij и qi , 13,266 k1 + 7 = 0, откуда k1 = - 0,528.

Шаг 5. Составляем условное уравнение поправок

ν 1 + ν 2 + ν 3 + ν 4 + 7′′ = 0

(16.96)

и формулы для вычисления поправок (с вычислением их значений):

ν1 = q1a11k1 = 4,520 (+1) (0,528) = −2,387′′ = −2,4′′

ν2

= q2 a21k1

= 2,181 (+1) (0,528) = −1,152′′ = −1,2′′

ν3

= q3a31k1

= 2,113 (+1) (0,528)

= −1,116

′′

 

′′

 

 

= −1,1

 

ν4

= q4 a41k1

= 4,452 (+1) (0,528)

 

′′

= −2,4

′′

= −2,351

 

 

Контроль по формуле (15.96): условие выполнено! Отступление при округлениях значений поправок на 0,1" является допустимым.

Шаг 6. Вычисляем уравненные значения углов:

β1' = 80° 16' 44,3" – 2,4" = 80° 16' 41,9"

β2' = 91°

45'

00,7"

– 1,1"

= 91° 44' 59,6"

β3' = 69° 25' 56,8" – 1,2" = 69° 25' 55,6"

β4' = 118°

32'

25,2"

– 2,4"

= 181° 32' 22,8"

Контроль: подстановка уравненных значений углов в уравнение (16.93) – условие выполнено!

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

153.2. Уравнивание системы нивелирных ходов с несколькими узловыми точками

На местности пройдена система нивелирных ходов с четырьмя узловыми точками 1, 2, 3 и 4 (рис. 16.8). В результате измерений образовано 9 секций, превышения в которых по указанному направлению приведены непосредственно на схеме. Указаны также высоты исходных реперов Р10, Р20 и Р30. В табл. 16.8 приведены длины ходов в секциях и значения весов и обратных весов превышений в секциях, вычисленные по формулам:

pi =

sе

;qi =

si

=

1

,

(16.97)

 

 

 

 

si

sе

pi

 

где sе = [sni ] ≈ 2км,n = 9.

Требуется определить уравненные значения высот узловых точек.

Шаг 1. Общее число измерений n = 9, число необходимых измерений k = 4, число избыточных измерений r = 5.

Шаг 2. Составим r = 5 условных уравнений:

1.h1' + h3' – h2' = 0

2.h4' + h5' – h3' = 0

3.h6' + h7' + h5' = 0

(16.98)

4.h7' + h8' – h9' – (HP30 – HP20) = 0

5.h1' + h4' + h8' – (HP30 – HP10) = 0

Шаг 3. Приведем условные уравнения к линейному виду, продифференцировав их по аргументам hi . Получим коэффициенты aij условных уравнений поправок:

а11 = +1 ;

а31 = +1;

а21 = - 1 ;

 

а42 = +1 ;

а52 = +1 ;

а32 = - 1 ;

 

а63 = +1 ;

а73 = +1 ;

а53 = +1 ;

417

а74 = +1 ;

а84 = +1 ;

а94 = - 1 ;

а15 = +1 ;

а45 = +1 ;

а85 = +1 .

 

Рис. 16.8. Схема нивелирных ходов с четырьмя узловыми точками

 

 

 

 

 

 

 

 

 

Таблица 16.8

 

 

Превышение

Длина хода

Вес p пре-

Обратный

 

 

секции

 

h,

s

в секции,

вышения

вес q пре-

 

 

 

 

мм

 

км

 

вышения

 

 

1

 

+3586

 

0,84

2,38

 

0,42

 

 

2

 

+2841

 

1,36

1,47

 

0,68

 

 

3

 

-752

 

2,15

0,93

 

1,08

 

 

4

 

-1243

 

0,78

2,56

 

0,39

 

 

5

 

+509

 

2,63

0,76

 

1,32

 

 

6

 

+5338

 

2,05

0,98

 

1,03

 

 

7

 

-5863

 

3,02

0,66

 

1,51

 

 

8

 

+4639

 

3,44

0,58

 

1,72

 

 

9

 

-3024

 

2,38

0,84

 

1,19

 

 

 

 

 

se = 2 км

 

 

 

 

Составим матрицу коэффициентов aij со строкой обратных весов qi

(табл. 16.9).

 

 

 

 

 

 

Таблица 16.9

 

 

 

 

 

 

 

 

1

2

3

4

5

6

7

8

9

i

 

 

 

 

 

 

 

 

 

j↓

 

 

 

 

 

 

 

 

 

1

+1

-1

+1

0

0

0

0

0

0

2

0

0

-1

+1

+1

0

0

0

0

3

0

0

0

0

+1

+1

+1

0

0

4

0

0

0

0

0

0

+1

+1

-1

5

+1

0

0

+1

0

0

0

+1

0

qi

0,42

0,68

1,08

0,39

1,32

1,03

1,51

1,72

1,19

Вычислим свободные члены (в мм), подставив в уравнения (16.98)

418

измеренные значения hi в секциях:

W1 = h1 + h3 – h2 = 3586 – 752 – 2841 = - 7 мм W2 = h4 + h5 – h3 = -1243 + 509 – (-752) = +18 мм

W3 = h6 + h7 + h5 = 5338 – 5863 + 509 = - 16 мм

W4 = h7 + h8 – h9 – (HP30 – HP20) = - 5863 + 4639 – (-3024) – 1794 = +6 мм W5 = h1 + h4 + h8 – (HP30 – HP10) = 3586 – 1243 + 4639 – 6965 = +17 мм

Шаг 4. Найдем по формулам (16.90) коэффициенты bjj нормальных уравнений коррелат:

1.(q1a11

2 + q2 a21

2 + q3a31

2 )k1 +

q3a31a32 k2 +

q1a11a15 k5 + W1

=

 

0

 

 

 

 

 

2.q3a32 a31k1 + (q3a32 2 + q4 a42 2 + q5a52 2 )k2 + q5a52 a53 k3 + q4 a42 a45 k5 + W2 = 0

3.q5a53 a52 k2

+ (q5a53

2 +

q6 a63

2

+

q7 a73

2 )k3 +

q7 a73 a74 k4 + W3

 

=

0

 

 

 

(16.99)

4.q7 a74 a73 k3 + (q7 a74 2 + q8a84 2 + q9 a94 2 )k4 + q8 a84 a85 k5 + W4

= 0

 

 

 

 

5.q a a k

1

+ q

4

a

45

a

42

k

2

+ q

8

a

85

a

84

k

4

+ (q a

2 + q

4

a

2 + q

8

a

85

2 )k

5

+ W

5

= 0

1

15

11

 

 

 

 

 

 

 

 

 

 

1

15

 

 

45

 

 

 

 

 

После подстановки значений aij и qi в уравнения (16.99) получим исходные нормальные уравнения коррелат:

1.

2,18k1 – 1,08 k2 + 0,42 k5 – 7 = 0

 

2.

-1,08 k1

+ 2,79 k2

+ 1,32 k3

+ 0,39 k5

+ 18 = 0

3. 1,32 k2

+ 3,86 k3

+ 1,51 k4

– 16 = 0

(16.100)

4.1,51 k3 + 4,42 k4 + 1,72 k5 + 6 = 0

5.0,42 k1 + 0,39 k2 + 1,72 k4 + 2,53 k5 + 17 = 0

Из решения системы уравнений (16.100) одним из способов получим:

k1 = - 2,137; k2 = - 11,552; k3 = +9,606; k4 = -3,882; k5 = -1,945. Контроль вычисления коррелат выполняем подстановкой их значений в

исходные уравнения (16.100):

1.2,18 (-2,137) – 1,08 (-11,552) + 0,42 (-1,945) – 7 = + 0,001

2.-1,08(-2,137) + 2,79(-11,552) + 1,32(+9,606) + 0,39(-1,945) + 18 = -0,001

3.1,32(-11,552) + 3,86(+9,606) + 1,51(-3,882) – 16 = - 0,031

4.1,51(+9,606) + 4,42(-3,882) + 1,72(-1,945) + 6 = +0,001

5.0,42(-2,137) + 0,39(-11,552) + 1,72(-3,882) + 2,53(-1,945) + 17 = -0,001

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

Шаг 5. Составляем условные уравнения поправок vi , пользуясь формулами (16.88) и табл. 16.9:

1.ν 1 = q1a11 k1 + q1a15 k5

6.ν 6 = q6 a63 k3

 

 

 

 

 

2.

 

q

2

a

21

k

1

7.ν

7

= q

7

a

73

k

3

+ q

7

a

74

k

4

ν 2 =

 

 

 

 

 

 

 

 

 

 

(16.101)

= q3 a31 k1 + q3 a32 k2

8.ν 8 = q8 a84 k4 + q8 a85 k5

3.ν 3

4.ν 4

= q4 a42 k2 + q4 a45 k5

9.ν 9 = q9 a94 k4

 

 

 

 

 

419