СОДЕРЖАНИЕ
Математическое моделирование элементов сложных экологических систем
Лекция 1. Введение в моделирование. Исторический экскурс.
1. Основы моделирования в экологии 1.1. Общие принципы построения моделей в экологии
2.2. Этапы построения математической модели
1.4. Элементы теории подобия, применяемые в моделировании
3.2.2. Выборочный метод в экологометрике.
Зависимость числа интервалов от объема выборки
Статистический ряд по интервалам
Лекция 4. Статистические оценки параметров распределения случайных величин по выборкам
4.4. Статистические оценки гипотез об экологических моделях
Выборка из генеральной совокупности
Статистическая таблица эксперимента
Пример преобразования членов уравнения регрессии
Вычисление данных для линеаризации уравнения регрессии
Нормальные уравнения мнк для некоторых функций
Статистическое оценивание уравнения регрессии и парной корреляции.
Обработка результатов наблюдений
Рекомендации по выбору вида функции
3.4. Динамические статистические модели
Данные по объему сброса качественных сточных вод
Данные по объему сброса сточных вод за 5-летие
Пример расчета 5-летних средних
Расчетные значения для определения уравнения динамики
Ряд динамики для определения сезонных колебаний
Эксперименталъный материал исследования
8.1. Анализ влияния отдельных факторов в экологической модели.
Лекция 9. Методы оптимизации. Метод Лагранжа
Лекция 10. Метод линейного программирования.
Лекция 11. Функциональные модели.
Лекция 12. Модели процессов содержащие обыкновенные дифференциальные уравнения.
Численные ошибки использованных для вычисления данных
Лекция 13. Статистические модели динамики.
Лекция 16. Информационные технологии в экологии. Экологические информационные системы.
1 6.1. Экологические информационные системы
1. Какова область значения для числовых характеристик?
Лекция 17. Использование информационных технологий для решения задач экологии.
Критические значения коэффициента корреляции rk;α
2.1. Подобие физических явлений и его признаки
2.4. Применение методов подобия в математическом
11.3. Численные методы решения дифференциальных уравнений
11.3.2. Процесс численного решения
11.3.4. Модифицированный метод Эйлера
11.3.6. Метод Рунге – Кутта для систем дифференциальных уравнений
11.3.7. Общая характеристика одношаговых методов
11.3.9. Методы прогноза и коррекции
11.3.10. Краткая характеристика методов прогноза и коррекции.
11.3.11. Выбор шага и погрешность решения.
11.4. Имитационное моделирование систем
11.4.1. Принципы имитационного моделирования
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.1. Сущность метода электроаналогий.
11.6.2. Электромеханические аналогии
11.6.3. Построение имитационных моделей методом электроаналогий
11.6.4. Плоское прямолинейное движение звеньев
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
Требуется построить имитационную модель теплового процесса, который происходит при нагреве штампа для горячего прессования пластмасс.

Рис. 6.41. Конструктивная схема штампа
Конструкция штампа содержит встроенный электронагревательный элемент (рис. 6.41). Рассмотрим процесс передачи тепла через тело штампа (рис.6.42).

Рис. 6.42. Процесс передачи тепла
Тепловой поток от нагревателя Ф1 проходит последовательно через цепь выделенных элементов. Поперечное сечение элементов заданно. Количество элементов в данном случае - 4. Проходя через теплопроводящую цепь элементов, тепловой поток преодолевает тепловое сопротивление этих элементов. Кроме того, часть теплового потока идет на потери тепла, связанные с теплоемкостью металла, из которого изготовлен штамп. Поэтому тепловой поток, проходящий через наружную поверхность Ф2< Ф1. Эквивалентная электрическая схема замещения представляет собой линию задержки, состоящую из последовательно включенных RC-цепочек (рис. 6.43). Активное сопротивление цепочки имитирует тепловое сопротивление, а электрическая емкость – теплоемкость элемента.

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

Рис. 6.44. Структурная схема
Таким образом, температура поверхности штампа не остается постоянной. Она периодически изменяется от минимального значения до максимального значения. Чтобы построить структурную схему математической модели, надо определить тепловой поток Ф1, проходящий через поперечное сечение элементов. Весь тепловой поток, излучаемый электрическим нагревателем, равен мощности этого нагревателя, а часть теплового потока, приходящаяся на поперечное сечение элемента, во столько раз меньше, во сколько раз поперечное сечение элемента меньше площади всей поверхности штампа.
Электрической схеме замещения соответствует структурная схема математической модели (рис. 6.44).
Непосредственно из структурной схемы следуют уравнения связи.


где
![]()
где, в свою очередь,
![]()
а T - полупериод рабочего цикла. Переходные процессы в интегрирующих звеньях определяются с помощью следующих дифференциальных уравнений.

Коэффициенты модели определим исходя из того, что материал штампа - железо. Для железа:

Примем
,
,![]()
Графики переходных процессов (рис. 6.45) показывают, что при заданном режиме работы температура поверхности штампа изменяется в пределах от 198оС до 144оС. Время нагрева штампа до заданной температуры составляет ок. 200 c .

Рис. 6.45
Рассмотренная методика синтеза имитационных моделей имеет как преимущества, так и недостатки. Преимущество этой методики в том, что синтез имитационной модели сопровождается построением структурной схемы динамической системы. Поскольку структурная схема содержит передаточные функции всех звеньев, то появляется возможность исследовать динамическую систему с привлечением методов теории управления. Однако такое моделирование оправдано лишь для сравнительно простых систем, содержащих небольшое количество звеньев. С повышением сложности системы количество передаточных функций и уравнений связи, которые рассчитываются вручную, возрастает, и работа становится громоздкой и трудоемкой. Чтобы облегчить синтез сложных моделей, применяют компьютерное моделирование, основанное на структурно - матричных представлениях метода электроаналогии, либо специальные программные средства моделирования, такие как Matlab, Mathcad, Anylogic и др.