Файл: Математическое моделирование в экологии.doc

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 31.03.2024

Просмотров: 1077

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

СОДЕРЖАНИЕ

Конспект лекций

Владикавказ

Математическое моделирование элементов сложных экологических систем

Лекция 1. Введение в моделирование. Исторический экскурс.

1. Основы моделирования в экологии 1.1. Общие принципы построения моделей в экологии

Лекция 2.

2.1. Элементы моделирования

2.2. Этапы построения математической модели

1.4. Элементы теории подобия, применяемые в моделировании

Лекция 3

3.2. Экологические модели

3.2.1. Основы экологометрики

3.2.2. Выборочный метод в экологометрике.

Зависимость числа интервалов от объема выборки

Статистический ряд по интервалам

Лекция 4. Статистические оценки параметров распределения случайных величин по выборкам

4.4. Статистические оценки гипотез об экологических моделях

Определение вариантов выборок

Выборка из генеральной совокупности

Статистическая таблица

Лекция 5.

Результаты эксперимента

Статистическая таблица эксперимента

Пример преобразования членов уравнения регрессии

Вычисление данных для линеаризации уравнения регрессии

Нормальные уравнения мнк для некоторых функций

Статистическое оценивание уравнения регрессии и парной корреляции.

Обработка результатов наблюдений

Лекция 6.

Рекомендации по выбору вида функции

3.4. Динамические статистические модели

Посадка леса

Данные по объему сброса качественных сточных вод

Данные по объему сброса сточных вод за 5-летие

Пример расчета 5-летних средних

Условное обозначение времени

Расчетные значения для определения уравнения динамики

Ряд динамики для определения сезонных колебаний

Лекция 7. Многофакторные эколого-математические модели. Анализ влияния отдельных факторов в экологической модели.

Эксперименталъный материал исследования

Результаты проведенных опытов

8.1. Анализ влияния отдельных факторов в экологической модели.

Лекция 9. Методы оптимизации. Метод Лагранжа

Лекция 10. Метод линейного программирования.

Лекция 11. Функциональные модели.

Лекция 12. Модели процессов содержащие обыкновенные дифференциальные уравнения.

Численные ошибки использованных для вычисления данных

Лекция 13. Статистические модели динамики.

Лекция 14. Балансовые модели.

Лекция 15.

Лекция 16. Информационные технологии в экологии. Экологические информационные системы.

1 6.1. Экологические информационные системы

1. Какова область значения для числовых характеристик?

Лекция 17. Использование информационных технологий для решения задач экологии.

Специальные приложения.

Значение функции

Значение критерия

Значение критерия

Критические значения коэффициента корреляции rk;α

2. Основы теории подобия

2.1. Подобие физических явлений и его признаки

2.2. Анализ размерностей

2.3. Первая теорема подобия

2.4. Применение методов подобия в математическом

11.3. Численные методы решения дифференциальных уравнений

11.3.1. Постановка задачи

11.3.2. Процесс численного решения

11.3.3. Метод Эйлера

11.3.4. Модифицированный метод Эйлера

11.3.5. Метод Рунге – Кутта

11.3.6. Метод Рунге – Кутта для систем дифференциальных уравнений

11.3.7. Общая характеристика одношаговых методов

3.8. Многошаговые методы

11.3.9. Методы прогноза и коррекции

11.3.10. Краткая характеристика методов прогноза и коррекции.

11.3.11. Выбор шага и погрешность решения.

11.3.12. Жесткие задачи

11.4. Имитационное моделирование систем

11.4.1. Принципы имитационного моделирования

11.4.2. Объекты моделирования

11.4.3. Динамическая модель исследуемого объекта

11.4.4. Построение имитационных моделей динамических систем

11.4.5. Преобразование передаточных функций звеньев в дифференциальные уравнения в форме Коши

11.4.6. Синтез имитационной модели на основе структурной схемы

11.5. Теоретические основы построения математических моделей систем

11.5.1. Компонентные и топологические уравнения

11.5.2. Компонентные и топологические уравнения механической системы

11.5.3. Компонентные и топологические уравнения электрической системы

11.5.4. Компонентные и топологические уравнения гидравлической системы

11.5.5. Компонентные и топологические уравнения тепловой системы

11.6. Метод электроаналогий

11.6.1. Сущность метода электроаналогий.

11.6.2. Электромеханические аналогии

