Известия РАН. Механика твердого тела, 2020, № 4, стр. 158-164

Моделирование ползучести и длительной прочности стержней при переменных напряжениях

Ю. Г. Басалов a, А. М. Локощенко a*, В. В. Терауд a

a Научно-исследовательский институт механики Московского государственного университета им. М.В. Ломоносова
Москва, Россия

* E-mail: loko@imec.msu.ru

Поступила в редакцию 03.10.2019
После доработки 08.10.2019
Принята к публикации 24.10.2019

Полный текст (PDF)

Аннотация

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

Ключевые слова: ползучесть, длительная прочность, моделирование, переменные напряжения, растяжение, кручение, поврежденность

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

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

3. Описание результатов испытаний никелевого сплава ХН55МВЦ. В [2] приведены результаты испытаний никелевого сплава ХН55МВЦ (ЧС57) при растяжении в условиях циклического приложения кусочно-постоянного напряжения и разгрузки при температуре 900°С. Выбор для исследования данного структурно стабильного сплава объясняется тем обстоятельством, что период формирования равновесной структуры в нем на порядок меньше, чем его долговечность. На рис. 1 приведены кривые ползучести, полученные либо в случае стационарного нагружения при σ0 = 20 МПа (кривая 1), либо в случае циклического нагружения, осуществляющегося в следующей последовательности: мгновенно приложенное напряжение σ0 = 20 МПа, действующее в течение ${{t}_{1}} = 25$ час, затем мгновенная полная разгрузка (σ = 0) и последующий отдых (σ = 0) в течение ${{t}_{2}} = 50$ час (кривая 2). Испытания как при постоянном напряжении, так и в течение 10 последовательных циклов нестационарного нагружения проводились на цилиндрических образцах диаметром 7 мм и длиной рабочей части 70 мм. Эксперименты показали, что в процессе отдыха наблюдается обратная ползучесть. Кривые ползучести как при σ0 = 20 МПа, так и при σ = 0 практически не зависят от номера цикла, т.е. фактически повторяют соответствующие кривые первого цикла.

Рис. 1

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

$\dot {p} = G{\text{sign}}[{{(\sigma - \chi )}^{n}}],\quad \dot {\chi } = A\left( \sigma \right)\dot {p} - D\left( {{\text{sign}}\chi } \right){{\chi }^{n}},\quad A\left( \sigma \right) = C{{\sigma }^{{ - m}}}$
в которых p – деформация ползучести, $\chi $ – микронапряжение, G, $n$, D, C, m – материальные константы, точка означает дифференцирование по времени t. В [3] предложен метод вычисления пяти приведенных материальных констант. Применение этого метода основано на использовании результатов испытаний при двух постоянных уровнях положительного осевого напряжения σ0 = 20 и 14 МПа, однако в [3] приведены экспериментальные данные только при σ0 = 20 МПа. На рис. 1 в качестве примера моделирования экспериментальных результатов приведены расчетные кривые ползучести при стационарном и нестационарном нагружении (кривые 3 и 4 соответственно) при учете либо изотропного, либо анизотропного упрочнения. Однако в тексте статьи [3] учет анизотропного упрочнения не проводится, значения используемых материальных констант G, $n$, $D$, $C$, m также не приводятся.

Авторы данной статьи для описания ползучести сплава ХН55МВЦ в рассматриваемых условиях используют вариант кинетической теории [5]. В этом варианте учитываются неустановившаяся стадия ползучести и явление последействия.

Рассмотрим систему определяющего и кинетического уравнений с тремя материальными константами в виде

(3.1)
$\dot {p} = A\left( {\sigma - q} \right),\quad p\left( {t = 0} \right) = 0$
(3.2)
$\dot {q} = B\sigma - Cq,\quad q\left( {t = 0} \right) = 0$
здесь и далее $q(t)$ – кинетический параметр. Система дифференциальных уравнений (3.1), (3.2) при $\sigma = {{\sigma }_{0}}$ приводит к следующим зависимостям $q(t)$ и $\dot {p}(t)$:

(3.3)
$q\left( t \right) = \frac{B}{C}{{\sigma }_{0}}\left[ {1 - \exp \left( { - {\kern 1pt} Ct} \right)} \right],\quad 0 \leqslant t \leqslant {{t}_{1}}$
(3.4)
$\dot {p}\left( t \right) = A\sigma _{0}^{{}}\left[ {1 - \frac{B}{C}\left( {1 - \exp \left( { - {\kern 1pt} Ct} \right)} \right)} \right],\quad 0 \leqslant t \leqslant {{t}_{1}}$

Интегрирование (3.4) приводит к следующей зависимости $p(t)$:

(3.5)
$p\left( t \right) = \frac{A}{B}\sigma _{0}^{{}} \cdot \int\limits_0^t {\left[ {1 - \frac{B}{C}\left( {1 - \exp ( - {\kern 1pt} Ct)} \right)} \right]dt} ,\quad 0 \leqslant t \leqslant {{t}_{1}}$

Зависимости (3.3) и (3.5) при постоянном растягивающем напряжении $\sigma (t) = {{\sigma }_{0}}$ описывают постепенное увеличение параметра q(t) во времени от $q(0) = 0$ до $\mathop {\lim }\limits_{t \to + \infty } q(t)$ = = $\frac{B}{C}{{\sigma }_{0}}$ и постепенное увеличение деформации ползучести во времени в соответствии с уравнением (3.5).

Непосредственно перед разгрузкой ($t = {{t}_{1}} = 25$ час) q(t) и p(t) принимают следующие значения:

(3.6)
$q({{t}_{1}}) = {{q}_{1}} = \frac{B}{C}{{\sigma }_{0}}\left[ {1 - \exp \left( { - {\kern 1pt} C{{t}_{1}}} \right)} \right],\quad {{p}_{1}} = p\left( {{{t}_{1}}} \right)$

При разгрузке ($t > {{t}_{1}} = 25$ час) из (3.2), (3.6) следует:

$q({{t}_{1}}) = {{q}_{1}},\quad \frac{{dq}}{{dt}} = - {\kern 1pt} {\kern 1pt} Cq,\quad \frac{{dq}}{q} = - Cdt,\quad \ln \frac{q}{{{{q}_{1}}}} = - {\kern 1pt} {\kern 1pt} C\left( {t - {{t}_{1}}} \right)$
(3.7)
$q\left( t \right) = {{q}_{1}}\exp \left[ { - {\kern 1pt} C\left( {t - {{t}_{1}}} \right)} \right],\quad {{t}_{1}} \leqslant t \leqslant \left( {{{t}_{1}} + {{t}_{2}}} \right) = 75\,\,[час]$

Подставляя (3.7) в (3.1) при σ = 0, получаем

$\dot {p} = - Aq = - Aq_{1}^{{}}\left\{ {\exp \left[ { - {\kern 1pt} С(t - {{t}_{1}})} \right]} \right\}$

Отсюда уравнение кривой ползучести при разгрузке принимает такой вид:

(3.8)
$p(t) = {{p}_{1}} - Aq_{1}^{{}}\int\limits_{{{t}_{1}}}^t {\left\{ {\exp \left[ { - {\kern 1pt} С(t - {{t}_{1}})} \right]} \right\}dt} ,\quad {{t}_{1}} \leqslant t \leqslant \left( {{{t}_{1}} + {{t}_{2}}} \right)$
${{p}_{2}} = p({{t}_{2}} = 75\;час) = {{p}_{1}} - Aq_{1}^{{}}\int\limits_{{{t}_{1}}}^{{{t}_{2}}} {\left\{ {\exp \left[ { - {\kern 1pt} С(t - {{t}_{1}})} \right]} \right\}dt} $

На рис. 1 наряду с экспериментальной кривой при постоянном нагружении (кри-вая 1) приведена теоретическая кривая 5, соответствующая уравнению (3.5). Экспериментальной кривой ползучести при циклическом нагружении (кривая 2) соответствует теоретическая кривая 6, построенная согласно уравнениям (3.5), (3.8). Все теоретические кривые получены при $A = 7.2 \times {{10}^{{ - 6}}}$ (МПа ⋅ час)–1, $B = 0.178$ (час)–1, $C = 0.271$ (час)–1.

4. Описание результатов испытаний жаропрочного сплава ВЖ-159. В [6] npиведены результаты экспериментального исследования процесса высокотемпературной ползучести жаропрочного сплава ВЖ-159, которые были получены в НИИ механики Нижегородского государственного университета им. Н.И. Лобачевского.

Опыты были выполнены на цилиндрических образцах с длиной рабочей части l = 50 мм и диаметром d = 8 мм в условиях одноосного растяжения для различных уровней нормального напряжения $\sigma $ и температуры T. Результаты этих испытаний при T = 750°С приведены на рис. 2 в виде кривых ползучести 1 и 2 при постоянных напряжениях $\sigma = 350$ МПа и $\sigma = 450$ МПа соответственно. Аналогично результаты испытаний при T = 850°С приведены на рис. 3 в виде кривых ползучести 3 и 4 при постоянных напряжениях $\sigma = 120$ МПа и $\sigma = 150$ МПа соответственно.

Рис. 2.

Ползучесть сплава ВЖ-159 при T = 750°C

Рис. 3.

Ползучесть сплава ВЖ-159 при T = 850°C

