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

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2017. Т. 21, № 2. С. 326-361 ISSN: 2310-7081 (online), 1991-8615 (print) d http://doi.org/10

УДК 539.372

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

А. В. Хохлов

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

Аннотация

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

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

На основе сравнения обнаруженных свойств теоретических кривых с типичными свойствами экспериментальных кривых ползучести вязко-упругопластичных материалов при ступенчатых нагружениях изучены возможности линейной теории вязкоупругости по описанию различных

Статья

3 © Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования

Хохлов А. В. Анализ общих свойств кривых ползучести при циклических ступенчатых нагружениях, порождаемых линейной теорией наследственности // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2017. Т. 21, № 2. С. 326-361. doi: 10.14498/vsgtu1533. Сведения об авторе

Андрей Владимирович Хохлов& http://orcid.org/0000-0002-9212-2579

кандидат технических наук; старший научный сотрудник; лаб. упругости и пластичности;

e-mail: andrey-khokhlov@ya.ru

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

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

Получение: 14 марта 2017 г. / Исправление: 17 мая 2017 г. / Принятие: 12 июня 2017 г. / Публикация онлайн: 10 июля 2017 г.

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

Эксперименты на ползучесть при ступенчатых (с кусочно-постоянным напряжением) и циклических нагружениях — важные виды квазистатических испытаний, позволяющие обследовать разные аспекты поведения материала и уловить детали реализации многих эффектов [1-40], например, исследовать:

- зависимость кривых ползучести от уровня напряжения и предыстории нагружения, восстановление и остаточную деформацию при полной и частичной разгрузке [1-31, 36-38];

- влияние параметров начальной стадии нагружения и перестановки ступеней нагружения [2,4,5,12,29,33-39];

- скорость затухания памяти, эффекты, сопровождающие скачок напряжения вниз или его быстрое убывание за малый промежуток времени (dip tests, rate reversal in creep tests, non-monotonic creep behavior) [24,31,33-36];

- влияние скачков напряжения на длительную прочность и отклонение от правила линейного суммирования поврежденности (Miner's rule, linear

damage rule, cumulative damage theory) [3,4,11,14,21-23,29,30,40,41];

- процессы накопления или стабилизации деформации при циклических нагружениях (cyclic softening, cyclic hardening, ratcheting, shakedown), влияние задержек напряжения, коэффициента асимметрии (коэффициента амплитуд), частоты и других параметров цикла на скорость рэт-четинга [3,9,10,13,16,17,19,26,38,41-44];

- возможные проявления виброползучести (эффекты ускорения ползучести по сравнению с ползучестью при среднем или максимальном напряжении цикла, понижения деформации разрушения и долговечности) [3,11,16,19,26,29,45-47] и т. п.

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

Однако для расшифровки результатов испытаний (и их планирования), правильной интерпретации, верного выбора определяющего соотношения (ОС), его идентификации, настройки, грамотного применения и численной реализации необходимо точное знание общих свойств теоретических кривых ползучести при ступенчатом (и квазиступенчатом) нагружении (а также — кривых релаксации, деформирования и др.), порождаемых рассматриваемыми ОС, полученное в результате аналитического изучения их уравнений. Необходимо изучение влияния характеристик материальных функций и параметров программы нагружения (длительностей ступеней, скачков напряжения, скоростей нагружения на переходных участках, параметров цикла и частоты нагружения и т. п.) на общие качественные свойства теоретических КП (интервалы монотонности, выпуклости, точки экстремума и излома, скачки деформации и скорости деформации, асимптотику при t ^ те, отклонение от обычной КП при мгновенном нагружении до характерного значения напряжения), которые порождает применяемое ОС с произвольными материальными функциями (а не только для их конкретных видов). Такое исследование не было выполнено даже для одномерного линейного интегрального ОС вязкоупругости (или ОС «нелинейной теории наследственности» Ю. Н. Ра-ботнова [4]), хотя оно весьма полезно для уточнения арсенала его возможностей и области (не)применимости и выявления индикаторов (не)линейности поведения материалов по результатам испытаний.

Данная работа — очередной этап качественного анализа линейного ОС вязкоупругости

с произвольными функциями ползучести и релаксации П(£) и Я(Ь) и обобщающего его нелинейного ОС Ю. Н. Работнова [36,37,39,48-50]. Цель работы — аналитическое изучение общих свойств КП ОС (1) при циклических ступенчатых нагружениях (их интервалов монотонности и выпуклости, поведения последовательностей максимальных и минимальных значений деформации в каждом цикле и их полусуммы, отклонения от обычной КП при среднем напряжении, скорости рэтчетинга, условий циклической стабильности и т. п.) и их зависимости от параметров цикла нагружения и характеристик функции

ползучести.

В статье используются следующие сокращения и обозначения:

- ФР и ФП — функции релаксации и ползучести;

- КП — кривая ползучести;

- h(t) — функция Хевисайда;

- ¿(t) —дельта-функция Дирака;

- РеМ — регулярные модели (с ФП: П(0) = 0);

- СиМ— сингулярные модели (ФР содержит слагаемое ^¿(t), П > 0);

- y(t*) := y(t* + 0) — y(t* — 0) — скачок функции y(t) в точке t*;

- y(0) := y(0+) — предел y(t) справа в точке t = 0.

1. Некоторые свойства линейного определяющего соотношения вязкоупругости. ОС (1) описывает изотермические одноосные процессы в структурно стабильных вязкоупругопластичных материалах, характеризуемые историями напряжения a(t) и деформации e(t) в точке тела (входные процессы предполагаются кусочно-непрерывными и кусочно-гладкими при t ^ 0). Операторами вида (1) часто задаются и зависимости между историями нагрузок и перемещений при испытаниях образцов на растяжение, кручение, изгиб, индентирование и т. п. Функции ползучести и релаксации (ФП и ФР) в (1) предполагаются положительными и дифференцируемыми на (0; те), n(t) — возрастающей и выпуклой вверх на (0; те) [37], а R(t) — убывающей и выпуклой вниз на (0; те), ФР может иметь интегрируемую особенность или ¿-сингулярность в точке t = 0. Из монотонности и положительности ФП и ФР на (0; те) следует, что в точке t = 0 существуют пределы справа П(0+) = infn(t) ^ 0 и R(0+) = sup R(t) > 0 (R(0+) = +те, если ФР не ограничена) и предел Я(+те) = inf R(t) ^ 0.

Операторы (1) определены на линейном пространстве кусочно-непрерывных и кусочно-гладких функций. Они представимы в виде

y(t) = yr (t)+ ys(t), где yr(t) — регулярная часть (непрерывная функция, yr(0) = 0),

ys

i= 1

— ступенчатая функция

и оператор (1) переводит такой процесс в деформацию ^ п

е(£) = П(£ - т)<тг(т)^т + ст(0)П(£) + У" <т(^)П(£ - - ¿¿), £ ^ 0. ./о 1

Если П(0) = 0 (модель регулярна), то Е(0) = 1/П(0) < те, и на линеале непрерывных кусочно-гладких при £ ^ 0 функций операторы (1) представимы в виде

= Д(0)е(г) + - т)ф)^т,

е(£) = П(0)ст(£) + 1П(£ - т)ст(т)^т, £ ^ 0.

о

Операторы (1) взаимно обратны и поэтому ФП и ФР связаны уравнением

[ n(t - т)R(t)dr = t, или f n(t - т)R(t)dr + n(0)R(t) = 1, t> 0. (3) J 0 Jo

Зная ФР, можно найти ФП из уравнения (3), и наоборот. При П(0) = 0 уравнения (3) и (2) — уравнения Вольтерры второго рода с ограниченными (если П(0+) < те) ядрами, и поэтому они однозначно разрешимы в пространствах Li [0,b]. Случай П(0) = 0 приводит к уравнению Вольтерры первого рода, некорректным задачам, нерегулярным моделям с особенностью в нуле у функции и кривых релаксации и касательного модуля диаграмм деформирования и т. п. [37,48].

Свойства семейств основных теоретических квазистатических кривых (диаграмм деформирования при постоянных скоростях деформации или на-гружения, кривых ползучести при ступенчатом нагружении, кривых ползучести и релаксации с произвольной начальной стадией нагружения и др.), порождаемых ОС (1) с произвольной ФП, необходимые математические и феноменологические ограничения на ФП и ФР проанализированы в цикле работ [36, 37, 39, 48, 49] и др. Анализ, в частности, показал, что среди моделей, описываемых ОС (1) с различными ФР и ФП, необходимо выделять три основных класса, поскольку качественные свойства базовых теоретических кривых моделей этих классов (а также особенности постановки и решения краевых задач) заметно отличаются:

1) регулярные модели (РеМ) — у которых П(0) = 0 (тогда мгновенный модуль E = R(0) = 1/П(0) диаграмм деформирования с постоянной скоростью конечен, П(0)/П(0) = — R(0)/R(0) [37], а ОС (1) и первое уравнение (3) сводятся к уравнениям Вольтерры второго рода (2) и (3));

2) сингулярные (СиМ) — модели с ФР, содержащей слагаемое n^(t), П > 0 (ФР R = n^(t) задает ньютоновскую жидкость с ОС а = цё и входит слагаемым в ФР половины реологических моделей из пружин и демпферов), тогда П(0) = 0 и П(0) = п-1;

3) нерегулярные модели с неограниченной ФР, не содержащей слагаемого n^(t), но имеющей интегрируемую особенность в точке t = 0 (R(0+) = = +те).

При малых временах (и больших скоростях деформации) РеМ ведут себя как твердые тела, а СиМ — как жидкости. Третий класс занимает промежуточное положение между первыми двумя. К нему относится, например, ФР

R(t) = At-U, u е (0; 1), A > 0,

задающая так называемый «фрактальный» элемент «фрактальных» моделей ("fractional models") с оператором дробного дифференцирования; соответствующая (в силу (3)) ФП имеет вид n(t) = A-1C(u) tu и обладает не только свойством П(0) = 0, как и СиМ, но и свойством П(0) = те, переходным к П(0) = 0, характеризующему РеМ.

Все структурные реологические модели, собранные из линейных пружин и демпферов посредством последовательных и параллельных соединений, описываются ОС (1). Схемы и названия всех двух-, трех- и четырехзвенных моделей (в терминологии нет единства) приведены в [37] на рис. 1. Можно

доказать, что множество всех несократимых n-звенных моделей распадается ровно на два класса эквивалентности: РеМ-n и СиМ-n (структурно различные модели мы называем эквивалентными, если они задаются одинаковыми семействами ФП или ФР). В частности:

1) эквивалентны трехзвенные РеМ Пойнтинга—Томсона и Кельвина [37, рис. 1, а];

2) все четыре РеМ-4 [37, рис. 1, в] эквивалентны модели стандартного тела (последовательному соединению моделей Максвелла и Фойгта, т. е. РеМ-2 и СиМ-2);

3) все РеМ-2к эквивалентны параллельному соединению к моделей Максвелла с разными временами релаксации;

4) все СиМ-2к эквивалентны последовательному соединению k моделей Фойгта с разными временами ползучести (retardation time);

5) РеМ-(2к + 1) получается из СиМ-2к последовательным присоединением упругого элемента (РеМ-1), а СиМ-(2к + 1) — из РеМ-2к параллельным подключением вязкого элемента (СиМ-1).

Например, семейство ФП

n

порождает все РеМ-4 при y е (0;в), а,в > 0, а при a = 0 — РеМ-3 (модель Кельвина). Так как П(0) = в — Y, ФП (4) порождает СиМ, когда y = в: при Ав = 0 — ньютоновскую жидкость, при a = 0 — модель Фойгта, при a > 0 получаются (все) СиМ-3. При y = 0 (4) дает модель Максвелла. Эти классические модели будут использованы для иллюстрации общих свойств циклических кривых ползучести ОС (1).

2. Кривые обратной ползучести и требование выпуклости функции ползучести. Кривые ползучести (КП), порождаемые ОС (1) при на-гружении o(t) = Oh(t), имеют вид

e

КП (5) возрастает по t (так как n(t) возрастает) и линейно зависит от уровня напряжений.

Отклик ОС (1) на прямоугольный импульс нагрузки o(t) = <r[h(t) — h(t — T)], O > 0, длительностью T > 0, имеет вид

e(t) = OS(t;T), S(t;T) := n(t)h(t) — n(t — T)h(t — T), t ^ 0.

