7.3 Блок-схема алгоритма оптимизации энергетического режима ДСП-180

Краткое описание используемых переменных:

- t[i] - текущее время;

- X[i] - текущее значение входной величины;

- Х[0] - начальное значение входной величины в момент времени t[0];

- sigma - значение переключающей функции;

- Sigma - переменная хранящая значение переключающей функции до остановки ИМ;

- tk - время выдержки после остановки ИМ, по достижении которого осуществляется поверочный шаг;

- tv - время выдержки сигнум-реле;

- Y[i] - значение выходной величины после статического звена в момент времени t[i];

- Z1[i] - значение выходной величины после инерционного звена 1 порядка в момент времени t[i];

- Z[i] - значение выходной величины после звена запаздывания в момент времени t[i];

- f(x,t) - уравнение статической характеристики с учетом дрейфа;

- dZ1[i] - скорость изменения выходной величины после инерционного звена в момент времени t[i];

- dZ[i] - скорость изменения контролируемой выходной величины после звена запаздывания в момент времени t[i];

- dZmax - запомненное максимальное значение скорости dZ;

- Кu - скорость ИМ;

- T1 - постоянная времени инерционного звена;

- Т2 - время запаздывания;

- Zn - зона нечувствительности САО;

- U - значение управляющей функции;

Блок-схема алгоритма оптимизации энергетического режима ДСП-180 представлена на рисунке 7.1.

7.4 Программа расчета переходного процесса в САО с запоминанием максимума изменения выходной величины

Программа max_remember.pas предназначена для расчета переходного поискового режима работы САО с запоминанием максимальной скорости выходного параметра и остановкой ИМ в момент достижения максимальной скорости. В программе реализована возможность учета сложного дрейфа статической характеристики. Для предотвращения ложного реверса предусмотрено время выдержки сигнум-реле равное 0,2Т_о в течении которого вводится запрет на реверс ИМ. Для повышения устойчивости поиска в САО, реализующей используемый ОАУ в программе предусмотрен поверочный принудительный рабочий шаг или шаг с реверсом через время, 3T_о после остановки ИМ.

Рисунок 7.1 - Блок схема алгоритма оптимизации энергетического режима ДСП-180



Для предотвращения колебательного режима в области экстремума оптимизируемого параметра предусмотрена остановка САО в точке предполагаемого максимума[12].

Результатом выполнения программы является:

- построение таблицы рассчитанных значений величин Х(ф),Y(ф), Z(ф) и dZ(ф)/dф;

- построение графика изменения во времени входного X(t), выходного Z(t) и установившегося Y(t) значений параметров в текущий момент времени;

- построение траектории поискового режима САО в осях X(ф) - Y(ф), Z(ф).

Текст рабочей программы расчета переходного процесса в САО с запоминанием максимальной скорости выходного параметра приведен в приложении А.

8. ИССЛЕДОВАНИЕ ВЛИЯНИЯ ПАРАМЕТРОВ ДИНАМИЧЕСКОЙ НАСТРОЙКИ САО НА ПОКАЗАТЕЛИ КАЧЕСТВА ПЕРЕХОДНОГО ПРОЦЕССА

САО с запоминанием максимума скорости выходной величины и с остановкой ИМ в момент достижения максимума для управления энергетическим режимом имеет два параметра динамической настройки:

- Скорость изменения входной величины _ К_и, (кА) /с;

- Зона нечувствительности СЭР - ДZ_н, Вт/c.

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

- Потеря на поиск - Z_п, Вт;

- Время выхода СЭР на экстремум (время поиска) - Т_п, с;

- Размах поисковых колебаний - А_х, кА.

Величина потери на поиск Z_п определяется как разность между максимальным значением выходного параметра в установившемся режиме Y_max и средним значением выходного параметра Z(x) в пределах одного предельного цикла [12]. Графики исследования переходных процессов, фазовые портреты и графики зависимости показателей качества от величины скорости ИМ и от зоны нечувствительности представлены ниже.

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

