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

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

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

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

“Авиационные двигатели”

Харьков 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. Давление в конце такта сжатия:

2. Температура в конце такта сжатия:

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

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

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

,

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

;

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

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

Действительное количество воздуха для сгорания 1кг топлива будет составлять:

.

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

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

.

Тогда:

.

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

Тогда

.

Подставим все известные величины в исходное уравнение:

Решим данное уравнение с помощью программного пакета MathCAD 13.

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

В результате получим: ,

2. Определим максимальное давление сгорания

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

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

1. Находим давление в конце такта расширения:

2.Находим температуру в конце такта расширения:

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

1. Индикаторное давление

,

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

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

Тогда:

2. Определяем индикаторный КПД.

3. Удельный индикаторный расход топлива равен

.

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

1. Среднее эффективное давление

,

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

Эффективный КПД нагнетателя:

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

Тогда

.

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

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

,

где

Среднее эффективное давление:

2. Механический КПД

3. Значение эффективного КПД

4. Удельный эффективный расход топлива

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

1. Рабочий объем цилиндра двигателя

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

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

.

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

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

5. Проверяем правильность расчетов основных размеров двигателя

.

 

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

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

Выполнение динамического расчета авиационного поршневого двигателя связано с довольно большим объемом расчетной работы, поэтому целесообразно проводить его на ЭВМ. Особенность такого расчета – учет в нем главного динамического эффекта, создаваемого прицепными механизмами, - сил второго порядка. Динамический расчет звездообразного двигателя без учета этих сил неприемлем, поскольку при этом создается ложное впечатление об уравновешенности механизма и о запасах прочности коленчатого вала, редуктора и воздушного винта.

2.1 Допущения

1. Учитываем только силы избыточного давления газов на поршень и силы инерции КШМ.

2. Индикаторные диаграммы во всех цилиндрах считаем одинаковыми. Теоретические диаграммы корректируем только в точке, соответствующей концу сгорания.

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

3. Предполагаем геометрическое подобие деталей КШМ проектируемого двигателя и прототипа.

4. Для расчета сил инерции реальное распределение масс в КШМ приводим к расчетной схеме, в которой все массы считаем точечными, сосредоточенными на осях поршневых пальцев и оси шатунной шейки коленчатого вала.

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

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

2.2 Определение основных размеров КШМ

Схема кривошипно-шатунного механизма с прицепными шатунами показана на рисунке 2.

Рисунок 2 Схема кривошипно-шатунного механизма с прицепными шатунами.

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

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

.

Угол прицепа:

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

.

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

,

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

.

Радиус остальных прицепов находим по формуле:

, где .

Результаты расчета сводим в таблицу 1

Таблица 1 Радиусы прицепов шатунов

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

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

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

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

,

.

Рисунок 3 – Схема приведения масс главного шатуна.

3. Приведенная масса поступательно-движущихся частей.

Эта масса различна в цилиндрах с главным шатуном и с прицепным.

В цилиндре с прицепным шатуном

,

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

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

В цилиндре с главным шатуном

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

,

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

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

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

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

 

Силы инерции поступательно-движущихся масс переменны по величине и направлению и действуют по осям цилиндров. Силу инерции в цилиндре с главным шатуном находят из уравнения:

,

а силу инерции в цилиндре с прицепным шатуном – из уравнения:

,

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

Силы инерции вращательно-движущихся масс находят по формулам:

,

.

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

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

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

Таким образом, расчету подлежат только давления в ходе расширения и сжатия, определяемые по формуле:

,

где  - давление в НМТ;

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

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

,

где

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

.

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

Таблица 1 – Давления и объемы в ходе расширения и сжатия

Сжатие			Расширение
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 Силы, действующие в центральном КШМ

В центральном КЩМ действуют силы, показанные на рисунке 4. Причем они имеют положительные значения. При направлении, противоположном указанному, силы считают отрицательными. Они равны:

,

,

,

,

.

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

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

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

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

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

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

Динамический расчет выполнен на ЭВМ. Это позволило значительно упростить процедуру расчета, а вместе с тем и повысить его точность. Результаты вычисления суммарных окружных и радиальных сил отображены в таблице 2; полных сил, действующих на шатунную шейку в окружном и в радиальном направлениях – в таблице 3.

Результаты расчета избыточного (Р) и абсолютного (PR) давления газов на поршень, силы инерции поступательно-движущихся частей (F), силы, действующей по оси цилиндра (PS), нормально к оси цилиндра (NS), окружной силы от одного цилиндра (Т) и радиальной силы от одного цилиндра (ZS) помещены в таблице 4.