В точке t = T e(t) имеет место скачок — <гП(0), а ¿(t) —скачок —<гП(0). На интервале t > T уравнение кривой обратной ползучести (КОП) имеет вид

e

У всех реономных материалов (кроме вязких жидкостей) после снятия нагрузки наблюдается постепенное убывание деформации до некоторого уровня

¿те ^ 0 (для сетчатых полимеров в высокоэластичном состоянии, как правило, = 0). Это явление называется последействием, упругим восстановлением, обратной ползучестью [3-8, 20-29]. Из требования (нестрогого) убывания КОП (6) (с любым T) следует невозрастание n(t) [37]. Поэтому на ФП в ОС (1) следует накладывать ограничение: n(t) не имеет участков выпуклости вниз. Можно доказать, что ограничение n(t) ^ 0 не является следствием остальных ограничений на ФР и ФП (и тождества (3)): существует гладкая убывающая ФР с R(t) ^ 0, такая, что соответствующая ФП строго монотонна, но имеет участок с n(t) > 0 (хотя на всем этом участке R(t) > 0).

Итак, n(t) убывает на луче t > 0 и n(t) > 0, значит, существует предел v := П(те) ^ 0. Так как КОП (6) убывает и положительна (ограничена снизу), она имеет при t ^ те предел ^ 0, и = OS(ro,T). Легко доказать, что S(ro,T) = vT [37], т. е. в случае v > 0 после полной разгрузки (бесконечно долго) сохраняется остаточная деформация = OvT. Это означает, что при v > 0 память интегрального оператора П, задающего ОС (1), не затухает.

Если ФП ограничена, то v = 0 и = 0. Нетрудно доказать, что ФП всех моделей СиМ-2к и РеМ-(2к +1) ограничены, а у РеМ-2к и СиМ-(2к +1), к е N, v > 0. Может быть v = 0 и для неограниченной ФП, например, для степенной ФП

П = ctu, u е (0; 1), c > 0.

У классических моделей (4) v = a и

S(t; T) = aT + ce-Ai, c := Y(eAT — 1).

В частности, v > 0 при a > 0 (для РеМ-2 и РеМ-4); у модели Максвелла и ньютоновской вязкой жидкости (СиМ-1) y = 0 и

e(t) = O aT = const

при t > T, т. е. упругое восстановление не происходит. У моделей Фойгта (СиМ-2) и Кельвина (РеМ-3) a = 0 и поэтому КОП (6) имеет вид

e(t) = <rce-Ai

и = 0. Случай y < 0 в (4) приводит к нарушению ограничения n(t) ^ 0, что влечет за собой возрастание кривой обратной ползучести (противоречие с экспериментальными данными).

Из ограничения n(t) ^ 0 следует, что КП (5) всегда выпуклы вверх, и потому ОС (1) не способно описывать ускоряющуюся ползучесть (третью стадию).

Отметим, что КОП (6) с O > 0 выпукла вниз для любых t > T > 0 (как в испытаниях материалов) тогда и только тогда, когда n(t) возрастает (S(t; T) = n(t) — n(t — T) > 0). Это условие выполняется, в частности, для всех ФП, представимых суммами экспонент с отрицательными показателями (моделей СиМ-2к, РеМ-(2к + 1) и др.) и суммами функций ctu.

3. Кривые ползучести при произвольных ступенчатых нагруже-

ниях. Кусочно-постоянной программе нагружения с n ступенями

n— 1

o

i=i

(полагаем, что ¿0 = 0 и ¿г > ^1) ОС (1) ставит в соответствие деформацию [36]

га—1

еф =

г=1

или

г

е(г) = ^1П(г) + ^ (кп(г - гк), г е (гг; ¿г+1), к=1

где Гг := гг+1 -гг — скачки напряжения в точках ¿г, г0 := 0. Скачки в(г) и ¿т(г) в точках ¿г выражаются формулами [36]:

= (Г ¿П(0), ^ = П(0)( + П(0)Гг = 11(0) (Г..

Они не зависят от момента времени ¿г и предыстории нагружения, а только от скачка напряжения гг+1 - гг; скачок деформации при разгрузке на Дг всегда равен по модулю скачку при догрузке на Дг. Эти свойства можно использовать как индикатор применимости линейного ОС (1). Линейное ОС вязко-упругости (1) с произвольной ФП неспособно моделировать явление «дрейфа мгновенно-упругой деформации вследствие деформации ползучести» [51-54, 31], проявляющееся, в частности, в том, что скачок деформации при разгрузке отличается по модулю от скачка при догрузке (с тем же модулем скачка напряжения), так как зависит от предыстории нагружения. В [50] установлено, что нелинейное ОС наследственности Ю. Н. Работнова [4] способно описывать явление дрейфа.

При г ^ те из (7) (в силу существования предела Б(те, Т) = -Т) получим

е

т

Рт := -^Гг

г=1

Перестановка ступеней нагружения с г < п не влияет на рп-1 [39] и асимптотику (8).

Если гп = 0, то е^) ^ рп-1 при г ^ те, т. е. КП (7) имеет горизонтальную асимптоту е = рп-1, а рп-1 приобретает смысл пластической (необратимой) деформации.

Если гп = 0, но V = 0, то рп-1 = 0 и поэтому при г ^ те имеем

е(г) - ГпП^ - ¿п-1) ^ 0, е^) - ГпПед ^ 0, п(г) - п(г - ^-1) ^ 0.

В частности, V = 0, если П^) ограничена (например, у всех СиМ-2п и РеМ-(2п + 1)); в этом случае рп-1 = 0 и е(г) ^ при г ^ те. Если П^) не

ограничена и гп = 0, то

е(г) ~ ГпП^ - ¿п-1),

и | неограниченно возрастает.

Если V = 0 (как у РеМ-2к и СиМ-(2к + 1), к е М) и Гп = 0, то е(г) ^ те при г ^ те,

^(П^) - П(г - ¿п-1)) ^ ^-¿п- = 0,

¿(t) - a„n(t - in-i) ^ Pn-i

(значение pn-i зависит от программы {(ti, ai)} и может быть любым, но если все ai ^ 0 и хотя бы одно a^ > 0, то и pn-i > 0), а ¿(t) — ann(t), вообще говоря, не стремится к нулю при t ^ те. Отметим, что если все ai ^ 0 для i < n, то отклонение ¿(t) — ann(t — tn-i ) убывает на луче t > tn-i, т. к. все слагаемые в (7), кроме последнего, убывают по t при t > ti.

В статье [36] изучены и проиллюстрированы общие свойства КП (7), в частности, найдены условия существования интервала убывания и точки минимума у кривых ползучести с неполной разгрузкой. В [39] исследовано влияние перестановки ступеней нагружения и обнаружено свойство асимптотической коммутативности ОС (1) при ступенчатых нагружениях. Ниже будет изучено влияние материального параметра v на скорость накопления пластической деформации при циклических нагружениях.

4. Деформация при циклическом импульсном нагружении с полной разгрузкой. Пусть программа нагружения состоит из одинаковых прямоугольных «отнулевых» циклов с отдыхом между ними (ti = iT, a2k-i = a, a2k = 0, k = 1, 2,... ):

те те

a(t) = a ^ (—1)ih(t — ti), или a(t) = ^ a[h(t — t2k) — h(t — t2fc+i)]- (10) i=0 k=0

Тогда деформация определяется по формуле

те

¿(t) = a ^ (—1)in(t — ti)h(t — ti), или

¿(t) = a ^ [n(t — t2k)h(t — t2k) — n(t — t2k+i)h(t — t2k+i)].

k=0

Ряды (10), (11) сходятся равномерно на любом отрезке, ибо на нем отличны от нуля лишь конечное число членов ряда. Очевидно,

0 < ¿

при всех t > 0 в силу возрастания ФП. Скачкам напряжения ai := (—1)ia в точках ti, i = 0,1,..., соответствуют скачки ¿(t) и ¿(t):

li = ain(0) = (—1)in(0)a, ¿i = n(0)a + n(0)ai = (—1)iii(0)a,

т. к. в данном случае ai = 0.

При t G (t2i-i; t2i), i = 1, 2,... (i — номер цикла) из (11) имеем

ii ¿(t) = a ^ [n(t — t2k-2) — n(t — t2k-i)] = a ^ S(t — t2k-2), k=i k=i т. е. ¿(t) — сумма сдвигов кривой обратной ползучести (6) вдоль оси времени (S(t) — краткое обозначение для S (t, T ) из (6)). Каждое слагаемое S (t—t2k-2) положительно (если n(t) = const) и убывает на луче t > t2k-i, а при t ^ те

S(t — i2fc—2) ^ vT, где v := П(те). Значит, e(t) > 0 и монотонно убывает на (¿2г—1; t2i) (лишь для упругого и вязкого элемента и модели Максвелла

e(t) = const = ¿a0T при t G (t2i-1; t2i),

т. к. S(t—t2k-2) = aT; для упругого элемента S(t, T) = 0 при t > T и e(t) = 0). Нижняя грань e(t) на интервале (t2i-i; t2i) (наименьшее значение в ¿-том цикле):

£ = inf e(t)= e(t2i — 0) = ffV S(t2i — t2fc-2) = ffV S(2(i — k + 1)T),

fe-ife) fc=i fc=i

т. е.

i

^ = 0 £ S

j=i

Так как S(2jT) > 0, то ei > 0 (только для РеМ-1 ei = 0) и последовательность ei возрастает. Так как S(2jT) ^ S(те) = vT, справедлива оценка снизу £i ^ ¿0р, где р := vT остаточная деформация от единичного импульса напряжения (равенство = ¿0р достигается лишь для вязкого элемента и модели Максвелла, для которых v = a и S(t — t2k-2) = aT = const).

При t G (t2i-2; t2i-1), i = 1,..., из (11) имеем

i-1 i-1 e(t) = ¿n(t) — 0 £ [n(t — t2k-1) — n(t — t2k)] = ¿n(t) — 0 £ S(t — t2k-1), k=1 k=1 каждое слагаемое S(t — t2k-1) положительно и убывает, следовательно, e(t) монотонно возрастает. Наибольшее значение деформации на интервале (t2i-2; t2i-1) выражается через предел в его правом конце:

£i = sup e(t) = e(t2i-1 — 0) =

= ¿n(t2i-1) — [n(t2i-1 — t2fc-1) — n(t2i-1 — t2fc)],

(t2i-2 ;t2i-1) ^ 1

k=1

т. е.

i

^i = Ш(0) + 0 £ [n((2j — 1)T) — n((2j — 2)T)] =

j=1 i

= 0П(0)+ 0 £ S((2j — 1)T). (13) j=1

Последовательность (ег} возрастает (ибо ФП возрастает и S((2j — 1)T) > 0). Так как S((2j — 1)T) > S(те) = vT, справедлива оценка снизу

ёг ^ стП(0) + ¿стр.

Из (12) и (13) с учетом (6) находим

ёг + ё = стП(0) + Ш(2гТ) — стП(0) = стП(2гТ),

т. е. последовательность средних деформаций ег := 0.5(ёг+ёг) целиком лежит на КП ё = 0.5<rn(i) для среднего напряжения цикла 0.5ст. Поэтому последовательность (ег} ограничена тогда и только тогда, когда ФП ограничена. В этом случае v = 0, р = 0 и ё ^ 0.5стП(те) при i ^ те.

Если же v > 0, то р = vT > 0, пластическая деформация (9) за i циклов равна ¿стр, и

ёг ^ ¿стр, ё ^ стП

т. е. деформация неограниченно возрастает, пластическая деформация неограниченно накапливается со скоростью р за цикл, и ОС моделирует рэтче-тинг.

В силу формул (12) и (13) размах «колебаний» e(t) на ¿-том цикле выражается формулой

2г г

ё — ёг = стП

m=1 fc=1

A

Так как ФП выпукла вверх, A(t, T) ^ 0 при всех t ^ 2T и поэтому последовательность ёг — ёг возрастает (не убывает). Лишь у модели Максвелла (с П = at + в, а, в > 0) и упругого и вязкого элементов (с а = 0 или в = 0) A(t, T) = 0 и поэтому ёг—ёг = стП(0) = ств = const (не зависит от T). При этом ёг = ¿CTaT, ёг = ств + ¿aaT и обе последовательности растут неограниченно при a > 0 (и возрастают по T).

Пример 1. Для моделей Кельвина

ёг = —ст £ Y(e-2fcTA — e-(2k-1)TA) = fc=i

= —ст7(е-лт — 1)e-AT (1 — e-2ATi)(1 — e-2AT )-1 =

