Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского

М А Т И» - Р О С С И Й С К И Й Г О С У Д А Р С Т В Е Н Н Ы Й

Т Е Х Н О Л О Г И Ч Е С К И Й У Н И В Е Р С И Т Е Т

Имени К. Э. Ц И О Л К О В С К О Г О

_______________________________________________________________

Кафедра «Двигатели летательных аппаратов и теплотехника»

Утверждено

Редакционно-издательским

Советом института

РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ ЦИКЛОВ

ТЕПЛОВЫХ МАШИН

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

по курсу «Термодинамика»

Составили: Попов В.Г.

Колесников С.П.

Хахин В.Н.

Москва 2005

1. ЦЕЛЬ РАБОТЫ

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

2. СОДЕРЖАНИЕ ЗАДАНИЯ И ИСХОДНЫЕ ДАННЫЕ

При выполнении курсовой работы по термодинамике студенты рассчитывают замкнутый термодинамический процесс (цикл) по четырём составляющим его политропным участкам. По результатам расчета строятся графики цикла в р-v и Т-s координатах.

Задание на курсовую работу выдается каждому студенту индивидуально на специальном бланке (Таблица 2.1), в котором приводятся необходимые для расчета исходные данные и общий вид рассчитываемого цикла в р-v координатах.

В процессе выполнения курсовой работы по исходным данным требуется определить:

1) Параметры состояния рабочего тела (воздуха): давление р, удельный объём v, температуру Т и энтропию s во всех четырех узловых точках 1,2,3 и 4 цикла.

2) Параметры состояния воздуха (p,v,T и s) в 2-3-х промежуточных точках для каждого из 4-х процессов цикла.

3) Показатель политропы n для каждого процесса цикла, теплоёмкость cn и работу расширения (сжатия) l воздуха в рассчитываемом процессе, количество подведенной (отведенной) в нём теплоты q, а также изменение внутренней энергии Δu, изменение энтальпии Δi воздуха, коэффициент ψ распределения тепла между внутренней энергией и совершаемой воздухом работой в рассматриваемом процессе данного цикла.

4) Работу lц воздуха за цикл, количество теплоты qподв, подведенной к воздуху извне за цикл, и теплоты qотв, отведенной от воздуха в охладитель в ходе совершения цикла, а также количество теплоты qц, превращенной в работу, изменения внутренней энергии Δuц, энтальпии Δiц и энтропии Δsц воздуха за цикл, термический к.п.д. ηt рассчитываемого цикла и термический к.п.д. ηtк, соответствующего ему цикла Карно.

Результаты расчета вносятся в бланк задания и по ним строятся в соответствующем масштабе рабочая р-v (рис.2.1) и тепловая Т-s (рис.2.2) диаграммы рассчитанного термодинамического цикла.

ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮТаблица 2.1

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского

Кафедра «Двигатели летательных аппаратов и теплотехника»

ЗАДАНИЕ ПО РАЗДЕЛУ «ТЕРМОДИНАМИКА»

Студент_________________________Группа_______________________________ФакультетВариант № .

(фамилия, и., о.)

  Узловые точки Промежуточные точки Диа-грамма
Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru x2 ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru x3 ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru x4 ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru x5 ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru x6 ; Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru
р МПа 0,08     0,392             p - v точки х1, х2…..
V м3/кг                    
T К                   T - s точки х11, х12
s кДж/кг·К                    

Произвести расчёт и анализ термодинамического цикла

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru

Для процесса n с кДж/кг·К l кДж/кг q кДж/кг Du кДж/кг Di кДж/кг Ds кДж/кг·К y = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru
1 - 2 к              
2 - 3            
3 - 4                
4 - 1 +∞              

Рабочее тело – воздух, в количестве 1 кг.

Для воздуха:

R = 0,287 Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ;

сv = 0,712 Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru .

Определить величины: а) р;v;Ти sдля

узловых точек цикла; б) n; с; l; q; Du;

Для цикла q подв кДж/кг q отв кДж/кг l ц кДж/кг q ц кДж/кг Δu ц кДж/кг Δi ц кДж/кг Δs ц кДж/кг·К ht
               

Di; Ds и yдля процессов цикла; в) q подв.;

q отв.; lц; qц; Duц; Diц; Dsц и htдля цикла в целом.

Дополнительные данные: v1 / v2 = 12

Выдано «_____»_____________ 200 г. Принято «_____» ________________________200 г. Преподаватель______________________

