Курсовая Расчет авиационного поршневого двигателя

Поможем написать учебную работу

Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.

Предоплата всего

от 25%

Подписываем

договор

Имя

Выберите тип работы:

Курсовая работа Дипломная работа Контрольная работа Реферат Отчет по практике Эссе Доклад Решение задач Научно-исследовательская работа Монография Аспирантский реферат Магистерская работа Научная статья Публикация статьи в ВАК Публикация статьи в Scopus Дипломная работа MBA Бизнес-план Тест/экзамен online Чертеж Рецензия

Принимаю Политику конфиденциальности

Продолжить

Скидка 25% при заказе до 28.2.2025

Министерство образования и науки Украины

Национальный аэрокосмический университет

им. Н.Е. Жуковского “ХАИ”


Расчетно-пояснительная записка к курсовому проекту по дисциплине

“Авиационные двигатели”
Харьков 2008


Содержание

1. Тепловой расчет

1.1 Исходные данные

1.2 Выбор дополнительных исходных данных

1.3 Расчет процесса наполнения

1.4 Расчет процесса сжатия

1.5 Расчет процесса сгорания

1.6 Расчёт процесса расширения

1.7 Определение индикаторных параметров двигателя

1.8 Определение эффективных параметров двигателя

1.9 Определение геометрических параметров двигателя

2. Динамический расчет

2.1 Допущения

2.3 Разнос масс КШМ с прицепными шатунами

2.4 Силы инерции

2.5 Построение верхней петли индикаторной диаграммы

2.6 Суммарная сила, действующая на поршень

2.7 Силы, действующие в центральном КШМ

2.8 Суммарные радиальные и окружные силы действующие на шатунную шейку

3. Уравновешивание двигателя

4. Расчет на прочность коленчатого вала

4.1 Силы, действующие на колено коленчатого вала

4.2 Определение запаса прочности в шатунной шейке

4.3 Определение запаса прочности в коренной шейке

5. Расчет поршневого пальца

6. Расчет поршневых колец

Список использованной литературы

1. Тепловой расчет

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

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

1.1 Исходные данные

1) Эффективная мощность на расчетной высоте –.

2) Частота вращения коленчатого вала – .

3) Число цилиндров – .

4) Степень сжатия – .

5) Давление наддува – .

6) Расчетная высота – .

7) Прототип двигателя – АИ-26.

1.2 Выбор дополнительных исходных данных

1. Коэффициент избытка воздуха принимаем равным .

2. Топливо. Сорт применяемого топлива зависит от степени сжатия и давления наддува. В нашем случае (относительно невысокая степень сжатия) наиболее подходящим является бензин Б-91/115.

Низшую теплотворную способность топлива определим по формуле
.
Элементарный состав авиабензина Б-91/115 включает в себя: ; ; , средняя молекулярная масса – . Тогда:

3. Параметры воздуха на расчетной высоте.

1.3 Расчет процесса наполнения

Цель расчета процесса наполнения – определение давления  и температуры  свежего заряда в конце хода выпуска.

1. Согласно заданию давление наддува . Находим температуру воздуха в конце такта наполнения
,


где  - повышение температуры воздуха в нагнетателе. Адиабатическая работа сжатия 1кг воздуха равна:
.
Адиабатический КПД центробежного нагнетателя примем равным .

Тогда:
;

.
2. Определяем коэффициент наполнения двигателя с наддувом на расчетной высоте
,
где  - приведенный коэффициент наполнения.

Примем , тогда
.


3. Находим давление в конце такта наполнения.
,
где p_r – давление остаточных газов в конце такта наполнения. Принимаем:
.
Степень подогрева свежей смеси в процессе наполнения  условно характеризует результат суммарного теплообмена смеси со стенками цилиндра и донышком поршня, а также понижение температуры за счет испарения топлива.

При ,. Тогда:
.
После подстановки найденных и полученных величин получим
.
4. Определяем коэффициент остаточных газов
,




где  - температура остаточных газов.

Примем , тогда:
.
5. Находим температуру газов в конце такта впуска
.

1.4 Расчет процесса сжатия

Цель расчета процесса сжатия – определение давления  и температуры  газов в конце этого процесса.

1.5 Расчет процесса сгорания

Цель расчета процесса сгорания – определение максимальных значений давления  и температуры  газов при сгорании топлива.

1. Температура  газов определим из уравнения сгорания, полученного на основании первого принципа термодинамики:

,

 - коэффициент эффективного выделения теплоты. Примем ;

 – теоретически необходимое количество воздуха для сгорания  топлива.