= ст7(1 — e-2iAT )(eAT + 1)-1,

ёг = ст(в — Y) + CTY(1 — e-2iAT )(1 + e-AT )-1. Очевидно, ёг и ёг связаны линейной зависимостью

ёг = eAT ёг + ст

(для модели Фойгта в — 7 = 0), обе последовательности возрастают и ограничены, а при i ^ те сходятся к пределам

е := a7(eAT + 1)-1, е := (1 + e-AT)-1 + е(в — y).

Величины е, е,, е и е зависят от отношения длительности полуцикла T к времени ретардации модели 1/А: предел е убывает по T (по AT), а е и (е — е) возрастают; при T ^ 0 е ^ 0.5е7, е ^ ев — 0.5е7, е — е ^ е (в — Y), а при AT ^ те е ^ 0, е ^ ев, е — е ^ ев.

Модель стандартного тела (РеМ-4) получается из модели Кельвина (РеМ-3) последовательным подключением вязкого элемента (П = at). При последовательном соединении моделей их ФП складываются, поэтому складываются и функции (16) и последовательности {е}, {е,}, {е — е,}. Поэтому для РеМ-4

v = a > 0, е — ie aT = eAT(е — ieaT) + e (в — y),

или

е = еЛте + (1 — eAT )eaTi + е (в — y).

Примечательно, что е, ^ те, но последовательность {е — е,} ограничена (ибо для вязкого элемента е — = 0), не зависит от a и совпадает с {е — е,} для РеМ-3:

ее — е = е(в — y) + 07(1 — e-2^ )[(1 + e-^ )-1 — (e^ + 1)-1] =

= е (в — 7) + ёу^ — 1)(e^ + 1)-1(1 — e-2^).

При i ^ те е — ^ е(в — 7) + <JY(e^ — 1)(e^ + 1)-1 = L; для любого i е(в — 7) < е — < L < ев; предел L возрастает с ростом AT, при AT ^ 0 L ^ е(в — 7), а при AT ^ те L ^ ев.

Пример 2. Для модели со степенной ФП П = ctu, u G (0; 1), имеем

A(t, T) = ctu — 2c(t — T )u + c(t — 2T )u =

= ctu[1 — 2

При t ^ те, очевидно, A(t, T) = —2cu(1 — u)T2tu-2 + O(tu-3). Последовательности

í í = <rcTu ^ [(2k)u — (2k — 1)u], е = <rcTu ^ [(2k — 1)u — (2k — 2)u] fc=1 fc=1

не ограничены (ряды с неотрицательными членами расходится по признаку сравнения); для любого i функции £¿(T), £¿(T) возрастают по T и е(0+) = 0,

е(0+) = 0.

На рис. 11 приведены кривые ползучести (11) четырех моделей для на-гружения с T = 5, е = 1:

ХВ онлайн-версии статьи все рисунки выполнены в цвете.

О

£

3

2

1

О

10

20

30

40 г

Рис. 1.

1) КП модели со степенной ФП П = сЬи, с = 1, и = 0.5 (черные КП);

2) КП РеМ-4 с П = аЬ+в—7е—Л4, Л = 0.1 (время ретардации т = 1/Л = 10), а = 0.1, в = 1.5, 7 = 1 (голубые КП);

3) КП модели Фойгта с П = 7 — 7в-Л4 с Л = 0.1, 7 = 1 (красные);

4) КП модели Максвелла (синие) с П = аЬ + в, а = 0.1, в = 0.5 (время релаксации т = в/а = 5).

РеМ-4 и модель Максвелла (РеМ-2) регулярны и поэтому их КП имеют в точках Ь = кТ разрывы первого рода со скачками ±<гП(0); свойство V > 0 вызывает накопление пластической деформации. У модели Фойгта и степенной модели V = 0 и П(0) = 0 и поэтому КП непрерывны при всех Ь > 0. Штрихпунктирные кривые 1-4 — обычные КП этих моделей при постоянном напряжении а = 0.5а (среднем за цикл). При достаточно больших Ь КП РеМ-4 и РеМ-2 лежат выше КП 1 степенной модели. КП вязкого элемента (не показана) не имеет разрывов и совпадает с КП модели Максвелла на интервалах Ь € (Ь2г— 1; Ь2г) (где КП — горизонтальные отрезки). Аналогично, КП СиМ-3 не имеют разрывов и отличаются от КП РеМ-4 только на интервалах Ь € (Ь2г—2;Ь2г—^ г = 1, ...,т (они получаются сдвигом дуги КП РеМ-4 вниз на в, а убывающие участки КП на интервалах Ь € (Ь2г—2; Ь2^—1) у СиМ-3 и РеМ-4 совпадают). Для модели Фойгта и степенной модели дополнительно приведены КП с удвоенной длительностью цикла, т. е. Т = 10 (штриховые кривые 5, б).

Исследуем условия ограниченности (существования пределов) последовательностей максимальных и минимальных деформаций (12) и (13), в частности докажем, что последовательность {ёг — ег} (см. (15)) ограничена для любой допустимой ФП. Будем использовать краткое обозначение Д(Ь) для функции (16).

Лемма 1. Пусть Т > 0, а > 0, а ФП П(Ь) дифференцируема, положительна, возрастает и выпукла вверх при Ь > 0. Тогда:

1) ряд ^Д(кТ) сходится, а его сумма выражается формулой

(£ = 0 лишь для упругого и вязкого элементов и модели Максвелла);

£ = vT + П

2) ряды ^3=1 Д((2? + 1)Т) и ^3=1 Д(2^Т) сходятся, их частичные суммы связаны формулой

БП + БП = П

суммы рядов £1, £2 —формулой

£1 + £2 = vT + П(0) — П(Т),

и справедливы неравенства БП ^ БП ^ 0 и £2 ^ £1 ^ 0. Доказательство.

1) По формуле

Д(кТ) = П((к — 2)Т) — 2П((к — 1)Т) + П(кТ),

поэтому каждое слагаемое П(кТ) входит ровно в три члена ряда, а при их суммировании пропадает. Значит частичная сумма ряда содержит лишь две пары крайних слагаемых:

Бп = П(0) — П(Т) + П(пТ) — П((п — 1)Т).

По теореме Лагранжа

П(пТ) — П((п — 1)Т) = П(£)Т,

где £ € ((п — 1)Т;пТ). При п ^ те имеем £ ^ те, П(£) ^ V (в силу убывания П(Ь)) и Бп ^ vT + П(0) — П(Т). Оценка £ ^ 0 следует из Д(Ь,Т) ^ 0 при Ь ^ 2Т, а £ ^ П(0) — П(Т) — из V ^ 0 (напомним, что П(Ь) ^ 0, и поэтому П(те) ^ 0).

2) Последовательность частичных сумм каждого из двух подрядов сходящегося знакопостоянного ряда £Д

Д

и каждое слагаемое с четным индексом, кроме первого и последнего, входит в два члена ряда, а слагаемое с нечетным — в один. Поэтому

п

БП = ^ Д(2;Т) = П(0) + П(2пТ)+

3 = 1

п— 1 п

+ 2 ^ П(2;Т) — 2 ^ П((2; — 1)Т) =

3=1 3=1

2п—1

= П(0) + П(2пТ) +2 ^ (—1)тП(тТ).

т=1

Аналогично

п 2п

БП = £ Д((2; + 1)Т) = = П(Т) + П((2п + 1)Т) — 2 £ (—1)тП(тТ),

3=1 т=2

и формула (18) доказана.

Итак, обе последовательности (е»}, (е»} возрастают, но (а — е»} ограничена (и имеет предел аП(0) — аХ2 в силу леммы 1 и формулы (15)). Поэтому последовательности (а} и (е»} могут быть ограничены только одновременно. Так как е» + а = аП(2гТ), обе последовательности (е»}, (е} ограничены (и сходятся) тогда и только тогда, когда ФП ограничена (в этом случае V = 0, р = 0 и е» := 0.5(а + е») ^ 0.5аП(те)). В этом случае ряды в формулах (13) и (12) сходятся, а пределы последовательностей (а»} при г ^ те выражаются через их суммы $1 и формулами

а = аП(0) + а е = а82.

Таким образом, доказана следующая теорема о свойствах отклика (11) на 2Т-периодическую программу нагружения (10) из прямоугольных импульсов.

Теорема 1. Пусть Т > 0, а > 0, а ФП П(£) дифференцируема, положительна, возрастает и выпукла вверх при £ > 0. Тогда кривая ползучести (11) обладает следующими свойствами:

1) на интервалах £ € (2(г — 1)Т, (2г — 1)Т), г = 1,..., кривая ползучести (11) возрастает и, выпукла вверх, а на интервалах ((2г — 1)Т, 2гТ) е(£) убывает (а если П(£) ^ 0, то и выпукла вниз);

2) максимальное и минимальное значения а и е» деформации е(£) в г-том цикле выражаются формулами (13) и (12); а ^ е» ^ 0, последовательности (а} и (е»} возрастают, е» := 0.5(а + е») = 0.5аП(2гТ), и справедливы оценки снизу (14);

3) размах колебаний е(£) на г-том цикле а — е» выражается формулой (15), последовательность а — е» возрастает, ограничена и сходится; ее предел равен аП(0) — аХ2;

4) последовательности (а»} и (е»} ограничены тогда и только тогда, когда ФП ограничена; в этом случае их пределы выражаются формулами:

а = аП(0) + = аП(те) — а$2, е =

(где в1 = ^((2; — 1)Т), «2 = 5(2,Т));

е» ^ 0.5аП(те), в1 + $2 = П(те) — П(0);

для деформации (11) верна двусторонняя оценка 0 ^ е(£) ^ а при £ > 0, приращение пластической деформации за цикл ар = <гоТ равно нулю, а ОС (1) моделирует циклическое упрочнение материала;

5) если V = П(те) > 0, то р = vT > 0 и справедливы оценки

е ^ гар, а» ^ аП(0) + гар и е» ^ 0.5аП(0) + гар;

в этом случае пластическая деформация после г циклов равна гар, а ОС (1) моделирует рэтчетинг с постоянной скоростью накопления пластической деформации ар за цикл (а средняя полная деформация е» = 0.5аП(2гТ) растет еще быстрее);

6) при Т ^ 0 (для любого фиксированного г ^ 1) а ^ аП(0), е» ^ 0, а если ФП ограниченна, то е ^ аП(0) и е ^ 0.

Замечание 1. При последовательном соединении моделей их ФП складываются, поэтому суммируются и функции Д(Ь,Т), и последовательности {ёг}, {ёг}, {ёг — ёг}. Линейная комбинация (с положительными коэффициентами) ФП с V = 0 наследует свойство V = 0; наличие хотя бы одного слагаемого с V > 0 обеспечивает свойство V > 0 у суммы ФП.

Замечание 2. С ростом Т величины ёг и ёг не обязаны всегда возрастать или убывать. Достаточные условия монотонности:

1) если функция ЬП(Ь) возрастает, т. е. ЬП(Ь) + П(Ь) ^ 0 (это выполняется, в частности, для степенных ФП), то ¿¿(Т) и ёг(Т) возрастают, ибо производная по Т каждого члена Б(2^'Т) = П(2^Т) — П((2^ — 1)Т) суммы (13) и (12) неотрицательна;

2) если ЬП(Ь) убывает, то ёг(Т) и ёг(Т) убывают.

У моделей Кельвина и Фойгта ЬП(Ь) = 7ЛЬв—Л сначала возрастает, а потом убывает, ёг(Т) убывает, ёг(Т) возрастает.

5. Произвольная программа двухступенчатого циклического на-гружения. В силу п. 2 теоремы 1 последовательность средних деформаций ёг целиком лежит на КП ё = 0.5аП(Ь) при среднем напряжении цикла 0.5а (совпадает с ее значениями в точках Ь = 2гТ). Это связано с тем, что программа импульсного нагружения (10) разлагается в сумму постоянного (среднего) напряжения 0.5а и симметричного нагружения с тем же периодом и амплитудой 0.5<г , т. е.

а+(Ь, Т) = 0.5а + 0.5аа+(Ь, Т), и поэтому (в силу линейности ОС (1)) отклик представляется в форме

ё+ (Ь, Т) = 0.5аП(Ь) + 0.5аё+(Ь, Т),

где ё+(Ь,Т) —отклик на симметричное периодическое нагружение а+(Ь,Т) с единичной амплитудой.

Произвольная программа (равнодольного) двухступенчатого циклического нагружения

