Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
учебники / Клюшин Е.Б., Михелев Д.Ш., Барков Д.П. - Практикум по прикладной геодезии - 1993.pdf
Скачиваний:
379
Добавлен:
26.04.2015
Размер:
13.91 Mб
Скачать

любой заданной функции с помощью простых вычислительных средств. Однако при вычислении самой матрицы (120) необходимо применять ЭВМ. При этом теряется оперативность оценки проекта сети. Поэтому приближенные методы оценки проектов разбивочных сетей также используются в практике работ.

§ 18. Приближенные способы вычисления обратного веса функции при оценке проекта

Рассмотрим наиболее распространенные методы получения при­ ближенных формул для оценки проектов PC.

1.В матрице обратных весов координат (120) вычисляют только

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

2.Не учитывают малые по абсолютному значению коэффициенты исходной матрицы А или В и коэффициенты заданной весовой функции / (пренебрегают тем самым некоторыми «незначимыми»

избыточными измерениями, т. е.

упрощают

схему

сети).

3. При вычислении матрицы

(120) по

формуле

(118) или (119)

принимают матрицу весов измерений равной Е, т. е. пренебрегают

разновесностью

результатов измерений (считают их равноточными).

4. Матрицу

(120), полученную для типовой фигуры (т. е. для

сети правильной геометрической формы) и ориентированной по осям координат, используют для оценки реальной сети, произвольно

ориентированной и

отличающейся

по форме от типовой фигуры,

т. е. пренебрегают

геометрической

формой сети

(влияние ориен­

тирования

устраняется умножением

полученных

оценок координат

на матрицу

поворота).

 

 

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

6.Выделяют в сети типовое построение, включающее в себя все пункты сети и по известным для данного типового построения формулам оценивают заданную функцию Q'f . Все избыточные измере­

ния Qf учитывают из приближенного соотношения

е / = £ е > .

