2 МОДЕЛИРОВАНИЕ НА МАКРОУРОВНЕ

2.1 Исходные данные

Дана схема гидравлической системы, представленная на рисунке 4. В системе используется в качестве рабочей жидкости вверенное масло АУ. Материал трубопровода – сталь. Основные параметры системы и жидкости приведены в таблице 1. Параметры трубопровода приведены в таблице 2.

P_B1

P_B2

Qn

P_B3

P_B4

Рисунок 4 – Схема гидравлической системы

1, 2, 3, 4 - магистрали потребителей; P_B1, P_B2, P_B3, P_B4- давление потребителей; Q_H – насос.

Таблица 1 – Параметры системы и жидкости

Обозначение	Основные параметры	Значение
	Плотность рабочей жидкости	860 кг/м^3
	Вязкость	0,15·10^-4 м^2/с
ЕС	Модуль упругости системы	1,7·10^8 Па
Етр	Модуль упругости трубопровода	9·10^10 Па
	Коэффициент потерь на трение при турбулентном потоке	0,028
	Толщина стенки трубопровода	2,2·10^-^3 м

^{Таблица 2 – Параметры трубопровода}

^Параметр	^Обозначение	^Номер трубопровода
		^1	^2	^3	^4	^5
^Диаметр трубопровода	^dтр^, м	^0,015	^0,015	^0,01	^0,01	^0,015
^Длина трубопровода	^l^,м	^1,5	^2	^2	^2	^1
^Коэффициент местных сопротивлений	^ξ	^4,5	^5	^5	^5	^3
^Давление потребителей и насосов	^p^, ×^10^6^Па	^0,1	^0,23	^0,17	^0,15	^-

Q_Н1 = 50×10^-6 м³/с;

Q_Н2 = 300×10^-6 м³/с.

2.2 Графические формы математической модели гидросистемы

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

2.2.2 Орграф. На основании динамической модели построен ориентированный граф, являющийся графической формой модели гидравлической системы. Узлы орграфа соответствуют сосредоточенным массам, а ветви – компонентам математической модели.

µ₁

µ₃₃

µ₅

m₂

m₁

m₃

m₄

C₁

Q_н

Q₃

Q₁

Q₂

Q₄

P_B4

P_B3

P_B2

m₅

µ₄

µ₄

Рисунок 5 – Динамическая модель гидравлической системы

Базовый узел (с нулевым номером) отображает инерциальную систему отсчета фазовых координат типа потока. Источник обеспечивает возрастание потоковой переменной узла, поэтому сигналы направляют от базы к узлу. В магистралях потребителей – наоборот. Во всех ветвях инерционных и диссипативных элементов направление сигналов от узла к базе. Такое направление характеризует затраты энергии источников на увеличение кинетической энергии и на трение. В ветвях упругих компонентов стрелки указывают направление передачи энергии от источников к потребителям. В ветвях всех элементов кроме направления записывается параметр каждого элемента. На рисунке 6 представлен полученный орграф.

2.2.3 Матрица инциденций. Для формирования полной математической модели на основе компонентных и топологических уравнений широкое применение получил узловой метод, для него необходимо сформировать матрицу инциденций, отражающую структуру связей всех элементов системы. Матрица инциденций формируется на основании ориентированного графа. Число строк матрицы соответствует числу узлов орграфа, число столбцов – числу ветвей. Отсутствие связи между узлом и ветвью обозначается «0», если ветвь в ходит в узел – «1», если выходит – «-1».

Рисунок 6 – Ориентированный граф гидравлической системы

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

А=[A_И, А_Д, А_У, А_В] (23)

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

Узлы	Ветви
	Диссипативные					Упругие	Источники потенциалов
№	μ1	μ 2	μ3	μ 4	μ 5	С1	Рв1	Рв2	Рв3	Рв4	Q5
1	-1	0	0	0	0	1	-1	0	0	0	0
2	0	-1	0	0	0	1	0	-1	0	0	0
3	0	0	-1	0	0	1	0	0	-1	0	0
4	0	0	0	-1	0	1	0	0	0	-1	0
5	0	0	0	0	-1	-1	0	0	0	0	1
Подматрица	АД					АУ	АВ

2.3 Узловой метод формирования математической модели гидросистемы

Из матрицы инциденций можно получить систему равнений (24), математически описывающие функционирование гидравлической системы:

(24)

где ;

А_Д, А_У, А_В – подматрицы инциденций;

- векторы давлений;

- векторы расходов,

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

Для нашего случая система будет иметь вид:

(25)

Так как в исходной системе насос постоянной производительности, то = 0 и пятое уравнение (25) преобразуется к виду:

P_H = P_Д5 + P_У1 (26)

С учетом возможного турбулентного режима течения рабочей жидкости, а так же возможных местных потерь, связанных с изгибом трубопровода, зависимость для диссипативных элементов будет носить более сложный характер:

(27)

где μ_л - коэффициент гидравлического сопротивления, характеризующий

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

μ_н - коэффициент гидравлического сопротивления, характеризующий нелинейные потери при турбулентном режиме, по длине и местные.

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

Вычисление параметров трубопровода гидросистемы

Значения коэффициентов линейных и нелинейных потерь для конкретной магистрали находят по формулам:

- площадь сечения трубопровода, м²:

; (28)

- коэффициент линейных потерь, H·с/м⁵;

; (29)

- коэффициент нелинейных потерь, H·с/м⁵.

; (30)

- коэффициент жесткости участка:

; (31)

где Ψ - доля объема трубопровода;

V_тр - объем трубопровода, м³:

V_тр=S_тр·l . (32)

Доля объема трубопровода рассчитывается как отношение объема отдельного участка к сумме объемов всех n соединенных между собой участков:

(32)

где - объем трубопровода i-ого участка, м³.

Коэффициент жесткости упругого элемента:

. (34)

По исходным данным и полученным результатам находим жесткость упругого элемента c₁= 6,755·10¹¹ Н/м⁵.

Коэффициент массы:

. (35)

Полученные результаты для отдельных участков трубопровода приведены в таблице 4.

Таблица 4 – Параметры трубопровода гидросистемы

Параметр	Номер магистрали
	1	2	3	4	5
Площадь сечения трубопровода, Sтр, ×10^-4м^2	1,767	1,767	0,785	0,785	1,767
Объем трубопровода, Vтр, ×10^-4м^3	2,651	3,534	1,571	1,571	1,767
Доля объема трубопровода,	0,239	0,319	0,142	0,142	0,159
Коэффициент массы, mг, ×10^6кг/м^4	7,3	9,733	21,9	21,9	4,867
Коэффициент линейных потерь,, ×10^7H·с/м^5	1,561	2,082	10,54	10,54	1,041
Коэффициент нелинейных потерь,, ×10^11H·с/м^5	1,005	1,202	7,388	7,388	0,670
Коэффициент жесткости участка, cг, ×10^12Н/м^5	2,649	1,490	7,576	7,576	5,960