CC BY
65
ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1512 УДК 539.378
КРИВЫЕ ДЛИТЕЛЬНОЙ ПРОЧНОСТИ НЕЛИНЕЙНОЙ МОДЕЛИ ВЯЗКОУПРУГОПЛАСТИЧНОСТИ ТИПА МАКСВЕЛЛА И ПРАВИЛО СУММИРОВАНИЯ ПОВРЕЖДЕННОСТИ ПРИ СТУПЕНЧАТЫХ НАГРУЖЕНИЯХ
А. В. Хохлов
Московский государственный университет имени М. В. Ломоносова, Научно-исследовательский институт механики, Россия, 119192, Москва, Мичуринский проспект, 1.
Аннотация
Исследуется нелинейное определяющее соотношение типа Максвелла для вязкоупругопластичных материалов с двумя произвольными материальными функциями одного аргумента с целью выявления арсенала его возможностей и области применимости, в частности, возможности подключения к нему критериев для описания разрушения при ползучести с постоянным или кусочно-постоянным напряжением. При минимальных первичных ограничениях на материальные функции выведены уравнения теоретических кривых длительной прочности, порожденных соотношением типа Максвелла в сочетании с деформационными и энергетическим (диссипативным) критериями разрушения, а также с интегральными критериями, учитывающими историю деформирования. Аналитически изучены и сопоставлены общие свойства этих кривых. Показано, что они адекватно описывают типичные свойства экспериментальных кривых длительной прочности вязкоупругопластиче-ских материалов. Для найденных зависимостей времен разрушения от напряжения проверено выполнение правила линейного накопления по-врежденности при ступенчатом нагружении и выведены формулы для величины (и знака) отклонений от него. В частности доказано, что для нелинейного соотношения типа Максвелла в сочетании с диссипативным критерием разрушения всегда в точности выполняется правило линейного накопления поврежденности при любом ступенчатом нагружении, а для деформационного критерия оно, наоборот, никогда не выполняется, но дает оценку сверху или снизу для времени разрушения в зависимости от знака скачка напряжения.
Ключевые слова: вязкоупругопластичность, кривые ползучести, повре-жденность, критерии разрушения, диссипация, время разрушения, кривые длительной прочности, сверхпластичность.
© 2016 Самарский государственный технический университет. Образец для цитирования
Хохлов А. В. Кривые длительной прочности нелинейной модели вязкоупругопластично-сти типа Максвелла и правило суммирования поврежденности при ступенчатых нагруже-ниях // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. Т. 20, № 3. С. 524-543. аок 10.14498/vsgtu1512. Сведения об авторах
Андрей Владимирович Хохлов (к.т.н; andrey-khokhlov@ya.ru), старший научный сотрудник, лаб. упругости и пластичности.
Введение. Надежное моделирование кривых длительной прочности (КДП) материала или элемента конструкции по результатам испытаний на ползучесть при высоких напряжениях необходимо для оценки срока безопасной эксплуатации при низких напряжениях, когда время разрушения слишком велико, чтобы можно было получить экспериментальную информацию о нем.
Для прогнозирования длительной прочности материалов к выбранному (построенному) определяющему соотношению (ОС) необходимо добавить критерий разрушения (КР), задающий момент разрушения £* (обрыва кривых ползучести) по достижению критических значений некоторой меры повре-жденности (скалярной, векторной или тензорной [1-7]). Желательно, чтобы КР и ОС хорошо взаимодействовали, то есть позволяли вывести общее уравнение теоретической кривой ползучести е = е(£; а, Т) и КДП £* = /(а, Т) (или а = Г(£*,Т)), где £* — время разрушения при данном уровне напряжения а и температуре Т, аналитически исследовать в общем виде зависимость свойств КДП от материальных функций и параметров ОС и КР и указать ограничения на них, обеспечивающие совпадение качественных свойств теоретических КДП с типичными свойствами экспериментальных КДП.
Эта (идеальная) программа реализована в работах [8-11 и др.] для линейного ОС вязкоупругости, нелинейного ОС Ю. Н. Работнова и двух нелинейных ОС, построенных и исследованных в [8, 9] (каждое из перечисленных нелинейных ОС содержит две материальные функции одного аргумента в одноосном случае и несколько параметров). Эти ОС нацелены на описание комплекса основных механических эффектов, типичных для вязкоупругопла-стичных материалов, обладающих выраженной наследственностью и высокой чувствительностью к скорости деформирования (полимеры, твердое топливо, асфальтобетон, композиты, пены, титановые и алюминиевые сплавы, углеродные и керамические материалы при высоких температурах, связки, сухожилия и другие биологические ткани и т.п.).
В данной работе эта программа реализована для нелинейного ОС типа Максвелла [12]:
е(£) = ее + е^, ее = Г(а)/Е, е^ = V(а)/п, t> 0. (1)
Оно связывает истории напряжения а(£) и деформации е(£) в точке тела и содержит две материальные функции (МФ) Г(х) и V(х) и две постоянные Е > 0, п > 0 (см. подробнее п. 1). Ниже будут выведены и аналитически изучены уравнения теоретических кривых длительной прочности, порожденных ОС (1) в сочетании с деформационными и энергетическим КР, а также — с интегральными КР, учитывающими историю деформирования [10].
Классический деформационный критерий разрушения (ДКР) постулирует, что разрушение происходит в момент когда деформация е(£, а) достигает критического значения е*:
е(£*,а) = е*. (2)
Он может описывать две разные ситуации: разрушение материала или элемента конструкции. Соответственно, е* обозначает либо материальную константу (зависящую только от температуры), либо предельную допустимую (по конструктивным соображениям) деформацию, рассчитанную по предель-
ному перемещению или напряжению (классический пример — разрушение турбины вследствие ползучести ее лопаток). В первом случае ДКР с постоянным е* = е* (Т) не всегда адекватен данным испытаний материалов на ползучесть: у многих материалов е* зависит от напряжения [1-7] (этот факт можно учесть, задавая е* = е*(а) в (2)), причем эта зависимость может быть немонотонной [3,4,13,14]. Во втором случае ДКР (2) адекватен «по определению», если конструктивный критерий срабатывает раньше, чем происходит разрушение материала.
В работе [10] предложены и исследованы два параметрических семейства скалярных критериев разрушения при монотонном одноосном деформировании, обобщающих ДКР (2), но учитывающих историю нарастания деформации и зависящих от характеристик кривых ползучести (они хранят информацию об условиях развития деформации и процессах, приводящих к разрушению). Их строение мотивировано желанием иметь арсенал более чутких и допускающих регулировку КР, позволяющих точнее описывать экспериментальные КДП различных материалов (или элементов конструкций) и учитывающих информацию, которой пренебрегает ДКР. Каждое из двух построенных в [10] семейств КР образует монотонную и непрерывную шкалу критериев (монотонно и непрерывно зависящих от параметра и > 0), включающую ДКР как предельный случай. При этом критерии первого семейства при всех и > 0 дают большее время разрушения, чем ДКР, критерии второго семейства при и > 1 — меньшее, а различие можно сделать сколь угодно малым за счет выбора значений управляющего параметра у края шкалы (и ^ <х или и = 1).
Помимо КР, выраженных через меры поврежденности деформационного типа и их критические значения, уже полвека разрабатывается класс энергетических мер поврежденности и критериев разрушения, — например, в работах О. В. Соснина, А. Ф. Никитенко, Б. В. Горева, В. В. Федорова, В. П. Рад-ченко [4, 14, 15-21 и др.] В частности в цикле нескольких десятков работ О. В. Соснина и его коллег предложена и экспериментально исследована (для титановых и алюминиевых сплавов) мера поврежденности, равная «удельной работе рассеяния» (то есть, удельной диссипации в точке тела), и КР, постулирующий, что разрушение материала наступает при достижении диссипацией Ш(£, а) критического значения Ш* = Ш*(Т) [15-18]:
W (и,а) = Ш*. (3)
Основная цель данной статьи — аналитически исследовать общие свойства КДП ¿*(а), порожденных ОС (1) в сочетании с КР (2) и (3), и доказать, что они совпадают с качественными свойствами типичных экспериментальных КДП вязкоупругопластических материалов. А именно: данные испытаний показывают, что КДП ¿*(а) всегда убывают, ¿*(а) ^ при а ^ а0 и ¿*(а) ^ 0 при а ^ а*, где ао ^ 0 — порог ползучести, а* > 0 — предел мгновенной прочности а*. КДП многих материалов в координатах «^а — ^¿*» представляются отрезком прямой линии с угловым коэффициентом ^ < 0 или двухзвенной ломаной. Это означает, что регистрируемые значения ¿*(а&) хорошо аппроксимируются (при достаточно больших напряжениях) степенной функцией
¿* = С а!1,
где ^ < 0, С > 0, а > ао, или двумя степенными функциями с различными показателями. Именно такими, как будет доказано, получаются и КДП, порожденные ОС (1) с указанными КР, при найденных ограничениях на МФ.
Вторая задача анализа — проверить, выполняется ли для теоретических КДП ОС (1) (для каких КР, при каких условиях на МФ) «закон» линейного накопления поврежденности при ступенчатых нагружениях, т.е. условие разрушения
5 = 1,
где 5 — сумма парциальных времен:
^__t*(í, а) £г-1 + ^^ £т- 1 (4)
' £* (аг) £*(ат) ,
\ / т=1 4 '
г — число ступеней, а = {(ат,£т) | т = 1,...,г}, ат — уровни напряжений на интервалах (£т-1; £т) (£0 = 0, > £г-1)), £*(г, а) —время разрушения на ступени с а = аг; £*(аг) — время разрушения при ползучести с постоянным напряжением а^. В многочисленных испытаниях разных материалов (сталей, алюминиевых сплавов, композитов и т.п.) на ползучесть до разрушения при двухступенчатом нагружении (г = 2) обнаружены значительные отклонения от условия 5 = 1 и два основных типа поведения материала [3, с. 117-119]:
1) Б > 1 для любых а1, а2 > 0;
2) Б > 1 при а2 < а1 и Б < 1 при а2 > а1.
В данной работе выведены общие формулы для отклонения Б — 1, порожденные ОС (1) в сочетании с разными КР, и исследован его знак в зависимости от знака разности а2 — а1 .
1. Нелинейная модель типа Максвелла. Будем рассматривать изотермические одномерные процессы, характеризуемые в точке тела историей напряжения а(£) и (логарифмической) деформации е(£), £ > 0. Связь между процессами а(£) и е(£) зададим по аналогии с реологической моделью Максвелла (последовательным соединением линейных упругого и вязкого элементов), т.е. постулируем, что деформация е(£) представима суммой упругой и вяз-копластической компонент, каждая из которых зависит от (безразмерного) напряжения а(£), но нелинейно. Это приводит к ОС (1); его можно записать в интегральной и дифференциальной формах:
е(£) = Е-1Г(а(£)) + п-1 / V(а(т))Жг.
Jо
г
V(а(т)) ¿т,
0 (5)
или ( )
е = Е-1(Г'(а)а + тг-^(а)), £> 0,
где тг = п/Е — время релаксации линейной модели Максвелла. ОС (5) содержит две МФ Г(х) и V(х), х € (ш_; ), ш_ < 0, > 0, и две постоянные: «модуль упругости» Е > 0 и коэффициент вязкости п > 0. Параметры ш+ и могут быть интерпретированы как пределы прочности при растяжении и сжатии. Параметры Е и п выделены из МФ Г и V для удобства сопоставления с линейной моделью Максвелла (получающейся при V(х) = Г(х) = х) и учета влияния температуры в форме Е = Е(Т), п = п(Т) [22,23]. Процессы а(£) предполагаются кусочно-непрерывными и кусочно-гладкими (при
t < 0 считаем, что а = 0). Обезразмеривание напряжения можно производить делением на cE, c > 0, или на характерное напряжение материала (предел упругости, ползучести, текучести, прочности и т.п.). Безразмерное время вводится делением на характерное время, например, на тг при фиксированной температуре (при таком масштабировании тг = 1 и n = E во всех последующих формулах при данной T).
МФ F определяет в (1) упругую деформацию £е(а). Поэтому первичные ограничения на F(ж) (естественные с точки зрения феноменологии и минимальные математические [23]) таковы: F(ж) —непрерывная кусочно-гладкая строго возрастающая функция на интервале (w_; ш+), такая, что F(0) = 0 (тогда xF(ж) > 0, т.е. sgnF(ж) = sgnж). Последние два условия обеспечивают совпадение знаков упругой деформации ее(а) и а и соблюдение условия £e(0) = 0. Из строгого возрастания F следует возрастание £е(а) и существование обратной функции f = F-1 на промежутке (у;у), где
у := sup F(ж) = F(w+ - 0), у := inf F(ж) = F(w_ + 0).
Функция вязкости V(ж)/п в ОС (1) регулирует чувствительность напряжения к скорости деформации, наследственные свойства, скорость диссипации, ползучести и накопления пластической деформации [22,23]. Чем больше |V(а)|/п, тем меньше вязкое сопротивление (больше |ev| и |ё| при том же а) и тем ближе поведение моделируемого материала к поведению жидкости. Минимальные первичные ограничения на V(ж): V(ж) — непрерывная (нестрого) возрастающая функция на интервале (w_; ш+), такая, что V(0) = 0 (тогда жV(ж) ^ 0).
Исследование кривых релаксации, ползучести и деформирования, порожденных ОС (5), показывает [22,23], что следует различать два основных случая, в которых ОС (5) (моделируемый материал) ведет себя по-разному:
1) |V(ж) | > 0 при ж = 0;
2) V(ж) = 0 на некотором отрезке [а_; а+ ] С (w_; ш+), а_ ^ 0, а+ ^ 0, а+ = а_.
Во втором случае при а £ [а_; а+ ] ОС (1) моделирует (нелинейно) упругое поведение материала, а- , а+ — пределы упругости материала при сжатии и растяжении (и пороги ползучести), а при а > а+ (или а < а-) начинают проявляться диссипативные и вязкопластические свойства. Для материалов с одинаковым поведением при растяжении и сжатии обе МФ в (5) должны быть нечетными.
Наложенные на МФ ОС (5) первичные ограничения обеспечивают, в частности, термодинамическую согласованность модели, т.е. положительность работы напряжений в любом процессе деформирования и неотрицательность диссипации (соблюдение диссипативного неравенства). Действительно, работа напряжений а(т) в процессе деформирования е(т), связанном с а(т) соотношением (5), выражается следующей формулой:
A = f а(т)ё(т)dr = f E-1aF'(a)adr + f kV(a)adr = U + W, J 0 J 0 J 0
U := E-1 f ж^(ж)^ж, W = n-1 f а(т)V(а(т))dr, (6)
00
где и = и (а) — энергия упругой деформации, Ш = Ш [£; а(т)] — диссипация. Энергия и (а) положительна и возрастает при а = 0 в силу ограничения Г'(х) > 0 (и'(а) = Е-1аГ'(а) > 0 при а > 0). Скорость диссипации выражается формулой
Ж(*) = п-1а(£) V (а(£)).
В силу непрерывности МФ V(х) и ограничения xV(х) ^ 0 при х = 0 справедливы неравенства ^ 0 и Ш(£) ^ 0 при всех £ > 0. Равенство Ш(£0) = 0 возможно лишь в случае, когда V(х) = 0 на некотором отрезке [а_; а+ ] С (ш_; ш+) и а(т) € [а_; а+ ] для всех т < £0 (т.е. модуль напряжения не превышает пределов упругости).
Базовые свойства нелинейного оператора
П : а(£) ^ е(£)
из ОС (5), вытекающие из общих ограничений на МФ Г и V, и результаты исследования свойств теоретических квазистатических кривых (ползучести, релаксации, ступенчатой ползучести, деформирования при постоянных и кусочно-постоянных скоростях нагружения и деформирования и др.), порождаемых ОС (5) с произвольными МФ Г, V, приведены в [12,22,23].
Нелинейное ОС (5) нацелено на описание комплекса основных реологических эффектов, типичных для материалов, обладающих наследственностью и высокой чувствительностью к скорости деформирования (таких как полимеры, их расплавы и растворы, композиты, твердое топливо, асфальтобетон, титановые и алюминиевые сплавы, углеродные и керамические материалы при высоких температурах и др.); материалов, имеющих выраженную стадию установившейся ползучести, «площадку текучести» на диаграмме деформирования и предел текучести, зависящий от скорости деформирования; материалов, проявляющих (в определенных состояниях и при определенных режимах деформирования) свойства как твердого тела, так и жидкости. В частности, ОС (5) (и его модификации) могут быть полезны для моделирования материалов в режимах сверхпластического деформирования [23].
Общая тензорная формулировка нелинейных ОС максвелловского типа для (больших деформаций) вязкоупругих сред, родственных ОС (5), описание кинематики, термодинамические аспекты и способы конкретизации ОС изучались в работах [24-29]. Внимание авторов было сосредоточено на описании поведения жидкостей, они обсуждали эксперименты и эффекты, присущие жидкостям (расплавам и растворам полимеров и т.п.), не рассматривали кривые ползучести, релаксации и деформирования, порождаемые ОС, не задавали многих вопросов, специфичных для механики деформируемого твердого тела.
Зафиксировав Г(х) = х и V(х) = х|х|п-1, п > 1, в (5), получим трехпара-метрическую модель с линейной упругостью и степенной вязкостью (будем ее называть полулинейной):
е = Е-1а + п-1а|а|п-1. (7)
Модель (7) применялась, например, в работах [30,31]. В статье [32] исследовались кривые ползучести и обратной ползучести параллельного соединения двух таких моделей с различными показателями п. Более подробный обзор
литературы и родственных моделей (в теории ползучести, сверхпластичности и реологии полимеров) приведен в [22,23].
2. Кривые ползучести при постоянном и кусочно-постоянном напряжении. При a(t) = а = const, t > 0, ОС (5) дает уравнение семейства кривых ползучести (КП):
¿(t,a) = kV(a)t + E-1F(а) или ¿(t, а) = E-1 (V(a)rr-1t + F(а)) , (8)
где k := n-1, Tr := n/E. Если V(ж) = 0 на некотором отрезке [а_; а+ ], а_ < 0, а+ > 0, то при а £ [а_; а+ ] ОС (1) моделирует (нелинейно) упругое поведение материала и при а £ [а_; а+ ] ползучесть отсутствует (пределы упругости на сжатие и растяжение совпадают с пределами ползучести). Если V(а) = 0, то все КП линейны по времени при t > 0, т.е. при любых МФ ОС (5) моделирует только ползучесть с постоянной скоростью ¿(t) = kV(а). ОС (5) неспособно описывать стадии замедленной и ускоренной ползучести, а также ограниченную ползучесть (свойственную, например, многим полимерам). Так как V(ж) > 0 при ж > 0 и V(ж) возрастает, то КП (8) возрастает по t (при <г > 0) и по а, что совпадает с типичными качественными свойствами экспериментальных КП структурно стабильных однородных материалов.
Выраженная стадия установившейся ползучести характерна для многих пластичных металлов (в частности, в состоянии сверхпластичности) и полимеров в вязкотекучем состоянии.
Любую программу нагружения из N ступеней
N-1
a(t) = ^ а [h(t - ti-1) - h(t - ti)] + gnh(t - tN-1),
i=1
где ai £ (w_; w+) и tj > 0 (to = 0, tj > ti-1), оператор П, задающий ОС (5), переводит в сумму откликов на каждую ступеньку [22]:
¿(t) = E-1 [V(<Tn)rr-1(t - tN-1) + F(<n)] h(t - tN-1) +
N-1
+ ^ S(t - ti-1; ti - ti-1,<i). (9)
¿=1
Здесь
S(t; ts, а) := E-1 [V(a)rr-1t + F(а)] [h(t) - h(t - ts)] + E-1V(a)tsTr-1h(t - ts)
— отклик на импульс напряжения а с носителем [0; ts] (задающий кривую обратной ползучести [22]: N = 2, а1 = а, а2 = 0), S(t - ti-1; ti - ti-1,ai) —
его сдвиги по времени на ti_1. Так как все слагаемые кусочно линейны, то
КП (9) кусочно-линейна и имеет на интервалах (ti_1; ti), i = 1,..., N, вид
¿(t) = E-1 [V(<i)rr-1(t - ti-1) + F(<i)] + pi-1, t £ (ti-1; ti), (10) i- 1
Pi-1 := E-Ч-1 ^ V(ffm)(tm - tm-1), i = 2,..., N. (11)
m=1
Накопленная пластическая (необратимая) деформация (11) выражает влияние предыдущих (г — 1) ступеней нагружения (ро := 0). Отклонение КП (10) при £ > _1 от обычной КП ст(£) = ам — _1) при мгновенном нагру-жении в момент £ = ¿м_1 (с нулевой предысторией) равно постоянной _1. В общем случае грм_1 = 0, следовательно, ОС (5) не обладает свойством затухающей памяти. Скачки деформации в точках £ = порожденные скачками напряжения ст^+1—стч, равны [Е(04+1) — Е(ст^)] /Е, скачки скорости ползучести е(£) равны кУ(стг+1) — (ст^). Общие свойства семейства КП (10) и условия моделирования рэтчетинга изучены в [22].
3. Кривые длительной прочности для деформационных критериев разрушения. Если к ОС (5) с произвольным пределом упругости ст+ ^ 0 добавить деформационный КР (2), то время разрушения ¿* при ползучести с постоянным напряжением (ниже будем для краткости обозначать постоянный уровень напряжения СТ через ст), таким, что ст > ст+ и Е(ст) < Ее*, находится из (8) и (2):
П_1У (ст)** + Е _1Е (ст) = е*,
т.е. КДП имеет вид
¿*(ст) = гг [Ее* — Е(ст)]/V(ст), ст е (ст+; /(Ее*)) , (12)
где гг = п/Е, / = Е _1.
При ст е [0; ст+ ] ползучести нет, при Е (ст) ^ Ее* разрушение происходит в момент £ = 0. КДП (12) определена при ст е (ст+; /(Ее*)) и убывает (так как Е и V возрастают), ¿* ^ при ст ^ ст+ +0 и ¿* ^ 0 при ст ^ /(Ее*) — 0 (т.е. свойства КДП (12) совпадают с типичными качественными свойствами экспериментальных КДП).
Время разрушения ¿*(г, ст) при ползучести с кусочно-постоянным напряжением, если оно происходит на г-той ступени нагружения (тогда Стi > ст+, иначе деформация постоянна и разрушения быть не может), находится из (10) и (2):
кУ(ст^ (¿*(г, ст) — ^_1) + Е(стi)E 1 + = е*, ¿*(г, ст) — £»_1 = т_1 [Ее* — Е(стi) — Ер^] /V(стi). (13)
Проверим соблюдение «закона» линейного накопления поврежденности при ступенчатом нагружении, т.е. условия разрушения 5 = 1, где 5 — сумма парциальных времен (4) с учетом (12). В многочисленных испытаниях различных материалов на ползучесть до разрушения при двухступенчатом нагружении обнаружены значительные отклонения от условия 5 = 1 и два основных типа поведения материала [3, с. 117-119]:
1) Б > 1 для любых ст1, ст2 > 0;
2) Б > 1 при ст2 < ст1 и Б < 1 при ст2 > ст1.
Если разрушение происходит на второй ступени нагружения (тогда ст2 > ст+), то при г = 2 из (13) находим
¿*(ст) — ¿1 = п [е* — Е_1Е(ст2) — Р1] /V(ст2), Р1 = кУ(ст^ь а по (4) —
5 := ¿1/£*(ст 1) + (¿*(ст) — ¿1)/£*(ст2),
т.е.
_ = + П [в* - Е-1Е(02) - кУ(01^1] =
= П [в* - Е-1Е(01)] + п [в* - Е-1Е(02)] =
= -,+ ¿1У(01) [Е(01) - Е(02)] ( )
+ Тг [Ев* - Е(01)] [Ев* - Е(02)]' ( )
Отклонение 5 - 1 пропорционально разности Е(01) - Е(02) и пластической деформации р1, накопленной к моменту £ = ¿1. Очевидно, 5 = 1 тогда и только тогда, когда V(01) = 0, т.е. это возможно лишь для ОС с 0+ > 0, когда 01 < 0+ (на первой ступени отсутствуют неупругие деформации). Так как Ев* - Е(0г) > 0 (в противном случае разрушение происходит сразу в момент приложения нагрузки: при £ = 0 или £ = ¿1), из (14) следует, что 5 > 1 при 02 <01 и 5 < 1 при 02 >01 (и 01 > 0+). Т.е. правило линейного суммирования поврежденностей 5 = 1 для ОС (5) в сочетании с ДКР (2) никогда не соблюдается (за исключением тривиального случая 01 < 0+), отклонение от него возможно в обе стороны в зависимости от знака 02 - 01 (при 02 < 01 оно дает для времени разрушения оценку снизу, а при 02 > 01 — сверху). Обнаруженное свойство — полезный индикатор применимости ОС (5) в сочетании с ДКР для описания длительной прочности: такая модель не может описать разрушение тех материалов, у которых всегда 5 > 1 при любом знаке 02 - 01.
Если, применяя ДКР (2), пренебречь упругой деформацией и взять в качестве меры поврежденности только вязкопластическую компоненту деформации (деформацию ползучести), то (положив Е = 0 в (12)) получим время разрушения
¿*(0) = ТгЕв*/У(0), 0 е (0+;), 0+ ^ 0' (15)
Очевидно,
¿*(0)/£*(0) = 1 - Е(0)(Ев*)-1, £*(0) > ¿*(0) и ¿*(0) > 0,
т.е. ¿*(0) = 0 при конечном напряжении (в отличие от (12)). Формула (14) при Е = 0 превратится в
5 = 1 + Е-2в-2тг-1£1 V (01),
т.е. 5 > 1 для любых 01, 02 > 0+ (тогда V(01) > 0), и знак отклонения 5 - 1 от правила линейного суммирования в этом случае уже не зависит от знака разности 02 - 01.
4. Время разрушения для критериев, учитывающих историю деформирования. В работе [10] предложены и исследованы два параметрических семейства критериев разрушения при монотонном одноосном деформировании, обобщающих ДКР (2):
и) = в*, и > 0. (16)
В качестве меры поврежденности и) предлагается использовать одно из двух семейств средних значений деформации за время ¿, зависящих от параметра и > 0:
г г \ 1/и
*-1у |в(т, ¿> 0, (17)
ft \ 1/u t-1j( ^(т. (18)
В основе КР (16), (18) лежит идея определения средней деформации через среднюю скорость vu(t) на интервале (0; t). Целесообразность использования критерия (16), (18) связана с тем, что он обобщает и формализует эмпирический факт [2]: в испытаниях на ползучесть для разных материалов произведение t* на скорость установившейся ползучести v(a) — величина постоянная, не зависящая (или слабо зависящая) от напряжения: vt* = const. КР (16), (18) обобщает это наблюдение на неустановившуюся ползучесть, а при ^ = v = const дает ¿u(t) = ¿(t) и vt* = const. Критическое значение поврежденности w* = ¿* в (16) может быть константой (зависящей только от температуры) или функцией напряжения.
Детальное аналитическое исследование КР (16), (17) и (16), (18) позволило обнаружить ряд полезных свойств, подтверждающих их пригодность для описания разрушения материалов при ползучести и прогнозирования длительной прочности. В работе [10] доказано, что для любого неубывающего процесса деформации ¿(t) справедливы следующие утверждения:
1) средние ¿u(t), u > 0, и ¿u(t), u ^ 1, возрастают по t и потому их можно использовать как меру поврежденности в КР для монотонных процессов деформирования (требование возрастания w(t) формализует представление о том, что процессы разрушения и нарастания деформации ползучести, связанные общими структурными механизмами, «сонаправ-лены»);
2) скорости роста поврежденностей (17) и (18) выражаются следующими формулами:
4(t) = ¿u(t)(ut)-1 [(¿(t)/£u(t))u - 1] , ¿u(t) = vu(t)u-1 [u - 1 + (¿(t)/vu(t))u] , t > 0;
3) семейства функций ¿u (t) и ¿u(t) непрерывны и возрастают по параметру u > 0;
4) для всех t > 0 ¿u(t) < ¿(t) при u > 0 и ¿u(t) ^ ¿(t) при u > 1 (¿u(t) = ¿(t) лишь при ^ = const);
5) ¿u(t) ^ ¿(t) при u ^ то, т.е. КР (16), (17) при достаточно большом u дает значение t*, сколь угодно близкое к ДКР (2); при u = 1 ¿1(t) = ¿(t), t > 0, и КР (16), (18) совпадает с ДКР (2).
Эти свойства показывают, что критерии (16), (17) образуют монотонную и непрерывную шкалу критериев (монотонно и непрерывно зависящих от параметра u > 0), более мягких, чем ДКР, но сколь угодно мало отличающихся от него при больших значениях u, а КР (16), (18) с u > 1 — шкалу критериев, более жестких, чем ДКР и включающую ДКР как предельный случай при u = 1. КР семейства (16), (17) при всех u > 0 дают большее время разрушения, чем ДКР, семейства (16), (18) при u > 1 — меньшее, а различие можно сделать сколь угодно малым за счет выбора значений параметра u у края шкалы (u ^ то или u = 1). Это свойство очень полезно для точной настройки модели на имеющиеся экспериментальные данные о зависимости времени разрушения от уровня напряжения, температуры и других факторов: если
w(t, u) = ¿u(t) := ¿(0) + Vu(t)t, Vu(t) :=
эти данные плохо описывает ДКР, то можно выбрать более подходящий КР из построенных семейств (16), (17) или (16), (18), плавно и монотонно изменяя значение и.
Так как скорость ползучести в(£) = кУ(0) не зависит от времени, средняя скорость равна = кУ(0), и (18) дает поврежденность
ш = ви(£) := в(0) +
совпадающую с полной деформацией (8). Таким образом, при любом и КР (16), (18) в сочетании с ОС (5) порождает такую же КДП, как и деформационный КР (2), т.е. КДП (12).
Подстановка КП (8) в (17) дает для поврежденности (17) выражение
£„(*) = Е-1 Е(0) (д(0)*/тг)-1/и [(д(0)*/тг + 1)ад+1 - 1]1/и (и + 1)-1/и,
д(ж) := V(ж)/Е(ж), а КР (16), (17), т.е. (£«(£*))" = в^, дает для времени разрушения ¿*(0) = тгТ*(0) уравнение:
(д(0)Т* + 1)"+1 - 1] (?(0)Т*)-1 = (и + 1)(Ев*)и(Е(0))-", 0> 0. (19)
Если пренебречь упругой составляющей деформации и применять КР (16), (17) к вязкопластической компоненте деформации (деформации ползучести) р(£) = ^(0)£, то формулы для поврежденности и времени разрушения упростятся:
ри (¿) = ¿V (0)(и + 1)-1/^,
г* = в*п(и + 1)1/^-1(0) или г* = тгЕв*(и + 1)^^(0)) 1, 0 е (0+;ш+).
(20)
Свойства КДП (20) совпадают со свойствами КДП (15), так как г*(0)Д*(0) = (и + 1)1/ад = с(и), и > 0; функция с(и) убывает, с(0+) = е, с(+то) = И с(и) = 1, и потому
1 < г* (0)Д*(0) < е.
5. Кривая длительной прочности для энергетического КР и сравнение
КДП. В цикле работ О. В. Соснина и его коллег [15-18] предложена и экспериментально исследована мера поврежденности, равная «удельной работе рассеяния» (то есть удельной диссипации в точке тела), и энергетический КР (3), постулирующий, что разрушение материала наступает при достижении диссипацией Ш(¿, 0) критического значения Ш* = Ш*(Т).
В силу (6) скорость диссипации при ползучести постоянна:
= fc0V (0) = 0г,
где г = fcV(0) — скорость установившейся ползучести. Так как
Ш (¿) = fc0V (0)£,
и
энергетический КР (ЭКР) (3) дает для времени разрушения уравнение
кстУ (ст)£* = Ж*,
т.е. порождает КДП
= (ст)] или £*(ст) = ТгЕЖ»/[стУ(ст)], ст е (ст+; ). (21)
КДП (21) убывает, как и (12), и £*(ст+ + 0) = (как в случае ст+ = 0, так и при ст+ > 0). КДП (21) отличается от КДП (12), порожденной ДКР, тем, что на нее не влияет МФ Е (упругая деформация), и тем, что £ *(ст) = 0 (£*(ст) > 0 и £*(ст) ^ 0 при ст ^ если = +то). Их поведение в правой окрестности предела упругости-ползучести ст = ст+ однотипно: при ст ^ ст+ +0 ¿*(ст) ^ и для (21), и для (12) (см. рис.), так как V(ст+) = 0.
Отношение времен разрушения (12) и (21) для ДКР и ЭКР следующее:
**(ст)Д*(ст) = [Ее* - Е(ст)] ст/(Ж*Е), ст е (ст+; /(Ее*)).
Если ст+ = 0, то оно обращается в нуль при ст = 0 и ст = /(Ее*), т.е. в окрестностях этих точек время разрушения по ДКР (2) значительно меньше, чем по ЭКР (3), а при ст ^ 0 абсолютное и относительное отклонения £ *(ст) — ¿*(ст) и [£*(ст) — ¿*(ст)]Д*(ст) стремятся к бесконечности. Отношение времен разрушения (15) (ДКР без учета упругой деформации) и (21) задается так:
¿*(ст)Д*(ст) = е*ст/Ж*, ст е (ст+;).
Если V = жга, то уравнения (21), (15) и (20) дают степенные зависимости времени разрушения от напряжения:
= пЖ* ст-п-1, ** = пе*ст-п, ** = пе*(и + 1)1/ист-п,
т.е. в логарифмических координатах ст—^ ¿*» эти КДП изображаются прямыми линиями с угловыми коэффициентами — (п + 1) и — п (см. правый рисунок), что типично для экспериментальных КДП многих материалов. Если МФ V склеена из двух степенных функций, определенных на [0; ж1 и [ж1; ж2], то зависимости £ *(ст), ¿*(ст) и ¿*(ст) тоже будут склеены из двух степенных функций вида с отрицательными показателями, и в координатах «^ ст — ^¿*» КДП будет изображаться двузвенными ломаными. Если ст+ > 0, V(ж) = 0 при ж < ст+ и V (ж) = (ж — ст+ )п при ж > ст+, то КДП (21) превращается в
= тгЕЖ*(ст — ст+ )-пст-1, ст е (ст+;).
На левом рисунке приведены КДП полулинейной модели (7) с ст+ = 0, V(ж) = ж|ж|, Е(ж) = ж, Е = 100, п = 100 (тогда тг = 1), е* = 0.02, Ж* = 0.02. Черная кривая —КДП (12) для ДКР, голубая —КДП (15) для ДКР без учета упругой деформации, красная —КДП (21) для ЭКР (3) (эта КДП растет быстрее остальных при ст ^ 0+), желтые кривые —КДП (20) с и = 0.1; 1; 10; 30 для интегрального КР (16), (17) без учета упругой деформации (всегда ¿*(ст, и) > ¿*(ст) > ¿*(ст)). Голубая штриховая линия — предельная кривая
для семейства КДП (20) при и ^ 0 (она получается из КДП (15) умножением на е); при и ^ семейство КДП (20) сходится к КДП (15). Синяя штрихпунктирная линия — КДП (19) (порожденная КР (16), (17) с учетом полной деформации): для выбранных МФ q(а) = а, и уравнение КДП (19), принимает вид
[(аТ* + 1)ад+1 - 1] (аТ*)-1 = (и + 1)(Ее*)иа-и;
для и = 1
¿* = 2тг(Ее* — а)а-2, а е (0; Ее*).
Все КДП имеют вертикальную асимптоту а = 0, так как а+ = 0. Четыре штриховые линии — КДП модели с положительным пределом упругости (ползучести) а+ = 1 и МФ V(х) = 0 при х < а+, V = (х — а+ )2 при х > а+ (по-прежнему Р(х) = х, х > 0). Цвета КДП соответствуют тем же КР, желтая штриховая — КДП (20) при и = 1. Все эти КДП имеют вертикальную асимптоту а = а+.
На правом рисунке приведены те же самые КДП, но в координатах а — ^ ¿*» (степенные КДП (15), (20) и (21) в них становятся прямолинейными). У КДП, которые обращались в нуль в точке а = Ее* = 2, появляется вертикальная асимптота х = ^(Ее*). Отметим, что так как
W (*) = kаV (а%
для ОС (1) в сочетании с ЭКР выполняется гипотеза единой кривой ползучести Никитенко—Соснина [15,16] в координатах «W — t/t*»: по (21)
х := = tаV(а)/(№*п),
и КП (8) имеет вид
W = W*ж,
Кривые длительной прочности полулинейной модели (7) с а+ =0, V = x2 и с а+ > 0, V = (x — а+ )2, x > а+ для разных критериев разрушения (онлайн в цвете) [Long-term strength curves of the semilinear model (7) for different failure criteria when a+ =0, V = x2, and a+ > 0, V = (x — a+ )2, x > a+ (color online)]
т.е. не зависит от напряжения. Для ДКР (2) и других рассмотренных КР эта гипотеза не выполняется, так как ¿*(ст) и х(ст) := ¿/¿* имеют иной вид: например, для (15) получим
х := ¿Д* = ¿V(а)/(пе*)
W = е*стх.
При ползучести с кусочно-постоянным напряжением скорость диссипации постоянна на каждом интервале (¿¿-1; ¿¿):
= А^У(ст) > 0, поэтому диссипация выражается следующей формулой:
¿-1
Ж(*)= V (^(¿т - ¿т-1) + кст^У- ¿¿-1), * е (¿¿-1; ¿¿). (22)
т=1
Время разрушения £*(г, ст) на ¿-той ступени нагружения находится из (22) и ЭКР (3):
**(г, ст) - ¿¿-1 = ЖпМУ(ст) - V ^(¿т - ¿т-1). (23)
Проверим, соблюдается ли «закон» линейного накопления поврежденно-сти при ступенчатом нагружении, т.е. условие разрушения 5 = 1, где сумма парциальных времен вычисляется по (4), а ¿*(стт) = ¿*(<тт) — по (21):
= V (стт)].
Подставим (23) в (4):
5 = __у"1 СТтУ (стт) (^ - ^ ) + > (¿т - ¿т-1)
т.е. в точности Б = 1 (ведь £ *(^)Д*(стт) = (^¿)). Таким обра-
зом, для ОС (1) (с произвольными МФ) в сочетании с ЭКР (3) всегда выполняется правило линейного накопления поврежденности при ступенчатых нагружениях, и критерий 5 = 1 дает в точности такое же время разрушения (23).
Заключение. Исследовано нелинейное определяющее соотношение типа Максвелла (1) для вязкоупругопластичных материалов с целью выявления его возможностей и области применимости. Одна материальная функция соотношения (1) управляет упругими свойствами, вторая — вязкопластиче-скими: регулирует наследственные свойства, вязкость, скорость диссипации, чувствительность напряжения к скорости деформации, зависимость скорости ползучести и релаксации от величины и знака напряжения, величину и скорость накопления необратимой деформации.
и
При минимальных первичных ограничениях на материальные функции выведены уравнения теоретических кривых длительной прочности (КДП), порожденных соотношением (1) в сочетании с деформационными и энергетическим критериями разрушения (2) и (3), а также — с интегральными критериями (16)—(18), учитывающими историю деформирования; аналитически изучены и сопоставлены общие свойства этих кривых, доказано, что они адекватно описывают типичные свойства экспериментальных КДП вязкоупруго-пластических материалов. Для найденных зависимостей времен разрушения от напряжения проверено выполнение правила линейного накопления повре-жденности при ступенчатом нагружении и выведены формулы для отклонения от него. В частности, доказано, что для определяющего соотношения (1) в сочетании с энергетическим КР (3) всегда выполняется правило линейного накопления поврежденности при ступенчатом нагружении, и условие разрушения S = 1 (см. (4)) дает точное время разрушения (23), а для соотношения (1) в сочетании с деформационным критерием разрушения (2) оно никогда не выполняется, но дает оценку сверху или снизу в зависимости от соотношения между уровнями напряжения на ступенях нагружения (для произвольных материальных функций). Если же, применяя деформационный критерий разрушения (2), пренебречь упругой деформацией, то знак отклонения S — 1 уже не будет зависеть от знака разности — о\ : будет выполняться неравенство S > 1 для любых о2 и а\ (превосходящих предел упругости а+ ^ 0).
Создание арсенала критериев разрушения, хорошо согласованных с построенным (или выбранным) определяющим соотношением, аналитическое исследование уравнений порожденных ими КДП в общем виде и их настройка за счет параметров, изменяющих КДП известным (в силу результатов анализа) образом, полезны для аттестации определяющего соотношения, обработки результатов испытаний, моделирования и прогнозирования длительной прочности вязкоупругопластических материалов и конструкций.
Декларация о финансовых и других взаимоотношениях. Исследование не имело спонсорской поддержки. Автор несет полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи была одобрена автором. Автор не получал гонорар за статью.
ORCID
Андрей Владимирович Хохлов: http://orcid.org/0000-0002-9212-2579
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Качанов Л. М. Теория ползучести. М.: Физматгиз, 1960. 456 с.
2. Работнов Ю. Н. Ползучесть элементов конструкций. М.: Наука, 1966. 752 с.
3. Локощенко А. М. Ползучесть и длительная прочность металлов. М.: Физматлит, 2016. 504 с.
4. Радченко В. П., Еремин Ю. А. Реологическое деформирование и разрушение материалов и элементов конструкции. М.: Машиностроение-1, 2004. 265 с.
5. Betten J. Creep Mechanics. Berlin, Heidelberg: Springer-Verlag, 2008. 367 pp.
6. Bergstrom J. S. Mechanics of Solid Polymers. Theory and Computational Modeling. Amsterdam: William Andrew is an imprint of Elsevier, 2015. 520 pp.
7. Махутов Н. А. Деформационные критерии разрушения и расчет элементов конструкций на прочность. М.: Машиностроение, 1981. 272 с.
8. Хохлов А. В. Определяющее соотношение для реологических процессов: свойства теоретических кривых ползучести и моделирование затухания памяти // Известия РАН. МТТ, 2007. №2. С. 147-166.
9. Хохлов А. В. Определяющее соотношение для реологических процессов c известной историей нагружения. Кривые ползучести и длительной прочности // Изв. РАН. МТТ, 2008. №2. С. 140-160.
10. Хохлов А. В. Критерии разрушения при ползучести, учитывающие историю деформирования, и моделирование длительной прочности // Изв. РАН. МТТ, 2009. №4. С. 121135.
11. Хохлов А. В. Общие свойства кривых ползучести и длительной прочности, порождаемых нелинейной теорией наследственности Ю. Н. Работнова: Отчет о НИР № 5288. М.: НИИ механики МГУ, 2015. 74 с.
12. Хохлов А. В. Свойства нелинейной модели вязкоупругопластичности типа Максвелла с двумя материальными функциями // Вестник МГУ. Сер. 1. Математика, механика, 2016. №6. С. 36-41.
13. Локощенко A. М., Шестериков С. А. К проблеме оценки длительной прочности при ступенчатом нагружении// ПМТФ, 1982. №2. С. 139-143.
14. Радченко В. П., Кичаев П. Е. Энергетическая концепция ползучести и виброползучести металлов. Самара: СамГТУ, 2011. 157 с.
15. Никитенко А. Ф., Соснин О. В. О разрушении при ползучести// ПМТФ, 1967. №3. С. 74-75.
16. Соснин О. В., Горев Б. В., Никитенко А. Ф. Энергетический вариант теории ползучести. Новосибирск: ИГиЛ СО АН СССР, 1986. 95 с.
17. Никитенко А. Ф. Ползучесть и длительная прочность металлических материалов. Новосибирск: НГАСУ, 1997. 278 с.
18. Соснин О. В., Никитенко А. Ф., Горев Б. В. К обоснованию энергетического варианта теории ползучести и длительной прочности металлов // ПМТФ, 2010. С. 188-197. 4 с.
19. Федоров В.В. Термодинамический метод оценки длительной прочности // Пробл. прочности, 1972. №9. С. 45-47.
20. Радченко В.П. Энергетический вариант одноосной теории ползучести и длительной прочности// ПМТФ, 1991. №4. С. 172-179.
21. Радченко В. П., Кичаев Е. К., Симонов А. В. Энергетический вариант модели реологического деформирования и разрушения металлов при совместном действии статических и циклических нагрузок // ПМТФ, 2000. №3. С. 169-175.
22. Хохлов А. В. Нелинейная модель вязкоупругопластичности типа Максвелла: свойства кривых ползучести при ступенчатых нагружениях и условия накопления пластической деформации // Машиностроение и инженерное образование, 2016. №3. С. 35-48.
23. Хохлов А. В. Нелинейные модели вязкоупругости типа Максвелла. Особенности их поведения, скоростная чувствительность и возможность использования для описания ползучести и сверхпластичности материалов: Отчет о НИР № 5193. М.: НИИ механики МГУ, 2013. 108 с.
24. Городцов В. А., Леонов А. И. О кинематике, неравновесной термодинамике и реологических соотношениях в нелинейной теории вязкоупругости // ПММ, 1968. Т. 32, №1. С. 70-94.
25. Leonov A. I., Lipkina E. Ch., Paskhin E. D., Prokunin A. N. Theoretical and experimental investigations of shearing in elastic polymer liquids// Rheol. Acta, 1976. vol. 15, no. 7/8. pp. 411-426. doi: 10.1007/BF01574496.
26. Пальмов В. А. Реологические модели в нелинейной механике деформируемых тел // Успехи механики, 1980. Т. 3, №3. С. 75-115.
27. Прокунин А. Н. О нелинейных определяющих соотношениях максвелловского типа для описания движения полимерных жидкостей // ПММ, 1984. Т. 48, №6. С. 957-965.
28. Leonov A. I., Prokunin A. N. Non-linear Phenomena in Flows of Viscoelastic Polymer Fluids. London: Chapman and Hall, 1994. xvii+475 pp. doi: 10.1007/978-94-011-1258-1.
29. Leonov A. I. Constitutive equations for viscoelastic liquids: Formulation, analysis and comparison with data// Rheology Series, 1999. vol. 8. pp. 519-575. doi: 10.1016/ S0169-3107(99)80040-9.
30. Cao Y. Determination of the creep exponent of a power-law creep solid using indentation tests// Mech. Time-Depend. Mater, 2007. vol.11, no. 2. pp. 159-172. doi: 10.1007/ s11043-007-9033-6.
31. Naumenko K., Altenbach H., Gorash Y. Creep Analysis with a Stress Range Dependent Constitutive Model// Arch. Appl. Mech., 2009. vol.79, no. 6. pp. 619-630. doi: 10.1007/ s00419-008-0287-5.
32. Радченко В. П., Шапиевский Д. В. Анализ нелинейной обобщенной модели Максвелла// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2005. №38. С. 55-64. doi:10.14498/vsgtu372.
Поступила в редакцию 15/VII/2016; в окончательном варианте — 04/IX/2016; принята в печать — 09/IX/2016.
Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki
[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 3, pp. 524-543
ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1512
MSC: 74R20, 74D99, 74C99
LONG-TERM STRENGTH CURVES GENERATED BY THE NONLINEAR MAXWELL-TYPE MODEL FOR
VISCOELASTOPLASTIC MATERIALS AND THE LINEAR DAMAGE RULE UNDER STEP LOADING
A. V. Khokhlov
Lomonosov Moscow State University,
Institute of Mechanics,
1, Michurinsky prospekt, Moscow, 119192, Russian Federation.
Abstract
The nonlinear Maxwell-type constitutive relation with two arbitrary material functions is formulated for viscoelastoplastic materials and studied analytically in uni-axial case to reveal capabilities of the model and its applicability scope. Its coupling with a number of fracture criteria is analyzed in order to simulate creep rupture under constant and piecewise-constant loading and to compare creep life estimates arising as a result. The limit strain criterion, the critical dissipation criterion and two proposed new families of failure criteria taking into account a strain history (i.e. a whole creep curve) are considered. Long-term strength curves equations generated by each one of the four chosen failure criteria are derived. Their general qualitative properties are analyzed and compared to each other under minimal restrictions on material
© 2016 Samara State Technical University. Please cite this article in press as:
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. doi: 10.14498/vsgtu1512. (In Russian) Authors Details:
Andrew V. Khokhlov (Cand. Techn. Sci.; andrey-khokhlov@ya.ru), Senior Researcher, Lab. of Elasticity and Plasticity.
functions of the constitutive relation. It is proved that qualitative properties of all theoretic long-term strength curves coincide with basic properties of typical test long-term strength curves of viscoelastoplastic materials. For every failure criteria considered herein, rapture time under step-wise loading is evaluated for arbitrary material functions and compared to the lifetime yielding from the linear damage accumulation rule (i.e. "Miner's rule"). General formulas for cumulative damage ("Miner's sum") deviations from unity are obtained for all failure criteria coupled with the nonlinear Maxwell-type constitutive relation. Their dependences on material functions and loading program parameters are examined. In particular, it is proved that the linear damage rule is exactly valid for the critical dissipation criterion whatever material functions, number of loading steps and stress levels are chosen. On the contrary, for the limit strain criterion, the linear damage rule is never valid for two-step loading and cumulative damage at rapture instant is greater or less than unity depending on the sign of stress jump.
Keywords: viscoelastoplasticity, creep curves, damage, failure criteria, dissipation, creep rupture, creep lifetime, long-term strength curves, superplas-ticity.
Declaration of Financial and Other Relationships. The research has not had any sponsorship. The author is absolutely responsible for submitting the final manuscript in print. The author has approved the final version of manuscript. The author has not received any fee for the article.
ORCID
Andrew V. Khokhlov: http://orcid.org/0000-0002-9212-2579
REFERENCES
1. Kachanov L. M. Teoriia polzuchesti [Creep theory]. Moscow, Fizmatgiz, 1960, 456 pp. (In Russian)
2. Rabotnov Yu. N. Creep problems in structural members. Amsterdam, London, North-Holland Publ. Co., 1969, xiv+822 pp.
3. Lokoshchenko A. M. Polzuchest' i dlitel'naia prochnost' metallov [Creep and Long-Term Strength of Metals]. Moscow, Fizmatlit, 2016, 504 pp. (In Russian)
4. Radchenko V. P., Eremin Yu. A. Reologicheskoe deformirovanie i razrushenie materialov i elementov konstruktsii [Rheological deformation and fracture of materials and structural elements]. Moscow, Mashinostroenie-1, 2004, 265 pp. (In Russian)
5. Betten J. Creep Mechanics. Berlin, Heidelberg, Springer-Verlag, 2008, 367 pp.
6. Bergstrom J. S. Mechanics of Solid Polymers. Theory and Computational Modeling. Amsterdam, William Andrew is an imprint of Elsevier, 2015, 520 pp.
7. Makhutov N. A. Deformatsionnye kriterii razrusheniia i raschet elementov konstruktsii na prochnost' [Deformation and fracture criteria for the calculation of structural elements of strength]. Moscow, Mashinostroenie, 1981, 272 pp. (In Russian)
8. Khokhlov A. V. Constitutive relation for rheological processes: Properties of theoretic creep curves and simulation of memory decay, Mech. Solids, 2007, vol.42, no. 2, pp. 291-306. doi: 10.3103/S0025654407020148.
9. Khokhlov A. V. Constitutive relation for rheological processes with known loading history. Creep and long-term strength curves, Mech. Solids, 2008, vol. 43, no. 2, pp. 283-299. doi: 10. 3103/S0025654408020155.
10. Fracture criteria under creep with strain history taken into account, and long-term strength modelling, Mech. Solids, 2009, vol.44, no. 4, pp. 596-607. doi: 10.3103/ S0025654409040104.
XOX.OB A. B.
11. Khokhlov A. V. The general properties of curves and creep rupture strength generated by the nonlinear theory of heredity by Yu. N. Rabotnov, Research report no. 5288. Moscow, Lomonosov Moscow State University, Institute of Mechanics, 2015, 74 pp. (In Russian)
12. Khokhlov A. V. Properties of nonlinear Maxwell-type models of the visco material with two functions, Vestnik MGU. Ser. 1. Matematika, mekhanika, 2016, no. 6, pp. 36-41 (In Russian).
13. Lokoshchenko A. M., Shesterikov S. A. Problem of estimating the creep strength under step loading, J. Appl. Mech. Tech. Phys., 1982, vol.23, no. 2, pp. 289-292. doi: 10.1007/ BF00911015.
14. Radchenko V. P., Kichaev P. E. Energeticheskaia kontseptsiia polzuchesti i vibropolzuch-esti metallov [Energy Concept of Creep and Vibrocreep of Metals]. Samara, Samara State Technical Univ., 2011, 157 pp. (In Russian)
15. Nikitenko A. F., Sosnin O. V. Fracture in creep, J. Appl. Mech. Tech. Phys., 1967, vol.8, no. 3, pp. 50-51. doi: 10.1007/BF00913209.
16. Sosnin O. V., Gorev B. V., Nikitenko A. F. Energeticheskii variant teorii polzuchesti [Energetic variant of the creep theory]. Novosibirsk, Inst. of Hydrodynamics, USSR Acad. of Sci., 1986, 95 pp. (In Russian)
17. Nikitenko A. F. Polzuchest' i dlitel'naia prochnost' metallicheskikh materialov [Creep and long-term strength of metallic materials]. Novosibirsk, Novosibirsk State University of architecture and Civil Engineering Publ., 1997, 278 pp. (In Russian)
18. Sosnin O. V., Nikitenko A. F., Gorev B. V. Justification of the energy variant of the theory of creep and long-term strength of metals, J. Appl. Mech. Tech. Phys., 2010, vol.51, no. 4, pp. 608-614. doi: 10.1007/s10808-010-0078-y.
19. Fedorov V. V. Thermodynamic method of estimating long-term strength, Strength Mater., 1972, vol.4, no. 9, pp. 1076-1079. doi: 10.1007/BF01529581.
20. Radchenko V. P. Energy variant of the uniaxial theory of creep and rupture strength, J. Appl. Mech. Tech. Phys., 1991, vol.32, no. 4, pp. 632-640. doi: 10.1007/BF00851575.
21. Radchenko V. P., Kichaev E. K., Simonov A. V. Energetic variant of the model of rheological deformation and destruction of metals 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.
22. Khokhlov A. V. Nonlinear Maxwell-type model of viscoelastoplasticity: the properties of creep curves at step loadings and conditions for the accumulation of plastic deformation, Mashinostroenie i inzhenernoe obrazovanie, 2016, no. 3, pp. 35-48 (In Russian).
23. Khokhlov A. V. Nonlinear Maxwell-type models of viscoelasticity. Features of their behavior, rate sensitivity, and the possibility of using to describe the creep and superplastic materials, Research report no. 5193. Moscow, Lomonosov Moscow State University, Institute of Mechanics, 2013, 108 pp. (In Russian)
24. Gorodtsov V. A., Leonov A. I. On the kinematics, nonequilibrium thermodynamics, and rheological relationships in the nonlinear theory of viscoelasticity, J. Appl. Math. Mech., 1968, vol.32, no. 1, pp. 62-84. doi: 10.1016/0021-8928(68)90148-2.
25. Leonov A. I., Lipkina E. Ch., Paskhin E. D., Prokunin A. N. Theoretical and experimental investigations of shearing in elastic polymer liquids, Rheol. Acta, 1976, vol. 15, no. 7/8, pp. 411-426. doi: 10.1007/BF01574496.
26. Pal'mov V. A. Rheological models in the nonlinear mechanics of deformable bodies, Uspekhi mekhaniki, 1980, vol. 3, no. 3, pp. 75-115 (In Russian).
27. Prokunin A. N. On the non-linear Maxwell-type defining equations for describing the motions of polymer liquids, J. Appl. Math. Mech., 1984, vol.48, no. 6, pp. 699-706. doi:10.1016/0021-8928(84)90037-6.
28. Leonov A. I., Prokunin A. N. Non-linear Phenomena in Flows of Viscoelastic Polymer Fluids. London, Chapman and Hall, 1994, xvii+475 pp. doi: 10.1007/978-94-011-1258-1.
29. Leonov A. I. Constitutive equations for viscoelastic liquids: Formulation, analysis and comparison with data, Rheology Series, 1999, vol. 8, pp. 519-575. doi: 10.1016/S0169-3107(99) 80040-9.
30. Cao Y. Determination of the creep exponent of a power-law creep solid using indentation tests, Mech. Time-Depend. Mater, 2007, vol.11, no. 2, pp. 159-172. doi: 10.1007/ s11043-007-9033-6.
31. Naumenko K., Altenbach H., Gorash Y. Creep Analysis with a Stress Range Dependent Constitutive Model, Arch. Appl. Mech., 2009, vol.79, no. 6, pp. 619-630. doi: 10.1007/ s00419-008-0287-5.
32. Radchenko V. P., Shapievskii D. V. Analysis of a nonlinear generalized Maxwell model, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2005, no. 38, pp. 55-64 (In Russian). doi: 10.14498/vsgtu372.
Received 15/VII/2016;
received in revised form 04/IX/2016;
accepted 09/IX/2016.