8.1 Исследование переходных процессов САО при отсутствии дрейфа статической характеристики

Определим оптимальную скорость изменения входной величины К_и при постоянной зоне нечувствительности ДZ_н = 0,075 Вт/c. Расчетные траектории поискового процесса СЭР для зоны нечувствительности ДZ_н=0,075Вт/c и для разных значений скорости изменения входной величины К_и приведены на рисунках 8.1 - 8.9.

Рисунок 8.1 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при ДZ_н = 0,075Вт/с и К_и = 0,33•10³ кА/с

Рисунок 8.2 - Фазовый портрет САО при ДZ_н = 0,05Вт/c и К_и = 0,1•10³ кА/с



Рисунок 8.3 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при ДZ_н = 0,075Вт/с и К_и = 0,36•10³ кА/с

Рисунок 8.4 - Фазовый портрет САО при ДZ_н = 0,075Вт/c и К_и = 0,36•10³ кА/с



Рисунок 8.5 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при ДZ_н = 0,075Вт/с и К_и = 0,385•10³ кА/с

Рисунок 8.6 - Фазовый портрет САО при ДZ_н = 0,075Вт/c и К_и =0,385•10³кА/с



Рисунок 8.7 - Увеличенная часть фазового портрета при ДZ_н = 0,075Вт/c и К_и = 0,33•10³кА/с

Рисунок 8.8 - Увеличенная часть фазового портрета при ДZ_н = 0,075Вт/c и К_и = 0,36•10³кА/с



Рисунок 8.9 - Увеличенная часть фазового портрета при ДZ_н = 0,075Вт/c и К_и=0,385•10³кА/с

Для более точного определения показателей Z_п и A_x приведены графики увеличенной верхней части фазового портрета. Сведем показатели качества работы СЭР для ДZ_н=0,075Вт/c и различных значений К_и при постоянной зоне нечувствительности в таблицу 8.1.

Таблица 8.1 - Зависимость показателей качества от изменения скорости ИМ

Ки	ДZн	Тп	Zп	Ах
103кА/с	Вт/c	c	Вт	103•кА
0,33	0,075	595	1,1	6,9
0,36	0,075	545	1,1	7,1
0,385	0,075	577	2,1	8,1

Из рисунка 8.10 видно, что при увеличении скорости ИМ время выхода на экстремум Т_п и величина потери на поиск Z_n уменьшается, а затем увеличивается, это объясняется тем, что при значительной скорости выходная величина в области экстремума прежде чем остановиться делает три реверса, а при меньшей скорости _ два.



Рисунок 8.10 - Зависимость показателей качества работы САО Z_n, Т_п и A_x от скорости изменения входного параметра К_и

Для зоны нечувствительности ДZ_н = 0,075Вт/c оптимальной скоростью ИМ является скорость К_и=0,36•10³кА/с.

Теперь определим влияние зоны нечувствительности на показатели качества переходного процесса и покажем на рисунках 8.11 - 8.14.

Рисунок 8.11 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при ДZ_н = 0,11Вт/с и К_и = 0,36•10³кА/с



Рисунок 8.12 - Увеличенная часть фазового портрета при ДZ_н = 0,11Вт/c и К_и = 0,36•10³кА/с

Рисунок 8.13 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при ДZ_н = 0,05Вт/с и К_и = 0,36•10³ кА/с



Рисунок 8.14 - Увеличенная часть фазового портрета при ДZ_н = 0,05°C/c и К_и = 0,36•10³кА/с

Из графиков, см. рисунки 8.11 - 8.14 определяем показатели качества работы системы и сводим их в таблицу 8.2. График зависимости показателей качества от зоны нечувствительности САО при постоянной скорости ИМ представлен на рисунке 8.15.

Таблица 8.2 - Зависимость показателей качества от изменения зоны нечувствительности при постоянной скорости ИМ

Ки	ДZн	Тп	Zп	Ах
103 кА/с	Вт/c	c	Вт	103кА
0,36	0,05	522	0,8	6,1
0,36	0,075	545	1,1	7,1
0,36	0,11	563	1,7	8,2

