Анализ пространственно-временного состояния объекта методами математического моделирования

Министерство  образования и науки Российской Федерации

Федеральное агентство по образованию

Государственное образовательное учреждение высшего  профессионального образования

«СИБИРСКАЯ  ГОСУДАРСТВЕННАЯ ГЕОДЕЗИЧЕСКАЯ  АКАДЕМИЯ»

(ГОУ  ВПО «СГГА») 

Институт  геодезии и менеджмента

Кафедра прикладной информатики 

Курсовой  проект 

   Анализ  пространственно - временного состояния объекта методами математического моделирования по геодезическим данным 
 
 
 
 
 
 

Выполнили:                                                                              Проверила:

ст.гр. КГ-41                                                                                Бугакова Т.Ю.

Крупин Д.С.                                                                               Оценка__________

Радевич М.И.

(Вариант 14) 
 
 

Новосибирск 2010

 

1. Введение 

2. Теоретическая  часть……………………………………………………………………….…...

3. Практическая часть 

      3.1. Разработка  имитационной модели изменения пространственно-временного состояния объекта в трехмерном пространстве относительно неподвижной системы координат. 

      3.2. Построение концептуальной модели изменения пространственно-временного состояния объекта в трехмерном пространстве. 

      3.3 Построение модели  изменения состояния объекта в фазовом и Гильбертовом пространствах. 

      3.4 Оценка математической модели пространственно-временного состояния объекта. 

      3.5 Оценка и анализ результатов моделирования. 

      3.6 Статистический метод оценки изменения пространственно-временного состояния объекта. 

      3.7 Прогнозирование функции отклика объекта на изменение его геометрических свойств. 

      4. Заключение..........................................................................................

5. Список литературы................... 

      Приложение 1………………………………………………………………………………… 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

              1.Введение

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

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

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

   Задачами курсовой работы является:

      1. Разработка  имитационной модели изменения пространственно-временного состояния объекта в трехмерном пространстве относительно неподвижной системы координат

      2. Построение концептуальной модели изменения пространственно-временного состояния объекта в трехмерном пространстве

      3. Построение модели  изменения состояния объекта в фазовом и Гильбертовом пространствах

      4.Оценка математической модели пространственно-временного состояния объекта

      5. Оценка и анализ результатов моделирования.

      6. Статистический метод оценки изменения пространственно-временного состояния объекта

      7. Прогнозирование функции отклика объекта на изменение его геометрических свойств

3. ПРАКТИЧЕСКАЯ ЧАСТЬ

3.1. Разработка  имитационной  модели изменения  пространственно-временного  состояния объекта  в трехмерном пространстве  относительно неподвижной  системы координат.

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

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

При имитационном моделировании на компьютере выделяют следующие основные этапы:

  1. формулировка проблемы;
  2. построение математической модели системы;
  3. составление и отладка программы;
  4. планирование имитационных экспериментов;
  5. проведение экспериментов и обработка результатов.

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

Прямой  путь решения данной задачи весьма прост. Интервал времени [0;T], в течение  которого рассматривается работа системы, разбивается на интервалы длиной дельтаt, из-за чего данный способ решения  получил название принципа Δt. В пределах каждого интервала последовательно  вычисляются приращения всех процессов  в модели, и производится, если это  необходимо, изменение состояния  отдельных элементов модели. При  достаточно малых Δt получаем хорошее  приближение имитируемых процессов  к процессам в реальной системе  с параллельным выполнением операций. Очевидно, что точность моделирования  при этом достигается ценой больших  затрат времени. Принцип Δt является наиболее универсальным принципом  построения моделирующих алгоритмов, хотя и наименее экономичным с  точки зрения вычислительных ресурсов. Чаще всего он применяется для  моделирования непрерывных динамических систем. Однако данный способ мало пригоден для имитации систем, динамика которых  состоит в переходе из состояния  в состояние, причём в промежутках  между переходами состояние системы  остаётся неизменным. Каждый такой  переход связан с наступлением некоторого события в системе, например, приход входного или управляющего дискретного  сигнала, отказ элемента, достижение некоторой характеристикой системы  заданного порогового значения и  другие. Для такого класса систем наиболее часто используется принцип особых состояний. При построении алгоритма  имитации в соответствии с данным принципом функционирование системы, формализованное в математической модели, рассматривается как совокупность параллельно протекающих процессов, причём каждый процесс есть некоторая  последовательность событий, с каждым процессом связано изменение  состояния системы. Событие, возникающее  в системе, определяется как особое состояние. Процессы в общем случае не являются независимыми, а взаимодействуют  между собой.

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

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