( н о

где п — число всех измерений в типовом построении; N — число всех измерений в сети (способ проф. А. И. Дурнева).

7. На основе численного анализа (с использованием строгих формул и применением ЭВМ) строят приближенные уравнения регрессии или полиномальные уравнения, связывающие обратный вес

68

оцениваемой функции с числом типовых фигур в сети или их продвигом, которые используют затем при оценке реальных постро­

ений аналогичного вида и формы.

 

 

Рассмотренные методы или их различные

комбинации

приводят

к формулам, которые могут использоваться

в полевых

условиях

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

Наиболее распространенный метод при приближенной оценке проектов PC — метод эквивалентной замены. Технология вычислений при этом заключается в следующем. Намечают ходовую линию Z от исходного пункта О до оцениваемого пункта сети К. Вычисляют обратный вес приращений координат пунктов ходовой линии по формулам

v2

Q ax :j = Qd cos2aij + Q, sin2a0-( у

2

(142)

 

QAY,. = Qdsin2ay + Qacos2aiy I —

 

где QD и Qa— обратные веса длин сторон ij хода Z и их дирекционных углов соответственно.

Вычисляют обратный вес координат конечной точки К ходовой линии Z как сумму обратных весов соответствующих приращений

координат предыдущих точек:

 

£б**,;

Q rK.z = I Сду,,

(143)

i=l

 

i=l

 

где г — число сторон

ходовой линии

Z.

Затем находят веса координат точки К по данной Z -й ходовой

линии:

 

 

 

Рхкл = 7,— ; р г

KZ

= —

р гк.г = — -(144)

fix*.,

бгх.,'

 

Таким же образом рассчитывают веса координат К по всем другим ходовым линиям и находят общий вес как сумму весов по отдельным ходам:

i

РХк,г' Р *к= i РуК.г

Z = 1

Z = 1

где п — число ходовых линий, сходящихся в точке Наконец, вычисляют обратный вес координат X

( 145>

К.

и Y точки К:

Заметим,

что

формулы

(143)— (146) справедливы при

оценке

функции

любого

вида.

Так,

например,

если

оценивается

длина

D стороны сети (или ее

направление а), то в этих формулах вместо

индексов

X и

У необходимо

использовать

D

или

ос.

 

Способом

эквивалентной замены целесообразно

оценивать сети

с небольшим числом избыточных измерений. Для сложных сетей лучше воспользоваться способом, изложенным в п. 6.

При строгом методе оценки ошибки исходных данных учитываются автоматически (структурой исходных матриц А или В, а также матриц Р и /). Если в этих матрицах часть коэффициентов, относящихся к исходным данным, выделить в отдельный блок, то в итоге расчетов по формулам (117), (118), (119) определятся две

матрицы

обратных

весов

оцениваемых

функций:

одна — зависящая

только от

ошибок

измерений, другая — зависящая только

от ошибок

исходных

данных.

Из их

суммы будет

получен

общий

обратный

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

§ 19. Оценка проекта триангуляции

Строгая оценка проекта сети триангуляции была рассмотрена ранее (см. §§ 16, 17). Применяя ее к малым PC правильной формы, можно получать сравнительно простые рабочие формулы для вычис­ ления обратных весов отдельных элементов (см. § 18). Так, например, исходя из коррелатного способа уравнивания по направлениям и приняв за единицу веса измеренное направление |х = тн, для PC триангуляции в виде геодезического четырехугольника прямоугольной формы справедливы следующие формулы [3]. Для продвига фигуры К можно записать

К= — =~. АВ Ъ

Если сторона А В — базис (рис. 25,а), то

(147)

Если исходное направление совпадает с базисом, то

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

Если

исходное направление совпадает со

стороной ВС, то

 

 

 

 

 

 

 

 

 

04 + 4?2 + 4

(149)

Q*AB

Q*DC

1,2;

Q . = 1,5;

Q.AC = Q.DB= 1,25-

 

 

 

 

 

 

 

 

 

i(q* +q2+\)'

 

]

 

 

 

 

 

 

 

 

 

где Ч=К

Ь— диагональ

АС

 

 

 

Если

базис

(рис. 25,6),

то

 

 

 

3

S \ D

 

ЪКл+2К2+2

 

 

QAD — QBC~ Z

р2

К 2 (КА+2К2+\)’

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

(150)

о = 0

° ьс( 1 25

4* 4+4* 2+|

за:2

 

 

 

Q D C

Q A B

„ 2

\

»

8 (л :4 + к 2 + 1 )

2(л:4+я 2+1)‘

Если

исходным направлением

будет сторона ВС, то

 

 

, .

^

 

^

 

Л

Л

4 Я 4+4Л:2+1

(151)

Q * A D ~ ^ 5 ;

Q < * d c ~ Q * a b ^

Q * a c ~ Q * b d -

8 ( t f 4 + J £ 2 + l ) ’

 

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

проекта). Отметим, что формулы

получены при \i = mH.

Если за единицу веса принять

не направление, а угол ц= т р,

то вычисленные по формулам (147)— (151)

обратные веса необходимо

уменьшить в два

раза, так как

[1 н = т р/2.

В

приведенных

формулах

величину S целесообразно выражать в

километрах,

а

р следует

принять равным

0,206.

 

 

Qs и

Qa,

 

Вычислив по

формулам

(147)— (151)

величины

ms и та находят по выражению (102) (считая исходные данные безошибочными):

при

|х = ти

 

 

 

 

 

ms = m„ S/ Q S\ тл= ти J~Q^

 

 

При

Ц= /Ир

 

 

 

 

 

ms ~-

=mt>'/QTs’ m°=mt>J ~ ^ = mt s / o i -

 

 

!=mp J ^

 

 

Если же необходимо учитывать ошибки исходных данных (базиса

и направления),

то (|х= тн )

 

 

m l = m l K 2 +m%Qs,

m l = m 2Kcx+ m b Q a.

(153)

 

Выполним

оценку

разбивочной сети

микротриангуляции

(см.

рис. 25, а), в которой

исходным базисом

служит сторона А В,

а ис­

ходным направлением — сторона ВС.

За единицу веса принята средняя

квадратическая ошибка

измерения

угла шр = 2".

Расчеты

проведем

по формулам (147), (149) и (152) при

К= 2; ( 7 = 0,5; р = 0,206 (табл. 16).

 

Из расчетов следует, что относительная погрешность в длинах

сторон

1: Т= ms : S, вызванная

только ошибками

угловых

измерений

и

геометрией

фигуры

(форма,

размер),

 

составляет

mDC: DC =3,5:259 000= 1:74 ООО. На

полученную

точность

не

влияют

ошибки

в измерении

базиса.

Поэтому базис

должен

быть

измерен

с относительной ошибкой не более

mb: b = 1:150 000,

что

в

данном

примере

соответствует

абсолютной

ошибке

mb= 259 000: 150 000 =

= 1,5 мм. Если базис измерен грубее, то необходимо воспользоваться

формулой (153). Пусть, например,

базис

измерен

светодальномером

с точностью

т ь = 5мм. Тогда по

формуле

(153)

для стороны DC

получим:

mDC =

• 22 + 3,52 = 10,6 мм

(1:24 400).

 

 

Зная ошибки сторон и их направлений, можно оценить ошибки

взаимного положения

смежных пунктов

по

формуле

 

 

 

 

S 2

 

 

 

 

 

 

 

(154)

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т аб л и ц а

16.

Резул ьтаты

оценки точности

сторон

разбивочной сети

 

Обозначение

ВС

AD

 

DC

 

АВ

BD

АС

S,

км

 

0,518

0,518

 

0,259

0,259

0,579

0,579

Qs,

мм2/с2

24,9

24,9

 

6,32

 

24,3

24,3

 

Q* 7

7

 

1,5

 

1,2

 

1,2

0,77

0,77

Q 's,

мм

12,4

12.4

 

3,16

 

12,2

12,2

<2;

 

 

0,75

 

0,60

 

0,60

0,38

0,38

ms, мм

7

7

 

3,5

 

7

7

 

с

 

 

\ Л

 

1,5

 

1,5

1,2

1,2

Для данного примера (см. табл. 16) для пунктов D и С имеем (без учета ошибки измерения базиса):

 

0,259 х2

М = / 3,52Н-1,52( —— =4,0 мм.

V

\0’206/

Из формулы (154) следует, что обратный вес взаимного положения вершин треугольников триангуляции (соединенных измеренным на­ правлением)

QbX,AY = QDC + Qa •

(155)

Для пунктов D и С при [i = mH будем иметь

Qax.ay 6,32+ 1,2 • 1,58 = 8,22 мм2/с2.

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

при mf = mBC= 10 мм,

Ю

Wp= — = = 2 ,8 "; У*2А

при mf = m^AB = 5",

= —Т===6,4";

>Дбо

при mf =

ау = 5 мм,

т^ = — = 2,5 .

V 4

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

Q = £

 

Qi-

 

(156)

i= 1

 

 

 

Определим обратный вес стороны EF сети (рис. 26).

Ранее

 

было получено QBC = 24,9 мм2/с2. Считая для стороны EF

исходным

 

базисом

сторону ВС (АГ=0,5), по формуле (147) получаем

С ледовательно,

общ ий вес длины

стороны E F составит

Q EF = Q BC + Q EF = 24,9 +1,58 = 26,5.

 

Если

задана ош ибка w£F=10mm,

то по формуле (103) найдем

73

Рис. 26. Разбивочная сеть три­ ангуляции в виде двух смежных геодезических четырехугольни­ ков

Рис. 27. Ряд триангуляции

М

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

e i = 3,54S2X /?i; Q'a = - K ,

 

 

 

 

(157)

где S

и

а — соответственно длина

(в км)

и

направление

связующей

стороны

К треугольника ряда.

 

 

 

 

 

Формулы (157) получены для случая уравнивания ряда по углам,

т. е. при

|!= тр.

 

 

 

 

 

 

Величина /?, для каждого i треугольника ряда выбирается из

таблиц

по

значениям

связующих

углов

А

и В.

По

значениям

Qs и Qa для связующих сторон ряда

и

формулам

(152)— (155)

выполняется

его оценка.

 

 

 

 

 

С

учетом формул

(157) приближенная оценка

обратных весов

координат точек ряда триангуляции может быть выполнена по выражениям (142), (143).

При оценке проекта сети микротриангуляции (рис. 28), состоящей

А В В

рис. 28. Разбивочная сеть триангуляции из трех прямоугольных треугольников

из прямоугольных треугольников (с углами 30, 60°), в качестве весовой функции выбрано положение точки Е.

Обратные веса координат точки Е определим приближенно по формулам (142), (143) с учетом выражений (157). Для этого выберем

ходовую

линию

от

В

к Е — она идет по

связующим

сторонам ВС,

CD

и DE треугольников

ABC, BCD и DCE. (табл. 17).

 

По данным табл. 17 составляем весовую диагональную матрицу

(121)

координат

точки

Е

для

данной сети (см. рис. 28):

 

е * . . = ( 31

 

23)

 

 

 

 

 

 

 

 

Т абл и ц а

17. Обратные

веса элементов триангуляционной сети

 

Связую­

 

 

 

 

Номер

 

а,

Q*

 

 

 

щая

сто­

Связующий

 

 

тре­

S , км

Qa

QAX,

QAY*

рона

хо­

угол,

градус

Ri

уголь­

градус м м 22

мм 22

мм22

довой

 

 

 

 

 

ника

 

 

 

 

 

 

линии

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

3

4

5

6

7

8

9

10

вс

 

30,

60

 

19

1

0,4

30

10,8

0,67

8,73

4,59

C D

 

60, 90

 

1

2

0,5

180

17,7

1,34

17,70

7,89

I

 

 

 

20

26,43

12,48

D E

 

90, 30

 

13

3

0,3

90

10,51

2,00

4,24

10,51

I

 

 

 

33

30,67

22,99

1I

Примечание. Данные граф 1, 2, 4, 5, 6 выбирают со схемы сети. Значение Ri (графа 3) определяют с учетом соответствующих связующих углов в единицах 6-го знака логарифмов.

= = 30,67*31; QY^ Q A, ^ 22,99*23,0; 6 *= 31.0+23,0 = 54,0.

без учета

корреляции координат X

и У. И обратный

вес положения

Е (точечная оценка) будет равен

Q E = 54 м м 2/ с 2.

 

Если

|л= шр = 2", то

прямая оценка точности по

формуле (111)

приводит

к следующим

результатам:

 

«1 *£= 2"ч/зТ = П мм;

шу1.= 2"%/23=10мм; М Е= 2 "^5 4 = 15 мм,

а обратная оценка точности по выражению (103) при заданном значении ошибок функций тх = т у= М Е= 10 мм — к следующим значе­ ниям ошибок единицы веса т р:

т>‘ =Ж 2 У - " ь - ж ' *

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

Допустим, что по условиям объекта и строительства разбивочную сеть (см. рис. 28) стало возможно трансформировать (рис. 29). Тре­

угольник ВСЕ— равносторонний

(5=400 м), а ходовая линия от

В к Е включает в себя стороны

ВС и BE. Расчеты обратных весов

координат точки Е для данного случая аналогичны предыдущим

(табл. 18).

 

 

 

 

Следовательно, при прямой оценке точности (гар = 2") по

формуле

(111) получаем

 

 

 

 

w*£= 2"v/!4 = 7,5 мм; Wy£= 2"4/l8 = 8,5 мм;

М Е= 2"^/зТ = 11 мм,

а при обратной

оценке точности (тХЕ = туЕ = М Е = \ 0 мм) по

выраже­

нию (103) находим

 

 

 

Ю ~

Ю -

Ю

t

 

т° * ~ Ж

’ т,' ~ ж

■ ■ тн' Ж

'

 

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

Если в типовом ряде триангуляции имеется два и более базиса (или исходного направления), то оценку обратного веса заданной стороны (направления, координат пункта) можно выполнить методом эквивалентной замены по формулам (144)— (146).

Рассмотрим сеть (рис. 30) с двумя измеренными базисами и со­ стоящую из двух однотипных геодезических четырехугольников прямо­ угольной формы. Пусть в качестве весовой функции задана длина стороны ВС.

Рис. 29. Разбивочная сеть триангуляции из прямоугольного и равнобедренного треугольников

Рис. 30. Разбивочная сеть триангуляции со сдвоенными геодезическими прямо­ угольниками

Т а б л и ц а

18. Обратные

веса элементов

трансформированной триангуляционной

сети

 

 

 

 

 

 

 

 

 

 

 

Связую­

 

 

 

Номер

 

 

 

 

 

 

щая

сто­

Связую­

 

тре­

 

а,

Qs,

 

 

<2дь

рона

хо­

щий угол,

Ri

уголь­

5, км градус мм2/с2

Q ;

мм2/с2

мм2/с2

довой

градус

 

ника

 

 

 

 

 

 

линии

 

 

 

 

 

 

 

 

 

 

В С

 

30,

60

19

1

0,4

30

10,8

0,67

8,73

4,59

B E

 

60,

60

4

2

0,4

90

13,0

1,34

5,00

13,00

I

 

 

23

 

 

 

 

 

13,78

17,59

 

Qxt — 13,8 мм2/с2;

QYk= 17,6 м м 2/ с 2; Qe = 13,8+ 17,6 = 31 мм2/с2.

 

 

Определим обратный вес длины стороны ВС относительно базиса А В из первого прямоугольника. Ранее было получено Q вс = 24,9. Найдем обратный вес этой же стороны ВС относительно базиса СЕ из второго прямоугольника. Поскольку по условию задачи оба прямоу­ гольника одинаковы и имеют один и тот же продвиг К= 2, то, очевидно, получим тот же результат, что и в первом случае, т. е. Р"ъс = 24,9. Теперь по формулам (144)— (146) последовательно вычисляем

/ > ', с

=2^9 =0’°4; П с =2^9 =°’04; Р *с = Р'вс + Р'кс =0,08;

QHC= — =^4s= 12,5 мм2/с2.

 

 

Рве

0,08

1

 

 

При

заданном т$ = 2" по

формуле (102)

получаем

твс = 2 "

Г\2^5

 

 

 

——

=5 мм,

 

 

а при

заданной

швс= Ю м м

по выражению

(103) находим

Сравнив эти результаты с ранее полученными значениями (твс = 7 мм, т р = 2,8"), можно сделать вывод, что измерение второго базиса улучшило

результаты, как и следовало ожидать, в y / l раза. Следует заметить, что в такой же сети, но без второго измеренного базиса СЕ, дополнительные измерения, выполненные во втором четырехугольнике BCEF (висячий), практически не влияют на точность определения стороны ВС.

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

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

1. Подсчитывают число N всех измеренных направлений в данной (сложной) сети триангуляции, включающую г пунктов.

2.Выделяют простую цепочку (ряд) треугольников, включающую

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

3.Определяют обратный вес Q' оцениваемого элемента из простой

цепочки

треугольников по известным рабочим формулам.

4. По

формуле (150) вычисляют окончательный обратный вес

оцениваемого элемента.

Описанный способ позволяет получить удовлетворительные резуль­ таты при N н ^Ю , причем конечные (наиболее удаленные) пункты сети должны быть связаны с исходными непосредственно наблюден­ ными направлениями. Такие сложные триангуляционные сети созда­ ются при строительстве тоннелей и гидростанций.

Рассмотрим сложную сеть гидротехнической триангуляции (рис. 31) (Видуев И. Г., Ракитов Д. И. Приложение геодезии в инженерно-стро­ ительном деле — М., 1964), в которой для определения координат четырех пунктов измерено N= 30 направлений. Простая цепочка из четырех треугольников включает и=18 направлений. Пусть требуется оценить длину стороны EF= 300 м (ось плотины). Воспользовавшись формулами (147), найдем из простой цепочки двух треугольников (связующие углы 63е и 80°): £ Л = 2 + 2 = 4 и 05=1,3.

Затем по формуле (141) получим окончательное значение обратного веса Qs (с учетом всех остальных избыточных измерений, но без учета ошибок измерения базиса):

18

Qs = — • 1,3 = 0,78 мм2/с2.

Соседние файлы в папке учебники