Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Криогенные поршневые детандеры

..pdf
Скачиваний:
45
Добавлен:
24.10.2023
Размер:
10.95 Mб
Скачать

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

Теплоприток к рабочему газу определяется при прочих рав­ ных условиях коэффициентом теплопередачи от окружающей среды к газу через стенку цилиндра и зависит в общем случае от а 0.с и а газа, но так как а0.с «С а газа, то определяющим, оче­ видно, будет Оо.с- Поэтому можно предположить, что удельный тепловой поток (Вт/м2) будет постоянной величиной.

Тогда

 

Q o . c

(169)

где

— удельный теплоприток к газу

(Вт/кг).

Влияние увеличения оборотности на потери от регенератив­ ного теплообмена можно проследить по формуле (97).

До настоящего времени нет достоверных данных о взаимо­ связи а со скоростью движения поршня. Экспериментально полученные значения а существенно превышают значения, под­ считанные по средней скорости движения поршня. Поэтому в данном расчете можно допустить независимость а от скорости

движения поршня.

_____

 

Тогда

 

 

Ят*= ?,• - { /- (^гЛ Вт/кг-

(17°)

УV п ’ )

(повторяется зависимость для q0.с).

уплотнении,

Тепло трения, выделяющееся в поршневом

можно представить в виде

 

QTP = frPnDb • 2Sp,-60 • ni,

(171)

где fTp— коэффициент трения; S, D — ход и диаметр поршня;

b — рабочая ширина уплотнительного элемента; п — число обо­ ротов вала детандера; i — количество уплотнительных элемен­ тов; pi — среднее давление за цикл по уплотнительным элементам.

Давление pi зависит в общем случае как от газовых усилий, так и от упругости уплотнительного элемента. Полагая, что по­ следняя не меняется, найдем, что pi = const и не зависит от числа оборотов вала машины. Полагая из конструктивных со­

ображений, что Ь'

= 6

 

f тр

 

Q r p — Q r p

 

f тр

, получим при i = const

(172)

103

В случае независимости коэффициента трения от скорости поршня = Qrр и с ростом числа оборотов машины потери

от трения не изменяются. В случае же, когда скорость4 поршня возрастает до величин, когда она оказывает влияние на fTP, то <2тр становится функцией числа оборотов. Однако в диапазоне реальных скоростей поршня расширительных поршневых ма­ шин такая зависимость маловероятна и с достаточной степенью точности можно полагать /тр, а следовательно, и QTP независи­ мыми от скорости вращения вала машины.

Таким образом, при увеличении числа оборотов до величин, при которых удается сохранить неизменными а0 и сопротивле­ ния в клапанах, увеличение числа оборотов приводит к умень­ шению потерь от теплопритоков (из окружающей среды и от регенеративного теплообмена) и к увеличению термодинамиче­ ской эффективности детандеров. Этот случай наиболее харак­ терен для прямоточных машин с одним клапаном впуска и машин с внутренним приводом клапанов (где имеется достаточ­ ный запас места для расположения органов газораспределе­ ния), а также для бесклапанных машин, в которых мертвое пространство должно быть существенно увеличено.

В проведенном анализе допускалось постоянство и2 и v5, причем если постоянство и2 вполне вероятно, то vs должно изменяться; и так как с ростом числа оборотов цад должен рас­ ти, то естественно допустить, что Т5 будет падать, а и5 — расти, что в конечном итоге приведет к возможности еще больше уменьшить У0, чем по соотношению

а это еще больше позволит почувствовать положительный

эффект от увеличения п.

Другой случай наблюдается, когда дальнейшее увеличение числа оборотов машины не позволяет разместить в цилиндре органы газораспределения без уменьшения их проходных сечений. Предположим, что уменьшение сечений органов газо­ распределения происходит пропорционально уменьшению диаметра цилиндра (поршня) D. В первом приближении при неизменной конструкции органов газораспределения

1™ = 1кл= const.

Тогда путем несложных преобразований получим

2

(173)

104

Предел повышения числа оборо­

ffr

 

 

 

 

 

тов вала должен

определяться для

 

 

 

 

 

 

каждого конкретного

типа

маши­

8

 

 

 

 

 

ны ‘. Изменение потерь в воздушном

 

 

 

 

 

 

детандере высокого давления с па­

 

>

 

 

 

раметрами:

рвх =

19,6

МН/м21

 

 

/

 

(200 КГС/СМ2) ,

Гвх = 293

К,

Р в ы х

=

 

 

 

 

 

= 1,47 МН/м2

(15 кгс/см2),

Гвы*

=

 

 

 

* —

 

_== 170 К, расход воздуха 64,6 кг/ч,

л » .

 

 

 

 

п = 285 об/мин,

в зависимости

от

 

 

 

 

 

 

 

 

 

скорости вращения

вала

машины,

Опн

 

2п„

 

 

 

подсчитанными

по

вышеприведен­

 

Зпн

пн,°й /н и н

ным

зависимостям,

 

показано

на

 

 

 

 

 

 

рис.

35.

 

 

 

 

 

 

 

Рис.

35.

Изменение

прира­

Проведенный анализ показывает,

щения температуры

возду­

что при малых значениях

удельных

ха

за

детандером

в

ре­

потерь от гидросопротивлений и вы­

зультате

дросселирования

(бГд), внутреннего тепло­

соких от регенеративного

теплооб­

обмена

(6ГТ) и теплоприто-

мена и теплопритока извне увеличе­

ков

(67^ ) при

изменении

ние числа оборотов до определенно­

числа оборотов вала в ми­

го предела даже

при

сокращении

нуту; бГ2 — суммарное

из­

проходных сечений органов газорас­

менение

потерь;

пя — номи­

пределения приводит к повышению

нальное число оборотов

 

термодинамической

эффективности

 

 

 

 

 

 

детандера.

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

мендованных для машин сп — 200 -г- 400 об/мин.

а Вып на

Влияние углов опережения впуска

а вп и

выпуска

к. п. д. детандера значительно слабее,

чем

влияние

скорости

вращения вала машины. Экспериментальная зависимость к. п. д. от углов опережения была получена для детандеров с внутрен­ ним приводом клапанов Грачевым А. Б. [31] и Бумагиным Г. И. Результаты обработки зависимости адиабатического к. п. д. т]ад от углов опережения показали, что правильный выбор зна­ чений а вп и а Вып могут повысить Т1ад машины на 3 5 % . Одно­ временно проведенное осциллографирование по давлению и температуре позволило объяснить падение т)ад при превышении

