Моделирование простой гидравлической системы в стационарном режиме. Рекомендации к выполнению лабораторной работы №2

1. Методические рекомендации к выполнению лабораторной работы № 2 по теме «Моделирование простой гидравлической системы в

динамическом режиме»
Преподаватели:
1. Советин Филипп Сергеевич, кандидат технических
наук, доцент;
2. Царёва Елена Владимировна, кандидат технических
наук, доцент;
3. Павлов Александр Сергеевич, старший
преподаватель.

2. Создать папкуC:\Users\HP\Documents\О-33

3. Разместить в своей папке программы для выполнения ЛР

4. Запустить пакет МАТЛАБ Помнить, что МАТЛАБовские работы открывать только в среде МАТЛАБ

5.

В среде МАТЛАБ открыть программный код для
выполнения ЛР № 2
Главная
Программа и
м-файл
Функция
должны
находиться в
одной папке
Пока что необходимо открыть только два
выделенных файла

6.

для переключения на текущую директорию
необходимо запустить программу

7.

Для этого курсор должен быть в программе главной программе
gidrstat и запустить следующим образом

8.

Программа требует сменить директорию, сменить её

9.

Был проведён запуск программы моделирования ПГС, той, ПГС,
которая была разобрана на семинаре.
Результаты выводятся в командное окно пакета МАТЛАБ

10. Изменение программы для моделирования своей ПГС, в соотв. своему варианту

ФРАГМЕНТ ПРОГРАММЫ (гл. программы)
%np-количество давлений (всего);nk-количество вентилей;
%nv-кол-во ур. без учёта дифференциальных
np=8; nk=5; nv=11; s=[1,1]; g=9.815;
disp ('Высота емкостей'); hg=[10,10];
disp ('плотность (кг/м3)'); ro=1000;
disp ('Начальное давление (Па)'); pn=100000;
disp ('Площадь внутреннего проходного сечения трубопровода (м^2)'); S=0.01;
disp ('Давление (1-4 5-8)'); p=[1000000, 1000000, 100000, 100000, 0, 0, 0, 0];
disp ('Коэф. пропускной способности (1-5)'); k=[0.05, 0.01, 0.01, 0.01, 0.03];
Синим цветом выделен кусок программы, который изменяется

11.

%ВНИМАНИЕ - НИ В КОЕМ СЛУЧАЕ НАЗВАНИЕ м-ФАЙЛА-ФУНКЦИИ НЕ МЕНЯТЬ!!!
function F=fpr(t,h)
global ro pn p ak vm hg g v s
%Давления газов над жидк. в закрытых емк.
p(7) = pn * hg(1) / (hg(1) - h(1));
p(8) = pn * hg(2) / (hg(2) - h(2));
%Давления жидк. в закрытых емк.
p(5) = p(7) + ro * g * h(1);
p(6) = p(8)+ ro * g * h(2);
%Об. расх. ж-ти через вентили
v(1) = ak(1) * sign(p(1) - p(5)) * sqrt(abs(p(1) - p(5)));
v(3) = ak(3) * sign(p(5) - p(3)) * sqrt(abs(p(5) - p(3)));
v(5) = ak(5) * sign(p(5) - p(6)) * sqrt(abs(p(5) - p(6)));
v(2) = ak(2) * sign(p(2) - p(6)) * sqrt(abs(p(2) - p(6)));
v(4) = ak(4) * sign(p(6) - p(4)) * sqrt(abs(p(6) - p(4)));
%Правые части д/у (не забывать учитывать знаки, это важно)
F=[(v(1)-v(3)-v(5))/s(1); (v(2)+v(5)-v(4))/s(2)];
vm= ro*v
Мои комментарии читать внимательно. После исправления
необходимо сохранить изменения в м-файле функции

12.

ФРАГМЕНТ ГЛАВНОЙ ПРОГРАММЫ
Начальные условия по высотам
H10=0;H20=0 (обе высоты равны нулю в
режиме «пуска остановки»)
Y0=[H10;H20];
H10=0;H20=0 – в режиме пуска установки.
Прошу обратить внимание, что при моделировании
динамического режима в режиме слива или
«остановки системы» необходимо задать начальные
условия по высотам столбов жидкости в емкостях
H10 H20
Далее необходим запуск изменённой программы для моделирования
своей гидравлической системы (след. слайд)

13.

14. Результаты компьютерного моделирования простой гидравлической системы в динамическом режиме