Параметры ввода

Число цилиндров 7

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

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

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

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

Ход поршня, .

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

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

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

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

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

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

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

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

Таблица 2 – Суммарные окружные и радиальные силы.

-----------------------------------------------

| TC | ZC | ZCD |

-----------------------------------------------

| 24334.54000| -28539.68000| -28539.68000|

| 16694.21000| -34419.33000| -34419.33000|

| 10324.36000| -34938.66000| -28069.61000|

| 8007.59900| -25581.32000| -25581.32000|

| 9021.97900| -24514.77000| 27376.85000|

| 26920.21000| 28033.34000| 28033.34000|

| 33954.51000| 4737.10500| 4564.03100|

| 31262.81000| -15762.70000| -15762.70000|

-----------------------------------------------

Таблица 3 - Полные силы действующие на шатунную шейку

в окружном и в радиальном направлениях.

| Угол| TSI | ZSI | ZSID |

--------------------------------------------------------

|.00|.24335E+05 | -.20928E+05 | -.20928E+05 |

| 12.86|.15000E+05 | -.26998E+05 | -.26998E+05 |

| 25.71|.70217E+04 | -.28081E+05 | -.21212E+05 |

| 38.57|.32616E+04 | -.19630E+05 | -.19630E+05 |

| 51.43|.30707E+04 | -.19769E+05 |.32123E+05 |

| 64.29|.20062E+05 |.31336E+05 |.31336E+05 |

| 77.14|.26533E+05 |.64309E+04 |.62579E+04 |

| 90.00|.23651E+05 | -.15763E+05 | -.15763E+05 |

|102.86|.16913E+05 | -.30233E+05 | -.30233E+05 |

|115.71|.98361E+04 | -.37722E+05 | -.37722E+05 |

|128.57|.43731E+04 | -.39685E+05 | -.32816E+05 |

|141.43|.32616E+04 | -.31533E+05 | -.31533E+05 |

|154.29|.57193E+04 | -.31373E+05 |.20519E+05 |

|167.14|.25226E+05 |.20612E+05 |.20612E+05 |

|180.00|.33954E+05 | -.28748E+04 | -.30479E+04 |

|192.86|.32957E+05 | -.23184E+05 | -.23184E+05 |

|205.71|.27637E+05 | -.35398E+05 | -.35398E+05 |

|218.57|.21440E+05 | -.40371E+05 | -.40371E+05 |

|231.43|.16276E+05 | -.39685E+05 | -.32816E+05 |

|244.29|.14866E+05 | -.28884E+05 | -.28884E+05 |

|257.14|.16443E+05 | -.26209E+05 |.25683E+05 |

|270.00|.34532E+05 |.28033E+05 |.28033E+05 |

|282.86|.41376E+05 |.64309E+04 |.62578E+04 |

|295.71|.38121E+05 | -.12460E+05 | -.12460E+05 |

|308.57|.30286E+05 | -.23794E+05 | -.23794E+05 |

|321.43|.21440E+05 | -.28468E+05 | -.28468E+05 |

|334.29|.13627E+05 | -.28081E+05 | -.21212E+05 |

|347.14|.97015E+04 | -.18160E+05 | -.18160E+05 |

|360.00|.90220E+04 | -.16903E+05 |.34989E+05 |

|372.86|.25226E+05 |.35454E+05 |.35454E+05 |

|385.71|.30652E+05 |.11595E+05 |.11422E+05 |

|398.57|.26517E+05 | -.98114E+04 | -.98114E+04 |

|411.43|.18383E+05 | -.23794E+05 | -.23794E+05 |

|424.29|.98361E+04 | -.31117E+05 | -.31117E+05 |

|437.14|.29033E+04 | -.33245E+05 | -.26376E+05 |

|450.00|.39568E+03 | -.25581E+05 | -.25581E+05 |

|462.86|.16009E+04 | -.26209E+05 |.25683E+05 |

|475.71|.20062E+05 |.24731E+05 |.24731E+05 |

|488.57|.28003E+05 | -.87905E+01 | -.18186E+03 |

|501.43|.26517E+05 | -.21714E+05 | -.21714E+05 |

|514.29|.21032E+05 | -.35398E+05 | -.35398E+05 |

|527.14|.15000E+05 | -.41840E+05 | -.41840E+05 |

|540.00|.10324E+05 | -.42551E+05 | -.35682E+05 |

|552.86|.97013E+04 | -.33002E+05 | -.33002E+05 |