.

 - действительный коэффициент молекулярного изменения, где

 -

Для случая  определяем

.

.

 

- средняя молярная теплоемкость газов в интервале температур от 0 до .

.

Расчет температуры  в пакете MathCAD 13

,

1.6 Расчёт процесса расширения

Цель расчёта процесса расширения – определение давления  и температуры  в конце такта расширения.

1.7 Определение индикаторных параметров двигателя

,

где  - коэффициент полноты (скругления) индикаторной диаграммы. Примем .

 - степень повышения давления.

.

1.8 Определение эффективных параметров двигателя

,

где  - коэффициент, оценивающий долю индикаторной мощности, затраченной на привод нагнетателя.

.

.

Среднее давление механических потерь  характеризует мощность, затраченную на преодоление сил трения, на привод вспомогательных механизмов и агрегатов и на “насосные” потери.

Для определения  пользуются эмпирическими уравнениями, полученными на основании экспериментальных данных.

,

1.9 Определение геометрических параметров двигателя

2. Определяем диаметр цилиндра  и ход поршня . Обозначим отношение . Тогда , откуда .

Значение m принимаем по прототипу .

.

3. Ход поршня .

4. Общий рабочий объем двигателя

.

2. Динамический расчет

2.1 Допущения

Ход поршня  и радиус кривошипа  найдены в тепловом расчете.

Основные размеры центрального КШМ вполне определяются радиусом  и длиной шатуна . Отношение  принимаем таким же как и у прототипа, . Тогда длина шатуна:

.

Радиусы прицепов  прицепных шатунов в различных цилиндрах неодинаковы. Из условия геометрического подобия следует, что

.

,

где  - угол между плоскостью симметрии главного шатуна и r_min.

.

, где .

	1	2	3
	0,078	0,79	0,076

2.3 Разнос масс КШМ с прицепными шатунами

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

2. Под “приведенным” главным шатуном (рисунок 3) понимают собственно главный шатун плюс массы пальцев прицепных шатунов  и массы , сосредоточенные на осях этих пальцев. Обозначим; . Приведенный главный шатун заменяем массами М_ПШ, сосредоточенной на оси поршневого пальца, и , сосредоточенной на оси шатунной шейки. Величины  и  определяем из формул:

,

.

,

где  – масса комплекта поршня;

 – часть массы прицепного шатуна, отнесенная к оси поршневого пальца.

,

где - масса вращательно-движущихся частей;

 - часть массы шатуна;

– приведенная масса кривошипа.

2.4 Силы инерции

,

,

где  - ускорения масс

 

.

,

.

Силы , постоянные по модулю, приложены к оси шатунной шейки и направлены по радиусу кривошипа.

2.5 Построение верхней петли индикаторной диаграммы

Согласно принятым ранее допущениям считаем, что в такте наполнения и выхлопа разность абсолютных давлений в цилиндре и картере равна нулю. Абсолютные давления в тактах сжатия и расширения меняются по политропам. Сгорание происходит при постоянном объеме. В конце сгорания давление составляет 0.85 от расчетного. Расширение заканчивается скачкообразным падением давления в НМТ от расчетного  до давления в картере .

,

 

 – полный объем цилиндра,

 – текущий объем над поршнем,

,

где

 – показатель политропы (в процессе сжатия , в процессе расширения ). После подстановки получим:

.

Сжатие			Расширение
180	89267	0,00312	360	1119735	0,00048
190	86819	0,00318	370	862225	0,00058
200	80365	0,00337	380	488627	0,00089
210	72102	0,00365	390	273784	0,00136
220	64068	0,00399	400	167575	0,00196
230	57477	0,00432	410	113280	0,00261
240	52823	0,00460	420	83752	0,00327
250	50251	0,00477	430	66934	0,00386
260	49855	0,00480	440	57259	0,00433
270	51896	0,00466	450	52023	0,00465
280	57000	0,00435	460	49884	0,00480
290	66478	0,00388	470	50197	0,00478
300	82961	0,00329	480	52692	0,00461
310	111869	0,00264	490	57271	0,00433
320	164913	0,00198	500	63797	0,00400
330	268448	0,00138	510	71794	0,00367
340	477973	0,00090	520	80081	0,00338
350	847071	0,00059	530	86642	0,00319
360	1119735	0,00048	540	89265	0,00312

2.6 Суммарная сила, действующая на поршень

,

 

 -

 -

2.7 Силы, действующие в центральном КШМ

,

,

,

,

.

2.8 Суммарные радиальные и окружные силы действующие на шатунную шейку

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

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