3. ПОДГОТОВКА К ВЫПОЛНЕНИЮ КУРСОВОЙ РАБОТЫ

Приступая к выполнению курсовой работы, студенты должны проработать разделы, посвященные первому и второму законам термодинамики, а также исследованию политропных процессов и циклов, т.е. должны знать следующие основные соотношения термодинамики, справедливые для идеальных газов [1]:

1. Уравнение состояния идеального газа

рv = RT, (1)

где R – газовая постоянная.

2. Уравнение политропного термодинамического процесса

рvn = const, (2)

где n – показатель политропы, равный

n = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (3)

Здесь и далее индексами «1» и «2» обозначены параметры состояния газа соответственно в начале и в конце рассматриваемого процесса.

3. Уравнение первого закона термодинамики

q = Δu + l. (4)

Здесь Δu – изменение внутренней энергии, которое для всех процессов идеального газа равно

Δu = u2 – u1 = cv(T2 – T1), (5)

где cv – теплоемкость газа в процессе при постоянном объёме, а

l – работа расширения (сжатия) газа, которая во всех процессах, кроме изотермического, может быть определена из выражения

l = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru2 – Т1). (6)

Работа газа в изотермическом процессе равна

l = RT·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = RT·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (7)

Для определения работы газа в адиабатном процессе, кроме выражения (6) можно использовать следующую формулу

l = -Δu = cv(T1 – T2). (8)

Изменение энтальпии для всех процессов идеального газа равно

Δi = i2 – i1 =cp(T2 – T1), (9)

где в соответствии с формулой Майера теплоёмкость газа в процессе при постоянном давлении равна

cp = cv + R. (10)

4. Уравнение для определения количества теплоты

q = cn(T2 – T1), (11)

где cn – теплоемкость газа в политропном процессе, равная

cn = cv + Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , (12)

или

cn = cv Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (13)

Здесь k – показатель адиабаты, равный k = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , для воздуха k=1,4.

5. Аналитическое выражение второго закона термодинамики:

ds = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , (14)

или Δs = s2 – s1 = cn·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , (15)

и Δs = s2 – s1 = cp·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru R·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (16)

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru

Рис.2.1. Рабочая диаграмма цикла.

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru

Рис. 2.2. Тепловая диаграмма цикла

4. ПОРЯДОК ВЫПОЛНЕНИЯ РАБОТЫ

Расчет термодинамического цикла по четырем составляющим его политропным процессам для идеальных газов (воздуха) состоит в применении к ним рассмотренных выше основных соотношений термодинамики.

Выполнение курсовой работы следует начинать с составления плана проведения расчетов. В зависимости от исходных данных производят либо последовательный расчет параметров воздуха от точки 1 до точки 4 цикла, либо сначала определяют параметры воздуха в точках 1, 2, и 4, и лишь затем в точке 3 цикла. В остальном порядок проведения расчетов для различных вариантов задания не имеет существенных различий. Он включает в себя пять основных этапов.

1. Расчет параметров воздуха (p, v, T, s) в узловых точках цикла.

Порядок расчета следующий:

а) определяют начальное состояние воздуха в точке 1, при этом давление, удельный объём или температура вычисляются по формуле (1), а энтропия воздуха – по формуле (16). При этом считают, что для нормальных физических условий (Т0=273 К, P0=105 н/м2), s0 =0. Тогда формула (16) имеет вид:

s1 – s0 = s1 = cP·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru - R·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ; (17)

б) определяют состояние воздуха в точке 2, при этом давление, удельный объём или температура воздуха находятся из уравнений (1) и (2).

Расчетные зависимости имеют следующий вид:

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ; (18)

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ; (19)

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (20)

Энтропия воздуха в точке 2 определяется из уравнения (15) или (16):

s2 = s1 + cn·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru ,

или

s2 = s1 + cР·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru - R·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru .

Аналогично рассчитывают параметры состояния воздуха и в других узловых точках цикла. При этом, если заданными являются величины Δu2-3, Δi2-3 или q2-3, то сначала из уравнений (5), (9) или (11) определяют температуру воздуха в точке 3, а затем по заданным величинам давления или удельного объёма, используя уравнения (1), (15) или (16), вычисляют остальные параметры в точке 3 цикла.