|565.71|.12325E+05 | -.31373E+05 |.20519E+05 |

|578.57|.31666E+05 |.22082E+05 |.22082E+05 |

|591.43|.39906E+05 | -.89038E+01 | -.18198E+03 |

|604.29|.38121E+05 | -.19065E+05 | -.19065E+05 |

|617.14|.31756E+05 | -.30234E+05 | -.30234E+05 |

|630.00|.24306E+05 | -.34419E+05 | -.34419E+05 |

|642.86|.17745E+05 | -.33245E+05 | -.26376E+05 |

|655.71|.14866E+05 | -.22279E+05 | -.22279E+05 |

|668.57|.14973E+05 | -.19769E+05 |.32123E+05 |

|681.43|.31666E+05 |.33985E+05 |.33985E+05 |

|694.29|.37257E+05 |.11595E+05 |.11422E+05 |

|707.14|.32957E+05 | -.83416E+04 | -.83416E+04 |

|720.00|.24335E+05 | -.20928E+05 | -.20928E+05 |

Таблица 4 – Результаты расчета давлений и сил.

Такт впуска

-------------------------------------------------------------------------------------

| Угол |Давлен. Р | Сила PR |Сила ин. F|Сум.силаPS| Сила NS |Окp.сила T|Pад.сила Z|

-------------------------------------------------------------------------------------

|.00|.0000E+00|.0000E+00|-.1295E+05|-.1295E+05|.0000E+00|.0000E+00|-.1295E+05|

| 12.86|.0000E+00|.0000E+00|-.1242E+05|-.1242E+05|-.7559E+03|-.3501E+04|-.1194E+05|

| 25.71|.0000E+00|.0000E+00|-.1090E+05|-.1090E+05|-.1300E+04|-.5900E+04|-.9254E+04|

| 38.57|.0000E+00|.0000E+00|-.8572E+04|-.8572E+04|-.1481E+04|-.6502E+04|-.5779E+04|

| 51.43|.0000E+00|.0000E+00|-.5725E+04|-.5725E+04|-.1251E+04|-.5256E+04|-.2592E+04|

| 64.29|.0000E+00|.0000E+00|-.2682E+04|-.2682E+04|-.6807E+03|-.2712E+04|-.5506E+03|

| 77.14|.0000E+00|.0000E+00|.2385E+03|.2385E+03|.6585E+02|.2472E+03|-.1113E+02|

| 90.00|.0000E+00|.0000E+00|.2777E+04|.2777E+04|.7882E+03|.2777E+04|-.7882E+03|

|102.86|.0000E+00|.0000E+00|.4766E+04|.4766E+04|.1316E+04|.4354E+04|-.2344E+04|

|115.71|.0000E+00|.0000E+00|.6146E+04|.6146E+04|.1560E+04|.4860E+04|-.4072E+04|

|128.57|.0000E+00|.0000E+00|.6961E+04|.6961E+04|.1521E+04|.4494E+04|-.5529E+04|

|141.43|.0000E+00|.0000E+00|.7336E+04|.7336E+04|.1267E+04|.3583E+04|-.6526E+04|

|154.29|.0000E+00|.0000E+00|.7434E+04|.7434E+04|.8868E+03|.2427E+04|-.7083E+04|

|167.14|.0000E+00|.0000E+00|.7416E+04|.7416E+04|.4514E+03|.1210E+04|-.7331E+04|

|180.00|.0000E+00|.0000E+00|.7396E+04|.7396E+04|.4638E-02|.1170E-01|-.7396E+04|

-------------------------------------------------------------------------------------

Такт сжатия

-------------------------------------------------------------------------------------

| Угол |Давлен. Р | Сила PR |Сила ин. F|Сум.силаPS| Сила NS |Окp.сила T|Pад.сила Z|

-------------------------------------------------------------------------------------

|180.00|.8927E+05|.1731E+03|.7396E+04|.7569E+04|.5239E-02|.1378E-01|-.7569E+04|

|192.86|.9021E+05|.1897E+03|.7416E+04|.7606E+04|-.4629E+03|-.1241E+04|-.7518E+04|

|205.71|.9315E+05|.2416E+03|.7434E+04|.7676E+04|-.9157E+03|-.2505E+04|-.7313E+04|

|218.57|.9843E+05|.3351E+03|.7336E+04|.7671E+04|-.1325E+04|-.3747E+04|-.6824E+04|

|231.43|.1067E+06|.4817E+03|.6961E+04|.7443E+04|-.1626E+04|-.4805E+04|-.5912E+04|

