РАЗВИТИЕ ТРЕЩИН В РОТОРНОЙ СТАЛИ Р2М ПРИ ПОВЫШЕННОЙ ТЕМПЕРАТУРЕ

Шлянников В.Н., Косов Д.А., Федоренков Д.И., Zhang X.C., Tu S.T. Развитие трещин в роторной стали Р2М при повышенной температуре // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. -№ 3. С. 199-210. DOI: 10.15593/peim.mech/2021.3.19

Shlyannikov V.N., Kosov D.A., Fedorenkov D.I., Zhang X.C., Tu S.T. Crack Growth in Rotor P2M Steel at Elevated Temperature. PNRPUMechanics Bulletin, 2021, no. 3, pp. 199-210. DOI: 10.15593/perm.mech/2021.3.19

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 3,2021 PNRPU MECHANICS BULLETIN

https://ered.pstu.ru/index.php/mechanics/index

DOI: 10.15593/perm.mech/2021.3.19 УДК 539.4

РАЗВИТИЕ ТРЕЩИН В РОТОРНОЙ СТАЛИ Р2М ПРИ ПОВЫШЕННОЙ ТЕМПЕРАТУРЕ

В.Н. Шлянников1, Д.А. Косов1, Д.И. Федоренков1, X.C. Zhang2, S.T. Tu2

Институт энергетики и перспективных технологий ФИЦ «Казанский научный центр РАН», Казань, Россия

2

East China University of Science and Technology

О СТАТЬЕ

АННОТАЦИЯ

Получена: 13 апреля 2021 г. Принята: 30 июля 2021 г. Опубликована: 22 октября 2021 г.

Ключевые слова:

скорость развития трещины, взаимодействие усталости и ползучести, повышенная температура, сталь Р2М.

Проведено комплексное расчетно-экспериментальное исследование скорости развития трещин при взаимодействии усталости и ползучести на компактных образцах из стали Р2М при температуре 550Х. Теоретическая часть исследования состояла в формулировке параметров сопротивления разрушению через классические и новые конституционные уравнения состояния среды с учетом накопления повреждений. Численные расчеты включали определение полей параметров напряженно-деформированного состояния для условий упругости, пластичности и ползучести, а также распределений нелинейных коэффициентов интенсивности напряжений и С*-интеграла по длине и фронту трещины для каждого испытанного компактного образца. Интерпретация экспериментальных результатов для одинаковых по форме в плане, но различной толщины образцов, дана в терминах упругих и нелинейных коэффициентов интенсивности напряжений с учетом накопления повреждений при ползучести. Получено, что скорость роста трещины при совместном действии усталости и ползучести по мере увеличения размеров трещины монотонно возрастает по сравнению с гармонической усталостью на порядок и более на образцах одинаковой геометрии. Учет поврежденности через коэффициенты интенсивности напряжений ползучести обусловливает различия в диаграммах циклического разрушения. Суперпозиция вкладов усталости и ползучести в терминах времени выдержки под нагрузкой показывает увеличение суммарной скорости роста трещины на порядок по сравнению с интерпретацией экспериментальных данных в терминах чистой ползучести.

©ПНИПУ

© Шлянников Валерий Николаевич - д.т.н., проф., e-mail: shlyannikov@mail.ru, : 0000-0003-2468-9300. Косов Дмитрий Александрович - асп., м.н.с., e-mail: dima45001@gmail.ru, : 0000-0003-1510-4884. Федоренков Дмитрий Игоревич - асп., м.н.с., e-mail: dif-96@mail.ru, : 0000-0002-7087-549Х. Xian-Cheng Zhang - Professor, e-mail: fzxuan@ecust.edu.cn.

Shan-Tung Tu - Vice-Rector, Professor, e-mail: sttu@ecust.edu.cn, ¡D: 0000-0002-8328-1666.

Valery N. Shlyannikov - Doctor of Technical Sciences, Professor, e-mail: shlyannikov@mail.ru, : 0000-0003-2468-9300. Dmitry A. Kosov - PhD Student, Junior Researcher, e-mail: dima45001@gmail.ru, ¡D: 0000-0003-1510-4884. Dmitry I. Fedorenkov - PhD Student, Junior Researcher, e-mail: dif-96@mail.ru, ¡D: 0000-0002-7087-549X. Xian-Cheng Zhang - Professor, e-mail: fzxuan@ecust.edu.cn.

Shan-Tung Tu - Vice-Rector, Professor, e-mail: sttu@ecust.edu.cn, ¡D: 0000-0002-8328-1666.

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License

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

CRACK GROWTH IN ROTOR P2M STEEL AT ELEVATED TEMPERATURE V.N. Shlyannikov1, D.A. Kosov1, D.I. Fedorenkov1, X.C.Zhang2, S.T.Tu2

Institute of Power Engineering and Advanced Technologies FRC Kazan Scientific Center, Russian Academy of Sciences, Kazan, Russian Federation East China University of Science and Technology

ARTICLE INFO

Received: 13 April 2021 Accepted: 30 July 2021 Published: 22 October 2021

Keywords:

crack growth rate, fatigue and creep interaction, elevated temperature, P2M steel.

ABSTRACT

The paper presents a comprehensive computational and experimental study of the crack growth rate during the creep-fatigue interaction in compact specimens of P2M steel at a temperature of 550 °C. The theoretical part of the study consisted in the formulation of the fracture resistance parameters through the classical and new constitutive equations of the cracked body state taking into account the accumulation of damage. Numerical calculations included the determination of the stress-strain state fields for the conditions of elasticity, plasticity, and creep, as well as the distributions of nonlinear stress intensity factors and the C*-integral along the crack length and along the crack front for each tested compact specimen. The interpretation of the experimental results for specimens of the same shape in plan but different thicknesses is given in terms of elastic and nonlinear stress intensity factors, taking into account the accumulation of creep damage. It was found that the crack growth rate under the fatigue and creep interaction increases monotonically as the crack size increases in comparison with harmonic fatigue by an order of magnitude or more on specimens of the same geometry. Taking into account the damage through the creep stress intensity factors causes differences in the cyclic fracture diagrams. The superposition of fatigue and creep contributions in terms of load holding time shows an order of magnitude increase in the total crack growth rate compared to the interpretation of the experimental data in terms of pure creep.

©PNRPU

Введение