По зависимости показателей качества от величины зоны нечувствительности, изображенной на рисунке 8.15 видно, что при уменьшении зоны нечувствительности время выхода на экстремум Т_п, величина потери на поиск Z_n и размах поисковых колебаний А_х уменьшаются. Однако уменьшение зоны нечувствительности приводит к уменьшению диапазона поисковых колебаний входного параметра и к уменьшению времени выхода на максимум.

Рисунок 8.15 - Зависимость показателей качества от величины зоны нечувствительности ДZ_н

В результате исследования переходных процессов в САО при различных значениях зоны нечувствительности и различных значениях скорости ИМ были получены оптимальные параметры настройки:

ДZ_н =0,05 Вт/c;

К_и=0,36•10³(кА)/с.

Переходный процесс при оптимальных параметрах настройки представлен на схеме Д.А.220200.005.БР.10.ВО рисунок 2.

Фазовый портрет при оптимальных параметрах настройки представлен на схеме Д.А.220200.005.БР.10.ВО рисунок 3.

Полученные результаты характеризуют эффективность работы САО в условиях отсутствия дрейфа статической характеристики оптимизируемого процесса. Однако для окончательного выбора параметров настройки САО необходимо исследовать поведение системы в условиях дрейфа статической характеристики оптимизируемого процесса.

8.2 Исследование переходных процессов САО при наличии дрейфа статической характеристики

В реальных производственных условиях наличие различных технологических возмущений приводит к смещению статической характеристики оптимизируемого процесса, что оказывает заметное влияние на поисковый режим[13].

Ниже на рисунках приведены расчетные траектории изменения Z(t) X(t) и Y(t) при наличии сложного дрейфа статической характеристики для оптимальных параметров САО: ДZ_н=0,05Вт/c К_и=0,36•10³(кА)/с.

Дрейф статической характеристики оптимизируемого процесса вправо-вниз соответствует тепловому режиму рабочего пространства при повышении производительности нагревательной печи и изображен на рисунках 8.16. и 8.17.

Рисунок 8.16 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при дрейфе статической характеристики вправо - вниз



Рисунок 8.17 - Фазовый портрет САО при дрейфе статической характеристики вправо - вниз

Дрейф статической характеристики оптимизируемого процесса вправо-вверх, соответствует энергетическому режиму дуговой печи при увеличении рабочего тока, изображенный на рисунках 8.18 и 8.19.

Рисунок 8.18 - Расчётные траектории изменения во времени входного Х(ф), выходного Z(ф) и установившегося Y(ф) значений параметров САО при дрейфе статической характеристики вправо - вверх



Рисунок 8.19 - Фазовый портрет САО при дрейфе статической характеристики вправо - вверх

Из представленных графиков видно, что выбранные ранее настройки системы экстремального регулирования позволяют с высокой точностью отслеживать экстремум выходной величины в условиях сложного монотонного дрейфа статической характеристики. Это доказывает высокую эффективность работы САО, в условиях, близких к реальным производственным[14].

ЗАКЛЮЧЕНИЕ

Основываясь на исследовании переходных процессов в САО, с запоминанием максимума выходного параметра, см. рисунки 8.1 - 8.15, при переборе различных значений зоны нечувствительности ДZ_н и скорости К_и ИМ, можно сделать следующие выводы:

Каждому технологическому процессу соответствует определенная скорость изменения входного управляющего воздействия, при которой осуществляется максимальное приближение момента остановки исполнительного механизма к моменту достижения оптимального значения выходного параметра, для рассматриваемого процесса такая скорость К_и = 0,36•10³(кА)/с, любое отклонение значения этой скорости неблагоприятно влияет на показатели качества работы САО;