|244.29|.1192E+06|.7014E+03|.6146E+04|.6847E+04|-.1738E+04|-.5415E+04|-.4536E+04|

|257.14|.1375E+06|.1026E+04|.4766E+04|.5792E+04|-.1599E+04|-.5291E+04|-.2848E+04|

|270.00|.1649E+06|.1510E+04|.2777E+04|.4287E+04|-.1217E+04|-.4287E+04|-.1217E+04|

|282.86|.2063E+06|.2241E+04|.2386E+03|.2480E+04|-.6847E+03|-.2570E+04|-.1157E+03|

|295.71|.2702E+06|.3371E+04|-.2682E+04|.6883E+03|-.1747E+03|-.6960E+03|.1413E+03|

|308.57|.3705E+06|.5143E+04|-.5725E+04|-.5817E+03|.1271E+03|.5340E+03|-.2633E+03|

|321.43|.5265E+06|.7900E+04|-.8572E+04|-.6723E+03|.1161E+03|.5100E+03|-.4532E+03|

|334.29|.7503E+06|.1185E+05|-.1090E+05|.9563E+03|-.1141E+03|-.5177E+03|.8121E+03|

|347.14|.9975E+06|.1622E+05|-.1242E+05|.3803E+04|-.2314E+03|-.1072E+04|.3656E+04|

|360.00|.1117E+07|.1834E+05|-.1295E+05|.5387E+04|-.1026E-01|-.4786E-01|.5387E+04|

Такт сгорания

-------------------------------------------------------------------------------------

| Угол |Давлен. Р | Сила PR |Сила ин. F|Сум.силаPS| Сила NS |Окp.сила T|Pад.сила Z|

-------------------------------------------------------------------------------------

|360.00|.4054E+07|.7023E+05|-.1295E+05|.5728E+05|-.7929E-01|-.3724E+00|.5728E+05|

|372.86|.4298E+07|.7454E+05|-.1242E+05|.6212E+05|.3781E+04|.1751E+05|.5972E+05|

|385.71|.3308E+07|.5706E+05|-.1090E+05|.4616E+05|.5507E+04|.2499E+05|.3920E+05|

|398.57|.2390E+07|.4082E+05|-.8572E+04|.3225E+05|.5571E+04|.2446E+05|.2174E+05|

|411.43|.1731E+07|.2918E+05|-.5725E+04|.2345E+05|.5124E+04|.2153E+05|.1062E+05|

|424.29|.1295E+07|.2148E+05|-.2683E+04|.1880E+05|.4770E+04|.1900E+05|.3858E+04|

|437.14|.1011E+07|.1645E+05|.2384E+03|.1669E+05|.4609E+04|.1730E+05|-.7789E+03|

|450.00|.8227E+06|.1313E+05|.2777E+04|.1591E+05|.4515E+04|.1591E+05|-.4515E+04|

|462.86|.6964E+06|.1090E+05|.4766E+04|.1567E+05|.4326E+04|.1431E+05|-.7704E+04|

|475.71|.6104E+06|.9383E+04|.6146E+04|.1553E+05|.3941E+04|.1228E+05|-.1029E+05|

|488.57|.5517E+06|.8345E+04|.6961E+04|.1531E+05|.3344E+04|.9882E+04|-.1216E+05|

|501.43|.5122E+06|.7646E+04|.7336E+04|.1498E+05|.2588E+04|.7318E+04|-.1333E+05|

|514.29|.4868E+06|.7199E+04|.7434E+04|.1463E+05|.1746E+04|.4776E+04|-.1394E+05|

|527.14|.4727E+06|.6949E+04|.7416E+04|.1437E+05|.8743E+03|.2344E+04|-.1420E+05|

|540.00|.4682E+06|.6869E+04|.7396E+04|.1427E+05|.2962E-01|.8129E-01|-.1427E+05|

-------------------------------------------------------------------------------------

Такт выхлопа

-------------------------------------------------------------------------------------

| Угол |Давлен. Р | Сила PR |Сила ин. F|Сум.силаPS| Сила NS |Окp.сила T|Pад.сила Z|

-------------------------------------------------------------------------------------

|540.00|.0000E+00|.0000E+00|.7396E+04|.7396E+04|.1536E-01|.4214E-01|-.7396E+04|

|552.86|.0000E+00|.0000E+00|.7416E+04|.7416E+04|-.4513E+03|-.1210E+04|-.7331E+04|

|565.71|.0000E+00|.0000E+00|.7434E+04|.7434E+04|-.8868E+03|-.2427E+04|-.7083E+04|

