Full Text

Введение Опыт эксплуатации разнообразных конструкционных элементов показывает, что условия циклического нагружения и виды напряженного состояния зачастую отличаются от условий лабораторных испытаний, таких как одноосное растяжение-сжатие, изгиб или кручение. Как правило, в эксплуатации конструкционные элементы подвержены сложным (многоосным) напряженным состояниям, поэтому для оценки долговечности реальных элементов необходимо использовать многоосные усталостные критерии. а б Рис. 1. Поверхность излома в титановом (ВТ3-1) образце с критической плоскостью: а - изображение, полученное с помощью сканирующего электронного микроскопа [23]; б - схематическое изображение поверхности излома с критической плоскостью Fig. 1. The fracture surface in a titanium (BT3-1) specimen with a critical plane (a) an image obtained using the scanning electron microscope [23]; (b) a schematic image of the fracture surface with a critical plane Для оценки усталостной долговечности элементов конструкций существуют различные модели и критерии усталостного разрушения с учетом многоосного напряженного состояния [1-7]. Современные многоосные усталостные критерии позволяют оценить количество циклов N до разрушения образца или конструкционного элемента (усталостная долговечность) с учетом ориентации так называемой критической плоскости развития усталостных повреждений [8-17]. Такую площадку (элемент критической плоскости) можно наблюдать при фрактографическом исследовании поверхностей изломов разрушенных при циклических испытаниях образцов в зоне зарождении усталостной микротрещины (рис. 1). В данной работе предлагается процедура определения ориентации критической плоскости и, следовательно, усталостной долговечности образца, подверженного многоосному циклическому нагружению с произвольным сдвигом фаз для двух классических диапазонов усталости - малоциклового (МЦУ) и многоциклового (МНЦУ). Предлагаемая процедура основана на хорошо апробированном критерии [9, 18]. Отметим, что ранее в [22] была предпринята попытка определения ориентации критической плоскости для критерия [8], в котором при определении ее ориентации кроме размаха касательных напряжений дополнительно учитывается величина нормального к ней напряжения. В силу этого обстоятельства попытка удалась только для случая синфазного и противофазного циклического нагружения. Можно ожидать, что предложенная в настоящей работе методика позволит аналитически исследовать иные критерии, в формулировке которых с критической плоскостью связаны только касательные напряжения и не участвуют нормальные. 1. Усталостный критерий с определением критической плоскости Рассмотрим частицу, подверженную многоосному напряженному состоянию, которое описывается тензором напряжений , который зависит от времени. Выберем систему координат, связанную с главными напряжениями , , . Выберем площадку, ориентированную единичной нормалью . Для определения долговечности одноосного циклического нагружения образца вплоть до его усталостного разрушения существует соотношение Баскина [19], которое аналитически представляет усталостную кривую при различных коэффициентах асимметрии цикла (левая ветвь бимодальной усталостной кривой [20]): Выберем для анализа хорошо апробированный критерий многоосного усталостного разрушения, учитывающий ориентацию критической плоскости развития повреждений [9]. Этот критерий, обобщающий (как и многие другие критерии многоосного усталостного разрушения) соотношение Баскина для одноосного нагружения, выглядит следующим образом: где коэффициенты и определяются из одноосных усталостных испытаний на растяжение-сжатие с двумя разными коэффициентами асимметрии и [21]; - количество циклов до разрушения. Величина есть максимальное значение размаха сдвигового напряжения за цикл нагружения среди всех возможных плоскостей с нормальным вектором , проходящих через рассматриваемую частицу. Выражение для величины имеет следующий вид: Величина - эквивалентное пиковое гидростатическое напряжение в рассматриваемой частице за цикл нагружения: Сдвиговое напряжение на заданной плоскости определяется следующим образом: Выражение для размаха сдвигового напряжения, как следует из [22], можно привести к виду 2. Частные случаи напряженного состояния 2.1. Трехосное нагружение «растяжение-сжатие» Выберем гармонический закон изменения главных напряжений в цикле с учетом произвольного сдвига фаз. Запишем главные напряжения следующим образом: где - частота нагружения; и - произвольные сдвиги фаз; , и - средние значения, а, и - амплитуды компонент тензора напряжений . Размах главных компонент тензора напряжений за цикл нагружения можно записать в виде где моменты времени и определяются из условия максимума по времени размаха касательного напряжения на площадке с нормалью; , . Кратко величину сдвигового напряжения можно записать как , , , , ,, , со следующими дополнительными обозначениями: Моменты времени, для которых достигаются экстремальные значения размаха касательных напряжений внутри цикла при произвольном определяются так: Далее, из условия экстремума функции , i = 1,2,3, при ограничении для моментов времени и определяются компоненты вектора нормали , : Отбор компонент, соответствующих максимальным и минимальным значениям размаха касательных напряжений, проводился явным образом, так как отделить их путем исследования квадратичной формы вторых производных в данном случае не удается, в отличие от случая синфазного и антифазного нагружения [22]. Если функция не имеет максимума внутри области, то с помощью упрощенного анализа определяются значения компонент , соответствующие экстремумамна границах этой области: либо либо После вычисления функции при найденных значениях выбирается тот случай из найденных экстремумов, где эта функция имеет максимальное значение. Аналогичным образом получается выражение для определения момента времени, когда эквивалентное гидростатическое напряжение максимально. С учетом полученных значений , и компонент нормали из критерия многоосного усталостного разрушения можно определить количество циклов N для произвольных величин амплитуд напряжений при циклическом нагружении в зависимости от сдвига фаз. Проведенные расчеты показывают значительное влияние сдвига фаз на усталостную долговечность, что видно на рис. 2, где показана зависимость для конкретных значений амплитуд главных напряжений. Так, например, на рис. 2, а приведены результаты для гидростатического напряженного состояния, которое без сдвига фаз вообще не должно приводить к усталостному разрушению в связи с отсутствием касательных напряжений на любой площадке. а б Рис. 2. Усталостная долговечность, представленная в виде количества циклов N, как функция от сдвигов фаз и Fig. 2. Fatigue life in terms of N cycles as a function of shifts of phases , 2.2. Двухосное нагружение вида «кручение с изгибом» Аналогичный анализ был проведен для широко распространенного вида нагружения - кручения с изгибом. При этом также учитывается произвольный сдвиг фаз. В данном случае в циклическом процессе нагружения отличны от нуля компоненты напряжений и , которые меняются по гармоническому закону: где - частота нагружения; - произвольный сдвиг фаз; и - средние значения, а и - амплитуды компонент тензора напряжений . Размах главных компонент тензора напряжений за цикл нагружения можно записать в виде где моменты времени и определяются из условия максимума по времени размаха касательного напряжения на площадке с нормалью; , . Кратко величину сдвигового напряжения можно записать как , , , со следующими дополнительными обозначениями: Были определены моменты времени, для которых достигаются экстремальные значения размаха касательных напряжений внутри цикла при произвольном: Из условия экстремума функции , i = 1,2,3, при ограничении были определены соответствующие данным моментам времени и - компоненты вектора нормали : , Отбор компонент , соответствующих максимальным и минимальным значениям размаха касательных напряжений, также проводился явным образом. Если функция не имеет максимума внутри области , то с помощью упрощенного анализа были определены значения компонент , соответствующие экстремумам на границах этой области. Аналогичным образом было получено выражение для определения момента времени, когда эквивалентное гидростатическое напряжение максимально. Как и в предыдущем случае нагружения, в данном варианте можно определить количество циклов N для произвольных величин амплитуд напряжений многоосного циклического нагружения в зависимости от сдвига фаз. Проведенные расчеты показывают значительное влияние сдвига фаз на усталостную долговечность, что видно на рис. 3. а б Рис. 3. Усталостная долговечность, представленная в виде количества циклов N, как функция от сдвига фазы Fig. 3. Fatigue life in terms of cycles N as a function of a shift of phase 3. Сравнение рассматриваемого критерия с экспериментальными данными Для оценки корректности определения ориентации критической плоскости было проведено сравнение с немногочисленными экспериментальными результатами, с учетом сдвига фаз для углеродистой стали 0,51 % при усталостных испытаниях на изгиб с кручением [24], а также с другим известным многоосным критерием усталостного разрушения [11]. Приведем сравнительную таблицу для нагружения вида «изгиб с кручением». Первые два столбца представляют собой амплитуды изгибающих и крутильных напряжений. Третий столбец - сдвиг фаз между ними. В четвертом по седьмой приведены данные для угла ориентации, который является углом между продольной осью образца и нормалью к критической плоскости (см. рис. 1). Численные результаты были получены путем прямого перебора для рассматриваемого критерия [9], а именно, с малым шагом были рассчитаны возможные комбинации дискретных значений времени и сдвига фаз. Ориентация плоскости с наибольшим найденным значением размаха касательного напряжения была выбрана в качестве критической; аналитические значения для угла α получены с применением формул, приведенных выше. Предлагаемая аналитическая процедура приводит к результату в 200 раз быстрее, чем прямой числовой перебор. Из приведенной таблицы видно, что аналитическая процедура дает результаты, полностью совпадающие с трудоемким численным определением ориентации критической плоскости для выбранного критерия усталостного разрушения Пападопулоса [9, 18] и относительно близкие к экспериментальным и полученным на основе иного критерия усталостного разрушения с определением критической плоскости [11]. Экспериментальные и расчетные данные Experimental and calculated data , MPa , MPa экспер. [24] Carpinteri [11] численно [9] аналит. 0,00 201,11 0 0,249 0,250 0,250 0,250 162,85 195,69 0 0,194 0,193 0,193 0,193 274,68 137,34 0 0,128 0,125 0,125 0,125 141,95 171,18 0,177 0,197 0,193 0,193 255,06 127,53 0,090 0,119 0,125 0,125 147,15 177,56 0,120 0,206 0,213 0,213 255,06 127,53 0,045 0,094 0,125 0,125 152,45 184,23 0,158 0,207 0,000 0,000 264,87 132,44 0,000 0,076 0,000 0,000 308,03 63,86 0,000 0,055 0,000 0,000 4. Численный подход Рассмотрена задача усталостного разрушения диска компрессора газотурбинного двигателя в полетных циклах нагружения. Аналитическим и численным методам расчета напряженно-деформированного состояния дисков и лопаток ГТД на различные воздействия посвящена обширная литература [25-31]. Для примера был рассмотрен диск переменной толщины (модельный вариант диска компрессора газотурбинного двигателя) под воздействием центробежных нагрузок. Напряженное состояние диска рассчитывалось численно-аналитическим методом [30-31], обобщающим на неосесимметричный случай метод, разработанный в [26]. Геометрия половины поперечного сечения диска представлена на рис. 4, а, распределение по радиальной координате некоторых компонент тензора напряжений - на рис. 4, б. Значения параметров частоты вращения = 600 1/c, модулей упругости и плотности 78 МПа, 44 МПа, 4370 (титановый сплав ВТ3-1). Для данного диска была проведена оценка количества циклов N до разрушения с выбранным пределом усталости для реверсного цикла 350 МПа. Было получено несколько схожих результатов для разных частей диска, все они лежат в диапазоне от до циклов, что соответствует экспериментальным данным [30]. а б Рис. 4. Распределение напряжений в диске рассматриваемой геометрии без учета v-образных углублений Fig. 4. Stress distribution in a simplified disc of the considered geometry without v-shaped notches for blades Также был рассмотрен более реалистичный элемент конструкции - сегмент диска компрессора низкого давления газотурбинного двигателя с лопаткой, представленные на рис. 5, a. Лопатки подвержены аэродинамическим нагрузкам, лопатки и диск подвержены центробежным нагрузкам. Рассматриваются следующие типы контактных условий на границе диска и лопаток - полное сцепление и проскальзывание с трением с возможностью отлипания. Расчеты показали, что наиболее опасными с точки зрения зарождения усталостных трещин являются окрестности зон контакта диска и лопаток, имеющих вид соединения типа «ласточкин хвост». Также показано, что наилучшее соответствие расчетных и экспериментально наблюдаемых зон концентрации напряжений достигается при учете возможности отлипания и проскальзывания контактных границ диска и лопаток. Концентрация напряжений в окрестности соединения «ласточкин хвост» приводит к необходимости значительного сгущения расчетной сетки в этой области. Оценки показывают, что расчет полномасштабной модели диска с полным набором лопаток и с достаточной степенью сгущения, а также с учетом значительного времени обработки каждого контактного условия неэкономичен по времени. В то же время при рассмотрении сектора диска с единственной лопаткой остаются неизвестными граничные условия на его боковых поверхностях, так как в силу разворота лопаток относительно оси вращения и самого вращения задача не является осесимметричной. Поэтому расчет проводился в два этапа. Сначала была создана сетка с умеренным сгущением для полной модели. Количество элементов ~. В дальнейшем для задачи на секторе диска с сильным сгущением сетки в качестве условий на боковых поверхностях использовались перемещения, полученные из расчета на полной модели. Сетка значительно сгущена в окрестности ожидаемых концентраторов напряжений. Общее количество элементов не превосходит 100 000, что вполне приемлемо для проведения расчетов на персональном компьютере. По результатам конечно-элементного расчета на рис. 5, б показано распределение напряжений в зоне контакта диска и лопатки. На рис. 5, в показано количество циклов N до разрушения в наиболее нагруженной области диска - в углу паза для закрепления лопатки для трех критериев (слева направо): классических критериев [6] и [7], сформулированных для инвариантов тензора напряжений и не использующих понятие критической плоскости, и для критерия [9], учитывающего ориентацию критической плоскости. Критерии [6] и [9] дали оценку долговечности дисков ГТД на уровне 35 000-50 000 циклов. Критерий [7] предсказал возможность усталостного разрушения в 20 000 полетных циклов. В целом все эти критерии дали сходное расположение зон усталостного разрушения. a б в Рис. 5. Оценка долговечности диска компрессора; а - сегмент диска компрессора и лопатка; б - напряженное состояние в зоне контакта диска и лопатки; в - количество циклов N до разрушения в наиболее нагруженной области диска для критериев [6, 7, 9] Fig. 5. Fatigue life assessment of the disc; (a) is a segment of a compressor disc and a blade, (b) is the stress state within the contact zone between the disc and the blade, (c) is the amount of cycles N to breakdown within the most vulnerable zone of the disc for criteria [6, 7, 9] Заключение Предложена процедура аналитического вычисления ориентации критической плоскости для многоосного циклического нагружения при произвольном сдвиге фаз в случаях классических усталостных диапазонов - малоцикловой и многоцикловой усталости. Проведено сравнение исследуемого критерия с экспериментальными данными и численными расчетами, произведенными с использованием иного критерия. Рассмотрен диск компрессора газотурбинного двигателя, подверженный полетным циклам нагружения. Путем распределения напряжений, рассчитанного приближенным численно-аналитическим методом, а также распределения напряжений, рассчитанного методом конечных элементов, определены зоны наибольшей концентрации напряжений, ориентация критической плоскости в зонах концентрации и получены оценки долговечности (число циклов до разрушения N ~ 35 000 - 50 000) эксплуатации диска с использованием выбранного многоосного критерия усталостного разрушения.

About the authors

N G Burago

Ishlinsky Institute for Problems in Mechanics of the Russian Academy of Sciences

I S Nikitin

Institute for Computer Aided Design of the Russian Academy of Sciences

A D Nikitin

Institute for Computer Aided Design of the Russian Academy of Sciences

B A Stratula

Institute for Computer Aided Design of the Russian Academy of Sciences

References

Statistics

Views

Abstract - 251

PDF (Russian) - 202