CC BY
85
Расчет и конструирование
621.432.98
ПОРШНЕВАЯ ТЕПЛОВАЯ МАШИНА ДЛЯ ПРЕОБРАЗОВАНИЯ ЭНЕРГИИ И ПРОИЗВОДСТВА РАБОТЫ НА ТРАНСПОРТЕ
Б.А. Шароглазов, А.В. Гофман
На базе положений равновесной термодинамики, применительно к конструктивной схеме машины-двигателя по патенту РФ № 2460890, даётся характеристика метода расчёта энергетических параметров такой машины, отличительной особенностью которой является непосредственное преобразование возвратно поступательного движения поршня во вращательное грузового вала. Анализируется численная характеристика достоинств реализуемого в ней рабочего цикла.
Ключевые слова: поршневая тепловая машина, рабочее тело, рабочий цикл, параметры цикла, степень сжатия, кинематика двигателя, индикаторная диаграмма.
Введение
По состоянию на сегодня одной из важных составляющих энергетических ресурсов на транспорте являются жидкие углеводородные топлива, производимые из нефти, природные запасы которой ограничены. В балансе потребления топливно-энергетических ресурсов (ТЭР) на долю транспорта и транспортных систем приходится 25-30 %. В этой связи важной проблемой современной энергетики, наряду с решением задач по альтернативным источникам энергии, является проблема режима расходования ТЭР, относящихся к числу исчерпаемых. Сегодня основным средством преобразования заключенной в топливе энергии и производства механической работы, в том числе и на транспорте, является поршневая тепловая машина с кривошипно-шатунным механизмом (КШМ). Однако двигатели с КШМ наряду с достоинствами обладают и рядом недостатков. К одним из них относятся сложность и высокая стоимость этих механизмов (20-30 % от стоимости двигателя в целом).
Названные обстоятельства побуждают исследователей к отысканию новых схем и типов машин-двигателей, свободных от указанных недостатков. Примером таких решений может служить ряд конструктивных схем машин, предложенных различными авторами [1-4]. Важной особенностью таких машин-двигателей является непосредственное преобразование возвратно-поступательного движения поршня во вращательное грузового вала. Использование механизма преобразования движения позволяет исключить из конструктивной схемы двигателя шатун и кривошип. Благодаря этому в схеме сил, действующих на поршень, удаётся исключить так называемую нормальную силу, действующую перпендикулярно поверхности зеркала цилиндра и обуславливающую рост затрат на преодоление трения в сопряжении «поршень-цилиндр». Это обстоятельство (снижение трения) обуславливает повышение механического и эффективного КПД бескри-вошипной машины (БКМ).
Поскольку закон движения поршня в таких двигателях отличен от закона, обеспечиваемого КШМ, при их создании необходимо использование иного расчётного инструментария.
Конструктивные особенности и рабочий цикл БКМ
Особенности расчета БКМ вытекают из анализа конструктивной схемы двигателя, в котором механизм преобразования движения вынесен за пределы поршня (рис. 1). В нём возвратнопоступательное движение поршня 1 посредством штока 2 передаётся каретке 3. В каретке установлена ось 4, цапфы которой перемещаются по криволинейному замкнутому (например синусоидальному) пазу (беговой дорожке), образованному криволинейными торцевыми поверхностями втулок 5 и 6. Движение цапф по беговой дорожке обуславливает вращательное (вокруг оси двигателя) движение оси 4 и каретки 3. Каретка в плоскости, перпендикулярной оси 4, охвачена
стержнями двухрожкового грузового вала 7. Вращательное движение каретки, благодаря такой кинематической связи, передается грузовому валу 7.
Беговая дорожка, образованная криволинейными профилями торцевых поверхностей втулок 5 и 6, в общем случае может быть многопериодной. Например, может иметь профиль двухпериодной синусоиды. В этом случае одному полному обороту грузового вала соответствует четыре хода поршня (два возвратно-поступательных). Таким образом, если иметь ввиду схему реализации 2-тактного цикла, то за один оборот грузового вала таких циклов реализуется два. Следовательно, применительно к прочим равным условиям, работа, производимая тепловой машиной (например, транспортным двигателем) выполненной по приводимой схеме, в два раза может превысить работу, производимую за такой же отрезок времени, традиционной машиной с КШМ.
Рис. 1. Конструктивная схема бескривошипной тепловой машины-двигателя (применительно к случаю реализации 2-тактного рабочего цикла)
Рабочий цикл бескривошипной тепловой машины-двигателя
Характер изменения параметров состояния рабочего тела (РТ) в цилиндре двигателя в течение цикла определяется многими факторами. В частности: особенностями подвода (или отвода) теплоты к РТ, динамикой выгорания топлива, кинематикой движения поршня, режимом нагрузки.
Кинематика движения поршня массово распространённых машин с КШМ описывается широко известными соотношениями (по этой причине здесь эти соотношения не приводятся).
Применительно к условиям реализации рабочего цикла (РЦ) в бескривошипной машине рассматриваемой схемы с двухпериодной беговой дорожкой кинематические особенности движения поршня могут быть отображены уравнениями:
5 = • (1- о^2а) 1
V = ы • 5П • sin2а К (1)
_/ = 2 • ы2 • 5П • cos2а)
где 5, V, ] - текущие значения перемещения, скорости и ускорения поршня соответственно; а -текущий угол поворота грузового вала; «п - полное перемещение (ход) поршня; ю - угловая частота вращения.
Приведённые кинематические соотношения (1) для БКМ и соответствующие широко известные соотношения для машин с КШМ позволяют сделать (применительно к числено равным значениям полных перемещений «п поршня) сопоставление текущих параметров, характеризующих кинематику этих машин (рис. 2).
Ещё раз отметим то обстоятельство, что в БКМ с двухпериодной беговой дорожкой за один оборот грузового вала совершается четыре хода поршня, а в двигателе с КШМ - два. Поэтому сопоставление кинематических параметров машин целесообразно выполнить, учитывая период изменения этих параметров. Для двигателя с КШМ он составляет 360 град., для БКМ - 180. Иными словами 5, V, ] для названной цели удобно представить функциями относительного времени £ = — = —, где гц, ац - соответственно временная и фазовая (по углу поворота вала) продол-
£ц я ц
жительность цикла. Графики, приведённые на рис. 2, выполнены с учётом названных обстоятельств. Анализ кинематических параметров машин (см. рис. 2), позволяет отметить, что характер их протекания подобен. Экстремальные значения 5, V,] приведены в табл. 1.
Рис. 2. Характер протекания кинематических параметров двигателей в функции относительного времени при {ц = 6,3 • 10_3, с
Таблица 1
Численные значения параметров кинематики в экстремальных точках
s, м v, м/с j, м/с2
t 0 0,5 0,25 0,208 0,5 0 0,5
КШМ 0,034 0 -1б,91 -17,43 - макс. 0 -21199 12450
БКМ 0,034 0 -1б,91 - макс. - 0 -1б824,93 1б824
Моделирование параметров состояния РТ в процессе сжатия
Функцию для текущего перемещения поршня 5 в (1) перепишем в виде:
7' (1- о^2а) = ^ • а(а), (2)
здесь а(а) - кинематическая функция хода поршня (численное значение её определяется, в нашем случае, только текущим а).
По аналогии с "!еЬе [5] используя (2), запишем функцию, характеризующую изменение объёма РТ в цилиндре:
ф(а) = 1 + ■ а (а), (3)
в которой е - степень сжатия.
Тогда система уравнений для определения параметров состояния рабочего тела (давления р и температуры Т) в процессе сжатия запишется в виде:
1«1
«1-і
P(a)=M^gf]
T(a) = Та • [^]"‘
ф(а) = 1 + ^ • о(а), о(а) = 1 — cos2a,
(4)
где n1 - показатель политропического процесса.
Исходными условиями для решения системы являются численные значения параметров в начальной точке (точка «а» индикаторной диаграммы - положение в НМТ): p = pa; T = Ta; y(a) = y(aa); o(a) = o(aa). Параметры состояния РТ в процессе сжатия применительно к двигателю с КШМ отображаются системой, аналогичной (4) [6]. Но в ней для определения значений a(a) используется уравнение
a(a) = 1+ ^ — (cosa + ^ • Vl—A2 • sin2aj, где A = r/Lm (г - радиус кривошипа, 1ш - длина шатуна).
Моделирование параметров состояния РТ в процессе сгорания
Параметры состояния РТ в процессе сгорания (прежде всего р и Т) являются функцией многих факторов. В этой связи определение их текущих значений в этот период РЦ выполняется пошагово на базе методов численного интегрирования. Нами использован метод расчёта, предложенный И.И. Вибе [5], достаточно подробно изложенный в [6].
Важной составляющей метода является расчёт количеств теплоты, сообщаемых РТ вследствие горения топлива на каждом достаточно малом шаге расчёта. При этом для описания характера выделения теплоты используется уравнение выгорания Вибе:
в котором /, 4 - текущее время и продолжительность сгорания соответственно; т - показатель характера сгорания.
Для целей сопоставительных исследований (оценки качеств РЦ двигателя с КШМ и БКМ) необходимы численные значения параметров т и ф2, а также 0 (угол опережения воспламенения). Это делалось на основании следующего. Полагалось, что характер горения углеводородных смесей определенного качества (например, характеризуемых одинаковыми значениями коэффициента избытка воздуха а) при близких численных значениях е в 2- и 4-тактных двигателях во времени одинаков. Тогда значения т и 4, определённые для 4-тактных двигателей, можно распространить и на 2-тактные. Так, по результатам обобщения экспериментального материала по бензиновым 4-тактным двигателям, работающим на частотах вращения 3500...5600 мин-1 с а = 0,85...0,95, установлено, что ф2 находится в пределах 35.50 град. ПКВ, т = 3.3,5. Учитывая связь временного интервала и соответствующей ему фазы поворота вала (ф = 6- п • £), определили, что продолжительность горения смеси названного качества находится в пределах 1,5 • 10_3 ... 1,6 • 10_3 с. Соответствующие значения 0 отвечают временным интервалам горения, находящихся в диапазоне 0,7-10_3...0,35-10_3 с. Сказанное даёт основания для оценки численных значений 0 и фг применительно к двигателям, работающим по 2-тактному циклу. Это, в свою очередь, позволяет расчётным путём выявить численные значения параметров, характеризующих состояние рабочего тела в процессе горения.
Отметим, что начальными условиями для расчёта горения являются параметры РТ, соответствующие завершению сжатия (момент, определяемый значением 0: точка у с соответствующими ру, Ту, уу).
Процесс расширения
Моделирование параметров в процессе расширения базируется на условии, что этот процесс является политропическим. Начальными условиями для расчёта являются численные значения параметров РТ в момент окончания сгорания (точка ъ на индикаторной диаграмме): р = рг; Т = Тг; у(а) = у(аъ); а(а) = о(аъ). Таким образом, система уравнений для расчёта текущих параметров имеет вид (применительно к БКМ):
а(а) = 1 — cos2а. 7
Для моделирования процесса в двигателе с КШМ используется система уравнений аналогичная (6). Но в ней о(а) рассчитывается с учётом замечания, сделанного по поводу этой функции в тексте, касающегося моделирования параметров при сжатии.
Изменение параметров состояния РТ в течение РЦ для БКМ и для двигателя с КШМ иллюстрируется «р-у» - диаграммой, рис. 3. На рис. 4 приведены графики изменения давлений и температур РТ в функции относительного времени. Сопоставительная оценка параметров по некоторым показателям цикла приводится в табл. 2.
(5)
(6)
ро I—I—І------------------------------------г--Т"
0,2 0,4 0,6 0,8 1 1,2 V, м3/кг
Рис. 3. Индикаторные «р-у» диаграммы двигателей с КШМ и БКМ
0 0.25 0,5 0,75 I
Рис. 4. Характер протекания давлений и температур РТ в функции относительного времени
Сравнение основных термодинамических параметров двигателей
Таблица 2
Рмакс, МПа ^рмак^ град. пкв Т макс ^Тмак^ град. пкв Р„ МПа Пі gi, г/(кВт ч)
КШМ 4,19 10 2565 30 0,686 0,37 220
БКМ 4,24 5 2604 15 0,67 0,36 225
Расчётная оценка показателей РЦ бескривошипной поршневой тепловой машины и машины-двигателя с КШМ применительно к одинаковым условиям реализации РЦ (по величине е, а, т, tz, Ь, Би) даёт основания утверждать, что бескривошипный двигатель по таким параметрам цикла как рмакс, Тмакс, ри п не уступает поршневому двигателю с КШМ. При этом производимая им работа, приходящаяся на один оборот грузового вала, превышает (почти в два раза) работу, производимую соответствующим двигателем с КШМ.
Литература
1. Daisaku Odawara, No. 520 Ueno Shiba-cho 5-cho, Sakai-shi, Japan. Reciprocating Machine. US3598094, Appl. No 18.794; Filed Mar. 16, 1970; Patented Aug. 10, 1971; Priority Apr. 28, 1967, June 10, 1967; Japan; 42/26889 and 42/36781 Continuation of application Ser. No. 717,518, Apr. 1, 1968. now abandoned. U.S. Cl. 123/58 R, 123/43 A, 123/58 B, 123/58 C, 74/60, 91/205, 91/411; Int. Cl. F02b 75/26, F01b 13/06, F16h 33/00.
2. Hiromassa Kitaguchi, Chi-l0l, Simamachi, Komatsu-shi, Isikawa-ken, 923-03, Japan. Crancless Engine. US5007385, Apr. 16, 1991, Appl. No.: 501,712; Filed: Mar. 30. 1990 Int. Cl. F02B 75/04. US. Cl. 123/48 B; 123/58 BA; 123/78 E.
3. Пат. 2156871 C1 Российская Федерация, МПК 7 F 02 B 75/26, F 01 B 9/06. Бескривошип-ный двигатель внутреннего сгорания / Б.А. Шароглазов, В.Г. Ваганов. - № 99110007/06; заявл. 18.05.1999; опубл. 27.09.2000, Бюл. № 27.
4. Пат. 2460890 C1 Российская Федерация, МПКF02B 75/26 (2006.01), F02B 75/32 (2006.01), F01B 9/06 (2006.01). Бескривошипная поршневая тепловая машина-двигатель /Б.А. Шароглазов, С.Ю. Терехов, И.И. Колбин. - № 2011128944/06; заявл. 12.07.2011; опубл. 10.09.2012, Бюл. № 25.
5. Wiebe, I.I. Brennverlauf und Kreisprozefi von Verbrennungsmotoren / I.I. Wiebe. - Berlin: Veb-verlagtechnik, 1970. - 275p.
6. Шароглазов, Б.А. Поршневые двигатели: теория, моделирование и расчёт процессов / Б.А. Шароглазов, В.В. Шишков. - Челябинск. Издат. центр ЮУрГУ, 2011. - 525 с.
Шароглазов Борис Александрович. Доктор технических наук, профессор, заслуженный деятель науки Российской Федерации, профессор кафедры «Двигатели внутреннего сгорания», Южно-Уральский государственный университет (Челябинск), [email protected].
Гофман Александр Викторович. Аспирант кафедры «Двигатели внутреннего сгорания», Южно-Уральский государственный университет (Челябинск), [email protected].
Поступила в редакцию 21 января 2014 г.
Bulletin of the South Ural State University Series “Mechanical Engineering Industry” _____________2014, vol. 14, no. 1, pp. 5-11
PISTON HEAT ENGINE FOR ENERGY CONVERSION AND TRANSPORT WORK PERFORMANCE
B.A. Sharoglazov, South Ural State University, Chelyabinsk, Russian Federation, [email protected],
A.V. Gofman, South Ural State University, Chelyabinsk, Russian Federation, [email protected]
The report based on principles of equilibrium thermodynamics contains specification of method of energy characteristic calculation for piston heat engine (patent № 2460890) that is featured with direct conversion of reciprocating motion to rotational motion of the main shaft. Also, the report provides numerical characteristic of working cycle advantages of such engine.
Keywords: piston heat engine, working medium, working cycle, cycle characteristics, compression ratio, engine kinematics, indicator diagram.
References
1. Daisaku Odawara, No. 520 Ueno Shiba-cho 5-cho, Sakai-shi, Japan. Reciprocating Machine. US3598094, Appl. 18.794; Filed Mar. 16, 1970; Patented Aug. 10, 1971; Priority Apr. 28, 1967, June 10, 1967; Japan; 42/26889 and 42/36781 Continuation of application Ser. 717,518, Apr. 1, 1968. now
abandoned. U.S. Cl. 123/58 R, 123/43 A, 123/58 B, 123/58 C, 74/60, 91/205, 91/411; Int. Cl. F02b 75/26, F01b 13/06, F16h 33/00.
2. Hiromassa Kitaguchi, Chi-l0l, Simamachi, Komatsu-shi, Isikawa-ken, 923-03, Japan. Cranc-less Engine. US5007385, Apr. 16, 1991, Appl. 501,712; Filed: Mar. 30. 1990 Int. Cl. F02B 75/04. US. Cl. 123/48 B; 123/58 BA; 123/78 E.
3. Sharoglazov B.A., Vaganov V.G. Beskrivoshipnyy dvigatel' vnutrennego sgoraniya [Crankless Internal Combustion Engine]. Patent RF, № 2156871 C1 7 F 02 B 75/26, F 01 B 9/06, Application: 99110007/06, 18.05.1999; Date of publication: 27.09.2000. Bull. 27.
4. Sharoglazov B.A., Terekhov S.Yu., Kolbin I.I. Beskrivoshipnaya porshnevaya teplovaya mashi-na-dvigatel' [Axial Piston Thermal Engine]. Patent RF, № 2460890 C1 Int. Cl. F02B 75/26 (2006.01), F02B 75/32 (2006.01), F01B 9/06 (2006.01), Application: 2011128944/06, 12.07.2011.; Date of publication: 10.09.2012. Bull. 25.
5. Wiebe I.I. Brennverlauf und Kreisprozefi von Verbrennungsmotoren. Berlin, Vebverlagtechnik, 1970. 275 p.
6. Sharoglazov B.A., Shishkov V.V. Porshnevye dvigateli: teoriya, modelirovanie i raschetprotses-sov [Piston Engine: Theory, Simulation and Calculation Process]. Chelyabinsk, South Ural St. Univ. Publ., 2011. 525 p.
Received 21 January 2014