Для уменьшения величины потерь на поиск (?Z_п), и размаха поисковых колебаний (А_х) можно уменьшить зону нечувствительности (?Z_н), но чрезмерное уменьшение величины зоны нечувствительности приведет к необходимости приобретения более дорогостоящих приборов, а также к чрезмерной чувствительности системы (ложным срабатываниям на неконтролируемые воздействия) и, как следствие, ухудшению показателей качества.

Исследование системы в условиях сложного монотонного дрейфа, см. рисунки 8.16 - 8.19 доказывает, что выбранные оптимальные параметры настройки позволяют САУ с высокой точностью отслеживать экстремум выходной величины в условиях, приближенных к производственным.

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1 Электрические промышленные печи: Дуговые печи и установки специального нагрева: Учебник для вузов/ А. Д. Свенчанский и др.; Под ред. А. Д. Свенчанского. - 2-е изд., перераб. и доп. - М.: Энергоиздат, 2001. - 296 с., ил.

2 Ишметьев Е. Н. и др. Автоматизация и оптимизация управления технологическими процессами внепечной доводки стали: Монография. - Магнитогорск: ГОУ ВПО "МГТУ", 2008. - 312 с.

3 Автоматическое управление электротермическими установками: Учебник для вузов/ А.М. Кручинин и др.; Под ред. А. Д. Свенчанского. - М.: Энергоатомиздат, 2000. - 416 с.: ил.

4 Л. В. Лапшин. Автоматизация технологических процессов дуговой сталеплавильной печи. - М.: ООО "Квадратум", 2002. - 157 с.

5 Окороков Н. В. Электроплавильные печи черной металлургии. - М.: Металлургия, 2005. - 220 с.

6 Н. А. Марков. Электрические печи и режимы дуговых электропечных установок. - М.: Энергия, 2003. - 204 с.

7 Б. Н. Парсункин, М. В. Бушманова. Расчет переходных процессов в системах экстремального регулирования с запоминанием экстремума: Учебное пособие. - Магнитогорск: МГТУ, 2001. - 164 с.

8 Электрические печи сопротивления и дуговые печи./ Под ред. Гутмана М. Б. - М.: Энергоатомиздат, 2003. - 320 с.

9 Глинков Г. М., Маковский В. А. АСУ ТП в агломерационных и сталеплавильных печах. - М.: Металлургия, 2002. - 360 с.

10 Промышленные установки электродугового нагрева и их параметры./ Под ред. Никольского Л. Е. - М.: Энергия, 2001. - 270 с.

11 Лапшин И. В. Автоматизация дуговых печей. Учебник для вузов. - М: Издательство "МИСиС", 2004. - 165 с.

12 Бигеев А. М., Бигеев В. А. Металлургия стали. Теория и технология плавки стали. Учебник для вузов, 3-е изд. перераб. и доп. - Магнитогорск: МГТУ, 2000. - 544 с.

13 Ефроймович Ю. Е. Оптимальные электрические режимы дуговых сталеплавильных печей. - М: Металлургиздат, 2005. - 98 с.

14 Идентификация элементов систем управления и оптимизации контуров управления технологическими процессами./ Парсункин Б. Н. - Магнитогорск, 2006. - 148 с.

Приложение

Текст рабочей программы расчета переходного процесса в САО с запоминанием максимальной скорости выходного параметра

program max_remember;

uses crt,graph;

const

a0=1210326.5354538; {коэффициенты статической характеристики}

a1=574.107323836;

a2=0.005478383;

a3=-0.000000096;

a4=0;

T1=0,05; {постоянная времени объекта }

T2=0,02; {время запаздывания }

ku=412.5; {скорость исполнительного механизма }

zn=0.1; {зона нечувствительности }

n=1500; {количество расчетных точек }

dt=1; {период дискретизации }

xn=12;

xk=42;

at=0; {скорость горизонтального дрейфа }

bt=0; {скорость вертикального дрейфа }

var

i,j,gd,ga,x_gr,z_gr,y_gr,mx,my,tk,tv,dreif:integer;

gr:string;

dzmax,U, sigma, sigma0,ht,hte,hx,maxx,maxy, minx, x1, x2:real;