В последние десять лет активно разрабатываются модели прогнозирования скорости развития трещин и остаточной долговечности при взаимодействии малоцикловой усталости и ползучести. Актуальность этих исследований обусловлена необходимостью обеспечения безопасной эксплуатации при соответствующей экономической эффективности ответственных элементов конструкций тепловой и ядерной энергетики, авиации и общего машиностроения. Условия эксплуатации подобных изделий относятся к одному из самых трудных в реализации теоретического, экспериментального и вычислительного плана вариантов термомеханического нагружения.

Первые предложения возможности учета взаимодействия малоцикловой усталости и ползучести в расчете общей долговечности за счет суперпозиции отдельных повреждений были высказаны в работах [1] и [2], дальнейшее развитие этого подхода представлено Skelton & Gandy [3]. Последовательный переход в наукоемких отраслях машиностроения на эксплуатацию ответственных элементов конструкций на основе принципов допускаемой повреждаемости с учетом диагностики оперативного индивидуального технического состояния потребовал детального рассмотрения развития трещин при взаимодействии усталости и ползучести. Первоначально отдельно рассматривались модели прогнозирования развития трещин усталости [4] и ползучести [5-9]. В дальнейшем были уже представлены

простые модели скорости развития трещин, относящиеся именно к условиям взаимодействии усталости и ползучести [10-14].

Рассмотрение стадии развития трещин при высокой температуре в условиях, моделирующих эксплуатационные, потребовало формулировки новых определяющих конституционных уравнений поведения среды с учетом образования, накопления и развития повреждений. Теоретические основы инженерных приложений при оценке несущей способности конструкции с учетом накопления и развития повреждений при высокой температуре заложены фундаментальными работами [15] и [16]. Дальнейшее развитие механики повреждений в основном состояло в обобщении и развитии введенных авторами функций сплошности и поврежденности с учетом многоосности исходного номинального напряженного состояния и доминирующих механизмов разрушения. При этом различают силовой (stress) и деформационный (ductile) доминирующие механизмы разрушения по телу зерна или по межзеренным границам. Наиболее развит в отечественной и зарубежной литературе механизм кавитации по границам зерен.

Дальнейшее развитие формулировок уравнений состояния на стадии формирования и роста пор в твердом теле с начальным макродефектом в виде трещины представлено в обобщающей работе [17]. В этой работе сформулированы конституционные уравнения для механизма кавитации по границам зерен с учетом плотности распределения повреждений в пористой среде

и введены основы построения полей напряжений и скоростей деформаций в локальной области вершины трещины. Более детально вопрос учета многоосного напряженного состояния в вершине трещины на стадии развития повреждений при ползучести рассмотрен в работе [18]. Данная работа получила широкое распространение в вычислительных комплексах при анализе напряженно-деформированного состояния ряда модельных задач. Основополагающими работами по деформационному доминирующему механизму разрушения в поврежденной области вершины трещины при ползучести является статья [19]. Приложение деформационной и силовой трактовки поведения материала для оценки эффектов стеснения по стадиям накопления и развития повреждений при ползучести содержится в работе [20]. Совсем недавно опубликован наиболее обстоятельный обзор [21] по моделям накопления повреждений при ползучести.

Те или иные конституционные модели накопления и развития повреждений при ползучести в первую очередь используются для анализа полей напряжений и скоростей деформации в области вершины трещины. Наиболее простое представление в форме закона Бей-ли-Нортона относится к категории неповрежденных полей, когда кинетические процессы накопления и развития повреждений по стадиям ползучести не участвуют в системе разрешающих уравнений того или иного вычислительного комплекса. Большинство работ в открытой литературе по механике повреждений относится именно к этой категории постановки численных исследований. На рубеже 1990-х годов получили свое развитие основополагающие формулировки поведения поврежденной среды [17] при исследовании НДС в области вершины стационарных трещин. Обобщение работ данного направления, доведенных до соответствующих стандартов испытаний по определению скорости роста трещин при ползучести, представлено в монографии [22].

Аналитические исследования влияния накопления и развития повреждений в формулировках конституционных уравнений состояния в области вершины трещины и, как следствие, кинетики полей напряжений и скоростей деформаций ползучести представлены в работах [23-25]. Существенным результатом работ этого направления явились установленные тренды изменения локальных поврежденных полей параметров НДС по возрастающим стадиям ползучести. Однако должного учета установленной кинетики повреждений в базовых параметрах, используемых при интерпретации экспериментальных данных по скорости развития трещин при ползучести и взаимодействии усталости и ползучести, не получилось. Это обстоятельство накладывает ограничение на использование традиционных моделей прогнозирования роста дефектов при высокой температуре в силу их консерватизма. Более того, известные модели прогнозирования скорости роста трещин при ползучести и усталости [5, 20, 6, 8] основаны исключительно на упрощенных полях напряжений и скоростей

деформаций для неповрежденных состояний, в которых отсутствует учет влияния кинетики накопления и развития повреждений в локальной зоне процесса разрушения как функций времени выдержки под нагрузкой.

На основе приведенного выше обзора определены задачи настоящей работы, которые состояли в формулировке нелинейных параметров сопротивления разрушению, их численном определении и последующем приложении к интерпретации экспериментального исследования скорости роста трещин при взаимодействии усталости и ползучести. Объектом исследований при повышенной температуре выступали компактные образцы из роторной стали Р2М.

1. Модели состояния в вершине трещины

В настоящей работе интерпретация экспериментальных результатов по скорости развития трещин в компактных образцах из стали Р2М при повышенной температуре будет представлена в терминах сопротивления разрушению для условий упругости, пластичности и ползучести. В связи с этим кратко напомним конституционные уравнения состояния среды для соответствующих условий деформирования.

Упругие поля напряжений и перемещений в вершине трещины для плоской задачи описываются в соответствии с решением путем удержания первых двух

-1/2

членов разложения по степеням г 1 в полярных координатах (г, 8), центрированных на вершину трещины, следующим образом [26]:

СТП CTJ2 Кг 'in (0) fn (0)" 1 T 0"

_СТ21 СТ22 _ <j2nr _ f21 (0) f22 (0)_ + 0 0

u К П~ "b (0,к)" T • L + +-- 8 G

v G rfj

\

где V и а - коэффициент Пуассона и угол приложения нагрузки; к = 3 - ^ и к = (3 -v)/(1 + v) для плоской деформации и плоского напряженного состояния соответственно; Т - несингулярный член и /¡¡(8), ^ (8),

