CC BY
52
УДК 539.376
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВЯЗКОУПРУГОГО РАЗУПРОЧНЯЮЩЕГОСЯ МАТЕРИАЛА С ЭКСПОНЕНЦИАЛЬНЫМ ЯДРОМ ПОЛЗУЧЕСТИ
С. В. Горбунов
Самарский государственный технический университет,
443100, Россия, Самара, ул. Молодогвардейская, 244.
E-mail: 0gorbunov0@gmail. com
Предложен вариант математической модели одноосного деформирования вязкоупругого материала с ядром ползучести экспоненциального вида. Исследована устойчивость решений модели по Ляпунову при постоянных напряжениях. Установлена область устойчивости решений системы дифференциальных уравнений математической модели, соответствующая асимптотически ограниченной ползучести материала. Область неустойчивости решений соответствует появлению третьей стадии ползучести. Установлена связь между устойчивостью решений по Ляпунову и устойчивостью вычислительного алгоритма при численном решении системы уравнений. В качестве иллюстрации приведено исследование модельной задачи.
Ключевые слова: вязкоупругий материал, устойчивость решений по Ляпунову, ядро ползучести экспоненциального вида, область устойчивости решений, третья стадия ползучести, устойчивость вычислительного алгоритма.
1. Проблема построения моделей неупругого реологического деформирования имеет давние корни [1-11] и в настоящее время далека от завершения. Разработанные модели ориентированы в основном на конструкционные материалы. Однако имеется класс вязкоупругих материалов (биоматериалы, полимеры, органопластики), деформирование которых в рамки имеющихся моделей не укладывается. Так, для указанных материалов существует напряжение <т = сгПДс, называемое пределом длительного сопротивления, такое, что если а < <тПДс) то деформация ползучести р = p(t) при постоянном напряжении er = do = const является асимптотически ограниченной, т. е. lim p(t) = р°° (<то), где р°° (<то) — предельное значение деформации пол-
£—>■ со
зучести, причём разрушения образца не происходит, а при полной разгрузке образца при t = t* (<то = 0, t > t*) вязкоупругая деформация полностью обратима. Если же а > <тПДс, то вязкоупругая деформация является неограниченной, появляется третья стадия ползучести, и процесс деформирования заканчивается разрушением образца. Схематически реологическое деформирование одноосных образцов показано на рисунке. Цель настоящей работы — разработка модели реологического деформирования этого класса материалов, при этом основной задачей является установление связи разупрочнения материала с процессом неустойчивого поведения решения предложенной ниже системы дифференциальных уравнений.
2. Для описания асимптотически ограниченной ползучести вязкоупругих материалов можно использовать экспоненциальные ядра ползучести [3,12]. Поэтому в настоящей работе рассматривается частный случай теории пол-
Сергей Владимирович Горбунов, аспирант, каф. прикладной математики и информатики.
Схематические кривые ползучести при постоянных напряжениях:
О 1 < (72 "С Опдс ^ <73 < О4
зучести [3,12] применительно к исследуемому классу вязкоупругих материалов, но при этом вводятся дополнительные кинетические уравнения, описывающие процесс накопления повреждённости. Предлагаемый вариант модели имеет следующий вид:
p(t) =pi(t) +p2(t), (1) Pi(t) = Ai (2)
P2{t) = A2 (fl2(^) -P2(i)), (3)
a = a0(l +uj), (4)
oj = aap, (5)
с начальными условиями p(0) = 0, w(0) = 0. В системе (l)-(5) p(t) — деформация ползучести; pi(t), P2(t) — составляющие деформации ползучести; do — номинальное напряжение; <т — истинное напряжение; со — параметр повреждённости; Ai, А2, ai, d2i rn, a, a* — константы модели.
Если повреждённость незначительна и величиной со можно пренебречь {со ¡=3 0), то <т = do и решение (1)-(5) при <то = const задаётся функцией
¡>(0 = £а* (l-e->“)(^)m, (6)
k= 1
т. е. имеем асимптотически ограниченную (при t —> 00) ползучесть с двумя экспоненциальными слагаемыми. Покажем, что при существенных значениях параметра повреждённости со система (1)—(5) описывает третью стадию ползучести, которая связана с неустойчивостью решения данной системы при do = const.
Найдём первый интеграл системы, разрешив совместно (4) и (5):
= aaodp] In (1 + со) = aaop + In С.
1 + со
Используя начальные условия р(0) = 0, w(0) = 0, находим С = 1. Таким образом, после элементарных преобразований первый интеграл системы
принимает вид
1+и = еааор. (7)
Запишем уравнения (2), (3) с учётом (7):
Ш = Ьк (ак^уеа°°рт-Рк(1)) , к = 1,2. (*
к т
-а"'
Введём следующие обозначения:
т
-4 ) =Ъ, аа0т = г]. (9)
С учётом (9) уравнения (8) принимают вид:
Рк&) = Хк(акЬеГ]Р-рк^)), к = 1,2. (10)
Таким образом, задача свелась к исследованию на устойчивость системы (10) с учётом (1) при различных постоянных значениях напряжения <то. Для этой цели применим метод Ляпунова исследования устойчивости по первому приближению. Введём возмущённые переменные
Ук=Рк~Р°к, к = 1,2, (11)
где р^—установившиеся движения. Согласно (1), р° = р\ + р2—также установившееся движение. Найдём его из системы (10), положив возмущения равными нулю (у 1 = 0, 2/2 = 0):
р\ = афёГ1Р\ р% = а2Ъё^\ (12)
отсюда
р° = (а1 + а2)ЬеГ1Р°. (13)
Уравнение (13) относительно р° при заданном значении <то можно разрешить численными методами.
Составим систему уравнений возмущённого движения, учитывая, что согласно (11), рк = ук'-
Г у\ = А1а16е’7(р°+?/1+да) - А! (р\ + у\) ,
I У2 = Х2а2Ье^р0+У^ - Х2 (р°2 + у2). [ 1
Разложим экспоненту в правых частях уравнений системы (14) как функцию двух переменных у\ и у2 в ряд Тейлора в окрестности точки (0; 0):
У1 = АхахЪе^0 + ХхсцЩё™0 {у\ + у2) - Xхр^ - Хт +о(ух + у2), ^
у2 = Х2а2Ье11Р° + А(уг + у2) - X2р\ - Х2у2 + о(ух + у2).
Исключим из уравнений (15) слагаемые более высокого порядка малости, чем (у\ +у2), тем самым получим систему уравнений первого приближения. С учётом равенства (12) окончательно получим
Г У1 = А1 (щ>\ - 1)2/1 + Х1цр\у2, \ у2 = ^2Г]Р2У1 + А2 (г]р°2 — 1) г/2-
Исследуем на устойчивость положение равновесия у\ = 0, у2 = 0 системы уравнений первого приближения (16). Представим систему уравнений в матричном виде:
Ч = Лу,
(М
где у = [уі у2\т и
А =
Лі (г]р\ - і) ХіщРі Мт>% Л2 (т>% -1)
Найдём корни характеристического уравнения ёе! (А — гЕ) = 0, где Е — единичная матрица:
г2 + г (А1 (1 - г]Р\) + А2 (1 - Г]Р°2)) + А1А2 (1-Г] (р5 +р°2)) = 0.
В данном случае воспользуемся необходимым и достаточным условием асимптотической устойчивости положения равновесия — условием Стодолы, которое здесь принимает вид
Г Аі (1 — Г]р°і) + А2 (1 - т%) > 0, , ■
\ АіА2 (1 - г] (рЧ + р%)) >0.
Таким образом, при выполнении условия (17), отвечающего заданному значению сто, положение равновесия (13) системы (1)—(5) будет асимптотически устойчивым.
Конкретизируем условие (17) для рассматриваемого частного случая, в котором (из физических предпосылок)
А1 >0, А2 > 0, а\ >0, я>2 > 0, р\ > 0, р2 > 0. (18)
Согласно (18) А1А2 > 0, и второе условие в (17) упрощается:
А1 + А 2-г] (А1 р\ + А2р°) > 0,
V ^ ~о | о >
Р1+Р2
отсюда
А1 + А2
Г] <
АіРі + Х2Р2 ’
Л<~о-
ри
(19)
Рассмотрим подробнее первое неравенство в (19), преобразовав знаменатель с учётом (12):
А1 + А2 , Л
Л (А^! + А2а2) ЬёЧР°'
Записывая (13) в виде 6е’7Р° = р°/(а\ + а2) и подставляя это выражение в (20), получим
1 (^ | А1Я.2 + А20,1 ^ р°
Л< ч°\1+Хкн + Х2а2 У ' (21)
Из условий (18) следует, что
\\ti2 + Лгй! р Л101 + \2Ci2
поэтому правая часть неравенства (21) будет больше правой части второго неравенства (19). Значит, с учётом обозначения (9), окончательно получаем условие асимптотической устойчивости в виде
где р° находится из (13).
3. В качестве модельного примера использовался материал со следующими характеристиками: а\ = 2,3 • 10_3; = 0,47 • 10_3; Ai = 0,015; Л2 = 0,056;
т = 0,51; а = 0,22; а* = 150 МПа. В результате применения критерия (22) получено значение <тПДс = 587 МПа, которое отделяет область устойчивого деформирования (асимптотически ограниченная ползучесть при <то = const) для gq ^ <7ПДС от области неустойчивого деформирования (появление третьей стадии ползучести) для Gq > (Тпдс-
Следует отметить, что систему (1)—(5) можно разрешить численно, например методом Эйлера («шагами» по времени). Установлено, что значение g = (Тпдс является граничным значением для номинального напряжения <то, при котором численный метод либо сходится (<то ^ <тПдо асимптотически ограниченная ползучесть), либо расходится (<то > <тПДс, неограниченная ускоренная ползучесть), при этом получено, что значение (Тпдс, расчитанное на основе численного решения задачи (1)—(5) методом Эйлера, практически совпадает со значением <тПДс по критерию (22).
Таким образом, установлена прямая связь между устойчивостью решений системы (1)—(5) и устойчивостью численных методов решения этой системы, а появление стадии ускоренной ползучести связано с нарушением условий устойчивого деформирования (22) и потерей устойчивости (расходимостью) численного алгоритма. Аналогичные результаты для пластически разупроч-няющихся материалов получены В. В. Стружановым [10,11,13]. В заключение следует отметить, что никаких принципиальных трудностей с исследованием системы уравнений типа (1)—(5), но с большим количеством экспоненциальных слагаемых в (6) и, соответственно, с большим количеством слагаемых деформации Pi(t) в (1)—(5), не возникает.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Работное Ю. Н. Ползучесть элементов конструкций. М.: Наука, 1966. 752 с. [Rabotnov Yu. N. Creep of Structural Elements. Moscow: Nauka, 1966. 752 pp.]
2. Радченко В. П., Павлова Г. A., Горбунов C.B. Устойчивость по Ляпунову решений эндохронной теории пластичности без поверхности текучести в условиях плоского напряженного состояния// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008. №2(17). С. 143-151. [Radchenko V.P., Pavlova G. A., Gorbunov S.V. Stability by Lyapunov of solutions in endochronic plasticity theory without fluidity surface in flat tension conditions// Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2008. no. 2(17). Pp. 143-151].
3. Радченко В. П., Еремин Ю. А. Реологическое деформирование материалов и элементов конструкций. М.: Машиностроение-1, 2004. 264 с. [Radchenko V.P., Eremin Yu. А. Rheological Deformation and Failure of Materials and Structural Elements. Moscow: Mashinostroenie-1, 2004. 264 pp.]
4. Ибрагимов В. А., Клюшников В. Д. Некоторые задачи для сред с падающей диаграммой// Изв. АН СССР. МТТ, 1971. №4. С. 116-121. [Ibragimov V.A., Klyushnikov V.D. Some problems for media with an incident diagram// Izv. AN SSSR. MTT, 1971. no. 4. Pp. 116-121].
5. Кадашевич Ю. И. Теория пластичности и ползучести, учитывающая микроразрушения // Докл. АН СССР, 1982. Т. 266, №6. С. 1341-1344. [Kadashevich Yu. I. The theory of plasticity and creep taking into account the microfracture // Dokl. AN SSSR, 1982. Vol. 266, no. 6. Pp. 1341-1344].
6. Лебедев А. А., Чаусов H. Г. Феноменологические основы оценки трещиностойкости материалов по параметрам спадающих участков диаграммы деформаций // Пробл. прочности, 1983. №2. С. 6-10; англ. пер.: Lebedev A. A., Chausov N. С. Phenomenological fundamentals of the evaluation of crack resistance of materials on the basis of parameters of falling portions of strain diagrams // Strength of Materials. Vol. 15, no. 2. Pp. 155-160.
7. Новожилов В. В., Кадашевич Ю. И. Микронапряжения в конструкционных материалах. JL: Машиностроение, 1990. 223 с. [Novozhilov V. V., Kadashevich Yu. I. Microstresses in Structural Materials. Leningrad: Mashinostroenie, 1990. 223 pp.]
8. Кадашевич Ю.И., Помыткин С. П. Эндохронная теория неупругости для разупроч-няющихся материалов с учётом больших деформаций / В сб.: Современные проблемы ресурса материалов и конструкций. М.: МАМИ, 2009. С. 158-165. [Kadashevich Yu.I., Pomytkin S. P. Endochronic theory of inelasticity for softening materials with regard to large deformations / In: Modem problems of resource materials and structures. Moscow: MAMI, 2009. Pp. 158-165].
9. Стружанов В. В. Свойства разупрочняющихся материалов и определяющие соотношения при одноосном напряженном состоянии // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2007. №2(15). С. 69-78. [Struzhanov V. V. The properties of softening materials and constitutive relations under uniaxial stress state // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2007. no. 2(15). Pp. 69-78].
10. Стружанов В. В., Миронов В. И. Деформационное разупрочнение материала в элементах конструкций. Екатеринбург: УрО РАН, 1995. 191 с. [Struzhanov V. V, Mironov V. I. Deformational Softening of Material in Structural Elements. Ekaterinburg: UrO RAN, 1995. 191 pp.]
11. Стружанов В. В., Вахарева Е. А. Итерационные процедуры расчёта параметров равновесия и устойчивость процесса чистого изгиба балок из пластических и хрупких разупрочняющихся материалов // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2010. №1(20). С. 84-95. [Struzhanov V. V., Bakhareva Е.А. Iterative procedures of equilibrium parameters estimation and process stability of pure bending of beams from soft and brittle weakening materials// Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2010. no. 1(20). Pp. 84-95].
12. Самарин Ю. П. Уравнения состояния материалов со сложными реологическими свойствами. Куйбышев: Куйбышевский госуниверситет, 1979. 84 с. [Samarin Yu. P. Equation of State for Materials with Complex Rheological Properties. Kuibyshev: Kuibyshev State Univ., 1979. 84 pp.]
13. Жижерин С. В., Стружанов В. В., Миронов В. И. Итерационные методы расчёта напряжений при чистом изгибе балки из повреждаемого материала // Вычисл. тех-нол., 2001. Т. 6, №5. С. 24-33. [Zgizherin S. V, Struzhanov V. V, Mironov V. I. Iterative computational methods of calculating stresses for a pure bending of a beam of a damageable material// Vychisl. Tekhnol, 2001. Vol. 6, no. 5. Pp. 24-33].
Поступила в редакцию 02/XI/2011; в окончательном варианте — 13/111/2012.
MSC: 74G55; 74Cxx
MATHEMATICAL MODEL OF VISCOELASTIC SOFTENING MATERIAL WITH EXPONENTIAL CREEP KERNEL
S. V. Gorbunov
Samara State Technical University,
244, Molodogvardeyskaya St., Samara, 443100, Russia.
E-mail: 0gorbunov0agmail.com
The variant of mathematical model of uniaxial strain for viscoelastic material with exponential creep kernel is proposed. Lyapunov stability of the solution of the model in case of permanent stress is investiga,ted,. The stability region of solutions of mathematical model’s differential equations, corresponding to asymptotically restricted creep of material, is established. Instability region of solutions is in accord with appearance of tertiary creep. Relation between stability of solutions by Lyapunov and, stability of iterative calculation for numerical solving the system of equations is established. As an illustration the investigation of model problem is quoted,.
Key words: viscoelastic material, Lyapunov stability of solutions, exponential creep kernel, stability region of solutions, tertiary creep, stability of numerical iterative calculation.
Original article submitted 02/XI/2011; revision submitted 13/111 /2012.
Sergey V. Gorbunov, Postgraduate Student, Dept, of Applied Mathematics & Computer Science.