miny,maxz,minz:real;

x,y,t,z1,dz1,z,dz: array [0..n] of real;

function f(x,t:real):real;

begin

f:=(a0+a1*(x+at*t)+a2*(x+at*t)*(x+at*t)+a3*(x+at*t)*(x+at*t)

*(x+at*t)+a4*(x+at*t)*(x+at*t)*(x+at*t)*(x+at*t))+bt*t;

end;

{расчет изменения во времени величин входного X(ф), }

{выходного Y(ф) и установившегося Z(ф) значений }

Begin

t[0]:=0; {ввод начальных условий}

x[0]:=xn;

y[0]:=f(x[0],t[0]);

z1[0]:=f(x[0],t[0]);

z[0]:=f(x[0],t[0]);

dz1[0]:=0;

dz[0]:=0;

dzmax:=0; U:=1;

sigma:=1 ;

For i:=1 to n do

BEGIN

t[i]:=t[i-1]+dt;

begin

x[i]:=x[i-1]+ku*sigma*dt;

if (x[i]<= 10) or (x[i]>=38) then break; {концевые ограничители}

y[i]:=f(x[i],t[i]);

z1[i]:=dz1[i-1]+z1[i-1];

dz1[i]:=(1/T1)*(y[i]-z1[i]);

z[i]:=dz[i-1]+z[i-1];

dz[i]:=(1/T2)*(z1[i]-z[i]);

if dz[i]>dzmax then dzmax:=dz[i];

if dz[i]>=0 then

begin

if dz[i]-dzmax+zn>0 then

begin U:=1; tv:=0; end;

if dz[i]-dzmax+zn<=0 then

begin U:=0;

end;

end;

if dz[i]<0 then

if (dz[i]+zn<0) then

begin

U:=(-1); dzmax:=0; tv:=tv+1; end;

if U=1 then

begin

sigma0:=sigma;

sigma:=sigma;

end;

if U=0 then

sigma:=0;

if U=(-1)then

begin

sigma:=sigma0*(-1); end;

if tv<=0.2*T1 then sigma:=sigma;

if (U=-1) and (sigma=1) then x1 :=x[i];

if (U=-1) and (sigma=-1) then x2:=x[i];

if x[i]=x1+(x2-x1)/2 then

begin

sigma:=0; x1:=0;

x2:=0; end;

if sigma=0 then tk:=tk+1; if (tk>3*T1) then

begin

sigma:=-1 ; tk:=0;

dzmax:=dz[i]; end; end;

write('t=',t[i]:3:0, ' y=',y[i]:5:2,' x=',x[i]:2:0,' z=',z[i]:3:2,' sigma=',sigma:1:0,' dz[i]= ',dz[i]:1:3,' ');

writeln(' dzmax=', dzmax:2:3);

END; readkey;

{Построение переходного процесса}

detectgraph(gd,ga);

initgraph(gd,ga,'c:\bp\bgi');

setcolor(15);

minx:=0;

maxx:=n ;

miny:=1250;

maxy:=1400;

mx:=getmaxx;

my:=getmaxy;

maxz:=45;

minz:=15;

x_gr:=15;

y_gr:=5;

z_gr:=10;

ht:=(mx-100)/(maxx-minx);

hte:=(my-100)/(maxy-miny);

hx:=(my-100)/(maxz-minz);

line(50,my-50,mx-50,my-50);

line(50,50,50,my-50);

line(mx-50,50,mx-50,my-50);

outtextxy(200,0,'perehod process');

outtextxy(460,my-10,'vremya v sekundah');

settextstyle(0,1,0);

outtextxy(10,110,'temperatura');

outtextxy(570,110,'rashod vozduha');

settextstyle(0,0,0);

setlinestyle(0,1,1);

for i:=0 to x_gr do

begin

line(round(i*(mx-100)/x_gr)+50,425,round(i*(mx-100)/x_gr)+50,my-45);

str(minx+i*(maxx-minx)/x_gr:5:0,gr);

