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

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

ФГАОУ ВПО ЮУрГУ (НИУ)
Политехнический институт

Аэрокосмический факультет
Кафедра «Летательных Аппаратов»

СЕМЕСТРОВАЯ РАБОТА

по дисциплине «Прочность конструкции» на тему

Расчет на прочность крыла ЛА

Руководитель
Овчинников А.М.
____________________ «___»____________2017 г. Автор работы
студент группы П-424
Иванов С.В.
____________________
« » 2017 г.

Работа защищена с оценкой
____________________ « » 2017 г.

Челябинск, 2017

Аннотация

Иванов М.В. Проектирование силового набора крыла самолета: семестровая работа по дисциплине «Прочность конструкций» - Челябинск: ЮУрГУ, 2017 - 25 с., 19 илл., 2 наименования литературы.

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

Проведен проверочный расчет спроектированного крыла в программном пакете Ansys.

Исходные данные. 2

Проектировочный расчет.. 3

1.1 Описание нагрузок. 3

1.2 Расчетная схема конструкции крыла. 7

Подбор сечения лонжеронного крыла.. 8

2.1 Подбор обшивки. 8

2.2 Подбор стрингеров нижней панели. 9

2.3 Расчет силовых элементов крыла на устойчивость. 10

2.4 Подбор поясов лонжеронов верхней панели. 12

2.5 Проверка нижней панели крыла на сжатие. 13

2.6 Подбор толщины стенок лонжеронов. 14

Проверочный расчет.. 16

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

В данной работе предлагается провести проектировочный расчет силового подкрепления крыла самолета, а затем выполнить проверочный расчет силового набора с помощью конечно-элементного пакета Ansys.

Принимаются следующие исходные данные для расчета:

1) длина крыла

2) хорда корневой кромки

3) хорда концевой кромки крыла

4) Масса самолета

5) Масса двигателей

7) Координаты подвески двигателей от конца крыла:

8) самолет движется на крейсерской;

9) материал обшивки, поясов лонжеронов, стенок лонжеронов, стрингеров - алюминиевый сплав АМг6: предел прочности модуль упругости

10) Аэродинамический профиль ЦАГИ-734.

Рисунок 1. Профиль крыла ЦАГИ-734.

1. Проектировочный расчет

1.1 Описание нагрузок

На крыло в полете действует распределенная подъемная сила распределенный вес крыла m и сосредоточенные массовые силы агрегатов - веса двигателей

Крыло длиной 8 [м] разбивается на 30 участков длиной [м] каждый. Разбиение показано на рисунке 2.

Подъемная сила на участках крыла и перерезывающая сила определятся по формулам:

Площадь i -го участка крыла; - коэффициент подъемной силы, для выбранного профиля =0,528; - плотность воздуха

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

Интегрирование выполним так же, как и при расчете перерезывающей силы, численным методом трапеций. Для участка крыла Δξi определяем приращение изгибающего момента:

Суммируя с нарастающим итогом приращения ΔMi от края крыла, получаем изгибающий момент в сечении:

Крутящий момент определим по формуле:

В таблице 1 приведены расчетные значения.

Таблица 1.

По данным таблицы 1 построим графики изменения перерезывающей силы и моментов.

Рисунок 2. Изменение подъемной силы по длине крыла.

Рисунок 3. Изменение перерезывающей силы по длине крыла.

Рисунок 4. Изменение изгибающего момента по длине крыла

Рисунок 5. Изменение крутящего момента по длине крыла

1.2 Расчетная схема конструкции крыла

При назначении силового набора крыла следует руководствоваться следующими рекомендациями:

1) передний лонжерон располагается на расстоянии от носка сечения, а задний - на, где - хорда сечения крыла;

2) расстояние между соседними стрингерами лежит в пределах от 120...300 мм для лонжеронного крыла;

3) расстояние между нервюрами в лонжеронном крыле обычно принимают 200...300 мм.

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

Назначение силового набора приведено на рисунке 7.

Рисунок 6. Назначение силового набора крыла.

2. Подбор сечения лонжеронного крыла

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

Здесь Н - плечо пары нормальных сил

где μ = 0,95 - коэффициент, показывающий насколько расстояние между центрами тяжести поясов лонжерона меньше габаритной высоты лонжерона; Н1 и Н2 - габаритные высоты лонжеронов. Под Н1 - понимается высота самого высокого лонжерона в сечении крыла.

2.1 Подбор обшивки

Минимальную потребную толщину обшивки рассчитываем из условия работы ее на сдвиг при кручении крыла по формуле

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

1.4.2 Подбор поясов лонжеронов нижней панели.

Минимальную потребную площадь поперечного сечения первого лонжерона находим по формуле

где к = 0,7...0,8 - коэффициент, определяющий долю нормальной силы N, воспринимаемой поясами лонжеронов; - разрушающее напряжение материала растянутого пояса.

Для второго лонжерона принимаем:

По потребным площадям и подбираем ближайшие большие по пло- щади стандартные прессованные профили, . Выбираем профили ПР 101 и ПР 111 - уголкового сечения, не равнополочные (ГОСТ 13738 - 91);

Рисунок 7. Профиль ПР 101.

Для первого лонжерона выбран профиль ПР101-47.

2.2 Подбор стрингеров нижней панели.

Задаёмся количеством стрингеров m, исходя из диапазона рекомендуемых расстояний между ними. Стрингеры в пределах межлонжеронной части крыла располагаем равномерно и находим фактическое расстояние между ними

где В - ширина межлонжеронной части крыла; m - число стрингеров в верхней (нижней) панели крыла.

Вычисляем нормальную силу в поясах лонжеронов

и в обшивке

где - редукционный коэффициент.

Оставшаяся часть растягивающей силы воспринимается стрингерами. Минимальную потребную площадь стрингера вычисляем по формуле

В формулах - разрушающие напряжения при растяжении пояса лонжерона, обшивки и стрингера соответственно.

По потребной величине подбираем стандартный ближайший по площади профиль. Выбираем профиль ПР 100- уголкового сечения, равнополочные (ГОСТ 13737-90);

