CC BY
35
УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА.
_ СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ
2017, Т. 159, кн. 4 С. 509-517
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)
УДК 539.3
ИССЛЕДОВАНИЕ КОНЕЧНЫХ УПРУГОПЛАСТИЧЕСКИХ ДЕФОРМАЦИЙ: АЛГОРИТМ РЕШЕНИЯ, ЧИСЛЕННЫЕ ПРИМЕРЫ
Л.У. Султанов
Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия
Аннотация
Работа посвящена разработке методики расчета упругопластических трехмерных тел с учетом конечных деформаций. Кинематика упругопластических деформаций основана на мультипликативном разложении полного градиента деформации на упругую и неупругую составляющие. Напряженное состояние характеризуется тензором напряжений Коши. Физические соотношения получены на основе уравнения второго закона термодинамики с введением функции свободной энергии. Функция свободной энергии записана в виде зависимости от инвариантов левого тензора упругой деформации Коши-Грина. Рассмотрена упругопластическая модель с изотропным упрочнением. На основе аналога ассоциированного закона пластического течения и критерия пластичности разработан метод проецирования напряжений на поверхность текучести с итерационным уточнением текущего напряженно-деформированного состояния. Итерационная процедура основана на введение в разрешающее уравнение мощности дополнительных напряжений. Построены определяющие соотношения для скоростей и приращений истинных напряжений Коши. В рамках метода последовательных нагружений получено вариационное уравнение, основанное на принципе возможных мощностей. Пространственная дискретизация основана на методе конечных элементов, использован восьмиузловой конечный элемент. Представлено решение задачи о растяжении стержня круглого поперечного сечения и дано сравнение с результатами других авторов.
Ключевые слова: нелинейная упругость, конечные деформации, пластичность
Введение
Настоящая работа является продолжением работы [1], в которых приведены основные кинематические соотношения и определяющие соотношения. Используется мультипликативное разложение градиента деформаций на упругую и неупругую составляющие. Для разделения упругих и неупругих деформаций применяется метод проецирования напряжений на поверхность текучести, получен аналог ассоциированного закона течения. Решение задачи ведется в рамках метода последовательных нагружений на базе уравнения принципа виртуальных мощностей. Численная дискретизация основана на методе конечных элементов.
1. Кинематика
Кинематику конечных упругопластических деформаций будем описывать с помощью градиента деформации Е, который представляется в виде мультипликативного разложения [1—8]:
Е = Ее ■ Ер,
где Ее - градиент упругих деформаций, Ер - градиент неупругих деформаций. Напряженное состояние описывается тензором напряжений Коши а [1, 9]:
2 В д'Ф 2ро _ д'ф
^ = 2рВе ■ Ж = ~ГВе ■ (1)
где Ве = Ее • Ее - левый тензор Коши-Грина упругих деформаций, J = ро/р -относительное изменение объема, ро - плотность среды в начальном состоянии, р -плотность среды в актуальном состоянии, ф - функция свободной энергии, которая определяется как скалярная функция от тензора меры упругих деформации Ве. В качестве условия пластического течения примем
Фр (а, х) = 0, (2)
где Фр - функция текучести, х - параметр упрочнения. Воспользуемся выра-
V
жением для обобщенной скорости Ве , которое получено из условия минимума обобщенного функционала [1]:
Ве = -2А• Ве. (3)
д а
Линеаризируя определяющие соотношение (1) и используя (3), получим тензор скорости напряжения Коши [1]:
а = Ле • •а + Ь • а + а • Ьт - 1ика-
- -р { дфР • Ве • £ + Ве • дВфтУ^Ве} =
где
Ле • а + Ь • а + а • Ьт - Л \ 2^ • а + Ле • • ^ \ = ае + &р, (4)
,дФр . дФр
ы
Ле • а + Ь • а + а • Ьт - Тыка, ар = -Л1 • а + Л
да да
= 4ро В д 2Ф В
Ле = J Ве ^ дВедВе ^ Ве'
2. Алгоритм интегрирования определяющих соотношений
Будем использовать метод проецирования напряжений на поверхность текучести [6, 10-14]. Для этого запишем скорость напряжений Коши (4) в к-м состоянии
ка = кЛе • ка + кЬ • ка + ка • кЬт - 11клка-
- кл {2+ -Л. ••д+а}. (5)
По известным параметрам к-го состояния, используя (5), определим (к + 1)-е по следующей формуле:
к+1 _ к _ , к
а = ка + к а Дt
ка+
кЛе • • ка + кЬ • ка + ка • кЬт - 1ыка - кЛ ( 2^ • ка + кЛе • • ^
а
дка е дк
Д,
к-+ [кЛе • • ка + кЬ • ка + ка • кЬт - 1ыка] Д-Л ^2дф • ка + кЛе •
а Д-(2^ а + *Ле • •да ) Д,
^ -Л (2^'ка + кЛе • • дфА Д, (6)
Введя в рассмотрении тензор пробных напряжений
к+1а = ка + [к+1Ле ■ ка + кЬ ■ ка + к а ■ кЬт - 1Хк/а] АЬ, из (6) получим уравнение для определения истинных напряжений
к+1 а = ^ - кл (2^ ■ а + ■ ■ АЬ. (7)
Алгоритм расчета состоит в нахождении тензора истинных напряжений и скорости пластических деформаций из критерия пластичности (2) и соотношения (7).
3. Интегрирование уравнения движения
Представим процесс деформирования в виде последовательности напряженно-деформированных состояний (НДС) исследуемого тела, реализующихся в определенные моменты времени (временные слои). Подобная стратегия решения нелинейных задач является в настоящее время доминирующей и с успехом применяется в задачах как статики (метод последовательных нагружений), так и динамики (метод пошагового интегрирования) [13-15]. В соответствии с этой методикой считаем,
к
что в некоторый момент времени кЬ известны все параметры процесса, включая конфигурацию, напряженное состояние, величины упругих и неупругих деформаций и т. д. Задача состоит в определении НДС в момент времени к+1Ь = кЬ + АЬ.
В качестве базового примем уравнение виртуальных мощностей в актуальном состоянии [9, 16, 17]. Запишем его в виде операторного уравнения О = 0. На к-м временном слое должно выполняться уравнение к О = 0. Аналогичное уравнение на следующем временном слое можно представить в виде к+1О = кО + к (ОАЬ = 0. В нашем случае уравнение виртуальных мощностей в актуальном состоянии (при Ь = к+1Ь) принимает вид
У ка ■ ■ёкавУ - I кгп5иё,Б - I кfёивУ I +
\Ук Б? Ук I
+ в < Iка ■ ёкавУ - I к^ ■ ёивБ - I кf ■ ёивУ I АЬ = 0. (8)
(Ук
Ук
Линеаризуя (8) и принимая во внимание соотношение (5), получим разрешающее уравнение
Ук
а ■ кЛе ■ ■¿а +2ка ■■ [¿Ьт-кЬ+кЬт ■ ёЬ] - [кV ■ ки] кf ■ ёи-
- л
2 дфр а + Л • дфр
2"я- ■ а + Ле • -
да да
вУ - {кК ■ кЬт - [кV ■ ки] kt4 ёивБ
Jk Аtn ■ ёивБ + ! к Аf ■ ёивУ-
Ук
1
АЬ
кк
• -ёкавУ - ки ■ ёивЬ1- к f ёивУ
а ■ -и — J
[Ук Б? Ук
Так как решаются квазистатические задачи, то можно перейти от скоростей к приращениям, например, полагая и = Ди/Д,.
Решая уравнение (9), получим вектор перемещений и, который определяет конфигурацию на следующем шаге нагружения
к+1х = кх + иД,
В задачах статики приращение времени является параметром нагружения, который обычно полагают равным единице. При возникновении пластических деформаций используется метод проецировании напряжений с итерационным уточнением, который заключается в следующем. Путем решения упругой задачи (9) находим приращение поля перемещений, по которым определяем поле пробных напряжений на следующем шаге нагружения
а = к а + [к+1Ле • • ка + кЬ • ка + ка • кЬт - /1Мка] .
Затем применяем метод проецирования напряжений на поверхность текучести, который заключается в решении системы нелинейных уравнений
к + 1,
Д 2
дФр
'д^а
Ф,
(к+1
сг
•а + к+1Ле--
х (ДЛ)) =0.
д Фр
дк+1 а
0,
(10) (11)
Решив систему нелинейных уравнений (10), (11) методом Ньютона, получим тензор истинных напряжений к+1а и приращение скорости пластических деформаций Д . В силу того, что полученное истинное напряженное состояние не удовлетворяет уравнению равновесия (9), строится итерационная процедура уточнения полученного НДС, которая основана на введении в разрешающую систему уравнений (9) мощности дополнительных напряжений на возможных деформациях скорости
Ук
га • • кЛе • • за + 2к а •• [ЗЬт • ктЬ + ктЬт • ЗЬ] - [кV • ки] kf • 6и-
-Дк Л
т
2 -дФ^ . к а+к Л .. дФр 2 дт а та+ Ле дт а
} dV - I {кЩ • тЬт - ^ • кти] кtn} *ис£!
= J Дк^ • SvdS + у Дкf • ёиЗУ-
як Ук
I ка • • ЗкadV - I к^^ SvdS - I к f• SvdV I + I т,аай • • ЗкadV,
Ук
Ук
Ук
где т+1 ааа=т+1а - а •
к
4. Численный пример
Рассмотрим задачу растяжения круглого стержня со следующими параметрами: радиус поперечного сечения К = 6.413 мм, длина стержня Ь = 53.334 мм. Отметим, что для конкретизации места образования шейки в центре стержня задается снижение радиуса на 1.8% [2-4]. Функция свободной энергии задается в следующем виде:
роФ = Л + (11В - 3)2 + - (11в - 3) - 2 (12В - 3),
Рис. 1. График зависимости силы от перемещения. Сплошная линия - зависимость, полученная с помощью разработанной нами методики, □ - [2], А - [3]
Рис. 2. Интенсивность пластических деформаций
где А, ¡л - параметры Ляме. В качестве критерия упругого деформирования примем условие Губера - Мизеса, которое для изотропной среды допускает обобщение
в виде
Ф = ^ - ат (х) < 0-
Здесь ^ = V 2 ^ • а'- и—ть напряжений, (х) - функция упрочнения, х - параметр упрочнения, а' - девиатор тензора напряжений.
Функция нелинейного изотропного упрочнения имеет вид [2-4]
ат(х) = ат + Нх + (аж - ат)(1 - е-&х).
Параметры материала следующие: Е = 206900 МПа, V = 0.29, аж = 715 МПа, ат = 450 МПа, Н = 0.129, 5 = 16.93. Численная реализация основана на методе конечных элементов. Используется трехмерный 8-узловой конечный элемент. На рис. 1 показан график зависимости силы, возникающей на торце стержня, от перемещения. Сплошной линией показана зависимость, полученная с помощью вышеописанной методики, □ - [2], ▲ - [3]. Видно, что наши результаты хорошо согласуются с результатами других авторов На рис. 2 изображено деформированное состояние стержня с полем интенсивности пластических деформаций. Интенсивность пластических деформаций принимает максимальное значение в области образования шейки.
Заключение
В работе представлена методика решения задач деформирования с учетом конечных упругопластических деформаций. Используется мультипликативное разло-
жение градиента деформации. Определяющее соотношение записано в актуальном состоянии в виде зависимости тензора напряжения Коши от левого тензора деформации Коши-Грина. Для учета пластических деформаций применяется метод проецирования напряжений на поверхность текучести с итерационным уточнением текущего НДС. В рамках метода последовательных нагружений получено разрешающее уравнение. Численная реализация основана на методе конечных элементов. В качестве демонстрации работоспособности разработанной методики решена задача о растяжении стержня с образованием шейки и проведено сравнение результатов решения с решением других авторов.
Благодарности. Работа выполнена при финансовой поддержке РНФ в рамках научного проекта 16-11-10299.
Литература
1. Султанов Л.У. Исследование конечных упругопластических деформаций. Кинематика среды и определяющие соотношения // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2015. - Т. 157, кн. 4. - С. 158-165.
2. Eidel B, Gruttmann F. Elastoplastic orthotropy at finite strains: multiplicative formulation and numerical implementation // Comput. Mater. Sci. - 2003. - V. 28, No 3-4. - P. 732-742. - doi: 10.1016/j.commatsci.2003.08.027.
3. Schroder J., Gruttmann F., Loblein J. A simple orthotropic finite elasto-plasticity model based on generalized stress-strain measures // Comput. Mech. - 2002. - V. 30, No 1. -P. 48-64. - doi: 10.1007/s00466-002-0366-3.
4. Simo J.S. A framework for finite strain elastoplasticity based on maximum plastic dissipation and the multiplicative decomposition: Part I. Continuum formulation // Comput. Meth. Appl. Mech. Eng. - 1988. - V. 66, No 2 - P. 199-219. - doi: 10.1016/0045-7825(88)90076-X.
5. Miehe C. A theory of large-strain isotropic thermoplasticity based on metric transformation tensors // Arch. Appl. Mech. - 1995. - V. 66, No 1-2. - P. 45-64. -doi: 10.1007/BF00786688.
6. Basar Y., Itskov M. Constitutive model and finite element formulation for large strain elasto-plastic analysis of shell // Comput. Mech. - 1999. - V. 23, No 5-6. - P. 466-481. -doi: 10.1007/s004660050426.
7. Meyers A., Schievbe P., Bruhns O.T. Some comments on objective rates of symmetric Eulerian tensors with application to Eulerian strain rates // Acta Mech. - 2000. - V. 139, No 1-4. - P. 91-103. - doi: 10.1007/BF01170184.
8. Xiao H., Bruhns O.T., Meyers A. A consistent finite elastoplasticity theory combining additive and multiplicative decomposition of the stretching and deformation gradient // Int. J. Plasticity. - 2000. - V. 16, No 2. - P. 143-177. - doi: 10.1016/S0749-6419(99)00045-5.
9. Bonet J., Wood R.D. Nonlinear continuum mechanics for finite element analysis. -Cambridge: Cambridge Univ. Press, 1997. - 283 p.
10. Rouainia M., Wood D.M. Computational aspects in finite strain plasticity analysis of geotechnical materials // Mech. Res. Commun. - 2006. - V. 33, No 2. - P. 123-133. - doi: 10.1016/j.mechrescom.2005.06.014.
11. Simo J.S., Ortiz M. A unified approach to finite deformation elastoplastic analysis lased on the use of hyperelastic constitutive equations // Comput. Methods. Appl. Mech. Engng. - 1985. - V. 49, No 2. - P. 221-245. - doi: 10.1016/0045-7825(85)90061-1.
12. Eterovic A.L., Bathe K.-J. A hyperelastic-based large strain elasto-plastic constitutive formulation with combined isotropic-kinematic hardening using the logarithmic stress and
strain measures // Int. J. Numer. Meth. Engng. - 1990. - V. 30, No 6. - P. 1099-1114. -doi: 10.1002/nme.1620300602.
13. Давыдов Р.Л., Султанов Л.У. Численный алгоритм решения задачи о больших упру-гопластических деформациях МКЭ // Вестн. ПНИПУ. Механика. - 2013. - № 1. -C. 81-93.
14. Davydov R.L., Sultanov L.U. Numerical algorithm for investigating large elasto-plastic deformations // J. Eng. Phys. Thermophys. - 2015. - V. 88, No 5. - P. 1280-1288. - doi: 10.1007/s10891-015-1310-7.
15. Golovanov A.I., Sultanov L.U. Numerical investigation of large elastoplastic strains of three-dimensional bodies // Int. Appl. Mech. - 2005 - V. 41, No 6. - P. 614-620. - doi: 10.1007/s10778-005-0129-x.
16. Abdrakhmanova A.I., Sultanov L.U. Numerical modelling of deformation of hyperelastic incompressible solids // Materials Phys. Mech. - 2016. - V. 26, No 1. - P. 30-32.
17. Голованов А.И., Коноплев Ю.Г., Султанов Л.У. Численное исследование конечных деформаций гиперупругих тел. IV. Конечноэлементная реализация. Примеры решения задач // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2010. - Т. 152, кн. 3. -С. 115-126.
Поступила в редакцию 29.08.17
Султанов Ленар Усманович, кандидат физико-математических наук, доцент кафедры
теоретической механики
Казанский (Приволжский) федеральный университет
ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: [email protected]
ISSN 2541-7746 (Print) ISSN 2500-2198 (Online) UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI (Proceedings of Kazan University. Physics and Mathematics Series)
2017, vol. 159, no. 4, pp. 509-517
Analysis of Large Elastic-Plastic Deformations: Integration Algorithm and Numerical Examples
L.U. Sultanov
Kazan Federal University, Kazan, 420008 Russia E-mail: [email protected]
Received August 29, 2017 Abstract
An algorithm for investigation of elastic-plastic solids with regard to finite deformations has been developed. The kinematics of elastic-plastic deformations is based on multiplicative decomposition of the deformation gradient into elastic and inelastic parts. The stress state is determined by the Cauchy stress tensor. The constitutive equations have been obtained from
the second law of thermodynamics with the introduction of the elastic free energy function. The elastic free energy function is formulated in an invariant form of the left Cauchy-Green tensor. The von Mises yield criterion with isotropic hardening has been used. The radial return method with an iterative refinement of the current state of deformation has been applied for dividing the elastic and plastic deformations. The principle of virtual work in terms of the virtual velocity has been used. The numerical implementation is based on the finite element method. The solution of the necking of a circular bar has been presented.
Keywords: nonlinear elasticity, finite deformations, plasticity
Acknowledgments. The study was supported by the Russian Science Foundation (project no. 16-11-10299).
Figure Captions
Fig. 1. Force-displacement graph. Solid line - a dependence obtained with the help of the developed algorithm, □ - [2], A - [3].
Fig. 2. Plastic deformation intensity.
References
1. Sultanov L.U. Analysis of finite elastoplastic deformations. Kinematics and constitutive equations. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2015, vol. 157, no. 4, pp. 158-165. (In Russian)
2. Eidel B., Gruttmann F. Elastoplastic orthotropy at finite strains: Multiplicative formulation and numerical implementation. Comput. Mater. Sci., 2003, vol. 28, nos. 3-4, pp. 732-742. doi: 10.1016/j.commatsci.2003.08.027.
3. Schroder J., Gruttmann F., Loblein J. A simple orthotropic finite elasto-plasticity model based on generalized stress-strain measures. Comput. Mech., 2002, vol. 30, no. 1, pp. 4864. doi: 10.1007/s00466-002-0366-3.
4. Simo J.S. A framework for finite strain elastoplasticity based on maximum plastic dissipation and the multiplicative decomposition: Part I. Continuum formulation. Com-put. Methods Appl. Mech. Eng., 1988, vol. 66, no. 2, pp. 199-219. doi: 10.1016/0045-7825(88)90076-X.
5. Miehe C. A theory of large-strain isotropic thermoplasticity based on metric transformation tensors. Arch. Appl. Mech., 1995, vol. 66, nos. 1-2, pp. 45-64. doi: 10.1007/BF00786688.
6. Basar Y., Itskov M. Constitutive model and finite element formulation for large strain elasto-plastic analysis of shell. Comput. Mech., 1999, vol. 23, nos. 5-6, pp. 466-481. doi: 10.1007/s004660050426.
7. Meyers A., Schievbe P., Bruhns O.T. Some comments on objective rates of symmetric Eulerian tensors with application to Eulerian strain rates. Acta Mech., 2000, vol. 139, nos. 1-4, pp. 91-103. doi: 10.1007/BF01170184.
8. Xiao H., Bruhns O.T., Meyers A. A consistent finite elastoplasticity theory combining additive and multiplicative decomposition of the stretching and deformation gradient. Int. J. Plast., 2000, vol. 16, no. 2, pp. 143-177. doi: 10.1016/S0749-6419(99)00045-5.
9. Bonet J., Wood R.D. Nonlinear Continuum Mechanics for Finite Element Analysis. Cambridge, Cambridge Univ. Press, 1997. 283 p.
10. Rouainia M., Wood D.M. Computational aspects in finite strain plasticity analysis of geotechnical materials. Mech. Res. Commun., 2006, vol. 33, no. 2, pp. 123-133. doi: 10.1016/j.mechrescom.2005.06.014.
11. Simo J.S., Ortiz M. A unified approach to finite deformation elastoplastic analysis lased on the use of hyperelastic constitutive equations. Comput. Methods Appl. Mech. Eng., 1985, vol. 49, no. 2, pp. 221-245. doi: 10.1016/0045-7825(85)90061-1.
12. Eterovic A.L., Bathe K.-J. A hyperelastic-based large strain elasto-plastic constitutive formulation with combined isotropic-kinematic hardening using the logarithmic stress and strain measures. Int. J. Numer. Meth. Eng., 1990, vol. 30, no. 6, pp. 1099-1114. doi: 10.1002/nme.1620300602.
13. Davydov R.L., Sultanov L.U. Numerical algorithm of solving the problem of large elastic-plastic deformations by FEM. Vestn. PNIPU. Mech., 2013, no. 1, pp. 81-93. (In Russian)
14. Davydov R.L., Sultanov L.U. Numerical algorithm for investigating large elasto-plastic deformations. J. Eng. Phys. Thermophys., 2015, vol. 88, no. 5, pp. 1280-1288. doi: 10.1007/s10891-015-1310-7.
15. Golovanov A.I., Sultanov L.U. Numerical investigation of large elastoplastic strains of three-dimensional bodies. Int. Appl. Mech., 2005, vol. 41, no. 6, pp. 614-620. doi: 10.1007/s10778-005-0129-x.
16. Abdrakhmanova A.I., Sultanov L.U. Numerical modelling of deformation of hyperelastic incompressible solids. Mater. Phys. Mech., 2016, vol. 26, no. 1, pp. 30-32.
17. Golovanov A.I., Konoplev Yu.G., Sultanov L.U. Numerical investigation of finite deformations of hyperelastic bodies. IV. Finite-element implementation. Examples of the solution of problems. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, 2010, vol. 152, no. 3, pp. 115-126. (In Russian)
I Для цитирования: Султанов Л.У. Исследование конечных упругопластических деформаций: алгоритм решения, численные примеры // Учен. зап. Казан. ун-та. \ Сер. Физ.-матем. науки. - 2017. - Т. 159, кн. 4. - С. 509-517.
/ For citation: Sultanov L.U. Analysis of large elastic-plastic deformations: Integration ( algorithm and numerical examples. Uchenye Zapiski Kazanskogo Universiteta. Seriya \ Fiziko-Matematicheskie Nauki, 2017, vol. 159, no. 4, pp. 509-517. (In Russian)