те

а

г=1

когда в моменты времени Ьг = гТ (Т —длительность полуцикла) происходит перескок с напряжения а1 на а2 и обратно, тоже разлагается в сумму постоянного напряжения ат = 0.5(а1 + а2) и симметричного нагружения с тем же периодом 2Т, амплитудой <г = ±0.5а и тем же модулем скачков напряжения. В силу линейности ОС (1) отклик на нее представим в виде суммы обычной КП и отклика на симметричное циклическое нагружение 0.5аа+(Ь,Т):

ё(Ь) = атП(Ь) + 0.5аё—(Ь,Т).

Отклик 0.5аё+(Ь, Т) равен отклонению ё(Ь) от КП атП(Ь), и анализ его свойств позволит исследовать, в частности, возможности моделирования с помощью ОС (1) некоторых эффектов наложения циклического возмущения на постоянное нагружение, например, ускорения или замедления ползучести по

сравнению с ползучестью при среднем напряжении (а не при максимальном, что принято называть виброползучестью) [3, 9,11,13,16,17,19, 26, 29, 45-47]. Нетрудно заметить, что, в силу линейности и инвариантности оператора (1) относительно сдвигов вдоль оси времени отклик ё+(Ь, Т) выражается через линейные комбинации функции (16) и ее сдвигов по времени, и поэтому для исследования деформации в случае произвольного циклического нагружения (19) следует детально изучить свойства функции Д(Ь, Т).

6. Функция влияния симметричного прямоугольного импульса растяжения—сжатия. Исследуем отклик ОС (1) на трехступенчатую программу нагружения

а«(Ь) = а[Ь(Ь) — И(Ь — Т)] — а[И(Ь — Т) — И(Ь — 2Т)], а> 0, Т> 0

(а(Ь) = 0 при Ь > 2Т), состоящую из двух прямоугольных импульсов растяжения и сжатия одинаковой продолжительности и последующей полной разгрузки. Оператор (1) переводит этот процесс в

ё(Ь) = а Д(Ь; Т),

где

Д(Ь; Т) := Б(Ь; Т) — Б(Ь — Т; Т),

т.е.

Д

Очевидно, Д(Ь, Т) = 0 при Ь < 0 (и поэтому Д(Ь — т, Т)И(Ь — т) = Д(Ь — т, Т)). В дальнейшем (при фиксированном Т) будем использовать сокращенное обозначение Д(Ь).

Если Ь € (0; Т), то Д(Ь) = П(Ь) —положительная возрастающая выпуклая вверх функция.

Если Ь € (Т; 2Т), то Д(Ь) = П(Ь) — 2П(Ь — Т) = Б(Ь; Т) — П(Ь — Т) — убывающая функция, причем Д(2Т — 0) = П(2Т) — 2П(Т) < 0, так как П(Ь) выпукла вверх. Скачок Д(Ь, Т) в точке Ь = Т равен —2П(0), а в точках Ь = 0 и Ь = 2Т равен П(0) (для сингулярных моделей П(0) =0 и поэтому Д(Ь, Т) непрерывна). Скачки производной Д(Ь) в этих трех точках равны 11(0), —211(0) и 11(0). Если П(Ь) ^ 0, то £(Ь; Т) ^ 0 и Д(Ь) = Б(Ь; Т) — П(Ь — Т) ^ 0, т. е. Д(Ь) выпукла вниз на интервале (Т; 2Т).

Если Ь ^ 2Т, то Д(Ь) имеет вид (16): Д(Ь,Т) := П(Ь) — 2П(Ь — Т)+П(Ь — 2Т). Для удобства Д(Ь) доопределена в точке Ь = 2Т по непрерывности справа: Д(2Т) := Д(2Т + 0). Так как

2/

для любой выпуклой вверх функции, Д(Ь) ^ 0 при Ь ^ 2Т; равенство Д(Ь) = 0 при некотором Ь > 2Т возможно только для упругого и вязкого элементов и модели Максвелла, для которых Д(Ь, Т) = 0 при Ь > 2Т. Предел Д(Ь) при Ь ^ те всегда равен нулю, так как Д(Ь) = Б(Ь — Т; Т) — Б(Ь; Т) и существуют пределы Б(Ь; Т) ^ V, Б(Ь — Т; Т) ^ V, где V = П(те). Так как

Д(Ь) = П(Ь) — 21Л(Ь — Т) + П(Ь — 2Т),

в случае n(t) ^ 0 будет A(t) ^ 0, т. е. A(t) возрастает на всем луче t ^ 2T (а |A(t| убывает).

Наибольшее значение A(t, T) достигается в точке t = T:

sup A(t) = A(T - 0) = П(Т) := Amax.

Наименьшее значение A(t) достигается в точке t = 2T:

inf A(t) = A(2T - 0) = П(2Т) - 2П(Т) := Amin.

Если n(t) = const, то Amin < 0 и |Amin| < Amax, так как

Л max -| Amini = П

в силу возрастания ФП (|Amin| = Amax только в случае П^) = const, т. е. для упругого элемента). С ростом T Amax возрастает, а Amin убывает:

A'min(T)=2(П(2Т) - П(Г)) < 0,

так как II(t) ^ 0.

Таким образом, доказана следующая лемма о свойствах функции (20).

Лемма 2. Пусть T > 0, <г > 0, а ФПП(^ дифференцируема, положительна, возрастает и строго выпукла вверх при t > 0. Тогда функция влияния (20) обладает следующими свойствами:

1) при t е (0; T) A(t) = П(^ — положительная возрастающая выпуклая вверх функция;

2) при t е (T; 2T) A(t) — убывающая функция и A(2T - 0) < 0;

3) на луче t ^ 2T A(t) выражается формулой (16), A(t) < 0 и предел A(+to) = 0;

4) ряд Y1 ь=2 A(kT) сходится, а его сумма £ равна vT + П(0) - П(Т); интеграл

I = A(t)dt

J2T

сходится, и для него справедлива оценка

[vT + П(0) - П(Т)]Т < I < [vT + П(Т) - П(2Т)]T;

5) у A(t) имеются ровно три точки разрыва на вещественной оси: t = 0, t = T и t = 2T, скачки в них равны П(0), -2П(0), П(0), а скачки производной— Й(0), -2П(0) и Й(0); у нерегулярных моделей это точки устранимого разрыва и A(t) непрерывна на всей оси;

6) наибольшее и наименьшее значения A(t) достигаются в точках t = T - 0 и t = 2T - 0:

Amax = A(T - 0) = П(Г) > 0,

Amin = A(2T - 0) = ^2T) - 2П(T) < -П(0) < 0, | Amini < Amax;

Amax и |Amin| возрастают с ростом T, а при T ^ 0 Amax ^ П(0), Amin ^ -П(0);

7) если П(Ь) ^ 0, то Д(Ь) выпукла вниз на интервале (Т; 2Т) и возрастает на всем луче Ь ^ 2Т.

Замечание 3. Пункт 4 вытекает из леммы 1 по интегральному признаку сходимости рядов.

Для моделей с ФП вида (4) (Максвелла, Фойгта, Кельвина и РеМ-4) Д(Ь) = —7е-Л + 27е-Л(4-Т) — 7в-Л(4-2Т) = —7(еЛТ — 1)2е-Л при Ь ^ 2Т, Д(Ь) не зависит от а, в и возрастает (П(Ь) = 7Л3е-Л4 > 0), £1 = е-лт £2, £2 = —7 (еЛТ — 1)/(еЛт + 1).

В самом деле,

£2 = —7(еЛТ — 1)2 £ е-2кЛТ =

= —7(еЛТ — 1)2е-2ЛТ (1 — е-2ЛТ )-1 = — 7 (еЛТ — 1)(еЛТ + 1)-1.

Для модели со степенной ФП П = сЬи, и € (0; 1), также П(Ь) > 0, и Д(Ь, Т) имеет вид (17).

На рис. 2 приведены графики функции (20) с Т = 5 для шести моделей разных классов:

1) модели со степенной ФП П = сЬи, с = 1, и = 0.5 (черная кривая 1) и с с = 0.5, и = 0.1 (штриховая кривая 1');

2) РеМ-4 с ФП (4) при Л = 0.1 (время ретардации т = 1/Л = 10), а = 0.1, в = 1.5, 7 = 1 (голубая кривая);

3) модели Фойгта с П = 7 — 7е-Л при Л = 0.1, 7 = 1 (красная кривая 3);

4) модели Максвелла (синяя) с П = аЬ + в, а = 0.1, в = 0.5 (время релаксации т = в/а = 5);

5) упругого элемента с в = 0.5 (зеленые штриховые «ступеньки» 5);

6) вязкого элемента с а = 0.1, в = 0 (синий «зуб»).

Нумерация кривых на рисунке совпадает с порядком перечисления моделей.

Упругий элемент, модель Максвелла и РеМ-4 регулярны (П(0) = 0), и поэтому их Д(Ь) имеют ненулевые скачки в точках Ь = 0; Т; 2Т; у остальных моделей П(0) = 0, и поэтому скачков нет, и Д(Ь) непрерывны при всех Ь ^ 0. Отклик (20) для РеМ-4 (и РеМ-3) на луче Ь > 2Т совпадает с откликом (20) модели Фойгта с теми же значениями Л и 7 (так как их ФП отличаются лишь линейным слагаемым):

Д(Ь) = —7 (еЛТ — 1)2е-Л4.

Для близкого к нулю показателя и реакция степенной модели (см. штриховую кривую 1' для и = 0.1) близка к реакции упругого элемента с П = с (Д(0) = +те; свойство Д(Т + 0) = —те создает резкий провал, визуально не отличимый от скачка вниз, а свойства Д(2Т + 0) = +те и Д(Ь) ~ 0 при Ь > 2Т хорошо приближают нулевой отклик упругого элемента). Для близких к единице и реакция степенной модели близка к реакции вязкого элемента (нет провала при Ь = Т и Л(Ь) ~ 0 на луче Ь > 2Т). Штрихпунктиром показан

Рис. 2.

для сравнения график функции £(£; Т) из (6) (т. е. кривая обратной ползучести с а = 1) для степенной модели: при Ь ^ те она стремится к 7г>Т = 0, но гораздо медленнее, чем Д(Ь), связанная с ней формулой (20).

7. Деформация в случае симметричного циклического растяжения-сжатия. Пусть периодическая программа нагружения состоит из одинаковых прямоугольных полуциклов растяжения и сжатия длительности Т и амплитуды <7:

те те

<

г=1 к=0

где Ьг = ¿Т. В силу линейности и инвариантности оператора (1) относительно сдвигов отклик

те

е

г=1

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

тете

е

¿=0 .7=0

Ряд (22) сходится при любом Ь (равномерно на любом отрезке), ибо в точке Ь отличны от нуля лишь конечное число членов ряда. На ¿-том периоде Ь € ((2г - 2)Т, 2гТ), ¿ = 1, 2,..., отклик (22) представляется конечной суммой с ] ^ ¿ - 1:

г-1

е(Ь) = 7 £ Д(Ь - 2Т?) при Ь € ((2¿ - 2)Т, 2¿T). ¿=0