Рисунок 8. Профиль ПР 100 (ГОСТ 13737-90).

Необходимое условие выполняет профиль ПР100-53.

2.3 Расчет силовых элементов крыла на устойчивость.

Работа обшивки на устойчивость зависит от работы её отдельных участков. Участок обшивки шириной и длиною а (а - расстояние между нервюрами) рассматривается как плоская пластинка, которая опирается по всему контуру на стрингеры и нервюры (рис. Д.1).

Рисунок 9. Фрагмент панели крыла.

Критическое напряжение пластинки при сжатии в направлении стрингерного набора определяют по формуле

где к - коэффициент, учитывающий характер закрепления пластинки по контру. При а ≥ коэффициент к= 4.

Стрингер

Расчет на местную потерю устойчивости

Критическое напряжение местной потери устойчивости для i-той полки стрингера (рис. Д1), рассматриваемой как пластинка шириной bi и толщиной δi, определяется по формуле:

где к= 0,46 - коэффициент для полок стрингера, имеющих один свободный край вдоль длинной стороны;

Введем поправку на пластичность материала:

Расчет на общую потерю устойчивости

Критические напряжения общей потери устойчивости стрингера определяются по формуле

Здесь m - коэффициент, зависящий от характера закрепления стрингера по концам (принято в крыле закрепление стрингера по концам принимать в виде так называемой приторцовки, для которой m = 2); F, Ix - площадь и момент инерции поперечного сечения стрингера относительно оси х, проходящей через центр тяжести стрингера и параллельной обшивке (в приближенном проектировочном расчете); а - расстояние между нервюрами.

Поправка на пластичность материала

Критическое напряжение потери устойчивости стрингера равно минимальному из двух напряжений

2.4 Подбор поясов лонжеронов верхней панели

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

Определяем эффективную площадь стрингера и присоединенной к нему обшивки

Требуемые площади сечений поясов лонжеронов рассчитываем по формулам

Здесь σкр - критическое напряжение местной потери устойчивости пояса самого высокого лонжерона. Этой величиной следует сначала задаться в пределах:

По вычисленным площадям подбираем стандартные профили с

По потребным площадям и подбираем ближайшие большие по пло- щади стандартные прессованные профили. Выбираем профили ПР 101 и ПР 111 - уголкового сечения, не равнополочные (ГОСТ 13738 - 91);

Рисунок 10. Профиль ПР 101.

Для первого лонжерона выбран профиль ПР111-40.

2.5 Проверка нижней панели крыла на сжатие

Критические напряжения потери устойчивости поясов первого и второго лонжеронов нижней панели, определятся по формулам

Нижняя панель крыла, подобранная из условия работы её на растяжение в расчетном случае А, будет работать на сжатие в расчетном случае D. Поэтому ее следует проверить на устойчивость по случаю D:

Осевая сила в панели в расчетном случае D.

2.6 Подбор толщины стенок лонжеронов.

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

Где и - расчетные разрушающие значения силовых факторов для случая А; - часть перерезывающей силы, воспринимаемой стенками лонжеронов; - перерезывающая сила, воспринимаемая стенкой первого лонжерона; - перерезывающая сила, воспринимаемая стенкой второго лонжерона; Н= 0,5(Н1 + Н2) - средняя высота лонжеронов в расчетном сечении; - угол сходимости поясов лонжеронов (в радианах)

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

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

3. Проверочный расчет

На основании проведенного проектировочного расчета была построена 3D-модель конструкции крыла с силовым набором (рисунок 11).

Рисунок 11. 3D-модель конструкции крыла с силовым набором.

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

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

Силовой набор и обшивка принимается оболочечными элементами Shell 181, каждой поверхности присваивается соответствующая толщина.

По координатам, указанным ранее, созданы элементы сосредоточенный массы (элемент Mass 21). Эти элементы соединены жестко (Rigid Region) с узлами, соответствующими нижним поясам лонжеронов. Эти элементы соответствуют сосредоточенной силе от агрегатов (двигателей).

Крыло считается закрепленным абсолютно жестко во всех направлениях (All DOF) по корневому торцу.

На рисунке 12 приведена конечно-элементная модель с сосредоточенными силами и закрепленной стороной.

Рисунок 12. Конечно-элементная модель для расчета.

На рисунках показан результат расчета напряжений (Nodal solution).

Рисунок 13. Распределение главных растягивающих напряжений.

Рисунок 14. Распределение главных сжимающих напряжений.

Для сравнения приведем расчеты (Element solution)

Рисунок 15. Распределение главных растягивающих напряжений.

Рисунок 16. Распределение главных сжимающих напряжений.

Рисунок 17. Распределение эквивалентных напряжений.

Далее проведен расчет потери устойчивости (Eigen Buckling) с учетом рассчитанных эффектов предварительного нагружения (Pre-Stress Effects). В этом расчете были вычислены первые 5 форм потерь устойчивости конструкции.

Все вычисленные формы потери устойчивости локализованы в растянутой зоне хвостовой части крыла, и отличаются друг от друга количеством возникающих волн. Первая форма потери устойчивости приведена на рисунке 18, пятая - на рисунке 19.

Рисунок 18. Первая форма потери устойчивости.

Рисунок 19. Пятая форма потери устойчивости.

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

Геометрические характеристики силового набора крыла и расчетные напряжения.

Толщина обшивки: ;

Стрингеры: Профиль ПР 100- уголкового сечения, равнополочные (ГОСТ 13737-90);

Рисунок 20. Профиль ПР 100 (ГОСТ 13737-90).

Профиль ПР100-53.

Для второго лонжерона выбран профиль ПР111-38.

Для второго лонжерона выбран профиль ПР101-47.

Численные результаты проверочного расчета:

Проверочные расчеты показали, что спроектированная конструкция неработоспособна, так как:

1) в силовом наборе возникают напряжения, большие предела прочности выбранного материала:

2) происходит потеря устойчивости обшивки (см. рисунки 18, 19).

На основании проверочного расчета сформулированы следующие рекомендации по изменению конструкции:

1) необходимо увеличить площади несущих элементов силового набора, выбрав при этом угловые профили с большей толщиной стенок и меньшей длиной.

2) Увеличить толщину стенок лонжеронов.

3) в проверочных расчетах необходимо учитывать подкрепление хвостовой части (выполняется в виде сотового наполнителя, а также силовых элементов механизации крыла);

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

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

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

1) Тарасов, Ю.Л., Лавров, Б.А. Расчет на прочность элементов конструкции самолета [Текст] / Ю.Л. Тарасов, Б.А. Лавров - Самара, Самарский государственный аэрокосмический университет, 2000 г. - 112 с.

2) Мехеда, В.А. Подбор сечений силовых элементов нестреловидных крыльев [Текст] / В. А. Мехеда - Самара, Самарский государственный аэрокосмический университет, 2008 г. - 48 с.

Скачать: У вас нет доступа к скачиванию файлов с нашего сервера.

Расчет аэродинамических характеристик крыла с использованием программного комплекса ANSYS CFX

Создание летательного аппарата нового поколения невозможно без анализа его аэродинамических характеристик еще на ранних стадиях проектирования. От глубины исследования формы несущих поверхностей и обводов планера напрямую зависят летно-технические характеристики разрабатываемого самолета. Развитие теоретических основ численных методик расчета аэродинамических характеристик летательных аппаратов можно разделить на несколько этапов:

линейная теория (60-е годы);
нелинейная теория полного потенциала скорости (70-е годы);
уравнения Эйлера (80-е годы);
уравнения Навье — Стокса, осредненные по Рейнольдсу (90-е годы).

Физику процесса обтекания тела произвольной формы потоком газа наилучшим образом отражают методики, основанные на решениях уравнений Навье — Стокса. С появлением программных средств, базирующихся на численных решениях уравнений Навье — Стокса, стало возможно получить расчетным путем ряд важных аэродинамических характеристик самолета, в частности вычислить максимальное значение коэффициента подъемной силы Cy max . При расчетах аэродинамических характеристик объектов сложной пространственной конфигурации с использованием такого подхода требуются большие объемы оперативной памяти компьютера, поскольку допустимые размеры расчетной сетки пропорциональны объему оперативной памяти компьютера. Рост возможностей вычислительной техники, наблюдаемый в последние годы, позволяет применять программы, основанные на численных решениях уравнений Навье — Стокса, для расчета характеристик обтекания таких объектов, как самолет. Одной из популярных коммерческих программ в этой области является ANSYS CFX (лицензия ЦАГИ № 501024).

Использование CFX в области авиастроения является рациональным, поскольку пакет ANSYS, помимо аэродинамического модуля CFX, содержит ряд других вычислительных модулей (STRUCTURAL, FATIQUE и д.р.), что обеспечивает возможность совместного решения задач аэродинамики, аэроупругости и прочности.

Рассмотрим особенности расчета обтекания прямого крыла бесконечного размаха с профилем GA(W)-1. Этот профиль был создан известным американским аэродинамиком Уиткомбом для применения на дозвуковых скоростях полета.

Комплекс ANSYS оснащен встроенными интерфейсами ряда основных CAD-программ. Геометрическая модель, созданная в программе трехмерного графического моделирования, считывается любой из программ комплекса. Твердотельная геометрическая модель отсека крыла, сохраненная в формате Parasolid, была импортирована в профессиональный сеточный генератор ANSYS ICEM, где методом Octree была построена неструктурированная расчетная сетка, состоящая из 3 млн объемных тетраэдрических элементов (рис. 1). Вблизи поверхности крыла параметры Tetra Size Ratio и Height Ratio были равны 1.2. Максимальный размер элементов на передней кромке крыла составил 1 мм. Для обеспечения нужной точности решения и сходимости расчета элементы расчетной сетки имели Aspect Ratio более 0.3 и Min Angle более 20°. Кроме того, необходимо, чтобы габаритные размеры расчетной области многократно превышали характерный размер исследуемого объекта. В данном случае использовалась прямоугольная расчетная область длиной 35 и высотой 30 м. Размах крыла равен 4 м, а хорда крыла — 3,3 м. Моделирование крыла бесконечного размаха осуществлялось путем задания в препроцессоре CFX-PRE справа и слева от крыла граничных условий типа Symmetry. Типы граничных условий, используемых в данной задаче, показаны на рис. 2.

В пристеночных областях при построении расчетной сетки для наилучшего моделирования пограничного слоя образованы слои призматических элементов (см. рис. 1). При решении задачи обтекания крыла (где одной из расчетных величин является касательное напряжение) очень важно контролировать величину Y+ . Значение Y+ характеризует относительную высоту первой ячейки пограничного слоя, которая задается в ICEM при построении призматических элементов. После окончания вычислений в среде постпроцессора CFX-POST можно визуализировать Y+ на расчетной модели (рис. 3).

При использовании методик, основанных на численных решениях уравнений Навье — Стокса, качество полученного результата во многом зависит от выбора модели турбулентности. В программном комплексе ANSYS CFX реализовано достаточно большое число моделей турбулентности. Однако ни одна из них не является универсальной для всех существующих классов задач. Из многообразия моделей турбулентности, используемых при расчетах аэродинамических характеристик, можно выделить известные модели турбулентности k -ε и k -ω. Они являются двупараметрическими моделями турбулентности, которые базируются на рассмотрении кинетической энергии турбулентных пульсаций k . В качестве второго уравнения применяют уравнение либо переноса скорости диссипации турбулентной энергии ε, либо удельной скорости диссипации энергии ω. Модель переноса касательных напряжений SST (двухслойная модель Ментера) использует модель k -ω в пристеночной области и преобразованную модель k -ε вдали от стенки. В новые версии программы CFX включен бета-вариант модели турбулентности Spalart-Allmaras (S-A). Эта модель является однопараметрической, использующей одно дифференциальное уравнение переноса.

Расчеты с применением программного комплекса ANSYS CFX проводились на сервере с 8-ядерным процессором Intel Xeon 2,83 ГГц и 16 Гбайт ОЗУ. Для получения стационарного решения в зависимости от типа модели турбулентности и угла атаки крыла потребовалось осуществить 40-60 итераций.