На шейку кривошипа дополнительно к  действует в радиальном направлении центробежная сила  вращательно-движущихся масс шатуна .

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

Число однорядных звезд, .

Число прицепных шатунов, .

Частота вращения коленчатого вала, .

Степень сжатия, .

Ход поршня, .

Диаметр цилиндра, .

Отношение радиуса к длине главного шатуна, .

Радиус прицепа шатунов, .

Давление в конце pасшиpения, .

Давление в конце наполнения, .

Атмосферное давление, .

Показатель политpопы сжатия, .

Показатель политpопы pасшиpения, .

3. Уравновешивание двигателя

,

- центробежная сила вращающихся частей равна:

Если исходить из положения, что все шатуны в двигателе центральные, то силы  и

 

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

,

Тогда . - боковой цилиндр.

 

,

,

Тогда в момент  

 

4. Расчет на прочность коленчатого вала

4.1 Силы, действующие на колено коленчатого вала

; ;

, (смотри рисунок 7)

где

Тогда

4) В плоскости колена действует сила инерции от массы шатунной шейки

5) В плоскости колена действуют силы инерции от масс щек _.

7) От сил, действующих в плоскости колена, на опорах колена возникают реакции . В случае симметричного колена

.

8) От силы  на опорах колена возникают реакции , действующие перпендикулярно к плоскости колена. В случае симметричного колена .

Силы Т, Z, реакции  и  переменны по величине и направлению и меняются в зависимости от угла поворота колена. Расчет этих сил реакции приведен в таблице 5.

4.2 Определение запаса прочности в шатунной шейке

;

α	T,Н		Z,Н		,Н		,Н				,		My		,				,
0	0		-2950		0		21119		395		0		197		15		0		0
12,86	-3501		-11940		-1751		21624		400		-228		3		0		-133		-5
25,71	-5900		-9254		-2950		22967		413		-384		-125		-9		-224		-8
38,57	-6502		-5779		-3251		24705		431		-423		-151		-11		-247		-9
51,43	-5256		-2592		-2628		26298		447		-342		-73		-5		-199		-7
64,29	-2712		-551		-1356		27319		457		-176		76		6		-103		-4
77,14	247		-11		124		27588		460		16		244		18		9		0
90,00	2777		-788		1389		27200		456		181		384		28		105		4
102,86	4354		-2344		2177		26422		448		283		469		35		165		6
115,71	4860		-4072		2430		25558		439		316		493		37		184		7
128,57	4494		-5529		2247		24830		432		292		469		35		171		6
141,43	3583		-6526		1792		24331		427		233		415		31		136		5
154,29	2427		-7083		1214		24053		424		158		349		26		92		3
167,14	1210		-7331		605		23929		423		79		280		21		46		2
180,00	0		-7396		0		23896		423		0		211		16		0		0
180,00	0		-7569		0		23810		422		0		211		16		0		0
192,86	-1241		-7518		-621		23835		422		-81		141		10		-47		-2
205,71	-2505		-7313		-1253		23938		423		-163		70		5		-95		-4
218,57	-3747		-6824		-1874		24182		425		-244		2		0		-142		-5
231,43	-4805		-5912		-2403		24638		430		-312		-55		-4		-182		-7
244,29	-5415		-4536		-2708		25326		437		-352		-86		-6		-205		-8
257,14	-5291		-2848		-2646		26170		445		-344		-75		-6		-201		-7
270,00	-4287		-1217		-2144		26986		453		-279		-15		-1		-163		-6
282,86	-2570		-116		-1285		27536		459		-167		85		6		-98		-4
295,71	-696		141		-348		27665		460		-45		191		14		-26		-1
308,57	534		-263		267		27462		458		35		259		19		20		1
321,43	510		-453		255		27367		457		33		257		19		19		1
334,29	-518		812		-259		28000		464		-34		203		15		-20		-1
347,14	-1072		3656		-536		29422		478		-70		179		13		-41		-2
360,00	0		5387		0		30288		486		0		243		18		0		0
360,00	0		57280		0		56234		746		0		373		28		0		0
372,86	17510		59720		8755		57454		758		1138		1365		101		665		25
385,71	24990		39200		12495		47194		656		1624		1734		128		948		35
398,57	24460		21740		12230		38464		568		1590		1661		123		928		34
411,43	21530		10620		10765		32904		513		1399		1468		109		817		30
424,29	19000		3858		9500		29523		479		1235		1309		97		721		27
437,14	17300		-779		8650		27205		456		1125		1202		89		657		24
450,00	15910		-4515		7955		25337		437		1034		1114		83		604		22
462,86	14310		-7704		7155		23742		421		930		1016		75		543		20
475,71	12280		-10290		6140		22449		408		798		895		66		466		17
488,57	9882		-12160		4941		21514		399		642		756		56		375		14
501,43	7318		-13330		3659		20929		393		476		608		45		278		10
514,29	4776		-13940		2388		20624		390		310		464		34		181		7
527,14	2344		-14200		1172		20494		389		152		326		24		89		3
540,00	0		-14270		0		20459		388		0		194		14		0		0
540,00	0		-7396		0		23896		423		0		211		16		0		0
552,86	-1210		-7331		-605		23929		423		-79		143		11		-46		-2
565,71	-2427		-7083		-1214		24053		424		-158		75		6		-92		-3
578,57	-3583		-6526		-1792		24331		427		-233		12		1		-136		-5
591,43	-4494		-5529		-2247		24830		432		-292		-37		-3		-171		-6
604,29	-4860		-4072		-2430		25558		439		-316		-54		-4		-184		-7
617,14		-4354		-2344		-2177		26422		448		-283		-21		-2		-165		-6
630,00		-2777		-788		-1389		27200		456		-181		71		5		-105		-4
642,86		-247		-11		-124		27588		460		-16		216		16		-9		0
655,71		2712		-551		1356		27319		457		176		381		28		103		4
668,57		5256		-2592		2628		26298		447		342		519		38		199		7
681,43		6502		-5779		3251		24705		431		423		581		43		247		9
694,29		5900		-9254		2950		22967		413		384		539		40		224		8
707,14		3501		-11940		1751		21624		400		228		397		29		133		5
720,00		0		-12950		0		21119		395		0		197		15		0		0

