Регулируемые параметры ГТУПолученные
характеристики ГТУ
|
|
|
1143.85428096 кВт
|
|
3065.16863501 кВт
|
|
3208.02703281 кВт
|
|
1000.00099588 кВт
|
|
0.326820272955
|
|
0.506133284302
|
|
736.703802982 К
|
|
1010.27777001 К
|
Данное исследование позволяет теоретически рассчитать ГТУ
необходимой нам мощности, при различных регулируемых начальных параметрах. На
практике регулируемые начальные параметры определяются исходя из технико-экономического
обоснования, характеристик имеющегося оборудования, а также характера
производственного процесса, в котором данная ГТУ будет применена.
3.4
Исследование влияния температур газа перед турбинами на КПД ГТУ
Исследовать влияние температуры газа перед турбинами на КПД ГТУ. Сравнение проведем
при заданном уровне мощности ГТУ.
Задаются ,, , , и два значения температуры (например, 650°С и 850
°С). Требуется найти параметры и характеристики цикла, соответствующие
максимуму КПД.
Для решения этой задачи так же, как и в исследовании 3.3
настоящей работы, определяем G, , , , для одного значения. После этого находим эти
же параметры и характеристики для другого значения и результаты сравниваем.
Сравниваются расходы воздуха и КПД. В ТS-диаграмме строим циклы
этих двух ГТУ. Недостающие для вычерчивания циклов данные рассчитываются по
(2.11) и (2.23).
Вводим начальные данные (стандартные регулируемые параметры)
в расчетную программу, работающую по алгоритму, представленному на рис.3.8:
Температура t3 = 650°С:
Макс КПД =0.249290676449 при p2=1.0 и p5= 0.520001
Расход G= 6.3445 кг\с
Температура t3 = 800°С:
Макс КПД =0.302055312208 при p2=1.4 и p5=0.700001
Расход G= 7.5291 кг\с
Анализируя полученные зависимости, видим, что при повышении
температуры после камеры сгорания в ГТУ c двухступенчатым
расширением рабочего тела в турбине, увеличивается её КПД на
(0.3021-0.249)/ 0.3021=17,58%.
Также приходится увеличивать давления p2 и p5, чтобы мощность всей
установки оставалось оптимальной.
Однако и расход рабочего агента увеличивается на 1.18 кг\с,
что также влияет на экономическое обоснование данной схемы.
Заключение
В ходе работы, используя математический расчет
термодинамических свойств рабочего тела (воздуха) было произведено
моделирование простой газотурбинной установки, а также ГТУ с одноступенчатым
сжатием в компрессоре и двухступенчатым расширением рабочего тела в турбине.
Это позволило наглядно выбирать оптимальное сочетание параметров ГТУ, а также
производить регулирование её КПД и удельной мощности и определять их
максимальное значение. Были получены следующие зависимости:
. Параболическая зависимость с четко выраженным
максимумом мощности и КПД ГТУ от роста давления в камере сгорания.
. Увеличение мощности и КПД ГТУ при увеличении
температуры в камере сгорания (уменьшении температуры окружающей среды), причем
возрастает давление в камере сгорания для их оптимальных значений.
. Увеличение мощности и КПД ГТУ при увеличении
внутренних относительных КПД турбины и компрессора.
Разработанная программа легка в использовании, листинги
команд разбиты на отдельные функции для каждой задачи.
Получены навыки программирования в среде разработки Python 2.7, а также навыки
графического представления полученных результатов (pyGlade и pyQT).
Библиографический список
1) Ривкин
С.Л. Термодинамические свойства газов. - 3-е изд. - М.: Энергия, 1973 - 288с.
2) Костюк
А.Г., Шерстюк А.Н. Газотурбинные установки. - М.: Высшая школа, 1979. - 254 с.
) Цанев
С.В., Буров В.Д., Ремезов А.Н. Газотурбинные и парогазовые установки тепловых
электростанций. - М.: МЭИ, 2002. - 584 с.
) Сычев
В.В., Вассерман А.А. Термодинамические свойства воздуха. - М.: Издательство
стандартов, 1978. - 276 с.
Приложение
Исходный код программы расчетов простой ГТУ
# -*- coding: cp1251 -*-
#Массивы с константами_c = (29.438265, -1.610822, -11.99174,
68.828384, -98.23993,
64.88351, -20.9094, 2.66524)_h = (-542,
29438.265, -805.411, -3997.2481, 17207.1, -19647.99, 10813.92, -2987.054,
.155)_s = (230.1763, -1.610822, -5.995872,
22.942794, -24.55998, 12.9767, -3.4849, 0.380749)_s = 29.43821
#Считаем энтальпиюmathpylabmatplotlib
import mlabental(x):= 0= float(s)i in xrange(len(a_h)):+= a_h[i]*((x/1000)**i)=
(s+488)/28.97r
#Считаем Энтропию
s = 0= float(s)i in xrange(len(a_s)):+=
a_s[i]*((x/1000)**i)(s + b_s*math.log(x/1000))/28.97
#Cчитаем теплоемкостьteploem(x):
s = 0= float(s)i in xrange(len(a_c)):+=
a_c[i]*((x/1000)**i)= s/28.97r
#Вспомогательная функция энтропии для построения TS диаграмм
def s_gr(x,y):= 0= float(s)=
1*math.log(x/273.15)-0.287*math.log(y/0.1)
return s
#Основная функция рассчитывающая параметры ГТУ
def osnova(x1,x2,x3,x4,x5,x6,x7,x8):= x2+273.15=
x3+273.15= X3/X2= (x7/x1)**((x8-1)/x8)= (x4*teploem(X2)*X2*tetta*(tau-1)*x5)/tau=
(x4*teploem(X2)*X2*(tau-1))/x6= Nt - Nk=
x4*teploem(X2)*X2*(tau-1)*((tetta/tau)-1)= (tau-1)/tau=
((tau-1)*((tetta*x5/tau)-(1/x6)))/(tetta-1-((tau-1)/x6))d =
X3*(1-((tau-1)*x5)/tau)d = X2*(1+(tau-1)/x6)(X2, X3, tetta, tau, Nt, Nk, Ngtu,
Ngtut, KPDt, KPDi, T4d, T2d)
#Вспомогательная функция для корректной печатиprent(x):
print " ",round(x[0], 3),"
",round(x[1], 3)," ",round(x[2], 3)," ",round(x[3],
3)," ",round(x[4], 3), \
" ",round(x[5], 3),"
",round(x[6], 3)," ", \(x[7], 3)," ",round(x[8],
3)," ",round(x[9], 3)," ",round(x[10], 3),"
",round(x[11], 3)
#Решение первой задачиzadacha1():=p1 #переменная х будет
содержать в себе изменяющиеся значения p2 начиная с p1= 0 #M=параметр
MM=значение p2. Содержат в себе координаты максимумов= 0
M3 = 0= 0= 0= 0x<3:" p2=",x=
osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)
if M1 < os1[6]: #Высчитывается точка максимума Nгту=
os1[6]= xM2 < os1[4]: #Высчитывается точка максимума Nт= os1[4]= xM3 <
os1[9]: #Высчитывается точка максимума относительного внутреннего кпд KPDi
M3 = os1[9]= x"# T1 T3 tetta tau nt nk ngtu
ngtut kpdt kpdi t4d t2d"(os1)+= 0.2" ""Max znach Ngtu=
",M1," pri davlenii P2=",MM1"Max znach Nt=",M2,"
pri davlenii P2=",MM2"Max znach KPDi =",M3," pri davlenii
P2=",MM3
print " "
#Рисуем графики: 3 графика мощностей+2графика кпд + 2 TS
диаграммы= p1= 3= 0.1
# -------------------Рисуем график Ngtu
xlist = mlab.frange (xmin, xmax, dx)=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)[6] for x in xlist]= mlab.frange (0, MM1,
dx)= [M1 for x in xlist4]= mlab.frange (0, M1, dx)= [MM1 for y in ylist5].plot
(xlist, ylist, color= 'green').plot (xlist4, ylist4, color= 'black').plot
(xlist5, ylist5, color = 'black').xlabel('Davlenie P2
MPa').title('Zavisimost1').ylabel('Ngtu-green Nt-blue Nk-red')
# -------------------Рисуем график Nt =
mlab.frange (xmin, xmax, dx)= [osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)[4] for x in
xlist2].plot (xlist2, ylist2)
# -------------------Рисуем график Nk=
mlab.frange (xmin, xmax, dx)= [osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)[5] for x in
xlist2].plot (xlist6, ylist6, color= 'red').show()
# -------------------Рисуем график KPDi=
mlab.frange (xmin, xmax, dx)= [osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)[9] for x in
xlist3]= mlab.frange (0, MM3, dx)= [M3 for x in xlist8]= mlab.frange (0, M3,
0.001)= [MM3 for y in ylist9].plot (xlist3, ylist3, color= 'red').plot (xlist8,
ylist8, color= 'black').plot (xlist9, ylist9, color= 'black').xlabel('Davlenie
P2 MPa').ylabel('KPDi-red KPDt-blue').title('Zavisimost 2')
# -------------------Рисуем график KPDt= mlab.frange
(xmin, xmax, dx)= [osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)[8] for x in
xlist7].plot (xlist7, ylist7).show()
# -------------------Рисуем TS диаграмму для максимальной
Ngtu
pylab.xlabel('Entropya S').ylabel('Temperatura
T').title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= 1= [ymin,
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],MM1)].plot
(xlist10, ylist10, color= 'red')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],
0.01)= [s_gr(y,MM1) for y in ylist11].plot (xlist11, ylist11, color= 'red')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],MM1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10],p1)].plot
(xlist12, ylist12, color= 'red')= mlab.frange
(ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'red')
# -------------------Рисуем TS диаграмму для максимального
КПД
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= [ymin,
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11]]= [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],MM3)].plot
(xlist14, ylist14, color= 'blue')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],
0.01)= [s_gr(y,MM3) for y in ylist15].plot (xlist15, ylist15, color= 'blue')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],MM3),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10],p1)].plot
(xlist16, ylist16, color= 'blue')= mlab.frange
(ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist17].plot (xlist17, ylist17, color= 'blue').show()()
#решаем 2ю задачуzadacha2():= 600= float(-40) = 0= [] #Массив bsys[[],[],[]...]
содержит y-координаты 10 графиков Ngtu. bsysx[x1,x2,x3,] содержит х-координаты
этих графиков.= []= [] #Массив csys[[],[],[]...] содержит y-координаты 10
графиков KPDi. csysx[x1,x2,x3,] содержит х-координаты этих графиков.= []= []
#Массив Nm[y1,y2,y3...] содержит y-координаты максимумов Ngtu. Nmx[x1,x2,x3,]
содержит х-координаты этих максимумов.= []= [] #Массив Nm2[y1,y2,y3...]
содержит y-координаты максимумов KPDi. Nmx2[x1,x2,x3,] содержит х-координаты
этих максимумов.= [] "Budem menyat temperaturu t3:" #В первой части
задачи меняем температуру после КС t3
while y<=1500:= 0= 0= 0= 0= []=p1x<3:=
osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)M1 < os1[6]:= os1[6]= x M3 < os1[9]:=
os1[9]= x+= 0.2" t3=",y," Ngtu максимальная=",round(M1,3),"
при оптимальном
p2=",round(MM1,3)," Kpdi макс=",round(M3,3)," при оптимальном
p2=",round(MM3,3)
y += 100
# Для их построения найдем максимумы поточнее чем необходимо.
y = 600y<=1500:= 0= 0= 0= 0= []= []=p1
x<3:M1<osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[6]:
M1 = osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[6] #По пути запоминаем
точки максимумов
MM1 =
xM3<osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[9]:
M3 = osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[9] #По пути
запоминаем точки максимумов= x= a +[osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[6]] #В
массив а записываем y координаты мощностей а в b кпд при определенном t3
b = b +[osnova(p1,t1,y,G,KPDoit,KPDoik,x,k)[9]]+=
0.01= Nm + [M1]= Nmx +[MM1]= Nm2 + [M3]= Nmx2 +[MM3]
bsys = bsys + [a] #переносим массивы с координатами в ячейку
другого массива, меняем t3 и повторяем пока t3<=1500= csys + [b]+= 100
# Создаем массив с координатами х. количество элементов в
массиве должно быть равно количеству итераций в предыдущем цикле.СМ НА ШАГ!
x=0.1x<3: = bsysx + [x]
x +=0.01
#Строим графики Nгту в одном окне прокручивая цикл
i = 0i <int(len(bsys)):.plot (bsysx, bsys[i],
color= 'blue')+= 1
#соединяем точки максимумов .plot (Nmx,
Nm, color= 'red')
pylab.show()
#Строим графики KPDi в одном окне прокручивая цикл
i = 0i <int(len(csys)):.plot (bsysx, csys[i],
color= 'blue')+= 1.plot (Nmx2, Nm2, color= 'red').show()
#строим TS диаграммы(t3=800)
t3 = 800=p1 #переменная х будет содержать в себе изменяющиеся
значения p2 начиная с p1= 0 #M=параметр MM=значение p2. Содержат в себе
координаты максимумов= 0
MM1 = 0= 0x<3:=
osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)
if M1 < os1[6]: #Высчитывается точка максимума Nгту=
os1[6]= xM3 < os1[9]: #Высчитывается точка максимума относительного
внутреннего кпд KPDi= os1[9]= x+= 0.2
# -------------------Рисуем TS диаграмму для максимальной
Ngtu
pylab.xlabel('Entropya S').ylabel('Temperatura
T').title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= 1= [ymin,
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],MM1)].plot
(xlist10, ylist10, color= 'red')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],
0.01)= [s_gr(y,MM1) for y in ylist11].plot (xlist11, ylist11, color= 'red')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],MM1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10],p1)].plot
(xlist12, ylist12, color= 'red')= mlab.frange
(ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'red')
# -------------------Рисуем TS диаграмму для максимального
КПД
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= [ymin, osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],MM3)].plot
(xlist14, ylist14, color= 'blue')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],
0.01)= [s_gr(y,MM3) for y in ylist15].plot (xlist15, ylist15, color= 'blue')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],MM3),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10],p1)].plot
(xlist16, ylist16, color= 'blue')= mlab.frange
(ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist17].plot (xlist17, ylist17, color= 'blue').show()()zadacha2_2():=
float(-40)= 0= [] #Массив bsys[[],[],[]...] содержит y-координаты 10 графиков Ngtu.
bsysx[x1,x2,x3,] содержит х-координаты этих графиков.= []= [] #Массив
csys[[],[],[]...] содержит y-координаты 10 графиков KPDi.
csysx[x1,x2,x3,] содержит х-координаты этих графиков.= []= [] #Массив
Nm[y1,y2,y3...] содержит y-координаты максимумов Ngtu. Nmx[x1,x2,x3,] содержит х-координаты этих
максимумов.= []= [] #Массив Nm2[y1,y2,y3...] содержит y-координаты максимумов
KPDi. Nmx2[x1,x2,x3,] содержит х-координаты этих максимумов.= [] "Budem
menyat temperaturu t1:" #В этой части задачи меняем температуру наружного
воздуха t1
while y<=20:= 0= 0= 0= 0= []=p1x<3:=
osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)M1 < os1[6]:= os1[6]= x M3 < os1[9]:=
os1[9]= x+= 0.2" t1=",y," Ngtu максимальная=",round(M1,3),"
при оптимальном
p2=",round(MM1,3)," Kpdi макс=",round(M3,3)," при оптимальном
p2=",round(MM3,3)
y += 20
# Для их построения найдем максимумы поточнее чем необходимо.
y = float(-40)y<=20:= 0= 0= 0= 0= []= []=p1
x<3:M1<osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[6]:
M1 = osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[6] #По пути
запоминаем точки максимумов
MM1 =
xM3<osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[9]:
M3 = osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[9] #По пути
запоминаем точки максимумов= x= a +[osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[6]] #В
массив а записываем y координаты мощностей а в b кпд при определенном t3
b = b +[osnova(p1,y,t3,G,KPDoit,KPDoik,x,k)[9]]+=
0.01= Nm + [M1]= Nmx +[MM1]= Nm2 + [M3]= Nmx2 +[MM3]
bsys = bsys + [a] #переносим массивы с координатами в ячейку
другого массива, меняем t3 и повторяем пока t3<=1500= csys + [b]+= 20
# Создаем массив с координатами х. количество элементов в
массиве должно быть равно количеству итераций в предыдущем цикле.СМ НА ШАГ!
x=0.1x<3: = bsysx + [x]
x +=0.01
#Строим графики Nгту в одном окне прокручивая цикл
i = 0i <int(len(bsys)):.plot (bsysx, bsys[i],
color= 'blue')+= 1
#соединяем точки максимумов .plot (Nmx,
Nm, color= 'red')
pylab.show()
i = 0i <int(len(csys)):.plot (bsysx, csys[i],
color= 'blue')+= 1.plot (Nmx2, Nm2, color= 'red')
pylab.show()= 0=p1 #переменная х будет содержать в себе
изменяющиеся значения p2 начиная с p1= 0 #M=параметр MM=значение p2. Содержат в
себе координаты максимумов= 0
MM1 = 0= 0x<3:=
osnova(p1,t1,t3,G,KPDoit,KPDoik,x,k)
if M1 < os1[6]: #Высчитывается точка максимума Nгту=
os1[6]= xM3 < os1[9]: #Высчитывается точка максимума относительного
внутреннего кпд KPDi= os1[9]= x+= 0.2
# -------------------Рисуем TS диаграмму для максимальной
Ngtu
pylab.xlabel('Entropya S').ylabel('Temperatura
T').title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= 1= [ymin,
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],MM1)].plot
(xlist10, ylist10, color= 'red')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],
0.01)= [s_gr(y,MM1) for y in ylist11].plot (xlist11, ylist11, color= 'red')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[1],MM1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10],p1)].plot
(xlist12, ylist12, color= 'red')= mlab.frange (ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM1,k)[10],
0.01)= [s_gr(y,p1) for y in ylist13].plot (xlist13, ylist13, color= 'red')
# -------------------Рисуем TS диаграмму для максимального
КПД
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)[1]= [ymin, osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],MM3)].plot
(xlist14, ylist14, color= 'blue')= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[11],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],
0.01)= [s_gr(y,MM3) for y in ylist15].plot (xlist15, ylist15, color= 'blue')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[1],MM3),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10],p1)].plot
(xlist16, ylist16, color= 'blue')= mlab.frange
(ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM3,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist17].plot (xlist17, ylist17, color= 'blue').show()()zadacha3():
t3 = 400t3 <= 1000:= 1 # S будет хранить в себе максимум
Ngtu= 0.95 #kk=KPDoiK задаем через переменные так как придется менять внутри
функции а глобальные желательно не трогать
tt = kk #tt=KPDoiT"t3=",t3S>0:= 0=
0= p1
while x<3: #как обычно прокручиваем давление p2 от 0 до 3
Мпа.
if maxn < osnova(p1,t1,t3,G,tt,kk,x,k)[6]:=
osnova(p1,t1,t3,G,tt,kk,x,k)[6]= x+= 0.05"KPDoiK=KPDoiT=", kk,
"MaxN=", maxn, "при P2опт=", maxp= maxn= kk-0.05= kk+= 200= 400= 1=
[]= []t3 <= 1000:= []= []= 0.95 = 1S>0:= a+ [kk]= b +
[osnova(p1,t1,t3,G,kk,kk,p2,k)[6]]= osnova(p1,t1,t3,G,kk,kk,p2,k)[6]= kk-0.001=
ass + [a]= bss + [b]+= 200= 0i <int(len(bss)):.plot (ass[i], bss[i], color=
'blue')+= 1.show()= 800= 0.9= 1.1.xlabel('Entropya S').ylabel('Temperatura
T').title('TSdiagramma pri kpd=0.9-red ngtu=0-blue')= osnova(p1,t1,t3,G,kk,kk,p2,k)[0]=
osnova(p1,t1,t3,G,kk,kk,p2,k)[1]= 1= [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)].plot (xlist10,
ylist10, color= 'red')= mlab.frange
(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)=
[s_gr(y,p2) for y in ylist11].plot (xlist11, ylist11, color= 'red')=
[osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[1],p2),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[10],p1)].plot
(xlist12, ylist12, color= 'red')= mlab.frange
(ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'red')= 0.9 = 1S>0:= kk-0.001=
osnova(p1,t1,t3,G,kk,kk,p2,k)[6]= 800= osnova(p1,t1,t3,G,kk,kk,p2,k)[0]=
osnova(p1,t1,t3,G,kk,kk,p2,k)[1]= 1= [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)].plot (xlist10,
ylist10, color= 'blue')= mlab.frange
(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)=
[s_gr(y,p2) for y in ylist11].plot (xlist11, ylist11, color= 'blue')=
[osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[1],p2),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[10],p1)].plot
(xlist12, ylist12, color= 'blue')= mlab.frange
(ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'blue').show() ()zadacha3_2():= 400= 1=
[]= []t3 <= 1000:'t3=', t3= []= []= 0.95 = 1S>0:= a+ [kk]= b + [osnova(p1,t1,t3,G,kk,kk,p2,k)[9]]=
osnova(p1,t1,t3,G,kk,kk,p2,k)[9]' kk=', kk, ' kpd gtu = ', S= kk-0.001= ass +
[a]= bss + [b]+= 200= 0i <int(len(bss)):.plot (ass[i], bss[i], color=
'blue')+= 1.show()= 800= 0.9= 1.1.xlabel('Entropya S').ylabel('Temperatura T').title('TSdiagramma
pri kpd=0.9-red ngtu=0-blue')= osnova(p1,t1,t3,G,kk,kk,p2,k)[0]=
osnova(p1,t1,t3,G,kk,kk,p2,k)[1]= 1= [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)].plot (xlist10,
ylist10, color= 'red')= mlab.frange
(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)=
[s_gr(y,p2) for y in ylist11].plot (xlist11, ylist11, color= 'red')=
[osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[1],p2),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[10],p1)].plot
(xlist12, ylist12, color= 'red')= mlab.frange
(ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'red')= 0.9 = 1S>0:= kk-0.001= osnova(p1,t1,t3,G,kk,kk,p2,k)[9]=
osnova(p1,t1,t3,G,kk,kk,p2,k)[0]= osnova(p1,t1,t3,G,kk,kk,p2,k)[1]= 1= [ymin,
osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)].plot (xlist10,
ylist10, color= 'blue')= mlab.frange
(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)=
[s_gr(y,p2) for y in ylist11].plot (xlist11, ylist11, color= 'blue')=
[osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]=
[s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[1],p2),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[10],p1)].plot
(xlist12, ylist12, color= 'blue')= mlab.frange
(ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)= [s_gr(y,p1) for y in
ylist13].plot (xlist13, ylist13, color= 'blue')
pylab.show()()main(): #главная функция она же тело
программы"выберите пункт меню""1- основа 2-энтальпия+энтропия
3-Завершение программы"
sch = int(raw_input(""))sch == 1:
"""p1 = float(raw_input("Введите давление
Р1: "))= float(raw_input("Введите температуру t1: "))=
float(raw_input("Введите температуру t3: "))=
float(raw_input("Введите Расход Воздуха G: "))=
float(raw_input("Введите внутренний относительный КПД турбины: "))=
float(raw_input("Введите внутренний относительный КПД компрессора:
"))= float(raw_input("Введите давление P2: "))= float(raw_input("Введите
показатель адиабыты воздуха: "))"""
os =
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)"T1=",os[0]"T3=",os[1]"tetta=",os[2]"tau=",os[3]"nt=",os[4]"nk=",
os[5]"ngtu=",os[6]"ngtut=",os[7]"kpdt=",os[8]"kpdi=",os[9]"t4d=",os[10]"t2d=",os[11]"Reshaem
zadachi""выберите пункт меню""1-
zadacha 1; 2-zadacha 2; 3-zadacha 2_2; 4-zadacha3 5-zadacha 3_2"=
int(raw_input(""))sch3==1:()sch3==2:()sch3==3:_2()sch3==4:()sch3==5:_2()sch
== 2:= int(raw_input("Введите температуру: "))= t
+ 273.15= float(T)= ental(T)"Entalpiya Kdj/Kg=", h=
entro(T)"Entropiya =", s= teploem(T)"Teploem =", c"1-график энтальпии 2- график энтропии 3 -
teploem"= int(raw_input(""))sch2 == 1:= 273.15= 1000.0= 1.0=
mlab.frange (tmin, tmax, dx)= [ental (x) for x in xlist].plot (xlist,
ylist).show()sch2 == 2:= 273.15= 2000.0= 1.0= mlab.frange (tmin, tmax, dx)=
[entro (x) for x in xlist].plot (xlist, ylist).show()sch2 == 3:= 273.15=
1000.0= 1.0= mlab.frange (tmin, tmax, dx)= [teploem (x) for x in xlist].plot
(xlist, ylist).show()sch == 3:= 0.1= 15= 900= 8.3803 = 0.85= 0.85= 1
k = 1.4()
Код программы расчетов ГТУ с перегревом газа
# -*- coding: cp1251 -*-
#Массивы с константами_c = (29.438265, -1.610822, -11.99174,
68.828384, -98.23993,
64.88351, -20.9094, 2.66524)_h = (-542,
29438.265, -805.411, -3997.2481, 17207.1, -19647.99, 10813.92, -2987.054,
.155)_s = (230.1763, -1.610822, -5.995872,
22.942794, -24.55998, 12.9767, -3.4849, 0.380749)_s = 29.43821
#Считаем энтальпиюmathpylab
matplotlib import mlabental(x):= 0= float(s)i in xrange(len(a_h)):+=
a_h[i]*((x/1000)**i)= (s+488)/28.97r
#Считаем Энтропию
def entro(x):
s = 0= float(s)i in xrange(len(a_s)):+=
a_s[i]*((x/1000)**i)(s + b_s*math.log(x/1000))/28.97
#Cчитаем теплоемкостьteploem(x):
s = 0= float(s)i in xrange(len(a_c)):+=
a_c[i]*((x/1000)**i)= s/28.97r
#Вспомогательная функция энтропии для построения TS диаграмм
def s_gr(x,y):= 0= float(s)=
1*math.log(x/273.15)-0.287*math.log(y/0.1) s
#Основная функция рассчитывающая параметры ГТУ
# osnova(p1,
t1,t3,g,kpdoiT,kpdoiK,p2,k,p5)
def osnova(x1,x2,x3,x4,x5,x6,x7,x8,x9):=
x2+273.15= x3+273.15= X3/X2= (x7/x1)**((x8-1)/x8)= (x7/x9)**((x8-1)/x8)=
(x9/x1)**((x8-1)/x8)= (x4*teploem(X2)*X2*(tau-1))/x6d = X2*(1+(tau-1)/x6)
= G*teploem(X2)*X3*(1-1/tau1)*x5=
G*teploem(X2)*X3*(1-1/tau2)*x5= Ntvd + Ntnd - Nk=
(tetta*x5*x6*(2-1/tau1-1/tau2)-(tau-1))/((tetta-1)*x6-(tau-1)+tetta*x5*x6*(1-1/tau1))=
(tetta*(2-1/tau1-1/tau2)-(tau-1))/((tetta-1)-(tau-1)+tetta*(1-1/tau1)) d =
X3*(1-(1-1/tau2)*x5)d = X3 - (Ntvd/Ntnd)*(X3-T6d)(X2, X3, tetta, tau, tau1,
tau2, Nk, Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)
#Вспомогательная функция для корректной печатиprent(x):
print " ",round(x[0], 3),"
",round(x[1], 3)," ",round(x[2], 3)," ",round(x[3],
3)," ",round(x[4], 3), \
" ",round(x[5], 3)," ",round(x[6],
3)," ", \(x[7], 3)," ",round(x[8], 3),"
",round(x[9], 3)," ",round(x[10], 3)," ",round(x[11],
3)s_gr(x,y):= 0= float(s)=
1*math.log(x/273.15)-0.287*math.log(y/0.1)szadacha1():= 0.5= 0= 0= 0 = 0= 0= 0p2<1.5:=
p1 + 0.000001= 0= 0= 0 = 0= 0= 0 p5< p2:= osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)M1
< os[9]:= os[9]= p2= p5N1 < os[11]:= os[11]= p2= p5+= 0.01M2 < M1:=
M1= MM1= MMM1N2 < N1:= N1= NN1= NNN1+= 0.1M2, MM2, MMM2N2, NN2, NNN2= 0.7=
mlab.frange (p1+0.05, 3, 0.01)= mlab.frange (p1+0.05, 3, 0.01)= [osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[9]
for p2 in xlist1]= [osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[11] for p2 in
xlist2].plot(xlist1, ylist1).show().plot(xlist2, ylist2, color='red').show()
.xlabel('Entropya S').ylabel('Temperatura T').title('TSdiagramma pri
Ngtu=max-red Kpd=max-blue')= osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]=1 =
[ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14]
,NN2)].plot(xlist12, ylist12, color = 'blue')=
mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14],(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],
0.01)= [s_gr(y,NN2) for y in ylist23].plot(xlist23,ylist23,color='blue')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[13]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],NN2),_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[13],NNN2)].plot(xlist34,ylist34,color='blue')
= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[13],
osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1], 0.01)= [s_gr(y,NNN2) for y in
ylist45].plot(xlist45,ylist45,color='blue')=
[osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12]]=
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],NNN2),_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12],p1)].plot(xlist56,ylist56,color='blue')=
mlab.frange (ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12], 0.01)=
[s_gr(y,p1) for y in ylist16].plot(xlist16,ylist16,color='blue')n =
[ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14]]n =
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14]
,MM2)].plot(xlist12n, ylist12n, color = 'red')n =
mlab.frange (osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14],(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],
0.01)n = [s_gr(y,MM2) for y in ylist23n].plot(xlist23n,ylist23n,color='red')n =
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13]]n
= [s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],MM2),_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13],MMM2)].plot(xlist34n,ylist34n,color='red')n
= mlab.frange
(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13],(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],
0.01)n = [s_gr(y,MMM2) for y in ylist45n].plot(xlist45n,ylist45n,color='red')
n =
[osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12]]n =
[s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],MMM2),_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12],p1)].plot(xlist56n,ylist56n,color='red')n
= mlab.frange (ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12], 0.01)n =
[s_gr(y,p1) for y in
ylist16n].plot(xlist16n,ylist16n,color='red').show()zadacha2():= 1000= 0.5= 0=
0= 0 p2<1.5:= p1 + 0.000001p5< p2:= osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)N1
< os[11]:= os[11]= p2= p5+= 0.01+= 0.1N1, NN1, NNN1=
0.5abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>5:+=
0.01abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>1:+= 0.0001gg=
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)"T1=",os[0]"T3=",os[1]"tetta=",os[2]"tau=",os[3]"tau1=",os[4]"tau2=",
os[5]"nk=",os[6]"ntvd=",os[7]"ntnd=",os[8]"ngtu=",os[9]"kpdt=",os[10]"kpdi=",os[11]"t6d=",os[12]"t4d=",os[13]"t2d=",os[14].xlabel('Entropya
S').ylabel('Temperatura T').title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]=1=
[ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]]=
[s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]
,NN1)].plot(xlist12, ylist12, color = 'blue')=
mlab.frange
(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
0.01)= [s_gr(y,NN1) for y in ylist23].plot(xlist23,ylist23,color='blue')=
[osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13]]=
[s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NN1),_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],NNN1)].plot(xlist34,ylist34,color='blue')=
mlab.frange (osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
0.01)= [s_gr(y,NNN1) for y in ylist45].plot(xlist45,ylist45,color='blue')=
[osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12]]=
[s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NNN1),_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12],p1)].plot(xlist56,ylist56,color='blue')
= mlab.frange
(ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12], 0.01)= [s_gr(y,p1) for
y in ylist16].plot(xlist16,ylist16,color='blue').show()
# return (X2, X3, tetta, tau, tau1, tau2, Nk,
Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)zadacha3():= 650= t3= 1000= 0.5= 0=
0= 0 p2<1.5:= p1 + 0.000001p5< p2:=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)N1 < os[11]:= os[11]= p2= p5+=
0.01+= 0.1N1, NN1, NNN1=
0.5abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>5:+=
0.01abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>1:+=
0.0001'rashod pri t3=650', gg= 800= t3= 0.5= 0= 0= 0 p2<1.5:= p1 + 0.000001p5<
p2:= osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)M1 < os[11]:= os[11]= p2= p5+=
0.01+= 0.1M1, MM1, MMM1=
0.5abs(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[9]-Nz)>5:+=
0.01abs(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[9]-Nz)>1:+=
0.0001'rashod pri t3=800', gg2= a.xlabel('Entropya S').ylabel('Temperatura
T').title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]=1=
[ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]]= [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]
,NN1)].plot(xlist12, ylist12, color = 'blue')
= mlab.frange
(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
0.01)= [s_gr(y,NN1) for y in ylist23].plot(xlist23,ylist23,color='blue')=
[osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13]]=
[s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NN1),_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],NNN1)].plot(xlist34,ylist34,color='blue')=
mlab.frange
(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
0.01)= [s_gr(y,NNN1) for y in ylist45].plot(xlist45,ylist45,color='blue')=
[osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12]]=
[s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NNN1),_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12],p1)].plot(xlist56,ylist56,color='blue')=
mlab.frange (ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12], 0.01)=
[s_gr(y,p1) for y in ylist16].plot(xlist16,ylist16,color='blue')
# return (X2, X3, tetta, tau, tau1, tau2, Nk,
Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)= b =
[ymin,osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14]]= [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14]
,MM1)].plot(xlist12, ylist12, color = 'red')=
mlab.frange
(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14],(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],
0.01)= [s_gr(y,MM1) for y in ylist23].plot(xlist23,ylist23,color='red')=
[osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13]]=
[s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],MM1),_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13],MMM1)].plot(xlist34,ylist34,color='red')=
mlab.frange
(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13],(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],
0.01)= [s_gr(y,MMM1) for y in ylist45].plot(xlist45,ylist45,color='red')=
[osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12]]=
[s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],MMM1),_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12],p1)].plot(xlist56,ylist56,color='red')=
mlab.frange (ymin,osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12], 0.01)=
[s_gr(y,p1) for y in ylist16].plot(xlist16,ylist16,color='red').show()= 0.1=
15= 900= t3= 7= 0.85= 0.85= 1= 0.6= 1.4=
osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
"T1=",os[0]"T3=",os[1]"tetta=",os[2]"tau=",os[3]"tau1=",os[4]"tau2=",
os[5]"nk=",os[6]"ntvd=",os[7]"ntnd=",os[8]"ngtu=",os[9]"kpdt=",os[10]"kpdi=",os[11]"t6d=",os[12]"t4d=",os[13]"t2d=",os[14]"1-
1 2-э2 3-задача3 4-выход"=
int(raw_input(""))sch == 1:()sch == 2:()sch == 3:()sch == 4:()