Моделирование систем. Елизаров И.А - 93 стр.

UptoLike

91
t_v(1)=173.9754;
t_k(1)=278.1349;
p(1)=8.9951e+005;
M_v(1)=2.5e+003;
pi(1)=p(1);
M_p(1)=32.6929;
% описание теплофизических характеристик
c_v=1*4.19e3; % удельная теплоемкость воды
c_k=1500; % удельная теплоемкость газовоздушной среды
M=18e-3; % молярная масса воды
R=8.31; % универсальная газовая постоянная
ro_p=8.4; % плотность пара
ro_k=1; % плотность газовоздушной среды
ro_v=1000; % плотность воды
r_t=3.6e7; % удельная теплота сгорания топлива
Lamda=1945e3; % удельная теплота парообразования
% описание характеристик парового котла и значений входов
V_0=10; % общий объём котла
K=1500; % коэффициент теплопередачи
F=10; % поверхность теплопередачи
G_v_n=0.6032; % расход воды в номинальном режиме
t_v_vx=20; % температура воды на входе в котел
t_t_vx=20; % температура газовоздушной среды на входе
p_t=2*9.81e4; % давление газовоздушной среды на входе
p_pot=5*9.81e4; % давление пара в линии потребления
p_k=1*9.81e4; % давление в камере сгорания
% проводимости клапанов
mu1=0.0000008;
mu2_n=0.0000002;
% характеристики насыщенного пара
xx=[0.98 0.99 1 1.01 1.013 1.02 2 3 4 5 6 7 8 9 10 20 40]*1e5;
% давление насыщенного пара
yy=[99.23 99.46 99.7 99.93 100 100.27 120 133 143 151 158 164 169 174 179
212 250]; % температура насыщенного пара
% переменные для определения интегральных составляющих
ПИ-регуляторов
sum=0;
sum2=0;
% настройки ПИ-регуляторов
Kp2=0.000000000002;
Ti2=100;