Моделирование газотурбинной установки

Методы расчета термодинамических свойств рабочих тел. Исследование циклов простых газотурбинных установок. Проектирование заданной установки с использованием математической модели. Изучение влияния температур газа перед турбинами на КПД газотурбины.

Рубрика Программирование, компьютеры и кибернетика
Вид курсовая работа
Язык русский
Дата добавления 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

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.