При Ь € (^ - 2)Т, (2¿ - 1)Т), ¿ = 1, 2,..., КП (22) монотонно возрастает и выпукла вверх, так как в представлении (22) все ненулевые слагаемые Д(Ь - 2Т?'), ] ^ ¿ - 1, возрастают и выпуклы вверх на интервалах

(2iT-2T, 2iT-T) (они совпадают с (2jT, 2jT+T) или лежат в (2jT+2T, те)). Наибольшее значение деформации на интервале (Í2¿—2; Í2i-i) (и в i-том цикле):

г—i

ёг = sup e(t) = e(Í2¿—i - 0) = а ^ Д((2г - 2j - 1)T - 0) =

(Í2i-2;Í2i-l) j=0

г—2

= оД(Т - 0) + а ^ Д((2г - 2j - 1)T), j=0

т. е.

г

0г = аП

m=2

г—1

Si—1 = £ Д((2к + 1)T). k=i

Здесь SjLi, i > 1, — частичная сумма ряда из леммы 1. Так как Д(^ ^ 0 при t ^ 2T (см. лемму 2), все Д((2к + 1)T) ^ 0, Sii_i ^ 0 и поэтому последовательность {ёг} убывает (ёг+ - ёг = стД(^ + 1)T) ^ 0) и ёг ^ ёi = аП(Т).

Докажем,, что ёг > <тП(0) ^ 0, а значит, последовательность {ёг} ограничена снизу и поэтому при i ^ те имеет предел ё = аМ (T), причем ё ^ оП(0) ^ 0.

С учетом (18) величину ёг в формуле (23) можно выразить через сумму

Si2—i = Ek—=i Д(2кТ):

ёг/а = П(Т) + Si—i = П(Т) + [П(0) - П(Т) + n((2i - 1)T) - n((2i - 2)T) - S2—i], т. е.

ёг = оП

Так как n((2i - 1)T) - n((2i - 2)T) > 0 (в силу возрастания ФП) и -S2_i ^ 0, величина ёг > аП(0). При v = 0 можно уточнить оценку по теореме Лагран-жа:

n((2i - 1)T) - n((2i - 2)T) > 0 = П(£)Т,

где £ G ((2i - 2)T; (2i - 1)T). Поскольку II(t) не возрастает, справедливы оценки

П(£) ^ IÍ((2i - 1)T) ^ П(+те) = v, n((2i - 1)T) - n((2i - 2)T) ^ vT

и поэтому ёг/а ^ П(0) +vT (равенство только в случае IT(t) = const, т. е. только для упругого и вязкого элементов и модели Максвелла). Таким образом, доказаны двусторонние оценки для последовательности ёг и ее предела ё :

аП

Отметим, что для всех регулярных моделей

7П(0) + 7vT > 0,

а для нерегулярных с v = 0

Ш(0) + OvT = 0.

Из (23) и (24) следует, что предел последовательности {е} при i ^ те выражается через пределы Xi и X2 частичных сумм S1 и S2 рядов из леммы 1 формулами

е = <гП

При t g (t2i-1; t2i), i = 1, 2,..., КП (22) убывает и выпукла вниз (нестрого), так как в сумме (22)) все ненулевые слагаемые A(t — 2Tj), j ^ i — 1, убывают и выпуклы вниз на интервалах ((2i — 1)T, 2iT) (они совпадают с интервалом (2jT + T; 2jT + 2T), на котором функция A(t — 2Tj) убывает и выпукла вниз). Только для упругого элемента e(t) = const = —o. Минимальное значение e(t) на (t2i-1; t2i) (и в i-том цикле) вычисляется по формуле (22):

i— 1

£i = inf e(t)= e(t2i — 0) = o V A(2(i — j )T — 0),

(t2i-i;t2i) j_Q

или

ei = 7A

A(2kT ) - "" C*2

fc=2

Так как A(t) ^ 0 при t ^ 2T, то

£ < стА(2Т - 0) = ст[П(2Т) - 2П(Т)] < —стП(0) < 0, г = 1, 2,...

(лишь для вязкого элемента £ = 0), | ^ <гП(0) и последовательность {е,} убывает. Только для упругого и вязкого элементов и модели Максвелла, для которых П = а£ + в и А(£, Т) = 0 при £ > 2Т, последовательность {е,} постоянна: £ = — стП(0) = —ств.

Чтобы доказать ограниченность последовательности {е,} снизу и найти оценку сверху для преобразуем (27), выразив через 51 по формуле (18):

£ = —стП

Так как —стб;1 ^ 0 и П((2г + 1)Т) — П(2гТ) > 0, то £ ^ — стП(Т). При V = 0 можно уточнить эту оценку, как и при доказательстве (25): из

П((2г + 1)Т) — П(2гТ) ^ vT

следует, что

£ ^ — стП(Т) + отТ (равенство лишь для упругого и вязкого элементов и модели Максвелла).

Таким образом, доказаны двусторонние оценки для е, и предела е:

— аП(Т) + ОиТ < е, < аП(2Т) — 2аП(Т) < —аП(0) < 0

и — аП

В силу (23) и (28) имеем

0.5(а + е,) = 0.5а [—Д((2г + 1)Т) + П((2г + 1)Т) — П(2гТ)] =

= 0.5а[П

Из убывания последовательности {е,} и ее ограниченности снизу следует, что при г ^ те она имеет предел е = ат(Т), а из (27) и (28) следует, что он выражается формулами

е = а £2 — аП

Отметим, что для сходимости последовательностей {е,} и {е,} в случае симметричного цикла (21) не требуется, чтобы ФП была ограниченной, как было в случае импульсных нагружений (10) (см. п. 4). В силу формул (26), (31) и леммы 2

е + е = аП(Т) + + а£2 — аП(0) = ^Т

и последовательность — = + е, сходится при г ^ те к числу ОиТ ^ 0.

Последовательность средних е, := 0.5(е, + е,) неотрицательна и убывает (так как {е} и {е,} убывают). Вместе с ограниченностью (сходимостью) {е,} и {е,} это означает, что при симметричном циклическом нагружении (21) рэт-четинг не происходит, деформация и петля гистерезиса стабилизируются, т. е. линейное ОС (1) верно воспроизводит эффект, наблюдаемый в испытаниях большинства реономных материалов при симметричных нагружениях.

У моделей с V = 0

а = — е = аП(0) — а£2 = аП(Т) + а £1

и с ростом числа циклов происходит симметризация деформации (22) относительно оси е = 0:

е, ^ 0, а — ^ 0.

Из (25) и (29) следует, в частности, что при Т ^ 0 е ^ аП(0), е ^ —аП(0) и а, ^ аП(0), е, ^ —аП(0) для любого фиксированного г, т. е. при стремлении частоты циклов к бесконечности, размах деформации уменьшается и в пределе стремится к мгновенно-упругой реакции (что и наблюдается в испытаниях реономных материалов до момента, пока не становится заметным влияние тепловыделения на механические свойства).

Подытожим обнаруженные свойства кривой ползучести (22) для произвольной периодической программы нагружения (21) с симметричным циклом.

Теорема 2. Пусть Т > 0, а > 0, а функция ползучести П(Ь) дифференцируема, положительна, возрастает и выпукла вверх при Ь > 0. Тогда:

1) на интервалах £ € ((2г — 2)Т, (2г — 1)Т), г — 1,..., циклическая КП (22) возрастает и выпукла вверх, а на интервалах ((2г — 1)Т, 2гТ) е(£) убывает (а если П(£) ^ 0, то и выпукла вниз);

2) в точках £ = кТ КП (22) нерегулярных моделей имеют разрывы первого рода со скачками 2стП(0)(—1)к и скачками скорости деформации £(£), равными 2стП(0)(—1)к; КП регулярных моделей непрерывны и в точках £ = кТ, а формула для скачков £(£) сохраняется;

3) максимальное и минимальное значения и деформации е(£) в г-том цикле выражаются формулами (23), (24) и (27), (28); для и справедливы двусторонние оценки (25) и (29) (не зависящие от г) и

+ £ > (гоТ ^ 0; о > |£к | при любых г, к ^ 1;

4) последовательности {£,} и {£,} убывают, ограничены и сходятся, для их пределов £ и £ справедливы формулы (26), (31) и £ + £ — отТ, а также двусторонние оценки (25) и (29);

5) КП (22) ограничена на луче £ ^ 0 и для нее справедливы оценки

£ — ст — стП(0) ^ £(£) ^ £1 — оП(Т) при £ > 0, £ ^ £(£) ^ £+1 — стП(Т) + при £ ^ 2гТ;

6) последовательность £ :— 0.5(£ + £,) (см. (30)) убывает, а при г — те £ — 0.5отТ;

7) если ФП обладает свойством V — 0, то £ — —£ , и с ростом числа циклов происходит симметризация деформации (22) относительно оси £ — 0 : £ — 0 и — || —^ 0;

8) при Т — 0 для любого г — стП(0), £ — — <тП(0) и £ — стП(0), £ — —оП(0).

Замечание. Функции £,(Т) и (Т) не обязаны быть монотонными. Достаточные условия их монотонности:

1) если функция Ш(£) выпукла вверх (т. е. 2П(£) + Ш(£) ^ 0, это выполняется, в частности, для степенных ФП), то £,(Т) и £,(Т) при симметричном нагружении убывают, а размах при импульсном нагружении (15), наоборот, возрастает (ибо производная по Т каждого слагаемого А(2^Т) сумм (23), (27) и (15) отрицательна);

2) если Ш(£) выпукла вниз (т. е. 2П(£) + Ш(£) ^ 0, убывает), то £,(Т) и о(Т) возрастают.

У моделей Кельвина и Фойгта Ш(£) — 7Л£в-Л: сначала выпукла вверх, а при больших £ — вниз. Поэтому при малых и больших значениях ^Т £,(Т) и £,(Т) могут вести себя по разному.

Из теоремы 2 следует, что ОС (1) с произвольной ФП моделирует отсутствие рэтчетинга при симметричных циклических нагружениях. Это и наблюдается в испытаниях структурно стабильных материалов, одинаково сопротивляющихся растяжению и сжатию [3, 7,9,10,13,17,32].

На рис. 3 приведены КП (22) пяти моделей (тех же, что и на рис. 2) для симметричного циклического нагружения (1) с Т — 5, ст — 1:

1) КП модели со степенной ФП П — с£и, с — 1, и — 0.5 (черная КП);

2) КП РеМ-4 с П — а£+в—7в-Л4 с Л — 0.1 (время ретардации т — 1/Л — 10), а — 0.1, в — 1.5, 7 — 1 (голубая КП);

3) КП модели Фойгта с П = 7 — 7в Л4, Л = 0.1, 7 = 1 (красная);

4) КП модели Максвелла (синяя) с П = а* + в, а = 0.1, в = 0.5 (время релаксации т = в/а = 5, Д(*) = 0 при * ^ 2Т, £1 = £2 = 0, = П(Т) = = в + аТ = е, е, = —ап(0) = — ав = е);

5) КП вязкого элемента с П = а*, а = 0.1 (синяя ломаная: она не имеет разрывов и периодична, так как Д(*) = 0 при * ^ 2Т, е(*) ^ 0 на луче * ^ 0, е, = аП(Т) = <гаТ, е, = 0).

КП модели Максвелла и РеМ-4 имеют в точках = ¿Т разрывы первого рода со скачками ±аП(0). У модели Фойгта и степенной модели П(0) =0 и поэтому их КП непрерывны при всех * > 0. У этих двух моделей V = 0 и поэтому происходит симметризация КП (22) относительно оси е = 0. У РеМ-4, РеМ-2 и вязкого элемента (СиМ-1) V = а > 0 и поэтому 0.5(е, + е,) ^ 0.5ааТ > 0 и е = |е| + а аТ > |е| (штриховые прямые на рис. 3 — е = 0.5а аТ = 0.25 и е = а аТ = 0.5). Две штрихпунктирные красные кривые — КП модели Фойгта при постоянном напряжении а = ±1.

Рис. 3.

8. О кривой циклической ползучести и ее отклонении от кривой ползучести при среднем напряжении. Произвольная программа двухступенчатого циклического нагружения (19) (в моменты времени = ¿Т происходит периодический перескок с напряжения Ст1 на Ст2 и обратно) разлагается в сумму постоянного напряжения ат = 0.5(а1 + а2) и симметричного ступенчатого нагружения (21) с тем же периодом 2Т, амплитудой а = 0.5а, а = Ст1 — а2, и теми же скачками напряжения в точках = кТ, что и у (19). Отклик на нее представим в виде суммы обычной КП для среднего напряжения ат и отклика г(*) на симметричное циклическое нагружение:

е

п

г(*)=0.5а ^ Д(* — 2Т(г — 1),Т), * е [0;2Тп), п = 1, 2,....

г=1

Свойства г(*) в случае а = 0.5а > 0 (когда а1 > а2) установлены в теореме 2. Если же Ст1 < Ст2, то а < 0, знаки всех производных КП (22) меняются на противоположные, интервалы возрастания и убывания, выпуклости вверх и

вниз, точки максимума и минимума меняются местами (но сохраняется их чередование). Поэтому при любых 0 1, 02 (вне зависимости от знаков 0 = 01 — 02 и и1 — 0m) и T > 0 для КП (32) справедлива следующая теорема.

Теорема 3. Пусть функция ползучести n(t) дифференцируема, положительна, возрастает и выпукла вверх при t > 0. Тогда циклическая кривая ползучести (32) и ее отклонение r(t) = e(t) — 0mn(t) от КП для среднего напряжения 0m = 0.5(и1 + и2) обладают следующими свойствами:

1) функция r(t) монотонна на интервалах (kT, kT + T), интервалы возрастания и убывания чередуются; ее значения f := r(t2j-1 — 0) и fj := r(t2j — 0) в точках экстремума t2j-1 = (2i — 1)T и t2j = 2iT (максимальные отклонения в i-том цикле) выражаются через n(t) формулами (23), (24) и (27), (28) с 0 = 0.50 = 0.5(и1 — и2) и заменой fj, f на а и r^

2) rr < 0, последовательности fj, r и r := 0.5(fj + rj монотонны, ограничены и сходятся, для их пределов f, r и 0.5(f + r) справедливы формулы (26), (31) и f + r = 0vT, где v = П(те);

3) rjfj > 0, rjrj < 0, последовательности |fj| и |rj| убывают, а |rj| возрастает и

4) отклонение r(t) ограничено на луче t ^ 0 :

—П(Т) + vT — E1 < r(t)/f < П(Т), где S1 — сумма ряда из леммы 1, S1 ^ 0;

последовательность норм pj = sup{|r(t)| | t G (2iT; те)} убывает и Pi ^ |f|, где

f = f-n(T) + 0X1 = f П(0) — f-^2 + fvT;

5) полусумма экстремальных значений КП (32) выражается формулой

f = 0.50m[(1 + A)n(2iT) + (1 — A)n(2iT — T)],

где A = <f/0m = (01 — и2)/(и1 + и2) — «коэффициент амплитуд»; последовательности отклонений

= f — 0mn(2iT — T) и 5+ = 0mn(2iT) — e

убывают по модулю и сходятся к 0.5ит(1 + A)vT и 0.5ит(1 — A)vT при i ^ те; отношение

= (1 — A)/(A + 1)= 02/01

не зависит от i, 0m и n(t);

6) в случае A G (—1;1) (т. е. в случае 02/01 > 0) для ej при любом i справедливы оценки

n(2iT — T) < n(2iT — T) + 0.5(1 + A)vT < £j/0m,

ej/0m < n

7) у моделей с V = 0 при г ^ те происходит асимптотическая симметризация отклонения циклической КП (32) от КП при среднем напряжении : г = —Г, |г,| — |г,| ^ 0, г, ^ 0 и е, — атП(2гТ) ^ 0; если П(*) ограничена, то е, ^ атП(те) и происходит симметризация циклической КП (32) относительно прямой е = атП(те) (асимптоты КП е = атП(*));

8) в случае ат = 0 КП (32) всегда ограничена на луче * ^ 0, а в случае ат = 0 КП (32) ограничена тогда и только тогда, когда ограничена П(*).

Таким образом, ОС (1) с произвольной ФП П(*) моделирует отсутствие рэтчетинга (стабильность) при симметричных циклических нагружениях, а при наложении симметричного (ступенчатого) циклического возмущения на постоянную нагрузку не происходит ускорения ползучести по сравнению с ползучестью при среднем напряжении. Если же есть асимметрия в циклической компоненте нагрузки (например, длительности полуциклов растяжения-сжатия не совпадают или а_ = —а+), то при V = 0 происходит ускорение или замедление ползучести в зависимости от знака приращения пластической деформации (см. (9)) за цикл. В случае двухступенчатого цикла нагружения она выражается формулой р = v(а1T1 + а2Т2), а в случае произвольного количества ступеней нагружения (а&, ) в цикле — формулой р = .

Этот аналитический результат показывает, в частности, что наблюдаемый в испытаниях различных материалов «эффект интенсификации ползучести по отношению к среднему напряжению» [3,16,19,26,32,42] может быть вызван не только физическими причинами, но и неидеальностью задания (регистрации) циклической компоненты напряжения: отклонения от симметричности цикла нагружения (Т = Т2 или а_ = — а+) могут дать заметный (накапливающийся) вклад в наблюдаемое смещение циклической КП.

Поведение КП (32) в случае произвольного циклического нагружения (19) и ее форма на участках разгрузки (т. е участках с меньшим модулем напряжения, где слагаемые в (32) имеют разные знаки) могут существенно меняться в зависимости от знака и величины отношения

А = а/ат =

амплитуды а к среднему ат (или от коэффициента асимметрии цикла а = = а2/а1 = (1 — А)/(А + 1)) и от свойств ФП. В частности, модуль е(*) не обязан убывать на участках разгрузки и может возрастать при а2/а1 > а+ или иметь точку минимума внутри интервалов (кТ, (к + 1)Т) (а не на концах), если а2/а1 е (а_; а+), где 0 ^ а_ ^ а+ ^ 1, а_,а+ зависят от ФП и Т (см. [36] и рис. 4). Значения КП (32) в концах полуциклов е, := е(*2,_1 — 0) и е, := е(*2, — 0) в общем случае уже не обязаны быть наибольшим и наименьшим значениями в г-том цикле. Детальный анализ общих свойств КП (32), описание возможных случаев и коррекция некоторых распространенных представлений о циклических КП — тема будущих работ автора. Здесь отметим только, что средняя деформация в цикле е, := 0.5(е, + е,) выражается формулой

е, = 0.5ат[(1 + А)П(2гТ) + (1 — А)П(2гТ — Т)],

а в случае A е (—1;1) (т. е. когда a = o2/o1 > 0) для нее справедлива универсальная двусторонняя оценка

Отклонения величин е от значений КП при среднем напряжении выражаются формулами

5+ := omn(2iT) - е = 0.5om(1 - A)S(2iT), 5- := е - omn(2iT - T) = 0.5om(1 + A)S(2iT),

где S(2iT) = П(2гТ) - n(2iT - T) > 0-значения функции (6). Последовательности |5+| и |5-| убывают, так как S(t) убывает (см. п. 2). Так как S(те) = vT, пределы последовательностей 5+ и 5- равны 0.5om(1 - A)vT и 0.5om(1 + A)vT и поэтому можно уточнить оценку (34) до (33)). Лишь у моделей с v = 0 5+ ^ 0 и 5- ^ 0 для любых om и A. С ростом A модуль |5+| убывает, |5-| возрастает при всех i ^ 1, а для импульсного нагружения (19) (с 02 = 0 или oi =0) A = 1 и 5+ = 0 или A = -1 и 5- =0. Примечательно, что отношение отклонений

5+/5- = (1 - A)/(A + 1) = a

зависит лишь от A (от a = o2/oi), но не зависит от i, om и функции ползучести. Это свойство циклических КП (32), порождаемых линейным ОС (1) при нагружениях вида (19), можно использовать при анализе данных испытаний как индикатор границы линейности поведения материала и (не)применимости линейного ОС (1).

На рис. 4 приведены КП (32) модели со степенной ФП П = ctu, c = 0.1, u = 0.5, для циклических нагружений (19) с разными коэффициентами амплитуд A = o/om = {1; 1/2; 1/3; 1/4}, т. е. с 02/01 = {0; 0.(3); 0.5; 0.6} (черные КП 1-4). Они показывают, что (в зависимости от величины 02/01) циклические КП на участках разгрузки могут не только убывать, но и возрастать

Рис. 4.

или иметь точку минимума внутри интервалов (кТ, (к + 1)Т) (см. [36]), и что убывание восстанавливается на последующих циклах. Для удобства совмещения КП на одном рисунке выбрана постоянная амплитуда — 1 и разные средние напряжения стт — {1; 2; 3; 4}. Обычные КП для стт показаны штрих-пунктирными голубыми линиями. Красные кривые 5 и 6 — КП (32) модели Фойгта (П — 7 —7в-Л4, Л — 0.1, 7 — 0.25) для стт — 1 и стт — 3 (ст2/ст1 — 0; 0.5).

Качественно схожими свойствами обладают экспериментальные КП различных реономных материалов при ступенчатых нагружениях [1-12, 24, 25, 27-29]. Общие свойства КП, доказанные в теоремах 1-3, позволяют приглядеться к экспериментальным кривым пристальнее и попытаться разглядеть дополнительные детали, на которых акцентируют внимание теоремы (например, с целью адекватного выбора типа и характеристик функции ползучести или поиска признаков, свидетельствующих об отклонении от линейности в поведении материала, или с целью поверки надежности регистрации и обработки данных измерений).

9. Заключение. В статье продолжен качественный анализ линейного определяющего соотношения (ОС) вязкоупругости (1) с произвольной (возрастающей выпуклой вверх) функцией ползучести. Аналитически изучены общие свойства и качественные особенности кривых ползучести (32), порождаемых ОС (1) при циклических двухступенчатых нагружениях (19) с произвольным коэффициентом асимметрии цикла. Для импульсных, симметричных и произвольных циклических ступенчатых нагружений получены общие формулы и точные двусторонние оценки для максимальных и минимальных значений деформации в каждом цикле и их пределов, для отклонения кривой циклической ползучести от обычной кривой ползучести при среднем напряжении цикла, для скорости накопления пластической (необратимой) деформации и рэтчетинга. Исследованы интервалы монотонности и выпуклости кривых циклической ползучести, поведение последовательностей максимальных и минимальных значений деформации и их полусуммы (условия их ограниченности, монотонности, сходимости), возможности моделирования циклического упрочнения и разупрочнения, а также зависимость всех обнаруженных свойств от характеристик функции ползучести и параметров цикла нагружения. Основные свойства кривых циклической ползучести ОС (1) собраны в теоремах 1-3 (пп. 4, 7, 8 статьи). В частности, выявлены ключевые роли выпуклости вверх функции ползучести и величины предела ее производной V — П(те): доказано, что равенство V нулю является критерием полного восстановления после снятия (произвольной) нагрузки и затухания памяти ОС (1), критерием отсутствия накопления пластической деформации при циклической ползучести и критерием асимптотической симметризации кривых циклической ползучести относительно КП при среднем напряжении.

Сравнение обнаруженных свойств теоретических кривых ползучести (32) с типичными свойствами экспериментальных кривых ползучести вязкоупру-гопластичных материалов при ступенчатых нагружениях позволило исследовать возможности ОС (1) по описанию различных эффектов при циклических ступенчатых нагружениях (в частности, выявлен ряд эффектов, которые линейное ОС (1) не способно описывать), сферы влияния качественных характеристик функции ползучести и способы идентификации и настройки ОС. В частности, доказано, что ОС (1) с произвольной ФП моделирует

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

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

Конкурирующие интересы. У меня нет конкурирующих интересов. Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.

Финансирование. Работа выполнена при поддержке РФФИ (грант № 17-08-01146_а).

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

1. Качанов Л. М. Теория ползучести. М.: Физматгиз, 1960. 456 с.

2. Наместников B. C., Хвостунков А. А. Ползучесть дюралюмина при постоянных и переменных нагрузках // ПМТФ, 1960. Т. 1, №4. С. 90-95.

3. Kennedy A. J. Processes of creep and fatigue in metals / Wiley series on the science and technology of materials. vol. 19. New York: Wiley, 1963. 480 pp.

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

5. Odqvist F. K. G. Mathematical Theory of Creep and Creep Rupture. Oxford: Clarendon ress, 1966. 170 pp.

6. Самарин Ю. П., Сорокин О. В. О ползучести поливинилхлоридного пластиката при переменных нагрузках // Докл. АН СССР, 1970. Т. 195, №2. С. 333-336.

7. Бугаков И. И. Ползучесть полимерных материалов. М.: Наука, 1973. 287 с.

8. Findley W. N., Lai J. S., Onaran K. Creep and Relaxation of Nonlinear Viscoelastic aterials. Amsterdam: North Holland, 1976. 368 pp.

9. Стрижало В. А. Циклическая прочность и ползучесть металлов при малоцикловом нагружении в условиях низких и высоких температур. Киев: Наукова думка, 1978. 238 с.

10. Kujawski D., Kallianpur V., Krempl E. An experimental study of uniaxial creep, cyclic creep and relaxation of aisi type 304 stainless steel at room temperature // J. Mech. Phys. Solids, 1980. vol.28, no. 2. pp. 129-148. doi: 10.1016/0022-5096(80)90018-6.

11. Шестериков С. А., Локощенко А. М. Ползучесть и длительная прочность металлов/ Итоги науки и техники. Сер. Механ. деформ. тверд. тела, Т. 13. М.: ВИНИТИ, 1980. С. 3-104.

12. Малинин Н. Н. Расчеты на ползучесть элементов машиностроительных конструкций. М.: Машиностроение, 1981. 221 с.

13. Москвитин В. В. Циклическое нагружение элементов конструкций. М.: Наука, 1981. 344 с.

14. Локощенко A. M., Наместникова И. В., Шестериков С. А. Описание длительной прочности при ступенчатом изменении напряжения// Проблемы прочности, 1981. №10. С. 47-51.