Необходимо  задать имитацию случайного движения, при условии, что каждая координата в период времени от 0 до 10,14 изменяется случайным образом в пределах [0;1] при h = 0,5  и получить плановые координаты марок X(м) и Y(м) (Таблицы 2 и 3).

  Объект (Рисунок 1) представляет жилой комплекс, состоящий из пяти блоков не жестко связанных друг с другом. Грунты основания образованы песчаными суглинками и глиной. Наличие просадочных грунтов и дополнительных нагрузок могло привести к деформированию основания здания, наблюдения за объектом выполнено нивелированием 2 класса

                                     У(м)   
 
 
 
 
 

                                       1                             15          14        13

                                                                                                                          12      11 
 

                                                       А                                   В                           С  

                                                                                                                                            10

                                       2                                                                                        

                    

                                                                                                                 7      8           9

                                                             3               4           5              6

Рn.(Н=23.4561) 
 
 
 

                                                                                                                                             Х(м) 

         

      Н(м)                                                                                                                                         

                   

       Рисунок 1. Схема расположения геодезических марок  фундамента             
 
 
 
 
 
 

                                                                                                                               

Таблица 1. Значение высотных отметок марок

Дата

(мес)

Отметки высот марок (м)
Номера  марок
1 2 3 4 5 6 7
0 141,3675 141,3555 141,3517 141,3584 141,3626 141,3587 141,3517
0,14 141,3686 141,3514 141,3556 141,3515 141,3654 141,3566 141,3526
1,02 141,3665 141,3553 141,3595 141,3556 141,3627 141,3515 141,3527
1,16 141,3684 141,3567 141,3554 141,3520 141,3668 141,3504 141,3556
2,09 141,3653 141,3598 141,3515 141,3565 141,3679 141,3563 141,3517
3,12 141,3614 141,3559 141,3534 141,3586 141,3663 141,3525 141,3595
4,17 141,3607 141,3580 141,3595 141,3577 141,3622 141,3556 141,3541
5 141,3656 141,3551 141,3526 141,3524 141,3613 141,3577 141,3518
6,1 141,3624 141,3522 141,3527 141,3505 141,3604 141,3598 141,3589
6,28 141,3697 141,3513 141,3566 141,3566 141,3625 141,3587 141,3595
7,18 141,3658 141,3505 141,3514 141,3557 141,3696 141,3516 141,3584
8,14 141,3688 141,3596 141,3593 141,3558 141,3707 141,3557 141,3535
10,21 141,3688 141,3597 141,3594 141,3552 141,3626 141,3518 141,3576
12 141,3668 141,3588 141,3585 141,3555 141,3684 141,3509 141,3577
 

Продолжение таблицы 1 

Отметки высот марок (м)
Номера  марок
8 9 10 11 12 13 14 15
141,3525 141,1225 141,351 141,331 141,123 141,3564 141,352 141,2194
141,3526 141,1254 141,355 141,336 141,126 141,3565 141,352 141,2175
141,3517 141,1265 141,36 141,339 141,123 141,3516 141,353 141,2126
141,3528 141,1284 141,356 141,333 141,127 141,3507 141,352 141,2114
141,3557 141,1295 141,359 141,336 141,123 141,3568 141,352 141,2125
141,3586 141,1226 141,355 141,333 141,127 141,3526 141,36 141,2156
141,3527 141,1227 141,351 141,332 141,126 141,3555 141,359 141,2165
141,3568 141,1238 141,351 141,333 141,122 141,3574 141,355 141,2183
141,3559 141,1257 141,352 141,33 141,121 141,3596 141,356 141,2115
141,3563 141,1216 141,356 141,332 141,123 141,3587 141,359 141,2126
141,3594 141,1205 141,351 141,338 141,121 141,3518 141,355 141,2165
141,3505 141,1299 141,359 141,335 141,124 141,3529 141,351 141,2112
141,3573 141,1290 141,359 141,335 141,123 141,351 141,359 141,2155
141,3582 141,1280 141,358 141,336 141,128 141,3502 141,354 141,2166
 

Ошибка  измерений  Т=1мм,

Допустимая  относительная разность осадок ±0,5мм, 
 
 
 
 
 
 
 
 

Таблица 2.  Значение плановой координаты Х марок

Дата

(мес)

