Численное исследование гиперупругих материалов
К.ф.-м.н, доцент Л.У. Султанов; аспирант Л.Р. Фахрутдинов,
ФГАОУ ВПО «Казанский (Приволжский) федеральный университет»
Аннотация. В настоящей работе рассматривается методика исследования гиперупругих тел, использующихся для определения больших деформаций нелинейных упругих материалов (различных полимеров, в том числе эластомеров, пены) с использованием меры деформации Фингера.
Дается кинематика движения среды, напряженное состояние описывается тензором истинных напряжений Коши - Эйлера. Большое внимание уделено алгоритму получения линеаризованных физических соотношений в напряжениях Коши - Эйлера. Рассмотрен пример получения линеаризованных физических соотношений в напряжениях Коши - Эйлера для материала, описываемого потенциалом Муни - Ривлина. Численная реализация основана на методе конечных элементов в рамках инкрементального метода.
Для проверки работоспособности методики решен ряд задач: задача о плоской деформации квадратной полосы и задача об упругом деформировании плиты под действием равномерного давления. Полученные в ходе решения результаты не противоречат данным, опубликованным ранее.
Ключевые слова: гиперупругие материалы; большие деформации; метод конечных элементов
Введение
В современном строительстве, наряду с традиционными стройматериалами, широко применяются новые, более технологичные типы материалов. К ним относят различные эластомеры и другие подобные полимеры, отличительной особенностью которых является то, что они допускают большие деформации, сохраняя при этом упругие свойства. Подобные материалы также называются гиперупругими. С точки зрения механики деформируемого твердого тела речь идет о нелинейно упругих телах, при деформировании которых необходимо учитывать геометрическую нелинейность в рамках больших деформаций. Подобного рода проблемы рассматриваются во многих статьях и обобщены в ряде монографий, среди которых можно выделить работы [1-9]. В них приводится обширная библиография по данному вопросу и изложены основные положения по построению определяющих соотношений, и в данной статье этот материал активно используется. Из числа журнальных публикаций можно отметить статьи [10-17], которые посвящены гиперупругим и термогиперупругим материалам, в которых также рассматриваются различные варианты построения определяющих соотношений для конкретных материалов и приводятся примеры решения задач. Методика, описанная в этой статье, в идейном плане близка к методам, использованным в работах [18, 19]. Соотношения приводятся в наиболее компактной форме прямого тензорного исчисления. Подробнее о методах и технологии работы с использованием такой формы записи можно ознакомиться в работах [1, 3, 6].
В настоящей работе рассмотрена методика исследования гиперупругих тел с использованием меры деформации Фингера. Приведена кинематика движения среды. Рассмотрен пример получения линеаризованных физических соотношений в напряжениях Коши - Эйлера для материала, описываемого потенциалом Муни - Ривлина. Численная реализация основана на методе конечных элементов в рамках инкрементального метода. Также решены некоторые задачи для проверки работоспособности методики.
Кинематика среды
Рисунок 1. Кинематика среды
В глобальной неподвижной системе координат с ортами в^, е\ , вз положение исследуемого деформируемого тела описывается с помощью:
г () = х' ())■ - радиус-вектор материальной частицы в недеформированном состоянии;
Я = у1 Щ - радиус-вектор материальной частицы в деформированном
состоянии;
и
и
(?) = Я (?)-г (?) = и' (?) -
)=Я)=у)е=и )в' -
вектор перемещения; вектор скорости.
Базовым тензором, играющим ключевую роль в кинематике конечных деформаций, является тензор градиента деформаций:
Т дУ
ЖМ) =№).
дх/
В качестве тензоров, описывающих деформацию и скорость деформации, используются:
Т
• левый тензор Коши - Грина (мера деформации Фингера): (В) = (Г)'(Г) ;
• тензор пространственного градиента скорости: (Ь) = (г )'(г 1);
тензор деформации скорости:
(<> )=2 [(* )+(а ))
)'(г-1) + (г-1 )Т ( )
Определяющие соотношения
Напряженное состояние описывается с помощью тензора истинных напряжений (Е) = (Ту (ëiëj), определенного в актуальном состоянии.
В качестве аргументов функции потенциальной энергии деформации примем компоненты тензора меры деформации Фингера, т. е. :
W = W (Bj),
где вводятся в рассмотрение меры деформации, которые не сопровождаются изменением объема:
) - 2/ (В) - 3 /з (В)
Здесь з - - относительное изменение объема.
Таким образом:
Ж - Ж0(3) + ж'(¡1В,12В). Тогда тензор напряжений Коши - Эйлера будет выражаться в следующем виде:
(Е)=3 с№ )т (£)-2 (В ).(£).
Скорости изменения напряжений Коши - Эйлера:
(Е )-(Ле)-(^ )+(и )-(Е)+(Е).(Л )т-(Е)/ы,
где введено обозначение:
( ~22ы \
4
(Л^)=т(в)" ^ '(B)•
J
d 2W
dB dB ,
В результате получаем физическое соотношение для производной Трусделла в виде линейного уравнения:
(еТг )-(Е )-(* )-(Е)-(Е)-(А )Т + 4, (Е)-(Ле)"(</).
Таким образом, для гиперупругого тела определяющие соотношения для скоростей имеют вид линейной зависимости при использовании обобщенной производной Трусделла для тензора истинных напряжений Коши - Эйлера и тензора деформации скорости.
Алгоритм расчета
Для решения задачи используется метод последовательных нагружений. Считается, что известно к -е состояние, по которому нужно найти (к +1) -е состояние.
В качестве базового уравнения используется уравнение виртуальных мощностей, записанного для (к +1) -го шага:
§( «е)^) -Л/ "+./ и + ( к+'?„ и Ук+1 Ук+1 ( где Ук+1 - текущий объем; $к+1 - часть его поверхности, на которой заданы усилия; / , Тп - векторы массовых и поверхностных сил.
Переходя к приращениям
(*+12) = ( * 2) + (л * 2),
получим разрешающее уравнение, решение которого дает вектор перемещении для текущего шага Л*и = Л, с помощью которого определяется конфигурация следующего шага:
*+1Я = кЯ + Ли.
Тогда напряженное состояние находится по соотношению: (*+12) = — (*+15^
дЖ
Рассмотрен пример построения физических соотношений для потенциала упругих деформаций Муни - Ривлина:
К
Ж = и1 (11В - 3) + и2 (12в - 3) + у / -1)2,
где Ж0/) = -1)2 , Ж'(Iв, 12В) = и1(11В -3) + и2(12В -3), и1,и2 - постоянные материала.
Тензор напряжений Коши - Эйлера для этого материала имеет вид:
5 г 7
(В)-311В (I)
+ 2и2/ 3
I
1
(в)-31\е (I)-2 (в2)
+К с -ое)
Линеаризуя полученное соотношение, получим выражение для скорости изменения напряжений Коши - Эйлера:
(Л2) = (А2)-(^ ) + (А )-(2) + (2)-(А )Т-(2)^,
где
(ЛЕ) = (ЛЕ') + (Л,0) Л2' = / (в)
Гд 2Ж1 ^ дВдВ
V J
°0 / >
^д 2Ж0 > дВдВ
V J
Таким образом, получена система уравнений для определения напряженного состояния, которая может быть дискретизирована методом конечных элементов для получения численных решений.
Численный пример
Для проверки работоспособности методики был решен ряд задач. Рассмотрим некоторые из них.
1. Задача о плоской деформации
квадратной полосы со сторонами 0.02х 0.02 м2, и1 = 0.36МРа, и2 = 0,22МРа ,
К = 2000МРа. По вертикальным краям полосы заданы перемещения. Так как полоса имеет две оси симметрии, то была рассмотрена четверть полосы с заданием соответствующих условий симметрии [20].
0.02 м
N 2
Рисунок 2. Плоское деформирование квадратной полосы
Рисунок 4. Конечно-элементная модель плиты
2. Задача об упругом
деформировании плиты под действием
равномерного давления ^ = 85 . Верхнее ребро плиты не имеет вертикального смещения. Плита квадратная со стороной а = 0.02 м и
толщиной И = 0.01 м, к = 2000 °ра . На рисунке 4 изображено недеформированное состояние плиты.
На рисунках 5, 6 изображено поле вертикальных перемещений и поле касательных напряжений.
I-
L
Рисунок 5. Поле вертикальных перемещений плиты
Рисунок 6. Поле касательных напряжений плиты
Заключение
Таким образом, была построена методика численного исследования гиперупругих слабосжимаемых материалов, для которых физические соотношения задаются с помощью
упругого потенциала. Получены линеаризованные определяющие соотношения и разрешающее уравнение для потенциала Муни - Ривлина. Данная методика после конечно-элементной дискретизации может быть использована для определения НДС, что видно из рассмотренных примеров. Решенные задачи демонстрируют эффективность полученной методики.
Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов №12-01-00955, №12-01-97026, №12-01-31212, №13-01-97059, №13-01-97058.
Литература
1. Голованов А.И., Султанов Л.У. Математические модели вычислительной нелинейной механики деформируемых сред. Казань: КГУ, 2009. 465 с.
2. Оден Д. Конечные элементы в нелинейной механике сплошных сред. М.: Мир, 1976. 464 с.
3. Черных К.Ф. Нелинейная теория упругости в машиностроительных расчетах. Л.: Машиностроение, 1986. 336 с.
4. Васидзу К. Вариационные методы в теории упругости и пластичности. М.: Мир, 1987. 542 с.
5. Грин А., Адкинс Д. Большие упругие деформации и нелинейная механика сплошной среды. М.: Мир, 1965. 455 с.
6. Коробейников С.Н. Нелинейное деформирование твердых тел. Новосибирск: Изд-во СО РАН, 2000. 262 с.
7. Трусделл К. Первоначальный курс рациональной механики сплошных сред. М.: Мир, 1975. 592 с.
8. Гузь А.Н. Основы трехмерной теории устойчивости деформируемых тел. Киев: Вища школа, 1986. 511 с.
9. Елисеев В.В. Механика упругих тел. СПб: СПбГТУ, 1999. 341 с.
10. Конюхов А.В., Коноплев Ю.Г. Модель термогиперупругости и ее применение к исследованию потери устойчивости раздуваемых пластин. I // Известия вузов. Авиационная техника. 2006. №3. С. 12-16.
11. Чернышов А. Д. Простые определяющие уравнения для упругой среды при конечных деформациях // Известия РАН. Механика твердого тела. 1993. №1. С. 75-81.
12. Simo J.S., Pister K.S. Remarks on rate constitutive equations for finite deformation problems: computational implications // Computer Methods in Applied Mechanics and Engineering. 1984. Vol. 46. Issue 2. Pp. 201-215.
13. Simo J.S., Ortiz M. A unified approach to finite deformation elastoplastic analysis lased on the use of hyperelastic constitutive equations // Computer Methods in Applied Mechanics and Engineering. 1985. Vol. 49. Issue 2. Pp. 221-245.
14. Панов А.Д. Теория определяющих соотношений при деформировании изотропного твердого тела // Известия РАН. Механика твердого тела. 2004. №6. С. 27-44.
15. Vujosevic L., Lubarda V.A. Finite strain thermoelasticity based on multiplicative decomposition of deformation gradient // Theoretical and applied mechanics. 2002. Vol. 28-29. Pp. 379-399.
16. Liu C.H., Wong J.Y., Mang H.A. Large strain finite element analysis of sand: model, algorithm and application to numerical simulation of tire-sand interaction // Computers & Structures. 2000. Vol. 74. Issue 3. Pp. 253-265.
17. Maniatty A.M., Liu Y., Klaas O., Shephard M. Higher order stabilized finite element method for hyperelastic finite deformation // Computer Methods in Applied Mechanics and Engineering. 2002. Vol. 191. Issues 13-14. Pp. 1491-1503.
18. Голованов А.И., Коноплев Ю.Г., Кузнецов С.А., Султанов Л.У. Численное моделирование больших деформаций неупругих трехмерных тел // Наукоемкие технологии. 2004. №4(5). С. 52-60.
19. Голованов А.И., Султанов Л.У. Численное исследование больших упругопластических деформаций трехмерных тел // Прикладная механика. 2005. №6(41). С. 36-43.
20. Bonet J., Wood R.D. Nonlinear continuum mechanics for finite element analysis. New York: Cambridge University Press, 1997. 283 p.
Ленар Усманович Султанов, г. Казань, Россия Тел. раб.: +7(843)233-71-85; эл. почта: ls561@mail.ru
Ленар Рустамович Фахрутдинов, г. Казань, Россия Тел. раб.: +7(843)231-51-85; эл. почта: bishchumbek@gmail.com
©Султанов Л.У., Фахрутдинов Л.Р., 2013
Султанов Л.У., Фахрутдинов Л.Р. Численное исследование гиперупругих материалов
doi: 10.5862/MCE.44.9
Numerical study of hyperelastic materials
L.U. Sultanov,
Kazan (Volga region) federal university, Kazan, Russia +7(843)2337185; e-mail: ls561@mail.ru L.R. Fakhrutdinov, Kazan (Volga region) federal university, Kazan, Russia +7(843)2315185; e-mail: bishchumbek@gmail.com
Key words
hyperelastic materials; large deformations; finite-element method
Abstract
In this paper, the research method of large deformations of hyperelastic solids using the Finger strain measure is considered.
The kinematics of continua motion is given; the stress state is described by the Cauchy - Euler tensor of the true stresses. Much attention is given to the algorithm of the linearized constitutive equations in terms of the Cauchy - Euler stresses. Also, there is an example of obtaining the linearized physical ratio in the Cauchy - Euler stresses for material, which is described by Mooney - Rivlin potential. The numerical implementation is based on a finite-element method within the framework of incremental methods.
To verify efficiency of the technique, several problems were solved: the problem of plane strain of the square strip and the problem of elastic deformation of the plate under uniform pressure. The results obtained do not contradict those that have been published before.
References
1. Golovanov A.I., Sultanov L.U. Matematicheskiye modeli vychislitel'noy nelineynoy mekhaninki deformiruyemykh sred [Mathematical models of computational nonlinear mechanics of deformable continua]. Kazan: Kazan state university, 2009. 465 p. (rus)
2. Oden D. Konechnye elementy v nelineynoy mekhanike sploshnykh sred [Finite elements in nonlinear continuum mechanics]. Moscow: Mir, 1976. 464 p. (rus)
3. Chernykh K.F. Nelineynaya teoriya uprugosti v mashynostroitelnykh raschetakh [Nonlinear theory of elasticity in engineering calculations]. Leningrad: Mashinostroyeniye, 1986. 336 p. (rus)
4. Vasidzu K. Variatsionnyye metody v teorii upugosti i plastichnosti [Variational methods in the theory of elasticity and plasticity]. Moscow: Mir, 1987. 542 p. (rus)
5. Grin A., Adkins D. Bolshiye uprugiye deformatsii i nelineynaya mekhanika sploshnoy sredy [Large elastic deformations and nonlinear continuum mechanics]. Moscow: Mir, 1965. 455 p. (rus)
6. Korobeynikov S.N. Nelineynoye deformirovaniye tverdykh tel [Nonlinear deformations of solids]. Novosibirsk: Sibirskoye otdeleniye RAN, 2000. 262 p. (rus)
7. Trusdell K. Pervonachalnyy kurs ratsionalnoy mekhaniki sploshnykh sred [Initial course of rational continuum mechanics]. Moscow: Mir, 1975. 592 p. (rus)
8. Guz A.N. Osnovy trekhmernoy teorii ustoychivosti deformiruyemykh tel [Basis of three-dimensional theory of stability of deformable solids]. Kiev: Vishcha shkola, 1986. 511 p. (rus)
9. Yeliseyev V.V. Mekhanika uprugikh tel [Mechanics of elastic solids]. Saint-Petersburg: SPbGTU, 1999. 341 p. (rus)
10. Konyukhov A.V., Konoplev U.G. Izv. VUZ. Aviatsionnaya Tekhnika. 2006. No. 3. Pp. 12-16. (rus)
11. Chernyshov A.D. Mechanics of Solids. 1993. No. 1. Pp. 75-81. (rus)
12. Simo J.S., Pister K.S. Remarks on rate constitutive equations for finite deformation problems: computational implications. Computer Methods in Applied Mechanics and Engineering. 1984. Vol. 46. Pp. 201-215.
13. Simo J.S., Ortiz M. A unified approach to finite deformation elastoplastic analysis based on the use of hyperelastic constitutive equations. Computer Methods in Applied Mechanics and Engineering. 1985. Vol. 49. Issue 2. Pp. 221-245.
15. Vujosevic L., Lubarda V.A. Finite strain thermoelasticity based on multiplicative decomposition of deformation gradient // Theoretical and applied mechanics. 2002. Vol. 28-29. Pp. 379-399.
16. Liu C.H., Wong J.Y., Mang H.A. Large strain finite element analysis of sand: model, algorithm and application to numerical simulation of tire-sand interaction. Computers & Structures. 2000. Vol. 74. Pp. 253-265.
17. Maniatty A.M., Liu Y., Klaas O., Shephard M. Higher order stabilized finite element method for hyperelastic finite deformation. Computer Methods in Applied Mechanics Engineering. 2002. Vol. 191. Pp. 1491-1503.
18. Golovanov A.I., Konoplev U.G., Kuznetsov S.A., Sultanov L.U. Science Intensive Technologies. 2004. No. 4(5). Pp. 52-60. (rus)
19. Golovanov A.I., Sultanov L.U. Applied Mechanics. 2005. No. 6(41). Pp. 36-43. (rus)
20. Bonet J., Wood R.D. Nonlinear continuum mechanics for finite element analysis. New York: Cambridge University Press, 1997. 283 p.
Full text of this article in English: pp. 69-74