Scientific journal
Scientific Review. Technical science
ISSN 2500-0799
ПИ №ФС77-57440

CONSTRUCTION AND VALIDATION OF A FINITE ELEMENT MODEL OF A TWO-BLADED AUTOROTATING MAIN ROTOR

Korbukov G.M. 1 Menzulskiy S.Yu. 1
1 Bauman Moscow State Technical University
1247 KB
The article presents the results of constructing a finite element model of a two-bladed autorotating main rotor, its validation and calculation for its own oscillations. The conducted experiment to determine the stiffness characteristics of the prototype main rotor is described, boundary conditions and loads are indicated, and the obtained values of deflection and deformations of the rotor blade are presented. A detailed description of the process of constructing and structuring the resulting finite element model is given. A static analysis of the main rotor model was carried out in order to determine the stiffness characteristics. Validation of the model was performed, in which the parameters of the properties of the material and connecting nodes were selected in order to compare the calculated and experimental data of the static calculation. During validation, errors in the values of deformations of the rotor blade between the calculated model and the experiment were obtained, from which conclusions were drawn about the convergence of the validated finite element model. The calculation of the finite element model of the main rotor for its own vibrations is performed, the natural frequencies and forms of vibrations of the finite element model of the main rotor are determined. The minimum permissible rotational speed of a two-bladed autorotating main rotor at the pre-spin stage has been determined. A conclusion was reached on the further use of the model in subsequent tests and calculations.
autorotating rotor
finite element model
validation
eigenfrequencies
eigenforms

С 2020 г. МГТУ им. Н.Э. Баумана проводятся работы по созданию наземно-воздушных вездеходов (НВВ), способных перемещаться по земле и в воздухе. В процессе предпроектного анализа НВВ было принято решение об использовании аэродинамической схемы автожира и, соответственно, установке на прототип авторотирующего несущего винта (НВ) с общим горизонтальным шарниром (ГШ). Динамика конструкции несущего винта всегда являлась одной из основных проблем винтокрылых аппаратов [1, с. 15–18; 2, с. 180–182]. В подобных задачах отправной точкой является определение собственных частот колебания конструкции и построение валидированной конечно-элементной модели для дальнейших расчетных исследований [3, с. 82–88; 4].

В рамках инициативной научно-исследовательской работы в МГТУ им. Н.Э. Баумана был изготовлен специализированный экспериментальный стенд на базе мотовездехода BRP Outlander 6x6 с установленным авторотирующим несущим винтом от автожира «Казачок».

Помимо ходовых испытаний специализированный стенд был использован для экспериментальных исследований напряженно-деформированного состояния несущего винта мотовездехода (рис. 1). Жесткостные испытания заключались в тензометрировании лопасти и комлевой части НВ, нагружении концевой хорды лопасти грузом весом 99,47 Н и измерении возникающих вдоль размаха лопасти линейных деформаций.

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

Результаты жесткостных испытаний, осредненные по трем измерениям, представлены в табл. 1.

Более подробно процесс проведения и результаты испытаний наглядно представлены в статье [5].

Построение модели

Построение конечно-элементной модели авторотирущего несущего винта должно выполняться в соответствии с упруго-массовыми свойствами его конструкции [6, с. 379–381; 7, 82–83].

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

missing image file

Рис. 1. Жесткостные испытания несущего винта специализированного стенда

missing image file

Рис. 2. Схема установки тензодатчиков

Таблица 1

Результаты прогибов и деформаций несущего винта специализированного стенда

№ ступени нагружения

Вес груза, Н

Прогиб лопасти, мм

Показания тензодатчика № 1

Показания тензодатчика № 2

Показания тензодатчика № 3

1

0,00

0

0,000000

0,000000

0,000000

2

34,43

215

-0,000201

-0,000134

-0,000053

3

70,04

449

-0,000405

-0,000268

-0,000105

4

99,47

643

-0,000571

-0,000375

-0,000147

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

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

Центральный болт установлен в центральном отверстии центральной части балки, моделируется как bush-элемент, соединенный с краями отверстия в центральной балке через MPC-элементы.

Лопасть моделируется трехмерными тетра-призматическими и гексаэдрическими конечными элементами первого порядка. Для создания использовалась плоская геометрическая модель сечения лопасти. На ней создавалась сетка из двумерных треугольных и квадратных конечных элементов первого порядка, а затем эта сетка вытягивалась по нормали к плоскости сечения, образуя трехмерную сетку. Каждая лопасть крепится к боковой части несущей балки с помощью 14-и MPC-элементов (по 7 на верхней и нижней сторонах лопасти), соединенных с отверстиями в балке.

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

Таким образом, для расчета динамики конструкции несущего винта была создана конечно-элементная модель размерностью 265 тысяч узлов.

Валидация модели

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

В качестве основного материала всей конструкции принимается алюминий. Исходные параметры: модуль упругости 1 рода 71 ГПа, коэффициент Пуассона 0,3, плотность материала 2700 кг/м3.

Для согласования результатов статического анализа полученной модели с результатами эксперимента к ней приложены следующие граничные условия. На центральный болт наложены ограничение перемещений по всем осям и ограничение вращения вокруг вертикальной и продольной осей. На одной лопасти наложено ограничение перемещения узлов нижней дальней кромки по вертикальной оси. На конце другой лопасти приложена распределенная по узлам нагрузка с суммарным значением 99,47 Н вдоль вертикальной оси (рис. 3).

missing image file

Рис. 3. Конечно-элементная модель, граничные условия и схема нагружения конечно-элементной модели

Таблица 2

Расчетные величины прогибов и линейных деформаций несущего винта специализированного стенда

Значения

Прилагаемая нагрузка

Прогиб лопасти

Линейные деформации в точке № 1

Линейные деформации в точке № 2

Линейные деформации в точке № 3

Экспериментальные

99,47, Н

645, мм

-0,000586

-0,000371

-0,000169

Расчетные

99,47, Н

837, мм

-0,000571

-0,000382

-0,000189

Погрешность

0 %

22,9 %

2,6 %

2,9 %

10,6 %

При проведении статического анализа значения деформации конструкции определялись в точке размещения тензодатчика № 1 (на расстоянии 2789 мм от дальнего конца и 157 мм от задней части лопасти). Для получения значений, совпадающих со значениями из эксперимента, необходимо менять значение модуля упругости материала. В результате изменения этих свойств были получены требуемые значения деформации, данные приведены в табл. 2.

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

В процессе валидации конечно-элементной модели было обнаружено, что в зависимости от вариации параметров материалов, модель показывает высокую сходимость с экспериментом либо по величинам деформации, либо по значению прогиба кончика лопасти. Это связано с тем, что при высоких значениях прогиба концевой хорды лопасти конструкция становится геометрически нелинейной. Так как собственные колебания – это малые колебания вокруг положения статического равновесия, рационально [9, с. 44] проводить валидацию по значениям экспериментально измеренных деформаций конструкции.

Расчет погрешности конечно-элементной модели выполнен по формуле

missing image file,

где missing image file – значение среднего отклонения, хэ – экспериментально измеренная линейная деформация, хр – расчетная линейная деформация, n – количество валидируемых точек. По результатам расчетов missing image file составляет 5,4 %.

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

На рис. 4 представлено деформированное состояние конечно-элементной модели.

Результаты расчета собственных колебаний несущего винта

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

Для примера на рис. 5 приведена пятая форма собственных колебаний несущего винта.

missing image file

Рис. 4. Результаты статического анализа конечно-элементной модели

Таблица 3

Результаты модального анализа несущего винта на горизонтальном шарнире

Форма колебания

Частота, Гц

Частота, об/мин

1

Колебания несущего винта как твердого целого на горизонтальном шарнире

0

0

2

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

1,43

85,8

3

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

2,75

165

4

Первая горизонтальная симметричная

3,14

188,4

5

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

5,52

331,2

6

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

8,69

521,4

7

Третья вертикальная антисимметричная

17,5

1050

missing image file

Рис. 5. Пятая форма колебаний несущего винта

Выводы

1. Разработана конечно-элементная модель двухлопастного авторотирующего несущего винта. Модель провалидирована по результатам жесткостных испытаний. Получена сходимость 94,6 % по величинам экспериментальных и расчетных значений деформации конструкции.

2. Проведен расчет собственных частот и форм колебаний двухлопастного авторотирующего несущего винта. Показано, что система предраскрутки несущего винта должна обеспечивать частоту вращения НВ не менее 165 об/мин.

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