методические указание по лабораторым работам / Узловой метод_методичка1

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

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

Саратовский государственный технический университет

Балаковский институт, техники, технологии и управления

УЗЛОВОЙ МЕТОД ФОРМИРОВАНИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ

Методические указания к выполнению практических работ

по курсу "Моделирование систем" для студентов специальности 210100 всех форм обучения

Одобрено

редакционно-издательским советом Балаковского института техники,

технологии и управления

Балаково 2007

ВВЕДЕНИЕ

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

ЦЕЛЬ РАБОТЫ

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

ОСНОВНЫЕ ПОНЯТИЯ

Информация о математической модели гидравлической системы, которую содержит орграф, можно представить в виде матрицы инцидентности А, сформированной по следующему правилу.

1. Размерность матрицы α х β:

а) α - число строк соответствует каждому узлу орграфа, за исключением базового;

б) β - число столбцов соответствует ветвям орграфа.

2. Элементы матрицы:

а) отсутствие связи между узлом и ветвью – "0";

б) если ветвь входит в узел – "1";

в) если ветвь выходит из узла – "-1".

При наличии внешнего источника типа потока Q^*_В он формально заменяется на источник типа потенциала P^*_B и добавляется условная масса m^*_B.

Матрицу инциденций А можно представить состоящей из подматриц инерционных А_И, диссипативных А_Д, упругих А_У ветвей и подматрицы ветвей источников потенциалов А_В.

. (1)

Подматрица инерционных элементов А_И всегда единичная диагональная и может быть не отражена при составлении матрицы.

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

1. Для инерционных элементов

. (2)

2. Для диссипативных элементов

. (3)

3. Для упругих элементов

. (4)

В формулах (2)-(4) использованы следующие обозначения:

- векторы потенциалов (давлений) соответствующих элементов;

- векторы потоков (расходов) соответствующих элементов;

- диагональные матрицы параметров элементов гидравлической системы.

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

1. Уравнение равновесия потенциалов ветвей орграфа, инцидентных i-ому узлу:

, . (5)

2. Уравнение непрерывности фазовой переменной типа потока l-ой ветви:

, . (6)

В формулах (5), (6) n – число узлов орграфа, за исключение базы, l – число ветвей.

Используя матрицу инциденций уравнения (5), (6), можно записать в матричной форме:

; (7)

; В=И, Д, У, (8)

где А – матрица инциденций;

А^Т – транспонированная матрица инциденций;

- вектор потоковых переменных ветвей;

- вектор потоковых переменных узлов.

Отдельный элемент вектора определяется по формуле:

, (9)

где - элемент матрицы инциденций, характеризующий соответствии i-ой ветви орграфа j-ому узлу.

Вектор потенциалов системы представим состоящим из-под векторов инерционных, диссипативных, упругих элементов и внешних воздействий на систему

. (10)

Тогда, используя (1), уравнение (7) запишется в виде

. (11)

С учетом (2)-(4) уравнение (11) принимает вид

. (12)

Учитывая топологическое уравнение (8), фазовые переменные , , можно выразить через узловые потоковые переменные . В результате получают матричное уравнение, являющееся классическим вариантом узлового метода:

. (13)

Для перехода к системе дифференциальных уравнений в нормальной форме Коши выполняют следующие преобразования:

1. Подставляют значения в уравнение (13) и разрешают его относительно производной.

2. Определяют производную по времени от вектора и сводят в единую систему:

, (14)

. (15)

Так как матрицы А_И и m одного и того же порядка, то матричное произведение

. (16)

Тогда система (14), (15) запишется в виде:

; (17)

, (18)

где . (19)

Уравнения (17)-(18) представляют собой модифицированный узловой метод, позволяющий получить динамическую модель в нормальной форме Коши.

ЗАДАНИЕ

1. Сформировать матрицу инциденций на основе ориентированного графа, построенного при выполнении практической работы "Графические формы формирования математической модели гидравлической системы" согласно вариантам заданий.

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

ПРИМЕР ВЫПОЛНЕНИЯ ЗАДАНИЯ

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

