Геодезические методы анализа высотных и плановых деформаций инженерных сооружений (стр. 4 из 4)
POINT 1 x = 6180000.000, m(x) = 0.000 ; y = 7400000.000, m(y) = 0.000
POINT 2 x = 6180866.261, m(x) = 0.002 ; y = 7399999.998, m(y) = 0.001
POINT 3 x = 6180866.257, m(x) = 0.002 ; y = 7400999.999, m(y) = 0.001
POINT 4 x = 6180866.258, m(x) = 0.003 ; y = 7401999.999, m(y) = 0.002
POINT 5 x = 6180000.000, m(x) = 0.003 ; y = 7401999.998, m(y) = 0.002
POINT 6 x = 6179999.996 , m(x) = 0.002 ; y = 7400999.997, m(y) = 0.001.
Во втором цикле в программу вводится общее имя файла, после чего указывается номер цикла (S = 2). Как только программа признает S>1,то она спрашивает, изменились ли координаты предыдущего цикла, с целью выявления возможных деформаций относительно предыдущего цикла. В данном случае были смоделированы деформации для точек 3 и 5, равные соответственно: - 0.015, - 0.010; - 0.017, - 0.012.
Далее, в программу вводятся деформации для выбранных точек, после чего вводятся те же ковариационные матрицы и корреляционные коэффициенты, что и в I цикле. Выполняется уравнивание рекуррентным способом для контроля грубых ошибок. Затем выполняется уравнивание параметрическим способом. Деформации получают на плоскости на эллипсоиде (WGS-84) в проекции Гаусса – Крюгера. Программа дает имена всех файлов, где хранится информация о втором цикле. Следующим пунктом программы выполняется анализ обнаруженных деформаций:
1. POINT3 Dx = -0.0154 Mx = 0.003; Dy = -0.0099 My = 0.002
2. POINT 5 Dx = -0.0162 Mx = 0.002 ;Dy = -0.0117 My = 0.002.
Обнаруженные программой деформации по величине близкие к модели.
Координаты, полученные в результате обработки данных второго цикла, отличаются от первого в тех точках, которые подвергались деформациям (в точках 3 и 5):
POINT 3 x = 6180866.245, m(x) = 0.003 ; y = 7400999.990, m(y) = 0.002
POINT 5 x = 6179999.986, m(x) = 0.004 ; y = 7401999.988, m(y) = 0.003.
При обработке данных третьего цикла, задача ставится так, чтобы те точки, которые подвергались деформациям во втором цикле, также должны деформироваться. Таким образом, проверим, выполняется ли объединение циклов с целью повышения точности определения плановых деформаций, что является одним из преимуществ данного алгоритма.
Для точек (3 и 5) смоделированы деформации равные:
0.0117, -0.0137; -0.0099, -0.0099.
В результате после уравнивания получены следующие координаты:
POINT 3 x = 6180866.233, m(x) = 0.003; y = 7400999.975, m(y) = 0.002
POINT 5 x = 6179999.975, m(x) = 0.004; y = 7401999.975, m(y) = 0.003.
Следует отметить, что координаты неподвижных точек не отличаются от предыдущих циклов.
Деформации реперов в третьем цикле относительно второго цикла:
POINT 3 Dx = - 0.0116 Mx = 0.002; Dy = - 0.0139 My = 0.001
POINT 5 Dx = - 0.0097 Mx = 0.002; Dy = - 0.0099 My = 0.002.
Проверку объединения циклов можно сделать с помощью информации о векторе суммарных деформаций, так как он представляет деформации относительно и первого цикла, и второго цикла.
POINT 3 CYKLE 1- 2 Dx =- 0.014, Dy = -0.011; CYKLE 2- 3 Dx =- 0.012, Dy =- 0.014
ВЕКТОР СУММАРНЫХ ДЕФОРМАЦИЙ:Dx =- 0.026, Dy =-0.025
POINT 5 CYKLE 1- 2 Dx =- 0.017, Dy =- 0.011; CYKLE 2- 3 Dx =-0.010, Dy =-0.010
ВЕКТОР СУММАРНЫХ ДЕФОРМАЦИЙ: Dx = -.028, Dy = -.021
На рисунке № 5 показаны не только деформации, которым подвергается точка в текущем цикле, но и деформации предыдущего цикла.
Рис 5
Достоверность работоспособности программы для анализа плановых деформаций по GPS- измерениям.
С целью проверки работоспособности составленной программы, в четвертой главе выполнен детальный анализ полученных результатов. Обработаны четыре цикла тех же данных проекта, но «традиционным методом». Здесь под «традиционным методом» понимается метод, когда деформации определенных точек определяют по формуле:
; 
, (9)где
– горизонтальное перемещение ( определяемая деформация), 
и 
– координаты второго и первого циклов.Средние квадратические ошибки определения этих координат находят следующим образом:
и 
, (10)где
– средняя квадратическая ошибка определения деформации; 
; 
– СКО определения деформаций текущего и предыдущего циклов.Тогда:
и 
. (11)Результаты обработки данных «традиционным методом» ІІ цикла:
Табл. 9
| № | Имя марки | Х уΙΙ цикла | СКОх,у(м) | Х уΙ цикла | СКОх,у(м) | D х,у(м) | СКО(D)(м) |
| 1 | moz-1 | 6180000.0007400000.000 | 0.0000.000 | 6180000.0007400000.000 | 0.0000.000 | 0.0000.000 | 0.0000.000 |
| 2 | moz-2 | 6180866.2657400000.003 | 0.0030.002 | 6180866.2617399999.998 | 0.0020.001 | 0.0040.005 | 0.0040.002 |
| 3 | moz-3 | 6180866.2487400999.987 | 0.0030.002 | 6180866.2577400999.999 | 0.0020.001 | -0.009-0.012 | 0.0040.002 |
| 4 | moz-4 | 6180866.2607401999.997 | 0.0040.003 | 6180866.2587401999.999 | 0.0030.002 | 0.002-0.002 | 0.0050.004 |
| 5 | moz-5 | 6179999.9817401999.982 | 0.0040.003 | 6180000.0007401999.998 | 0.0030.002 | -0.019-0.016 | 0.0050.004 |
| 6 | moz-6 | 6180000.0037401000.002 | 0.0030.002 | 6179999.9967400999.997 | 0.0020.001 | 0.0070.005 | 0.0040.002 |
Таблица № 9 в столбцах 3 и 4 содержит уравненные координаты реперов и СКО их определения во втором цикле. Аналогичная информация первого цикла содержится в столбцах 5 и 6. Последние два столбца заполнены деформациями второго цикла и СКО определения этих деформаций.
Далее, составлены несколько таблиц для сравнения вычисленных разными методами деформаций с теми величинами, которые вводились при их моделировании. Проанализируем результаты обработки данных второго цикла:
Табл. 11
| 1 | 2 | 3 | 4 | 5 |
| № марки | Имя марки | D(х,у)(мм)программа | D(х,у) (мм)смоделированные | Разности столбцов 3-4 |
| 1 | moz-3 | Х = -0.0154У =- 0.0099 | Х = -0.015У = -0.0100 | Х = -0.0004У = 0.0001 |
| 2 | moz-5 | Х = -0.0162У = -0.0117 | Х = -0.0170У = - 0.0120 | Х = 0.0008У = 0.0003 |
В этой таблице в столбце 3 отражены деформации, полученные при использовании программы. В четвертом столбце величины деформаций, которые вводились для деформирования точек, имена которых содержатся в столбце 2. Последний столбец заполнен разностью между смоделированными и вычисленными программой деформациями.
В таблице № 12 отражена информация аналогичная той, которая содержится в таблице № 11, но обработанная «традиционным методом».
Табл. 12
| 1 | 2 | 3 | 4 | 5 |
| № марки | Имя марки | D(х,у) мм)Трад. метод | D(х,у) (мм)Смоделирован-ные | Разности столбцов 3-4 |
| 1 | moz-3 | Х = -0.009У =- 0.012 | Х = -0.015У = -0.010 | Х = +0.006У = -0.002 |
| 2 | moz-5 | Х = -0.019У = -0.016 | Х = -0.017У = - 0.012 | Х = -0.002У = -0.004 |
Анализ этих двух таблиц показывает, что вычисленные деформации с использованием составленной программы ближе по величине к смоделированным. При обработке тех же данных «традиционным методом» разность между смоделированными и вычисленными деформациями значительно отличается от полученных по программе.
Таким образом, анализ данных по всем циклам даёт возможность с учётом их объединения получить желаемые результаты, учитывая все реально существующие деформации на данном объекте.
Заключение
1. До сих пор применяемый «традиционный метод» как вычисление разности между предыдущим и текущим циклами менее точный, чем вычисления с использованием алгоритма последовательного объединения циклов.
2. С целью экономии памяти компьютера удалось разделить процесс контроля грубых ошибок, учитывая только необходимые измерения, а для избыточных измерений вычисляются только их свободные члены уравнений поправок, и выполняется контроль грубых ошибок. При этом отпадает необходимость вычислять всю матрицу (5).
3.Объединение циклов выполняется параметрическим способом, при этом порядок матрицы коэффициентов нормальных уравнений определяется только числом пунктов, тогда как применение формулы (5) требует удвоенное число.
4.Экспериментальные исследования, выполненные в данной диссертации, направлены на применение современных геодезических методов для анализа высотных и плановых деформаций инженерных сооружений и земной поверхности.
5.Эти программы особо актуальны для республики Мозамбик, так как из-за периодических наводнений необходимо контролировать состояние всех сооружений, находящихся в зоне наводнений.
ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ
1. Ассане Антонио Алфредо. Анализ плановых деформаций по GPS – измерениям. Деп. в ОНИПР ЦНИИГАиК, 14.11.2006, № 886.
2. Ассане Антонио Алфредо. Методы определения деформации инженерных сооружений. Деп. в ОНИПР ЦНИИГАиК, 14.11.2006,№ 887
3. Ассане Антонио Алфредо, Ю.И. Маркузе, Е.П. Власенко. Анализ плановых деформаций инженерных сооружений и земной поверхности. Геодезия и картография, № 3 март 2007, с. 28-32.