outtextxy(round(i*(mx-100)/x_gr),my-30,gr); end;

for i:=0 to y_gr do

begin

line(mx-585,(my-50)-round(i*(my-100)/y_gr),45,(my-50)-round(i*(my-100)/y_gr));

str(miny+i*(maxy-miny)/y_gr:4:0,gr);

outtextxy(13,(my-50)-round(i*(my-100)/y_gr),gr);

end;

for i:=0 to z_gr do

begin

line(mx-46,(my-50)-round(i*(my-100)/z_gr),585,(my-50)-round(i*(my-100)/z_gr));

str(minz+i*(maxz-minz)/z_gr:5:0,gr); outtextxy(600, (my-50)-round(i*(my-100)/z_gr),gr);

end;

t[0]:=0;

x[0]:=xn;

y[0]:=f(x[0],t[0]);

z1[0]:=f(x[0],t[0]);

z[0]:=f(x[0],t[0]);

dz1[0]:=0;

dz[0]:=0;

dzmax:=0; U:=1; sigma:=1;

setlinestyle(0,1,1); moveto (50,getmaxy-50);

For i:=1 to n do

BEGIN

t[i]:=t[i-1]+dt; begin

x[i]:=x[i-1]+ku*sigma*dt;

if (x[i]<=9) or (x[i]>=40) then break;

y[i]:=f(x[i],t[i]);

z1[i]:=dz1[i-1]+z1[i-1];

dz1[i]:=(1/T1)*(y[i]-z1[i]);

z[i]:=dz[i-1]+z[i-1];

dz[i]:=(1/T2)*(z1[i]-z[i]);

if dz[i]>dzmax then dzmax:=dz[i]; if dz[i]>=0 then

begin

if dz[i]-dzmax+zn>0 then

begin

U:=1;

tv:=0;

end;

if dz[i]-dzmax+zn<=0 then

begin

U:=0;

end;

end;

if dz[i]<0 then

if (dz[i]+zn<0) then

begin

U:=(-1);

dzmax:=0;

tv:=tv+1;

end;

if U=1 then

begin

sigma0:=sigma; sigma:=sigma;

end;

if U=0 then sigma:=0 ; if U=(-1) then

begin

sigma:=sigma0*(-1);

end;

if tv<=0.1*T1 then

sigma:=sigma;

if (U=-1) and (sigma=1) then

x1 :=x[i]; if (U=-1) and (sigma=-1)

then x2:=x[i];

if x[i]=x1+(x2-x1)/2then begin

sigma:=0; x1:=0; x2:=0; end;

if sigma=0 then tk:=tk+1; if (tk>3*T1)then

begin

sigma:=-1;

tk:=0;

dzmax:=dz[i];

end;

end;

setcolor(15);