,

,

,

.

;

.

, ,

 

Для шеек валов у края смазочных отверстий

.

,

4.3 Определение запаса прочности в коренной шейке

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

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

0,0	0,0	0,0	0,0	0,0	0,0	0,0	0,0
-66,5	-1,1	-23,6	-0,4	332,7	5,7	-23,0	-0,4
-112,1	-1,9	-47,6	-0,8	474,8	8,1	-46,1	-0,8
-123,5	-2,1	-71,2	-1,2	464,7	7,9	-68,1	-1,2
-99,9	-1,7	-91,3	-1,6	409,1	7,0	-85,4	-1,5
-51,5	-0,9	-102,9	-1,8	361,0	6,2	-92,3	-1,6
4,7	0,1	-100,5	-1,7	328,7	5,6	-82,7	-1,4
52,8	0,9	-81,5	-1,4	302,3	5,2	-52,8	-0,9
82,7	1,4	-48,8	-0,8	271,9	4,6	-4,7	-0,1
92,3	1,6	-13,2	-0,2	233,3	4,0	51,5	0,9
85,4	1,5	10,1	0,2	187,8	3,2	99,9	1,7
68,1	1,2	9,7	0,2	139,0	2,4	123,5	2,1
46,1	0,8	-9,8	-0,2	90,7	1,5	112,1	1,9
23,0	0,4	-20,4	-0,3	44,5	0,8	66,5	1,1
0,0	0,0	0,0	0,0	0,0	0,0	0,0	0,0

,

,

5. Расчет поршневого пальца

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


Рисунок 8 – Модель поршневого пальца
При плавающем пальце для каждой точки его расчетного сечения возможен такой случай, когда в момент действия максимального усилия (таковым обычно является ) она окажется либо в зоне максимально растянутых, либо в зоне максимально сжатых волокон. Таким образом, в каждой точке расчетного сечения напряжение может колебаться в соответствии с изменением нагрузки от  до , где .


Рисунок 9 – Схема нагружения поршневого пальца


Этот случай симметричного цикла нагружения опасен с точки зрения прочности поршневого пальца и поэтому является расчетным.
Изгибающий момент
 