Вычисления проводились при числе Маха 0,2 и числе Рейнольдса 2,2Ѕ106. В препроцессоре ANSYS CFX отсутствует возможность напрямую задавать число Рейнольдса. В связи с этим число Рейнольдса вычислялось в CFX-PRE по величине статического давления, соответствующего определенному коэффициенту кинематической вязкости.

В результате проведенных расчетов были получены величины сил и моментов, действующих на отсек крыла на заданных углах атаки. Зависимость коэффициента подъемной силы Сy от угла атаки сравнивалась с аналогичными экспериментальными данными, полученными американскими специалистами NASA Венцем и Ситхарамом (SAE Paper 740365). На линейном участке все рассмотренные модели турбулентности продемонстрировали удовлетворительное совпадение расчетных и экспериментальных данных. В зоне Сy max максимальное соответствие с экспериментальными данными показала модель турбулентности SST (рис. 4). С использованием постпроцессора CFX-POST файл с результатами расчета позволяет визуализировать картину обтекания крыла. Линии тока и поле скоростей хорошо иллюстрируют отрывное течение, соответствующее углу атаки, при котором достигается Cy max крыла (рис. 5).

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

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

Аэродинамическая нагрузка распределяется по размаху крыла по закону, близкому к параболическому. Для упрощения заменим его трапециевидным законом (Рис. 2.2). Если принять допущение, что С y постоянен по размаху крыла, то закон изменения аэродинамической силы q az пропорционален хорде крыла b z:

где Y - подъемная сила создаваемая крылом;

S k - несущая площадь полукрыльев, равная S k = S - b 0d ф = 61;

d ф - диаметр фюзеляжа;

b 0 - хорда корневой нервюры;

b z - значение текущей хорды.

Значение текущей хорды крыла bz вычислим из предлагаемой формулы:

где b к - хорда концевой нервюры;

Длина полукрыла без центроплана, равная;

Подставив в (3.10) уравнение (3.11), получим:

Считаем, что топливо распределено по крылу равномерно, тогда распределенная нагрузка от массовых сил крыла (его собственного веса и топлива) изменяется по его размаху тоже пропорционально хорде b z:

где m k - масса конструкции полукрыльев, равная m k = m k m взл = 1890;

m Т - масса топлива, равная m Т = 0,85m Tmax = 3570 ;

g - ускорение свободного падения, равная g = 9,81.

Рис.

Произведем расчет распределенных аэродинамических q az и массовых нагрузок q крz в концевой, корневой части крыла и (к примеру) в районе элеронов:

1) Расчет распределенной нагрузки на конце крыла, т.е. при Z= 0:

2) Расчет распределенной нагрузки в корневом сечении, т.е. при Z== 13,23:

3) Расчет распределенной нагрузки в районе двигатели+шасси, т.е. при Z=l 1 =1,17

5665,94-2142,07=3523,87Н/м

Рис. 2.3. Схема возникновения крутящего момента в сечении крыла

Поэтому погонный крутящий момент от распределенных аэродинамических q az и массовых сил крыла q крz равен:

Нм/м (3.15)

Приводим подобные, и получим:

Нм/м (3.16)

Обычно топливо в крыле расположено в передней части крыла, поэтому ц.м. топлива совпадает с ц.м. крыла. С учетом этого предположения формула (3.15) будет иметь вид:

Нм/м (3.17)

Подставим известные величины в формулу (3.17), получим:

Нм/м (3.18)

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

1) Расчет крутящего момента на конце крыла, т.е. при Z= 0:

2) Расчет крутящего момента в корневой части крыла, т.е. при Z= 13,23:

3) Расчет крутящего момента в районе двигателя+шасси, т.е. при Z= 1,17:

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

Рис.

Из рисунка видно, что равен (знак «минус» означает, что момент направлен в противоположную сторону, против часовой стрелки):

(Нм ), (3.19)

где - расстояние от ц.м. двигателя до ц.ж. крыла.

Так как двигатели находятся на разном расстоянии от ц.ж. крыла, то и моменты они будут создавать разные. По известным данным найдем:

Базовым вариантом является региональный самолет Ан-148-100, обеспечивающий перевозку в одноклассной компоновке от 70 пассажиров с шагом кресел 864 мм (34‘’) до 80 пассажиров с шагом кресел 762 мм (30‘’). С целью обеспечения гибкости удовлетворения требований различных авиакомпаний, а также с целью снижения эксплуатационных затрат и повышения рентабельности перевозок предусматривается сертификация базового самолета в вариантах с максимальной дальностью полета от 2200 до 5100 км. Крейсерская скорость полета 820-870 км/ч. Проведенные маркетинговые исследования показали, что базовый самолет по своим технико-экономическим характеристикам отвечает требованиям большого количества авиакомпаний.

Самолет Ан-148-100 выполнен по схеме высокоплана с двигателями Д-436-148, размещенными на пилонах под крылом. Это позволяет повысить уровень защищенности двигателей и конструкции крыла от повреждений посторонними предметами. Наличие вспомогательной силовой установки, бортовой системы регистрации состояния самолета, а также высокий уровень эксплуатабельности и надежности систем позволяют использовать Ан-148-100 на сети технически слабооснащенных аэродромов.

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

Комфорт пассажирам обеспечивается на уровне комфорта на магистральных самолетах и достигнут рациональной компоновкой и составом сервисных помещений, глубокой эргономической оптимизацией общего и индивидуального пространства пассажирского салона, применением современных кресел, дизайна и материалов интерьера, а также созданием комфортных климатических условий и низкого уровня шума. Рационально выбранная длина пассажирского салона и размещение пассажиров в ряду по схеме 2+3 позволяют силами эксплуатанта получить различные одноклассные и смешанные компоновки в диапазоне 55-80 пассажиров с салонами экономического, бизнес и первого класса. Высокая степень преемственности конструктивно-технологических решений и эксплуатационной унификации Ан-148-100 с успешно эксплуатируемыми самолетами «Ан», использованием «Hi-Tech» компонентов оборудования и систем отечественного и зарубежного производств обеспечивают самолету Ан-148-100 высокий конкурентный уровень экономической эффективности, технического и эксплуатационного совершенства.

Техническое обслуживание самолета Ан-148-100 основано на удовлетворении требований международных стандартов (ICAO, MSG-3) и обеспечивает поддержание летной годности самолета в пределах жизненного цикла эксплуатации по состоянию с интенсивностью до 300 ч в месяц с коэффициентом готовности более 99,4%, при минимизации затрат на ТО (1,3 чел-ч на 1 час налета).

Семейство самолетов Ан-148 также включает следующие модификации:

пассажирский самолет, обеспечивающий перевозку 40-55 пассажиров на дальность до 7000 км; административный на 10 – 30 пасс. с дальностью до 8700 км;

грузовой вариант с боковой грузовой дверью для перевозок генеральных грузов на поддонах и в контейнерах;

грузо-пассажирский вариант для смешанных перевозок «пассажиры + груз».

Принципиальной особенностью создания семейства Ан-148 является использование максимальной унификации и преемственности агрегатов и компонентов базового самолета – крыла, оперения, фюзеляжа, силовой установки, пассажирского и самолетного оборудования.

Расчет крыла большого удлинения

Геометрические данные крыла

–площадь стреловидного крыла;

Удлинение стреловидного крыла;

Размах стреловидного крыла;

Сужение стреловидного крыла;

Корневая хорда крыла;

Концевая хорда крыла;

Угол стреловидности крыла по передней кромке.

Так как крыло данного самолета стреловидное и угол по передней кромке более 15° (рис. 1), вводим эквивалентное равновеликое по площади прямое крыло, и все расчеты проводим для этого эквивалентного крыла. Прямое крыло введем путем поворота стреловидного так, чтобы прямая проходящая по половине хорды прямого крыла была перпендикулярна оси фюзеляжа (рис. 2). При этом размах спрямленного крыла

Площадь спрямленного крыла:

причем в качестве параметра примем значение, равное расстоянию от конца консоли спрямленного крыла до оси самолёта, так как схема данного самолета – высокоплан (рис. 3)

. Тогда .

Найдем относительную координату линии центров давления. Для этого определим коэффициент подъемной силы для расчетного случая А.

Взлетный вес данного самолета;

- плотность воздуха на высоте Н = 0 км;

- крейсерская скорость самолета ( = кг),

Скорость пикирования,

Тогда: С х = 0,013; С д = 0,339; α 0 = 2 о

Лонжероны в крыле располагаем:

Передний лонжерон на расстоянии 15% хорды от носка крыла;

Задний лонжерон на расстоянии 75% хорды от носка крыла (рис. 5).

В расчетном сечении () высота переднего лонжерона , заднего- .

Определение нагрузок на крыло

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

Массы агрегатов находим через их относительные массы от взлетной массы самолета:

Масса крыла;

Масса силовой установки;

Так как на самолёте 2 двигателя, то массу одного двигателя примем равной

Распределение воздушной нагрузки по длине крыла.

По длине крыла нагрузка распределяется по закону относительной циркуляции:

где - относительная циркуляция,

В случае стреловидного крыла относительная циркуляция определяется по формуле:

, где - влияние стреловидности крыла, ( - угол стреловидности по четверти хорды).

Таблица – Распределение воздушной нагрузки по консоли крыла

zотн	0	0,1	0,2	0,3	0,4	0,5	0,6	0,7	0,8	0,9	1
 Г45	-0,235	-0,175	-0,123	-0,072	-0,025	0,025	0,073	0,111	0,135	0,14	0
Г пл	1,3859	1,3701	1,3245	1,2524	1,1601	1,0543	0,9419	0,8271	0,7051	0,5434	0
Г	1,27404	1,2868	1,265952	1,218128	1,1482	1,0662	0,976648	0,879936	0,76936	0,61004	0
qв,H/м	36430,7	36795,5	36199,4	34831,9	32832,3	30487,6	27926,9	25161,4	21999,5	17443,9	0,0

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

, где - хорда крыла.

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

, где - удельный вес топлива.

где - вес топлива (для самолёта АН 148 ).

Суммарная погонная нагрузка на крыло находится по формуле:

Начало координат поместим в корне крыла, сечения нумеруем от корня в направлении конца крыла, начиная с .

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

	z, м	b(z), м				, кг/м	, кг/м	, кг/м	, кг/м
0	0	4,93	1,3435	-0,060421	1,283079	4048,02	505,33	2187,441	1355,25
0,1	1,462	4,559	1,3298	-0,044994	1,284806	4053,46	467,30	1870,603	1715,56
0,2	2,924	4,188	1,2908	-0,031625	1,259175	3972,60	429,27	1578,541	1964,79
0,2	2,924	4,188	1,2908	-0,031625	1,259175	3972,60	429,27	0	3543,33
0,3	4,386	3,817	1,2228	-0,018512	1,204288	3799,44	391,24	0	3408,20
0,4	5,848	3,446	1,1484	1,141972	3602,84	353,22	0	3249,62
0,4	5,848	3,446	1,1484	1,141972	3602,84	353,22	1068,742	2180,88
0,5	7,31	3,075	1,057	0,006428	1,063428	3355,03	315,19	851,0063	2188,84
0,6	8,772	2,704	0,9571	0,018769	0,975869	3078,79	277,16	658,0454	2143,59
0,7	10,234	2,333	0,8538	0,028539	0,882339	2783,71	239,13	489,86	2054,72
0,8	11,696	1,962	0,743	0,03471	0,77771	2453,62	201,11	346,45	1906,06
0,9	13,158	1,591	0,6091	0,035996	0,645096	2035,23	163,08	227,8153	1644,34
0,95	13,889	1,4055	0,4593	0,032139	0,491439	1550,45	144,06	177,7887	1228,60
1	14,62	1,22	0	0	0	0,00	0,00	0	0

Строим эпюры функций , и (рис. 7)

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

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

, ,

