Моделирование вязкоупругости полиуретана при умеренно высоких скоростях деформирования Текст научной статьи по специальности «Физика»

УДК 539.3

Моделирование вязкоупругости полиуретана при умеренно высоких скоростях деформирования

1 2 1 © А.Е. Белкин , И.З. Даштиев , Б.В. Лонкин

1 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

2 ЦНИИСМ, г. Хотьково Московской обл., 141372, Россия

Представлена математическая модель вязкоупругого поведения полиуретана СКУ-ПФЛ-100 для диапазона деформаций 0...30 % и умеренно высоких скоростей деформирования, не превышающих значения 10-1. Для определения вязкой составляющей деформации применена реологическая модель Бергстрема - Бойс. Связь напряжения с упругой составляющей деформации описана с помощью потенциала Арруды - Бойс. Для определения параметров модели использовались экспериментальные диаграммы сжатия полиуретана, полученные на машине 1т1тп Е1вс1го-риЬ 1000 при различных скоростях деформирования. Приведены значения параметров модели, найденные путем минимизации функции отклонений расчетных величин от результатов эксперимента. Показано, что в рассмотренном диапазоне деформаций и их скоростей модель позволяет описать поведение полиуретана с достаточной для практических целей точностью. Модель предназначена для расчета полиуретановых деталей амортизаторов, поглощающих аппаратов, буферов и других конструкций, испытывающих динамические нагрузки.

Ключевые слова: полиуретан, вязкоупругость, математическая модель Бергстрема - Бойс, упругий потенциал Арруды - Бойс, диаграмма сжатия, скорость деформации, гистерезис, определение параметров модели.

Введение. При анализе работы полиуретановых деталей амортизаторов, поглощающих аппаратов, буферов и других конструкций, испытывающих динамические нагрузки, важное значение приобретает учет зависимости механических свойств полиуретана от режима нагружения, в частности скорости нагружения. Для расчета таких деталей должна привлекаться теория вязкоупругих деформаций эластомеров [1, 2]. Среди формулировок законов вязкоупругости для ре-зиноподобных материалов весьма популярным у специалистов является предложение Бергстрема - Бойс [3]. Исследования [3-5] показывают, что модель Бергстрема - Бойс позволяет достаточно точно описать экспериментальные данные для резины, включая гистерезис при циклическом нагружении и зависимость диаграмм материала от скорости деформирования. Применимость этой модели к полиуретану для процессов с низкими скоростями деформирования изучалась в работе [6], где сделан вывод о ее хорошей предсказательной способности. В настоящей статье анализируется применение модели Берг-стрема - Бойс для описания вязкоупругого поведения полиуретана СКУ-ПФЛ-100 при умеренно высоких скоростях деформирования.

В

Числовые значения параметров модели определяются по результатам испытаний образцов на сжатие с различными скоростями нагружения.

^ Математическая модель вяз-коупругости полиуретана. Согласно [3], поведение вязкоупругого ма-_2 териала подобно поведению условной механической системы (рис. 1), состоящей из двух параллельно соединенных звеньев A, B разной природы. Упругое звено A определяет Рис. 1. Условная модель вязкоупруго- поведение при медленном статиче-го материала ском нагружении, когда демпфер

вязкоупругого звена B не сопротивляется деформации и не воспринимает напряжения. В случае конечной скорости деформации работают оба звена системы, причем сопротивление демпфера зависит от скорости деформации. При очень быстром динамическом (мгновенном) нагружении демпфер «запирается», в нем не происходит деформации и напряжение распределяется между упругими элементами звеньев A и B. В этом случае реакция системы является упругой, однако жесткость повышается по сравнению с медленным упругим деформированием.

Выбор потенциала для упругих элементов звеньев, а также закона изменения вязких деформаций демпфера определяет поведение модели.

При параллельном соединении структур (звеньев) их деформации одинаковы:

FA = FB = F, (1)

здесь F — градиент вектора места (градиент деформации).

Напряжения в материале складываются из откликов структур (звеньев):

о = о A + о B. (2)

Статическая составляющая напряжения оA определяется соотношениями гиперупругости на основе выражения для удельной энергии деформации материала

WA = "2 k (J -1)2 + WA (7Ю, 72 c ), (3)

где k — модуль объемного сжатия; J = det F — мера объемной деформации; WA (71c , 7 2c ) — энергия изохорической деформации;

7 ic = tr C, 7 2 c = 1 (liC - C : C j — соответственно первый и второй

инварианты тензора меры изохорической деформации C = J -2/3 Fт F.

В соответствии с выражением (3) второй тензор напряжений

дЖА

Пиолы-Кирхгофа 8а = 2-имеет вид [7]

дС

8 А = Ы(ы - 1)С"1 + 2фА (I - 3 11с С"1) + 2фА (11с I - С - 212с С-1), (4)

дШа дША

где фА = Ы-2/3; ФА = Ы-4/3; С =ГтГ; 11с = гг С; 12с =

д 11С д12С

= 1 (/1с - С : С).

Переходя к истинным напряжениям Коши по формуле о а = Ы-1 Г8А Гт, можно получить

оА = к(ы -1) I + 2ы -1(фА + /1С фА) йеу В - 2ы -1 ф^^ВВ), (5)

где В = ГГт — тензор меры деформации Фингера [7], символ ёеу обозначает девиатор тензора.

Конкретизируем вид функции удельной энергии деформации полиуретана. Как известно, для эластомеров предложено большое число вариантов определения этой функции, выбор сложен. Весьма успешным, способным описать экспериментальные данные при различных видах напряженного состояния, признается потенциал Арру-ды - Бойс [8] для цепочечной модели эластомера

WA = CA

XlU ß

V

(6)

где CA — коэффициент, пропорциональный начальному значению модуля сдвига эластомера; X chain = V!\с / 3 — кратность усредненного удлинения макромолекулярной цепи эластомера; X Alock — предельное удлинение молекулярной цепи; ß = L-1(X chain / XAock) —

функция, обратная к функции Ланжевена L(ß) = cth(ß) -ß-1.

Модель материала Арруды - Бойс включена в конечно-элементные комплексы ANSYS, ABAQUS, MARC, что способствует ее популярности.

Воспользуемся потенциалом (6) для определения напряжений. При этом учтем, что для обратной к функции Ланжевена существуют различные аппроксимации [9,10]; в настоящей работе применим простейшую из них:

L-1(x) = x

С 3 - x2 ^

V1 - x2 у

x = X chain / X 'lock.

При потенциале (6) соотношения упругости принимают вид

Л k(Л 1)тн GA L~l(Xchain/ XAock) d B (7)

oA = k (J -1) T н-----dev B , (7)

J X chain L-1( 1 / X Ажк)

CA -

где GA =-L-1(1 / X Aock) — постоянная материала, имеющая зна-

3 X iock

чение начального модуля сдвига; B = J -2/3 B — тензор меры Финге-ра для изохорической деформации.

Закон упругости Арруды - Бойс (7) содержит параметры GA, XAock, k, определяемые экспериментально.

Определение динамической составляющей напряжения 0B основывается на принципе мультипликативного разложения градиента полной деформации на упругую Fe и вязкую Fv части

F = Fe Fv (8)

с последующей формулировкой закона течения для скоростей вязких деформаций.

Пользуясь представлением (8) и составляя выражение для пространственного градиента скорости [2, 7], можно получить следующее кинетическое уравнение:

F F -1 = Fe Fe"1 + Fe Dv Fe"1, (9)

где Dv — тензор скоростей вязких деформаций.

Тензор Dv выражается по закону течения через девиатор напряжений

3 S

Dv =- dev оB, (10)

2 Gv

где s v, av — интенсивности скоростей деформаций и напряжений в вязкой среде.

Для установления зависимости sv = f (av,Xvchain) в работе [3] использована физическая концепция движения макромолекул эластомера в так называемом полимерном расплаве [11]. Исходя из этой концепции и привлекая результаты обширных экспериментов на резинах различных марок, Бергстрем и Бойс предложили следующую аппроксимацию закона:

A Gvm

(XVchain - 1 + S0)n

(11)

где А, т , п — параметры закона деформирования; 50 — малая постоянная деформация, добавляемая, чтобы описать скорость ползучести при нулевой деформации.

Заметим, что соотношение (11) — это уравнение теории упрочнения [12], записанное для полимерных цепочек. Используем это соотношение для расчета полиуретана.

Из кинетического уравнения (9) с учетом закона течения (10), (11) можно выразить скорость упругой деформации

К = Г Г -1 К -

А К)

В\т-1

2(№сЫп -1 + §0)п

-К ёеу ов.

(12)

Если предположить, что градиент деформации Я и его скорость Г фиксируются в ходе эксперимента и являются известными функциями времени, то соотношение (12) можно рассматривать как дифференциальное уравнение относительно градиента упругой деформации Ге. При этом динамическая составляющая напряжения ов выражается через тензор Ге.

Примем, что напряжение ов подчиняется тому же закону упругости Арруды-Бойс, что и напряжение оА равновесного состояния

^ 1(Хе сИсап 7 ХВоск ), -

Ов = к(Л -1)I +•

0В

^е ХВ скат ^ 1(1 7 Х Воск )

•ёеу Ве.

(13)

Интегрирование дифференциального уравнения (12) должно выполняться с начальным условием

и=0

I.

(14)

Одноосное напряженное состояние. Для определения числовых значений параметров модели, как правило, используются результаты испытаний материала при одноосном напряженном состоянии. Деформации материала при одноосном растяжении-сжатии характеризуется значениями кратностей удлинения-укорочения в продольном и поперечном направлениях, обозначаемых ниже Х1 и X2 соответственно. Для этого состояния составим матрицы из компонент тензоров В, В, ёеу В:

0 0 "(V Х2 )4 7 3 0 0

[В] = 0 Х 22 0 , [ В]= 0 (V X 2 )-2 7 3 0

0 0 Х22 0 0 (1/ X 2 )-2 7 3

[devB]= (V X2 )4/3-(VX2)

\-2 / 3

3 о

о - 3

о о

о

0 -1

(15)

Из соотношений упругости (7) с учетом определений (15) можно получить выражение для действующего напряжения растяжения-сжатия

aA =-

ga L-1(X

chain

X Aock )

J X

chain

L-1(1 / X U )

4 / 3

-2 / 3

(16)

При использовании формулы (16) следует считать, что кратность удлинения в поперечном направлении X2 является величиной, зависимой от кратности удлинения в продольном направлении Х1, и определяется через Х1, исходя из условия о А = 03 = 0. Такой подход к определению X 2 предполагает известным значение модуля объемного сжатия материала. К сожалению, этот параметр для эластомеров, в частности полиуретанов, редко известен с достаточной точностью. Однако при одноосном напряженном состоянии малая объемная сжимаемость материала практически не влияет на податливость образца. Поэтому при описании экспериментальных диаграмм растяжения-сжатия будем рассматривать полиуретан как несжимаемый материал, полагая, что объемная деформация отсутствует, т.е.

2 _1 / 2

3 = Х1Х2 =1. Подставляя X2 = Х1 в соотношение упругости (16), получим

aA =

ga L-1(X

chain

X iock )

X

chain

L-1(1 / X U)

(X12 -Xf1),

(17)

где

Xchain ( Xf2 + 2Xf 1).

Закон упругости (17) является двухпараметрическим. Для определения динамической составляющей напряжения оВ используем тот же закон упругости (17), в котором полная кратность удлинения материала Х1 заменена на кратность упругого удлинения

Xв ветви В (см. рис. 1) и параметры имеют новые значения Ов,

Х в

¡век"

В случае одноосного напряженного состояния уравнение (12) приводится к виду

i i B A oB\m

T 1 e A G1 ■ s B\ /1 o\

7--TT =-1-П sign(^B). (18)

i1 i

1e

Kokain -1 + 50

Путем численного интегрирования уравнения (18), дополненного соотношениями упругости о А (к^), оВ (к ^), можно воспроизводить

диаграммы растяжения-сжатия материала.

Если испытания образцов проводятся с заданным законом нагру-жения, то условное напряжение в материале (отношение нагрузки на образец к начальной площади его поперечного сечения) является известной функцией времени. В несжимаемом материале при одноосном напряженном состоянии истинное напряжение а1 и условное напряжение а* связаны соотношением

01 = к а*. (19)

Подставляя в (19) а1 =аА (к1) + аВ (к^) и дифференцируя полу-

ченное уравнение по времени, находим

дХ1

i 1 +

дК7

i Be =i 1 a*l+i1 ä*. (20)

При установленных связях oA (A4), gB (i fe) система дифференциальных уравнений (18), (20) может быть численно проинтегрирована по времени. Таким образом, для заданного закона нагружения а* (t) определяются кратности удлинений i1, iBe как функции времени.

Результаты испытаний образцов полиуретана СКУ-ПФЛ-100 на сжатие. Для анализа работы полиуретановых элементов амортизаторов важное значение имеют исследования зависимости механических свойств полиуретана от скорости нагружения. Авторами проведена серия испытаний образцов полиуретана СКУ-ПФЛ-100 при низких скоростях деформирования на машине Zwick 100 (статические испытания) и при умеренно высоких скоростях деформирования на машине Instron Electropuls 1000 (динамические испытания). В испытаниях использовались короткие цилиндрические образцы с отношением высоты к диаметру 1,6.

В статических испытаниях для получения стабильных характеристик осуществлялось циклическое нагружение. Форма цикла нагру-жения была треугольная, т. е. нагрузка и разгрузка происходили с постоянными, причем равными скоростями. Фаза отдыха между циклами нагрузка-разгрузка отсутствовала. Средняя скорость деформирования составляла 0,0022 с . Установление диаграммы сжатия наблюдалось фактически после двух-трех циклов.

На рис. 2 приведены диаграммы пяти циклов сжатия в координатах «условное напряжение - относительное укорочение».

20 30 40

Относительное укорочение, %

Рис. 2. Диаграммы сжатия полиуретана СКУ-ПФЛ-100 на первых пяти циклах

Для каждого цикла определены значения удельной накопленной, возвращенной и рассеянной энергии, а также коэффициент механических потерь (относительный гистерезис) (табл. 1). Как можно видеть, уже начиная со второго цикла, энергетические характеристики циклов различаются весьма незначительно, что свидетельствует об установлении процесса деформирования.

Таблица 1

Энергетические характеристики циклов нагружения

Цикл Накопленная Возвращенная Рассеянная Относительный

нагружения энергия, Дж/см3 энергия, Дж/см3 энергия, Дж/см3 гистерезис

1 3,353 2,490 0,863 0,26

2 3,097 2,469 0,628 0,20

3 3,061 2,459 0,602 0,20

4 3,041 2,452 0,589 0,19

5 3,028 2,447 0,581 0,19

По диаграммам сжатия (см. рис. 2) определены значения начального модуля упругости полиуретана при медленном деформировании. При третьем-пятом циклах нагружения начальный модуль упругости составляет 20...30 МПа для диапазона относительных укорочений 0...10 %.

В динамических испытаниях образцы подвергались периодической импульсной нагрузке. Во всех испытаниях номинальное значе-

ние максимальной сжимающей силы составляло 1 кН (максимальная сила для машины 1п81хоп Е1ес1хоРи18 1000). В использованных образцах диаметром 9,7 мм такая сила вызывает напряжение 13,5 МПа. Скорость нагружения изменялась от 1 кН/с (от 13,5 МПа/с) до 50 кН/с (до 677 МПа/с). В табл. 2 приведены основные параметры динамических испытаний.

Таблица 2

Параметры динамических испытаний на машине 1ш1гоп ЕкйгоРиЬ 1000

Испытание Номинальная скорость нагружения, кН/с (МПа/с) Средняя скорость дефор--1 мирования, с Номинальная максимальная нагрузка, кН Время выдержки под нагрузкой, с Время выдержки после разгрузки в цикле, с Количество циклов

1 1 (13,5) 0,29 0,1 10 10

2 10 (135) 2,26 1

3 35 (474) 6,01 0,01 1 6

4 50 (677) 8,61

В качестве примера одного из реализованных испытаний на рис. 3 представлены графики изменения сжимающей силы и относительного укорочения образца в течение одного цикла нагружения со скоростью 10 кН/с, записанные в процессе испытания.

юоо

500

40

-4=

N №

и &

20

о

В о

0,6 0,8 Время, с

1,0

1,2

Рис. 3. Графики изменения во времени силы сжатия (1) и относительного укорочения (2) в течение одного цикла нагружения

образца

На рис. 4 приведены диаграммы динамического сжатия с различными скоростями. Диаграммы соответствуют установившемуся состоянию, когда каждый последующий цикл нагружения практически повторяет предыдущий.

Относительное укорочение, %

Рис. 4. Диаграммы сжатия полиуретана СКУ-ПФЛ-100 при различных скоростях нагружения и разгрузки, МПа/с: 1 — 13,5; 2 — 135; 3 — 474; 4 — 677

Из диаграмм следует, что с увеличением скорости нагружения возрастает как начальный модуль упругости, так и среднее значение модуля во всем интервале деформаций. В табл. 3 приведены характеристики гистерезисных потерь.

Таблица 3

Характеристики гистерезисных петель, полученных при динамических испытаниях

Испытание Средняя скорость деформиро--1 вания, с Достигнутое условное напряжение, МПа Относительное укорочение, % Накопленная энергия, Дж/см3 Поглощенная энергия, Дж/см3 Относительный гистерезис

1 0,29 13,4 30,5 1,913 0,698 0,37

2 2,26 12,6 26,4 1,689 0,808 0,48

3 6,01 12,4 25,0 1,647 0,876 0,53

4 8,61 13,4 25,0 1,788 0,924 0,52

Можно видеть, что с ростом скорости нагружения в рассмотренном диапазоне наблюдается некоторый рост относительного гистерезиса. Минимальное значение гистерезиса наблюдается при статическом нагружении (см. табл. 1). Полученные результаты испытаний использованы для определения значений параметров модели Берг-стрема - Бойс.

Результаты моделирования. Значения параметров оа , ХАоск

модели Арруды - Бойс (17) подбирались по кривой статического нагружения, показанной на рис. 2. Относительные укорочения ограничивались значением 35 %. Для поиска указанных параметров ис-

пользовался алгоритм минимизации функции отклонений расчетных величин а* (Х1) от результатов эксперимента методом Нелдера - Мида. Независимо от этого для повышения достоверности осуществлялся поиск значений с помощью встроенных средств программы АКБУБ. В результате независимых решений установлены весьма широкие диапазоны значений параметров оа = 4,3... 5,5 МПа, Xаоск = 1,5...2,9. Конкретные значения оа, Xаоск зависят от способа формирования минимизируемой функции отклонений: числа и расположения точек сравнения расчетных и экспериментальных данных, выбора разных весовых коэффициентов для областей малых и больших деформаций.

На рис. 5 показаны расчетные и экспериментальная статические диаграммы сжатия полиуретана. Расчеты выполнены при четырех комбинациях параметров оа , Xаоск. В целом полученные результаты свидетельствуют о применимости модели Арруды - Бойс для анализа упругих деформаций полиуретана. Однако диаграммы на рис. 5 позволяют увидеть и некоторый дефект модели. На начальном участке нагружения, где согласно результатам испытаний наблюдается снижение жесткости полиуретана, расчетные кривые уклоняются от экспериментальных.

12

Ю

№ о ц о

и 6

о о И а о

Е5 л

X

2 1 \

// /г 7 5

/ /Л'

10 20 30

Относительное укорочение, %

40

Рис. 5. Результат описания экспериментальной диаграммы сжатия полиуретана (кривая 1) при помощи модели Арруды - Бойс:

2 - Оа = 5,3 МПа, хАоск = 2,5; 3 — Оа = 5,3 МПа, X Аоск = 2,9; 4 —Оа = 4,3 МПа, XАоск = 2,5; 5 — вА = 4,3 МПа,

X А

= 2,9

Задача калибровки вязкоупругой модели Бергстрема - Бойс, состоящая в нахождении шести параметров: Ов, Xвоск, А, т, п, 50,

существенно сложнее, чем подбор коэффициентов для соотношения упругости (17). Эта задача требует отдельного рассмотрения, в настоящей статье ограничимся представлением некоторых полученных результатов.

Для каждого динамического режима испытаний, указанного в табл. 2, 3, параметры модели подбирались отдельно. Выбранные значения приведены в табл. 4.

Таблица 4

Параметры модели Бергстрема - Бойс и рассчитанные значения относительного гистерезиса

Испытание Ов, МПа Х в 1оск A, МПа - с-1 П1 п 50 Расчетное значение относительного гистерезиса

1 4,8 7,0 0,007 2 0,3 10 -6 0,27

2 8,8 7,0 0,04 2 0,3 10 -6 0,40

3 11,0 7,0 0,09 2 0,3 10 -6 0,45

4 11,8 7,0 0,1 2 0,3 10 -6 0,46

Для оценки точности моделирования при подобранных параметрах восстановлены диаграммы установившегося циклического сжатия полиуретана путем интегрирования уравнений (18), (20). На рис. 6 сопоставлены расчетные и экспериментальные гистерезисные петли при сжатии с различными скоростями нагружения и разгрузки.

5 10 15 20 25 Относительное укорочение, % а

)

//<

/ /

/ /

У. У/ // //

//

5 10 15 20 25 Относительное укорочение, % б

30

Рис. 6 (начало). Сравнение расчетных и экспериментальных диаграмм сжатия полиуретана при различных скоростях нагружения:

а - 13,5 МПа/с; б - 135 МПа/с

о № а

0

1

Л? / /

> г / / / /

✓ ^ / /

// // /

// У /

As у /

5 10 15 20 25 Относительное укорочение, % в

30

Рис. 6 (окончание). Сравнение расчетных и экспериментальных диаграмм сжатия полиуретана при различных скоростях нагру-жения:

в - 677 МПа/с

Как можно видеть, формы гистерезисных петель весьма заметно различаются, наибольшие отличия между расчетом и экспериментом наблюдаются в кривых разгрузки. Вместе с тем величина рассеяния энергии, предсказываемая с помощью модели, хорошо соответствует эксперименту. Показателем этого служит близость экспериментальных и расчетных значений гистерезиса, приведенных в последних столбцах таблиц 3 и 4.

Заключение. Реологическая модель, предложенная Бергстремом и Бойс, позволяет описать поведение полиуретана при сжатии в рассмотренном диапазоне деформаций и их скоростей с достаточной для практических целей точностью. Для отдельного динамического режима удается определить параметры модели, обеспечивающие приближенное согласование с экспериментом. При изменении режима нагружения параметры приходится переопределять. Процедура определения единого набора значений параметров модели, пригодных для описания экспериментальных данных при разных скоростях нагружения, требует дополнительного исследования.

ЛИТЕРАТУРА

[1] Димитриенко Ю.И., Даштиев И.З. Модели вязкоупругого поведения эластомеров при конечных деформациях. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2001, № 1, с. 21-41.

[2] Димитриенко Ю.И. Нелинейная механика сплошной среды. Москва, Физ-матлит, 2009, 610 с.

[3] Bergström J.S., Boyce M.C. Constitutive modeling of the large strain time-dependent behavior of elastomers. Journal of Mechanic Physics Solids, 1998, vol. 46, pp. 931-954.

[4] Bergström J. S., Boyce M. C. Mechanical behavior of particle filled elastomers. Rubber Chem. Technol, 1999, vol. 72, pp. 633-656.

[5] Quintavalla S.J., Johnson S.H. Extension of the Bergström-Boyce model to high strain rates. Rubber Chem. Technol, 2004, vol. 77, pp. 972-981.

[6] Qi H.J., Boyce M.C. Stress-strain behavior of thermoplastic polyurethane. Mechanics of Materials, 2005, vol. 37, issue 8, pp. 817-839.

[7] Голованов А.И., Султанов Л.У. Математические модели вычислительной нелинейной механики деформируемых тел. Казань, Казанский государственный университет, 2009, 465 с.

[8] Arruda E.M., Boyce M.C. A three-dimensional constitutive model for the large stretch behavior of rubber elastic materials. Journal of Mechanic Physics Solids, 1993, vol. 41, no. 2, pp. 389-412.

[9] Cohen A. A Pade approximant to the inverse Langevin function. Rheologica Acta, 1991, vol. 30, pp. 270-273.

[10] Norgan C.O., Saccomandi G. A Molecular-statistical basis for the gent constitutive model of rubber elasticity. Journal of Elasticity, 2002, vol. 68, pp. 167-176.

[11] Дой М., Эдвардс С. Динамическая теория полимеров. Москва, Мир, 1998, 440 с.

[12] Качанов Л.М. Основы теории пластичности. Москва, Наука, 1969, 420 с.

Статья поступила в редакцию 05.11.2014

Ссылку на эту статью просим оформлять следующим образом:

Белкин А.Е., Даштиев И.З., Лонкин Б.В. Моделирование вязкоупруго-сти полиуретана при умеренно высоких скоростях деформирования. Математическое моделирование и численные методы, 2014, № 3, с. 39-54.

Белкин Александр Ефимович окончил МВТУ им. Н.Э. Баумана. Д-р техн. наук, профессор кафедры «Прикладная механика» МГТУ им. Н.Э. Баумана. Автор более 70 научных работ в области численных методов расчета конструкций, механики ре-зинокордных композитов и пневматических шин. e-mail: [email protected]

Даштиев Идрис Зилфикарович окончил МВТУ им. Н.Э. Баумана. Д-р техн. наук, зам. главного конструктора ОАО «ЦНИИСМ». Автор более 125 научных работ и патентов в области механики и технологии полимерных композиционных материалов и конструкций на их основе.

Лонкин Борис Викторович — аспирант кафедры «Прикладная механика» МГТУ им. Н.Э. Баумана. Специализируется в области численных методов расчета конструкций и механики эластомерных материалов. e-mail: [email protected]

Modeling polyurethane viscoelasticity at moderately high strain rates

1 2 1 © A.E. Belkin1, I.Z. Dashtiev2, B.V. Lonkin1

1 Bauman Moscow State Technical University, Moscow, 105005, Russia

2

2 TSNIISM, Khotkovo-town, Moscow Region, 141350, Russia

The article presents a mathematical model of the viscoelastic behavior of polyurethane SKU-PFL-100 for strain range of 0...30 % and moderately high strain rates up to 10 . To determine the viscous component of the deformation Bergstrom - Boyce rheological

model has been applied. Relationship between stress and the elastic component of deformation is described by an Arruda - Boyce potential. We determined the model parameters using experimental compression diagrams of polyurethane obtained from Instron Electropuls 1000 machine at different strain rates. The model parameter values obtained by minimizing a function of the calculated value deviations from the experimental results are given. It is shown that in the considered range of deformations and strain rates model allows describing the polyurethane behavior with sufficient accuracy for practical purposes. The model is designed for calculating polyurethane shock-absorber parts, cushions, buffers and other structures subjected to dynamic loading.

Keywords: polyurethane, viscoelasticity, Bergstrom-Boyce mathematical model, Arruda-Boyce elastic potential, compression diagrams, strain rate, determination of model parameters.

REFERENCES

[1] Dimitrienko Yu.I., Dashtiev I.Z. Vestnic MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2001, no. 1, pp. 21-41.

[2] Dimitrienko Yu.I. Nelineynaya mekhanika sploshnoy sredy [Nonlinear continuum mechanics]. Moscow, Fizmatlit Publ., 2009, 610 p.

[3] Bergström J.S., Boyce M.C. Constitutive modeling of the large strain time-dependent behavior of elastomers. Journal of Mechanic Physics Solids, 1998, vol. 46, pp. 931-954.

[4] Bergström J. S., Boyce M. C. Mechanical behavior of particle filled elastomers. Rubber Chem. Technol., 1999, vol. 72, pp. 633-656.

[5] Quintavalla S.J., Johnson S.H. Extension of the Bergström-Boyce model to high strain rates. Rubber Chem. Technol, 2004, vol. 77, pp. 972-981.

[6] Qi H.J., Boyce M.C. Stress-strain behavior of thermoplastic polyurethane. Mechanics of Materials, 2005, vol. 37, issue 8, pp. 817-839.

[7] Golovanov A.I., Sultanov L.U. Matematicheskie modeli vychislitelnoy nelineynoy mekhaniki deformiruemykh tel [Mathematical models of computational nonlinear mechanics of deformable bodies]. Kazan, Kazan State University Publ., 2009, 465 p.

[8] Arruda E.M., Boyce M.C. A Three-dimensional constitutive model for the large stretch behavior of rubber elastic materials. Journal of Mechanic Physics Solids, 1993, vol. 41, no. 2, pp. 389-412.

[9] Cohen A. A Pade approximant to the inverse Langevin function. Rheologica Acta, 1991, vol. 30, pp. 270-273.

[10] Norgan C.O., Saccomandi G. A molecular-statistical basis for the gent constitutive model of rubber elasticity. Journal of Elasticity, 2002, vol. 68, pp. 167-176.

[11] Doi M., Edwards S.F. The Theory of Polymer Dynamics. Oxford University Press, Oxford, 1986. [Russian edition: Doi M., Edvards S.F. Dinamicheskaya teoriyapolimerov. Moscow, Mir Publ., 1998, 440 p.].

[12] Kachanov L.M. Osnovy teorii plastichnosti [Fundamentals of the theory of plasticity]. Moscow, Nauka Publ., 1969, 420 p.

Belkin A.E. graduated from Bauman Moscow Higher Technical School. Dr. Sci. (Eng.), professor of the Applied Mechanics Department at Bauman Moscow State Technical University. Author of more than 70 publications on numerical methods of design calculation, mechanics of rubber-cord composites and pneumatic tires. e-mail: [email protected]

А.Е. EeßKun, M.3. flawmuee, E.B. HOHKUH

Dashtiev I.Z. graduated from Bauman Moscow Higher Technical School. Dr. Sci. (Eng.), deputy chief designer of the «Central research institute for special machinery» Author of more than 125 publications and patents in the field of mechanics and technology of polymer composite materials and constructions on their base.

Lonkin B.V., a postgraduate of the Applied Mechanics Department at Bauman Moscow State Technical University. Specializes in the field of numerical methods for the calculation of structures and mechanics of elastomeric materials. e-mail: [email protected]