Описание этих результатов приведено в [6] с помощью варианта теории Ю.Г. Коротких [7]. При описании экспериментальных данных в [6] применяются системы уравнений с использованием 11 материальных констант.

Авторы данной статьи для описания экспериментальных данных, приведенных на рис. 2 и 3, использовали вариант кинетической теории Ю.Н. Работнова [1]:

(4.1)
$\dot {p} = B{{\left[ {\frac{{\sigma \left( {1 + p} \right)}}{{1 - \omega }}} \right]}^{n}},\quad p\left( {t = 0} \right) = 0$
(4.2)
$\dot {\omega } = C{{\left[ {\frac{{\sigma \left( {1 + p} \right)}}{{1 - \omega }}} \right]}^{m}},\quad \omega \left( {t = 0} \right) = 0$

Кинетический параметр ω характеризует поврежденность материала, накапливаемую в процессе ползучести. Теоретические кривые ползучести, определенные с помощью уравнений (4.1), (4.2) при использовании значений материальных констант n = 9.782, $m = 4.542$, $B = 3.065 \times {{10}^{{ - 3}}}$ (МПа)n час–1, $C = 6.760 \times {{10}^{3}}$ (МПа)n час–1, позволяет хорошо описывать результаты испытаний при постоянных напряжениях, они приведены на рис. 2 (кривые 5 и 6) и 3 (кривые 7 и 8).

При T = 750°С в случае скачкообразного увеличения напряжения с $\sigma = 350$ МПа до $\sigma = 450$ МПа экспериментальные данные представлены на рис. 2 кривой 9. При T = 850°С в случае скачкообразного увеличения напряжения с $\sigma = 120$ МПа до $\sigma = 150$ МПа экспериментальные данные представлены кривой 11. В случае ступенчатого увеличения растягивающего напряжения при T = 850°С система уравнений (4.1)–(4.1) хорошо описывает экспериментальные данные (теоретическая кривая 12). При T = 750°С, однако, к моменту скачкообразного увеличения напряжения (${{t}_{1}} = 30$ мин) уровень поврежденности $\omega $ достиг достаточно большого значения, в связи с этим сразу после догрузки кривая ползучести стала почти вертикальной. Для исправления этого результата, т.е. для уменьшения скорости ползучести при $t > {{t}_{1}}$, следует заменить уравнение ползучести (4.1) на

(4.3)
$\dot {p} = B{{\left[ {\frac{{\sigma \left( {1 + p} \right)}}{{1 - \omega }}} \right]}^{n}} \cdot \exp q$
где q – второй (кроме ω) структурный параметр, удовлетворяющий кинетическому уравнению

(4.4)
$dq = cpd\sigma ,$$c = - 9.6$

Теоретическая кривая ползучести при T = 750°С после внезапного увеличения напряжения (кривая 10), соответствующая уравнениям (4.3), (4.4) хорошо соответствуют экспериментальной кривой 9.

5. Описание ползучести при прямом и обратном кручении. Подавляющее большинство экспериментальных исследований процесса ползучести проводится при растягивающих напряжениях, при этом осевое напряжение, естественно, не меняет знак. В отличие от таких испытаний в [8] приведены экспериментальные данные по ползучести тонкостенных трубчатых образцов при кручении, при этом касательное напряжение в процессе испытания меняет знак. В качестве материала для исследования был выбран дюралюмин Д16Т в состоянии поставки, материал был поставлен в виде прутка. Температура испытаний 150°С. Трубчатые образцы имели внешний диаметр $15.5 \pm 0.01$ мм, внутренний диаметр $14.0 \pm 0.01$ мм и рабочую часть длиной $100 \pm 0.1$ мм.

В закрученном тонкостенном образце напряженное состояние однородно и оно описывается касательным напряжением τ, действующим в плоскости поперечного сечения образца. В начальный момент времени (t = 0) к образцу прикладывался крутящий момент, создающий напряжение $\tau (t) = {{\tau }_{0}} = const$, при этом напряжении образец испытывался в течение ${{t}_{1}} = 50$ час. По прошествии 50 час образец разгружался и затем снова нагружался таким же по величине, но обратным по направлению крутящим моментом, создающим в образце напряжение $\tau (t) = - {{\tau }_{0}}$, при этом он испытывался еще 50 час.

На рис. 4 результаты испытаний при $\tau = \pm 101.1$ МПа представлены маркером ×, при $\tau = \pm 115.6$ МПа – ⚫ и ▲, при $\tau = \pm 130$ МПа – +, при $\tau = \pm 140$ МПа – ◼ и ◆ؚ, при $\tau = \pm 147.4$ МПа – ★. Для обратного кручения время отсчитывалось с момента, при котором напряжение меняло знак. В случае прямого кручения ${{\gamma }^{c}}$ означает сдвиговую деформацию ползучести, накопленную после приложения крутящего момента, а в случае обратного кручения γ c означает абсолютную величину сдвиговой деформации ползучести, накопленной после перемены знака крутящего момента.

