Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Анализ статических и динамических характеристик простой гидравлической системы. Методические указания к расчётной работе по дисциплинам «Моделирование
.pdf
Приложение 3. Распечатка программы динамического моделирования
//Ввод известных значений давления в МПа (p1, p2 – давление входа, p3, p4 –
//давления выхода):
p1=1e6*input('p1=');
p2=1e6*input('p2=');
p3=1e6*input('p3=');
p4=1e6*input('p4=');
//Присвоение значений другим параметрам (pN – атмосферное давление
//в МПа, g – ускорение свободного падения в м/с^2,
//h1G – высота первого бака в м, h2G – высота второго бака в м, S – площадь
//основания баков в м^2, ro – плотность воды при данной температуре
//в кг/м^3, k(i) – сопротивление клапана):
g=9.81;
pN=0.1*1e6;
h1G=10;
h2G=10;
S1=1;
S2=1;
ro=1000;
for i=1:5
k(i)=0.0001*i;
end;
//Подпрограмма, реализующая правую часть системы уравнений
function dy=ff(t, y)
dy=zeros(2,1);
p7=pN*h1G/(h1G-y(1));
p5=p7+(ro*g)*y(1);
Q1=k(1)*sign(p1-p5)*sqrt(abs(p1-p5));
Q3=k(3)*sign(p5-p3)*sqrt(abs(p5-p3));
p8=pN*h2G/(h2G-y(2));
p6=p8+(ro*g)*y(2);
Q2=k(2)*sign(p2-p6)*sqrt(abs(p2-p6));
Q4=k(4)*sign(p6-p4)*sqrt(abs(p6-p4));
Q5=k(5)*sign(p5-p6)*sqrt(abs(p5-p6));
dy(1)=(Q1-Q5-Q3)/S1;
dy(2)=(Q2+Q5-Q4)/S2;
endfunction
//Расчёт заполнения ёмкостей
y0=[0;0];
t0=0;
t=0:0.001:20;
y=ode(y0,t0,t,ff);
subplot(1,2,1);
plot(t,y);
41
Окончание прил. 3
xgrid();
xtitle('Изменение высот жидкости в ёмкостях при их заполнении','t, c','h, м');
legend('h1','h2',4,%f);
subplot(1,2,2);
//Расчёт изменения высот жидкости при наличии возмущения в системе
//До появления возмущения
//Новые начальные условия
h1G=15;
y0=[14.5;1.50];
t0=0; t=0:0.1:15;
y=ode(y0,t0,t,ff);
//После появления возмущения
y0=[10;9.9]; //Параметры возмущения
t0=15; t1=15:0.1:25;
y1=ode(y0,t0,t1,ff);
t2=[t t1];
y2=[y y1];
plot(t2,y2);
xgrid();
xtitle('Изменение высот жидкости в ёмкостях при их заполнении при
возмущении в системе','t, с','h, м');
legend('h1','h2', 4,%f);
42
АНАЛИЗ СТАТИЧЕСКИХ И ДИНАМИЧЕСКИХ
ХАРАКТЕРИСТИК ПРОСТОЙ ГИДРАВЛИЧЕСКОЙ СИСТЕМЫ
Методические указания к расчётной работе по дисциплинам:
«Моделирование гидравлических систем», «Моделирование гидрогазодинамических
процессов»
Составитель Севостьянов Александр Владимирович.
Редактор М.Ю. Болгова
Подписано к печати 20.10.2014. Формат 60x84 1/16. Бумага офсетная.
Ризография. Объём 2,5 п.л. Тираж 50 экз. Заказ № .
Издательство Липецкого государственного технического университета.
Полиграфическое подразделение Издательства ЛГТУ.
398600, Липецк, ул. Московская, 30.
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]