11.6.3. Построение имитационных моделей методом электроаналогий

11.6.4. Плоское прямолинейное движение звеньев

11.6.5. Электрогидравлические аналогии

11.6.6. Электротепловые аналогии

Литература

Составим уравнение равновесия сил:

(6.13)

где

Составим уравнение равновесия моментов:

MT1+ MT2 - MT3 - Mu = 0, (6.14)

где MТ1 = T1 × l; MТ2 = T2 × l/2; MТ3 = T3 × l.

На основании данных уравнений составляем эквивалентную электрическую схему (рис. 6.30).

Рис. 6.30

Электрическая схема содержит 5 контуров, 3 узла и 3 трансформатора. Уравнения Кирхгофа для контуров 1-5 совпадают с уравнениями равновесия сил и с уравнениями равновесия моментов (6.13-6.14). Записываем уравнения Кирхгофа для узлов 1, 2 и 3.

Записываем уравнения трансформаторов.

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

Рис. 6.31

Рассматриваемое твердое тело имеет две степени свободы, так как оно может совершать прямолинейное и одновременно вращательное движения. Поэтому чтобы построить сеть связей динамической системы, необходимо в двух уравнениях Кирхгофа для контуров вынести вправо инерционные составляющие. Такими уравнениями являются уравнения (6.13) и (6.14). Используя дополнительно уравнения Кирхгофа для узлов 1–3, а также уравнения трансформаторов, строим сеть связей динамической системы (рис.6.31).

Таблица 6.9

Передаточные функции элементарных звеньев

Теперь можно построить структурную схему динамической системы. С этой целью каждую ветвь сети связей заменяем элементарным звеном. Передаточные функции элементарных звеньев приведены в табл. 6.9, а структурная схема на рис. 6.32.


Рис. 6.32

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

На рис.6.33 представлена программа моделирования рассмотренного механического объекта в математической среде MathCAD.

Вал на трех опорах

Рис. 6.33. Код программы моделирования в системе MathCAD

Рис. 6.33 а. График процесса при h1=h2=h3=10.

Результаты моделирования представлены в виде графиков линейных (рис. 6.34, а)

Рис.6.34

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


11.6.5. Электрогидравлические аналогии

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

Рис. 6.35. Конструктивная схема гидростатической опоры

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

сопротивлением Rδ в карман-камеру. Из кармана масло под давлением вытесняется через зазор δ , гидравлическое сопротивление которого

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

При оптимальных соотношениях параметров гидростатических направляющих изменение толщины масляного слоя при изменяющейся нагрузке оказывается достаточно малым, что позволяет применять эти направляющие в высокоточных станках. Расход масла через дроссель – Q1. Через щель вытесняется количество масла Q3. При постоянном давлении насоса и сопротивлении дросселя расход Q1 постоянен, но при движении платформы вверх-вниз расход масла через щель изменяется. Это связано с изменением объема под платформой. При уменьшении

объема излишек масла вытесняется из-под платформы, а при его увеличении под платформой создается вакуум и масло «втягивается» в карман. Это приводит к появлению мнимого расхода масла Q2, который, складываясь с расходом Q1, дает расход масла на выходе Q3. Наличие этого мнимого расхода обусловливает уникальное свойство гидравлических систем – демпфирование колебаний. Распределение давления в зазоре показано на рис. 6.36.


Рис. 6.36. Распределение давления в зазоре

В соответствии с эпюрой результирующая гидростатическая сила в зазоре равна:

где SЭФ - эффективная площадь.

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

Рис. 6.37. Эквивалентная электрическая схема

Запишем уравнения Кирхгофа для контуров 1 и 2, а также для узла 1:

Уравнения трансформатора:

На основе составленных уравнений строим сеть связей физических величин (рис. 6.38).

Рис. 6.38. Сеть связей

Таблица 6.10

Передаточные функции элементарных звеньев

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

Рис. 6.39. Структурная схема

Запишем дифференциальные уравнения интегрирующих звеньев.

Далее формируем уравнения связи. С этой целью сначала определим давление в кармане PK:

Откуда:

Расчёт исходных параметров модели.

Дано:

Примем для установившегося режима


Тогда

Результаты моделирования гидростатических направляющих представлены в виде графиков, которые показывают реакцию системы на внешнее усилие FH= 3×104H (рис. 6.40)

Рис. 6.40