Evaluation of service life of heat-resistant alloys under cyclic thermomechanical loading

Abstract


The problem of evaluating strength and service life of critical engineering facilities, whose operating properties are characterized by multi-parametric non-stationary thermomechanical effects is discussed in the paper. The main degradation mechanisms of structural materials (metals and alloys) specific to these objects are considered. The main requirements to the mathematical models of fatigue damage accumulation are formulated. The basic physical laws of complex thermoplastic deformation and fatigue damage accumulation in structural materials (metals and their alloys) for various modes of combined thermomechanical loading and their main difference from isothermal fatigue processes are considered. In the modern mechanics of damaged media (MDM), a mathematical model is developed which describes the processes of thermoplastic cyclic deformation and fatigue damage accumulation in structural alloys under multiaxial disproportionate modes of combined thermomechanical loading. A MDM model consists of three interrelated parts: defining relations of thermal plasticity accounting for their dependence on failure process, evolution equations of fatigue damage accumulation and a strength criterion of damaged material. It is shown that for certain parameters of equations of cyclic thermal plasticity, using a single point on the experimental fatigue curve, the parameters of evolution equations of damage accumulation can be determined which are used for a high-accuracy computational reconstruction of low-cycle fatigue curves for various complex deformation trajectories. The results of numerical modeling of thermoplastic cyclic deformation and fatigue damage accumulation in heat-resistant alloys (Haynes188) under combined thermomechanical loading are presented. Particular attention is paid to the issues of modeling the process of thermoplastic cyclic deformation and fatigue damage accumulation for complex processes of deformation accompanied by the rotation of the main sites of stress and strain tensors.

Full Text