Основные параметры в точке 4 определяют или из уравнений (18-20), описывающих связи между параметрами газа в точках 3 и 4 (когда в точке 4 задан лишь один параметр из трёх и задан показатель политропы n3-4), или непосредственно из уравнения состояния (когда в точке 4 заданы два параметра из трёх). При этом, если в последнем случае задаются также показатель политропы n3-4 и лишь один параметр воздуха в точке 3, то состояние воздуха в этой точке цикла определяют после точки 4, используя уравнения (18-20).

2. Расчет параметров воздуха (p, v, Т, s) в промежуточных точках процессов цикла.

Промежуточные точки назначаются студентами самостоятельно и раздельно для диаграмм p-v и T-s цикла таким образом, чтобы на каждом криволинейном участке диаграммы было не менее 2-х точек.

Принятые промежуточные точки обозначаются через Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru и т.д. для рабочей (p-v) и Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru и т.д. для тепловой (Т-s) диаграмм цикла. Параметры состояния воздуха в промежуточных точках цикла вычисляются по формулам (18) и (15), которые при принятых обозначениях принимают следующий вид:

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = P1 Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , (21)

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = s1 + cnln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , (22)

где значениями удельного объёма vx в точке x и температуры Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru в точке х1 цикла следует задаваться при расчете таким образом, чтобы промежуточные точки равномерно располагались на исследуемом участке диаграммы.

3. Построение рабочей (p-v) и тепловой (T-s) диаграмм цикла. Графики строят по узловым и промежуточным точкам, имея при этом в виду, что:

- цикл как в рабочей, так и в тепловой диаграммах, изображается замкнутой кривой, причем направление процессов цикла в этих диаграммах одинаково – по ходу движения часовой стрелки;

- все политропные процессы в тепловой диаграмме представляются логарифмическими кривыми, направленными выпуклостью вниз;

- расположение политропных процессов в рабочей диаграмме зависит от величины показателя политропы: при ∞>n>0, как это следует из уравнения (18), политропа имеет гиперболический характер; при 0>n>- ∞ политропы представляют собой кривые, проходящие через начало координат, причем при n > -1 кривая процесса направлена выпуклостью вверх, а при n < -1 направлена выпуклостью вниз, в частности при n = -1 политропа является прямой линией.

4. Определение величин, характеризующих процессы цикла:

а) показатель политропы вычисляется по формуле (3);

б) теплоёмкость воздуха в рассматриваемом процессе определяется по формуле (12) или (13), причем теплоёмкость отрицательна, если 1<n<k;

в) работа, производимая воздухом в политропном процессе, определяется по одному из выражений (6), (7) или (8); полученная величина должна быть численно равна площади под кривой этого процесса в (р-v) диаграмме; работа считается положительной при расширении газа и отрицательной при сжатии газа;

г) количество подводимой (или отводимой) теплоты в политропном процессе вычисляется по формуле (4) или (11); найденная величина должна быть численно равна площади под кривой процесса в (T-s) диаграмме; согласно уравнению (14), теплота подводится к газу извне (т.е. она положительная), если энтропия увеличивается (ds>0) и, наоборот, теплота отводится от газа в окружающую среду (т.е. она отрицательная), если энтропия уменьшается (ds<0).

д) изменение внутренней энергии воздуха для процесса определяется по формуле (5);

е) изменение энтальпии воздуха определяется по формуле (9);

ж) изменение энтропии определяется по формулам (15) или (16);

з) коэффициент распределения тепла между внутренней энергией и совершаемой работой воздуха за процесс находится по формуле:

ψ = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (23)

5. Определение величин, характеризующих цикл в целом:

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

lц = l1-2 + l2-3 + l3-4 + l4-1 = ∑ l ; (24)

в соответствии с этим выражением за один цикл воздух совершает полезную работу, равную площади цикла, представленного в (p-v) диаграмме;

б) количество теплоты, превращенной в работу цикла, равно:

qц = qподв – qотв, (25)

где qподв – теплота, подведенная к газу извне за цикл, а qотв – теплота, отведенная от газа в охладитель в ходе совершения цикла.

Из уравнения (25) видно, что теплота цикла qц, превращенная в работу цикла lц, определяется площадью цикла в (T-s) диаграмме. Следовательно, площади цикла в (p-v) и (T-s) диаграммах имеют одинаковые значения;

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

Δuц = Δu1-2 + Δu2-3 + Δu3-4 + Δu4-1 = 0, (26)

Δiц = Δi1-2 + Δi2-3 + Δi3-4 + Δi4-1 = 0, (27)