1 Отметим, что решающим фактором, сдерживаю­ щим повышение числа оборотов вала детандера, является снижение работоспособности органов газораспределения. Это относится прежде всего к машинам с клапанами.

105

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

Аналогичная картина наблюдалась и при исследовании зависимости т)ад от угла опережения выпуска аВьш- С увеличени­ ем аВЬш возрастали потери от неполноты расширения, а при малых значениях авьш процесс выхлопа не успевал закончиться к началу обратного хода поршня. Оптимальные значения углов

опережения авп и авып для детандера с рвх — 1,2 Мн/м2 и п =

=

500 об/мин лежали

в пределах:

авп = 12-=- 16° и аВЫл =

=

16

-г- 18°. Для детандера с рвх = 2 -^ 7 Мн/м2 и л = 385 об/мин

авп =

15 -=- 20° и аВып =

18 25°.

 

 

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

опережения на другие детандеры

невозможно. Приведенные

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

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

к изменению других параметров машины:

вредного простран­

ства, углов опережения и пр. В этом случае требуется

кон­

кретный анализ в каждом частном случае,

подобно тому,

как

это было сделано при изменении числа оборотов.

 

Глава V. МЕТОДЫ ТЕРМОДИНАМИЧЕСКОГО РАСЧЕТА КРИОГЕННЫХ ПОРШНЕВЫХ ДЕТАНДЕРОВ

Исходные данные. Конструктору машины должны быть заданы следующие данные: давления на входе и выходе (рвх и рвых) ; температура газа перед детандером (Гвх); род газа и массо­ вый часовой расход газа через детандер (G4)- Это необходимый и достаточный объем исходных данных. В некоторых случаях вместо Гвх задается температура газа за детандером ГвыхИногда в числе заданных величин указывается необходимое га­ рантированное значение адиабатного к. п. д.

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

Этапы расчета. Термодинамический расчет детандеров под­ разделяется на конструктивный и поверочный. Конструктивный расчет выполняют в начальной стадии проектирования. При этом выбирают тип детандера, тип уплотнения, скорость пор­ шня. Выбирается тип органов газораспределения и их привода. Значения параметров индикаторной диаграммы, гидравлические потери в органах газораспределения, вероятное значение адиа­ батного к. п. д. и другие параметры выбирают исходя либо из опытных данных, либо на основании теоретического исследо­ вания.

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

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