^ (8, к, а), 4 (8, к, а) - известные безразмерные функции [27]; ^ - упругий коэффициент интенсивности напряжений (КИН); G - модуль сдвига. Упругий КИН для компактного образца при внецентренном растяжении, который использован в экспериментальных исследованиях, определяется стандартной формулой

K =

P

1 b\fw

2 + a/w 1 -( a/wf

0,866 + 4,64 1-13,32

V w) V w

a I2 -14,72faY -5,6faY

где Р - приложенная нагрузка; Ь и w - толщина и ширина образца; а - длина трещины.

X

Поведение материала в упругопластической области при нелинейном деформировании описывалось по классической модели:

е =

™ f _ Y

ст а —I—ст„

E E

ст<ст„

ст>стг

где Е и а7 - модуль Юнга и предел текучести при простом одноосном растяжении; п и а - показатель и коэффициент деформационного упрочнения соответственно.

Поля напряжений, деформаций и перемещений в пластической области вершины трещины с сингулярностью, зависящей от пластических свойств материала

--1

г , представлены в соответствии с классическим решением [28-30]:

Щ (r')= КР7"+1ёу (9);

-п

Щ(г,в) = аК;г^Щ( 9);

п

и? (г,9) = aKp^Hf (9),

где Кр - пластический коэффициент интенсивности

напряжений [28]. В уравнениях (5) компоненты напряжений, деформаций и перемещений приведены к безразмерному виду за счет традиционной нормировки, принятой в нелинейном анализе: а =<а/а7 ,

щ = щЕ/аЬ, г = г/Ь, где а - предел текучести мате-

т « ~р ~р

риала; /. характерный размер тела с трещинои; гг.. ¡:

нормализованы так, чтобы а, пи.: (0) = {?I|2S¡¡S¡¡^ 2 =1, где Б. - девиатор тензора напряжений. Амплитудный коэффициент или пластический КИН К в формулах (5) соответствует первоначальной формулировке [31]

Kp =

aIpw

i

n+1

Авторы [31] разработали основанный на МКЭ численный метод определения /и - интеграла, входящего

в выражение для пластического КИН К (6), для общего трехмерного случая деформирования образцов и элементов конструкций путем учета углов девиации

наклонной трещины 9*:

л

1р

Qp(6,«,(a/w)) = —(ст^)"+ cosG-

-сг

к +-

~dQ

sin6 -

—Ц-^А" )cos0 •

Для состояния упруговязкопластичного материала в условиях ползучести поля напряженно-деформированного состояния в вершине трещины определяются с учетом безразмерного параметра повреждения ю:

ст,

и* (г,9,/,со) = = < (г, 9, со)1(ВЬ) = \К1 (со)]" (9, п, и со);

Кг = Ксг1а0,

где p и q - показатели сингулярности; ^ - нелинейный КИН для ползучести, который нормализован напряжением второй устойчивой стадии ползучести а0 и имеет смысл амплитудного коэффициента. Параметр повреждения ю в уравнениях

и и' (9,л,/,со) отличаются от таковых для условий пластичности (5). Авторы [32] ввели нелинейный КИН для ползучести, который определяется через С*-интеграл [33]:

Kr = — ст,

C

о v bi:L ,

где

C* = f Br^ f +1 1 ^ + 1 (в '

cos 9 d 9-

-f

БШ 9

ст„

59

дй[ ~дв

-rcos9

_сг ""г , ^сг дйв

дй? дг

дг

d 9;

i: (9, t,ncr, ro)= fa (9, t, ncr, ffl) d9 ;

a

(9,/,«сг,со) = ^-(аГ)"" cos 9-w, +1

du

~сг \

dQ

-ст.,

к +-

<Ш,

~сг \

dQ

sin9-

ю„

J+1 (

^)cos9.

В формулах (10)-(12) B и псг - константы ползучести по конституционному уравнению Нортона. Скорость накопления деформаций ползучести в уравнении (9) авторами [15, 33] обобщена на состояние многоосного деформирования через ./¿-теорию течения следующим образом:

d е_

dt

= - Ba":,

1

1 — ю

В работе [15] введен параметр сплошности, а в [16] предложено одно из первых уравнений накопления повреждений при ползучести и введено как функция приложенных номинальных напряжений ст:

d ю _ dt

где С и т - экспериментальные константы материала. Впоследствии это уравнение было модифицировано и введено несколько более сложных моделей накопления и развития повреждений [23-25], обзор которых представлен в работе [21]. На основе критерия Писаренко-Лебедева для сложного напряженного состояния авторы [34] ввели обобщенную деформационную модель накопления повреждений при ползучести для многоосного нагружения, которая использована в настоящей работе.

dm = dbi _(1 — х)/п,- + Х_" dt

tt

1

(1 — )g

где - время до разрушения; d и g - экспериментальные константы материала; Ь, и - функции вида напряженного состояния; % - отношение напряжений разрушения при одноосном сжатии и растяжении.

Результаты экспериментов, которые представлены в настоящей работе, относятся к гармоническому циклическому нагружению и сочетанию малоцикловой усталости и ползучести. Подобные ситуации предполагают учет каждой из составляющих общего процесса накопления и развития повреждений. В литературе известны несколько моделей для описания накопления и роста повреждений при взаимодействии ползучести и усталости. Ряд авторов предполагают, что помимо хорошо известного влияния повреждений при ползучести на усталостную выносливость существует возможность влияния усталостных повреждений на поведение при ползучести. В работе [3] предложены обобщенные диаграммы, включающие нелинейное взаимодействие ползучести и усталости в виде следующего комбинированного уравнения:

ю = -

1—m 1—ю.

где ю - общее повреждение; юсг - повреждение при ползучести; к/ - повреждение при усталости. При достижении равенства ю = 1 соотношение между ними определяется выражением

Ю = 0,5

( 2 — Ю W4Ю cr — 3ю2

Уравнения (13), (15)-(17) в дальнейшем будут использованы в численных МКЭ расчетах полей параметров напряженно-деформированного состояния в испытательных образцах при взаимодействии усталости и ползучести при повышенной температуре.

2. Определение параметров сопротивления разрушению

Объектом экспериментальных исследований в настоящей работе являлись цилиндрические и компактные образцы из роторной стали Р2М для испытаний при температуре 550 °C. Цилиндрические образцы испыты-вались при одноосном статическом растяжении, малоцикловой усталости и длительном статическом нагру-жении в порядке определения констант уравнений Рам-берга-Осгуда, типа Мэнсона-Коффина и Нортона. Вторая основная часть экспериментальных работ состояла в исследовании скорости роста трещины в компактных образцах при гармоническом циклическом на-гружении и совместном действии усталости и ползучести в трапециевидном цикле нагружения с временем выдержки под нагрузкой 60 с при повышенной температуре 550 °C. Во всех случаях коэффициент асимметрии цикла нагружения R = 0,1. В таблице приведены основные свойства пластичности и ползучести стали Р2М при температуре 550 °C, идентификация которых дана в пояснениях к соответствующим вышеприведенным формулам.

Основные механические свойства стали Р2М при температуре 550 °С

Main mechanical properties of P2M steel at the temperature of 550 °C

E ГПа V Gy Mna а np B 1/(Mna)n-hr) Пег C 1/(Mna)m-hr) m

180 0,3 460 1,20 6,67 1,41 x10-12 5,07 2x10-11 2,48

Испытания на скорость роста трещин проведены при температуре 550 °С на компактных образцах одинаковой формы в плане шириной w = 50 мм. Эксперименты при гармоническом нагружении выполнены на образцах толщиной Ь = 6,2 мм и 12,5 мм, а при взаимодействии усталости и ползучести на образцах трех толщин Ь = 6,2, 12,5 и 25 мм (рис. 1, а, б).

Каждый цикл нагружения при взаимодействии усталости и ползучести составлял 70 с, в котором по 5 с

ю

n

m

занимали периоды нагрузки и разгрузки и время выдержки под нагрузкой составляло 60 с (рис. 1, в). Испытания по определению скорости роста трещин выполнены на специализированной установке УТС-10 с непрерывной атоматизированной фиксацией длины трещины методом разности электрических потенциалов и раскрытия берегов трещины керамическим датчиком перемещений (рис. 2).

dwell

/ -- loading \

! unloading' Л

time

Рис. 1. Компактный образец (a, б) и профиль его нагружения (в) при усталости +ползучесть

Fig. 1. C(T) specimen (a, b) and its loading profile (c) under fatigue+creep

Рис. 2. Установка для экспериментов при взаимодействии усталости и ползучести

Fig. 2. Creep-fatigue test equipment

б

Рис. 3. Изломы образцов для усталости

Fig. 3. Fracture surfaces for fatigue

В результате испытаний серии образцов из стали Р2М были получены поверхности разрушения каждого образца, типовые примеры которых показаны на рис. 3, а и б для гармонической усталости и взаимодействия усталости и ползучести соответственно. В процессе испытаний каждого образца при высокой температуре формировались три маркера текущего положения фронта развивающейся трещины на изломе образца за счет изменения асимметрии цикла нагружения с R = 0,1 до R = 0,5 в течение одной минуты. При этом максимальные напряжения цикла оставались неизменными, а повышалось значение минимального номинального напряжения. Таким способом обеспечивалось отсутствие эффектов торможения роста трещины из-за возможной перегрузки в цикле нагружения. Типовые экспериментальные фронты трещин показаны на рис. 4.

Информация об экспериментально полученных положениях фронтов трещин положена в основу численных расчетов для определения упругих и нелинейных параметров трещиностойкости по уравнениям (3), (6), (10), (11) в компактных образцах с использованием метода конечных элементов. Для каждого испытанного образца формировалось по три расчетных схемы МКЭ для моделирования положения прямолинейных и криволинейных фронтов трещин (рис. 4, а) для соответствующих условий нагружения. Сетки конечных элементов строились с использованием объемных 20-узловых конечных элементов второго порядка. Вершина трещины моделировалась надрезом с конечным малым радиусом кривизны p/a = 1,3-10-4 (рис. 4, б, в), где р - радиус кривизны; а - длина трещины. Минимальный размер элемента составлял 3 мкм, а общее количество узлов расчетной схемы МКЭ превышало 3 млн. Сетка конечных элементов сгущалась по направлению к вершине трещины и по мере выхода фронта трещины на свободную поверхность образца.

а

б

а

в

В

б в

Рис. 4. Моделирование фронтов в компактном образце (а) и вершины трещины (б, в) с минимальным размером элемента 3 мкм

Fig. 4. Finite element mesh for the C

Программа численных расчетов состояла в определении полей параметров напряженно-деформированн-ного состояния и дальнейшем определении упругих и пластических коэффициентов интенсивности напряжений (3), (6), а также C^-интеграла и КИН при ползучести (10), (11) в каждом испытанном образце из стали Р2М. Для учета свойств поврежденности при ползучести в программный комплекс ANSYS [35] через пользовательскую процедуру USERMAT введены новые конституционные уравнения (13)-(15) состояния среды в соответствии с силовой моделью накопления и развития повреждений при сложном напряженном состоянии. Постпроцессорная обработка результатов расчетов осуществлялась по компьютерному коду, разработанному авторами [35, 36], который включает численное определение контурных интегралов (7), (12) и собственно расчет упругих и нелинейных параметров сопротивления разрушению в условиях пластичности и ползучести с учетом накопления и развития повреждений. В порядке интерпретации экспериментальных результатов по развитию трещин необходимо располагать зависимостями параметров сопротивления разрушению от длины трещины в испытательном образце. В результате выполненных численных расчетов подобные зависимости найдены для различных точек фронта в каждом испытанном компактном образце.

На рис. 5, а приведены расчетные зависимости пластического КИН КР и КИН при ползучести ^ от относительной длины трещины a/w в компактных образцах различной толщины при совместном нагружении усталостью и ползучестью.

5

3

ь

* 2

п.

U)

1

0.5

0.3 0.4 0.S 0.6 0.7

a/w

а

b=12.5mm

1-undamage у/ '

-2-damage /Г

P=12.8kN , /J?'

/s ' >P=10.6kN

creep+

fatigue

0.4 0.5 О.е 0.7

a/w

б

Рис. 5. Зависимости нелинейных КИН от длины трещины в образцах различной толщины

Fig. 5. Nonlinear SIF dependencies on crack length in specimens of different thicknesses

На этом рисунке КИН Kcr соответствует классическому (без учета повреждений) конституционному уравнению ползучести типа Norton. Различия в нелинейных КИН КР и Kcr от вкладов пластичности при малоцикловом деформировании и ползучести при выдержке под нагрузкой составляют почти половину порядка с доминированием коэффициента интенсивности напряжений при ползучести Kcr. Из рис. 5, б, следует, что эффект поврежденности через соответствующие конституционные уравнения ползучести (13), (15) приводит к увеличению значений КИН Kcr, который зависит от уровня приложенной к образцу нагрузки.

Изломы компактных образцов на рис. 3 после испытаний демонстрируют существенное изменение кривизны фронта трещины по толщине образца. Следуя подобным экспериментальным фронтам трещин, на рис. 6

и 7 приведены распределения С*-интеграла и КИН при ползучести Kcr в зависимости от приложенной нагрузки и выбора уравнений состояния среды. Видно, что различия в значениях этих параметров в середине толщины образца z/t = 0,5 и на свободной поверхности z/t = 0,0 могут достигать более 20 %.

Рис. 6. Распределения С*-интеграла вдоль фронта трещины по толщине образца

Fig. 6. C*-integral distributions along the crack front along the specimen thickness

X.

LJ.

(Л

cieep+fatigue b=12,5mm 1-undamage  2-damage front 3

"1 2

12 front 2

./-----

/V2 front 1

.............

0 0.1 0.2 0.3 0.4 0.5

z/b

б

Рис. 7. Распределения КИН Kcr вдоль фронта трещины по толщине образца

Fig. 7. SIF Kcr distributions along the crack front along the specimen thickness

На рис. 7, б представлено сравнение распределений КИН при ползучести для классического закона Norton и силовой модели накопления повреждений (13), (15). Очевидно, что на начальной стадии развития трещины (фронт 1) при взаимодействии усталости и ползучести, как и ожидалось, различий в поведении КИН Kcr почти не наблюдается. Однако на стадии устойчивого развития трещины (фронт 2) различия в распределениях КИН по фронту трещины из-за учета накопления повреждений становятся уже существенными. На стадии окончательного разрушения образца (фронт 3) эффекты повреждений от ползучести становятся менее значимыми из-за общей текучести образца.

3. Скорость развития трещин при повышенной температуре

Одним из основных вопросов исследований развития трещин при повышенной температуре является оценка влияния вида нагружения. На рис. 8 приведены диаграммы скорости развития трещин в стали Р2М при температуре 550 °C в терминах традиционных упругих коэффициентов интенсивности напряжений К1 (3). На рис. 8, а представлено сравнение диаграмм для гармонического малоциклового нагружения с частотой 7 Гц и программного трапецевидного нагружения (см. рис. 1, в) при взаимодействии усталости и ползучести. Сравнение проведено для образцов толщиной 6,2 и 12,5 мм. Из этих данных следует, что наложение времени выдержки ползучести при максимальной нагрузке цикла приводит к ускорению роста трещины по мере увеличения ее длины. Для ситуации взаимодействия усталости и ползучести скорость развития трещин может быть интерпретирована в терминах количества циклов нагружения N или времени выдержки под нагрузкой t. На рис. 8, б, приведены диаграммы циклического разрушения для совместного действия усталости и ползучести в координатах скорости роста трещин по времени и упругого коэффициента интенсивности напряжений для компактных образцов трех разных толщин, но одинаковой формы в плане. Как и ожидалось, в терминах упругих КИН, которые имеют одинаковые значения для условий плоского напряженного состояния и плоской деформации, толщина образцов не оказывает существенного влияния на скорость роста трещин.

Солгласно стандарту ASTM E2760-10 [37], скорость роста трещин при взаимодействии усталости и ползучести должна быть интерпретирована в терминах С*-интеграла. Для этих целей в предыдущем разделе по уравнению (11) численно на основе МКЭ были рассчитаны распределения С*-интеграла по фронтам трещин в каждом испытанном образце. На рис. 9, а, показаны диаграммы циклического разрушения с использованием этих данных для компактных образцов различной толщины. Можно заметить, что изменение толщины несколько увеличивает полосу разброса экспериментальных данных с сохранением общего тренда результатов.

а

0.01

0.001

те 0-0001

о

1Е-005

: A - b=6.2mm

: < - b=12.5mm £

о - b=25.0mm

А

4P ¥

А) ¿V

г oí

У Ш

А creep*

: А fatigue

20 ДО 60 80

К, [МРа*т0-5]

б

Рис. 8. Влияние вида нагружения (a) и толщины образца (б) на скорость роста трещин

Fig. 8. Influence of loading conditions (a) and specimen thickness (b) on the crack growth rate

На рис. 9, б, представлено сравнение диаграмм скорости роста трещин при совместном действии усталости и ползучести в терминах КИН Kcr (10) для двух типов конституционных уравнений состояния среды. В отличие от традиционных подходов на основе упругого КИН и С*-интеграла, интерпретация экспериментальных данных посредством КИН при ползучести Kcr позволяет выявить эффект влияния накопления и развития повреждений, который возрастает по мере увеличения размеров трещины и длительности времени приложения нагрузки.

Особенностью интерпретации результатов испытаний при совместном действии малоцикловой усталости и ползучести является необходимость построения диаграмм циклического разрушения в терминах количества накопленных циклов нагружения или времени ползучести, которые соответствуют трапецевидному профилю программного нагружения. Выбор той или иной формы представления данных в итоге приведет к определению собственных для данного вида нагружения констант трещиностойкости материала. Однако последующее использование этих констант в задачах прогнозирова-

ния циклического разрушения элементов конструкций приведет к различным несовпадающим между собой оценкам остаточной долговечности при одном и том же виде нагружения. В связи с этим при определении скорости роста трещин целесообразно использовать комбинированные модели взаимодействия ползучести и усталости, которые учитывают эффекты времени выдержки и накопленное количество циклов нагружения.

б

Рис. 9. Влияние толщины образца (a) и накопления повреждений (б) на скорость роста трещин

Fig. 9. The influence of the specimen thickness (a) and damage accumulation (b) on the crack growth rate

В этом плане наиболее простыми являются модели линейного суммирования повреждений, которые широко используются в оценках развития трещин при взаимодействии ползучести и усталости для прогнозирования срока службы элементов конструкций. Авторы [3] предложили простое правило суперпозиции для определения скорости роста трещин при совместном действии усталости и ползучести:

da dN '

da dt

da dN

da dt

da dN

( da + l —

fatigue V уcreep

da

1

da

dN у fatigie 3600,f V dt da

- 3600f

( da Л

dN > fatigue V dt )creep

а

а

+

+

где da/dN - скорость роста трещины по циклам нагру-жения (m/cycle); da/dt - скорость роста трещины по времени выдержки под нагрузкой (m/hour) и f - частота нагружени в Гц. На рис. 10 в качестве примера приведены диаграммы циклического разрушения при расчете суммарной скорости роста трещин в терминах количества циклов нагружения по уравнению (20) и времени ползучести по уравнению (21).

б

Рис. 10. Скорость роста трещин в терминах циклов нагружения

Fig. 10. Crack growth rate in terms of number of loading cycles

Из рис. 10, а, следует, что добавление к скорости роста трещин усталости составляющей от ползучести не приводит к существенному изменению суммарных значений da/dN. В отличие от этого наложение на ползучесть циклического деформирования (рис. 10, в) приводит к увеличению суммарной скорости роста трещин

почти на порядок. Подобный эффект влияния коррелирует с результатами численного расчета нелинейных коэффициентов интенсивности напряжений, представленных на рис. 5, а, где КИН ползучести Kcr существенно превышает КИН пластичности КР. На рис. 10 диаграммы циклического разрушения представлены как функции упругого коэффициента интенсивности напряжений К1. Тем не менее установленные различия в скоростях роста трещин при взаимодействии усталости и ползучести сохранятся и при интерпретации экспериментальных данных через нелинейные параметры сопротивления разрушению в форме С*-интеграла и КИН при ползучести Kcr.

Выводы

Проведено комплексное расчетно-эксперимен-тальное исследование скорости развития трещин при взаимодействии усталости и ползучести на компактных образцах из стали Р2М при температуре 550 °C. Формулировка параметров сопротивления разрушению представлена через классические и новые конституционные уравнения состояния среды с учетом накопления повреждений. Интерпретация экспериментальных результатов дана в терминах упругих, пластических коэффициентов интенсивности напряжений и КИН при ползучести. Получено, что скорость роста трещины при совместном действии усталости и ползучести по мере увеличения размеров трещины монотонно возрастает по сравнению с гармонической усталостью на порядок и более на образцах одинаковой геометрии. Учет повреж-денности через КИН ползучести обусловливает различия в диаграммах циклического разрушения. Суперпозиция вкладов усталости и ползучести в терминах времени выдержки под нагрузкой показывает увеличение суммарной скорости роста трещины на порядок по сравнению с интерпретацией экспериментальных данных в терминах чистой ползучести.

Благодарность

Работа выполнена при финансовой поддержке РФФИ (грант №20-58-53018\20_ГФЕН).

Acknowledgments

This work was supported by the Russian Foundation for Basic Research (Grant No. 20-58-53018 \ 20_GFEN).

Библиографический список

1. Lagneborg R., Attermo R. The effect of combined low-cycle fatigue and creep of the life of austenitic stainless steels // Metal. Trans. - 1971. - Vol. 2. - P. 1821-1827.

2. Lemaitre J., Plumtree A. Application of damage concepts to predict creep-fatigue failures // Trans. ASME. - 1979. -Vol. 101. - P. 284-292.

3. Skelton R.P., Gandy D. Creep-fatigue damage accumulation and interaction diagram based on metallographic interpretation of mechanisms // Mater. High Temp. - 2008. - Vol. 25. - P. 27-54.

4. Ellyin F. Crack growth rate under cyclic loading and effect of different singularity fields // Engineering Fracture Mechanics. - 1986. -Vol. 25. - P. 463-473. DOI: 10.1016/0013-7944(86)90260-2

а

5. Nikbin K.M., Smith D.J., Webster G.A. An engineering approach to the prediction of creep crack growth // Trans. ASME. - 1986. - Vol. 108. - P. 186-191.

6. Wasmer K., Nikbin K.M., Webster G.A. Creep crack initiation and growth in thick section steel pipes under internal pressure // International Journal of Pressure Vessels and Piping. -2003. - Vol. 80. - P. 489-498. DOI: 10.1016/S0308-0161(03)00103-0

7. Yatomi M., Nikbin K.M., O'Dowd N.P. Creep crack growth prediction using a damage based approach // International Journal of Pressure Vessels and Piping. - 2003. - Vol. 80. -P. 573-583. DOI: 10.1016/S0308-0161(03)00110-8

8. Theoretical and numerical modelling of creep crack growth in a carbon-manganese steel / M. Yatomi, N.P. O'Dowd, K.M. Nikbin, G.A. Webster // Engineering Fracture Mechanics. -2006. - Vol. 73. - P. 1158-1175. DOI: 10.1016/j. engfracmech.2005.12.012

9. Fookes A.J., Smith D.J. The influence of plasticity in creep crack growth in steels // International Journal of Pressure Vessels and Piping. - 2003. - Vol. 80. - P. 453-463. DOI: 10.1016/S0308-0161(03)00100-5

10. Djavanroodi F. Creep-fatigue crack growth interaction in nickel base super alloy // American Journal of Applied Sciences. -2008. - Vol. 5. - P. 454-460. DOI: 10.3844/ajassp.2008.454.460

11. Narasimhachary S.B., Saxena A. Crack growth behavior of 9Cr-1Mo (P91) steel under creep-fatigue conditions // International Journal of Fatigue. - 2013. - Vol. 56. - P. 106-113. DOI: 10.1016/j.ijfatigue.2013.07.006

12. A novel creep-fatigue interaction damage model with the stress effect to simulate the creep-fatigue crack growth behavior / L. Xu, L. Zhao, Z. Gao, Y. Han // International Journal of Mechanical Sciences. - 2017. - Vol. 130. - P. 143-153. DOI: 10.1016/j.ijmecsci.2017.05.036

13. Zhao L., Xu L., Nikbin K. Predicting failure modes in creep and creep-fatigue crack growth using a random grain/grain boundary idealised microstructure meshing system // Materials Science and Engineering: A. - 2017. - Vol. 704. - P. 274-286. DOI: 10.1016/j.msea.2017.08.035

14. Finite element simulation of creep-fatigue crack growth behavior for P91 steel at 625 °C considering creep-fatigue interaction / H. Jing, D. Su, L. Xu, L. Zhao, Y. Han, R. Sun // International Journal of Fatigue. - 2017. - Vol. 98. - P. 41-52. DOI: 10.1016/j.ijfatigue.2017.01.004

15. Качанов Л.М. Основы механики разрушения. - М.: Наука, 1974. - 312 с.

16. Работнов Ю.Н. Ползучесть элементов конструкций. -М.: Наука, 1966. - 752 с.

17. Hutchinson J.W. Constitutive behavior and crack tip fields for materials undergoing creep-constrained grain boundary cavitation // Acta Metall. - 1983. - Vol. 31. - P. 1079-1088.

18. Hayhurst D.R. Creep rupture under multi-axial states of stress // J. Mech. Phys. Solids. - 1972. - Vol. 20(6). - P. 381-390.

19. Cocks A.C.F., Ashby M.F. Intergranular fracture during power-law creep under multiaxial stresses // Metal Science. - 1980. -Vol. 14, no. 8-9. - P. 395-402. DOI: 10.1179/030634580790441187

20. Budden P.J., Ainsworth R.A. The effect of constraint on creep fracture assessments // International Journal of Fracture. -1997. - Vol. 87. - P. 139-149. DOI: 10.1023/A:1007416926604

References

1. Lagneborg R., Attermo R. The effect of combined low-cycle fatigue and creep of the life of austenitic stainless steels. Metal. Trans, 1971, vol. 2, pp.1821-1827.

21. Meng Q., Wang Z. Creep damage models and applications for crack growth analysis in pipes: a review // Engineering Fracture Mechanics. - 2019. - Vol. 205. - P. 547576. DOI: 10.1016/j.engfracmech.2015.09.055

22. Saxena A. Nonlinear fracture mechanics for engineers. -CRC Press LCC, 1998. - 472 p.

23. Bassani J.L., Hawk D.E. Influence of damage on crack-tip fields under small-scale-creep conditions // International Journal of Fracture. - 1990. - Vol. 42. - P. 157-172. DOI: 10.1007/BF00018384

24. Meng Q., Wang Z. Asymptotic solutions of mode I steady growth crack in materials under creep conditions // Acta Mechanica Solida Sinica. - 2015. - Vol. 28(5). - P. 578-591. DOI: 10.1016/S0894-9166(15)30051-3

25. Murakami S., Hirano T., Liu Y. Asymptotic fields of stress and damage of a mode I creep crack in steady-state growth // International Journal of Solids and Structures. - 2000. -Vol. 37 (43). - P. 6203-6220. DOI: 10.1016/S0020-7683(99)00267-X

26. Williams M.L. On the stress distribution at the base of a stationary crack // J. Appl. Mech. - 1957. - Vol. 24. - P. 109-114.

27. Eftis J., Subramonian N. The inclined crack under biaxial load // Engineering Fracture Mechanics. - 1978. - Vol. 10. -P. 43-67. DOI: 10.1016/0013-7944(78)90049-8

28. Hutchinson J.W. Singular behaviour at the end of a tensile crack in a hardening material // Journal of the Mechanics and Physics of Solids. - 1968. - Vol. 16. - P. 13-31. DOI: 10.1016/0022-5096(68)90014-8

29. Hutchinson J.W. Plastic stress and strain fields at a crack tip // Journal of the Mechanics and Physics of Solids. - 1968. -Vol. 16. - P. 337-347. DOI: 10.1016/0022-5096(68)90021-5

30. Rice J.R, Rosengren G.F. Plane strain deformation near a crack tip in a power-law hardening material // Journal of the Mechanics and Physics of Solids. - 1968. - Vol. 16. - P. 1-12. DOI: 10.1016/0022-5096(68)90013-6

31. Shlyannikov V.N, Tumanov A.V. Characterization of crack tip stress fields in test specimens using mode mixity parameters // International Journal of Fracture. - 2014. -Vol. 185. - P. 49-76. DOI: 10.1007/s10704-013-9898-0

32. Shlyannikov V.N., Tumanov A.V., Boychenko N.V. A creep stress intensity factor approach to creep-fatigue crack growth // Engineering Fracture Mechanics. - 2015. - Vol. 142. -P. 201-219. DOI: 10.1016/j.engfracmech.2015.05.056

33. Riedel H., Rice J.R. Tensile crack in creeping solids // ASTM STP 700, Fracture Mechanics. - 1980. - P. 112-130.

34. Shlyannikov V.N., Tumanov A.V. Creep damage and stress intensity factor assessment for plane multi-axial and three-dimensional problems // International Journal of Solids and Structures. - 2018. -Vol. 150. - P. 166-183. DOI: 10.1016/j.ijsolstr.2018.06.009

35. ANSYS Mechanical APDL Theory Reference Release 14.5. - ANSYS, Inc. Southpointe, 2012 - 275 Technology Drive, CanonBurg, PA.

36. Shlyannikov V.N., Tumanov A.V. Creep-fracture resistance parameters determination based on stress and ductility damage models // Fatigue & Fracture of Engineering Materials & Structures. -2018. - Vol. 41. - P. 2110-2129. DOI: 10.1111/ffe. 12766

37. ASTM E2760-10. Standard test method for creep-fatigue crack growth testing. Annual book of ASTM standards. Philadelphia (PA): American Society for Testing and Materials, 2010.

2. Lemaitre J., Plumtree A. Application of damage concepts to predict creep-fatigue failures. Trans. ASME, 1979, vol. 101, pp. 284-292.

3. Skelton R.P., Gandy D. Creep-fatigue damage accumulation and interaction diagram based on metallographic interpretation of mechanisms. Mater. High Temp, 2008, vol. 25, pp. 27-54.

4. Ellyin, F. Crack growth rate under cyclic loading and effect of different singularity fields. Engineering Fracture Mechanics, 1986, vol. 25, pp. 463-473. D01:10.1016/0013-7944(86)90260-2

5. Nikbin K.M., Smith D.J., Webster G.A. An engineering approach to the prediction of creep crack growth. Trans. ASME, 1986, vol. 108, pp. 186-191.

6. Wasmer, K., Nikbin, K.M., Webster, G.A. Creep crack initiation and growth in thick section steel pipes under internal pressure. International Journal of Pressure Vessels and Piping, 2003, vol. 80, pp. 489-498. DOI: 10.1016/S0308-0161(03)00103-0

7. Yatomi, M,. Nikbin, K.M., O'Dowd, N.P. Creep crack growth prediction using a damage based approach. International Journal of Pressure Vessels and Piping, 2003, vol. 80, pp. 573583. DOI: 10.1016/S0308-0161(03)00110-8

8. Yatomi, M., O'Dowd, N.P., Nikbin, K.M., Webster, G.A. Theoretical and numerical modelling of creep crack growth in a carbon-manganese steel. Engineering Fracture Mechanics, 200, vol. 73, pp. 1158-1175. DOI: 10.1016/j.engfracmech.2005.12.012

9. Fookes, A.J., Smith, D.J. The influence of plasticity in creep crack growth in steels. International Journal of Pressure Vessels and Piping, 2003, vol. 80, pp. 453-463. DOI: 10.1016/S0308-0161(03)00100-5

10. Djavanroodi, F. Creep-fatigue crack growth interaction in nickel base super alloy. American Journal of Applied Sciences. 2008, vol. 5, pp. 454-460. doi:10.3844/ajassp.2008.454.460

11. Narasimhachary, S.B., Saxena, A. Crack growth behavior of 9Cr-1Mo (P91) steel under creep-fatigue conditions. International Journal of Fatigue, 2013, vol. 56, pp. 106-113. DOI: 10.1016/j.ijfatigue.2013.07.006

12. Xu, L., Zhao, L., Gao, Z., Han, Y. A novel creep-fatigue interaction damage model with the stress effect to simulate the creep-fatigue crack growth behavior. International Journal of Mechanical Sciences, 2017, vol. 130, pp. 143-153. DOI: 10.1016/j.ijmecsci.2017.05.036

13. Zhao, L., Xu, L., Nikbin, K. Predicting failure modes in creep and creep-fatigue crack growth using a random grain/grain boundary idealised microstructure meshing system. Materials Science and Engineering: A, 2017, vol. 704, pp. 274-286. DOI: 10.1016/j.msea.2017.08.035

14. Jing, H., Su, D., Xu, L., Zhao, L., Han, Y., Sun, R. Finite element simulation of creep-fatigue crack growth behavior for P91 steel at 625°C considering creep-fatigue interaction. International Journal of Fatigue, 2017, vol. 98, pp. 41-52. DOI: 10.1016/j.ijfatigue.2017.01.004

15. Kachanov L.M. Osnovy mekhaniki razrusheniia [Fundamentals of fracture mechanics]. Moscow, Nauka, 1974, 312 p.

16. Rabotnov Iu.N. Polzuchest' elementov konstruktsii [Creep of structural elements]. Moscow, Nauka, 1966, 752 p.

17. Hutchinson J.W. Constitutive behavior and crack tip fields for materials undergoing creep-constrained grain boundary cavitation. Acta Metall, 1983, vol. 31, pp. 1079-1088.

18. Hayhurst D.R. Creep rupture under multi-axial states of stress. J. Mech. Phys. Solids, 1972, vol. 20(6), pp. 381-390.

19. Cocks, A.C.F., Ashby, M.F. Intergranular fracture during power-law creep under multiaxial stresses. Metal Science, 1980, vol. 14, no. 8-9, pp. 395-402. DOI: 10.1179/030634580790441187

20. Budden, P.J., Ainsworth, R.A. The effect of constraint on creep fracture assessments. International Journal of Fracture, 1997, vol. 87, pp. 139-149. DOI: 10.1023/A:1007416926604

21. Meng, Q., Wang, Z. Creep damage models and applications for crack growth analysis in pipes: a review. Engineering Fracture Mechanics, 2019,vol. 205, pp. 547-576. DOI: 10.1016/j. engfracmech.2015.09.055

22. Saxena A. Nonlinear fracture mechanics for engineers. CRC Press LCC, 1998, 472p.

23. Bassani, J.L., Hawk, D.E. Influence of damage on crack-tip fields under small-scale-creep conditions. International Journal of Fracture, 1990, vol. 42, pp. 157-172. DOI: 10.1007/BF00018384

24. Meng, Q., Wang, Z. Asymptotic solutions of mode I steady growth crack in materials under creep conditions. Acta Mechanica Solida Sinica, 2015, vol. 28(5), pp. 578-591. DOI: 10.1016/S0894-9166(15)30051-3

25. Murakami, S., Hirano, T., Liu, Y. Asymptotic fields of stress and damage of a mode I creep crack in steady-state growth. International Journal of Solids and Structures, 2000, vol. 37 (43), pp. 6203-6220. DOI: 10.1016/S0020-7683(99)00267-X

26. Williams M.L. On the stress distribution at the base of a stationary crack. J. Appl. Mech., 1957, vol. 24, pp. 109-114.

27. Eftis, J., Subramonian, N. The inclined crack under biaxial load. Engineering Fracture Mechanics, 1978, vol. 10, pp. 43-67. DOI: 10.1016/0013-7944(78)90049-8

28. Hutchinson, J.W. Singular behaviour at the end of a tensile crack in a hardening material. Journal of the Mechanics and Physics of Solids, 1968, vol. 16, pp. 13-31. DOI: 10.1016/0022-5096(68)90014-8

29. Hutchinson, J.W. Plastic stress and strain fields at a crack tip. Journal of the Mechanics and Physics of Solids, 1968, vol. 16, pp. 337-347. DOI: 10.1016/0022-5096(68)90021-5

30. Rice, J.R, Rosengren, G.F. Plane strain deformation near a crack tip in a power-law hardening material. Journal of the Mechanics and Physics of Solids, 1968, vol. 16, pp. 1-12. DOI: 10.1016/0022-5096(68)90013-6

31. Shlyannikov, V.N, Tumanov, A.V. Characterization of crack tip stress fields in test specimens using mode mixity parameters. International Journal of Fracture, 2014, vol. 185, pp. 49-76. DOI: 10.1007/s10704-013-9898-0

32. Shlyannikov, V.N, Tumanov, A.V, Boychenko, N.V. A creep stress intensity factor approach to creep-fatigue crack growth. Engineering Fracture Mechanics, 2015, vol. 142, pp. 201219. DOI: 10.1016/j.engfracmech.2015.05.056

33. Riedel H, Rice JR. Tensile crack in creeping solids. ASTMSTP 700, Fracture Mechanics, 1980, pp. 112-130.

34. Shlyannikov, V.N., Tumanov, A.V. Creep damage and stress intensity factor assessment for plane multi-axial and three-dimensional problems // International Journal of Solids and Structures, 2018, vol. 150, pp. 166-183. DOI: 10.1016/j.ijsolstr.2018.06.009

35. ANSYS Mechanical APDL Theory Reference Release 14.5. - ANSYS, Inc. Southpointe, 2012 - 275 Technology Drive, CanonBurg, PA.

36. Shlyannikov, V.N., Tumanov, A.V. Creep-fracture resistance parameters determination based on stress and ductility damage models. Fatigue & Fracture of Engineering Materials & Structures, 2018, vol. 41, pp. 2110-2129. DOI: 10.1111/ffe.12766

37. ASTM E2760-10. Standard test method for creep-fatigue crack growth testing. Annual book of ASTM standards. Philadelphia (PA): American Society for Testing and Materials; 2010.