close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, СИСТЕМНЫЙ АНАЛИЗ
УДК 51-74:621
ББК 22.311
О.В. РЕПЕЦКИЙ
проректор по международной деятельности
Байкальского государственного университета экономики и права,
доктор технических наук, профессор, г. Иркутск
e-mail: repetskiy@isea.ru
ФАН ВАН ТУАН
аспирант Байкальского государственного университета
экономики и права, г. Иркутск
e-mail: yeubeconlam@yahoo.com
ПОСТРОЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
ДЛЯ АНАЛИЗА ВЛИЯНИЯ ФРИКЦИОННЫХ ДЕМПФЕРОВ
НА КОЛЕБАНИЯ ЛОПАТОК ГАЗОТУРБИННЫХ ДВИГАТЕЛЕЙ
Проведено построение стандартных фрикционных демпферных элементов для
анализа колебаний лопаток турбин с фрикционным демпфером. Выполнен расчет
для двух видов демпферов ? круглого и трапециевидного ? с применением
метода прямого численного интегрирования.
Ключевые слова: математическая модель, трение, колебание, демпферы.
O.V. REPETSKIY
Vice Rector for International Relations, Doctor of Science in Engineering,
Professor, Baikal State University of Economics and Law, Irkutsk
e-mail: repetskiy@isea.ru
PHAN VAN TUAN
post-graduate student,
Baikal State University of Economics and Law, Irkutsk
e-mail: yeubeconlam @yahoo.com
CONSTRUCTION OF MATHEMATICAL MODEL
FOR ANALYSIS OF FRICTION DAMPERS INFLUENCES
ON VIBRATION OF GAS TURBINE ENGINES? BLADES
The authors construct standard friction damper elements for the vibration analysis
of turbine blades with friction damper and make the calculation for two types of dampers ? circle and trapezoidal ? with the use of direct numerical integration method.
Keywords: mathematical model, friction, vibration, damper.
ных элементов (КЭ). Тогда в каждый момент
необходимо решить две задачи: статическую
контактную задачу в зоне контакта между лопаткой и ФД и динамическую задачу
системы. Результаты решения контактной
задачи являются входными параметрами динамической задачи. Однако эта модель требует очень большого количества вычислений,
причем возникает несходимость численных
решений и можно говорить, что это неэффективная модель. Поэтому можно построить более простую модель решения данной
задачи ? на основе стандартных элементов
Для устранения колебаний лопаток существует большое количество способов. Одним
из них является использование фрикционных
демпферов (ФД). Один из принципов работы ФД показан на рис. 1. Под действием
центробежной силы FД демпферы давят на
полки, создают силу трения и приводят к
уменьшению колебаний лопаток.
Самое обобщенное решение поставленной задачи может быть получено путем исследования системы лопатки?ФД как единой
механической системы при моделировании
конструкции с помощью трехмерных конеч© ?.?. ????????, ??? ??? ????, 2011
200
О.В. РЕПЕЦКИЙ, ФАН ВАН ТУАН
фрикционных демпферов (ФДЭ). Требования для построения стандартных ФДЭ:
? точное отображение действия ФД на
колебание лопаток;
? обеспечение простоты решения динамических уравнений системы.
= [xД yД zД]T; перемещения контактной
точки на полке лопатки {uП} = [xП yП zП]T;
относительную скорость между ФД и полкой в момент t:
?x? ? ? x? ? ?
?
?
{u??? }= {u? ? }? {u? ? }= ?y? ? ? y? ? ?;
?
?
?z? ? ? z? ? ?
относительные перемещения между ФД и
полкой в момент t {?u} = {u??? }dt ; действующие силы от ФД на полку {FД} = [Tx Ty N]T.
Работа ФДЭ включает три режима:
1. Режим залипания. В этом режиме ФД
и полка являются одним целым ? их относительная скорость равна нулю:
?x? ? ? x? ? ? ?0 ?
?
?
{u??? }= ?y? ? ? y? ? ? = ??0 ??
?
? ? ?
?z? ? ? z? ? ? ?0 ?
????????????
???????
???????????
???????
?????
???????????
Рис. 1. Структура лопаток
с фрикционными демпферами
и виды демпферов
и zД = zП, а действующие силы от ФД на
полку определяются выражением
?Tx ?
{F? }= ??Ty ?? = ? ??K? ?? .{u? }+ {P}=
? ?
?N ?
??
?? модель
На рис. 2а показана трехмерная
ФДЭ.
Здесь uД ? степени свободы демпфе?
ра, uП ? степени свободы полки. Согласно
этой модели, ФД характеризуется следующими параметрами: матрица жесткости [KД],
постоянная сила давления {РД} = [0 0 P]T,
коэффициент трения µ.
?kxx kxy kxz ?
?
?
??K? ?? = ?kyx kyy kyz ? .
?
?
??kzx kzy kzz ??
Контактные параметры включают: перемещения контактной точки на ФД {uД} =
?kxx .x ? + kxy .y ? + kxz .z ? ? ?0 ?
?
? ? ?
= ? ?kyx .x ? + kyy .y ? + kyz .z ? ? + ?0 ?.
?
? ?P ?
?kzx .x ? + kzy .y ? + kzz .z ? ? ? ?
2. Режим скольжения, когда ФД и полка
скользят относительно друг друга. Динамические уравнения этого режима:
??Tx ??
µ N ??x? ? ? x? ? ??
z ? = z ? , ? ? = ? xyz ?
?,
u??? ??y? ? ? y? ? ??
??Ty ??
xy
= (x? ? ? x? ? )2 + (y? ? ? y? ? )2 .
где u???
?
??????????
???
?????
?
??
??????? ??
?
??
?
?
??
?
?
?
??
?
?????????
??
??
??
?
?????
??
??????????
?
?
??
?
????????? µ? ??
?
?µ?
?
Рис. 2. Модели ФДЭ: а ? трехмерная; б ? одномерная
201
???????? ????. 2011. ? 1 (75)
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, СИСТЕМНЫЙ АНАЛИЗ
3. Режим открытия. ФД и полка не контактируют. Динамические уравнения имеют
вид xД = 0, yД = 0, zД = P / kzz, Tx = 0,
Ty = 0, Nz = 0.
Во время работы ФДЭ часто имеет место переход от одного режима к другому.
Значения матрицы жесткости [KД] и коэффициента трения ? зависят от характеристик
материалов, геометрии ФД и полки в зоне
контакта. Определение значений этих параметров показано в [1?4].
Трехмерная модель является обобщенной моделью, изображающей действие ФД,
и корректна для трехмерных задач. Однако
для отдельных моделей ФД лучше использовать одномерную модель, чтобы уменьшить
затраты компьютерного времени. Одномерная модель ФДЭ изображена на рис. 2б.
Эта модель характеризуется параметрами:
жесткость KД, постоянная сила давления N и
коэффициент трения µ.
Режимы работы ФДЭ показаны в табл. 1.
Таблица 1
Режимы работы ФДЭ
Режим залипания
Режим скольжения
?T = ?k.x ?
?
? T ? µ.N
??
?x ? = x? ?
?T = ?sign(x? ? ).µ.N
?
µ.N
?
?x ? = ?sign(x? ? )
k
?
??x? ? = x? ?
Таблица 2
Моменты перехода
между режимами работы ФДЭ
Залипание к скольжению Скольжение к залипанию
µN
µN
?
?
?x ? =
?x ? = ?
k или ?
k
?
??x? ? > 0
??x? ? < 0
µN
µN
?
?
?x ? = k
?x ? = ? k
?
?
или ?x? ? = 0
?x? ? = 0
?x?? < 0
?x?? > 0
? ?
? ?
??
??
На практике мы часто встречаем ситуацию,
в которой контакты происходят не по точкам,
а по плоскостям или линиям. Тогда, согласно
[4], контакты по плоскостям или линиям моделируются контактами дискретных точек и в
каждой точке моделируются одним ФДЭ.
С данной моделью ФД влияние ФД на
колебание лопатки осуществляется с учетом двух факторов: изменение собственных
частот механической системы лопатки?ФД
(в режиме залипания); изменение колебаний механической системы лопатки?ФД под
действием внешних сил.
Динамические уравнения лопаток в МКЭ
имеют вид
[M]{???} + [C]{??} + [K]{? } = {F (t )} + {F? }, (1)
где [M], [C], [K] ? глобальные матрицы
масс, вязкого демпфирования и жесткости;
{???}, {??}, {? } ? глобальные векторы узловых ускорений, скоростей и перемещений;
{F(t)} ? глобальный вектор внешней динамической нагрузки; {FД} ? глобальный вектор
действующих сил на полку от ФД.
Моменты перехода между режимами
определены соотношениями, приведенными
в табл. 2.
?
? ?
????????????????
?
?
Рис. 3. Моделирование контактов дискретными ФДЭ:
а ? круглый демпфер; б ? трапециевидный демпфер
???????? ????. 2011. ? 1 (75)
202
О.В. РЕПЕЦКИЙ, ФАН ВАН ТУАН
Уравнение (1) является нелинейным
уравнением. Оно не имеет аналитического
решения и может быть решено только с
использованием численных методов, таких,
например, как метод прямого численного
интегрирования (Ньюмарка, Вилсона, Рунге?Кутта и др.), метод гармонического
баланса и др. Мы выбрали метод прямого
численного интегрирования Ньюмарка и
пакет программ ANSYS. Расчет выполняется для двух видов демпферов: круглого
(с 6 ФДЭ) и трапециевидного (с 18 ФДЭ)
(рис. 3). Материал лопатки: модуль Юнга
E ? 2,01.e5 МПа; коэффициент Пуассона
µ ? 0,3; плотность ? ? 7 700 кг/м3; жес-
ткость ФДЭ KД = 5.e7 Н/м; сила давления
P = 6 Н (согласно [4]).
Первый шаг: анализ влияния ФД на свободные колебания лопаток. Модель анализа
показана на рис. 4, результаты ? на рис. 5.
Второй шаг: анализ влияния ФД на вынужденные колебания лопаток. Выражение
вынужденной силы: Fx = 0,1sin(?o .t ), где
?o ~ 515 рад/с ? угловая собственная частота лопатки, при которой происходит резонанс. На рис. 6а показана схема действия
вынужденной силы, на рис. 6б ? смещение
вершины лопатки при резонансе без ФД. На
рис. 7 показаны вынужденные колебания лопаток с трапециевидным ФД при резонансе.
??
??
?
?
??
? ??? ?
????? ??????????????
?? ?
?
?
?
~ ???????
~
???
?????
?
??????
?
???
??
??
??
??
??? ?
?
????? ?
??
????
?
?
?
Рис. 4. Модель анализа влияния ФД на свободные колебания лопаток:
а ? размеры модели; б ? модель в ANSYS; в ? схема действующей силы
???
???
???
???
???????????
???
???
???
???
???
???
?
?
????
????
????
????
????
????
????
????
????
????
???????????????????????????????????????????????????????????????????
?????
???????????
???????????????????????????????????????????????????????????????????
?????
?
Рис. 5. Смещение вершины лопатки при свободных колебаниях с ФД:
а ? с круглым демпфером; б ? с трапециевидным демпфером, м
203
?
???????? ????. 2011. ? 1 (75)
?
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, СИСТЕМНЫЙ АНАЛИЗ
?????
?????
?????
?????
?
?
?
?
?????
?
???????????
?
?
??????
??????
??
??
??????
??????
??????
?????????????????????????????????????????????????????????????????????
?????
??
??
??
?????????????????????????????????????????????????????????????????????
?????
?
?
Рис. 6. Вынужденные колебания лопаток без ФД при резонансе:
а ? схема действия вынужденной силы Fx, Н; б ? смещение вершины лопатки, м
???????????
?
?
?
?
???????????
???
???
???
???
?
???
?
???
??
??
?
????
??
??
??
????
????
????
?????????????????????????????????????????????????????????????????????
?????
?
?????????????????????????????????????????????????????????????????????
?????
?
Рис. 7. Вынужденные колебания лопаток с трапециевидным ФД при резонансе:
а ? сила трения одного ФДЭ, Н; б ? смещение вершины лопатки, м
???
???
???????????
???
???
???
???
???
???
???
???
???
???
?
????
?
????
???????????
????
????
????
????
????
????
?????????????????????????????????????????????????????????????????????
?????????????????????????????????????????????????????????????????????
?????
?????
?
?
Рис. 8. Смещение вершины лопатки при вынужденных колебаниях лопаток
с ФД при отсутствии резонанса: а ? с круглым демпфером;
б ? с трапециевидным демпфером, м
???????? ????. 2011. ? 1 (75)
204
О.В. РЕПЕЦКИЙ, ФАН ВАН ТУАН
Следующим шагом выполнен расчет
вынужденной силы Fx = 0,1sin(?o .t ), где
?o = 300 рад/с, при отсутствии резонанса.
На рис. 8а показано смещение вершины лопатки с круглым ФД, на рис. 8б ? с трапециевидным ФД.
Мы можем сделать следующие выводы:
? использование стандартных ФДЭ и
дискретизаций упрощает расчет колебаний
лопаток с ФД;
? трехмерный ФДЭ универсальнее одномерного ФДЭ, но в большинстве случаев достаточно использования одномерного ФДЭ;
? метод прямого численного интегрирования является эффективным методом решения систем нелинейных дифференциальных
уравнений;
? при одинаковых значениях центробежной силы затухание колебаний при использовании круглого демпфера наступает быстрее,
чем при использовании трапециевидного.
Список использованной литературы
1. Иванов А.С. Контактная жесткость неподвижных соединений деталей машин: дис. ... д-ра техн. наук.
М., 2006.
2. Матвеев В.В. Исследование демпфирующей способности попарно бандажированных турбинных лопаток
в зависимости от условий сопряжения их бандажных полок // Проблемы прочности. 1978. № 8.
3. Пыхалов А.А. Контактная задача статического и динамического анализа сборных роторов турбомашин:
дис. ... д-ра техн. наук. М., 2006.
4. Dipl.-Ing. Lars Panning. Auslegung von Reiblementen zur Schwingungs-dampfung von Turbinenschaufeln /
Institut fur Dynamik und Schwingungen. Universitat Hannover, 2005.
Bibliography (transliterated)
1. Ivanov A.S. Kontaktnaya zhestkost? nepodvizhnykh soedinenii detalei mashin: dis. ... d-ra tekhn. nauk. M.,
2006.
2. Matveev V.V. Issledovanie dempfiruiushchei sposobnosti poparno bandazhirovannykh turbinnykh lopatok v
zavisimosti ot uslovii sopryazheniya ikh bandazhnykh polok // Problemy prochnosti. 1978. № 8.
3. Pykhalov A.A. Kontaktnaya zadacha staticheskogo i dinamicheskogo analiza sbornykh rotorov turbomashin:
dis. ... d-ra tekhn. nauk. M., 2006.
4. Dipl.-Ing. Lars Panning. Auslegung von Reiblementen zur Schwingungs-dampfung von Turbinenschaufeln /
Institut fur Dynamik und Schwingungen. Universitat Hannover, 2005.
205
???????? ????. 2011. ? 1 (75)
1/--страниц
Пожаловаться на содержимое документа