Full Text

Введение В настоящее время композиционные материалы активно внедряются в конструкцию силовых деталей и узлов современных авиационных двигательных установок [1]. Лопатки вентилятора либо лопатки спрямляющего аппарата перспективных зарубежных двигателей нового поколения: GEnx, Rolls-Royce семейства Advance и UltraFan, LEAP и HF120 выполнены из полимерных композиционных материалов (ПКМ). Разработка аналогичной лопатки спрямляющего аппарата (ЛСА) для российских двигателей является чрезвычайно актуальной. Ранее в работах [2-4] представлены результаты выбора материала и оценки прочности композитной ЛСА. Показано, что лопатка, выполненная из углепластика, обеспечивает снижение веса на до 40 %, имеет необходимую жесткость и обладает более высоким запасом статической прочности по сравнению с металлическим аналогом. В процессе эксплуатации ЛСА подвергается воздействию интенсивных аэродинамических нагрузок, которые могут приводить к колебаниям, усиливающимся при появлении резонанса, и как следствие - к разрушению. В связи с этим при проектировании ЛСА, так же как и при проектировании подобных конструкций, в том числе из композиционных материалов, [5-10] возникает необходимость в проведении модального анализа для определения собственных частот и форм колебаний, что и являются целью расчетно-экспериментальных исследований, представленных в настоящей работе. Кроме того, модальный анализ часто используется для верификации математических моделей оценки состояния и мониторинга конструкций [11-16], поскольку повреждения и трещины меняют их спектр собственных частот колебаний. а б Рис. 1. Общий вид образца ЛСА из ПКМ в исходном состоянии (а) и с нанесенным покрытием для лазерного модального анализа (б) Fig. 1. General view of OGV sample in the initial state (a) and with the coating for laser modal analysis (b) Объектом настоящих исследований является лопатка спрямляющего аппарата (рис. 1), выполненная из углепластика по препреговой технологии с автоклавным формованием (изготовленная на уникальной установке «Исследовательский комплекс для проведения научно-технологических исследований в области создания изделий из полимерных композиционных материалов»). Перо, верхняя полка, нижняя полка создаются за один технологический процесс, что обеспечивает наиболее высокие механические характеристики композиционного материала, реализуемые в конструкции [3, 4]. Для оценки влияния возможных технологических отклонений на амплитудно-частотные характеристики (АЧХ) ЛСА из ПКМ экспериментальные исследования проводились на трех натурных образцах композитной лопатки, полностью соответствующих конструктивным требованиям и изготовленных по одной технологии. Для численного расчета частот и форм ЛСА из ПКМ в настоящей работе была разработана математическая модель, которая позволяет в последующих исследованиях оценить изменение АЧХ ЛСА при вариации её схемы армирования и конструктивных особенностей её исполнения. 1. Методика экспериментального определения собственных частот и форм колебаний ЛСА Теоретической основой экспериментального модального анализа является возможность описания динамического поведения исследуемой конструкции компактным набором модальных параметров: собственных частот, форм и декрементов (или коэффициентов демпфирования) колебаний [17, 18]. При этом объект исследования представляют как конечное число n независимых колебательных систем с одной степенью свободы. Динамические свойства колебательной системы с n степенями свободы полностью определяются набором модальных параметров , σr, w0r, (r = 1, 2 … n) где ω0r - собственные частоты недемпфированной системы; - собственные формы; σr - коэффициенты демпфирования. Экспериментальное определение этих параметров базируется на анализе матрицы передаточных функций, которая вводится соотношением {X} = [H] {F}, (1) где {X} - вектор спектров механических колебаний, соответствующих всем степеням свободы системы; {F} - вектор спектров сил возбуждения для тех же самых степеней свободы. Каждый элемент Нij(ω) матрицы [H] представляет собой результат измерений отдельной частотной характеристики как отношения (2) где Xi(ω) - частотная функция реакции (в виде перемещения, скорости или ускорения), соответствующая степени свободы i, на воздействие Fj(ω), соответствующее степени свободы j. При экспериментальных исследованиях, когда измеряются скорости или ускорения, матрица [Н] обычно является матрицей подвижностей или ускоряемостей соответственно [18, 19, 20]. Собственные частоты можно определить по любой замеренной частотной характеристике Нij(ω). Для определения формы колебаний необходимо принять одну из степеней свободы опорной и провести серию измерений для всех степеней свободы. Для определения вектора формы конкретной моды необходимо определить модули замеренных частотных характеристик при соответствующем значении модальной частоты. Они представляют собой соответствующие модам модальные перемещения. Число степеней свободы при экспериментальном анализе конструкции должно быть выбрано с учетом надежного представления форм мод. Именно сложность форм мод в большей степени, чем количество ожидаемых мод, определяет необходимое число степеней свободы [20]. Описанный выше подход к экспериментальному определению динамических характеристик конструкции тем более точен, чем более точно выполняется условие линейности конструкции, т.е. пропорциональности реакции силе возбуждения. Свойство линейности имеет важные для экспериментальной методики следствия [18]: - измеряемые частотные характеристики не зависят от типа и временного закона сигнала возбуждения; возбуждение синусоидальной силой, изменяющейся во времени по синусоидальному закону, с разверткой частоты теоретически дает те же результаты, что и возбуждение широкополосной случайной силой; - измеряемые частотные характеристики не зависят от уровня возбуждения; - измеряемые частотные характеристики не зависят от того, какая точка используется для возбуждения, а какая - для измерения реакции (теорема взаимности Максвелла). Таким образом, экспериментальный модальный анализ сводится к определению элементов матрицы передаточных функций Нij(ω) одним из двух методов: - возбуждение колебаний Fj(ω) в одной из точек объекта j и регистрации отклика (виброускорения или виброскорости) Xi(ω) во множестве точек i; - возбуждение колебаний Fj(ω) во множестве точек объекта j и регистрации отклика (виброускорения или виброскорости) Xi(ω) в одной из точек i. В методе сканирующей лазерной виброметрии, используемом в настоящей работе, реализуется первый из этих подходов. К одной из точек объекта прикладывается воздействие в виде внешней сосредоточенной силы, изменяющейся во времени по гармоническому закону с переменной частотой или по более сложному закону, содержащему сумму гармоник с плотным спектром частот (удар, белый шум). Реакция в виде виброскорости измеряется во множестве точек при их последовательном сканировании [17, 18]. Далее по соотношению (2) определяются компоненты матрицы передаточных функций Нij(ω). 2. Описание экспериментального комплекса модального анализа Установка для экспериментального модального анализа (рис. 2) включает трехкомпонентный сканирующий лазерный виброметр PSV-400-3D, предназначенный для трехмерного анализа колебаний конструкций сложной формы бесконтактным методом. Виброметр представляет собой исследовательский лазерно-цифровой измерительный комплекс, состоящий из трех оптических сканирующих лазерных головок PSV-I-400, модуля сканирования геометрии PSV-А-420, размещенных на штативах либо на стапеле, и системы управления. В систему управления входят три управляющих модуля оптических головок (контроллеры) OFV-5000, блок сопряжения управляющих модулей с промышленным компьютером PSV-E-408. 2 1 Рис. 2. Общий вид экспериментальной установки: 1 - система крепления исследуемой ЛСА на свободном подвесе; 2 - динамик для возбуждения колебаний Fig. 2. General view of experimental setup: 1 - OGVfixingwith unrestrained suspension; 2 - device for vibrations initiating Оптические головки PSV-I-400 представляют собой устройства, в которых находятся гелий-неоновые источники лазерных лучей и приемники отраженного сигнала. Работа каждой оптической головки и связь с компьютером обеспечивается управляющим модулем (контроллером) OFV-5000 и блоком сопряжения. В верхней оптической головке дополнительно встроена видеокамера и лазерный модуль сканирования геометрии. Также для улучшения качества изображения используется портативная цифровая камера со сменными объективами, позволяющая проводить видео триангуляцию для прецизионного сведения лазерных лучей, что повышает качество проводимых экспериментов. Изображение с видеокамеры передается в компьютер и отображается на мониторе. Проведение модального анализа методом трехкомпонентной сканирующей лазерной виброметрии осуществляется с помощью программного комплекса POLYTECPSV9.0. Он позволяет производить анализ собственных форм колебаний, спектральный анализ, быстрое преобразование. 3. Подготовка и проведение экспериментального модального анализа Все лабораторные эксперименты для проведения модального анализа ЛСА из ПКМ методом трехкомпонентной сканирующей виброметрии проводились при закреплении лопатки в свободном подвесе с помощью специальных струбцин и эластичных жгутов (рис. 2) и бесконтактном возбуждении вибраций акустическим динамиком, который позволяет исследовать формы колебаний лопаток в диапазоне частот до 6400 Гц. Динамик располагался со стороны спинки лопатки в её центральной зоне, удаленный от лопатки на расстояние около 100 мм (см. рис. 2). При подготовке лопатки к исследованиям ее поверхность покрывается краской ARDROX 9D1B, обеспечивающей необходимые для лазерной виброметрии оптические свойства поверхности. В проводимых исследованиях в качестве возбуждающего воздействия использован гармонический сигнал с возрастающей частотой (PeriodicChirp) с постоянной амплитудной и возрастающей частотой. Сканирующие лазерные головки размещались на специальном стапеле, на расстоянии 1,5-2 м от исследуемой лопатки. При ориентировании головок выдерживался угол 30° между лазерным лучом и нормалью к поверхности. Построение сетки сканирования проводилось в автоматическом режиме. При этом учитывалось, что большая густота точек сканирования позволяет более точно определять высокие формы, вместе с тем шаг точек не должен быть менее установленного качества выполнения 3-D-выравнивания. Количество точек сканирования для исследованной лопатки в различных экспериментах составляет до 557. Сканирование поверхности вибрирующей лопатки лазерными лучами осуществлялось автоматически. Количество повторений измерений в каждой точке составляло 3-5. Обработка результатов лазерной виброметрии проводилась с применением программного обеспечения PolytecPSV9.0. В результате обработки экспериментальных данных получены амплитудно-частотные характеристики каждой из трех компонент виброскорости в любой точке сетки сканирования и осредненные по объекту в целом, а также анимационные представления форм колебаний. Оценка достоверности результатов и погрешности определения собственных частот выполнялась на основе совместного анализа амплитудно-частотных характеристик и форм колебаний, полученных в различных экспериментах. При этом эксперименты различались частотным диапазоном, разрешением по частоте и параметрами сетки сканирования. Для соответствующих форм колебаний сравнивались собственные частоты, и по коэффициенту вариации оценивался их разброс в различных экспериментах. При наличии значительных расхождений проводилась дополнительная проверка и уточнение собственных форм. 4. Результаты экспериментального определения собственных частот и форм колебаний лопатки спрямляющего аппарата из ПКМ Первичными результатами экспериментального модального анализа лопатки являются амплитудно-частотные характеристики, полученные по результатам каждого отдельного эксперимента. Все результаты экспериментальных исследований усреднялись по совокупности точек сетки сканирования. На рис. 3 в качестве примера приведена экспериментально полученная АЧХ лопатки № 3. В процессе исследований проведено сравнение амплитудно-частотных характеристик по экспериментам, полученным для трех различных лопаток. Значения частот, соответствующих пикам на АЧХ, в разных экспериментах хорошо согласуются между собой. Рис. 3. Амплитудно-частотные характеристики ЛСА Fig. 3. Amplitude-frequency characteristics of OGV Формы колебаний в каждом из экспериментов определены в результате анализа анимационных картин динамического поведения лопатки в процессе колебаний. На рис. 4 представлены первые четыре формы колебаний ЛСА из ПКМ. В диапазоне частот до 6400 Гц у исследованных лопаток выявлено до 22 собственных частот колебаний. Собственные частоты приведены в табл. 1. Рис. 4. Собственные формы колебаний ЛСА из ПКМ, полученные в рамках лабораторных испытаний Fig. 4. Experimental natural frequencies of OGV Для оценки воспроизводимости полученных собственных частот определены коэффициенты вариации (отношение среднего квадратического отклонения к среднему значению) каждой из собственных частот (табл. 1). Рассеяние собственных частот по коэффициенту вариации достигает 0,332 % для седьмой формы колебаний ЛСА, для остальных форм колебаний не превышает 0,25 %. Таким образом, получена хорошая воспроизводимость собственных частот для данной группы экспериментов для различных лопаток. Таблица 1 Собственные частоты для ЛСА Table 1 Naturalfrequenciesoutletguidevane Номер собств. формы Собственная частота, Гц Коэффициент вариации, % Номер эксперимента Среднее значение 1 2 3 1 205,1 204,5 204,5 204,7 0,169 2 277,3 276,9 276,9 277,0 0,083 3 484,3 482,8 483,8 483,6 0,158 4 715,9 716,5 716,3 716,2 0,043 5 736,8 736,7 735,6 736,4 0,090 6 811,2 810,9 812,7 811,6 0,119 7 839,5 834,2 838,3 837,3 0,332 8 941,2 936,7 937,5 938,5 0,25 9 1072,1 1070,2 1070,3 1070,9 0,100 10 1376 1372 1372,2 1373,4 0,164 11 1490,5 1494,9 1494,9 1493,4 0,170 12 - 1760,9 1760,9 1760,9 0,000 13 - 2002,6 2003,1 2002,9 0,018 14 - 2382,2 2382,8 2382,5 0,018 5. Сравнение с результатом численного расчета собственных частот и форм колебаний лопатки спрямляющего аппарата из ПКМ В данном разделе проводится численный модальный анализ с целью определения частот и форм (мод) собственных колебаний ЛСА из ПКМ. Для расчета собственных частот и определения форм колебаний ЛСА использовалась трехмерная модель, форма которой была построена с помощью CAD-системы SolidWorks (рис. 5, а), а структура армирующих слоев задавалась в программном комплексе ANSYSWorkbench. Конструктивно в ЛСА можно выделить перо и полки, посредством которых происходит закрепление лопатки. Для крепления верхней полки к обечайке корпуса используется болтовое соединение. На входной кромке ЛСА предусмотрена выемка, в которую устанавливается титановая накладка. Лопатка состояла из слоев различной конфигурации толщиной 0,21 мм, соответствующих набору технологических выкроек эпоксидного углепрепрега, со схемой армирования [0°/±45°]. Исследования по выбору материала и схемы армирования для изготовления лопатки описаны в работах [2-4]. В модели задавалась полная технологическая схема укладки слоев, каждый слой рассматривался в виде ортотропного однородного материала, главные оси симметрии которого, соответствующие локальной системе координат слоя, ориентированы под углами 0 либо ±45° относительно глобальной системы координат конструкции лопатки. Механические свойства материала слоев, использованные в расчетах, приведены в работах [3, 4]. Математическая постановка решаемой задачи динамической теории упругости в вариационной формулировке заключается в отыскании минимума функционала Лагранжа [21] с учетом сил инерции, распределенных по объему тела. Предполагалось, что задача является линейной, все виды нелинейности - нелинейное поведение материала, контактные граничные условия, конечные перемещения - не рассматривались. Для определения собственных частоти форм собственных колебаний рассматриваемой конструкции принималось, что внешние силы и демпфирование равны нулю. Решение поставленной задачи осуществлялось методом конечных элементов (МКЭ) с использованием многопроцессорного программного комплекса ANSYS Workbench. Для лучшей сходимости решения и снижения погрешностей получаемых результатов генерировалась расчетная сетка МКЭ, ячейки которой имели призматическую форму, с использованием конечного элемента Solid 186. При построении и локальном измельчении сетки не допускались резкие отличия геометрических размеров соседних элементов (более чем в 2 раза). Максимальный размер элемента составлял 2 мм, минимальный - 1 мм. Общее количество конечных элементов составило порядка 500 тысяч. При этом в рамках одного элемента объединялось от четырех до восьми слоев углепластика, ориентированных различным образом. Эффективные упругие свойства материала внутри каждого конечного элемента вычислялись на основе заданных свойств отдельных слоев автоматически процедурой ANSYS Workbench, реализованной в конечном элементе Solid 186. На рис. 5, б, в представлена конечно-элементная сетка для конструкции, использованная в расчетах. По результатам численного анализа были определены формы колебаний и собственные частоты. На рис. 6, а представлены первые четыре формы колебаний ЛСА из ПКМ, полученные в рамках численного расчета. Качественное сравнение результатов численного расчета (рис. 6) и экспериментальных данных (см. рис. 4) показало, что формы колебаний совпадают. а б в Рис. 5. Геометрическая (а) и сеточная (б, в) модель ЛСА Fig. 5. Three-dimensional model of OGV (a) and FE mesh (b, c) Рис. 6. Формы колебаний ЛСА из ПКМ, полученные в рамках численного расчета Fig. 6. Calculated mode shapes of OGV Количественное сравнение собственных частот ЛСА из ПКМ показал, что отклонения расчетных и экспериментальных значений (табл. 2) наблюдаются в диапазоне от 1,08 % (седьмая собственная частота) до 5,64 % (первая собственная частота). При этом наибольшие значения отклонений (5,64, 5,33, 4,33 %) соответствуют первой, третьей и десятой собственным частотам, для остальных частот отклонение не превышает 2,98 %. Таблица 2 Сравнительный анализ собственных частот ЛСА Table 2 Comparative analysis of natural frequencies of OGV Номер собств. формы Собственная частота, Гц Отличие, % Эксперимент Расчет 1 204,7 216,93 5,64 2 277,0 284,15 2,52 3 483,6 505,51 4,33 4 716,2 705,16 1,54 5 837,3 867,01 3,43 6 1070,9 1039,6 2,98 7 1373,4 1388,4 1,08 8 1493,4 1472,5 1,94 9 1760,9 1733,9 1,53 10 2002,9 1894,4 5,42 Таким образом, результаты численного моделирования, выполненного в трехмерной постановке с учетом полной технологической схемы укладки и анизотропии свойств армирующих слоев в конструкции, хорошо согласуются с результатами механических испытаний. Совпадение собственных частот и форм колебаний свидетельствует о том, что жесткость натурной ЛСА при наличии возможных отклонений и разбросов параметров технологического процесса, геометрии выкроек слоев и механических свойств материалов соответствует проектируемой, а лабораторный модальный анализ может использоваться для контроля стабильности геометрических размеров и механических свойств материала слоев. Заключение В результате проведенных работ предложена методика экспериментального исследования собственных частот и форм колебаний полноразмерных композитных лопаток спрямляющего аппарата, базирующаяся на трехкомпонентной сканирующей лазерной виброметрии с использованием программно-аппаратного комплекса PSV-400-3D. Методика позволяет находить собственные частоты и формы колебаний ЛСА из ПКМ в диапазоне частот до 6,4 кГц. На трех натурных образцах ЛСА, полностью соответствующих конструктивным требованиям и изготовленных по одной технологии, проведено 15 экспериментов с различными параметрами сетки, разрешения по частоте и частотного диапазона. В результате экспериментов получены средние значения и коэффициенты вариации собственных частот ЛСА из ПКМ и проиллюстрированы основные собственные формы их колебаний. Рассеяние по коэффициенту вариации экспериментальных значений собственных частот не превышает 1,6 %. Разработана математическая модель ЛСА из ПКМ для проведения численного модального анализа на основе трехмерной компьютерной модели, построенной с помощью CAD-системы SolidWorks. Разработанная численная модель полностью описывает геометрию, схему укладки и анизотропию армирующих слоев, решение задачи модального анализа осуществлялось МКЭ в программном пакете ANSYS Workbench. Результаты численного моделирования хорошо согласуются с экспериментальными данными. Максимальное отличие оценок собственных частот, полученных экспериментальными и численными методами, не превышает 5,64 %. Совпадение расчетных и экспериментальных собственных частот и форм колебаний свидетельствует о том, что жесткость натурной ЛСА при наличии возможных отклонений и разбросов параметров технологического процесса, геометрии выкроек слоев и механических свойств материалов, соответствует проектируемой. В последующих исследованиях разработанная модель позволит оценить изменение АЧХ ЛСА при вариации её схемы армирования и конструктивных особенностей исполнения, а лабораторный экспериментальный модальный анализ может использоваться для контроля стабильности геометрических размеров и механических свойств материала слоев.

About the authors

M A Grinev

OAO «Aviadvigatel»

A N Anoshkin

Perm National Research Polytechnic University

P V Pisarev

Perm National Research Polytechnic University

G S Shipunov

Perm National Research Polytechnic University

M Sh Nikhamkin

Perm National Research Polytechnic University

A A Balakirev

Perm National Research Polytechnic University

I P Konev

Perm National Research Polytechnic University

A Yu Golovkin

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 266

PDF (Russian) - 189