Введение Тенденция развития конструкций и аппаратов современного машиностроения характеризуется усложнением рабочих режимов, снижением металлоемкости за счет оптимального проектирования и применения новых высокопрочных материалов, а также значительным ростом удельного веса нестационарного нагружения. Возрастают требования к надежности и длительности безаварийной эксплуатации как конструкций в целом, так и отдельных её элементов. Все это приводит к актуальности соответствующих задач надежной оценки ресурса. Особенно эта задача актуальна для объектов, срок службы которых составляет несколько десятков лет (атомные энергетические установки, авиационные газотурбинные двигатели, ракетные двигатели и др.). Как правило, эксплуатационные условия работы таких конструкций и аппаратов современной техники характеризуются многопараметрическими нестационарными термосиловыми нагрузками, воздействиями внешних полей, приводящими к деградации начальных прочностных свойств конструкционных материалов и в конечном счете исчерпанию ресурса материала конструктивных узлов объекта [1-6]. Для элементов конструкций, работающих при нестационарном термомеханическом нагружении, основным деградационным механизмом конструкционных материалов (металлов и их сплавов) является механизм накопления усталостных повреждений, включающий: многоцкловую усталость (МнЦУ) при квазиупругой работе материала, соответствующую долговечностям циклов при симметричном одноосном циклическом нагружении; малоцикловую усталость (МЦУ) при нестационарном термопластическом деформировании материала, соответствующую долговечностям, меньшим циклов при симметричном одноосном циклическом нагружении; переходную область, соответствующую долговечностям циклов, где одновременно действуют оба механизма деградации материала [1-3]. Особенностью работы элементов конструкций, разрушающихся в результате МЦУ, является нестационарность теплового и силового нагружения, которая определяет в местах концентрации напряжений циклический характер процесса термопластического деформирования материала, сопровождающийся выдержками различной длительности при различных уровнях напряжений и температуры [1, 2]. Указанный циклический характер термопластической деформации приводит к интенсивному накоплению повреждений, обусловленному зарождением, ростом и слиянием микродефектов, приводящих материал к разрушению за ограниченное ( ) число циклов нагружения. Эффекты деформирования при термомеханическом нагружении определяются зависимостью физико-механических характеристик конструкционных материалов от температуры, наличием ограничений температурных деформаций и соотношением скоростей механической и температурной деформаций. В связи с тем что частоты и фазы циклов изменения температуры и компонент тензора механической деформации, как правило не совпадают, неизотермическое циклическое нагружение, как правило, является нерегулярным, непропорциональным, многоосным и сопровождается вращением главных площадок тензоров напряжений и деформаций. Циклическое неизотермическое деформирование материала в отличие от изотермического деформирования представляет собой несбалансированный циклический процесс. Степень дисбаланса является функцией диапазона температурного цикла и сдвига фаз изменения температуры и деформации. Даже в упругой области процесс деформирования зависит от фазы изменения температуры и механической деформации и является жестким в конце холодного температурного цикла и более податливым на высокотемпературном конце (цикл растяжения-сжатия вследствие зависимости модуля упругости от температуры - нелинейная убывающая функция растущей температуры) [7]. При неизотермическом упругопластическом деформировании процесс деформирования будет иметь более сложный характер вследствие зависимости упругопластического поведения материала от температуры. Объемная температурная деформация приводит, как правило, к многоосному процессу деформирования, сопровождающемуся при сочетании с механическим деформированием вращением главных площадок тензоров напряжений и деформаций за счет непропорциональности процесса деформирования. Это приводит к новым источникам дисбаланса: напряженно-деформированное состояние в любой точке цикла отлично от других точек. Фазирование изменения температуры и механической деформации является одним из важнейших факторов, связанных и с температурной усталостью. На основании исследований закономерностей усталости при повышенных температурах в изотермических условиях выполнен достаточно большой объем работ [3, 4, 7-15], в то время как термоусталость изучена недостаточно. Основная проблема заключается в том, что физико-механические характеристики конструкционного материала существенно зависят от температуры, в то время как экспериментальные исследования усталости в основном проводятся в изотермических условиях при ожидаемой максимальной температуре цикла. Двумя простейшими и наиболее распространенными техническими предпосылками для переноса изотермического поведения на неизотермические циклические условия являются: сопротивление изотермической усталости при максимальной температуре цикла, которое эквивалентно сопротивлению температурной неизотермической циклической усталости; минимальное сопротивление изотермической усталости в диапазоне температур термического циклирования, которое эквивалентно неизотермической усталости. В качестве других предпосылок иногда используют алгебраическую среднюю или интегральную эквивалентную температуру между двумя крайними значениями. Ни одно из этих допущений не имеет глубокого физического и экспериментального обоснования. Сопоставление экспериментальных результатов неизотермической термомеханической усталости с изотермической часто дает неожиданные расхождения: реальная долговечность при неизотермической термомеханической усталости может быть значительно ниже, чем предполагаемая на основании изотермических испытаний [7]. Неизотермическое деформирование конструкционных материалов по сравнению с изотермическим характеризуется рядом специфических особенностей [7]: нелинейной зависимостью напряжение-деформация даже в упругой области вследствие зависимости модулей упругости от температуры; неуравновешенностью диаграмм при симметричных циклах изменения механической деформации и сильной зависимостью формы петли гистерезиса напряжение-деформация от сочетания фаз изменения температуры и механической деформации. Термомеханическое усталостное циклирование может инициировать механизмы повреждения, которые не развиваются при изотермической усталости. Основными параметрами, определяющими темпы процессов термомеханической усталости являются: перепад температур (температурные градиенты), максимальная температура цикла, геометрические характеристики и условия закрепления элементов конструкций, физико-механические характеристики материалов (особенно коэффициент линейного расширения), неоднородность (анизотропность) конструктивной композиции материалов. Важным выводом из сравнительного анализа экспериментальных результатов по термомеханической усталости и изотермическим усталостным испытаниям является то, что оценка долговечности по результатам изотермических испытаний при максимальной температуре цикла не является консервативной. Снижение усталостной долговечности при термомеханическом нагружении происходит в результате действия дополнительных механизмов деградации, которые отсутствуют при изотермических нагружениях (металлургические изменения структуры материала, дополнительное деформационное старение, усиление влияния окружающей среды). Таким образом, несмотря на большое внешнее сходство явлений изотермической и неизотермической усталости, использовать результаты изотермической малоцикловой усталости для оценки неизотермической усталостной долговечности следует с большой осторожностью. Существующие на сегодняшний день подходы по предсказанию усталостной долговечности при помощи полуэмпирических формул (правил), основанные на стабилизированном анализе процесса деформирования и связывающие параметры петель упругопластического деформирования с количеством циклов до разрушения, требуют большого количества экспериментальной информации и справедливы только для узкого класса режимов нагружения [2, 8]. Повреждение и разрушение материалов в основном обусловлено зарождением микродефектов, их ростом и слиянием в макроскопические трещины. Задачи описания механического поведения этих микродефектов не менее важны, чем задачи развития макродефектов (трещин), и в последние годы для решения таких задач успешно развивается новое научное направление - механика поврежденной среды (МПС) [2, 6, 16, 17-23]. Поскольку процессы накопления повреждений тесно связаны с кинетикой НДС, точность расчетных оценок прочности и ресурса конструктивных элементов будет зависеть от того, насколько определяющие соотношения достоверно описывают процессы деформирования опасных зон элементов конструкций в заданных условиях эксплуатации. Такие параметры процесса термопластического деформирования, как длина и вид траектории деформирования, вид напряжённого состояния, история его изменения и др., существенно влияют на скорости протекания процессов накопления повреждений [3, 4, 9, 10]. В работах [2, 16, 17, 24, 25] с современных позиций МПС развита математическая модель, описывающая процессы циклического вязкопластического деформирования и накопления усталостных повреждений в конструкционных материалах (металлах и их сплавах) при многоосных непропорциональных путях комбинированного термомеханического нагружения. Ниже в целях оценки модели проведено исследование влияния законов изменения механической деформации и температуры (вида траектории деформирования) на усталостную долговечность жаропрочных сплавов при непропорциональных режимах термомеханического нагружения. Аттестация проводилась путём сопоставления результатов расчёта с экспериментальными данными по изотермическому и неизотермическому (синфазному и противофазному) циклическому деформированию лабораторных образцов. Проведено исследование влияния законов изменения механической деформации и температуры (вида траектории деформирования) на усталостную долговечность жаропрочного сплава Haynes188 [26] при непропорциональных режимах термомеханического нагружения [27]. 1. Математическая модель механики повреждённой среды Модель повреждённой среды, предложенная Ю.Г. Коротких и развитая в работах его учеников [2, 16, 17, 24, 25], состоит из трёх взаимосвязанных частей: соотношений, определяющих вязкопластическое поведение материала с учётом зависимости от процесса разрушения; уравнений, описывающих кинетику накопления повреждений; критерия прочности повреждённого материала. 1.1. Соотношения термопластичности Определяющие соотношения термопластичности базируются на следующих основных положениях: тензоры деформаций и скоростей деформаций состоят из упругих , , не зависящих от истории нагружения и определяющихся конечным состоянием процесса, и пластических - , , зависящих от истории процесса нагружения; начальная поверхность текучести для различных температур описывается поверхностью в форме Мизеса; эволюция изменения поверхности текучести описывается изменением ее радиуса и перемещением ее центра ; справедлив принцип градиентальности вектора скорости пластических деформаций к поверхности текучести в точке нагружения; изменение объема тела упруго ; рассматриваются изначально изотропные среды; учитывается только анизотропия, вызванная процессами необратимого деформирования. В упругой области связь между шаровыми и девиаторными составляющими тензоров напряжений и деформаций устанавливается с помощью закона Гука: (1) где - шаровые, а - девиаторные составляющие тензоров напряжений , деформаций и их скоростей , соответственно; - температура; - начальная температура; - модуль объемного сжатия; - модуль сдвига; - коэффициент линейного температурного расширения материала. Эффекты монотонного и циклического деформирования в пространстве напряжений учитываются с помощью поверхности текучести, уравнение которой имеет вид , . (2) Для описания сложных циклических режимов деформирования в пространстве напряжений вводится поверхность циклической «памяти». Уравнение поверхности «памяти», позволяющее при расчетах отделить монотонные процессы деформирования от циклических, имеет вид (3) где - максимальный за историю нагружения модуль переменной . Эволюционное уравнение для радиуса поверхности текучести принимается в виде [2] (4) , , , (5) , , ( ), (6) , , , (7) , , (8) где - модули изотропного упрочнения, соответствующие монотонным лучевым путям нагружения (q1), излому траектории деформирования на 90º (q2), температурному изменению радиуса поверхности текучести (q3); - постоянная, определяющая скорость процесса стабилизации формы петли гистерезиса циклического деформирования материала; - стационарное значение радиуса поверхности текучести при данных и ; и - длины траекторий пластического деформирования материала при циклическом и монотонном нагружениях; - начальное значение радиуса поверхности текучести. Первый член уравнения (4) описывает изотропное упрочнение в результате монотонного пластического деформирования ( и ), второй член - циклическое упрочнение материала ( и ), а третий - изменение радиуса поверхности текучести при изменении температуры. В целом уравнение (4) описывает локальную анизотропию пластического упрочнения в зависимости от параметра , характеризующего отклонение вектора догрузки от нормали к поверхности текучести в точке нагружения. Операторы и позволяют автоматически провести разделение процессов монотонного и циклического деформирования. Уравнение для смещения поверхности текучести основано на гипотезе А.А. Ильюшина: упрочнение зависит от истории деформирования лишь на некоторой ближайшей части траектории (запаздывание векторных свойств). Внутренняя переменная , учитывающая анизотропию упрочнения пластического деформирования, удовлетворяет эволюционному уравнению [2]: , , (9) где и - модули анизотропного упрочнения. Первый и второй члены уравнения (9) отвечают за анизотропную часть деформационного упрочнения, а третий - за изменение в результате воздействия температуры . Уравнение (9) описывает известный пространственный эффект Баушингера и анизотропию векторных свойств при изменении направления деформирования (изломе траектории деформирования). Введение второго члена в это соотношение моделирует исчезающую память внутренней переменной (скорость изменения является разностью между двумя составляющими и ). Для характеристики поведения поверхности «памяти» сформулировано эволюционное уравнение для : , (10) где и далее для любой величины , заключенной в скобки , выполняются условия Компоненты тензора скоростей пластических деформаций подчиняются закону градиентальности вектора скорости пластических деформаций к поверхности текучести в точке нагружения: , (11) где - коэффициент пропорциональности, определяемый из условия прохождения новой поверхности текучести через конец вектора девиатора напряжений в конце этапа нагружения. Материальные параметры, входящие в выражения (4)-(10), находятся из различных экспериментов [2]. На стадии коллективного взаимодействия рассеянных по объему усталостных повреждений наблюдается влияние величины поврежденности на физико-механические характеристики материала. В первом приближении это влияние может быть учтено на базе концепции деградирующего континуума введения эффективных напряжений (либо эффективных модулей упругости) [2, 19, 20]: (12) Аналогичным образом определяется эффективная внутренняя переменная : . (13) 1.2. Эволюционные уравнения накопления повреждений Экспериментальный и теоретический анализ процессов повреждённости материала позволяет представить эволюционное уравнение накопления повреждений в элементарном объёме материалов в следующем общем виде [2, 18, 19]: , (14) где функция описывает влияние кривизны траектории деформирования; - вид (объёмность) напряжённого состояния; - уровень накопленной повреждённости; - накопленную относительную энергию, затраченную на образование дефектов. Конкретизация данного соотношения при усталости приводит к следующему виду уравнения накопления повреждений [2, 16, 17, 24, 25]: (15) где ; (16) ; , (17) , (18) где , , - материальные параметры, зависящие от температуры ; - функции параметра объёмности напряжённого состояния ; - интенсивность тензора напряжений; - энергия, идущая на образование рассеянных усталостных повреждений при МЦУ; - значение в конце фазы зарождения микродефектов при МЦУ, а - значение в момент образования макротрещины при МЦУ. 1.3. Критерий прочности повреждённого материала В качестве критерия окончания фазы развития рассеянных микроповреждений (стадии образования макротрещины) принимается условие достижения величины повреждённости своего критического значения: (19) 1.4. Алгоритм интегрирования определяющих соотношений МПС Определение основных характеристик процесса вязкопластического деформирования повреждённых материалов (параметров состояния) может осуществляться двумя способами: интегрирование определяющих соотношений по времени одним из методов решения задачи Коши [28] при линеаризации алгоритма определения λ и записи определяющих соотношений МПС в приращениях, которые зависят от выбранного шага . Шаг по времени корректируется при прохождении сложных участков траектории деформирования в течение всего расчетного времени при условии устойчивости вычислений. В результате численных экспериментов [29, 30] установлено, что используемая явная схема Эйлера может успешно применяться, несмотря на то, что она является лишь условно устойчивой, при условии правильно выбранной величины шага интегрирования. Главное преимущество явной схемы - простота её реализации при минимизации времени вычисления. Успешное применение явной схемы в значительной степени определяется выбором шага интегрирования . Для этого необходима методика, позволяющая следить за скоростями изменения искомых параметров, определяя шаг интегрирования . Авторы [29, 30] предлагают выбирать шаг интегрирования в зависимости от скорости изменения искомой величины . На шаге интегрирования . Шаг интегрирования и ошибка интегрирования δ определяются по формулам [29] (20) где - верхняя граница Сглаживание численных результатов достигается за счёт прямого ограничения скоростей изменения вычисляемых величин . Такой подход [1, 2, 16, 17, 24, 25] наиболее удобен для решения задач механики деформируемого твёрдого тела, в которых параметры зависят от истории процесса нагружения и даёт результаты, необходимые для инженерных расчётов, как по точности и устойчивости расчёта процесса деформирования, так и по временным затратам. Данный подход и используется в настоящей работе. Данный алгоритм был реализован в программе «EXPMODEL» [31], которая позволяет по заданной истории изменения компонент тензора полных деформаций и температуры рассчитывать параметры процессов вязкопластического деформирования и накопления повреждений в конструкционных материалах (металлах и их сплавах) при комбинированных термомеханических воздействиях. 3. Численные результаты В примере представлены результаты численных и экспериментальных данных по термоциклическому деформированию тонкостенных трубчатых образцов из жаропрочного кобальтового сплава Haynes 188 в интервале температур . Изотермические и неизотермические эксперименты проводились при жёстком нагружении: амплитуда осевой деформации составляла 0,4 %, а амплитуда сдвиговой деформации - 0,9 % [27]. В экспериментах [27] варьировались: угол сдвига фаз между осевой деформацией и деформацией сдвига (при изотермическом деформировании при - пропорциональное нагружение; при осевая и сдвиговая деформации меняются в противофазе); угол сдвига фаз между амплитудами осевой деформации и температурой Т (при неизотермическом нагружении при осевая деформация и температура меняются в фазе, при - в противофазе). В работе [27] при сложном изотермическом и неизотермическом деформировании реализовано восемь типов испытаний: эксперимент И316МФ: механические деформации меняются в фазе ( ), а температура - постоянная (изотермические испытания); эксперимент И316МПФ: механические деформации меняются в противофазе ( ), а температура - постоянная (изотермические испытания); эксперимент И760МФ: механические деформации меняются в фазе ( ), а температура - постоянная (изотермические испытания); эксперимент И760МПФ: механические деформации меняются в противофазе ( ), а температура - постоянная (изотермические испытания); эксперимент МФТФ: механические деформации и температура меняются в фазе ( и ) - неизотермические испытания; эксперимент МФТПФ: механические деформации меняются в фазе, а температура в противофазе ( и ) - неизотермические испытания; эксперимент МПФТФ: механические деформации меняются в противофазе, а температура в фазе ( и ) - неизотермические испытания; эксперимент МПФТПФ: механические деформации и температура меняются в противофазе ( и ) - неизотермические испытания. Задавая различные законы изменения осевой деформации , деформации сдвига и температуры можно реализовать различные траектории деформирования (табл. 1). Таблица 1 Законы изменения осевой деформации , сдвиговой деформации и температуры Т (траектории деформирования) Laws of changing axial strain , shear strain and temperature T (deformation trajectories) И316МФ МФТФ Окончание табл. 1 И316МПФ МФТПФ И760МФ МПФТФ И760МПФ МПФТПФ Физико-механические характеристики жаропрочного кобальтового сплава Haynes 188 [26] и материальные параметры модели повреждённой среды [2, 16, 17, 24, 25] приведены в табл. 2-4. Таблица 2 Физико-механические характеристики и параметры модели кобальтового сплава Haynes 188 Physical and mechanical characteristics and parameters of cobalt alloy model Haynes 188 Характеристики Температура (°С) 316 760 К (МПа) 167916 138400 G (МПа) 77500 63680 a (1/град) 0,00001319 0,00001531 (МПа) 295 260 g1 (МПа) 18170 15770 g2 182,74 182,67 a 4 4 1481 201 На рис. 1, 2 представлены расчётные и экспериментальные зависимости амплитуд осевых и сдвиговых напряжений в зависимости от числа циклов нагружения. Сплошными линиями отмечены результаты численного моделирования, а маркерами соответствующие экспериментальные данные (МПФТПФ - треугольники, МПФТФ - окружности, И315ТПФ - ромбы, И760ТФ - квадраты). Видно, что степень циклического упрочнения при неизотермическом термоциклическом нагружении значительно больше, чем та, которая наблюдается в изотермических испытаниях при температурах 316 и 760 °С. Эти закономерности показывают, что механизм неизотермического термоциклического деформирования качественно отличается от изотермического. Таблица 3 Зависимость значения модуля монотонного изотропного упрочнения (МПа) от длины траектории пластического деформирования на монотонных участках для кобальтового сплава Haynes188 Value of monotonic isotropic hardening module (MPa) depending on the length of the plastic deformation trajectory on monotonous sections for cobalt alloy Haynes188 Т = 316 °С 0 0,0002 0,0004 0,0005 0,001 0,002 0,003 0,005 0,01 -2401 -2716 -3032 -3190 -3979 -5296 -6089 -6186 -2133 Т = 760 °С 0 0,0002 0,0004 0,0005 0,001 0,002 0,003 0,005 0,01 -2124 -2129 -2135 -2138 -2152 -2470 -3367 -5207 -1471 Таблица 4 Зависимость значения модуля циклического упрочнения (МПа) от величины максимального значения смещения центра поверхности текучести для кобальтового сплава Haynes 188 Value of cyclic hardening (MPa) depending on the value of the maximum displacement of the center of the yield surface for cobalt alloy Haynes 188 Т = 316 °С ρmax 24,35 46,81 49,16 Q1 520 520 500 Т = 760 °С ρmax 43,88 46,51 60,75 Q1 350 440 330 Анализ полученных результатов позволяет сделать вывод, что представленная в работе модель циклического термопластического деформирования достоверно описывает экспериментальные результаты и может в дальнейшем быть применена для расчета кинетики процесса накопления усталостных повреждений при циклических термомеханических воздействиях. Для исследования влияния траектории деформирования на малоцикловую усталость жаропрочных сплавов были проведены расчёты процессов накопления усталостных повреждений при различных заданных законах изменения осевой деформации , деформации сдвига и температуры (см. табл. 1). Расчётные результаты сравнивались с имеющимися опытными данными [27]. На рис. 3 показаны экспериментальные усталостные кривые жаропрочного сплава Haynes 188 при одноосном растяжении-сжатии для температур Т = 425,650 и 870 °С соответственно (пунктирные линии). Рис. 1. Расчётные и экспериментальные зависимости амплитуды осевых напряжений в зависимости от числа циклов нагружения N Fig. 1. Computational and experimental dependences between axial stresses and the number of loading cycles N Рис. 2. Расчётные и экспериментальные зависимости амплитуд сдвиговых σ12 напряжений в зависимости от числа циклов нагружения N Fig. 2 Computational and experimental dependences between shear σ12 stresses and the number of loading cycles N Изотермические усталостные кривые при температурах Т = 571 и 760 °С (сплошные линии) получены с использованием сплайн аппроксимации (квадратичный сплайн) усталостных кривых при температурах Т = 425,650 и 870 °С [27]. Для определения материальных параметров эволюционного уравнения накопления повреждений использовалась точка Р на усталостных кривых. Материальные параметры эволюционного уравнения накопления повреждений (15)-(18) жаропрочного сплава Haynes 188 приведены в табл. 2. Рис. 3 Экспериментальные усталостные кривые жаропрочного сплава Haynes 188 Fig. 3. Experimental fatigue curves of heat-resistant alloy Haynes 188 Результаты расчётов, их сравнение с имеющимися опытными данными приведены на рис. 4, 5 и в табл. 5. На рис. 4, 5 показана зависимость величины повреждённости от числа циклов нагружения N соответственно для двух вариантов нагружения изотермических и неизотермических соответственно (см. табл. 1). Рис. 4. Зависимость величины повреждённости ω от числа циклов нагружения N Fig. 4. Dependence of damage value ω on the number of loading cycles N Результаты исследований показывают следующее. Развитая в работах [2, 16, 17, 24, 25] модель повреждённой среды с достаточной для инженерных расчётов точностью описывает процессы циклического изотермического и неизотермического деформирования и накопления усталостных повреждений при комбинированном термомеханическом нагружении и позволяет рассчитывать малоцикловую усталостную долговечность жаропрочных сплавов при Рис. 5. Зависимость величины повреждённости ω от числа циклов нагружения N Fig. 5. Dependence of damage value ω on the number of cyclic loadings N Таблица 5 Число циклов до разрушения для различных режимов термоциклического нагружения Number of cycles before fracture for different modes of thermal cyclic loading Режим нагружения Число циклов нагружения, N эксперимент расчёт Изотермический И316МФ 9448 10512 И316МПФ 3410 2860 И760МФ 910 1205 И760МПФ 1089 785 Неизотермический МФТФ 282 143 МФТПФ 569 830 МПФТФ 270 341 МПФТПФ 1004 1307 произвольных траекториях деформирования (при сравнении опытных и расчётных данных следует иметь в виду, что информация в [27] представлена без учёта разброса экспериментальных данных в пределах каждого отдельного опыта). При изотермическом циклическом термопластическом деформировании Т = 316 °С механическое непропорциональное деформирование (тест И316МПФ) приводит к снижению усталостной долговечности по сравнению с пропорциональным (тест И316МФ) практически в три раза. С ростом температуры Т = 760 °С эта зависимость не работает. Неизотермическое циклическое термопластическое деформирование приводит к снижению усталостной долговечности по сравнению с изотермическим. Наименьшая долговечность наблюдается, когда максимальное значение температурного цикла достигается в полуцикле растяжения по механическим деформациям и , а минимальное значение температурного цикла - в полуцикле сжатия (тест МФТФ), что соответствует экспериментальным данным [11]. Наибольшая долговечность (тест МПФТПФ) наблюдается, когда механические деформации и температура меняются в противофазе ( и ). Заключение Проведена оценка достоверности модели повреждённой среды [2, 16, 17, 24, 25] путём сопоставления результатов численных расчётов с известными из литературы опытными данными по малоцикловой усталостной долговечности при изотермических и неизотермических режимах комбинированного темомеханического деформирования, которые подтверждают правильность моделирования процессов сложного циклического термомеханического деформирования и усталостной долговечности жаропрочных сплавов.