Расчет производим по следующим формулам:

;

; ,

, .

Аналогично рассчитываем величины изгибающих моментов:

Полученные результаты заносим в таблицу 2.

Таблица 2

	z,м	ΔQ, кг	Q, кг	ΔM, кгм	M, кгм
0	0	2244,77	20592,41	196758,3	1016728
0,1	1,462	2690,34	18347,64	172115,8	819969,8
0,2	2,924	2969,13	15657,30	152033,9	647854
0,3	4,386	3127,09	12688,17	130883,4	495820,1
0,4	5,848	3194,27	53414,20	121865,8	364936,7
0,5	7,31	3167,01	43712,46	87477,02	243070,9
0,6	8,772	3068,96	34081,88	66035,43	155593,9
0,7	10,234	2895,33	24644,21	57833,87	89558,46
0,8	11,696	2595,34	15538,14	24598,34	31724,59
0,9	13,158	1602,68	6337,4565	7126,248	7126,248
1	14,62	0	0	0	0

Необходимо учесть воздействие сосредоточенных массовых сил :

, ;

Построим эпюры , (рис. 8)

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

Для погонных моментов:

Расстояния от точек приложения нагрузок до оси приведения.

Момент считаем положительным, если он действует против часовой стрелки.

Интегрируя эпюру , получаем приведенные моменты от воздействия распределенных нагрузок. Схема расчета имеет вид:

Полученные результаты заносим в таблицу 3:

Таблица 3

qv	qkr	qt	av	akr	at	mz	dM	M
4027,11	502,72	2187,44	1,67127	2,2185	2,3664	438,75654	42399,48
4032,53	464,88	1870,60	1,69219	2,1982393	2,335009	1434,007	1368,9901	41030,49
3952,09	427,05	1578,54	1,713111	2,1779786	2,303619	2203,8936	2659,3053	38371,18
3952,09	427,05	1578,54	1,713111	2,1779786	2,303619	5840,2499	2659,3053	38371,18
3779,82	389,22	1311,25	1,734031	2,1577179	2,272228	6371,3749	3610,3448	34760,84
3584,23	351,39	1068,74	1,754951	2,1374572	2,240837	6780,5438	4297,6997	30463,14
3584,23	351,39	1068,74	1,754951	2,1374572	2,240837	3144,1876	4297,6997	30463,14
3337,71	313,56	851,01	1,775871	2,1171965	2,209446	3383,2196	4771,5346	25691,6
3062,89	275,73	658,05	1,796792	2,0969357	2,178056	3491,9366	5025,7392	20665,86
2769,34	237,90	489,86	1,817712	2,076675	2,146665	3488,2576	5102,522	15563,34
2440,94	200,07	346,45	1,838632	2,0564143	2,115274	3343,7442	4994,1933	10569,15
2024,72	162,24	227,82	1,859553	2,0361536	2,083884	2959,9915	4608,0307	5961,119
1542,45	143,32	177,79	1,870013	2,0260233	2,068188	2226,3231	3791,1959	2169,923
0,00	0,00	0,00	1,880473	2,0158929	2,052493	0	2169,9229	0

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

где - расстояние от цеyнтра тяжести -того бака до оси приведения.

Строим суммарную эпюру (рис. 9)

Проверка правильности построения эпюр нагрузок по крылу.

С эпюры =20592кг.

Определение точки положения поперечной силы в расчетном сечении

Зная поперечную силу и приведенный момент в расчетном сечении(=0.2), можно найти точку приложения поперечной силы по хорде крыла расчетного сечения:

Координату откладывают от оси приведения.

Проектировочный расчет сечения крыла

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

Таблица 4

Шаг стрингеров находят из условия получения волнистости поверхности крыла не выше определенного значения. Величина должна удовлетворять неравенству

Здесь и – давление в горизонтальном полете на нижней и верхней поверхностях крыла;

– коэффициент Пуансона, для дюраля ;

– модуль упругости первого рода материала обшивки.

Приближенно величины и считаем равными

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

Задаваясь шагом стрингеров, найдём толщину обшивки, удовлетворяя неравенство (табл. 5).

Таблица 5.

По соображениям прочности увеличим толщину обшивки, приняв

δ сж = 5(мм), δ р = 4(мм),

Определим количество стрингеров на верхней и на нижней частях поперечного сечения: . (рис. 10)

Нагрузки, воспринимаемые панелями будут равны

Нагрузка, воспринимаемая панелью может быть представлена

Подбор продольного силового набора в растянутой зоне

Усилие в растянутой зоне определяется равенством

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

– площадь поперечного сечения одного стрингера,

– толщина обшивки в растянутой зоне.

Так как панель цельнофрезерованная:

– коэффициент, учитывающий концентрацию напряжений и ослабление сечения отверстиями под заклепки или болты,

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

Тогда найдем потребную площадь стрингеров в растянутой панели: рис. 11

Зная потребную площадь стрингера, из сортамента профилей выберем стрингер с близкой площадью поперечного сечения. Выбираем угольник равностенный ПР100-22, , , (рис 11).

Определим площади поясов лонжерона

Площадь следует распределить между растянутыми полками переднего и заднего лонжеронов.

Подбор продольного силового набора в сжатой зоне

Усилие в сжатой зоне находят по формуле:

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

– расчетное разрушающее напряжение стрингера в сжатой зоне,

– площадь поперечного сечения одного стрингера в сжатой зоне,

Присоединенную площадь обшивки определим по формуле:

Тогда потребная площадь стрингера:

Зная потребную площадь стрингера, из сортамента профилей выберем стрингер с близкой площадью поперечного сечения (Рис. 12). Это бульбоугольник ПР102-23, , , . Рис. 12

Критические напряжения местной потери устойчивости выбранного стрингера определим по формуле:

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

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

для полки стрингера:

Так как >, их необходимо скорректировать по формулам:

, , ,

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

Площадь присоединенной обшивки:

Суммарная площадь полок лонжеронов:

Распределим площадь между сжатыми полками переднего и заднего лонжеронов пропорционально квадратам их высот:

Примем отношение ширины полки лонжерона к ее толщине , тогда