15. Cho U. W., Findley W. N. Creep and Plastic Strains of 304 Stainless Steel at 593°C Under Step Stress Changes, Considering Aging // J. Appl. Mech., 1982. vol. 49, no. 2. pp. 297-304. doi:10.1115/1.3162084.

16. Голуб В. П. Циклическая ползучесть жаропрочных никелевых сплавов. Киев: Наумова думка, 1983. 224 с.

17. Гохфельд Д. А., Садаков О. С. Пластичность и ползучесть элементов конструкций при повторных нагружениях. М.: Машиностроение, 1984. 256 с.

18. Малинин H Н. Ползучесть в обработке металлов давлением. М.: Машиностроение, 1986. 221 с.

19. Голуб В. П. Исследования в области циклической ползучести материалов (обзор) // Прикладная механика, 1987. Т. 23, №12. С. 3-19.

20. Tschoegl N. W. The Phenomenological Theory of Linear Viscoelastic Behavior. Berlin: Springer, 1989. 769 pp.

21. Никитенко А. Ф. Ползучесть и длительная прочность металлических материалов. Новосибирск: НГАСУ, 1997. 278 с.

22. Локощенко А. М. Ползучесть и длительная прочность металлов в агрессивных средах. М.: МГУ, 2000. 179 с.

23. Радченко В. П., Саушкин М. Н. Ползучесть и релаксация остаточных напряжений в упрочнённых конструкциях. М.: Машиностроение-1, 2005. 226 с.

24. Hamouda B. H., Laiarinandrasana L., Piques R. Viscoplastic behavior of a medium density polyethylene (MDPE): constitutive equations based on double nonlinear deformation model// Int. J. Plasticity, 2007. vol. 23, no. 8. pp. 1307-1327. doi: 10.1016/j.ijplas.2006. 11.007.

25. Betten J. Creep Mechanics. Berlin, Heidelberg: Springer-Verlag, 2008. 267 pp.

26. Радченко В. П., Кичаев П. Е. Энергетическая концепция ползучести и виброползучести металлов. Самара: Самар. гос. тех. ун-т, 2011. 157 с.

27. Darabi M. K, Al-Rub R. K. A., Masad E. A., Huang C.-W., Little D. N. A modified viscoplastic model to predict the permanent deformation of asphaltic materials under cyclic-compression loading at high temperatures // Int. J. Plasticity, 2012. vol. 35. pp. 100-134. doi:10.1016/j.ijplas.2012.03.001.

28. Bergstrom J. S. Mechanics ofSolid Polymers. Theory and Computational Modeling. Elsevier: William Andrew, 2015. 520 pp.

29. Локощенко А. М. Ползучесть и длительная прочность металлов. М.: Физматлит, 2016. 504 с.

30. Локощенко А. М., Фомин Л. В. Длительное разрушение пластин при переменных изгибающих моментах в присутствии агрессивной среды // Прикл. матем. и механ., 2016. Т. 80, №2. С. 276-284.

31. Dandrea J., Lakes R. S. Creep and creep recovery of cast aluminum alloys // Mech. Time-Depend. Mater., 2009. vol.13, no. 4. pp. 303-315. doi: 10.1007/s11043-009-9089-6.

32. Taleb L., Cailletaud G. Cyclic accumulation of the inelastic strain in the 304L SS under stress control at room temperature: Ratcheting or creep? // Int. J. Plasticity, 2011. vol. 27, no. 12. pp. 1936-1958. doi: 10.1016/j.ijplas.2011.02.001.

33. Khan F., Yeakle C. Experimental investigation and modeling of non-monotonic creep behavior in polymers// Int. J. Plasticity, 2011. vol.27, no. 4. pp. 512-521. doi: 10.1016/j. ijplas.2010.06.007.

34. Drozdov A. D. Time-dependent response of polypropylene after strain reversal // Int. J. of Solids and Structures, 2010. vol.47, no. 24. pp. 3221-3233. doi: 10.1016/j.ijsolstr.2010. 08.001.

35. Drozdov A. D., Dusunceli N. Unusual mechanical response of carbon black-filled thermoplastic elastomers // Mechanics of Materials, 2014. vol. 69, no. 1. pp. 116-131. doi: 10.1016/ j.mechmat.2013.09.019.

36. Хохлов А. В. Свойства семейств кривых ползучести при ступенчатом нагружении линейного определяющего соотношения вязкоупругости // Проблемы прочности и пластичности, 2G15. Т. 77, №4. С. 344-359.

37. Хохлов А. В. Качественный анализ общих свойств теоретических кривых линейного определяющего соотношения вязкоупругости // Наука и образование, 2G16. №5. С. 187245, http://technomag.bmstu.ru/doc/840650.html.

38. Хохлов А. В. Нелинейная модель вязкоупругопластичности типа Mаксвелла: свойства кривых ползучести при ступенчатых нагружениях и условия накопления пластической деформации // Машиностроение и инженерное образование, 2G16. №3. С. 55-68.

39. Хохлов А. В. Асимптотическая коммутативность кривых ползучести при ступенчатом нагружении в линейной теории наследственности // Машиностроение и инженерное образование, 2G16. №1. С. 7G-82.

4G. Хохлов А. В. Кривые длительной прочности нелинейной модели вязкоупругопластич-ности типа Mаксвелла и правило суммирования поврежденности при ступенчатых на-гружениях// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2G16. Т. 2G, №3. С. 524-543. doi:10.14498/vsgtu1512.

41. Fatemi A., Yang L. Cumulative fatigue damage and life prediction theories: A survey of the state of the art for homogeneous materials // Int. J. Fatigue, 1QQ8. vol. 20, no. 1. pp. 9-34. doi:10.1016/S0142-1123(97)00081-9.

42. Голуб В. П. О некоторых эффектах ползучести при циклических нагружениях // Проблемы прочности, 1987. №5. С. 2G-24.

43. Радченко В. П., Кичаев Е. К., Симонов А. В. Энергетический вариант модели реологического деформирования и разрушения металлов при совместном действии статических и циклических нагрузок// ПМТФ, 2GGG. Т. 41, №3. С. 169-175.

44. Zheng X.-T., Xuan F.-Z., Zhao P. Ratcheting-creep interaction of advanced 9-12% chromium ferrite steel with anelastic effect // Int. J. Fatigue, 2011. vol. 33, no. 9. pp. 12861291. doi: 10.1016/j.ijfatigue.2011.04.009.

45. Баренблатт Г. И., Козырев Ю. И., Mалинин Н. Н., Павлов Д. Я., Шестериков С. А. О виброползучести полимерных материалов// ПМТФ, 1965. №5. С. 68-75.

46. Локощенко А. M., Шестериков С. А. О виброползучести // Инженерный журнал. Механика твердого тела, 1966. №3. С. 141-143.

47. Локощенко А. M. Виброползучесть металлов при одноосном и сложном напряженных состояниях// Изв. РАН. МТТ, 2G14. №4. С. 111-12G.

48. Хохлов А. В. Характерные особенности семейств кривых деформирования линейных моделей вязкоупругости// Проблемы прочности и пластичности, 2G15. Т. 77, №2. С. 139-154.

49. Хохлов А. В. Кривые длительной прочности, порождаемые линейной теорией вязко-упругости в сочетании с критериями разрушения, учитывающими историю деформирования// Труды МАИ, 2G16. №91. С. 1-32, http://www.mai.ru/science/trudy/ published.php?ID=75559.

5G. Хохлов А. В. Анализ общих свойств кривых ползучести при ступенчатом нагружении, порождаемых нелинейным соотношением Работнова для вязкоупругопластичных материалов// Вестник МГТУ им. Н. Э. Баумана. Сер. Естественные науки, 2G17. №3. С. 93-123. doi: 10.18698/1812-3368-2017-3-93-123.

51. Радченко В. П., Самарин Ю. П. Влияние ползучести на величину упругой деформации слоистого композита// Механика композитных материалов, 1983. Т. 19, №2. С. 231237.

52. Радченко В. П., Шапиевский Д. В. О дрейфе упругой деформации для нелинейно-упругих материалов вследствие ползучести // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2GG6. №43. С. 99-1G6. doi: 10.14498/vsgtu458.

53. Радченко В. П., Шапиевский Д. В. Mатематическая модель ползучести микронеоднородного нелинейно-упругого материала// ПМТФ, 2GG8. Т. 49, №3. С. 157-163.

54. Mелнис А. Э., Лайзан Я. Б. Нелинейная ползучесть компактной костной ткани человека при растяжении // Механика полимеров, 1978. Т. 14, №1. С. 97-1GG.

Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki

[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 2, pp. 326-361

d http://doi.org/10.14498/vsgtu1533

ISSN: 2310-7081

MSC: 74D05

Analysis of creep curves produced by the linear viscoelasticity theory under cyclic stepwise loadings

A. V. Khokhlov

Lomonosov Moscow State University,

Institute of Mechanics,

1, Michurinsky prospekt, Moscow, 119192, Russian Federation.

Abstract

Basic qualitative properties of the creep curves generated by the linear integral constitutive relation of viscoelasticity (with an arbitrary creep compliance) under cyclic piecewise-constant uni-axial loadings (with an arbitrary asymmetry stress ratio) are studied analytically. General formulas and a number of exact two-sided bounds are obtained for maximal, minimal and ratcheting strain values during each cycle, for their sequences limits, for the rate of plastic (non-recoverable) strain accumulation and for cyclic creep curve deviation from the creep curve at constant stress which is equal to the cycle mean stress. Their dependence on loading cycle parameters and creep compliance properties are analyzed. Monotonicity and convexity intervals of cyclic creep curves, sequences of maximal and minimal strain values and ratcheting strain sequence, their evolution with cycle number growth and conditions for their boundedness, monotonicity and convergence are examined. The linear viscoelasticity theory abilities for simulation of ratcheting, creep acceleration, cyclic hardening or softening and cyclic stability under symmetric cyclic loadings are considered.

The analysis carried out revealed the importance of convexity restriction imposed on a creep compliance and the governing role of its derivative limit value at infinity. It is proved that the limit value equality to zero is the criterion for non-accumulation of plastic strain, for memory fading and for asymptotic symmetrization of cyclic creep curve deviation from the creep curve at the mean stress.

The qualitative features of theoretic cyclic creep curves are compared to basic properties of typical test creep curves of viscoelastoplastic materials under cyclic multi-step uni-axial loadings in order to elucidate the linear theory applicability scope, to reveal its abilities to provide an adequate description of basic rheological phenomena related to cyclic creep and to develop techniques of identification and tuning of the linear constitutive relation. In particular, it is proved that the linear constitutive relation with an arbitrary (increasing convex-up) creep compliance function provides the absence of ratcheting and cyclic softening under symmetric cyclic multi-step loadings and the absence of creep acceleration whenever a symmetric cyclic loading is added to a constant load.

Article

3 © The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:

Khokhlov A. V. Analysis of creep curves produced by the linear viscoelasticity theory under cyclic stepwise loadings, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 2, pp. 326-361. doi: 10.14498/vsgtu1533 (In Russian). Author's Details:

Andrew V. Khokhlov http://orcid.org/0000-0002-9212-2579 Cand. Techn. Sci.; Senior Researcher; Lab. of Elasticity and Plasticity; e-mail: andrey-khokhlov@ya.ru

Analysis of creep curves produced by the linear viscoelasticity theory.

Keywords: linear viscoelasticity, cyclic creep, creep curves at piecewise-constant loading, asymmetry stress ratio, mean stress, creep acceleration, plastic strain, ratcheting, cyclic stability.

Received: 14th March, 2017 / Revised: 17th May, 2017 / Accepted: 12th June, 2017 / First online: 10th July, 2017

Competing interests. I have no competing interests.

Author's Responsibilities. I take full responsibility for submitting the final manuscript

in print. I approved the final version of the manuscript.

Funding. This work was supported by the Russian Foundation for Basic Research (project no. 17-08-01146_a).

References

1. Kachanov L. M. Creep Theory. Moscow, Fizmatgiz, 1960, 456 pp.

2. Namestnikov B. C., Khvostunkov A. A. Creep of duralumin at constant and variable loads, Prikl. Mekh. Tekh. Fiz, 1960, vol. 1, no. 4, pp. 90-95 (In Russian).

3. Kennedy A. J. Processes of creep and fatigue in metals, Wiley series on the science and technology of materials, vol. 19. New York, Wiley, 1963, 480 pp.

4. Rabotnov Yu. N. Polzuchest' elementov konstruktsii [Creep of Structure Elements]. Moscow, Nauka, 1966, 752 pp.

5. Odqvist F. K. G. Mathematical Theory of Creep and Creep Rupture. Oxford, Clarendon ress, 1966, 170 pp.

6. Samarin Yu. P., Sorokin O. V. On the creep of masticated polyvinyl chloride rubber under variable loads, Doklady Akad. Nauk SSSR, 1970, vol. 195, no. 2, pp. 333-336 (In Russian).

7. Bugakov I. I. Polzuchest' polimernykh materialov [Creep of polymer materials]. Moscow, Nauka, 1973, 287 pp.

8. Findley W. N., Lai J. S., Onaran K. Creep and Relaxation of Nonlinear Viscoelastic aterials. Amsterdam, North Holland, 1976, 368 pp.

9. Strizhalo V. A. Tsiklicheskaia prochnost' i polzuchest' metallov pri malotsiklovom na-gruzhenii v usloviiakh nizkikh i vysokikh temperatur [Cyclic Strength and Creep of Metals under Low-Cycle Loading at Low and High Temperatures]. Kiev, Naukova Dumka, 1978, 238 pp.

10. Kujawski D., Kallianpur V., Krempl E. An experimental study of uniaxial creep, cyclic creep and relaxation of aisi type 304 stainless steel at room temperature, J. Mech. Phys. Solids, 1980, vol.28, no. 2, pp. 129-148. doi: 10.1016/0022-5096(80)90018-6.

11. Shesterikov S. A., Lokoshchenko A. M. Creep and Long-Term Strength of Metals, Itogi nauki i tekhniki. Ser. Mekhan. deform. tverd. tela [Results of Science and Engineering, Ser. Mechanics of Deformable Solids], 13. Moscow, VINITI, 1980, pp. 3-104 (In Russian).

12. Malinin N. N. Raschety na polzuchest' elementov mashinostroitel'nykh konstruktsii [Creep Calculations on Parts of Machine Constructions]. Moscow, Mashinostroenie, 1981, 221 pp.

13. Moskvitin V. V. Tsiklicheskoe nagruzhenie elementov konstruktsii [Cyclic Loading of Structure Elements]. Moscow, Nauka, 1981, 344 pp.

14. Lokoshchenko A. M., Namestnikova I. V., Shesterikov S. A. Description of long-term strength with stepwise stress changes, Strength of Materials, 1981, vol. 13, no. 10, pp. 1240-1244. doi: 10.1007/BF00775885.

15. Cho U. W., Findley W. N. Creep and Plastic Strains of 304 Stainless Steel at 593°C Under Step Stress Changes, Considering Aging, J. Appl. Mech., 1982, vol.49, no. 2, pp. 297-304. doi:10.1115/1.3162084.

16. Golub V. P. Tsiklicheskaia polzuchest' zharoprochnykh nikelevykh splavov [Cyclic Creep of Creep-Resisting Nickel Alloys]. Kiev, Naumova dumka, 1983, 224 pp.

17. Gokhfel'd D. A., Sadakov O. S. Plastichnost' i polzuchest' elementov konstruktsii pri pov-tornykh nagruzheniiakh [Plasticity and creep in structural elements under repeated loading]. Moscow, Mashinostroenie, 1984, 256 pp.

18. Malinin H N. Polzuchest' v obrabotke metallov davleniem [Creep theories in metal forming]. Moscow, Mashinostroenie, 1986, 221 pp.

19. Golub V. P. Investigations into cyclic creep of materials (review), Soviet Applied Mechanics, 1987, vol.23, no. 12, pp. 1107-1121. doi: 10.1007/BF00884887.

20. Tschoegl N. W. The Phenomenological Theory of Linear Viscoelastic Behavior. Berlin, Springer, 1989, 769 pp.

21. Nikitenko A. F. Polzuchest' i dlitel'naia prochnost' metallicheskikh materialov [Creep and Long-Time Strength of Metallic Materials]. Novosibirsk, Novosibirsk State Architectural Univ., 1997, 278 pp.

22. Lokoshchenko A. M. Polzuchest' i dlitel'naia prochnost' metallov v agressivnykh sredakh [Creep and Long-Term Strength of Metals in Corrosive Media]. Moscow, Moscow State Univ., 2000, 179 pp.

23. Radchenko V. P., Saushkin M. N. Polzuchest' i relaksatsiia ostatochnykh napriazhenii v up-rochnennykh konstruktsiiakh [Creep and Relaxation of Residual Stresses in Hardened Structures]. Moscow, Mashinostroenie-1, 2005, 226 pp.

24. Hamouda B. H., Laiarinandrasana L., Piques R. Viscoplastic behavior of a medium density polyethylene (MDPE): constitutive equations based on double nonlinear deformation model, Int. J. Plasticity, 2007, vol.23, no. 8, pp. 1307-1327. doi: 10.1016/j.ijplas.2006.11.007.

25. Betten J. Creep Mechanics. Berlin, Heidelberg, Springer-Verlag, 2008, 267 pp.

26. Radchenko V. P., Kichaev P. E. Energeticheskaia kontseptsiia polzuchesti i vibropolzuchesti metallov [Energy Concept of Creep and Vibrocreep of Metals]. Samara, Samara State Techn. Univ., 2011, 157 pp.

27. Darabi M. K, Al-Rub R. K. A., Masad E. A., Huang C.-W., Little D. N. A modified viscoplastic model to predict the permanent deformation of asphaltic materials under cyclic-compression loading at high temperatures, Int. J. Plasticity, 2012, vol.35, pp. 100-134. doi:10.1016/j.ijplas.2012.03.001.

28. Bergstrom J. S. Mechanics of Solid Polymers. Theory and Computational Modeling. Elsevier, William Andrew, 2015, 520 pp.

29. Lokoshchenko A. M. Polzuchest' i dlitel'naia prochnost' metallov [Creep and long-lasting strength of metals]. Moscow, Fizmatlit, 2016, 504 pp.

30. Lokoshchenko A. M., Fomin L. V. Creep fracture of plates with variable bending moments in the presence of an aggressive medium, J. Appl. Math. Mech., 2016, vol. 80, no. 2, pp. 198204. doi: 10.1016/j.jappmathmech.2016.06.006.

31. Dandrea J., Lakes R. S. Creep and creep recovery of cast aluminum alloys, Mech. Time-Depend. Mater., 2009, vol.13, no. 4, pp. 303-315. doi: 10.1007/s11043-009-9089-6.

32. Taleb L., Cailletaud G. Cyclic accumulation of the inelastic strain in the 304L SS under stress control at room temperature: Ratcheting or creep?, Int. J. Plasticity, 2011, vol. 27, no. 12, pp. 1936-1958. doi: 10.1016/j.ijplas.2011.02.001.

33. Khan F., Yeakle C. Experimental investigation and modeling of non-monotonic creep behavior in polymers, Int. J. Plasticity, 2011, vol.27, no. 4, pp. 512-521. doi: 10.1016/j. ijplas.2010.06.007.

34. Drozdov A. D. Time-dependent response of polypropylene after strain reversal, Int. J. of Solids and Structures, 2010, vol.47, no. 24, pp. 3221-3233. doi: 10.1016/j.ijsolstr.2010. 08.001.

35. Drozdov A. D., Dusunceli N. Unusual mechanical response of carbon black-filled thermoplastic elastomers, Mechanics of Materials, 2014, vol.69, no. 1, pp. 116-131. doi: 10.1016/ j.mechmat.2013.09.019.

36. Khohlov A.V. Properties of creep curves at piecewise-constant stress generated by the linear viscoelasticity theory, Problemy prochnosti i plastichnosti [Problems of strength and plasticity], 2015, vol.77, no. 4, pp. 344-359 (In Russian).

Analysis of creep curves produced by the linear viscoelasticity theory. . .

37. Khohlov A.V. The Qualitative Analysis of Theoretic Curves Generated by Linear Viscoelasticity Constitutive Equation, Science and Education, 2016, no. 5, pp. 187-245 (In Russian), http://technomag.bmstu.ru/doc/840650.html.

38. Khokhlov A. V. The nonlinear Maxwell-type viscoelastoplastic model: properties of creep curves at piecewise-constant stress and criterion for plastic strain accumulation, Mashinos-troenie i inzhenernoe obrazovanie, 2016, no. 3, pp. 55-68 (In Russian).

39. Khokhlov A. V. Asymptotic commutativity of creep curves at piecewise-constant stress produced by the linear viscoelasticity theory, Mashinostroenie i inzhenernoe obrazovanie, 2016, no. 1, pp. 70-82 (In Russian).

40. Khokhlov A. V. Long-term strength curves generated by the nonlinear Maxwell-type model for viscoelastoplastic materials and the linear damage rule under step loading, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2016, vol.20, no. 3, pp. 524-543 (In Russian). doi: 10.14498/vsgtu1512.

41. Fatemi A., Yang L. Cumulative fatigue damage and life prediction theories: A survey of the state of the art for homogeneous materials, Int. J. Fatigue, 1998, vol. 20, no. 1, pp. 9-34. doi:10.1016/S0142-1123(97)00081-9.

42. Golub V. P. Some creep effects in cyclic loading, Strength of Materials, 1987, vol. 19, no. 5, pp. 605-610. doi: 10.1007/bf01524289.

43. Radchenko V. P., Kichaev E. K., Simonov A. V. Energetic version of the model of rheological deformation and destruction of metal under a joint action of static and cyclic loads, J. Appl. Mech. Tech. Phys., 2000, vol.41, no. 3, pp. 531-537. doi: 10.1007/BF02465307.

44. Zheng X.-T., Xuan F.-Z., Zhao P. Ratcheting-creep interaction of advanced 9-12% chromium ferrite steel with anelastic effect, Int. J. Fatigue, 2011, vol.33, no. 9, pp. 12861291. doi: 10.1016/j.ijfatigue.2011.04.009.

45. Barenblatt G. I., Kozyrev Iu. I., Malinin N. N., Pavlov D. Ia., Shesterikov S. A. On the vibrocreep of polymeric materials, Prikl. Mekh. Tekh. Fiz., 1965, no. 5, pp. 68-75 (In Russian).

46. Lokoshchenko A. M., Shesterikov S. A. On Vibrocreep, Inzhenernyi zhurnal. Mekhanika tverdogo tela, 1966, no. 3, pp. 141-143 (In Russian).

47. Lokoshchenko A. M. Vibrocreep of metals in uniaxial and complex stress states, Mechanics of Solids, 2014, vol.49, no. 4, pp. 453-460. doi: 10.3103/S0025654414040104.

48. Khokhlov A. V. Specific features of stress-strain curves at constant stress rate or strain rate yielding from linear viscoelasticity, Problemy prochnosti i plastichnosti [Problems of strength and plasticity], 2015, vol. 77, no. 2, pp. 139-154 (In Russian).

49. Khokhlov A. V. Long-term strength curves produced by linear viscoelasticity theory combined with failure criteria accounting for strain history, Trudy MAI, 2016, no. 91, pp. 1-32 (In Russian), http://www.mai.ru/science/trudy/published.php?ID=75559.

50. Khokhlov A. V. Khokhlov A.V. Analysis of general properties of creep curves generated by the Rabotnov nonlinear hereditary relation under multi-step loadings, Vestn. Mosk. Gos. Tekh. Univ. im. N. E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Technical University. Series Natural Sciences], 2017, no. 3, pp. 93-123 (In Russian). doi:10.18698/1812-3368-2017-3-93-123.

51. Radchenko V. P., Samarin Yu. P. Effect of creep on the elastic deformation of a laminar composite, Mechanics of Composite Materials, 1983, vol.19, no. 2, pp. 162-168. doi: 10. 1007/BF00604218.

52. Radchenko V. P., Shapievskii D. V. Drift of elastic deformation due to creep for nonlinear elastic materials, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2006, no. 43, pp. 99-106 (In Russian). doi: 10.14498/ vsgtu458.

53. Radchenko V. P., Shapievskii D. V. Mathematical model of creep for a microinhomogeneous nonlinearly elastic material, J. Appl. Mech. Tech. Phys., 2008, vol.49, no. 3, pp. 478-483. doi:10.1007/s10808-008-0064-9.

54. Melnis A. É., Laizan Ya. B. Nonlinear creep of human compact bone tissue upon stretching, Polymer Mechanics, 1978, vol. 14, no. 1, pp. 82-84. doi: 10.1007/BF00859565.