ISSN 1992-6502 (P ri nt)_
2016. Т. 20, № 4 (74). С. 3-8
Ъьомшь QjrAQnQj
ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru
УДК 539.3
Термонапряженное состояние
центробежного колеса компрессора при наличии трещины в. с. Жернаков 1, и. н. Будилов 2
1 zhvs@mail.rb.ru, 2 budilov02@mail.ru
ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)
Поступила в редакцию 05.12.2016
Аннотация. Выполнен трехмерный упругий расчет напряженно-деформированного состояния (НДС) центробежного компрессора (ЦК) с учетом эксплуатационных нагрузок и температур. Произведена оценка влияния газодинамических и центробежных сил на НДС пера лопатки. Выявлены наиболее нагруженные участки ЦК, построены векторные поля главных напряжений. Получены качественные и количественные соотношения для КИН и ./-интеграла по фронту трещины с учетом ее размеров и формы фронта.
Ключевые слова: центробежное колесо; компрессор; титановый сплав; прочность; напряженно-деформированное состояние, коэффициент интенсивности напряжений, интеграл Черепанова-Райса.
Одним из наиболее напряженных частей вертолетных ГТД является центробежный компрессор (крыльчатка турбины), лопатки которой испытывают наибольшие центробежные и изгибающие нагрузки [1, 2]. В зонах сопряжения лопатки ЦК с основанием диска максимальные напряжения могут достигать предела текучести и переходить в область пластических деформаций. В этих высоконагруженных областях возникают и распространяются трещины, вызывающие усталостное разрушение и отказ турбины ЦК, что в свою очередь приводит к необходимости более точного расчета термонапряженного состояния лопатки ЦК и оценке силовых параметров разрушения в случае зарождения трещин.
ФОРМИРОВАНИЕ РАСЧЕТНОЙ МОДЕЛИ
Центробежный компрессор представляет собой монолитное тело, условно состоящее из основания диска, длинной и короткой лопаток в количестве 36 шт. и вала. На основании чертежей была разработана твердотельная модель ЦК с учетом ее симметрии в пакете SoHdWorks, показанная на рис. 1.
С помощью специально выбранной меридиональной поверхности, повторяющей поверхность спинки лопатки, было произведено отсечение фрагмента от исходной модели.
Рис. 1. Тельная модель центробежного компрессора
Поскольку модель обладает осевой симметрией, разумно с точки зрения оптимизации конечно-элементного рас-чета и использования более мелкой сетки, рассматривать только фрагмент исходной модели, включающий в себя только одну основную и одну короткую лопатку. Один из вариантов сетки КЭ для полноразмерного объекта ЦК приведен на рис. 2.
Рис. 2. Конечно-элементная сетка
Данная геометрическая модель позволила использовать в процессе расчетов более мелкую сетку конечных элементов для получения достоверных результатов расчета.
Центробежный компрессор показан на рис. 3.
Рис. 3. Центробежный компрессор
Все расчеты были выполнены в статической постановке в пакете ANSYS/Workbench [3, 4].
Расчетные параметры конечно-элементной модели. При выполнении расчетов использовались предварительно полученные данные о распределении давления и температуры по перу ко-
роткой и длинной лопаток. Для возможности учета центробежной нагрузки частота вращения была принята равной 44000 об/мин.
Материал диска - титановый сплав с механическими характеристиками: при температуре Т = 20 °С - модуль упругости при растяжении Е=1,1*105 МПа, коэффициент Пуассона V =0,35; при температуре Т = 300 °С -Е = 1,0* 105 МПа, V = 0,37. Механические свойства в промежутке температур от 20 до 300 °С определялись линейной аппроксимацией.
Давление на перо лопатки определялось при решении газодинамической задачи [5] в пакете CFX, а затем поля давлений передавались в расчетный модуль ANSYS/Workbench, тем самым решалась связанная мультифизическая задача с учетом всех сил, давлений и полей температур, действующих на центробежный компрессор. При решении термопрочностной задачи были определены поля температур, представленные в секторе диска в виде изополос, показанных на рис. 4, а, б.
Рис. 4. Поля температур в секторе характерного элемента (а) и секторе с фрагментом вала (б) центробежного компрессора
Граничные условия. Расчет полноразмерной модели предполагал ограничения переме-
а
B. C. Жернаков, И. Н. Будилов • Термонапряженное состояние центробежного колеса.,
5
щений в осевом направлении по внутренней поверхности вала. В радиальном и тангенциальном направлении поверхность вала была задана свободной командой Cylindrical Support.
При расчете сектора модели, наряду с ограничениями осевых перемещений в радиальном направлении по внутренней поверхности вала, задавались обобщенные условия симметрии по боковым поверхностям сектора, что позволяло учитывать угловую скорость, как и на полноразмерной модели. Модель трещины. Область с эллиптической трещиной описывалась изопараметрическими элементами второго порядка и задавалась отдельно от основной модели с помощью инструментария CRACK в модуле Static Stractural. В результате образовалась гибридная конечно-элементная модель ЦК с эллиптической трещиной. ЦК с трещиной с радиально-кольцевой структурой КЭ показан на рис. 5 с различным увеличением. Трещина моделировалась в форме эллипса с параметрами 2с=2,8 мм, а=0,4 мм, располагалась в зоне радиусного перехода от длинной лопатки к основанию диска ЦК. Ориентацию трещины выбирали т.о., чтобы ось Х располагалась вдоль короткого радиуса эллипса, а ось Z-вдоль длинной стороны фронта трещины.
#
- гг
щ
Рис. 5. Гибридная сетка КЭ в секторе ЦК с наведенной трещиной
РЕЗУЛЬТАТЫ РАСЧЕТОВ
Результаты выполненного расчета МКЭ по оценке силовых и энергетических параметров разрушения позволили выявить качественные и количественные закономерности с учетом размеров трещины и ее формы фронта. Размеры трещины варьировались 2С=1,8; 2,8; 3,2 мм, а=0,4; 0,8 мм.
На рис. 6 показаны поля эквивалентных напряжений в области наведенной трещины. Видно, что в зоне трещины наблюдается всплеск напряжений. Существенный рост напряжений, также проявляется в области всего фронта трещины, особенно в зонах выхода трещины на поверхность и в точке где трещины достигает максимальной глубины.
Рис. 6. КЭ сетка и напряжения по Мизесу в области вершины трещины, расположенной в основании длинной лопатки ЦК
На рис. 7 показаны действующие эквивалентные напряжения при одновременном действии газодинамических и инерционных сил.
Расчеты силовых параметров разрушения приведены на графиках рис. 8-10, здесь представлены зависимости КИН первого (К1), второго (К2), и третьего (Кз) рода по фронту трещи-
ны.
0.000
0.060 (т)
0.030
Л
б
Рис. 7. Напряжения по Мизесу в длинной лопатке при одновременном действии газодинамических и инерционных сил:
а - колесо; б - лопатки
Рис. 8. Значения силовых параметров разрушения К1, К2, Кз:
2С=1,8 мм; а=0,4 мм
Установлено, что наибольшее значение коэффициент интенсивности напряжений первого рода К1 достигает на максимальном удалении по фронту трещины, в то время как коэффициент интенсивности напряжений третьего рода К3 достигает наибольших значений при выходе фронта трещины на поверхность пера лопатки.
а
В. С. Жернаков, И. Н. Будилов • Термонапряженное состояние центробежного колеса.,
7
Получено, что для данного вида нагружения лопатки ЦК коэффициент интенсивности напряжений второго рода К2 по абсолютному значению меньше, чем КИН первого (К1) более чем в 3 раза.
Ъ
и е
«
¿2 714
11
|тш|
Рис. 9. Значения силовых параметров разрушения К1, К2, К3:
2С=3,2 мм; а=0,4 мм
Установлен характер поведения кривых силовых параметров разрушения от формы фронта трещины. Так, рост длины трещины при постоянной ее глубине не сказывается на виде кривых силовых параметров разрушения, в то время, как при изменении глубины трещины при постоянной ее длине происходит изменение характера КИН третьего рода. Максимальные значения силовых параметров разрушения составили К1 = 234,2 МПахмм1/2, К2=205,3 МПахмм1/2, Кз= -176,06 МПа*мм1/2. На рис.10 показаны зависимости интеграла Черепанова-Райса по фронту трещины, полученные по шести контурам интегрирования. Установлено, что для различных по форме трещины форма кривой интеграла Черепанова-Райса представляет собой кривую с максимумом в точке максимально удаленной по фронту трещины. Максимальное значение наблюдается для трещины 2С=3,2 мм; а=0,8 мм и составляет 0,93 М1/шш2.
Установлено, что характер полученных кривых типичен для силовых и энергетических параметров разрушения по фронту трещины и подтверждается авторами [2].
Рис. 10. Значения интеграла Черепанова-Райса:
а-2С=1,8 мм; б-2С=3,2 мм; а=0,4 мм
ЗАКЛЮЧЕНИЕ
В результате анализа термонапряженного состояния ЦК было установлено, что в результате действия центробежной радиальной нагрузки, в диске ЦК преобладают осевые перемещения, вызванные изгибной деформацией основания диска и пера лопатки.
Выполнен расчет силовых параметров разрушения, описаны новые зависимости КИН первого (К1), второго (К2), и третьего (К3), рода по фронту трещины. Установлены значения и диапазон изменения интеграла Черепанова-Райса по фронту трещины по шести контурам интегрирования.
Проведенное численное моделирование продемонстрировало эффективность компьютерного моделирования (МКЭ) при оценке прочности ЦК ГТД. Программный комплекс ANSYS/Workbench/CFX позволяет эффективно решать комплексные мультифизические задачи и проводить анализ прочности, а также другие виды расчетов в единой расчетной среде.
СПИСОК ЛИТЕРАТУРЫ
1. Фалалеев С. В. Современные проблемы создания двигателей летательных аппаратов. [Электронный ресурс]. Самара: Самар. гос. аэрокосм. ун-т им. С. П. Королева (нац. исслед. ун-т). 2012. 106 с. [Falaleev S. V. Current problems in the creation of aircraft engines. Samara: Samara State Aerospace University, 2012. 106 p.]
2. Раер Г. А. Динамика и прочность центробежных компрессорных машин. Л.: Машиностроение, 1968. 260с. [RAYER G. A. Dynamics and durability of centrifugal compressor cars. L.: Mechanical engineering, 1968. 260 P.]
3. ANSYS Release 15.0. Documentation. Canonsburg: ANSYS Inc., 2013.
4. Справочник по коэффициентам интенсивности напряжений. В 2 Т. / под ред. Ю. Мураками. М.: Мир, 1990. 1060 с. [ Reference book on coefficients of intensity of tension. In 2 T./under the editorship of Yu. Murakami. M.: World, 1990. 1060 p. ]
5. Кудринский В. З., Коваленко Н. В. Газодинамический расчет вертолетного турбовального газотурбинного двигателя. Рига: РВВАИУ, 1973. 50 с. [Kudrinskii V. Z., Kovalenko N. V. The gas-dynamic calculation of the turbine engine helicopter turboshaft. Riga: Riga Higher Military Aviation Engineering School, 1973. 50 p.]
ОБ АВТОРАХ
ЖЕРНАКОВ Владимир Сергеевич, зав. каф. сопр. материалов. Дипл. инж.- мех. (УАИ, 1967). Д-р техн. наук по тепловым двигателям ЛА (УГАТУ, 1992). Проф., засл. деят. науки РФ, чл.-кор. АН РБ. Иссл. в обл. механики деформируемого твердого тела.
БУДИЛОВ Игорь Николаевич, проф. каф. сопротивления материалов. Дипл. инженер по технологии и оборудованию сва-рочн. произв-ва (УАИ, 1978). Д-р техн. наук по тепловым двигателям летательн. аппаратов (УГАТУ 1998). Иссл. в обл. механики разрушения, прочности элементов ГТД и нефтехимических ТУ.
METADATA
Title: Thermal-Stress state of centrifugal compressor wheel witch of crack.
Authors: V. S. Zhernakov1, I. N. Budilov2
Affiliation:
1-2 Ufa State Aviation Technical University (UGATU), Russia.
Email: 1 zhvs@rb.ru, 2 budilov02@mail.ru.
Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 20, no. 4 (74), pp. 3-8, 2016. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).
Abstract: Completed a three-dimensional elastic analysis of state of stress of the centrifugal compressor (CC) with the operating loads and temperatures. Assessing impact of gas-dynamic and centrifugal forces on the stress-strain state of the blade. Determination of the most loaded sections, vector fields of principal stress of centrifugal compressor.
Key words: centrifugal wheel; compressor; titanium alloy; strength; stress-strain state.
About Authors:
ZHERNAKOV Vladimir Sergeyevich, head of the department of strength of materials. Degree in mechanical engineering. (UAI, 1967). Dr. of technical sciences in heat engines of aircraft (UGATU, 1992). Prof., Honored Scientist of Russia, Corresponding Member Academy of Sciences Rep. Bashkortostan. Research in the field of solid mechanics.
BUDILOV Igor Nikolaevich, prof. of the department of strength of materials. Dipl. engineering technology and welding equipment production (UAI, 1978). Dr. of technical sciences in heat engines of aircraft (UGATU, 1998). Studies of fracture mechanics, strength elements of turbine engines and petrochemical engineering construction.