Х
Номера  марок
1 2 3 4 5 6 7
0 25,0000 22,5000 5,0000 7,5000 25,0000 27,5000 35,0000
0,14 25,0300 22,5140 5,0260 7,5170 25,0350 27,5130 35,0490
1,02 25,0050 22,5180 5,0280 7,5190 25,0160 27,5130 35,0010
1,16 25,0120 22,5030 5,0490 7,5500 25,0340 27,5370 35,0200
2,09 25,0090 22,5250 5,0070 7,5330 25,0220 27,5250 35,0040
3,12 25,0270 22,5140 5,0000 7,5470 25,0040 27,5210 35,0080
4,17 25,0200 22,5230 5,0110 7,5500 25,0310 27,5380 35,0410
5 25,0460 22,5020 5,0050 7,5480 25,0080 27,5010 35,0210
6,1 25,0120 22,5270 5,0160 7,5210 25,0030 27,5170 35,0500
6,28 25,0170 22,5090 5,0100 7,5280 25,0120 27,5140 35,0440
7,18 25,0280 22,5220 5,0120 7,5380 25,0410 27,5330 35,0150
8,14 25,0400 22,5220 5,0100 7,5270 25,0110 27,5160 35,0180
10,21 25,0220 22,5300 5,0380 7,5160 25,0270 27,5350 35,0090
12 25,0500 22,5500 5,0500 7,5500 25,0500 27,5500 35,0500
 

Продолжение таблицы 2 

Х
Номера  марок
8 9 10 11 12 13 14 15
42,5000 42,5000 27,5000 47,5000 52,5000 60,0000 52,5000 47,5000
42,5410 42,5010 27,5060 47,5050 52,5190 60,0460 52,5020 47,5480
42,5110 42,5370 27,5020 47,5280 52,5430 60,0240 52,5430 47,5120
42,5190 42,5340 27,5420 47,5340 52,5180 60,0170 52,5350 47,5390
42,5340 42,5270 27,5370 47,5150 52,5450 60,0500 52,5440 47,5070
42,5270 42,5210 27,5390 47,5240 52,5050 60,0480 52,5300 47,5290
42,5250 42,5250 27,5030 47,5000 52,5500 60,0130 52,5490 47,5290
42,5150 42,5420 27,5450 47,5150 52,5010 60,0280 52,5390 47,5210
42,5440 42,5090 27,5070 47,5030 52,5500 60,0430 52,5020 47,5460
42,5400 42,5040 27,5410 47,5470 52,5220 60,0310 52,5060 47,5310
42,5260 42,5160 27,5210 47,5340 52,5350 60,0410 52,5090 47,5300
42,5010 42,5180 27,5140 47,5470 52,5030 60,0380 52,5200 47,5140
42,5490 42,5130 27,5500 47,5240 52,5200 60,0160 52,5050 47,5490
42,5500 42,5500 27,5500 47,5500 52,5500 60,0500 52,5500 47,5500
 
 
 
 
 
 
 
 
 
 
 
 
 

Таблица 3.  Значение плановой координаты Y марок

Дата

(мес)

Y
Номера  марок
1 2 3 4 5 6 7
0 7,5000 5,0000 7,5000 25,0000 22,5000 22,5000 25,0000
0,14 7,5310 5,0490 7,5420 25,0400 22,5430 22,5310 25,0050
1,02 7,5330 5,0230 7,5250 25,0030 22,5490 22,5030 25,0480
1,16 7,5230 5,0230 7,5040 25,0240 22,5190 22,5440 25,0100
2,09 7,5500 5,0030 7,5130 25,0040 22,5190 22,5200 25,0480
3,12 7,5020 5,0420 7,5460 25,0440 22,5150 22,5000 25,0160
4,17 7,5230 5,0360 7,5220 25,0340 22,5240 22,5470 25,0250
5 7,5150 5,0430 7,5200 25,0320 22,5160 22,5440 25,0230
6,1 7,5100 5,0260 7,5290 25,0170 22,5220 22,5370 25,0080
6,28 7,5250 5,0500 7,5170 25,0040 22,5270 22,5480 25,0500
7,18 7,5320 5,0120 7,5060 25,0300 22,5390 22,5500 25,0010
8,14 7,5360 5,0290 7,5160 25,0340 22,5110 22,5280 25,0150
10,21 7,5320 5,0310 7,5220 25,0090 22,5350 22,5350 25,0170
12 7,5500 5,0500 7,5500 25,0500 22,5500 22,5500 25,0500
Анализ пространственно-временного состояния объекта методами математического моделирования