line(trunc((t[i-1])*ht/dt)+50,getmaxy-50-trunc((f(x[i-1],t[i-1])-miny)*hte),

trunc(t[i]*ht/dt)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

setcolor(14);

line(trunc((t[i-1])*ht/dt)+50,getmaxy-50-trunc((z[i-1]-miny)*hte), trunc((t[i])*ht/dt)+50,getmaxy-50-trunc((z[i]-miny)*hte));

setcolor(15);

line(trunc((t[i-1])*ht/dt)+50,getmaxy-50-trunc((x[i-1]-minz)*hx),trunc((t[i])*ht/dt)+50,getmaxy-50-trunc((x[i]-minz)*hx));

END;

readkey;

begin

detectgraph(gd,ga);

initgraph(gd,ga,'c:\bp\bgi');

setcolor(15);

minx:=10;

maxx:=42;

miny:=1260;

maxy:=1380;

mx:=getmaxx;

my:=getmaxy;

x_gr:=8;

y_gr:=6;

hte:=(my-100)/(maxy-miny); {построение осей}

hx:=(mx-100)/(maxx-minx);

line(50,my-50,mx-50,my-50);

line(50,50 ,50,my-50);

outtextxy(200,0,'fazovy portret');

outtextxy(460,my-10,'rashod vozduha');

settextstyle(0,1,0);

outtextxy(10,100,'temperatura');

settextstyle(0,0,0);

setlinestyle(0,1,1);

for i:=0 to x_gr do

begin

line(round(i*(mx-100)/x_gr)+50,425,round(i*(mx-100)/x_gr)+50,my-45);

str(minx+i*(maxx-minx)/x_gr:5:0,gr);

outtextxy(round(i*(mx-100)/x_gr),my-30,gr);

end;

for i:=0 to y_gr do

begin

line(mx-585,(my-50)-round(i*(my-100)/y_gr),45,(my-50)-round(i*(my-100)/y_gr));

str(miny+i*(maxy-miny)/y_gr:4:0,gr);

outtextxy(15,(my-50)-round(i*(my-100)/y_gr),gr);

end;

t[0]:=0; {ввод начальных условий}

x[0]:=xn;

y[0]:=f(x[0],t[0]);

z1[0]:=f(x[0],t[0]);

z[0]:=f(x[0],t[0]);

dz1[0]:=0; dz[0]:=0;

dzmax:=0; U:=1; sigma:=1;

setlinestyle(0,1,1); moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=t[i-1]+dt;

begin

x[i]:=x[i-1]+ku*sigma*dt;

if (x[i]<=9) or (x[i]>=40) then break;

y[i]:=f(x[i],t[i]);

z1[i]:=dz1[i-1]+z1[i-1];

dz1[i]:=(1/T1)*(y[i]-z1[i]);

z[i]:=dz[i-1]+z[i-1];

dz[i]:=(1/T2)*(z1[i]-z[i]);

if dz[i]>dzmax then dzmax:=dz [i];

if dz[i]>=0 then

begin

if dz[i]-dzmax+zn>0 then

begin

U:=1;

tv:=0;

end;

if dz[i]-dzmax+zn<=0 then

begin

U:=0;

end;

end;

if dz[i]<0 then

if (dz[i]+zn<0) then

begin U:=(-1);

dzmax:=0;

tv:=tv+1;

end;

if U=1 then begin

sigma0:=sigma; sigma:=sigma; end;

if U=0 then

sigma:=0;

if U=(-1) then

begin

sigma:=sigma0*(-1);

end;

if tv<=0.2*T1 then sigma:=sigma;

if (U=-1) and (sigma=1) then x1 :=x[i];

if (U=-1) and (sigma=-1) then x2:=x[i];

if x[i]=x1+(x2-x1)/2 then

begin

sigma:=0;

x1:=0;

x2:=0;

end;

if sigma=0 then

tk:=tk+1;

if (tk>3*T1) then

begin

sigma:=-1; tk:=0;

dzmax:=dz[i];

end;

end;

setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((z[i-1]-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((z[i]-miny)*hte));

{setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1], t[i-1])-miny)*hte);

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));}

END;

x[0]:=xn; sigma:=1;

moveto (50,my-50); For i:=1 to n do

BEGIN t[i]:=0;

x[i]:=x[i-1]+ku*sigma*dt;

if x[i]>=xk then sigma:=0; setcolor(15);