1лонжерон:

, ; , ;

2лонжерон:

, ; , .

Подбор толщин стенок лонжеронов

Определим моменты инерции лонжеронов.

Перенося поперечную силу со статическим нулем в центр жесткости, замечаем, что эта сила эквивалентна двум силам:

и крутящему моменту

Эти силы вызывают потоки касательных усилий в стенках лонжеронов (рис. 13) .

Если предположить, что крутящий момент воспринимается только внешним контуром сечения крыла, то этот момент уравновешивается потоком касательных усилий

Тогда в зависимости от расположения поперечной силы (до или после центра жесткости)

Найдем толщину стенки:

, ,

. .

Определение расстояния между нервюрами

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

Критические напряжения потери устойчивости стрингера определяются по формуле:

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

– расстояние между нервюрами.

Проверочный расчет крыла

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

Для определения коэффициента редукции нулевого приближения построим диаграмму деформирования материалов обшивки, стрингеров и лонжеронов. Параметры деформирования приведены в таблице 4.

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

сжатой зоне :

Лонжерон : ,

Стрингер: .

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

Лонжерон: ,

Стрингер: .

Определим редуцированные площади элементов. Действительные площади элементов сечения:

Редуцированные площади:

Дальнейшие расчеты представлены в таблице 6.

Далее необходимо найти координаты центра тяжести редуцированного сечения. Определяем положение центральных осей редуцированного сечения. Исходные оси выбираем проходящими через носок профиля в соответствии с его геометрией (рис. 15).

Координаты центра тяжести редуцированного сечения определяем следующим образом:

где - число сосредоточенных площадей в сечении.

Координаты сосредоточенных элементов в центральных осях найдем так:

Определяем осевые и центробежные моменты инерции редуцированного сечения в центральных осях:

Вычислим координаты элементов в главных центральных осях

. (табл 6)

Определяем моменты инерции в главных центральных осях

Определяем проекции изгибающих моментов на главные центральные оси (рис. 17):

Определяем редуцированные напряжения в элементах сечения:

Определяем действительные напряжения в продольных элементах из условия равенства деформации действительных и редуцированных сечений по диаграмме деформирования (рис. 18).

После нахождения действительных напряжений определяем коэффициент редукции последующего приближения для каждого элемента конструкции:

Определение коэффициентов редукции последующих приближений для каждого элемента конструкции будет проведено с помощью ЭВМ. (приложение 1)

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

В растянутой зоне, - в сжатой зоне.

Таблица 5

Таблица 5 (продолжение)

Проверочный расчет на касательные напряжения

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

и нормальные напряжения , которые равны .(табл. 7)

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

Расстояние между нервюрами, - шаг стрингеров.

Если обшивка теряет устойчивость от сдвига () и работает как диагонально – растянутое поле (рис. 19), то в ней возникают дополнительные растягивающие нормальные напряжения, определяемые по формуле:

где – угол наклона диагональных волн.

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

. .

Условие прочности, соответствующее критерию энергии формообразования, имеет вид:

Коэффициент , характеризующий избыток прочности обшивки определяем по формуле:

Полученные результаты заносим в таблицу 7.

Строим эпюру касательных напряжений (рис. 20)

Таблица 7

Расчет центра жесткости сечения крыла

Центр жесткости – это точка, относительно которой происходит закручивание контура поперечного сечения, либо это точка, при приложении поперечной силы в которой закручивание контура не происходит. В соответствии с этими двумя определениями существуют 2 метода расчета положения центра жесткости: метод фиктивной силы метод фиктивного момента. Так как проверочный расчет на касательные напряжения проведен, и эпюра суммарных ПКУ построена, то для расчета центра жесткости сечения используем метод фиктивного момента.

Определяем относительный угол закручивания 1 го контура. Эпюра q S - известна.

В соответствии с формулой Мора к первому контуру прикладываем единичный момент:

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

Определяем относительный угол закручивания сечения крыла при приложении к нему момента М = 1 ко всему контуру. Неизвестными являются q 01 q 02 , для их определения запишем два уравнения: уравнение равновесия относительно т.А (нижний пояс переднего лонжерона) и уравнение равенства относительных углов закручивания первого и второго контуров (аналог ур-я совместности деформации).

где - удвоенные площади контуров.

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

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

Решая которые, находим

После нахождения `М 1 и`М 2 , определяем относительный угол закручивания первого контура, от приложения к сечению единичного момента:

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

Определяем расстояние от поперечной силы до центра жесткости (рис. 21).

м.

Эксплуатационная работа, поглощаемая амортизационной системой при посадке:

где - эксплуатационная вертикальная посадочная скорость, равная

Но так как , то принимаем м/с.

кДж.

Одна стойка воспринимает эксплуатационную работу

кДж.

Вычислив эксплуатационную работу, поглощенную пневматиками при посадке

найдем работу воспринимаемую амортизатором

Ход амортизатора вычисляем по формуле

Коэффициент полноты диаграммы обжатия амортизатора при восприятии работы .

φ э - передаточное число при ходе поршня S э.

Так как рассматривается телескопическая стойка и при этом предполагается, что в момент касания колесами земли ось стойки перпендикулярна поверхности земли, то η е =0,7 и φ э =1.

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

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

Зададимся значениями параметров:

МПа – начальное давление газа в амортизаторе;

– коэффициент предварительной затяжки амортизатора;

– передаточное число в момент начала обжатия амортизатора;

м 2 .

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

м.

Толщину уплотнительных колец полагаем .Тогда для внутреннего диаметра цилиндра

Начальный объем V 0 газовой камеры находим по формуле

Высота газовой камеры при необжатом амортизаторе

м.

Параметры и находим по следующему алгоритму.

Для нахождения неизвестных и используем уравнения

После некоторых преобразований

Здесь - передаточное число соответствующее ходу амортизатора

Коэффициент полноты диаграммы обжатия амортизатора при поглощении работы . Для телескопических стоек .

Первое из равенств (3) имеет вид квадратного уравнения

, 5

где , 6

из равенства (5)

