Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

методические указание по лабораторым работам / МУ_Практическая работа_2 2013

.pdf
Скачиваний:
62
Добавлен:
20.02.2014
Размер:
122.71 Кб
Скачать

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

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

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

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

всех форм обучения

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

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

Балаково 2007

ВВЕДЕНИЕ

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

нентных и топологических уравнений широкое применение получил узло-

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

циденций, отражающую структуру связей всех элементов системы. Мат-

рица инциденций формируется на основании графических моделей техни-

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

дические указания являются третьей работой по разработке и анализу ма-

тематических моделей.

ЦЕЛЬ РАБОТЫ

Получение навыков построения математических моделей техниче-

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

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

Информация о математической модели гидравлической системы,

которую содержит орграф, можно представить в виде матрицы инцидент-

ности А, сформированной по следующему правилу. 1. Размерность матрицы α х β:

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

нием базового;

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

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

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

2

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

При наличии внешнего источника типа потока Q*В он формально заменя-

ется на источник типа потенциала P*B и добавляется условная масса m*B.

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

ветвей источников потенциалов АВ.

А = [АИ , АД, АУ , АВ ].

(1)

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

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

ричной форме.

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

 

 

 

 

 

= m

dQ И

.

(2)

 

 

РИ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

2.

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

 

 

 

 

 

 

 

 

 

 

 

 

 

= m ×

 

 

Д .

 

 

 

РД

Q

(3)

3.

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

 

 

 

 

 

 

 

 

 

 

 

 

= с

 

У dt .

 

 

 

РУ

Q

(4)

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

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

Q И , Q Д , Q У - векторы потоков (расходов) соответствующих элементов; m, μ, c - диагональные матрицы параметров элементов гидравлической системы.

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

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

(

 

j ) i = 0 , i =

 

.

(5)

Р

1, n

j

 

3

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

(

 

k )l = 0 ,

l =

 

.

(6)

Q

1, k

k

 

 

 

 

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

число

ветвей.

 

 

 

 

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

матричной форме:

 

 

 

А ×

Р

= 0 ;

 

 

(7)

 

 

 

 

 

B + AT

 

= 0 ;

В=И, Д, У,

(8)

 

 

 

Q

Q

где А –

матрица инциденций;

 

 

АТ

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

 

 

 

 

 

 

 

Q

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

 

 

 

 

 

 

 

 

Q

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

 

 

 

 

Отдельный элемент QB i вектора

Q

B определяется по формуле:

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

Q B i = ИBij × Q j

,

(9)

 

 

 

 

 

 

 

 

 

j=1

 

 

где ИBij

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

ой ветви орграфа j-ому узлу.

 

 

Вектор потенциалов системы Р представим состоящим из-под век-

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

действий на систему

 

 

 

 

 

 

 

 

 

Т = [

 

И ,

 

 

Д ,

 

У ,

 

В ].

 

 

 

 

(10)

 

 

 

 

 

 

 

Р

Р

Р

Р

Р

 

 

 

 

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

 

 

 

 

 

 

АИ ×

 

 

+ АД ×

 

 

+ АУ ×

 

 

+ АВ ×

 

= 0 .

 

(11)

 

РИ

РД

РУ

РВ

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

× m

dQ

 

+ АД × m

 

+ АУ × c

 

dt + А

В ×

 

 

= 0 .

(12)

АИ

И

QД

QУ

Р

В

 

 

 

 

dt

 

 

 

 

 

Учитывая топологическое уравнение (8), фазовые переменные QИ , QД , QУ

4

можно выразить через узловые потоковые переменные Q . В результате

получают матричное уравнение, являющееся классическим вариантом уз-

лового метода:

АИ × m × АТИ

dQ

+ А

Д × m × АТД

 

+ А

У × c × АСТ

 

 

 

 

(13)

Q

Qdt + АВ × РВ = 0 .

 

 

dt

 

 

 

 

 

 

 

 

 

 

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

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

1.Подставляют значения РУ = -АУ × c × АТС Qdt в уравнение (13) и разре-

шают его относительно производной.

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

dQ

=

 

1

 

 

 

В

 

B + A

У

 

 

У - АД × m × АТД

 

 

(14)

 

 

 

 

Q

Q

 

Q) ,

 

 

Т

dt А

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

= −cAТQ .

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

дение

 

 

 

 

 

 

 

 

 

 

И

Т )−1 = m−1

.

 

 

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ

= m

 

 

 

 

 

 

 

 

 

 

 

 

 

(17)

 

в

Рв + А

у

Ру + А

д

Рд) ;

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

= −cAТQ ,

 

 

 

 

 

