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

Полный текст

Роторы турбомашин являются динамическими системами с конструктивной поворотной симметрией (циклически-симметричные системы (ЦСС)), которые позволяют исследовать их динамику на основе одного сектора. На практике в конструкциях всегда возникают отличия между лопатками (по массе, геометрии, материалу и т. д.) из-за изготовления, износа при эксплуатации и других факторов. Все эти малые отличия лопаток (так называемая расстройка параметров) нарушают циклическую симметрию. Причины расстройки параметров (неидентичность секторов) ЦСС могут быть разные. Чаще всего расстройка лопаток является маленькой величиной (т. е. собственная частота лопаток различна в нескольких процентах от номинального значения), но она может привести к возникновению резонансных напряжений, которые выше оптимальных на 20 % и более. Реальные увеличения максимальных напряжений на отдельных лопатках составляют от 0 до 20 % и зависят от размещения лопаток на диске. Можно предположить, что исследование расстройки будет играть большую роль при проектировании новых турбомашин, а также в правильной оценке и увеличении ресурса уже работающих [ 1; 2]. МУП используется, чтобы исследовать характеристики динамики ЦСС на основе одного сектора (рис. 1, б). В данной работе используется способ МУП, детально описанный в работах [3; 4], для исследования влияния расстройки лопаток на характеристики колебаний рабочих колес турбомашин, а именно: собственные частоты колебаний, формы собственных колебаний и максимум отклика вынужденных колебаний. Расчеты МУП проводятся МКЭ 60 Вестник СибГАУ. № 1(53). 2014 с помощью треугольных конечных элементов переменной толщины STI218, которые описаны в работах [5; 6]. Результаты данного исследования для рабочего колеса сравнены с экспериментальными данными и с результатами, которые получены с использованием МКЭ для ПМ (рис. 1, а) с использованием программного комплекса ANSYS. Погрешность результатов между методами МУП и ПМ и результатов метода МУП и результатов, полученных с помощью программы ANSYS, определены как £i = у-расс _ у-расс /муп - /пм /i расс ПМ .100 £2 = расс расс /МУП - J ANSYS расс /ANSYS .100 (1) Рис. 1. Модель конечных элементов рабочих колес: а - полный диск; б - один типичный сектор Общая схема моделирования уменьшенного порядка. Предположим, что степени свободы упорядочены так, чтобы получить блочно-диагональные формы при сборке матриц масс и жесткости всей конструкции. Они представлены в виде [4] M = K = Г Md 0 " ГI ® Md 0 " 0 Mb 0 I ® Mb Г Kd 0 " I ® K 0 " 0 Kb i 0 I ® Kb _ (4) где /МУП - собственная частота колебаний рабочих колес с расстройкой, полученная с использованием МУП; /ПМс, /ANsYS - собственная частота колебаний рабочих колес с расстройкой, полученная с использования МКЭ для полной модели и программы ANSYS. Реализация расстройки. Расстройка лопаток моделируется путем изменения эквивалентных жесткостей Л или изменения собственных частот лопаток, когда они стационарны (собственные частоты консольных лопаток). Эквивалентная расстроенная жесткость k-й формы n-й лопатки представлены в виде [7] Лb,n =p,n )2 =(1 + A/nk )лb , (2) где Ль - жесткость k-й настроенной формы колебаний консольной лопатки и А/П - расстройка параметров k-й формы n-й лопатки. В данной работе рассмотрим расстройку жесткости лопаток путем изменения модуля упругости лопаток. Тогда модуль упругости i-й лопатки определен как [3] Et =(1 + A/i) E0, (3) где E0 - модуль упругости лопатки без расстройки, А/ - параметр расстройки i-й лопатки. где I - единичная матрица размера N; N - число секторов; Md, Kd - матрицы масс и жесткости одного сектора диска; Mb, Kb - матрицы масс и жесткости одной лопатки; символ ® обозначает оператор Кро-некера, описанный в работе [4]. Перемещения одного сектора (рис. 1, б) разделены на две части: перемещение диска с безмассовыми лопатками, которое имеет циклическое свойство (часть форм колебаний диска с лопаткой принадлежат степеням свободы лопатки и обозначаются 5d, а часть диска - Sd), перемещение лопатки, закрепленной жестко с диском 5b. Перемещения всей структуры определены как [4] "Sd 1 Г0 1 (5) где а - вектор обобщенных координат диска для всех форм, а = [ад а1T ... aTP]" ; ап - один вектор обобщенных координат, соответствующих n узловым диаметрам диска; P - максимальное число гармоник или максимальное число узловых диаметров; b - вектор обобщенных координат для всех N лопаток, [т т T T b0 b1 ... bNJ ; bi - вектор обобщенных координат i-й лопатки. Уравнение движения системы «диск-лопатки» в виде матрицы определено как [4] "Sd" "0 " X = а + 43 MD | 5b _ Id + 5dTMb5d 5 dTMbbb . 5bTMb5d Ib +(1 + Yj) 0 0 Ki diag "2^k | k=1,...,mb 5b Kb5d b а g S3 ^ ^ j 05 it n di<2g (1 + A/nk ) k=1,..., mb 1 L ^ b 5dTF 5b F (6) где Id = S MdS , Kd = S KdS - соответственно, матрицы масс и жесткости диска, когда все формы + а б 61 Математика, механика, информатика с матрицей масс нормированы; Ib = 5 Mb5 , Kb = 5b Kb5b - соответственно, матрицы масс и жесткости N лопаток, когда все формы с матрицей масс нормированы; N - количество лопаток, Е,k - коэффициент вязкого демпфирования; у - коэффициент демпфирования структуры; F - соответственно, вектор возбуждающих силы степеней свобод всех лопаток при сборке. Все матрицы в уравнении (6) детально описаны в работах [3; 4]. Характеристики колебаний рабочих колес турбомашин. Собственные частоты колебаний. Метод МУП определяет собственные частоты колебаний рабочих колес без расстройки в зависимости от числа узловых диаметров и узловых окружностей. Число узловых диаметров одной формы колебаний нумеруется так, чтобы соответствовать одной фазе между смежными лопатками, определенной как [4] 2пп Фп =- (П = 0 ..., P) (7) 5- Nb Z 52,i j=1 Nb z 5: ,j=1 2 j, n (i = 1, ..., N), (8) ?jn (n = 1, ..., N), (9) где 5n - евклидова норма перемещений лопаток; Nb -число степеней одной лопатки; 5 j п - перемещение j-й степени n-й лопатки. Определение возбуждающей силы. Возбуждающие силы - это гармонические функции от времени, отличающиеся фазой между секторами. Фаза i-й лопатки определена как [4] 2nC(i -1) N (i = 1, ..., N) (10) где п - число узловых диаметров; N - количество лопаток рабочего колеса. Формы собственных колебаний. В данной работе используем евклидову норму перемещений, чтобы исследовать формы колебаний рабочих колес с расстройкой и без расстройки. Для рабочих колес без расстройки евклидова норма определена в зависимости от числа узловых диаметров и узловых окружностей и является гармонической функцией при присутствии узловых диаметров (см. рис. 6, а, в). Евклидова норма перемещений для собственных колебаний выражена скалярной величиной, называющейся формой относительных перемещений лопаток 5г- . Эта норма перемещений для i-й лопатки определена как [4] где C - порядок энергии возбуждений. Вектор внешней возбуждающей силы представлен в виде [4] F = [/е}Ф1 /г]Ф2 ... fe^N }" =4Nec+1 ® f , (11) где eC+1 - (С+1)-й столбец матрицы Фурье; f - вектор возбуждающей силы одной лопатки. На физической модели (рис. 2) математическое уравнение распределения вектора газодинамических сил {f} по длине лопатки может быть описано следующим видом [5]: [f } = [0 Py Рх mz 0 0]", (12) где Py = Py (z, t) и Рх = Рх (z, t) - возбуждающие силы по направлениям х и y; mz = mz (z, t) - крутящий момент (рис. 2). Определение Px, Py, mz детально описано в работе [5]. где 5 j t - перемещение j-й степени свободы i-й лопатки; Nb - число степеней свободы одной лопатки; N -число лопаток в системе. Максимум отклика вынужденных колебаний. В данной работе используем евклидову норму перемещений, чтобы исследовать максимальную амплитуду колебаний лопаток системы «диск-лопатки» с расстройкой и без расстройки. Для настроенной системы пик нормы перемещений всех лопаток одинаков, но это не так для расстроенной системы. Евклидова норма перемещений для вынужденных колебаний является скалярным выражением, оценивающим отклонение амплитуды лопаток, и определена как [4] Рис. 2. Физическая модель лопатки ротора турбины под действием газодинамических сил Результаты исследования. В качестве примера рассмотрим модельное рабочее колесо, содержащее 24 лопатки. Геометрические размеры и характеристики материала рассчитываемой конструкции: внутренний радиус 0,0135 м, внешний радиус 0,06 м, толщина диска и лопатки 0,002 м, длина лопатки 0,036 м, ширина лопатки 0,012 м, модуль упругости материала 210 ГПа, плотность 7850 кг/м3, коэффициент Пуассона 0,3. Возбуждающие силы приложены на все узлы пика каждой лопатки: [f }=[0 1 1 1 0 0]", коэффициент вязкого демпфирования ^ = 0, коэффициент демпфирования структуры у = 0,6 %. Параметры расстройки приведены в табл. 1. 62 Вестник СибГАУ. № 1(53). 2014 Конечно-элементная модель сектора «диск-лопатки» для МУП и системы «диск-лопатки» при использовании МКЭ для ПМ на основе треугольных конечных элементов STI218, соответственно, содержит 174 и 3 312 степени свободы. Результаты расчета собственных частот системы «диск-лопатки» без расстройки при использовании МУП, ПМ и программного комплекса ANSYS приведены в табл. 2 и на рис. 3. Подобные результаты исследований колебаний рабочего колеса с расстройкой параметров по собственным частотам, погрешностям результатов и евклидовым нормам перемещений лопаток методов МУП, ПМ и программы ANSYS приведены на рис. 4-6. Максимальная амплитуда перемещений лопаток, соответствующая порядкам энергии возмущения C = 5, приведена на рис. 7. Параметры расстройки Таблица 1 Лопатка 1 2 3 4 5 6 7 8 9 10 11 12 f (%) -0,52 -1,87 -1,82 -0,39 -5,01 -0,85 1,42 7,62 2,93 2,72 2,77 -4,92 Лопатка 13 14 15 16 17 18 19 20 21 22 23 24 f (%) -8,07 -4,91 5,93 -6,92 -0,41 0,43 1,84 -5,47 2,39 3,81 4,11 3,92 Сопоставление расчетных и экспериментальных значений частот собственных колебаний модельного диска [8] Таблица 2 Формы 1 2 3 4 5 6 7 8 9 10 11 12 Эксп. [1] 210 265 340 501 681 803 922 938 961 1008 1027 1030 МУП 255,1 260,1 321 492 669 808 911,9 987 1041 1079 1104 1119 ПМ 255,1 260,1 321 492 669 808 911,9 987 1041 1079 1104 1119 ANSYS 265 284 318 491 655 778 865 929 973 1006 1026 1039 Формы 13 14 15 16 17 18 19 20 21 22 23 24 Эксп. [1] 1032 1386 1362 1723 2109 2714 3452 4102 4738 5112 5513 5983 МУП 1123 1461 1521 1737 2147 2713 3354 4002 4613 5155 5605 5947 ПМ 1123 1461 1521 1737 2147 2713 3354 4002 4613 5155 5605 5947 ANSYS 1043 1481 1524 1743 2208 2815 3447 4040 4567 5154 5655 5700 Поворо тные зоны частот -- -<+-1 -«-< <►-0-“ -♦-i-*- I * * Число узловых диаметров п & 3000 £ 2000 о U 1000 -е- МУ п ....... 1М - ANSYS Число d 30 40 рм колебаний Рис. 3. График частот собственных колебаний «диск-лопатки» без расстройки в зависимости от числа узловых диаметров n и узловых окружностей m Рис. 4. График собственных частот рабочего колеса с расстройкой параметров 63 Математика, механика, информатика а б Рис. 5. Погрешность результатов собственных частот: а - при использовании МУП и ПМ; б - при использовании МУП и программы ANSYS О S 10 1S 20 25 Чпсло лопаток Число лопаток Рис. 6. График евклидовой нормы перемещений лопаток, соответствующей собственным расстроенным частотам 64 Вестник СибГАУ. № 1(53). 2014 а б Рис. 7. График максимальных значений отклика вынужденных колебаний рабочих колес при С = 5 В работе представлен способ МУП для исследования колебания рабочих колес турбомашин без расстройки и с расстройкой жесткости лопаток. В этом способе использован подход составляющей формы: один сектор системы «диск-лопатки» разделен на две части. Это диск с безмассовыми лопатками, жестко защемленный на внутреннем радиусе, и консольная лопатка, жестко защемленная на поверхности между диском и лопаткой [9]. При этом используются матрицы масс и жесткости конечных элементов одиночной лопатки. Поэтому расстройка жесткости лопатки легко моделируется в расчете. Дополнительно при использовании МУП число степеней свободы расчета значительно уменьшается (из 3 312 степеней свободы при использовании ПМ число степеней свободы уменьшено до 174 для МУП), что снижает трудоемкость и численные затраты времени на ЭВМ. Из рис. 3, 6 и 7 видно, что на графике собственных частот рабочих колес без расстройки в зависимости от числа узловых диаметров и окружностей часто появляются поворотные зоны частот. В этих зонах амплитуды собственных колебаний и максимальные амплитуды вынужденных колебаний рабочих колес даже с маленькой расстройкой жесткости сильно изменяются как качественно, так и количественно (рис. 6, г; 7, б). Из табл. 2 и рис. 5 видно, что численные результаты МУП для рабочих колес без расстройки и с расстройкой МКЭ на основе треугольных конечных элементов STI218 хорошо совпадают с экспериментальными данными и результатами программного комплекса ANSYS. Полученные результаты 2-х этих методов очень близки. Однако использование треугольных конечных элементов STI218 позволяет учитывать модели со сложной геометрией сечения. Здесь представлено влияние расстройки жесткости лопаток на характеристику колебаний рабочих колес турбомашин. Эти результаты совпадают с выводами работ, которые опубликованы авторами Gisli Ottarsson [3], J. R. Bladh [4], и позволяют решать задачи чувственности колебаний и повышения прочности рабочих колес с расстройкой при проектиро вании новых и эксплуатации старых конструкций турбомашин. Данные исследования являются базовыми для дальнейшего анализа влияния расстройки параметров на долговечность рабочих колес турбомашин.

Об авторах

Олег Владимирович Репецкий

Восточно-Сибирский институт экономики и права

Email: repetskiy@vsiep.ru
доктор технических наук, профессор, президент

- До Мань Тунг

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

Email: manhtungcvp@yahoo.com
аспирант

Список литературы

Дополнительные файлы