107

рого определяют действительный весовой расход рабочего тела.

Для приведения расчетной производительности детандера к заданной корректируют число оборотов детандера. Кроме того, определяют термодинамические параметры газа в харак­ терных точках индикаторной диаграммы и на выходе из маши­ ны. Оценивают различные потери холодопроизводительности. Вычисляют адиабатный к. п. д. и мощность детандера. В целом проверочный расчет проводят для оценки эффективности рабо­ ты машины, выявления необходимых реконструктивных меро­ приятий и оптимизационных расчетов, выбора типа тормозной системы, определения температур материала основных деталей, получения исходных данных для динамического и прочностного расчета машины и для выбора типа смазки.

Расчетное задание для поверочного термодинамического расчета должно дополнительно содержать следующие све­ дения:

1)чертежи цилиндровой группы детандера с данными о раз­ мерах и конструкции цилиндра и поршня детандера, типе поршневого уплотнения, числе оборотов детандера, мертвом объеме;

2)чертежи общего вида детандера с необходимыми сведе­ ниями о его компоновке и взаимном расположении холодной и

теплой зон машины;

наполнения

3) фазы газораспределения (величины отсечек

и выталкивания), тип и основные конструктивные

параметры

органов газораспределения.

 

Поверочный термодинамический расчет завершает проект и является его неотъемлемой частью.

Методы конструктивного расчета. Для конструктивного рас­ чета детандера необходим выбор параметров индикаторной диаграммы, средней скорости поршня и некоторых других вели­ чин. Выбор этих параметров находится в прямой связи с объемом накопленной информации об оптимальных их значе­ ниях для данного типа поршневого детандера. Для детандеров разных типов (см. гл. I, п. 1) объем такой информации разли­ чен. Наиболее просто можно рассчитать детандеры, по которым накоплены исчерпывающие экспериментальные сведения. Оче­ видно, что в этом случае расчет должен использовать методы подобия.

Есть типы детандеров, степень изученности которых недостаточна для применения методов моделирования, но для которых можно задаться вероятными значениями г)ад, а0, Ь0, с0. Наконец, всегда вероятна ситуация, когда конструктору пред­ стоит рассчитать машину такого типа, по которому эксперимен­ тальные данные практически отсутствуют. Таким образом, объем накопленной информации определяет выбор того или иного метода расчета.

108

1. Метод расчета детандеров, по которым накоплен обширный экспериментальный материал. (Метод моделирования)

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

6 = По “Ь Со •

(174)

На рис. 36, б даны значения б, обеспечивающие

максимальный к. п. д.

в зависимости от давления на входе [61]. Кроме того, для детандеров этого

типа

определена экспериментальным путем величина безразмерного

комп­

лекса:

 

 

 

G

(175)

 

— = 0,85- ■0 ,95(в среднем 0,9)

 

6

 

где

0 = — ---------- ;

 

 

" опчРвх

 

Уоп-ч — описанный часовой объем в м3/ч.

 

Диаметр цилиндра находится так:

 

 

мм,

(176)

где С„ — средняя скорость поршня в м/с.

Средняя скорость поршня может составлять величину 2—3 м/с_ для детандеров с металлическими кольцами со смазкой. Число оборотов п для де­

тендеров этого типа обычно находится в пределах

150—350 об/мин (редко

до 450 об/мин).

 

Ход поршня

 

30С„

(177)

S =

Отношение S/D может быть в пределах 1—4. Наивыгоднейшая степень

расширения в процессе выхлопа для этих машин:

а34 = - ^ - = 2 ч - 3 ,5 .

Р*

до

для детандеров высокого давления (а) и значения степени наполнения бо при максимальном значении к. п. д. (б) [61]

109

На практике

значения б

можно принимать

несколько

больше

значений

б0, определенных

из условия

максимального к. п.

д. (рис.

36). Это

объяс­

няется тем, что при увеличении б сокращаются размеры цилиндра и умень­

шаются

потери от

трения

и

теплопритоков.

В полном смысле оптимальные

значения

б можно

найти

при решении оптимизационной

задачи, постановка

которой дана в п. 4

гл II.

 

метода расчета,

основанного

на моделировании,

В связи с обсуждением

