CC BY
47
Б01: 10.15593/2224-9982/2017.51.12 УДК 620.17
М.Ш. Нихамкин, Д.Г. Соломонов
Пермский национальный исследовательский политехнический университет, Пермь, Россия
ПРИМЕНЕНИЕ ЭКСПЕРИМЕНТАЛЬНОГО МОДАЛЬНОГО АНАЛИЗА ДЛЯ ИДЕНТИФИКАЦИИ ПАРАМЕТРОВ МОДЕЛИ СЛОИСТОГО
УГЛЕПЛАСТИКА
Расчет модальных характеристик (собственных частот и форм), необходимых для того, чтобы исключить резонансные колебания изделий из полимерных композиционных материалов (ПКМ), требует надежных данных о механических характеристиках материала. Проблема состоит в большом (по сравнению с изотропными материалами) количестве характеристик упругости, а также в том, что эти параметры зависят от широкого круга структурных и технологических факторов. Данные о характеристиках материалов, приведенные в литературе, зачастую противоречивы и при проведении ответственных расчетов требуют дополнительной проверки.
Цель настоящей работы - разработка методики идентификации параметров модели упругого поведения ПКМ по результатам экспериментального модального анализа. Объект исследования - слоистый углепластик на основе равнопрочной углеродной ткани, широко применяемый в авиастроении.
Для экспериментального определения собственных частот и форм колебаний используется метод сканирующей лазерной виброметрии, для расчетного модального анализа - метод конечных элементов. Модель материала - слоистый композит с ортотропными линейно упругими слоями. Задача идентификации параметров модели материала рассматривается как задача минимизации расхождения расчетных собственных частот с экспериментальными. Для ее решения используется метод квазислучайного поиска.
Разработанная методика может быть рекомендована для определения параметров модели материала, необходимых для расчетов модальных характеристик изделий из ПКМ.
Ключевые слова: полимерные композиционные материалы, механические характеристики, модальный анализ, метод конечных элементов.
M.Sh. Nikhamkin, D.G. Solomonov
Perm National Research Polytechnic University, Perm, Russian Federation
APPLICATION OF EXPERIMENTAL MODAL ANALYSIS TO IDENTIFY THE PARAMETERS OF THE MODEL OF LAMINATED CARBON FIBER
REINFORCED PLASTIC
The calculation of modal characteristics (natural frequencies and vibration modes) of products made of polymer composite materials, necessary to exclude their resonance oscillations, requires reliable data on the mechanical characteristics of the material. The problem is a large (in comparison with isotropic materials) quantity of elasticity characteristics, and also in the fact that these parameters depend on a wide range of structural and technological factors. Data on the characteristics of materials, given in the literature, are often contradictory and require additional verification when carrying out responsible calculations.
The purpose of this work is to develop a methodology for identifying models of elastic behavior of polymeric composite materials based on the results of experimental modal analysis using the example of laminated carbon fiber reinforced plastic. The object of investigation is a layered carbon fiber reinforced plastic based on an equally strong carbon cloth widely used in aircraft building.
The method of scanning laser vibrometry is used for experimental determination of natural frequencies and modes of oscillations. For the numerical modal analysis, the finite element method is used. The material model is a layered composite with orthotropic linearly elastic layers. The task of identifying the material model is considered as the problem of minimizing the discrepancy between the calculated natural frequencies and the experimental ones. To solve it, the quasi-random search method is used.
The developed method can be recommended for determination of parameters of material models for calculating modal characteristics of polymer composite products.
Keywords: laminated carbon fiber reinforced plastic, mechanical characteristics, modal analysis, finite element method.
Введение
Полимерные композиционные материалы (ПКМ), в частности углепластики, широко используются для изготовления нагруженных деталей самолетов и авиационных двигателей благодаря возможности облегчения конструкции без снижения прочности [1-4]. Известно, например, что использование ПКМ в конструкции планера самолетов Boeing 787 достигает 50 % по весу [2, 4]. Широкое применение ПКМ предусматривается при создании отечественных самолетов и авиационных двигателей нового поколения [2, 3].
Одна из проблем применения ПКМ в авиастроении - обеспечение прочности в условиях вибраций, характерных для авиационной техники. При значительных амплитудах колебаний, которые могут возникнуть, например, в условиях резонанса, может произойти разрушение элементов конструкции из ПКМ, ведущее к опасным последствиям. Чтобы исключить резонансные колебания, надо еще на стадии проектирования рассчитывать модальные характеристики деталей: их собственные частоты и формы колебаний. В случае изделий из ПКМ такой расчет осложняется отсутствием надежных данных о механических характеристиках. Проблема состоит в большом (по сравнению с изотропными материалами) количестве параметров, входящих в модель материала, а также в том, что эти параметры зависят от материалов волокон и связующего, схемы армирования, технологических факторов [5, 6]. Так, например, модель ортотроп-ного упругого материала содержит девять параметров. Их определение - трудоемкая экспериментальная задача. Данные о характеристиках материалов, приведенные в литературе, зачастую противоречивы и при проведении ответственных расчетов требуют дополнительной проверки.
Для расчетного анализа напряженно-деформированного состояния и модального анализа конструкций, в том числе из ПКМ, эффективно используется метод конечных элементов (МКЭ) (см, например, работы [7-12]). Из-за неопределенностей и допущений, неизбежных при моделировании сложных конструкций, корректное применение МКЭ требует проведения верификации, а при необходимости - идентификации (настройки) конечно-элементных (КЭ) моделей по экспериментальным данным. Под идентификацией КЭ-модели понимается такое изменение ее параметров, которое сводит к минимуму отличия между расчетными и экспериментальными данными. В ряде работ последних лет [11-17] для этого используются результаты экспериментального модального анализа.
В работе [12] корректировка КЭ-модели рабочего колеса турбины проводилась за счет подбора условных значений плотности материала и модуля Юнга, обеспечивающего минимальное расхождение расчетных и экспериментальных значений собственных частот. Коррекция проводилась с помощью коммерческой программы LMS Virtual.Lab Optimization. В работах [13, 14] параметрами, за счет которых выполняется настройка КЭ-модели по собственным частотам, были условные жесткости отдельных конечных элементов. В работах [15, 16] идентификация КЭ-моделей рассматривается как задача многокритериальной многопараметрической оптимизации. Для ее решения предлагается алгоритм, основанный на анализе чувствительности. В работе [17] результаты экспериментального модального анализа использовались для оценки эффективных характеристик свойств сотового заполнителя. Анализ публикаций по проблеме идентификации КЭ-моделей по данным экспериментального модального анализа показывает, что опыт идентификации КЭ-моделей применительно к конструкциям из ПКМ недостаточен.
Цель настоящей работы - разработка методики идентификации параметров модели упругого поведения ПКМ по результатам экспериментального модального анализа на примере слоистого углепластика.
Объект исследования
Объект исследования - слоистый углепластик ВКУ-39 на основе равнопрочной углеродной ткани Porsher 3692 и связующего ВСЭ-1212, широко применяемый в авиастроении [3]. Исследование проводилось на образцах, представляющих собой пластину длиной 250 мм, шири-
ной 24,6 мм и толщиной 4,2 мм (рис. 1). Подобные образцы рекомендованы стандартом ЛБТМ для определения механических характеристик ПКМ при растяжении и усталостных испытаниях [18, 19]. Образец состоит из 20 слоев углепрепрега 0,21 мм со схемой армирования (90/45/0/-45/0)*4. Масса образца 40 г.
Рис. 1. Образец - пластина из углепластика
Методика и результаты экспериментального модального анализа
Задача экспериментального модального анализа (ЭМА) состояла в получении данных о собственных частотах и формах колебаний образца, необходимых для последующей идентификации расчетной модели. Современные средства ЭМА, в частности используемый в настоящей работе метод сканирующей лазерной виброметрии, позволяют не только получать с высокой точностью данные о собственных частотах, но и определять с высоким пространственным разрешением собственные формы колебаний.
Методика ЭМА основана на представлении исследуемого объекта как колебательной системы с конечным числом п степеней свободы. Экспериментальное определение собственных частот системы и соответствующих собственных форм базируется на анализе матрицы передаточных функций [Н], каждый элемент которой представляет собой результат измерений отдельной частотной характеристики как отношения [20, 21]
Ну(ш) = X,(ш)/^(ш), г ,] = 1, ..., п, (1)
где Хг(ю) - частотная функция реакции в виде скорости или ускорения для г-й степени свободы (СС) на воздействие Fj(ш), соответствующее ]-й СС; ш - частота.
При сканирующей лазерной виброметрии к одной СС объекта г прикладывается воздействие в виде внешней силы, а реакция в виде виброскорости измеряется во множестве СС j = 1, ..., п при их последовательным сканировании [21-24]. Далее по соотношению (1) определяются компоненты матрицы передаточных функций Н/(ю). Собственные частоты определяются по пикам на замеренной амплитудно-частотной характеристике (АЧХ) Ну(ш). Для определения собственной формы колебаний одна из степеней свободы принимается опорной и проводится серия измерений амплитуд виброскорости для всех СС при соответствующем значении модальной частоты. Число степеней свободы п при ЭМА выбирается так, чтобы обеспечить достаточное пространственное разрешение для надежного представления форм колебаний.
Принципиальное достоинство метода сканирующей лазерной вибромерии - бесконтактное измерение вибраций, благодаря чему колебательная система не искажается введением в нее дополнительных масс датчиков и жесткости кабелей, препарировка объекта датчиками не требуется. Другое достоинство - высокая пространственная разрешающая способность: плотность точек измерения отклика ограничена только точностью фокусировки лазера [21].
При проведении экспериментального модального анализа исследуемого образца использовалась методика, основанная на применении трехкомпонентного сканирующего лазерного виброметра Ро1у1ее РБУ400-3Б [23-25]. Он представляет собой исследовательский лазерно-цифровой измерительный комплекс, измерение виброскорости в котором основано на исполь-
зовании эффекта Доплера. PSV400-3D состоит из трех оптических сканирующих лазерных головок, модуля сканирования геометрии и системы управления на базе промышленного компьютера.
Исследуемый образец закреплялся на жесткой металлической раме в податливых упругих подвесах (рис. 2). Такая схема, часто используемая в экспериментальном модальном анализе, позволяет максимально приблизить условия закрепления к отсутствию ограничений на перемещения («свободный подвес»), удобному для воспроизведения в последующих расчетах. Возбуждение колебаний осуществлялось с помощью акустического динамика. На него подавался сигнал, изменяющийся во времени по гармоническому закону с постоянной амплитудой и возрастающей до 6400 Гц частотой (Periodic Chirp). Проведено пять экспериментов с различными параметрами сетки сканирования, которая включала в себя от 25 до 165 узлов.
На рис. 3 приведена в качестве примера АЧХ образца. В исследованном диапазоне частот выявлено восемь собственных форм колебаний. В табл. 1 приведены значения собственных частот колебаний, полученные в пяти экспериментах. Рассеяние экспериментальных данных в табл. 1 оценивается значениями коэффициента вариации, который лежит в пределах 0,57 %, что свидетельствует о высокой точности экспериментального определения собственных частот.
Рис. 3. Пример АЧХ образца
Таблица 1
Результаты определения собственных частот колебаний, Гц
№ Номер эксперимента / число точек сканирования Среднее Коэффициент
п/п 1 / 25 2 / 51 3 / 51 4 /165 5 /165 значение вариации, %
1 363 364 360 365 365 363,4 0,571
2 995 998 998 998 998 997,4 0,135
3 1760 1765 1764 1762 1762 1762,6 0,111
4 1935 1943 1941 1941 1944 1940,8 0,180
5 3178 3181 3177 3181 3181 3179,6 0,061
6 3540 3551 3556 3548 3549 3548,8 0,164
7 4694 4702 4690 4688 4700 4694,8 0,130
8 5342 5378 5371 5374 5381 5369,2 0,292
Рис. 2. Экспериментальная установка: 1 - образец; 2 - лазерный виброметр; 3 - акустический возбудитель колебаний
Собственные формы колебаний приведены в табл. 2. Среди них пять изгибных и три крутильных формы. Эти собственные формы используются для последующего сравнения с расчетными при оценке точности расчетов.
Таблица 2
Собственные частоты и формы колебаний пластины: эксперимент и расчет
№ п/п Эксперимент Расчет МКЭ Расхождение, %
Собственная форма колебаний Собственная частота, Гц Собственная форма колебаний Собственная частота, Гц
1-я изгибная 363 389,0 7,05
2-я изгибная 997 1063,2 6,60
1-я крутильная 1763 1744,7 -1,02
3-я изгибная 1940 2057,4 6,01
4-я изгибная V /X ' * 3180 3342,9 5,14
2-я крутильная 3548 3505,2 -1,23
- "" '
5-я изгибная / V /\уч 4695 4890,0 4,16
3-я крутильная 5360 5295,9 -1,37
Методика и результаты расчетного модального анализа
Для расчетного модального анализа использовался метод конечных элементов (МКЭ). При пренебрежении демпфированием собственные колебания конечно-элементной модели е п степенями свободы описываются в матричной форме уравнением [7, 26]
[ М ]{и} + [ К ]{и} = 0,
где [К] и [М] - матрицы жесткости и масс; {и} и {и} - векторы перемещений и ускорений в узлах конечно-элементной (КЭ) модели.
Это уравнение имеет вещественное периодическое решение вида {и} = {и0}ео8ш1, если выполняется условие
([ К ] -ш2[М ]))} = 0.
Задача расчета собственных частот и форм колебаний сводится к задаче о собственных значениях Шк и векторах {и0}к, которые обращают в ноль определитель:
К ] -ш2[ М ]| = 0.
Поставленная задача решалась методом конечных элементов в пакете ANSYS Workbench. Принималось, что образец-пластина свободна от ограничений на перемещения («свободный подвес»). Конечно-элементная
модель состоит из восьмиузловых призматиче- Рис. 4. Конечно-элементная модель образца ских элементов (рис. 4), при этом в рамках одного элемента объединяются все 20 слоев углепластика, ориентированных различным образом. Характерный размер элементов в плоскости пластины выбирался на основе оценки сходимости результатов расчета собственных частот. Серия расчетов с размерами элементов от 5 до 1,25 мм показала, что при размере элемента 2,5 мм (около 1000 элементов) имеет место сходимость значений первых восьми собственных частот с погрешностью в пределах 0,3 %.
Модель материала образца - слоистый композит с ортотропными линейно упругими слоями. Принималось, что различно ориентированные слои композита деформируются совместно без проскальзывания. Технические упругие постоянные углепластика принимались как по результатам испытаний стандартных образцов (для Ец, E22, Gi2), так и по оценочным данным из работ [3, 25, 27] (табл. 3).
Таблица 3
Характеристики* материала ВКУ-39
Источник р, кг/м3 Ец, ГПа Е22, ГПа Е33, ГПа G12, ГПа G13, ГПа G23, ГПа V21 V13 V32
По данным [3, 23, 25] 1580 63,9 63,9 20 4,8 2,7 2,7 0,04 0,3 0,3
После идентификации 1580 70,2 67,4 18,3 4,2 2,9 2,6 0,04 0,3 0,3
* Обозначения: Ец - модуль упругости вдоль основы; Е22 - модуль упругости поперек основы; Е33 - модуль упругости поперек плоскости слоя; 01Ъ 013, 023 - модули сдвига; v21) v1з, Vз2 - коэффициенты Пуассона, р - плотность.
В результате расчетного модального анализа в диапазоне до 6400 Гц выявлено пять изгибных собственных форм колебаний, три крутильных, а также две изгибные в плоскости пластины. В описанном выше экспериментальном исследовании две последние вследствие причин методического характера выявлены не были. Остальные собственные формы совпали с экспериментальными (см. табл. 2). Значения соответствующих собственных частот также приведены в табл. 2. Расхождение расчетных и экспериментальных значений собственных частот лежит в пределах 7 %. В некоторых практически важных случаях такая точность расчета собственных частот может оказаться недостаточной для отстройки от резонансных колебаний детали.
Идентификация параметров модели материала на основе результатов экспериментального модального анализа
Цель идентификации (настройки) модели материала состояла в уточнении характеристик материала: модулей упругости вдоль основы Е11 и поперек основы Е22, модуля упругости поперек плоскости слоя Е33, модулей сдвига С12, 013, 023 и коэффициентов Пуассона v21, v1з, Vз2. Для идентификации использовались результаты экспериментального модального анализа - собственные частоты и формы колебаний образца.
Задача идентификации модели материала рассматривается как задача оптимизации [28, 29] с целевой функцией:
( А - А Л2
А - А
т ^ ¡V Л к
1 = ^ а!
1=1
А
■ шш, (2)
V -Не J
где /¡р и /¡е - расчетное и экспериментальное значения ¡-й собственной частоты; а.. - весовые коэффициенты.
Параметрами управления считаются характеристики материала Е11, Е22, Е33, С12, 013, 023, v21, v13, v32. В качестве ограничений на параметры управления принималось отклонение от исходного значения на ±8 %.
В соотношении (2) важно, чтобы расчетные и экспериментальные значения собственных частот соответствовали одним и тем же собственным формам. В настоящей работе сравнение расчетных и экспериментальных собственных форм проводилось на основании анализа их анимационного представления. В более сложных случаях для этого может быть использован анализ МАС-матрицы [21].
Для решения задачи оптимизации использовался метод квазислучайного поиска. Весовые коэффициенты были приняты равными единице. Квазислучайный поиск включал в себя 100 реализаций.
В результате расчетов был получен набор параметров управления, при котором целевая функция (2) снизилась по сравнению с исходным значением в три раза. Эти значения приведены в табл. 3. Полученные для этих параметров значения собственных частот приведены в табл. 4. Расхождение с экспериментальными данными составляет 3,9 % для крутильных форм колебаний и 2,9 % для изгибных форм, что значительно лучше, чем до идентификации.
Таблица 4
Погрешность расчета собственных частот колебаний образца
№ п/п Эксперимент, собственная частота, Гц Расчет МКЭ без идентификации Расчет МКЭ после идентификации
Собственная частота, Гц Расхождение с экспериментом, % Собственная частота, Гц Расхождение с экспериментом, %
1-я изгибная 363 389,0 7,05 374 2,9
2-я изгибная 997 1063,2 6,60 1022 2,4
1-я крутильная 1763 1744,7 -1,02 1699 -3,8
3-я изгибная 1940 2057,4 6,01 1980 2,0
4-я изгибная 3180 3342,9 5,14 3220 1,2
2-я крутильная 3548 3505,2 -1,23 3414 -3,9
5-я изгибная 4695 4890,0 4,16 4717 0,5
3-я крутильная 5360 5295,9 -1,37 5158 -3,9
На рис. 5 приведена диаграмма, демонстрирующая чувствительность четырех низших собственных частот колебаний к параметрам материала. Видно, что чувствительность собственных частот, а значит, и целевой функции к параметрам материала неодинакова. Наибольшее влияние на все собственные частоты колебаний оказывает модуль упругости в направлении основы Е11, существенно меньшее - Е22. К модулям сдвига целевая функция чувствительна слабо. Влияние модулей сдвига и модуля упругости в направлении утка на крутильную форму колебаний несколько больше, чем на изгибные формы.
Этот вывод можно использовать для уменьшения трудоемкости расчетов по идентификации моделей подобных материалов.
Рис. 5. Оценка чувствительности расчетных собственных частот к параметрам материала:
1 - Еи; 2 - Е22, 3 - Оп; 4 - 023
Заключение
В результате проведенного исследования отработана методика идентификации параметров модели упругого поведения слоистых ПКМ на основании экспериментальных данных о модальных характеристиках (собственных частотах и формах колебаний) образцов. Основные элементы методики сводятся к следующему.
В качестве образца для исследований используется рекомендованная стандартами ЛБТМ для испытаний на растяжение и усталость пластина прямоугольной формы. Образец крепится в податливых эластичных подвесах таким образом, чтобы условия закрепления - отсутствие ограничений на перемещения - были максимально приближены к реализуемым в расчетах.
Для проведения экспериментального модального анализа рекомендуется трехкомпонент-ная сканирующая лазерная виброметрия, с помощью которой настоящей работе удалось получить собственные частоты колебаний образца из ПКМ в диапазоне частот до 6400 Гц с весьма незначительным рассеянием (в пределах 0,5 %), а собственные формы - с высоким пространственным разрешением.
Используемая для идентификации конечно-элементная модель образца должна обеспечивать сходимость результатов расчета собственных частот также в пределах 0,5 %. Для решения задачи оптимизации может быть эффективно использован метод квазислучайного поиска. В качестве параметров управления для слоистых ПКМ на тканой основе можно ограничиться модулями упругости Ец и Е22 .
Разработанная методика может быть рекомендована для определения параметров модели материала, необходимых для расчетов модальных характеристик изделий из ПКМ.
Библиографический список
1. Иноземцев А.А., Нихамкин М.Ш., Сандрацкий В.Л. Основы конструирования авиационных двигателей и энергетических установок. - М., Машиностроение, 2008. - Т. 2. - 368 с. - (Газотурбинные двигатели).
2. Савин С.П. Применение современных полимерных композиционных материалов в конструкции планера самолетов семейства МС-21 // Известия Самар. науч. центра РАН. - 2012. - Т. 14, № 4(2). -С. 686-693.
3. Технологии и задачи механики композиционных материалов для создания лопатки спрямляющего аппарата авиационного двигателя / А.Н. Аношкин, В.Ю. Зуйко, Г. С. Шипунов, А. А. Третьяков // Вестник Пермского национального исследовательского политехнического университета. Механика. -2014. - № 4. - С. 5-44.
4. Келли А. Инженерный триумф углеволокон // Композиты и наноструктуры. - 2009. - № 1. -С. 3S-49.
5. Стреляев Д.В., Умушкин Б.П., Никонов В.В. Перспективные композиционные материалы в конструкциях авиационной и космической техники / Моск. гос. техн. ун-т гражд. авиации. - М., 2012. - 73 с.
6. Скворцов Ю. В. Конспект лекций по дисциплине механика композиционных материалов / Са-мар. гос. аэрокосм. ун-т. - Самара, 2013. - 94 с.
7. Леонтьев Н.В. Применение системы ANSYS к решению задач модального и гармонического анализа / Нижегород. гос. ун-т им. Н.И. Лобачевского. - Н. Новгород, 2006. - 101 с.
S. Расчетно-экспериментальные исследования собственных частот и форм колебаний лопатки спрямляющего аппарата из полимерных композиционных материалов / М.А. Гринев, А.Н. Аношкин, П.В. Писарев, Г. С. Шипунов, М.Ш. Нихамкин, А. А. Балакирев, И.П. Конев, А. А. Головкин // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - № 4. -С. 106-119.
9. Нихамкин М.Ш., Семенова И.В. Концентрация напряжений в лопатках компрессора при повреждении их посторонними предметами // Известия вузов. Авиационная техника. - 2011. - № 4. - С. 15-1S.
10. Ефимик В.А. Применение метода конечных элементов к задаче обственных колебаний прямоугольных пластин и цилиндрических оболочек // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2014. - № 3S. - С. 72-92.
11. Mistuning identification and model updating of an industrial blisk / D. Laxalde, F. Thouverez, J.-J. Sinou, J.-P. Lombard, S. Baumhauer // International Journal of Rotating Machinery. - Vol. 2007. - Article ID 172S9. - 10 p. DOI: 10.1155/2007/172S9
12. Ткач В.В. Применение модального анализа в многодисциплинарном исследовании ЖРД // Электронный журнал «Труды МАИ». - 2010. - № 3S.
13. Безмозгий И.М., Софинский А.Н., Чернягин А.Г. Моделирование в задачах вибропрочности конструкций ракетно-космической техники // Космическая техника и технологии. - 2014. - № 3(6). -С. 71-S0.
14. Межин В.С., Обухов В.В. Практика применения модальных испытаний для целей верификации конечно-элементных моделей конструкции изделий ракетно-космической техники // Космическая техника и технологии. - 2014. - № 1(4). - С. S6-91.
15. Методика уточнения конечно-элементной модели механической системы с помощью анализа чувствительности / С.М. Николаев, В.А. Жулёв, И.А. Киселёв, П.С. Воронов // Наука и образование. Электрон. журн. МГТУ им. Н.Э. Баумана. - 2014. - № 12. - С. 12S-136.
16. Николаев С. М. Жулёв В. А. Киселёв И. А. Уточнение конечно-элементной модели лопатки ГТД на основе результатов вибрационных испытаний с учетом разброса модальных параметров // Наука и образование. Электрон. журнал. МГТУ им. Н.Э. Баумана. - 2015. - № 09. - С. 336-351.
17. Калиненкова А. О. Исследование свойств сотового заполнителя по результатам динамических испытаний // Молодеж. науч.-техн. сб.: электрон. журнал. - 2015. - № 12.
1S. ASTM D 3039/D 3039M-14. Standard Test Method for Tensile Properties of Polymer Matrix Composite Materials / ASTM International. - West Conshohocken, PA, 1942S-2959 USA, 2014.
19. ASTM D 3479. Standard Test Method for Tension-Tension Fatigue of Polymer Matrix Composite Materials / ASTM International. - West Conshohocken, PA, 1942S-2959 USA, 2012.
20. Ewins D.J. Modal testing: Theory, рractice and аpplication. - 2nd ed. - Baldock: Research Studies Press LTD, 2000.
21. Heylen W., Lamens S., Sas P. Modal аnalyses. Theory and testing / Leven. University. Leuven, 2003. - 325 p.
22. Методика экспериментального модального анализа лопаток и рабочих колес газотурбинных двигателей / А.А. Иноземцев, М.Ш. Нихамкин, Л.В. Воронов, А.Б. Сенкевич, А.Ю. Головкин, Б.П. Болотов // Тяжелое машиностроение. - 2010. - № 11. - C. 2-6.
23. Вибрационные измерения в аэрокосмической отрасли [Электронный ресурс] / Advancing Measurements by Ligh. Polytec. - URL: www.polytec.com (дата обращения: 20.09.2017).
24. Экспериментальный и расчетный модальный анализ лопаток вентиляторов полой конструкции / А.А. Иноземцев, М.Ш. Нихамкин, Л.В. Воронов, И.Л. Гладкий, А.Ю. Головкин, Б.П. Болотов // Авиационная промышленность. - 2010. - № 3. - C. S-11.
25. Компьютерное моделирование механического поведения композитной лопатки спрямляющего аппарата авиационного двигателя / М.А. Гринев, А.Н. Аношкин, П.В. Писарев, В.Ю. Зуйко, Г.С. Шипунов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 3. - С. 38-51.
26. Зенкевич О. Метод конечных элементов в технике. - М.: Мир. - 542 с.
27. Исследование НДС и оценка прочности композитной лопатки спрямляющего аппарата авиационного двигателя / М.А. Гринёв, А.Н. Аношкин, П.В. Писарев, В.Ю. Зуйко, Г.С. Шипунов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 4. -С. 293-307.
28. Колмогоров Г. Л., Лежнева А. А. Оптимальное проектирование конструкций: учеб. пособие / Перм. гос. техн. ун-т. - Пермь, 2005. - 168 с.
29. Алексеева Е.В., Кутненко О.А., Плясунов А.В. Численные методы оптимизации: учеб. пособие / Новосиб. ун-т. - Новосибирск, 2008. - 128 с
References
1. Inozemtsev A.A., Nikhamkin M.Sh., Sandratskiy V.L. Osnovy konstruirovaniya aviatsionnykh dvigateley i energeticheskikh ustanovok [Basics of designing aircraft engines and power plants]. Seriya: Gazo-turbinnye dvigateli t. 2. M., Mashinostroenie, 2008, 368 p.
2. Savin S.P. Primenenie sovremennykh polimernykh kompozitsionnykh materialov v konstruktsii planera samoletov semeystva MS-21 [The use of modern polymer composite materials in the design of the airframe of the MS-21 family aircraft]. Izvestiya Samarskogo nauchnogo tsentra Rossiyskoy akademii nauk, 2012, Vol. 14, no. 4(2), pp. 686-693.
3. Anoshkin A.N., Zuyko V.Yu., Shipunov G.S., Tretyakov A.A. Tekhnologii i zadachi mekhaniki kompozitsionnykh materialov dlya sozdaniya lopatki spryamlyayushchego apparata aviatsionnogo dvigatelya [Technologies and tasks of the mechanics of composite materials for the creation of a blade of a rectifying device for an aircraft engine]. Vestnik PNIPU. Mekhanika, 2014, no. 4, pp. 5-44.
4. Kelli A. Inzhenernyy triumf uglevolokon. Kompozity i nanostruktury [Engineering triumph of carbon fiber. Composites and nanostructures.]. 2009, no. 1, pp. 38-49.
5. Strelyaev D.V., Umushkin B.P., Nikonov V.V. Perspektivnye kompozitsionnye materialy v konstruktsiyakh aviatsionnoy i kosmicheskoy tekhniki [Perspective composite materials in the constructions of aviation and space technology]. Moscow, MGTU GA, 2012, 73 p.
6. Skvortsov Yu.V. Konspekt lektsiy po distsipline mekhanika kompozitsionnykh materialov [Lecture notes on the discipline of the mechanics of composite materials]. Samara. SGAU, 2013, 94 p.
7. Leontev N.V. Primenenie sistemy ANSYS k resheniyu zadach modalogo i garmonicheskogo analiza [Application of the ANSYS system to solving modal and harmonic analysis problems]. Nizhniy Novgorod, NNGU, 2006, 101 p.
8. Grinev M.A., Anoshkin A.N., Pisarev P.V., Shipunov G.S., Nikhamkin M.Sh, Balakirev A.A., Konev I.P., Golovkin A.A. Raschetno-eksperimentalnye issledovaniya sobstvennykh chastot i form kolebaniy lopatki spryamlyayushchego apparata iz polimernykh kompozitsionnykh materialov [Calculation and experimental studies of the natural frequencies and vibration modes of the blade of a rectifying apparatus made from polymer composite materials]. Vestnik PNRPU. Mekhanika, 2016, no. 4, pp. 106-119.
9. Nikhamkin M.Sh., Semenova I.V. Kontsentratsiya napryazheniy v lopatkakh kompressora pri povrezhdenii ikh postoronnimi predmetami [Concentration of stresses in the compressor blades if they are damaged by foreign objects]. Izvestiya vysshikh uchebnykh zavedeniy. Aviatsionnaya tekhnika, 2011, no. 4, pp. 15-18.
10. Efimik V.A. Primenenie metoda konechnykh elementov k zadache sobstvennykh kolebaniy prya-mougolnykh plastin i tsilindricheskikh obolochek [The application of the finite element method to the problem of general oscillations of rectangular plates and cylindrical shells]. Vestnik PNRPU. Aerokosmicheskaya tekhnika, 2014, no. 38, pp. 72-92.
11. Laxalde D., Thouverez F., Sinou J.-J., Lombard J.-P., Baumhauer S. Mistuning identification and model updating of an industrial blisk. International Journal of Rotating Machinery, vol. 2007, Article ID 17289, 10 p. DOI: 10.1155/2007/17289
12. Tkach V.V. Primenenie modalnogo analiza v mnogodistsiplinarnom issledovanii ZhRD [The use of modal analysis in a multidisciplinary study]. Elektronnyy zhurnal «Trudy MAI», 2010, no. 38.
13. Bezmozgiy I.M., Sofinskiy A.N., Chernyagin A.G. Modelirovanie v zadachakh vibroprochnosti kon-struktsiy raketno-kosmicheskoy tekhniki [Modeling in the problems of vibration resistance of rocket and space engineering structures]. Kosmicheskaya tekhnika i tekhnologii, 2014, no. 3 (6)/2014, pp.71-80.
14. Mezhin V.S., Obukhov V.V. Praktika primeneniya modalnykh ispytaniy dlya tseley verifikatsii konechno-elementnykh modeley konstruktsii izdeliy raketno-kosmicheskoy tekhniki [The practice of applying modal tests for the purposes of verification of finite element models of the design of rocket and space equipment]. Kosmicheskaya tekhnika i tekhnologii, 2014, no. 1(4)/2014, pp. 86-91.
15. Nikolaev S.M., Zhulev V.A., Kiselev I.A., Voronov P.S. Metodika utochneniya konechno-elementnoy modeli mekhanicheskoy sistemy s pomoshchyu analiza chuvstvitelnosti. [Method for refining the finite element model of a mechanical system using sensitivity analysis]. Nauka i Obrazovanie. MGTU im. N.E. Baumana. Elektron. Zhurn, 2014, no. 12, pp. 128-136.
16. Nikolaev S.M., Zhulev V.A., Kiselev I.A. Utochnenie konechno-elementnoy modeli lopatki GTD na osnove rezul'tatov vibratsionnykh ispytaniy s uchetom razbrosa modalnykh parametrov [Refinement of the finite-element model of the GTE blade based on the results of vibration tests taking into account the spread of the modal parameters]. Nauka i Obrazovanie. MGTU im. N.E. Baumana. Elektron. Zhurn, 2015, no. 09, pp. 336-351.
17. Kalinenkova A.O. Issledovanie svoystv sotovogo zapolnitelya po rezultatam dinamicheskikh ispy-taniy [Investigation of properties of honeycomb core by results of dynamic tests]. Molodezhnyy nauchno-tekhnicheskiy sbornik. Elektronnyy zhurnal, 2015, no. 12.
18. ASTM D 3039/D 3039M - 14. Standard Test Method for Tensile Properties of Polymer Matrix Composite Materials. ASTM International, West Conshohocken, PA, 19428-2959 USA, 2014.
19. ASTM D 3479. Standard Test Method for Tension-Tension Fatigue of Polymer Matrix Composite Materials. ASTM International, West Conshohocken, PA, 19428-2959 USA, 2012.
20. Ewins D.J. Modal Testing: Theory, Practice and Application// 2nd edition. Baldock, Research Studies Press LTD, 2000.
21. Heylen W., Lamens S., Sas P. Modal analyses. Theory and testing. Leven Univ. Belgium, 2003,
325 p.
22. Inozemtsev A.A., Nikhamkin M.Sh., Voronov L.V., Senkevich A.B., Golovkin A.Yu., Bolotov B.P. Metodika eksperimentalnogo modalnogo analiza lopatok i rabochikh koles gazoturbinnykh dvigateley [Technique of experimental modal analysis of blades and impellers of gas turbine engines]. Tyazheloe mashinostroe-nie, 2010, no. 11, pp. 2-6.
23. Vibratsionnye izmereniya v aerokosmicheskoy otrasli [Vibration measurements of the aerospace industry]. Advancing Measurements by Ligh. Polytec, available at: www.polytec.com (accessed 20 September 2017).
24. Inozemtsev A.A., Nikhamkin M.Sh., Voronov L.V., Gladkiy I.L., Golovkin A.Yu., Bolotov B.P. Eksperimentalnyy i raschetnyy modalnyy analiz lopatok ventilyatorov poloy konstruktsii [Experimental and calculated modal analysis of fan blades of hollow structure]. Aviatsionnayapromyshlennost, 2010, no. 3, pp. 8-11.
25. Grinev M.A., Anoshkin A.N., Pisarev P.V., Zuyko V.Yu., Shipunov G.S. Kompyuternoe modeliro-vanie mekhanicheskogo povedeniya kompozitnoy lopatki spryamlyayushchego apparata aviatsionnogo dvigate-lya [Computer modeling of the mechanical behavior of a composite blade of an aircraft engine straightener]. Vestnik PNIPU. Mekhanika, 2015, no. 3, pp. 38-51.
26. Zenkevich O. Metod konechnykh elementov v tekhnike [Finite element method in engineering]. Moscow: Mir, 542 p.
27. Grinev M.A., Anoshkin A.N., Pisarev P.V., Zuyko V.Yu., Shipunov G.S. Issledovanie NDS i ot-senka prochnosti kompozitnoy lopatki spryamlyayushchego apparata aviatsionnogo dvigatelya[Research of VAT and an estimation of durability of a composite blade of rectifying device of the aviation engine]. Vestnik PNIPU. Mekhanika, 2015, no. 4, pp. 293-307.
28. Kolmogorov G.L., Lezhneva A.A. Optimalnoe proektirovanie konstruktsiy [Optimal design of structures: Textbook]. PGTU. Perm, 2005, 168 p.
29. Alekseeva E.V., Kutnenko O.A., Plyasunov A.V. Chislennye metody optimizatsii [Numerical optimization methods: Textbook]. Novosib. un-t. Novosibirsk, 2008, 128 p.
Сведения об авторах
Нихамкин Михаил Шмерович (Пермь, Россия) - доктор технических наук, профессор кафедры «Авиационные двигатели» ФГБОУ ВО ПНИПУ (614990, Пермь, Комсомольский пр., д. 29, e-mail: nikhamkin @mail .ru).
Соломонов Данил Глебович - студент ФГБОУ ВО ПНИПУ (614990, Пермь, Комсомольский пр., д. 29, e-mail: solomonov1198@yandex.ru).
About the authors
Mikhail Sh. Nikhamkin (Perm, Russian Federation) - Doctor of Technical Sciences, Professor, Department of Aircraft Engines, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: nikhamkin@mail.ru).
Danil G. Solomonov (Perm, Russian Federation) - Student, Perm National Research Polytechnic University (29, Komsomolsky av., Perm, 614990, Russian Federation, e-mail: solomonov1198@yandex.ru).
Получено 18.11.2017