(где L – расстояние между серединами опорных участков в бобышках поршня, м; а – длина верхней головки шатуна, м; Р_с – сила, действующая на палец со стороны поршня, Н, изменяем в пределах от  до , где
, где
Масса поршня

Значение  в зависимости от угла поворота коленчатого вала приведены в таблица 7.
Таблица 7 – Значение  в зависимости от угла поворота коленчатого вала

0,00	-11831,50	51,43	-5235,48	102,86	4348,29	154,29	6791,70
12,86	-11347,48	64,29	-2456,71	115,71	5609,74	167,14	6775,33
25,71	-9958,18	77,14	209,64	128,57	6356,74	180,00	6756,88
38,57	-7835,03	90,00	2529,91	141,43	6700,84

Тогда
Среднее напряжение рассматриваемого цикла равно нулю, а амплитуда цикла . Момент сопротивления пальца изгибу
.


Запас прочности пальца на изгиб определяется по выражению
,
где коэффициент фактора размерности  для диаметра пальца 25мм составляет 0,92.

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

Величина  составляет 450 МПа. Тогда





Рисунок 10 – Схема овализации поршневого пальца
Так как по этой формуле запас прочности поршневого пальца значительно занижен по сравнению с действительным запасом, то полученный запас прочности является достаточно надежным критерием оценки прочности.

Под действием силы Р_С палец овализируется (рисунок 10). В направлении действия Р_С диаметр пальца уменьшается, в противоположном направлении – увеличивается. Как показали исследования Р.С. Кинасошвили, наибольшие напряжения при овализации возникают в точках  и  и их величина
,
где ξ=f(α), а α, в свою очередь равна




.
Графическая зависимость ξ=f(α) приведена в пособии [4]. По ней определяем .

Тогда

Так как допустимое напряжение овализации лежит в пределах 230-270МПа [4], то запас прочности по овализации обеспечивается.

Степень овализацци задается обычно как увеличение внешнего диаметра в направлении, перпендикулярном приложению нагрузки. Это увеличение определим по формуле Кинасошвили [4]:

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

Удельное давление пальца на бобышки поршня:
, где b – длина контактирующей с пальцем поверхности бобышки



Тогда как допустимое удельное давление пальца на бобышки поршня составляет 50 МПа [4].

6. Расчет поршневых колец

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

Если принять, что средний радиус кольца в свободном состоянии равен радиусу наружной поверхности поршня (рисунок 10), то формула
 
с учетом того, что  обычно в 20-25 раз меньше R принимает вид
.
Для чугунных колец Е=82500МПа. Тогда



Рисунок 10 – Схема расположения и геометрические размеры кольца


Допускаемые напряжения сжатия составляют 135-150МПа, поэтому запас прочности по напряжениям сжатия удовлетворяет требованиям прочности.

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

Величина зазора в свободном состоянии
. В частном случае при  






7. Расчет прицепного шатуна
Производим расчет на устойчивость прицепного шатуна.

Определим осевой момент инерции сечения прицепного шатуна.
,
где  - момент инерции прямоугольника, наиболее удаленного от оси х-х (рисунок 11)


Рисунок 11 – Расчетная схема прицепного шатуна





Полная площадь поперечного сечения (площадь брутто) составляет

Минимальный радиус инерции





Гибкость , где n - коэффициент приведения длины, зависящий от способа заделки стержня. Рассматриваем прицепной шатун как стержень с шарнирно закрепленными концами, при этом способе n=0,7.
Тогда .
Так как явление продольного изгиба существует не только в пределах упругости, но и за ними, то для определения критического напряжения нельзя пользоваться формулой Эйлера, так как она дает завышенные значения.

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

Определим критическое напряжение стержня по формуле Ф.С. Ясинского: , значения эмпирических коэффициентов a и b определяем по таблицам [5]:,. Тогда



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


Коэффициент запаса устойчивости:





8. Расчет поршня
Статический расчет твердотельной модели поршня, выполненный в пакете Solid Works, производим в пакете Cosmos Works.

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

Рисунок 13 – Результат расчета на статическую прочность
По результатам расчета видно, что максимальные напряжения, а следовательно и минимальные коэффициенты запаса, расположились у основания бобышек (с внутренней стороны поршня) и составляют 136,7МПа, а так же на проточках под поршневые кольца (95МПа) и на периферии маслоотводных отверстий (100МПа).

Максимальные деформации характерны для диаметрально-противоположных точек “дна” поршня, расположенных в плоскости, перпендикулярной к оси поршневого пальца и составляют м.

Список использованной литературы

1        . И.П. Пелепейченко, Н.И. Кормилов “Тепловые двигатели”, - Харьков: ХАИ, 1977. – 108с.

2. И.П. Пелепейченко, В.И. Крирченко “Динамический расчет авиационного однорядного звездообразного двигателя на ЭВМ”, - Харьков: ХАИ, 1982. – 56с.

3. В.И. Крирченко “ Динамический расчет поршневого звездообразного двигателя”, - Харьков: ХАИ, 1973. – 68с.

4. Ю.А. Гусев, С.В. Епифанов, А.В. Белогуб “Поршни двигателей внутреннего сгорания”, - Харьков: ХАИ, 1999. – 32с.

5. Г.С. Писаренко, В.Г. Попков “Сопротивление материалов”, - Киев: Вища школа, 1986. – 776с.