необходимо уточнить условия подобия и критерии подобия для поршневых

детандеров. Если исходить из

системы

уравнений,

описывающей

рабочий

процесс

в

поршневом

детандере (см.

гл.

IV),

и

составить определительные

уравнения

(по А. А. Гухману),

то можно

установить

следующие

Q

критерии

подобия

 

f __f

(аналог

адиабатного

к. п. д .);

а0, Ь0, с0,

 

(или

5------ ——

"3

_

 

 

т в х

 

 

 

 

 

 

 

 

 

 

/

у

в

 

ОТНОШеНИе S/D)', Ар^/Рвых'. Арбы/Двых!

 

k\ ffl\ tl\ бцикла/Еорвх

(или

G4/V оп.ч,

Рвх). Как видим, таких критериев

всего

двенадцать

и полное

по­

добие возможно только в том случае, когда

для

разных

машин

численные

значения этих критериев одинаковы.

Понятно,

что

на

практике

обеспечить

одинаковые значения ряда из этих величин

(например

m,

п, Ap W P bhx)

не­

возможно.

Поэтому подобие не

является

вполне

строгим

и

применимость

этого метода расчета ограничена. Ориентировочно можно считать применение такой методики достоверной, если размеры проектируемой машины и разме­ ры экспериментально изученных машин, при исследовании которых получены значения безразмерных параметров, различаются не более чем в 5—7 раз.

Пример. Рассчитать азотный двухклапанный детандер классического типа.

Известно:

рв1 = 19,6

МН/м2 (200 кгс/см2); р вЫх

= 0,147 МН/м2 (1,5 кгс/см2);

G , = 600

кг/ч; Твх =

300 К. Расчет. Принимаем

значение б = 0,3 (несколько

 

 

 

G

больше, чем значение бо, соответствующее максимуму к. п. д.). Значение ~

принимаем равным 0,9.

Тогда безразмерная величина расхода:

 

 

 

 

 

~G =

0,9

б =

0,27.

 

 

По

диаграмме

pv р для азота

находим удельный

объем при

парамет­

рах на входе и определяем значение плотности рВх-

 

 

 

 

 

 

Рвх —

 

1-

= 2 20

кг/м3.

 

 

 

 

 

 

 

 

 

 

Из

формулы

(175)

находим

величину

описанного часового объема:

 

 

 

 

G„

 

 

600

10 м3/ч.

 

 

 

 

 

GpBX

=

 

 

 

 

 

 

 

0,27-220

 

 

По

формуле

(176)

вычисляем

диаметр

цилиндра,

принимая

среднюю

скорость поршня Сп =

 

2 м/с:

 

 

 

 

 

 

 

D = 2 6 ,6

V

= 26,6 V ю = 59,5

мм.

 

Округляем значение D до 60 мм. Исходя из практических данных о чис­

ле оборотов вала детандера в минуту, при котором обеспечивается надежная работа клапанного привода и допустимая величина депрессии в клапанах, принимаем:

п = 300 об/мин.

Тогда ход поршня по формуле (177):

ЗОС,

30-2

= 0 , 2 м = 200 мм.

S = — =-*

300

п

 

ПО

S

200

e. лежит в допустимых пределах 1—4.

Отношение ~ ~ =

„„ = 3 ,3 4 ,

D

60

 

Фазы газораспределения для детандера этого типа определяются следующим

образом. Углами

опережения

впуска

и

выпуска

задаются:

а вп ~ 5-н8°;

ввып = 8 -Н 10°. Углы

наполнения

и выталкивания

вычисляются

из следую­

щих соотношений для кривошипно-шатунного механизма:

 

 

1

/ ,

 

 

г

sin2 а .

 

 

с„ = —

1 — cos а,. + - —

 

 

0

2

V

 

0

21

 

с

 

(178)

,

,

1

Л

 

 

 

г

 

— cos а 1 „ +

*1 - ь ,) >

 

1

ь0 = —

1 1

21

 

г

где I — отношение радиуса кривошипа к длине шатуна.

2. Метод расчета детандеров при известных наилучших значениях параметров Оо; &о; Со и известном вероятном к. п. Д. Т}ад

Описанный объем цилиндра может быть найден по одной из формул:

G4v3 ■106

Ко = -

 

см3;

(179)

60ni

1 + < J0 —

^5vs ( ° o 4 - bo)

 

 

G,ZV 106