line (trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

x[0]:=9; sigma:=1;

moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=300;

x[i]:=x[i-1]+ku*sigma*dt;

if x[i]>=40 then sigma:=0; setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

x[0]:=10;

sigma:=1;

moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=600;

x[i]:=x[i-1]+ku*sigma*dt; if x[i]>=10 then sigma:=0; setcolor(15);

line (trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

readkey;

end;

begin

detectgraph(gd,ga);

initgraph(gd,ga, 'c:\bp\bgi');

setcolor(15);

minx:=23;

maxx:=38;

miny:=1360;

maxy:=1375;

mx:=getmaxx;

my:=getmaxy;

x_gr:=10;

y_gr:=6;

hte:=(my-100)/(maxy-miny); {построение осей}

hx:=(mx-100)/(maxx-minx);

line(50,my-50,mx-50,my-50); line(50,50,50,my-50);

outtextxy(50,0,'uvelich fozoviy portret');

outtextxy(460,my-10,'rashod vozduha');

settextstyle(0,1,0);

outtextxy(10,100,'temperatura');

settextstyle(0,0,0);

setlinestyle(0,1,1);

for i:=0 to x_gr do begin

line(round(i*(mx-100)/x_gr)+50,425,round(i*(mx-100)/x_gr)+50,my-45);

str(minx+i*(maxx-minx)/x_gr:5:0,gr);

outtextxy(round(i*(mx-100)/x_gr),my-30,gr); end;

for i:=0 to y_gr do begin

line(mx-585,(my-50)-round(i*(my-100)/y_gr),45,(my-50)-round(i*(my -100)/y_gr));

str(miny+i*(maxy-miny)/y_gr:4:0,gr);

outtextxy(15,(my-50)-round(i*(my-100)/y_gr),gr);

end;

t[0]:=0; {ввод начальных условий}

x[0]:=xn;

y[0]:=f(x[0],t[0]);

z1[0]:=f(x[0],t[0]);

z[0]:=f(x[0],t[0]);

dz1[0]:=0;

,dz[0]:=0;

dzmax:=0;

U:=1;

sigma:=1;

setlinestyle(0,1,1);

moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=t[i-1]+dt;

begin

x[i]:=x[i-1]+ku*sigma*dt;

if(x[i]<=5) or (x[i]>=40) then break;

y[i]:=f(x[i],t[i]);

z1[i]:=dz1[i-1]+z1[i-1];

dz1[i]:=(1/T1)*(y[i]-z1[i]);

z[i]:=dz[i-1]+z[i-1];

dz[i]:=(1/T2)*(z1[i]-z[i]);

if dz[i]>dzmax then dzmax:=dz[i];

if dz[i]>=0 then

begin

if dz[i]-dzmax+zn>0 then begin U:=1; tv:=0;

end;

if dz[i]-dzmax+zn<=0 then

begin

U:=0;

end;

end;

if dz[i]<0 then

if (dz[i]+zn<0) then

begin

U:=(-1);

dzmax:=0; tv:=tv+1;

end;

if U=1 then

begin

sigma0:=sigma;

sigma:=sigma;

end;

if U=0 then

sigma:=0; if U=(-1) then

begin

sigma:=sigma0*(-1);

end;

if tv<=0.2*T1 then sigma:=sigma;

if (U=-1) and (sigma=1) then x1:=x[i];

if (U=-1) and (sigma=-1) then x2:=x[i];

if x[i]=x1+(x2-x1)/2 then

begin

sigma:=0; x1:=0;

x2:=0;

end;

if sigma=0 then

tk:=tk+1;

if (tk>3*T1) then

begin

sigma:=-1; tk:=0;

dzmax:=dz[i];

end;

end;

setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((z[i-1]-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((z[i]-miny)*hte));

setcolor(11);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i-1])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

x[0]:=xn; sigma:=1;

moveto (50,my-50);

for i:=1 to n do

BEGIN t[i]:=0;

x[i]:=x[i-1]+ku*sigma*dt; if x[i]>=xk then sigma:=0; setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

x[0]:=10; sigma:=1;

moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=300;

x[i]:=x[i-1]+ku*sigma*dt; if x[i]>=40 then sigma:=0; setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

x[0]:=10; sigma:=1;

moveto (50,my-50);

for i:=1 to n do

BEGIN

t[i]:=600;

x[i]:=x[i-1]+ku*sigma*dt; if x[i]>=10 then sigma:=0; setcolor(15);

line(trunc((x[i-1]-minx)*hx)+50,getmaxy-50-trunc((f(x[i-1],t[i])-miny)*hte),

trunc((x[i]-minx)*hx)+50,getmaxy-50-trunc((f(x[i],t[i])-miny)*hte));

END;

readkey;

end;

End.

Размещено на Allbest.ru