Подставляя из (8) во второе уравнение (3) получаем трансцендентное уравнение

корень которого есть искомая величина .

Вычисления сведены в табл. 8

Таблица 8.

Строим график в координатной системе (S max , f) (рис. 22).

Точка пересечения кривой с осью f = 0 дает значение S max =0,55.

Из зависимости (8) найдём

Давление газа в амортизаторе при его максимальном обжатии

МПа.

Высота уровня жидкости над верхней буксой

м.

При этом:

0,589 + 0,1045 = 0,6935 > 0,55 – условие выполняеться.

Задаваясь значениями параметров:

м - конструктивный ход амортизатора;

м - суммарная высота букс;

м - опорная база штока;

м - суммарный размер узлов крепления амортизатора;

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

Длина амортизатора при эксплуатационном обжатии

Определение нагрузок на стойку

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

Расчетная вертикальная и горизонтальная нагрузки на стойку равны:

Между колесами усилие распределяется в соотношении 316,87: 210,36, а усилие - 79,22: 52,81.

Построение эпюр изгибающих моментов

Стойка является комбинированной системой. Вначале методом сечений находим усилие в подкосе. Записываем для стойки уравнение равновесия относительно шарнира

Эпюра изгибающих моментов, действующих в плоскости движения самолёта, изображена на рисунке 23.

Максимальный момент, равный 489,57кНм, действует в точке навески шасси.

Эпюра изгибающих моментов, действующих в плоскости перпендикулярной плоскости движения самолёта, изображена на рисунке 24.

Скачек на эпюре в точке присоединения стержня к цилиндру, созданный эксцентриситетно приложенной силой (вертикальной проекцией усилия в стержне), равен кНм.

Крутящий момент равен величине

и нагружает только цилиндр.

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

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

где k – коэффициент пластичности, принимаем ;

W – момент сопротивления

, ;

МПа.

Из этого уравнения находим

Зная наружный диаметр штока получим внутренний

Тогда толщина стенки .

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

Построение эпюры осевой силы

Расчетное давление газа в амортизаторе

Газ давит на шток с силой

Несоответствие между силой Р ш и внешней нагрузкой 528,127 кН объясняется наличием сил трения в буксах. Таким образом, сила трения в одной буксе равна величине

кН.

На верхнем конце штока газ давит на шток с силой

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

ниже сечения нижней буксы – силой

На цилиндр газ воздействует через уплотнение с осевой силой

растягивающей цилиндр. При построении эпюры N ц, следует учесть также силы F тр и S z . Окончательный вид эпюр осевых сил N ц и N ш показан на рис. 25

К сожалению, я ненашел ни одной статьи по аэродинамики "для моделиста". Ни на форумах, ни в дневниках, ни в блогах- ни где нет нужной "выжимки" по этой теме. А вопросов возникает море, особенно у новичков, да и те, кто считает себя "уже не новичком", зачастую не утруждают себя изучением теории. Но мы это исправим!)))

Сразу скажу, сильно углубляться в эту тему не буду, иначе это получится, как минимум научный труд, с кучкой непонятных формул! И тем более я не стану пугать вас такими терминами, как "число Рейнольдса"- кому будет интересно- можете почитать на досуге.

Итак, договорились- только самое нужное для нас- моделистов.)))

Силы, действующие на самолет в полете.

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

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

Любое изменение любой из этих сил повлияет на характер полета самолета. Если бы подъемная сила, создаваемая крыльями, увеличилась по сравнению с силой тяжести, результатом оказался бы подъем самолета вверх. Наоборот, уменьшение подъемной силы против силы тяжести вызвало бы снижение самолета, т. е. потерю высоты.

Если равновесие сил не будет соблюдаться, то самолет будет искривлять траекторию полета в сторону преобладающей силы.

Про крыло.

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

Удлинение крыла - отношение размаха крыла к его средней аэродинамической хорде. Для непрямоугольного крыла удлинение = (квадрат размаха)/площадь. Это можно понять, если за основу возьмём прямоугольное крыло, формула будет проще: удлинение = размах/хорду. Т.е. если крылоимеет размах 10 метров а хорда = 1 метр, то удлинение будет = 10.

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

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

Аэродинамическим качеством крыла называется отношение подъемной силы к силе лобового сопротивления крыла.

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

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

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

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

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

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

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

Крыло обратной стреловидности (КОС) — крыло с отрицательной стреловидностью.

Преимущества:

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

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

Су-47 "Беркут" с обратной стреловидностью:

Чехословацкий планер LET L-13 с обратной стреловидностью крыла:

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

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

Средней аэродинамической хордой крыла (САХ) называется хорда такого прямоугольного крыла, которое имеет одинаковые с данным крылом площадь, величину полной аэродинамической силы и положение центра давления (ЦД) при равных углах атаки. Или проще- Хорда — отрезок прямой, соединяющей две наиболее удаленные друг от друга точки профиля.

Величина и координаты САХ для каждого самолета определяются в процессе проектирования и указываются в техническом описании.

Если величина и положение САХ данного самолета неизвестны, то их можно определить.

Для крыла, прямоугольного в плане, САХ равна хорде крыла.

Для трапециевидного крыла САХ определяется путем геометрического построения. Для этого крыло самолета вычерчивается в плане (и в определенном масштабе). На продолжении корневой хорды откладывается отрезок, равный по величине концевой хорде, а на продолжении концевой хорды (вперед) откладывается отрезок, равный корневой хорде. Концы отрезков соединяют прямой линией. Затем проводят среднюю линию крыла, соединяя прямой середины корневой и концевой хорд. Через точку пересечения этих двух линий и пройдет средняя аэродинамическая хорда (САХ).

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

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

Разработка аэродинамических профилей проводилась практически с начала истории авиации, проводится она и сейчас.Делается это в специализированных учреждениях. Ярчайшим представителем такого рода учреждений в России является ЦАГИ - Центральный аэрогидродинамический институт имени профессора Н.Е. Жуковского. А в США - такие функции выполняет Исследовательский центр в Лэнгли (подразделение NASA).

THE END?

Продолжение следует.....