Моделирование газотурбинной установки
Методы расчета термодинамических свойств рабочих тел. Исследование циклов простых газотурбинных установок. Проектирование заданной установки с использованием математической модели. Изучение влияния температур газа перед турбинами на КПД газотурбины.
| Рубрика | Программирование, компьютеры и кибернетика |
| Вид | курсовая работа |
| Язык | русский |
| Дата добавления | 30.11.2013 |
| Размер файла | 436,9 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
ymin = osnova(p1,t1,t3,G,kk,kk,p2,k)[0]
ymax = osnova(p1,t1,t3,G,kk,kk,p2,k)[1]
dy = 1
ylist10 = [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]
xlist10 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)]
pylab.plot (xlist10, ylist10, color= 'red')
ylist11 = mlab.frange (osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)
xlist11 = [s_gr(y,p2) for y in ylist11]
pylab.plot (xlist11, ylist11, color= 'red')
ylist12 = [osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]
xlist12 = [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)]
pylab.plot (xlist12, ylist12, color= 'red')
ylist13 = mlab.frange (ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)
xlist13 = [s_gr(y,p1) for y in ylist13]
pylab.plot (xlist13, ylist13, color= 'red')
kk = 0.9
S = 1
while S>0:
kk = kk-0.001
S = osnova(p1,t1,t3,G,kk,kk,p2,k)[6]
t3 = 800
ymin = osnova(p1,t1,t3,G,kk,kk,p2,k)[0]
ymax = osnova(p1,t1,t3,G,kk,kk,p2,k)[1]
dy = 1
ylist10 = [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]
xlist10 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)]
pylab.plot (xlist10, ylist10, color= 'blue')
ylist11 = mlab.frange (osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)
xlist11 = [s_gr(y,p2) for y in ylist11]
pylab.plot (xlist11, ylist11, color= 'blue')
ylist12 = [osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]
xlist12 = [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)]
pylab.plot (xlist12, ylist12, color= 'blue')
ylist13 = mlab.frange (ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)
xlist13 = [s_gr(y,p1) for y in ylist13]
pylab.plot (xlist13, ylist13, color= 'blue')
pylab.show()
main()
def zadacha3_2():
t3 = 400
p2 = 1
ass = []
bss = []
while t3 <= 1000:
print 't3=', t3
a = []
b = []
kk = 0.95
S = 1
while S>0:
a = a+ [kk]
b = b + [osnova(p1,t1,t3,G,kk,kk,p2,k)[9]]
S = osnova(p1,t1,t3,G,kk,kk,p2,k)[9]
print ' kk=', kk, ' kpd gtu = ', S
kk = kk-0.001
ass = ass + [a]
bss = bss + [b]
t3 += 200
i = 0
while i <int(len(bss)):
pylab.plot (ass[i], bss[i], color= 'blue')
i += 1
pylab.show()
t3 = 800
kk = 0.9
p2 = 1.1
pylab.xlabel('Entropya S')
pylab.ylabel('Temperatura T')
pylab.title('TSdiagramma pri kpd=0.9-red ngtu=0-blue')
ymin = osnova(p1,t1,t3,G,kk,kk,p2,k)[0]
ymax = osnova(p1,t1,t3,G,kk,kk,p2,k)[1]
dy = 1
ylist10 = [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]
xlist10 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)]
pylab.plot (xlist10, ylist10, color= 'red')
ylist11 = mlab.frange (osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)
xlist11 = [s_gr(y,p2) for y in ylist11]
pylab.plot (xlist11, ylist11, color= 'red')
ylist12 = [osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]
xlist12 = [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)]
pylab.plot (xlist12, ylist12, color= 'red')
ylist13 = mlab.frange (ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)
xlist13 = [s_gr(y,p1) for y in ylist13]
pylab.plot (xlist13, ylist13, color= 'red')
kk = 0.9
S = 1
while S>0:
kk = kk-0.001
S = osnova(p1,t1,t3,G,kk,kk,p2,k)[9]
ymin = osnova(p1,t1,t3,G,kk,kk,p2,k)[0]
ymax = osnova(p1,t1,t3,G,kk,kk,p2,k)[1]
dy = 1
ylist10 = [ymin, osnova(p1,t1,t3,G,kk,kk,p2,k)[11]]
xlist10 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,kk,kk,p2,k)[11],p2)]
pylab.plot (xlist10, ylist10, color= 'blue')
ylist11 = mlab.frange (osnova(p1,t1,t3,G,kk,kk,p2,k)[11],osnova(p1,t1,t3,G,kk,kk,p2,k)[1], 0.01)
xlist11 = [s_gr(y,p2) for y in ylist11]
pylab.plot (xlist11, ylist11, color= 'blue')
ylist12 = [osnova(p1,t1,t3,G,kk,kk,p2,k)[1],osnova(p1,t1,t3,G,kk,kk,p2,k)[10]]
xlist12 = [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)]
pylab.plot (xlist12, ylist12, color= 'blue')
ylist13 = mlab.frange (ymin,osnova(p1,t1,t3,G,kk,kk,p2,k)[10], 0.01)
xlist13 = [s_gr(y,p1) for y in ylist13]
pylab.plot (xlist13, ylist13, color= 'blue')
pylab.show()
main()
def main(): #главная функция она же тело программы
print "выберите пункт меню"
print "1- основа 2-энтальпия+энтропия 3-Завершение программы"
sch = int(raw_input(""))
if sch == 1:
"""p1 = float(raw_input("Введите давление Р1: "))
t1 = float(raw_input("Введите температуру t1: "))
t3 = float(raw_input("Введите температуру t3: "))
G = float(raw_input("Введите Расход Воздуха G: "))
KPDoit = float(raw_input("Введите внутренний относительный КПД турбины: "))
KPDoik = float(raw_input("Введите внутренний относительный КПД компрессора: "))
p2 = float(raw_input("Введите давление P2: "))
k = float(raw_input("Введите показатель адиабыты воздуха: "))"""
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k)
print "T1=",os[0]
print "T3=",os[1]
print "tetta=",os[2]
print "tau=",os[3]
print "nt=",os[4]
print "nk=", os[5]
print "ngtu=",os[6]
print "ngtut=",os[7]
print "kpdt=",os[8]
print "kpdi=",os[9]
print "t4d=",os[10]
print "t2d=",os[11]
print "Reshaem zadachi"
print "выберите пункт меню"
print "1- zadacha 1; 2-zadacha 2; 3-zadacha 2_2; 4-zadacha3 5-zadacha 3_2"
sch3 = int(raw_input(""))
if sch3==1:
zadacha1()
elif sch3==2:
zadacha2()
elif sch3==3:
zadacha2_2()
elif sch3==4:
zadacha3()
elif sch3==5:
zadacha3_2()
elif sch == 2:
t = int(raw_input("Введите температуру: "))
T = t + 273.15
T = float(T)
h = ental(T)
print "Entalpiya Kdj/Kg=", h
s = entro(T)
print "Entropiya =", s
c = teploem(T)
print "Teploem =", c
print "1-график энтальпии 2- график энтропии 3 - teploem"
sch2 = int(raw_input(""))
if sch2 == 1:
tmin = 273.15
tmax = 1000.0
dx = 1.0
xlist = mlab.frange (tmin, tmax, dx)
ylist = [ental (x) for x in xlist]
pylab.plot (xlist, ylist)
pylab.show()
elif sch2 == 2:
tmin = 273.15
tmax = 2000.0
dx = 1.0
xlist = mlab.frange (tmin, tmax, dx)
ylist = [entro (x) for x in xlist]
pylab.plot (xlist, ylist)
pylab.show()
elif sch2 == 3:
tmin = 273.15
tmax = 1000.0
dx = 1.0
xlist = mlab.frange (tmin, tmax, dx)
ylist = [teploem (x) for x in xlist]
pylab.plot (xlist, ylist)
pylab.show()
elif sch == 3:
exit
p1 = 0.1
t1 = 15
t3 = 900
G= 8.3803
KPDoit = 0.85
KPDoik = 0.85
p2 = 1
k = 1.4
main()
Код программы расчетов ГТУ с перегревом газа
# -*- coding: cp1251 -*-
#Массивы с константами
a_c = (29.438265, -1.610822, -11.99174, 68.828384, -98.23993,
64.88351, -20.9094, 2.66524)
a_h = (-542, 29438.265, -805.411, -3997.2481, 17207.1, -19647.99, 10813.92, -2987.054,
333.155)
a_s = (230.1763, -1.610822, -5.995872, 22.942794, -24.55998, 12.9767, -3.4849, 0.380749)
b_s = 29.43821
#Считаем энтальпию
import math
import pylab
from matplotlib import mlab
def ental(x):
s = 0
s = float(s)
for i in xrange(len(a_h)):
s += a_h[i]*((x/1000)**i)
r = (s+488)/28.97
return r
#Считаем Энтропию
def entro(x):
s = 0
s = float(s)
for i in xrange(len(a_s)):
s += a_s[i]*((x/1000)**i)
return (s + b_s*math.log(x/1000))/28.97
#Cчитаем теплоемкость
def teploem(x):
s = 0
s = float(s)
for i in xrange(len(a_c)):
s += a_c[i]*((x/1000)**i)
r = s/28.97
return r
#Вспомогательная функция энтропии для построения TS диаграмм
def s_gr(x,y):
s = 0
s = float(s)
s = 1*math.log(x/273.15)-0.287*math.log(y/0.1)
return s
#Основная функция рассчитывающая параметры ГТУ
# osnova(p1, t1,t3,g,kpdoiT,kpdoiK,p2,k,p5)
def osnova(x1,x2,x3,x4,x5,x6,x7,x8,x9):
X2 = x2+273.15
X3 = x3+273.15
tetta = X3/X2
tau = (x7/x1)**((x8-1)/x8)
tau1 = (x7/x9)**((x8-1)/x8)
tau2 = (x9/x1)**((x8-1)/x8)
Nk = (x4*teploem(X2)*X2*(tau-1))/x6
T2d = X2*(1+(tau-1)/x6)
Ntvd = G*teploem(X2)*X3*(1-1/tau1)*x5
Ntnd = G*teploem(X2)*X3*(1-1/tau2)*x5
Ngtu = Ntvd + Ntnd - Nk
KPDi = (tetta*x5*x6*(2-1/tau1-1/tau2)-(tau-1))/((tetta-1)*x6-(tau-1)+tetta*x5*x6*(1-1/tau1))
KPDt = (tetta*(2-1/tau1-1/tau2)-(tau-1))/((tetta-1)-(tau-1)+tetta*(1-1/tau1))
T6d = X3*(1-(1-1/tau2)*x5)
T4d = X3 - (Ntvd/Ntnd)*(X3-T6d)
return (X2, X3, tetta, tau, tau1, tau2, Nk, Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)
#Вспомогательная функция для корректной печати
def 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)," ", \
round(x[7], 3)," ",round(x[8], 3)," ",round(x[9], 3)," ",round(x[10], 3)," ",round(x[11], 3)
def s_gr(x,y):
s = 0
s = float(s)
s = 1*math.log(x/273.15)-0.287*math.log(y/0.1)
return s
def zadacha1():
p2 = 0.5
M2 = 0
MM2 = 0
MMM2 = 0
N2 = 0
NN2 = 0
NNN2 = 0
while p2<1.5:
p5 = p1 + 0.000001
M1 = 0
MM1 = 0
MMM1 = 0
N1 = 0
NN1 = 0
NNN1 = 0
while p5< p2:
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
if M1 < os[9]:
M1 = os[9]
MM1 = p2
MMM1 = p5
if N1 < os[11]:
N1 = os[11]
NN1 = p2
NNN1 = p5
p5 += 0.01
if M2 < M1:
M2 = M1
MM2 = MM1
MMM2 = MMM1
if N2 < N1:
N2 = N1
NN2 = NN1
NNN2 = NNN1
p2 += 0.1
print M2, MM2, MMM2
print N2, NN2, NNN2
p5 = 0.7
xlist1 = mlab.frange (p1+0.05, 3, 0.01)
xlist2 = mlab.frange (p1+0.05, 3, 0.01)
ylist1 = [osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[9] for p2 in xlist1]
ylist2 = [osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[11] for p2 in xlist2]
pylab.plot(xlist1, ylist1)
pylab.show()
pylab.plot(xlist2, ylist2, color='red')
pylab.show()
pylab.xlabel('Entropya S')
pylab.ylabel('Temperatura T')
pylab.title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]
ymax = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]
dy =1
ylist12 = [ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14]]
xlist12 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14]
,NN2)]
pylab.plot(xlist12, ylist12, color = 'blue')
ylist23 = mlab.frange (osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[14],
osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1], 0.01)
xlist23 = [s_gr(y,NN2) for y in ylist23]
pylab.plot(xlist23,ylist23,color='blue')
ylist34 = [osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],
osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[13]]
xlist34 = [s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],NN2),
s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[13],NNN2)]
pylab.plot(xlist34,ylist34,color='blue')
ylist45 = 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)
xlist45 = [s_gr(y,NNN2) for y in ylist45]
pylab.plot(xlist45,ylist45,color='blue')
ylist56 = [osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],
osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12]]
xlist56 = [s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[1],NNN2),
s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12],p1)]
pylab.plot(xlist56,ylist56,color='blue')
ylist16 = mlab.frange (ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,NN2,k,NNN2)[12], 0.01)
xlist16 = [s_gr(y,p1) for y in ylist16]
pylab.plot(xlist16,ylist16,color='blue')
ylist12n = [ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14]]
xlist12n = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14]
,MM2)]
pylab.plot(xlist12n, ylist12n, color = 'red')
ylist23n = mlab.frange (osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[14],
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1], 0.01)
xlist23n = [s_gr(y,MM2) for y in ylist23n]
pylab.plot(xlist23n,ylist23n,color='red')
ylist34n = [osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13]]
xlist34n = [s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],MM2),
s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13],MMM2)]
pylab.plot(xlist34n,ylist34n,color='red')
ylist45n = mlab.frange (osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[13],
osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1], 0.01)
xlist45n = [s_gr(y,MMM2) for y in ylist45n]
pylab.plot(xlist45n,ylist45n,color='red')
ylist56n = [osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1], osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12]]
xlist56n = [s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[1],MMM2),
s_gr(osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12],p1)]
pylab.plot(xlist56n,ylist56n,color='red')
ylist16n = mlab.frange (ymin,osnova(p1,t1,t3,G,KPDoit,KPDoik,MM2,k,MMM2)[12], 0.01)
xlist16n = [s_gr(y,p1) for y in ylist16n]
pylab.plot(xlist16n,ylist16n,color='red')
pylab.show()
def zadacha2():
Nz = 1000
p2 = 0.5
N1 = 0
NN1 = 0
NNN1 = 0
while p2<1.5:
p5 = p1 + 0.000001
while p5< p2:
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
if N1 < os[11]:
N1 = os[11]
NN1 = p2
NNN1 = p5
p5 += 0.01
p2 += 0.1
print N1, NN1, NNN1
gg = 0.5
while abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>5:
gg += 0.01
while abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>1:
gg += 0.0001
print gg
os = osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)
print "T1=",os[0]
print "T3=",os[1]
print "tetta=",os[2]
print "tau=",os[3]
print "tau1=",os[4]
print "tau2=", os[5]
print "nk=",os[6]
print "ntvd=",os[7]
print "ntnd=",os[8]
print "ngtu=",os[9]
print "kpdt=",os[10]
print "kpdi=",os[11]
print "t6d=",os[12]
print "t4d=",os[13]
print "t2d=",os[14]
pylab.xlabel('Entropya S')
pylab.ylabel('Temperatura T')
pylab.title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]
ymax = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]
dy =1
ylist12 = [ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]]
xlist12 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]
,NN1)]
pylab.plot(xlist12, ylist12, color = 'blue')
ylist23 = mlab.frange (osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1], 0.01)
xlist23 = [s_gr(y,NN1) for y in ylist23]
pylab.plot(xlist23,ylist23,color='blue')
ylist34 = [osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13]]
xlist34 = [s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NN1),
s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],NNN1)]
pylab.plot(xlist34,ylist34,color='blue')
ylist45 = mlab.frange (osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1], 0.01)
xlist45 = [s_gr(y,NNN1) for y in ylist45]
pylab.plot(xlist45,ylist45,color='blue')
ylist56 = [osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12]]
xlist56 = [s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NNN1),
s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12],p1)]
pylab.plot(xlist56,ylist56,color='blue')
ylist16 = mlab.frange (ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12], 0.01)
xlist16 = [s_gr(y,p1) for y in ylist16]
pylab.plot(xlist16,ylist16,color='blue')
pylab.show()
# return (X2, X3, tetta, tau, tau1, tau2, Nk, Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)
def zadacha3():
t3 = 650
a = t3
Nz = 1000
p2 = 0.5
N1 = 0
NN1 = 0
NNN1 = 0
while p2<1.5:
p5 = p1 + 0.000001
while p5< p2:
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
if N1 < os[11]:
N1 = os[11]
NN1 = p2
NNN1 = p5
p5 += 0.01
p2 += 0.1
print N1, NN1, NNN1
gg = 0.5
while abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>5:
gg += 0.01
while abs(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[9]-Nz)>1:
gg += 0.0001
print 'rashod pri t3=650', gg
t3 = 800
b = t3
p2 = 0.5
M1 = 0
MM1 = 0
MMM1 = 0
while p2<1.5:
p5 = p1 + 0.000001
while p5< p2:
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
if M1 < os[11]:
M1 = os[11]
MM1 = p2
MMM1 = p5
p5 += 0.01
p2 += 0.1
print M1, MM1, MMM1
gg2 = 0.5
while abs(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[9]-Nz)>5:
gg2 += 0.01
while abs(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[9]-Nz)>1:
gg2 += 0.0001
print 'rashod pri t3=800', gg2
t3 = a
pylab.xlabel('Entropya S')
pylab.ylabel('Temperatura T')
pylab.title('TSdiagramma pri Ngtu=max-red Kpd=max-blue')
ymin = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[0]
ymax = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)[1]
dy =1
ylist12 = [ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]]
xlist12 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14]
,NN1)]
pylab.plot(xlist12, ylist12, color = 'blue')
ylist23 = mlab.frange (osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[14],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1], 0.01)
xlist23 = [s_gr(y,NN1) for y in ylist23]
pylab.plot(xlist23,ylist23,color='blue')
ylist34 = [osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13]]
xlist34 = [s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NN1),
s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],NNN1)]
pylab.plot(xlist34,ylist34,color='blue')
ylist45 = mlab.frange (osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[13],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1], 0.01)
xlist45 = [s_gr(y,NNN1) for y in ylist45]
pylab.plot(xlist45,ylist45,color='blue')
ylist56 = [osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],
osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12]]
xlist56 = [s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[1],NNN1),
s_gr(osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12],p1)]
pylab.plot(xlist56,ylist56,color='blue')
ylist16 = mlab.frange (ymin,osnova(p1,t1,t3,gg,KPDoit,KPDoik,NN1,k,NNN1)[12], 0.01)
xlist16 = [s_gr(y,p1) for y in ylist16]
pylab.plot(xlist16,ylist16,color='blue')
# return (X2, X3, tetta, tau, tau1, tau2, Nk, Ntvd, Ntnd, Ngtu, KPDt, KPDi, T6d, T4d, T2d)
t3 = b
ylist12 = [ymin,osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14]]
xlist12 = [s_gr(ymin,p1),s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14]
,MM1)]
pylab.plot(xlist12, ylist12, color = 'red')
ylist23 = mlab.frange (osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[14],
osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1], 0.01)
xlist23 = [s_gr(y,MM1) for y in ylist23]
pylab.plot(xlist23,ylist23,color='red')
ylist34 = [osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],
osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13]]
xlist34 = [s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],MM1),
s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13],MMM1)]
pylab.plot(xlist34,ylist34,color='red')
ylist45 = mlab.frange (osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[13],
osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1], 0.01)
xlist45 = [s_gr(y,MMM1) for y in ylist45]
pylab.plot(xlist45,ylist45,color='red')
ylist56 = [osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],
osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12]]
xlist56 = [s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[1],MMM1),
s_gr(osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12],p1)]
pylab.plot(xlist56,ylist56,color='red')
ylist16 = mlab.frange (ymin,osnova(p1,t1,t3,gg2,KPDoit,KPDoik,MM1,k,MMM1)[12], 0.01)
xlist16 = [s_gr(y,p1) for y in ylist16]
pylab.plot(xlist16,ylist16,color='red')
pylab.show()
p1 = 0.1
t1 = 15
t3 = 900
t5 = t3
G= 7
KPDoit = 0.85
KPDoik = 0.85
p2 = 1
p5 = 0.6
k = 1.4
os = osnova(p1,t1,t3,G,KPDoit,KPDoik,p2,k,p5)
print "T1=",os[0]
print "T3=",os[1]
print "tetta=",os[2]
print "tau=",os[3]
print "tau1=",os[4]
print "tau2=", os[5]
print "nk=",os[6]
print "ntvd=",os[7]
print "ntnd=",os[8]
print "ngtu=",os[9]
print "kpdt=",os[10]
print "kpdi=",os[11]
print "t6d=",os[12]
print "t4d=",os[13]
print "t2d=",os[14]
print "1- 1 2-э2 3-задача3 4-выход"
sch = int(raw_input(""))
if sch == 1:
zadacha1()
if sch == 2:
zadacha2()
if sch == 3:
zadacha3()
if sch == 4:
exit()
Размещено на Allbest.ru
Подобные документы
Разработка системы расчета характеристик разомкнутых экспоненциальных сетевых моделей, выполняющая имитационное моделирование заданной сетевой модели. Построение модели на языке GPSS, анализ эффективности аналитической модели, выполняющей роль эталона.
курсовая работа [483,6 K], добавлен 01.12.2010Разработка математической модели и неявной конечно-разностной схемы для получения динамики изменения температур заготовки в период нагрева. Распределение температур по сечению сляба. Разработка алгоритма и блок-схемы, отладка прикладной программы для ЭВМ.
курсовая работа [658,5 K], добавлен 30.06.2011Проектирование установки, предназначенной для быстрого прототипирования (печати пластиковых моделей по готовой 3D-модели). Укрупнённая структурная схема системного проектирования. Разработка корпуса автоматизированной установки. Внешний вид контроллера.
дипломная работа [3,2 M], добавлен 10.01.2015Расчет параметров моделирования в системе Fortran. Описание алгоритма и математической модели системы, их составляющих. Моделирование шума с заданной плотностью распределения вероятностей. Выполнение моделирования работы системы при входном сигнале N(t).
курсовая работа [896,3 K], добавлен 20.06.2012Разработка математической модели заданной системы, ее внутреннее содержание и взаимосвязь отдельных компонентов. Формирование алгоритма решения задачи. Проектирование программы и ее листинг, порядок и основные этапы проведения необходимых расчетов.
контрольная работа [656,9 K], добавлен 02.02.2015Изучение элементов языка С++, программирование разветвлений и циклов с использованием операторов условного и перехода. Обработка одномерных массивов. Поиск максимального элемента массива с заданной размерностью. Листинги программы и результатов.
курсовая работа [647,7 K], добавлен 05.02.2013Исследование метода математического моделирования чрезвычайной ситуации. Модели макрокинетики трансформации веществ и потоков энергии. Имитационное моделирование. Процесс построения математической модели. Структура моделирования происшествий в техносфере.
реферат [240,5 K], добавлен 05.03.2017Табличный вывод значений суммы ряда и номера последнего элемента суммы в зависимости от значений величин входных параметров с применением операторов ветвления и циклов. Блок-схема алгоритма решения. Время работы программы для расчета одного значения.
контрольная работа [762,9 K], добавлен 14.05.2013Свойства и методы формирования криптопараметров и оценка стойкости. Криптографические хэш-функции. Методы и алгоритмы формирования рабочих ключей. Моделирование упрощенной модели электронной цифровой подписи файла с использованием метода Шнорра.
курсовая работа [47,9 K], добавлен 14.12.2012Создание web-страниц с использованием HTML. Работа с графикой в Adobe Photoshop и Flash. Создание динамических web-страниц с использованием JavaScript. Пример реализации "Эконометрической модели экономики России". Моделирование с использованием Powersim.
презентация [478,4 K], добавлен 25.09.2013