|578.57|.0000E+00|.0000E+00|.7336E+04|.7336E+04|-.1267E+04|-.3583E+04|-.6526E+04|

|591.43|.0000E+00|.0000E+00|.6961E+04|.6961E+04|-.1521E+04|-.4494E+04|-.5529E+04|

|604.29|.0000E+00|.0000E+00|.6146E+04|.6146E+04|-.1560E+04|-.4860E+04|-.4072E+04|

|617.14|.0000E+00|.0000E+00|.4766E+04|.4766E+04|-.1316E+04|-.4354E+04|-.2344E+04|

|630.00|.0000E+00|.0000E+00|.2777E+04|.2777E+04|-.7882E+03|-.2777E+04|-.7882E+03|

|642.86|.0000E+00|.0000E+00|.2386E+03|.2386E+03|-.6588E+02|-.2473E+03|-.1114E+02|

|655.71|.0000E+00|.0000E+00|-.2682E+04|-.2682E+04|.6807E+03|.2712E+04|-.5505E+03|

|668.57|.0000E+00|.0000E+00|-.5725E+04|-.5725E+04|.1251E+04|.5256E+04|-.2592E+04|

|681.43|.0000E+00|.0000E+00|-.8572E+04|-.8572E+04|.1481E+04|.6502E+04|-.5779E+04|

|694.29|.0000E+00|.0000E+00|-.1090E+05|-.1090E+05|.1300E+04|.5900E+04|-.9254E+04|

|707.14|.0000E+00|.0000E+00|-.1242E+05|-.1242E+05|.7560E+03|.3501E+04|-.1194E+05|

|720.00|.0000E+00|.0000E+00|-.1295E+05|-.1295E+05|.2237E-01|.1066E+00|-.1295E+05|

-------------------------------------------------------------------------------------

 

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

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

,

где:

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

- сила инерции от неуравновешенных частей равна:

Тогда получим, что

Рассмотрим вопрос уравновешивания сил инерции поступательно движущихся масс.

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

,

где  - поступательно движущаяся масса, относящаяся к одному цилиндру, ;

Z – число цилиндров в одной звезде.

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

Такую силу легко уравновесить, добавив к противовесам соответствующую массу.

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

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

,

где  - разность поступательно движущихся масс главного и бокового цилиндра: ,

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

Рисунок 5 – Результирующий вектор сил инерции первого порядка.

 

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

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

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

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

При работе двигателя колено вала нагружается следующими силами (рисунок 6)

Рисунок 6 – Силовое нагружение колена

1) В плоскости колена действует сила Z

2) Перпендикулярно к плоскости колена действует сила Т

3) В плоскости колена действует сила инерции от вращающихся масс шатуна:

; ;

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

где

Тогда

Рисунок 7 – Разнос масс шатуна

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

(Все необходимые объемы соответствующих частей деталей и радиусы их центров масс определены по твердотельной модели исполненной в пакете Solid Works)

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

6) В плоскости колена действуют силы инерции противовесов Р_ПР.

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

.

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

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

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

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

Рисунок 8 – Расчетная модель шатунной шейки

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

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

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

;

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

Рисунок 7 – Твердотельная модель коленчатого вала (фрагмент)

Опасные точки шатунной шейки расположены у масляного отверстия (смотри рисунок 7.). Если его ось составляет с плоскостью колена угол g (g=300º, так как при этом положении силы Т и Z минимальны), то изгибающий момент в плоскости, проходящей через ось масляного отверстия, равен:

Знак “-” показывает, что момент  вызывает у края отверстия напряжения сжатия.

Таблица 5 – Результаты расчета

α	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 Определение запаса прочности в коренной шейке

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

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

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

Таблица 6 – Величина крутящего момента и касательных напряжений

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с.

Министерство образования и науки Украины Национальный аэрокосмический университет им. Н.Е. Жуковского “ХАИ”Расчетно-пояснительная записка к курсовому проекту по дисциплине “Авиационные двигатели”

 

• 14:21 05 Июля 2016
  Российско-китайский медуниверситет откроет магистерскую программу по TCM
• 01:21 05 Июля 2016
  На форуме МГУ Россия и Китай обсудят научно-образовательное сотрудничество
• 11:00 04 Июля 2016
  Реморенко: "ЕГЭ-туризм" в России сходит на нет
• 10:52 04 Июля 2016
  Летняя российско-австрийская школа откроется в Москве
• 15:31 01 Июля 2016
  В СПбГУ могут воссоздать географический факультет

Заказать уникальную работу