Full Text

1. Определение критической плоскости для многоосного напряженного состояния в режимах МЦУ и МНЦУ Для оценки усталостной долговечности различных элементов конструкций в реальных условиях эксплуатации существуют различные модели и критерии усталостного разрушения с учетом многоосного напряженного состояния [1-5]. В современных подходах к построению многоосных критериев усталостного разрушения часто используется понятие критической плоскости [6-14]. Один из первых критериев с определением критической плоскости для классических режимов малоцикловой (МЦУ, число циклов до разрушения 103 < N < 105) и многоцикловой (МНЦУ, число циклов до разрушения 105 < N < 107) усталости был предложен Финдли [15]. Согласно этому критерию, развитие усталостных повреждений в циклическом процессе нагружения происходит по плоскости с нормалью n (с компонентами nk, k = 1, 2, 3), на которой максимум комбинации достигает определенного критического значения. Такая плоскость называется критической. В этой комбинации - размах касательного напряжения в цикле, - нормальное напряжение на этой плоскости, - тензор напряжений, который определяет напряженное состояние в материальной частице деформируемого тела, и , . Для определения долговечности одноосного циклического нагружения образца вплоть до его усталостного разрушения существует соотношение Баскина [16], которое аналитически представляет усталостную кривую при различных коэффициентах асимметрии цикла (левая ветвь бимодальной усталостной кривой [17] (рис. 1): где - предел усталости; - коэффициент усталостной прочности; - показатель усталостной прочности; N - число циклов до разрушения. На рис. 1 по оси ординат схематически показаны значения амплитуд циклического процесса нагружения при различных коэффициентах асимметрии цикла в зависимости от N. Рис. 1. Бимодальная усталостная кривая для режимов МЦУ и СВМУ Fig. 1. Bimodal fatigue curve for LCF and VHCF modes Обобщение соотношения Баскина на многоосный случай для критерия Финдли имеет вид [1] где < 0, , , - параметры, определяемые по данным эксперимента. Расчет усталостной долговечности по критерию Финдли требует определения ориентации плоскости, проходящей через заданную материальную точку, на которой внутри цикла достигается максимальное выражение функции Финдли Для многоосного напряженного состояния это непростая задача, которая при необходимости решается численно. Построим ее аналитическое решение для синфазного и противофазного циклического нагружения. 1.1. Синфазное многоосное циклическое нагружение Рассмотрим трехосное циклическое нагружение в системе координат, связанной с главными осями тензора напряжений. Предполагаем, что эти оси не изменяются в течение цикла, а главные значения тензора напряжений меняются по гармоническому закону без сдвига фаз относительно друг друга: , , , где дополнительным индексом m помечены средние за цикл значения напряжений, а индексом a их амплитуды. Размах главных напряжений в цикле , , . Выберем главные оси так, чтобы максимумы главных напряжений удовлетворяли неравенствам , где , , . Введем следующие обозначения: , , , , , . В главных осях тензора напряжений можно получить формулы (см. Приложение) , С учетом этих формул определим ориентацию критической плоскости с компонентами нормали: , , , . Можно показать, что задача определения критической плоскости для трехосного напряженного состояния сводится к определению максимума функции (см. Приложение): при ограничениях , , . Приведем результаты решения этой задачи для всевозможных значений максимумов и размахов главных напряжений. Вначале ищем экстремум функции внутри зоны ограничений , , : , . При этом возникают следующие возможные случаи соотношений между главными значениями напряжений. I. Случай Примем обозначения , , , I-а. Если , условия экстремума имеют вид , . I-б. Если , условия экстремума имеют вид . II. Если , необходимо , . Условия экстремума : . Найденные значения , должны удовлетворять условиям , , и условиям максимума которые имеют вид Эти неравенства следуют из условия отрицательности квадратичной формы Если таких значений не существует, то необходимо искать максимум на границах зоны ограничений: либо , либо либо III. В этом случае максимум достигается при значениях , Из этих трех пар следует выбрать ту, где функция принимает наибольшее значение. Таким образом, для режима синфазного многоосного циклического нагружения получено аналитическое решение задачи определения ориентации критической плоскости. Определив компоненты нормали к этой плоскости, можно вычислить значение функции Финдли и соответствующее ему число N, т.е. усталостную долговечность элемента конструкции, находящегося в многоосном напряженном состоянии [18]. 1.2. Противофазное многоосное циклическое нагружение Важным практическим примером противофазного циклического многоосного нагружения являются усталостные испытания на чистое кручение или изгиб с кручением. Нетрудно видеть, что в системе координат, связанной с главными осями тензора напряжений, противофазное циклическое нагружение с гармоническим законом изменения по времени всегда можно представить в виде В этом случае , , . Выражение для по сравнению с синфазным циклическим режимом не меняется: . Нормальное напряжение на площадке с нормалью n в противофазном режиме зависит от времени по закону Поэтому для определения , входящего в функцию Финдли, рассмотрим два случая. I. Если , то при =1, где , , . Перенумеруем значения , k = 1, 2, 3 так, чтобы они удовлетворяли неравенствам . После этого проводится такой же анализ и расчет ориентации критической плоскости для перенумерованных значений , , k = 1, 2, 3, как и для синфазного режима. Полученные значения компонент нормали проверяются на выполнение условия . Если это условие не выполняется, то рассматриваем следующий альтернативный случай. II. Если , то при = -1, где , , . И в этом случае перенумеруем значения , k = 1, 2, 3 так, чтобы они удовлетворяли неравенствам . После этого проводится такой же анализ и расчет ориентации критической плоскости для перенумерованных значений , , k = 1, 2, 3, как и для синфазного режима. Полученные значения компонент нормали проверяются на выполнение условия . Если выполняются оба условия, то следует выбрать такую ориентацию нормали, при которой значение функции Финдли имеет большее значение. 1.3. Определение параметров многоосного критерия усталостного разрушения Ранее в [18] была предложена процедура определения параметров классических многоосных критериев (Сайнса [19], Кроссланда [20]) в диапазоне МЦУ и МНЦУ по результатам одноосных испытаний с двумя разными коэффициентами асимметрии цикла (левая ветвь бимодальной одноосной усталостной кривой [17] (рис. 1). Приведем результаты этой процедуры для определения параметров многоосного критерия Финдли. Многоосный критерий Финдли для одноосного нагружения с коэффициентом асимметрии цикла будет выглядеть так: при Сравнивая с представлением Баскина одноосной усталостной кривой и имея данные испытаний при и , получаем где и - классические пределы усталости по усталостным кривым при коэффициентах асимметрии цикла -1 и соответственно; - статический предел прочности; . Анализ критерия Финдли для реверсного режима чистого кручения дает формулы После вычисления этих параметров по результатам двух одноосных испытаний с разными коэффициентами асимметрии цикла задача вычисления ориентации критической плоскости и долговечности для многоосного напряженного состояния решается с помощью полученных выше формул для компонент нормали и уравнения для числа циклов до разрушения N. 2. Обобщение многоосного критерия усталостного разрушения в режиме СВМУ В настоящее время установлено, что относительно небольшие циклические напряжения (меньшие классического предела усталости материала), действующие с высокой частотой (порядка 1 кГц и выше), могут приводить к разрушению конструкций [21-25]. Высокочастотное нагружение приводит к значительным наработкам (N ~ 109 - 1010 циклов) за время расчетного срока службы изделия. Указанный диапазон долговечностей N > 108 известен в литературе как сверхмногоцикловая усталость (СВМУ). В последние годы разработаны и реализованы методики СВМУ испытаний [25-27] для весьма ограниченного набора циклических нагружений, в первую очередь для реверсивного и пульсирующего растяжения-сжатия [26, 27], а также кручения образцов [28-30]. Однако для оценок долговечности различных элементов конструкций, подверженных высокочастотным длительным воздействиям в условиях сложного нагружения, необходимы критерии многоосного СВМУ-разрушения. В данной работе предлагается простая и естественная форма такого рода критериев, основанная на обобщении известных многоосных моделей для классических усталостных режимов. Опишем предлагаемую схему обобщения многоосного критерия разрушения на режим СВМУ (см. рис. 1, правая ветвь одноосной бимодальной усталостной кривой при N > 108). В основе обобщенного многоосного критерия типа Финдли лежит подобие поведения левой и правой ветвей бимодальной усталостной кривой. Это обобщение заключается в замене параметров левой ветви бимодальной усталостной кривой на параметры правой ветви: , , , , где и - «новые» пределы усталости на правой ветви усталостной кривой для коэффициентов асимметрии R = -1 и R = R0. Для режима СВМУ имеем (все параметры отмечены знаком ~) , , , Для реверсного режима чистого кручения получаем , Рис. 2. СВМУ экспериментальные данные для сплава ВТ3-1: а - растяжение-сжатие; б - кручение Fig. 2. VHCF experimental data for the alloy VT3-1: а - tension-compression; b - torsion Были определены параметры обобщенного критерия СВМУ разрушения по результатам испытаний титановых образцов: МПа, МПа, МПа, МПа, . На рис. 2, а приведены экспериментальные данные для R = -1 (растяжение-сжатие). Полученные значения параметров были использованы для расчетной оценки долговечности образцов при СВМУ-испытаниях на реверсное кручение. Сравнение показало близость расчетных и экспериментальных результатов (см. рис. 2, б, где сплошной линией проведена расчетная ветвь усталостной кривой). Таким образом, подтверждена возможность использования обобщенного критерия многоосного усталостного разрушения в режиме СВМУ и предложенной схемы определения его параметров для приближенной оценки долговечности элементов конструкций. Формулы для параметров обобщенного СВМУ-критерия и его формулировка на основе подобия левой и правой ветвей усталостной кривой позволяют определить ориентацию критической плоскости и для этого режима циклического нагружения. 3. Примеры расчета долговечности в МЦУ и СВМУ-режиме Рассмотрим задачу усталостного разрушения диска компрессора газотурбинного двигателя в полетных циклах нагружения в условиях малоцикловой усталости. Аналитическим и численным методам расчета напряженно-деформированного состояния дисков и лопаток ГТД на различные воздействия посвящена обширная литература [31-35]. В нашем случае предполагается, что циклом многоосного нагружения системы диск-лопатка является полетный цикл нагружения (ПЦН, взлет-полет-посадка). В качестве границ цикла примем нагрузки на крейсерской скорости полета воздушного судна и соответствующие угловые скорости вращения диска компрессора. Задача состоит в определении долговечности диска (число ПЦН до разрушения) из многоосного критерия Финдли для режима МЦУ. Для этого был проведен конечно-элементный расчет напряженно-деформированного состояния (НДС) системы диск-лопатка при учете центробежных нагрузок, распределенного аэродинамического давления на лопатку, а также сил нелинейного контактного взаимодействия диска, лопаток и иных дополнительных учитываемых элементов конструкции (рис. 3, а). Численные методы решения контактных задач описаны в работах [36, 37]. Параметры расчета принимались такими: угловая скорость вращения w = 314 рад/с (3000 об/мин), скоростной напор на бесконечности = 26000 Н/м2, что соответствует скорости потока 200 м/с при плотности 1,3 кг/м3. Свойства материалов принимались следующими: диск (титановый сплав): E = 116 ГПа, n = 0,32, r = 4370 кг/м3; лопатки (алюминиевый сплав): E = 69 ГПа, n = 0,33, r = 2700 кг/м3; фиксирующие штифты (сталь): E = 207 ГПа, n = 0,27, r = 7860 кг/м3. Расчеты показали, что наиболее опасными с точки зрения зарождения усталостных трещин являются окрестности зон контакта диска и лопаток, имеющих вид соединения типа «ласточкин хвост». На рис. 3, б представлена зона концентрации максимальных растягивающих напряжений в левом (скругленном) углу паза в диске, в который вставляется лопатка. а б Рис. 3. Сектор диска с лопаткой, нагружаемый центробежными силами (а); максимальное главное напряжение в соединении диска и лопатки (б) Fig. 3. Is the sector of the disk with the blade which is loaded with centrifugal forces (a), is the maximum main stress in the connections of the disk and the blade (b) Наименьшая долговечность получилась для окрестности левого угла контактного соединения диска и лопатки, которая выделена на рис. 4, а сплошными линиями. На рис. 4, б показаны рассчитанные величины числа полетных циклов до разрушения N для выбранного критерия многоосного усталостного разрушения. По горизонтальной оси отложена безразмерная координата закругления левого угла паза соединения, по вертикальной оси отложена безразмерная координата по глубине паза. Наименьшее значение N составило величину ~30 000 циклов. Если принять среднюю протяженность цикла за 2 ч, то в реальном времени долговечность элемента конструкции составит 60 000 ч. Дополнительно исследовался альтернативный механизм циклического нагружения, связанный с наблюдаемыми высокочастотными осевыми колебаниями бандажных полок. Амплитуды этих колебаний и вызываемые ими возмущения напряженного состояния в окрестности концентратора напряжений относительно невелики, однако число циклов при высокой частоте колебаний за время эксплуатации может достичь величин 109-1010, что приводит к необходимости исследования режима сверхмногоцикловой усталости (СВМУ) и возможности усталостного разрушения при уровне напряжений, меньшем классического предела усталости материала. На рис. 5, а показана схема приложения низкоамплитудных а б Рис. 4. Выделенная зона максимальной концентрации напряжений (а); уровни долговечности в выделенной зоне (б) Fig. 4. Is the maximum stress concentration zone (a); is the durability levels in the allocated zone (b) Рис. 5. Схема высокочастотных воздействий на диск и лопатку (а); уровни долговечности в выделенной зоне (б) Fig. 5. Is the scheme of high-frequency loading on the disk and blade (a), is the levels of durability in the allocated zone (b) осевых смещений. Таким образом, рассматриваются случаи, когда на левую границу приходится максимум или минимум амплитуды колебаний, равной ±1 мм, при частоте 3000 об/мин. Как и в случае режима МЦУ, зона концентрации напряжений сосредоточена в левом углу паза. Несмотря на относительно невысокий уровень амплитуд напряжений в цикле, и в этом случае возникают зоны возможного усталостного разрушения, расположенные в тыльной зоне левого угла паза диска, примерно там же, где и в случае полетных циклов нагружения. Величина N, при которой они появляются, имеет порядок 109-1010 (см. рис. 5, б), что в пересчете на реальное время процесса с периодом цикла 0,02 с (при выбранной частоте колебаний) дает величину 50 000 ч, вполне достижимую в процессе эксплуатации. Таким образом, значения долговечностей в реальном времени для двух различных режимов усталостного разрушения получились достаточно близкими. Хотя эти оценки весьма приблизительны, они указывают на принципиальную возможность развития усталостного разрушения в указанных зонах диска компрессора как по механизму МЦУ (полетный цикл нагружения), так и по механизму СВМУ (высокочастотные, низкоамплитудные вибрации элементов конструкции). Заключение Получено аналитическое решение задачи определения ориентации критической плоскости для многоосного напряженного состояния при синфазном и противофазном циклическом нагружении для классического диапазона МЦУ-МНЦУ (левая ветвь бимодальной усталостной кривой). Построено обобщение многоосного критерия усталостного разрушения на случай сверхмногоцикловой усталости, предложена процедура определения параметров этого обобщенного критерия по результатам одноосных экспериментов для двух коэффициентов асимметрии цикла, определена критическая плоскость развития усталостных повреждений для многоосного напряженного состояния (правая ветвь бимодальной усталостной кривой). Проведены расчеты напряженного состояния, и на этой основе даны сравнительные оценки усталостной долговечности ответственного элемента конструкции - диска компрессора газотурбинного двигателя для режимов нагружения МЦУ и СВМУ. Анализ показал, что с учетом характерных частот рассматриваемых циклических процессов значения долговечностей в реальном времени могут быть достаточно близки и, следовательно, при определении условий безопасной эксплуатации должны учитываться оба механизма усталостного разрушения. Приложение Выведем полезные формулы для касательного напряжения и его размаха на произвольно ориентированной площадке с нормалью n(n1, n2, n3). В системе координат, связанной с главными осями тензора напряжений, имеем для нормального напряжения на выбранной площадке и компонент касательного напряжения следующие формулы (по повторяющимся индексам k и l производится суммирование): , Для размаха компонент касательного напряжения и его модуля имеем , , Раскрывая скобки, получаем , так как , . Окончательно . Раскрывая скобки, эту формулу можно преобразовать к виду Введя обозначения , , , получим компактное выражение для размаха касательных напряжений . Исключая , получим формулу, используемую для вычисления максимума функции F:

About the authors

I S Nikitin

Institute for Computer Aided Design of the Russian Academy of Sciences; Moscow Aviation Institute (National Research University)

N G Burago

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

A D Nikitin

Institute for Computer Aided Design of the Russian Academy of Sciences; Moscow Aviation Institute (National Research University)

V L Yakushev

Institute for Computer Aided Design of the Russian Academy of Sciences

References

Statistics

Views

Abstract - 289

PDF (Russian) - 199