Рис. 1. Орграф гидравлической системы

Таблица 1

Матрица инциденций гидравлической системы

Узлы	Ветви
	Источники потенциалов					Упругие		Диссипативные
	рН1	рН2		рВ1	рВ2	с1	с2	μ1	μ2	μ3	μ4	μ5
1	1	0	0	0	0	-1	0	-1	0	0	0	0
2	0	1	0	0	0	-1	0	0	-1	0	0	0
3	0	0	0	0	0	1	-1	0	0	-1	0	0
4	0	0	0	-1	0	0	1	0	0	0	-1	0
5	0	0	0	0	-1	0	1	0	0	0	0	-1
_	АВ					АУ		АД

На основании матрицы инциденций записывают подматрицы упругих А_У, диссипативных А_Д элементов и подматрицу источников потенциалов А_В:

, , .

Матрицы параметров инерционных, упругих и диссипативных элементов гидравлической системы соответственно:

, , .

Матрица потенциалов источников Р_в, упругих Р_у и диссипативных

Р_д, элементов и матрица фазовых переменных типа потока Q:

, , , .

Вычислим матричное произведение слагаемых правой части уравнения (17):

,

Складывают полученные матрицы в соответствии с уравнением (17) и учитывают зависимость (18):

(20)

Давления диссипативных элементов на основании формулы (19):

(21)

С учетом численных значений параметров модели, найденных во второй работе, уравнения (20) приобретает вид:

(22)

А давления диссипативных элементов:

(23)

Для рассмотрения статического состояния система дифференциальных уравнений (20) приобретает вид:

(24)

Выполнив преобразования (24) с учетом (23), получают:

(25)

Таким образом, система уравнений (22), (23) является динамической моделью гидравлической системы, а (25) – статической моделью.

ВОПРОСЫ ДЛЯ САМОПРОВЕРКИ

1. Матрица инциденций. Алгоритм построения.

2. Классический вид уравнений узлового метода.

3. Уравнения модифицированного узлового метода.

4. Статическая и динамическая модель.

ВРЕМЯ, ОТВЕДЕННОЕ НА ВЫПОЛНЕНИЕ РАБОТЫ

1. Подготовка к работе – 1,0 акад.час.

2. Выполнение работы – 2,0 акад.час.

3. Оформление работы – 1,0 акад.час.

ПОРЯДОК ОФОРМЛЕНИЯ ОТЧЕТА

1. Название практической работы.

2. Цель работы.

3. Задание.

4. Краткие теоретические сведения.

5. Расчеты с выводами.

ЛИТЕРАТУРА

1. Тарасик В.П. Математическое моделирование технических систем: учебник для вузов / В.П.Тарасик. – Мн.: ДизайнПРО, 2004. – 640с.: ил.

2. Самарский А.А. Математическое моделирование: идеи, методы, примеры / А.А.Самарский, А.П.Михайлов. – М.: Физматлит, 2001. - 320с.: ил.

СОДЕРЖАНИЕ

1. Введение 2

2. Цель работы 2

3. Основные понятия 2

4. Задание 6

5. Пример выполнения задания 6

6. Вопросы для самопроверки 11

7. Время, отведенное на выполнение работы 11

8. Порядок оформления отчета 11

Литература

УЗЛОВОЙ МЕТОД ФОРМИРОВАНИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ

Методические указания к выполнению практических работ

по курсу "Моделирование систем" для студентов специальности 210100 всех форм обучения

Составили: Мефедова Юлия Александровна

Власов Вячеслав Викторович

Рецензент Т.Н. Скоробогатова

Редактор Л.В.Максимова

Корректор Н.Т.Мальчикова

Подписано в печать		Формат 60х84 1/16
Бумага тип.	Усл. печ. л. 0,75	Уч.-изд.л. 0,7
Тираж 100 экз.	Заказ	Бесплатно

Саратовский государственный технический университет

410054, г.Саратов, ул.Политехническая, 77

Копипринтер БИТТиУ, 413840, г.Балаково, ул. Чапаева, 140

13