Рис. 4

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

Таблица 1
№ кривой 1 2 3 4 5
|τ|, МПа 101.1 115.6 130.0 140.0 147.4

Для моделирования полученных результатов воспользуемся стандартной теорией упрочнения (γc – деформация ползучести при кручении, B, $\alpha $, $n$ – константы).

(5.1)
${{({{\gamma }^{c}})}^{\alpha }}{{\dot {\gamma }}^{c}} = B{{\tau }^{n}},\quad {\text{т}}{\text{.е}}{\text{.}}\quad {{\gamma }^{c}}(t) = {{[\left( {\alpha + 1} \right)B{{\tau }^{n}}t]}^{{1/(\alpha + 1)}}},\quad \alpha > 0,\quad n > 1$

Будем считать, что материальные константы $n$ и $\alpha $ не зависят от знака $\tau $, показатель степени $n$ выражается в виде отношения двух нечетных чисел, показатель $\alpha $ – в виде отношения четного числа к нечетному, а коэффициент $B$ принимает значения ${{B}_{1}}$ и ${{B}_{2}}$ при прямом и обратном кручении соответственно. В этом случае уравнение (5.1) можно записать в следующей форме:

${{({{\gamma }^{c}})}^{\alpha }}{{\dot {\gamma }}^{c}} = {{B}_{1}}{{\tau }^{n}}q{\text{,}}\quad q = \left\{ {\begin{array}{*{20}{l}} {1\quad при\quad \tau } \\ {{{B}_{2}}{\text{/}}{{B}_{{1\quad }}}при\quad - \tau } \end{array}} \right.$

На рис. 4 приведены теоретические кривые ползучести при $n = 799{\text{/}}99 \approx 8$, α = = 110/101 ≈ 1.1, ${{B}_{1}} = 0.8 \times {{10}^{{ - 14}}}$ (МПа)n час–1, ${{B}_{2}} = 2 \times {{10}^{{ - 14}}}$ (МПа)n час–1, при этом теоретические кривые ползучести при прямом кручении представлены сплошными линиями, а при обратном – штриховыми.

6. Заключение. Проведено моделирование результатов трех серий испытаний стержней на ползучесть при переменных нормальных или касательных напряжениях. Применение различных вариантов кинетической теории Ю.Н. Работнова при описании экспериментальных данных приводит к хорошему соответствию экспериментальных и теоретических результатов.

Соответствие номеров теоретических кривых и значений касательных напряжений представлено в табл. 1.

Работа выполнена при поддержке Российского научного фонда (проект № 19-19-00062).

Список литературы

  1. Работнов Ю.Н. Ползучесть элементов конструкций. М.: Наука, 1966. 752 с. (Перевод: Rabotnov Yu.N. Creep problems in structural members. Amsterdam, North Holland, 1969).

  2. Dushin Y.A., Gribov N.N., Ignatov V.A., Medvedev N.A. Structure-stable alloy for structural components of the high-temperature gas-cooled reactor with temperature at 950°C // Elsevier. 1991. V. 16. № 1/2. P. 317–326.

  3. Марголин Б.З., Душин Ю.А., Гуленко А.Г., Медведев Н.А. Ползучесть и разрушение структурно-стабильных материалов при нестационарном нагружении и объемном сжатии // Проблемы прочности. 1993. № 2. С. 13–25.

  4. Марголин Б.З., Швецова В.А. Влияние скорости деформирования на характер разрушения при статическом и циклическом нагружении. Сообщение 1. Формулировка общих подходов // Проблемы прочности. 1991. № 2. С. 3–14.

  5. Локощенко А.М. Ползучесть и длительная прочность металлов. М.: Физматлит, 2016. 504 с. (перевод: Lokoshchenko A.M. Creep and long-term strength of metals. CISP. CRC Press. Taylor & Francis Group. Boca. Raton. London. New York. 2018. 545 p.)

  6. Волков И.А., Игумнов Л.А., Казаков Д.А., Миронов А.А., Тарасов И.С., Шишулин Д.Н., Сметанин И.В. Модель поврежденной среды для описания длительной прочности конструкционных материалов (металлов и их сплавов) // Проблемы прочности и пластичности. 2017. Т. 79. № 3. С. 285–300.

  7. Волков И.А., Коротких Ю.Г. Уравнения состояния вязкоупругости сред с повреждениями. М.: ФИЗМАТЛИТ, 2008. 424 с.

  8. Наместников В.С. Прямое и обратное кручение в условиях ползучести. ПМТФ. 1960. № 1 С. 121–122.

Дополнительные материалы отсутствуют.