(18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д = −μАТ

 

.

 

 

 

 

 

(19)

где

 

 

P

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д

 

 

 

 

 

 

 

 

Уравнения (17)-(18) представляют собой модифицированный узло-

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

5

ЗАДАНИЕ

1. Сформировать матрицу инциденций на основе ориентированного графа,

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

2. Используя модифицированный узловой метод формирования моделей,

разработать математическую модель гидравлической системы в статиче-

ском и динамическом режимах.

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

На основании орграфа (рис.1) сформируем матрицу инциденций

(таблица 1), по правилам, изложенным в основных теоретических сведени-

ях.

с2

 

4

 

 

с2

μ4

рв1

 

μ3

m4

μ5

5

 

 

 

m5

3

 

 

 

 

 

m3

 

 

 

0

рв2

 

с1

m2

рн1

 

 

 

μ2

m1

 

 

 

 

 

рн2

μ1

 

с1

 

 

 

2

 

1

 

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

6

Таблица 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

 

 

 

АВ

 

 

 

 

АУ

 

 

АД

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании матрицы инциденций записывают подматрицы упру-

гих АУ, диссипативных АД элементов и подматрицу источников потенциа-

лов АВ:

 

 

1 0 0 0

0

 

 

−1 0

0

0

0

 

 

 

−1 0

 

 

0 1 0 0

0

 

 

0

−1 0

0

0

 

 

 

−1 0

A

в

= 0 0 0 0

0

 

,

А =

0

0

−1 0

0

 

,

А

у

=

1

−1 .

 

 

 

 

 

д

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0 −1

0

 

 

 

 

0

0

0

−1

0

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0 0 0

−1

 

 

 

0

0

0

0

−1

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

m

0

0

0

0

 

 

 

 

 

 

 

 

μ

0

0

0

0

 

 

 

 

 

 

 

 

 

1

μ

 

 

 

 

 

 

 

 

1

 

 

0

0

0

 

 

c 0

 

 

 

 

0

2

0

0

0

 

 

0

m

2

 

,

,

μ =

 

 

 

μ

 

 

 

 

.

 

0

 

m3

0

0

 

c =

1

 

 

0

0

3

0

0

m =

0

 

 

0 с

 

 

 

 

 

 

 

 

 

μ4

 

 

 

0

0

0

m4

0

 

 

 

2

 

 

 

0

0

0

0

 

0

0

0

0

m

5

 

 

 

 

 

 

 

 

 

0

0

0

0

μ

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

7

 

РН1

 

 

 

 

 

 

 

 

 

Р

Д1

 

 

 

 

Q1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

Н2

 

 

 

Р

 

 

 

 

 

Р

Д2

 

 

 

Q2

Рв =

0

 

 

 

 

 

РД = Р

 

 

 

Q = Q3

.

 

 

 

 

 

 

 

,

Ру =

У1

 

,

Д3

,

 

 

 

 

 

 

 

Р

У2

 

 

 

 

 

 

 

 

 

 

 

 

 

РВ1

 

 

 

 

 

 

 

 

 

Р

Д4

 

 

 

Q4

 

 

Р

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В2

 

 

 

 

 

 

 

 

 

Р

Д5

 

 

 

Q5

 

Вычислим матричное произведение слагаемых правой части урав-

нения (17):

 

 

 

 

 

 

 

РН1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

PН2

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

 

 

m 2

 

 

 

 

 

 

 

 

 

m

A

 

Р

=

 

 

 

0

 

 

 

 

В

В

 

 

− PВ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 4

 

 

 

 

 

 

 

 

 

 

 

 

− PВ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− Ру1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− Ру1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 2

 

 

 

 

 

 

 

 

 

−1

 

 

 

РУ1 − РУ2

m

 

A

Р

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

у

 

 

 

m3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РУ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 4

 

 

 

 

 

 

 

 

Р

У2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m5

 

 

 

 

 

 

 

 

,

 

 

 

 

 

− Р

д1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

− Р

1

 

 

 

 

 

 

 

 

д2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

2

 

 

 

 

 

 

 

 

 

 

−1

 

 

 

− Рд3

 

 

 

 

 

m

А

д

Р

=

 

 

 

 

 

 

 

 

 

 

 

д

 

m

3

 

 

 

 

 

 

 

− Р

 

 

 

 

 

 

 

 

д4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

4

 

 

 

 

 

 

− Р

 

 

 

 

 

 

 

д5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

сАТ Q = с1 (−Q1 − Q2 + Q3 )

Уc2 (−Q3 + Q4 + Q5 )

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

и учитывают зависимость (18):

8

 

 

 

dQ

1

 

 

=

 

Р

Н1

 

 

РД1

Р

У1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

m1

 

m1

 

m1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ

2

 

=

 

Р

Н2

 

РД2

 

 

 

Р

У1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

m 2

 

m2

 

m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= − РД3 + РУ1 РУ2

 

dQ3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

m3

 

 

 

 

 

m3

 

 

 

 

 

 

m3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ

4

 

 

 

 

 

 

 

 

Р

В1

 

 

 

 

РД4

 

 

 

 

 

 

Р

У2

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

m 4

 

 

 

 

 

m 4

 

 

 

 

 

 

m 4

(20)

 

 

 

 

 

 

 

 

= − РВ2 РД5 + РУ2

dQ5

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

m5

 

 

 

 

 

m5

 

 

 

 

 

 

m5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У1

= С1 (−Q1 − Q2 + Q3 )

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У2

 

 

= С2 (−Q3 + Q4 + Q5 )

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Р

Д1

= m× Q1 + m× Q12 × signQ1

 

 

= m× Q2

+ m× Q22 ×signQ2

Р

Д2

 

 

= m3 × Q3 2 ×signQ3

 

РД3

(21)

 

 

 

 

 

= m× Q4

+ m× Q42 ×signQ4

РД4

 

 

= m× Q5

2

× signQ5

РД5

+ m× Q5

С учетом численных значений параметров модели, найденных во

второй работе, уравнения (20) приобретает вид:

dQ

1

 

=

 

Р

Н1

 

 

 

-

 

 

РД1

 

 

-

 

Р

У1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

6

 

10,95 ×10

6

 

 

 

 

6

 

 

 

 

 

 

 

 

10,95 ×10

 

 

 

 

 

 

10,95 ×10

 

 

 

 

dQ 2

=

 

Р

Н2

 

 

 

-

 

РД2

 

 

-

 

Р

У1

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

×106

10,95 ×106

 

 

 

 

 

10,95 ×106

 

 

 

10,95

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

 

 

 

 

 

 

 

Р

 

 

 

 

 

 

 

Р

 

 

 

 

 

 

dQ

3

= -

 

Д3

 

 

 

 

+

У1

 

 

 

 

-

 

У2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(22)

 

dt

 

 

2,49

×10

6

 

2,49

×106

 

 

2,49

×106

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ

4

 

= -

 

Р

В1

 

 

 

 

-

Р

Д4

 

 

 

 

+

 

Р

У2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

21,9

×10

6

 

21,9

×10

6

 

 

21,9

×10

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ5

 

 

= -

 

Р

В2

 

 

 

 

-

Р

Д5

 

 

 

+

 

Р

У2

 

 

 

 

 

 

 

dt

 

 

 

32,9

×10

6

 

32,9

×10

6

 

 

32,9

×10

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

У1

 

= 2,08 ×1012 (-Q1

- Q 2

+ Q3 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У2

 

= 0,8 ×1012 (-Q 3

 

+ Q 4

+ Q5 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

РД1

= 5,27 ×107

 

 

× Q1

+ 26,48 ×1010

 

× Q12

 

× signQ1

 

 

 

= 5,27 ×10

7

 

 

× Q 2

+ 29,97 ×1010

 

× Q 22

× signQ2

 

 

РД2

 

 

 

 

 

 

= 430 ×1010

 

× Q3

2

× signQ3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(23)

РД3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 10,54 ×10

7

× Q 4

+ 52,97 ×10

10

× Q

2

× signQ4

 

 

РД4

 

 

 

 

 

4

 

 

 

= 15,81×10

7

 

× Q 5

+ 72,48

×10

10

× Q

2

× signQ5

 

 

РД5

 

 

 

 

 

 

5

 

 

Для рассмотрения статического состояния система дифференциаль-

ных уравнений (20) приобретает вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 = РН1 − РД1 − РУ1

 

 

 

 

 

 

 

 

 

 

 

 

 

0 = Р

Н

2

− Р

Д2

− Р

У1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

= −РД3 + РУ1 − РУ2

 

 

 

 

 

 

 

 

 

 

 

 

 

= −РВ1 − РД4 + РУ2

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

= −РВ2 − РД5 + Р

 

 

 

 

 

 

 

 

 

 

(24)

 

 

 

 

0

У2

 

 

 

 

 

 

 

 

 

 

 

 

0 = −Q1 − Q2 + Q3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 = −Q3 + Q4 + Q5

 

 

 

 

 

 

 

 

 

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

 

 

5,27 ×107

× Q1

 

+ 26,48 ×1010

× Q12

 

× signQ1

+ РУ1

= РН1

 

5,27 ×107

× Q 2

 

+ 29,97 ×1010

× Q 22

 

× signQ 2

+ РУ1

= РН2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×10

10

× Q3

2

 

× signQ 3 - РУ1

+ РУ2

= 0

 

 

 

 

430

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 52,97 ×1010

× Q 42 × signQ 4 - РУ2 = РВ1

(25)

10,54 ×107 × Q 4

 

 

 

 

 

7

× Q 5

 

+ 72,48 ×

 

 

10

 

2

 

× signQ 5 - РУ2

= РВ2

 

5,81 ×10

 

 

10

 

× Q 5

 

 

- Q1 - Q 2 + Q 3 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- Q

3 + Q 4 + Q5 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, система уравнений (22), (23) является динамической моде-

лью гидравлической системы, а (25) – статической моделью.

10