(180)

 

 

c m 3;

60tii

 

V2 ,

 

Co + c0— Vs (a0 + b0)

 

здесь i — число цилиндров.

 

v2, v3 и v3 необходимо

знать темпе­

Для определения удельных объемов

ратуры газа в характерных точках индикаторной диаграммы, которые можно

определить расчетным путем или принять исходя из опыта. Выражения

(179)

и (180) легко представить в виде:

 

 

 

 

 

Уо =

G4z3R - 106____________

(181)

1 + ар

г3

 

 

c m 3;

/

Qp + b0 \

 

60nip3

T3

z3T5

 

a34

)

 

\

 

 

V0 =

G4z2R- 106____________

(182)

a o 4" Co

z 2

 

 

CM3

60nip2

Co 4- b 3 \

 

T 2

2 ST

 

5a

/

 

 

 

 

где 2 — коэффициент сжимаемости.

 

 

от

точности определения

Точность вычислений зависит прежде всего

температур Т2 и Т3 или Т3 и Г5. Для оценки

того или иного метода опреде­

ления температур в табл. 13

приведены

значения ошибок в величине описан­

ного объема ДУ0 для гелиевого детандера при

определении температур

раз­

ными способами. Величины а0, bo, с0, G4, Твх, рвх и п являются значениями

действительных параметров реальной машины. Как видно из табл. 13, для определения описанного объема детандера при величине утечек до 4% можно рекомендовать формулу (182), с определением температур Т2 и Т3 по урав­

нениям (129), (120), (123) и (130) (графа 3 в табл. 13). Поскольку в этом случае расчетное значение Ко несколько меньше, то следует принять действи-

111

Т а б л и ц а 13

Сравнение расчетных значений описанных объемов гелиевого детандера1 (У0)р при определении температур в точках 2, 3 л 5 разными способами

О.

о

с

о

с

2

1

Способ определения температур и допущения

ТВ = Т4 — Твых; Твых— определяется по величине ве­

роятного

адиабатного к.

п. д.

с помощью

тепловых

диаграмм.

Р 5 — р вых>значением а 34 =

задаются.

Вы­

числяется

Дi34 = z3RT3

(/ —

1 /о 34),

при

этом

Т3

варьируется до тех пор, пока ее значения, найденные из тепловых диаграмм по величинам Ai34, не совпа­

дут с принятыми при вычислении Ai34

Номер

Расчетное

Величина

 

формулы

ошибки 6VV=

 

для

значение

(Vo)p~(Vо)д„

Примечания

опреде­

(Vo)p

ления

 

 

СМ3

(^о)д

 

0 '«)р

 

 

 

X 100%

 

(181)

80,6

+ 14

Метод К. С. Буткевича (под­

 

 

 

робнее см. [7, 46]); для опреде­

 

 

 

ления Ai34 следует пользовать­

 

 

 

ся формулой, приведенной в

 

 

 

данной таблице)

2 т 2 ~ Тъх' Т 5 = Твых; Твых-определяется по вели- (182) 57 - 1 7 ,7

чине вероятного значения г)ад

Р2 ~ Pr t Р 5 = ^вых

3

4

Т2 и Т5 определены

по формулам (129), (120), (123)

(182)

63,2

8 ,8

и (130)

при m = n — k. (В данном случае Т2 = 29 К;

 

 

 

Г5 = 1 1

К)

 

 

 

 

Р2 = Ръх— кРб12

 

 

 

 

Т2 и Т5 определены

по формулам (129), (120), (123)

(182)

73,6

+ 6 , 2

и (130)

при n = k и m = f( Т)ад)

 

 

 

P2 = Pm ~~P6i2

Вероятным значением к. п. д. для вычисления температур за­ даваться не надо

Вероятным значением к. п. д. задаваться не надо

1 Опытные данные для этого детандера: О

= 14 кг/ч;

(V0)„— 6 9 ,3 см»; л = 300 об/мин; 7 ^ = 2 8 К;

78%;

а0—0,079;

Ь0 - 0 ,3 4 4 ; Со - 0,249; а.щ тах

0 ,1 3 3 ; р 2 -

2 ,3 4 ; р3

~ 0 ,5 2 ;

р5 - 0 ,1 9 МН/м*; 7 /2 = 2 7 ,4 ; Т 3 ~ 17,45; T s - 11.45;

Т 6 -

3 5 ,8 К.