About the authors

I A Volkov

Volga State University of Water Transport

Email: pmptmvgavt@yandex.ru

L A Igumnov

Research Institute of Mechanics, National Research Lobachevsky State University of Nizhny Novgorod

Email: igumnov@mech.unn.ru

References

  1. Методы обоснования ресурса ЯЭУ / A.М. Митенков, В.Б. Кайдалов, Ю.Г. Коротких [и др.]. - М.: Машиностроение, 2007. - 445 с.
  2. Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругопластических сред с повреждениями. - М.: Физматлит, 2008. - 424 с.
  3. Коллинз Дж. Повреждение материалов в конструкциях. Анализ, предсказание, предотвращение. - М.: Мир, 1984. - 624 с.
  4. Корум, Сартори. Оценка современной методологии проектирования высокотемпературных элементов конструкций на основе экспериментов по их разрушению // Теоретические основы инженерных расчетов. - 1988. - № 1. - С. 104-118.
  5. Исследование термоциклической долговечности деталей с различными углами наклона охлаждаемых каналов / Н.Г. Бычков, А.Р. Лепёшкин, А.В. Першин, А.Д. Рекин, В.П. Лукаш // Авиационно-космическая техника и технологии. - 2009. - №. 10 (67). - С. 113-117.
  6. Капустин С.А., Чурилов Ю.А., Горохов В.А. Моделирование нелинейного деформирования и разрушения конструкций в условиях многофакторных воздействий на основе МКЭ: моногр. - Н. Новгород: Изд-во Нижегород. гос. ун-та им. Н.И. Лобачевского, 2015. - 347 с.
  7. Hulford G. R. Low cycle thermal fatigue. Mechanics and Mathematical Methods // F Thermal stress II. Chapter 6. Elsevier Science Publishers B.V., 1987. - P. 329-428.
  8. Исследование малоцикловой прочности при высоких температурах / под. ред. С.В. Серенсена. - М.: Наука, 1975.
  9. Chaboche J.L. Constitutive equation for cyclic plasticity and cyclic viscoplasticity // Int. J. Plasticity. - 1989. - Vol. 5, No. 3. - P. 247-302.
  10. Bernard-Connolly M., Biron A., Bue-Quic T. Low-cycle fatigue behaviour and cumulative dormage effect of SA-516-70 steel at room and high temperature // Journal Random Fatigue Life Predictions Asme Publ. - 1980. - Р. 297-302.
  11. Гусенков А.П., Казанцев А.Г. Прочность при малоцикловом и длительном циклическом нагружении и нагреве // Машиноведение. - 1979. - № 3. - С. 59-65.
  12. Сентоглу. Влияние ограничений на термомеханическую усталость // Теоретические основы инженерных расчетов. - 1985. - № 3. - С. 74-83.
  13. Liang Jin, Pellox R.M., Xie Xishan Thermomechanical fatique behavior of a nickel base superalloy // Chin. J. Met. Sci. Technol. - 1989. - Vol. 5. - P. 1-7.
  14. Temis Y.M., Azmetov K.K., Fakeev A.I. Numerical simulation of nonisothermal plasticity and thermomechanical fatigue of turbomachinery components // Computational Plasticity XII: Fundamentals and Applications - Proceedings of the 12th International Conference on Computational Plasticity - Fundamentals and Applications, COMPLAS, 2013. - С. 1130-1141.
  15. Экспериментально-теоретическое исследование усталости материалов и конструкций в условиях высокотемпературных многоцикловых нагружений / А.А. Антипов, А.Н. Горохов, В.А. Горохов, Д.А. Казаков, С.А. Капустин // Проблемы прочности и пластичности. - 2014. - Т. 76, № 1. - С. 26-38.
  16. Волков И.А, Коротких Ю.Г., Тарасов И.С. Моделирование сложного пластического деформирования и разрушения металлов при многоосном непропорциональном нагружении // ПМТФ. - 2009. - Т. 50, № 5. - С. 193-205.
  17. Коротких Ю.Г., Волков И.А., Игумнов Л.А., Шишулин Д.Н., Тарасов И.С. Моделирование процессов неупругого деформирования и разрушения жаропрочных сплавов при циклическом термомеханическом нагружении // Проблемы прочности и пластичности. - 2015. - Вып. 77, № 4. - С. 329-343.
  18. Боднер С.Р., Линдхолм У.С. Критерий приращения повреждения для зависящего от времени разрушения материалов // Тр. Амер. о-ва инж.-мех. Сер. Д. Теор. основы инж. расчетов. - 1976. - Т. 100, № 2. - С. 51-58.
  19. Леметр Ж. Континуальная модель повреждения, используемая для расчёта разрушения пластичных материалов // Тр. Амер. о-ва инж.-мех. Сер. Д. Теор. основы инж. расчетов. - 1985. - Т. 107, № 1. - С. 90-98.
  20. Murakami S., Imaizumi T. Mechanical description of creep damage and its experimental verification // J. Méc. Théor. Appl. - 1982. - Vol. 1. - P. 743-761.
  21. Bondar’ V.S., Danshin V.V. Mathematical simulation of cyclic thermoviscoplastic deformation and fracture of materials // Strength of Materials. - 2014. - Т. 46, № 5. - С. 613-618.
  22. Бондарь В.С., Даншин В.В., Макаров Д.А. Математическое моделирование процессов деформирования и накопления повреждений при циклических нагружениях // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 2. - С. 125-152.
  23. Бондарь В.С., Даншин В.В., Семенов П.В Нелинейные процессы накопления повреждений при циклических нагружениях // Проблемы прочности и пластичности. - 2013. - Т. 75, № 2. - С. 96-104.
  24. Волков И.А, Коротких Ю.Г. Моделирование процессов усталостной долговечности материалов и конструкций при малоцикловом нагружении // Механика твёрдого тела. Изв. РАН. - 2014. - № 3. - С. 66-78.
  25. Оценка ресурсных характеристик конструкционных сталей с использованием модели деградации, учитывающих усталость и ползучесть материала / И.А. Волков, В.В. Егунов, Л.А. Игумнов, Д.А. Казаков, Ю.Г. Коротких, Ф.Н. Митенков // ПМТФ. - 2015. - Т. 56, № 6. - С. 70-83.
  26. HAYNES® 188 ALLOY. STANDART PRODUCTS by Brand or Alloy Designation H-3001B / Global Headquarters 1020 West Park Avenue, P.O. Box 9013. Kokomo, Indiana 46904-9013 (USA).
  27. Kalluri S., Bonacuse P.J. An axial-torsional termomechanical fatique testing technique // Preparade for the Symposium on multiaxial fatique and deformation testing techniques. Denver, Colorado, may 15. 1995. - 25 p.
  28. Маковкин Г.А. Моделирование циклического упрочнения при блочном непропорциональном деформировании // Прикладные проблемы прочности и пластичности. Методы решения: межвуз. сб. - М., 1997. - С. 62-69.
  29. Бантхья, Мукерджи Об усовершенствованной схеме интегрирования по времени для системы определяющих соотношений неупругой деформации с нелинейностью жёсткого типа // Теоретические основы инженерных расчётов. - 1985. - № 4. - С. 54-60.
  30. Кумар, Морджариа, Мукерджи Численное интегрирование связанных определяющих моделей неупругой деформации // Journal of Eng. Mater. Technol. - 1980. - № 1. - P. 97-100.
  31. Программная реализация процессов вязкопластического деформирования и накопления повреждений в конструкционных сплавах при термомеханическом нагружении / И.А. Волков, Л.А. Игумнов, Ю.Г. Коротких, Д.А. Казаков, А.А. Емельянов, И.С. Тарасов, М.А. Гусева // Проблемы прочности и пластичности. - 2016. Вып. 78, № 2. - С. 119-134.

Statistics

Views

Abstract - 126

PDF (Russian) - 78

Cited-By


PlumX


Copyright (c) 2016 Volkov I.A., Igumnov L.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies