CC BY
85
УДК 539.3 ББК В251
М. С. Воронин, Л. А. Мержиевский
г. Новосибирск, Россия
Модель квазистатического и динамического деформирования эластомеров
В работе строится модель динамического деформирования эластомеров, основанная на максвелловских представлениях о механизмах необратимого деформирования. Построено уравнение состояния при нешаровом тензоре деформации и зависимость времени релаксации касательных напряжений от параметров, характеризующих состояние среды. В рамках сформулированной модели решен ряд задач квазистатического и ударноволнового деформирования резин.
Ключевые слова: динамическая деформация, эластомеры, уравнение состояния, тензор деформации, ударно-волновые процессы.
M. S. Voronin, L. A. Merzhievsky
Novosibirsk, Russia
Model of Quasi-static and Dynamic Deformation of Elastomers
We construct a model of dynamic deformation of elastomers, based on Maxwell’s views on the mechanisms of irreversible deformation. The equation of the state for the not spherical strain tensor and the dependence of relaxation time of shear stresses on the parameters characterizing the state of the environment is derived. In the framework of the formulated model some problems of quasi-static and shock-wave deformation of rubber is solved.
Keywords: dynamic deformation, elastomers, equation of state, the strain tensor, shockwave processes.
Среди органических полимеров в особую группу выделяются эластомеры, имеющие при обычных температурах ярко выраженные высокоэластические свойства. Температура стеклования их существенно ниже 0oC. Наиболее яркими представителями эластомеров являются каучуки и получаемые из них химически сшитые материалы, называемые резинами.
Исследованию химического строения, физических и механических свойств эластомеров, и построению моделей их деформирования и разрушения посвящены работы [1; 6; 7; 11; 17] и др. Развитие теории деформации и разрушения высокополимеров начиналось с описания высокоэластического состояния, которое связывалось гибкостью отдельных макромолекул, изменяющих свою длину под действием внешних сил [6]. Дальнейшее развитие привело к представлениям о сеточной структуре высокополимеров, связанной с образованием физических узлов - временных связей, находящихся в динамическом равновесии, и о существовании надмолекулярных образований, оказывающих существенное влияние на свойства и поведение полимеров [6; 12]. Построение моделей деформирования эластомеров базировалось на описании вязкоупругого характера их поведения, связанного с релаксационными явлениями [4; 12]. При этом развивались максвелловские представления [8], получившие широкое обобщение в работе [3].
В данной статье развитый в работах [9; 10] подход применяется к построению модели и решению ряда задач деформирования каучука и резин.
1. Формулировка модели. Полная механико-математическая формулировка уравнений используемой в работе модели упруговязкого тела максвелловского типа детально изложена в работе [3]. Основные уравнения модели, выражающие законы сохранения импульса, энергии и эволюции компонентов тензора эффективных деформаций, имеют вид:
du 1 да,п dE ал., Bui
— -o, -----——= 0,
dt p drj ’ dt p dr
j
© Воронин М. С., Мержиевский Л. А., 2011
53
^+си'+и»с+К6'+^]=" (1)
Здесь С = ||8^-1| - метрический тензор, ТУ = ||§^-||, ТУо = ^ аЦ ~ компоненты
тензора напряжений, щ - компоненты вектора скорости, р и ро - начальная и текущая плотности, Е - удельная упругая энергия, т - время релаксации касательных напряжений, t и - время и пространственные координаты; г,] = 1, 2, 3; I - единичный тензор,
(I д д п—- дЕ
л-т+щщ- р-^ыв, *,, = -2^—.
Уравнение неразрывности в данной модели следует из приведенных уравнений. Система замыкается уравнением, связывающим изменение удельной внутренней энергии Е с компонентами нешарового тензора деформаций и энтропией £:
Е = Е(^ ,Я) (2)
и зависимостью для времени релаксации касательных напряжений:
т = т (ёц,S), (3)
с помощью которой учитываются микроструктурные механизмы необратимого деформирования. Система соотношений (2), (3) является составным уравнением состояния среды.
Для решения поставленных задач достаточно использовать систему уравнений модели в нестационарном одномерном случае в главных осях координат (предполагается, что главные оси тензоров напряжений и деформаций совпадают), которая имеет вид:
д(ргр) д (ригр) д (риг") д [(ри? — а\) г"]
т + ~0’ а~г +^ а2~0’
д [/э (Е + и2/2) г1'] д [(ри (Е + и2/2) — а\и) г1'] ^
ді + д~г = °’
д^2 дЪ.2 и(3 — V) и ^2 дНз дНз V(V — 1) и ^з
“оГ + м_я----------о----=-------> ~7й~ + м_я---------о-----=------’ (4)
дЬ дг 2 г т дЬ дг 2 г т
Р = Ро ехр(—Н\ — Ь,2 — Нз), <1і = Ні — — 1п —,
3 ро
Е = Е(Ьі, 5і), т = т(Ьі, 5і), аі=Р^~Гі Т=~^Б^
дпі дЬ
где Ні - логарифмы коэффициентов растяжения элемента среды вдоль главных осей, Т - температура элемента среды, V - показатель симметрии: V = 0 - плоская, V =1 - цилиндрическая, V = 2 -сферическая.
2. Построение уравнения состояния. Наиболее удобной формой представления уравнения состояния для (4) является связь изменения удельной внутренней энергии Е с компонентами тензора деформаций и энтропией Ь:
Е = Е(Ні,Н2,Н3,Ь).
Поскольку энергия не зависима от системы координат, то она должна быть связана не с Ні, а с инвариантами тензора. Ограничимся зависимостью энергии только от первого и второго инвариантов тензора деформаций. Как показывает анализ, третий инвариант в рассматриваемых задачах оказывает малое влияние на характеристики среды и протекание процессов. Для выбранного тензора деформаций инварианты (6 - первый, Б - второй) запишутся в виде:
6 = ехр[—(Ні + Н-2 + Нз)], Б = (^2 + ^2 + ^3)-
Таким образом можно строить уравнение состояния в виде:
Е = Е(6, Б, Ь)- (5)
На практике чаще используется построение уравнения состояния в форме зависимости для свободной энергии:
Е = Е (5,В,Т). (6)
Переход от (5) к (6) и другим термодинамическим потенциалам осуществляется по формулам термодинамики:
Е(ё, Д Т) = Е(ё, д т) - т (10 ^ , в{5, А Т) = - ^
(7)
дЕ\ (дЕ
(Ш = —Б(1Т + —аы, (1Е = твв + —аы => = —•
р р \dhijj, \дЬ4) 3 р
По аналогии с классическим уравнением состояния Ми-Грюнайзена, полную энергию (6) представим в виде суммы упругой (или холодной), тепловой и девиаторной составляющих:
Е = Ех + Еаву + Ет. (8)
Конкретный вид зависимости каждого слагаемого в (8) выбирался с использованием имеющихся экспериментальных данных и уравнений состояния разных типов. Форма зависимости для холодной составляющей выбирается либо на основе потенциала взаимодействия атомов и молекул, либо с привлечением феноменологических представлений, базирующихся на косвенных экспериментальных данных. По аналогии с уравнением состояния РОСА [12] в качестве зависимости упругой составляющей Ех примем
/адт-1 Ып-1\
ад) = ^0 (
— 1 п — 1 )
Рх(6) = а6т + Ъ6п-Рхо, <5= А (9)
ро
где а, Ь, т, п - имеют смысл интерполяционных констант, а константы Ехо и Рхо необходимы для выполнения условий Ех (1) =0 и Рх (1) = 0. Все эти параметры находились из наилучшего соответствия с уравнением состояния [13].
В девиаторном слагаемом ограничимся линейной зависимостью внутренней энергии единицы массы вещества от второго инварианта В:
Еаеу (5,В) = 2е\5& В, (10)
где е± - поперечная скорость звука, £о - интерполяционная константа.
Тепловое слагаемое строилось на основе имеющейся в литературе информации о теплоемкости полимеров [14]. Анализ совокупности теоретических и экспериментальных работ, посвященных рассмотрению структурных механизмов (колебательных спектров), определяющих теплоемкость полимеров, и непосредственному измерению соответствующих параметров, позволил остановиться на следующей зависимости для тепловой составляющей уравнения состояния:
3 1
ЕТ(6,Т) = ДтТ]Гл^1п(1 -ех*) - -тт2,
i=1
Хг=в4^, в1=вЫРа*, №=Р о5-т% (И)
где Ят - универсальная газовая постоянная, делённая на молярную массу, 9оi - харакеристиче-ские температуры, определяемые экспериментально, Yоi - интерполяционные константы, имеющие смысл аналогов коэффициента Грюнайзена, для соответствующих колебательных мод, Ni - интерполяционные константы, удовлетворяющие условию ^ Ni = N, где N - полное число колебательных мод, во и Ге - интерполяционные константы, которые отличны от нуля в случае необходимости учета теплового возбуждения электронов. Входящие в приведенное соотношение параметры определяются с привлечением данных из [13; 18] и сравнения расчетной и экспериментальной ударной адиабаты.
3. Построение зависимости для времени релаксации касательных напряжений.
Впервые, по-видимому, идея о релаксационном характере необратимого деформирования была сформулирована Д. К. Максвеллом [8]. В качестве иллюстрации было выписано соотношение, получившее в дальнейшем широкую известность как определяющее уравнение реологической модели Максвелла. В рамках существовавших в то время представлений механизм необратимой деформации и сопутствующей релаксации напряжений любой конденсированной среды рассматривался как макроскопический результат перегруппировок её молекул (атомов) и перемещений разнообразных комплексов этих молекул к такому (стабильно-равновесному) расположению, при котором энергия взаимодействия совокупности молекул, образующих рассматриваемое тело, и, соответственно, его свободная энергия имеют минимальное значение. Это однозначно коррелирует с современными представлениями о микроструктурных механизмах высокоэластической деформации полимеров. Исходя из термоактивируемого характера микроструктурных механизмов [1], для времени релаксации примем классическое соотношение Больцмана-Аррениуса
т = то ехр(и/ЯТ),
где то - некоторое характерное время кинетической единицы, за счет движения которой происходит релаксация, и - энергия активации соответствующего механизма релаксации.
С учетом изменения энергии активации в процессе деформирования, а также с учетом возможности протекания нескольких механизмов релаксации это соотношение принимает вид
,цр,л-Г’»«»(Р”~ет(Т)Е)’
І ' '
ai(T) = aоi • (1 + ацАТ), АТ = Т — То, (12)
£ = У^ [(<71 - сгг)2 + (<71 - стз)2 + (<72 - 03)2],
где aоi, ali - интерполяционные константы, Я - интенсивность касательных напряжений.
Для определения входящих в зависимость (12) параметров воспользуемся разработанной в [9] методикой. Будем решать задачу об одноосном деформировании тонкого стержня. В предположении об однородности деформации все величины в исходной системе уравнений (4) являются функциями только времени. Поэтому в итоге система упрощается и принимает вид:
сШ1 . 2 hl — h2 <1Б 2 hl----h2 71
Л 3 г ’ Л 3 г рТ ’
Е = Е(3, В, Т), т = т(6,В,Т), 72 = 73 = 0, (13)
/12 = Ьз, 5 = ехр(-Ь 1-2/12), О =-(Ьх - к?)2,
Результатом решения этой задачи являются диаграммы деформирования. Значения параметров в (12) определяются на основе условия наилучшего совпадения расчетных и экспериментальных диаграмм деформирования.
4. Результаты. На рис. 1 линиями приведены результаты расчета диаграмм деформирования каучука в квазистатическом диапазоне. Здесь кривые 1-3 соответствуют е = 0.075, 0.225, 0.96 с-1. Символами показаны экспериментальные данные из [14]. Как следует из приведенных на рисунке данных, при е ~ 40% происходит резкое изменение наклона диаграмм. Это вызвано изменением механизма структурной релаксации полимера, связанной с изменением характера микроструктурного механизма необратимого деформирования - активацией молекулярных комплексов более высокого порядка. В построенной модели это передается изменением параметров в зависимости для времени релаксации касательных напряжений.
На рис. 2, 3 расчетные диаграммы деформирования (сплошные линии) сравниваются с экспериментальными данными [16] (символы) резин двух разных марок. На рис. 2 кривые 1-3 соответствуют скоростям деформации е = 1200, 1800, 3000 с-1, кривые 1-3 на рис. 3 получены для
ё = 900, 2500, 3000 с-1 соответственно. Здесь момент изменения механизма релаксации выражен менее четко, и структурный переход наблюдается в диапазоне деформаций ё ~ 20-40%. В модели данное обстоятельство учтено тем, что (12) представляется в форме двух слагаемых с соответствующими значениями параметров. Значения этих параметров обеспечивают постепенный переход от описания механизма деформации более низкого уровня к более высокому.
Для ряда практических задач необходимо описание поведения эластомеров при ударноволновом нагружении. Резиноподобные материалы зачастую используются для демпфирования импульсов ударного сжатия и в других практических случаях. С целью проверки применимости построенной модели для решения таких задач проведены расчеты ударной адиабаты и адиабат разгрузки резины. Результаты расчетов приведены на рис. 4 (зависимости в координатах гидростатического давления Р от массовой скорости и, 1 - ударная адиабата, 2-4 - адиабаты разгрузки, пунктир - расчет по уравнению состояния [13], символы - экспериментальные данные [15]). На рис.
5, 6 приведены эти же данные в координатах: скорость ударной волны Б - массовая скорость и давление - удельный объем V соответственно (использованы те же обозначения, что и на предыдущем рисунке).
Важной задачей, решение которой является тестом для проверки применимости моделей для описания реальных ударно-волновых процессов, является задача о взаимодействии ударной волны с догоняющей волной разрежения. Этот процесс часто реализуется в экспериментальных исследованиях и на практике. Следует указать, что численные расчеты показывают [10], что многие из моделей, традиционно применяющихся к описанию ударно-волновых процессов, неверно передают характер затухания ударных волн при их взаимодействии с догоняющими волнами разрежения. Для тестирования построенной модели на адекватность описания таких процессов решалась задача, постановка которой соответствовала экспериментальному исследованию [5]. Плоские ударные волны в образцах из вакуумной резины создавались ударом пластин из алюминия или оргстекла, или детонацией зарядов взрывчатого вещества, контактирующего с образцом. Ввиду недостатка в [5] информации о параметрах ударников при расчете в качестве начальных условий задавался экспериментально измеренный вблизи поверхности нагружения ударный импульс (на приведенных ниже рисунках показаны цифрами). Профили рассчитанных волн в нескольких лагранжевых точках (сплошные линии) сравниваются с экспериментальными (символы) на рис. 7, 8. Проведенное сравнение показывает, что построенная модель с хорошей точностью воспроизводит результаты экспериментальных исследований и может применяться для решения задач ударно-волнового деформирования.
Рис. 1.
Рис. 2.
Рис. 5.
Рис. 6.
Рис. 7.
Рис. 8.
Список литературы
1. Бартенев Г. М., Френкель С. Я. Физика полимеров. Ленинград: Химия, 1990. 432 с.
2. Глушак Б. Л., Гударенко Л. Ф., Стяжкин Ю. М. Полуэмпирическое уравнение состояние металлов с переменной теплоемкостью ядер и электронов // Вопросы атомной науки и техники. Сер. «Математического моделирования физических процессов». 1991. в. 2. C. 57-62.
3. Годунов С. К. Элементы механики сплошной среды. М.: Наука, 1978. 304 с.
4. Гуревич Г. И. О законе деформации аморфных и поликристаллических тел // ЖТФ. 1947. T. 17. № 12. C. 1491-1502.
5. Калмыков Ю. Б., Канель Г. И., Пархоменко И. П., Уткин А. В., Фортов В. Е. Поведение резины в ударных волнах и волнах разрежения // ПМТФ. 1990. № 1. C. 126-130.
6. Лебедев В. П. Структура аморфных полимеров // Успехи химии. 1978. T. 47. № 1. C. 127-151.
7. Лукомская А. И., Евстратов В. Ф. Основы прогнозирования механических свойств кау-чуков и резин. М.: Химия, 1975. 360 с.
8. Максвелл Д. К. Строение тел. // Д. К. Максвелл. Статьи и речи. М.: Наука, 1968. C. 182192.
9. Мержиевский Л. А., Шамонин С. А. Построение зависимости времени релаксации касательных напряжений от параметров состояния среды // ПМТФ. 1980. № 5. С. 170-179.
10. Мержиевский Л. А., Реснянский А. Д. Численное моделирование ударно-волновых процессов в металлах // ФГВ. 1984. Т. 20. № 5. C. 114-122.
11. Уорд И. Механические свойства твердых полимеров М.: Химия, 1975. 360 с.
12. Френкель Я. И., Образцов Ю. Н. Феноменологическая теория механических свойств аморфных тел и распространение колебаний в них // ЖЭТФ. 1939. T. 9. в. 9. C. 1081-1093.
13. URL: http://www.ficp.ac.ru/rusbank/ (дата обращения 30.04.2011).
14. Amin A. F., Alam M. S., Okui Y. An improved hyperelasticity relation in modeling viscoelasticity response of natural and high damping rubbers in compression: experiments, parameter identification and numerical verification // Mechanics of Materials. 2002. Vol. 34. P. 75-95.
15. S. P. Marsh (Ed.) LASL Shock Hugoniot Data. Univ. California Press. Berkeley, 1980.
16. Yang L. M., Shim V. P. W., Lim C. T. A visco-hyperelastic approach to modelling the constitutive behaviour of rubber // International Journal of Impact Engineering. 2000. Vol. 24. P. 545-560.
17. Treloar L. R. G. The mechanics of rubber elasticity // Proc. Roy. Soc. London. 1976. A 351. P. 295-406.
18. Wunderlich B. Thermal Analysis of Polymeric Materials. Springer Berlin Heidelberg, 2005.
Рукопись поступила в редакцию 10 мая 2011 г.
