Методические рекомендации к выполнению лабораторной работы № 2 по теме «Моделирование простой гидравлической системы в динамическом режиме» - Преподаватели:
- Советин Филипп Сергеевич, кандидат технических наук, доцент;
- Царёва Елена Владимировна, кандидат технических наук, доцент;
- Павлов Александр Сергеевич, старший преподаватель.
Создать папкуC:\Users\HP\Documents\О-33 Разместить в своей папке программы для выполнения ЛР Запустить пакет МАТЛАБ Помнить, что МАТЛАБовские работы открывать только в среде МАТЛАБ В среде МАТЛАБ открыть программный код для выполнения ЛР № 2 - В среде МАТЛАБ открыть программный код для выполнения ЛР № 2
- Главная
- Программа и
- м-файл
- Функция
- должны
- находиться в
- одной папке
- для переключения на текущую директорию
- необходимо запустить программу
- Для этого курсор должен быть в программе главной программе gidrstat и запустить следующим образом
- Был проведён запуск программы моделирования ПГС, той, ПГС,
- которая была разобрана на семинаре.
- Результаты выводятся в командное окно пакета МАТЛАБ
Изменение программы для моделирования своей ПГС, в соотв. своему варианту - ФРАГМЕНТ ПРОГРАММЫ (гл. программы)
- %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];
%ВНИМАНИЕ - НИ В КОЕМ СЛУЧАЕ НАЗВАНИЕ м-ФАЙЛА-ФУНКЦИИ НЕ МЕНЯТЬ!!! - %ВНИМАНИЕ - НИ В КОЕМ СЛУЧАЕ НАЗВАНИЕ м-ФАЙЛА-ФУНКЦИИ НЕ МЕНЯТЬ!!!
- 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
- ФРАГМЕНТ ГЛАВНОЙ ПРОГРАММЫ
- Начальные условия по высотам
- H10=0;H20=0 (обе высоты равны нулю в
- режиме «пуска остановки»)
- Y0=[H10;H20];
- H10=0;H20=0 – в режиме пуска установки.
- Прошу обратить внимание, что при моделировании
- динамического режима в режиме слива или
- «остановки системы» необходимо задать начальные
- условия по высотам столбов жидкости в емкостях
- H10 H20
- Далее необходим запуск изменённой программы для моделирования своей гидравлической системы (след. слайд)
Результаты компьютерного моделирования простой гидравлической системы в динамическом режиме
Do'stlaringiz bilan baham: |