Δsц = Δs1-2 + Δs2-3 + Δs3-4 + Δs4-1 = 0 (28)

г) термический к.п.д. цикла, который характеризует степень совершенства преобразования теплоты в работу, равен отношению полезной работы lц к подведенной теплоте qподв:

ηt = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (29)

Из этого выражения следует, что термический к.п.д. цикла, представленного в (T-s) диаграмме, определяется отношением площади цикла к площади, соответствующей количеству подводимой теплоты.

Совершенство произвольного обратимого цикла оценивается сравнением его термического к.п.д. с термическим к.п.д. обратимого цикла Карно, осуществляемого между крайними (Tmax и Tmin) температурами рассматриваемого цикла. Цикл Карно имеет максимальный термический к.п.д. для заданного интервала температур, равный:

ηtk = 1 - Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru . (30)

Если сравнить рассчитываемый цикл с описанным вокруг него циклом Карно в (Т-s) диаграмме, то площадь рассматриваемого цикла получается меньше площади цикла Карно. Отношение этих площадей показывает насколько близко приближается рассматриваемый цикл к циклу Карно, т.е. определяют насколько термический к.п.д. рассматриваемого цикла меньше термического к.п.д. цикла Карно:

Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru .

5. ПРИМЕР ВЫПОЛНЕНИЯ КУРСОВОЙ РАБОТЫ

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

1. Расчет параметров p,v,T,s в узловых точках цикла:

а) начальное состояние воздуха в точке 1 согласно исходным данным и формулам (1) и (17) составляет:

Р1 = 0,08 МПа = 0,8·105 н/м2; Т1 = 280 К;

v1 = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = 1 м3/кг;

s1 = cp·ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = 1,0∙ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru - 287∙ln Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru = 0,092 кДж/кг∙К,

где согласно формуле (10) сp = сv + R = 0,712 + 0,287 ≈ 1 кДж/кг∙К;

б) состояние воздуха в точке 2 можно определить, используя соотношение Российский Государственный Технологический Университет «МАТИ» им. К. Э. Циолковского - student2.ru , когда сначала по вычисленному значению v1 находится удельный объём v2, а затем по формуле (18) - давление Р2; температура Т2 находится из уравнения состояния (1); так как процесс 1-2 адиабатный (n = k), то изменение энтропии в этом процессе Δs1-2 = 0 (сn= 0) и s2 = s1;

в) учитывая далее, что процесс 2-3 является изобарным (n = 0), находим давление Р3 = Р2, а так как при этом теплоёмкость cn = cp, то по формуле (11) по заданному значению q2-3 и известной температуре Т2 вычисляем температуру воздуха Т3 в точке 3 цикла; энтропию s3 находим по формуле (15);

г) учитывая, что процесс 4-1 изохорный (n = + ∞), находим удельный объём воздуха v4 = v1, а далее по заданному значению Р4 и вычисленному значению v4 определяем температуру Т4 в точке 4 цикла, энтропию воздуха s4 вычисляем по формуле (15), учитывая при этом, что теплоёмкость газа в процессе 4-1 cn = cv.

2. Расчет параметров p,v,T,s в промежуточных точках процессов цикла (см. п.2, раздел 4).

3. Построение рабочей (p-v) и тепловой (T-s) диаграмм цикла. По вычисленным значениям параметров состояния воздуха в узловых и промежуточных точках цикла в масштабе строятся рабочая (рис.2.1.) и тепловая (рис.2.2.)диаграммы цикла.

4. Определение величин, характеризующих процессы цикла (см. п.4, раздел 4).

5. Определение величин, характеризующих цикл в целом (см. п.5, раздел 4).

Результаты расчетов сводятся в общую таблицу 5.1.

6. ОТЧЕТНОСТЬ ПО РАБОТЕ

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

Пояснительная записка должна содержать:

- оформленный и заполненный результатами расчетов бланк задания на курсовую работу (см. таблицу 5.1.);

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

- рабочую и тепловую диаграммы цикла (см. рис.2.1 и рис.2.2).

Все расчеты выполняются в системе СИ: давление измеряется в н/м2 или в Па; температура Т - в К; удельный объём v - в м3/кг; теплота q, внутренняя энергия u и энтальпия i - в кДж/кг; теплоёмкость cn, энтропия s и газовая постоянная R для 1 кг воздуха - в кДж/кг∙К.

ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮТаблица 5.1

Наши рекомендации