УДК 677.017.56
ОПИСАНИЕ ВЯЗКОУПРУГОСТИ ПОЛИМЕРНЫХ МАТЕРИАЛОВ СТАТИСТИЧЕСКИМ МЕТОДОМ
О.В. Рымкевич1, П.П. Рымкевич2, А.А. Романова3
1 Санкт-Петербургский государственнъшуниверситет технологии и дизайна (СПбГУТД),
191186, Санкт-Петербург, ул. БолъшаяМорская, 18 2,3Санкт-Петербургский государственнъшуниверситет сервиса и экономики (СПбГУСЭ),
191015, Санкт-Петербург, ул. Кавалергардская, 7, лит. А
Предложено описание вязкоупругости полимерных материалов, используемых для защиты автомобильных деталей, статистическим методом. Получено уравнение, позволяющее описать режимы деформации в различных температурных полях.
Ключевые слова: термоусаживающаяся трубка, эффект памяти формы, вязкоупругость, коэффициент продольной и поперечной усадки.
VISCOELASTICITY OF POLYMERIC MATERIALS DESCRIPTION BY STATISTICAL
METHOD
O. V. Rymkevich, P. P. Rymkevich, A.A.Romanova
Technology and design St. Petersburg State University, 191186, St. Petersburg, Bolshaya Morskaya St., 18 St. -Petersburg state university of service and economy (SPbSUSE), 191015, St.-Petersburg, streetKavalergardsky, 7, lit. A
Description of the polymeric materials viscoelasticity by the statistical method used for automobile parts defence is suggested. An equation for description of the deformation modes in different temperature fields is received.
Keywords: Heat-shrinkable tubing, the shape memory effect, viscoelasticity, ratio of the longitudinal and transverse shrinkage.
Автомобильные детали со сложным профилем в режиме эксплуатации подвергаются негативным воздействия окружающей среды. Поэтому очень важно для повышения их износостойкости, а также для улучшения изоляционных свойств защищать их специальными термоусаживающимися трубками. Эти трубки изготавливаются из, так называемых, "умных материалов” - материалов с эффектом памяти формы [1 - 2]. Такие материалы под действием температурного поля усаживаются, возвращаясь к своим первоначальным размерам, обеспечивая деталям со сложным профилем герметичную защиту. Выбор оптимального температурного воздействия для обеспечения термоусадки изделий является актуальной задачей.
Целью данной работы является разработка метода, позволяющего вывести уравнение для расчета коэффициент усадки термо-усаживающейся трубки в различных температурных режимах.
Разобьем всю структуру полимерной нити на множество сегментов макромолекул [3
- 7] Будем называть такие сегменты кластерами, которые могут находиться в двух основных энергетических состояниях. Состояние 1 соответствует нити перед ее вытяжкой при температуре, близкой к температуре плавления, и является наиболее энергетически выгодным, состояние 2 - состояние, которое занимает структура после охлаждения. На рисунке 1 изображена зависимость внутренней энергии и кластера от его линейных размеров
Рисунок 1. Зависимость внутренней энергии U кластера от его линейных размеров
26
НИИТТС
dN з = V13 dtN 1 - V зldtN з; (9)
dN 1 = IVзldtNз - IVlзdtN 1; (10)
Проведем исследование в обратимой области деформации. Найдем вероятность перехода кластера Ж1з в единицу времени в результате тепловых флуктуаций из состояния 1 в состояние з. Вероятность данного перехода
иБ
пропорциональна е кТ . Обозначим коэффициент пропорциональности за у0 и примем, что
иБ = и12 и ПБ - из = из2. С учетом этого вероятность Ж1з равна:
ип
V =уо в*. (1)
Аналогично, вероятность обратного перехода Жз1 можно представить как:
из2
Г31 = гов (2)
Перейдем к приведенным значениям энергий, т.е. к энергиям, отнесенным к абсолютной температуре. Приведенные величины будем обозначать надстрочным символом -звездочка. Тогда:
V =у0 в4; (з)
=у0в-и:. (4)
В условиях, когда система находится в деформированном состоянии, кластеры обладают еще и упругой энергией. Так как упругая
2
энергия пропорциональна а ,то механическое напряжение играет активирующую роль. Оно понижает барьер, т.е. способствует переходу 1
- з и повышает наоборот барьер в направлении з - 2 на величину упругой энергии. Все переходы в напряженном состоянии будем обозначать символом ”~”. С учетом вышесказанного:
и*2 = и; -х*X2; (5)
и з2 = и 3*2 + /' X2, (6)
, 1
где у------и является структурночувствитель-
Т
ным коэффициентом, зависящим от рода нити.
В напряженно-деформированном состоянии вероятности переходов 1 - з и з - 1:
^1з =пе-и”; (7)
{V з1 = Г0 в^2. (8)
Составим уравнение, отражающее закон сохранения общего числа кластеров с учетом двух возможных устойчивых состояний;
d (N з + N1) = dN^^ = 0, (11)
где: dNз - изменение числа кластеров в состоянии з за единицу времени; dN 1 - изменение числа кластеров в состоянии 1 за единицу времени; N0 - начальное число кластеров.
Рассмотрим образец синтетической нити с первоначальной длиной Ь0. В состоянии 1
находится число кластеров равной с длиной
а, в состоянии з - Nз с длиной Ь .
Назовем Ь - а = 5 - квантом деформации.
Введем число кластеров Ni с длиной с. N4 - число кластеров, которые дают остаточную деформацию.
Так как время жизни макромолекул на вершине барьера в состоянии 2 мало, их числом можно пренебречь. Тогда первоначальная длина образца равна:
Ь0 = N1a + NзЬ + Ni с . (12)
При деформировании (растяжении) образца число кластеров в состоянии один - N1, а в состоянии з - Nз. Число кластеров в состоянии N4 остается неизменным. Длина образца Ь при приложении нагрузки:
Ь = N 1а + N зЬ + N4 с + Ь0 х , (1з)
где Ь0х - упругая часть.
Найдем абсолютное удлинение М = - ^)а + (Nз - Nз)Ь + Ь0х (14)
Вычислим относительную деформацию, поделивуравнение (14) на Ь0:
£= х + а(N1 -N1) + ЬN -N3) (15)
Ь 1 Ь
Ь0 Ь0
а ~ Ь ~
Часть —(И 1 -N1) + — (Nз -N3) пред-
Ь 1 Ь
Ь0 Ь0
ставляет собой конформационную часть деформации £конф .
Продифференцируем уравнение (15) по времени, воспользуемся формулами (9) и (10) и
учтем, что (Vз1 Nз - V 1зN1) и (V 1зN1 - Vз1 Nз) отличаются только знаком. В итоге получим следующее выражение:
Ь - а і
— (є- х) =-IV 13 N1 - V 31 N3 | = ...
йі Ь0
... = — IV 13N - V31N3 |.
А1
Так как
(16)
Ж1 + N 3 = N 0;
(17)
[(IV13 + IV 31) ТУ 1 - IV 31N 0 ].
(18)
В стационарном состоянии левая часть уравнении (18) обращается в ноль. Состояния можно считать стационарными, так как вследствие малой скорости деформации переходы в кластерах успевают отрелаксировать быстрее, чем происходит растяжение. Тогда:
(V 1з + IV з1) N1 - V з^ = 0;
_ _ _ _ 0 (19)
(V 1з + V з1) N1 = V з1 N0.
Отсюда находим число кластеров, находящихся в состоянии 1, для образца в напряженно-деформированном состоянии.
Г ш \
N. = N
(20)
У V 1з + V з1 у
Аналогично, число кластеров, находящихся в состоянии 1 и з, для недеформирован-ного образца равно
Ґ
N = N
N = N„
Л
(21)
(22)
Перепишем уравнение (15) , учитывая уравнение (17):
SN1 ЬЫ 1
----=----1 -(е-х)---(аК + ЬN3) .
А А А
(23)
Подставим в уравнение (23) выражения (21) и (22):
(24)
N
Обозначим —- = т0, где т0 - число Ь0
кластеров на единицу длины, т.е. линейная
плотность кластеров. С учетом т0 уравнение
(24) привет вид:
8М1 т
---- = -(£- х) +-
N з = N 0 - N1,
то уравнение (16) преобразуется к следующему виду:
d 8 |— ~~
— (е-х) = —I V 1зN1 -VзМа + Vз1 N1 1 = ... dt Ь^ J
... = -(е - х) +
(25)
- +1
Подставим в уравнение (25) выражения для Ж1з и Ж31 из уравнений (з) и (4) и учтем, что
дN 1 от 5
= -(Є - X) Н--------;--
и (26)
Ь е-из+1
Подставим в уравнение (18) выражение
3
для — из уравнения (25):
йі
-(£- X) =
-(£-Х )'
от 8
-и,
Є 3 + 1J
(ж 13 + V 31)
5 ~
V13 + IV 31 I----------V 31^
А
(27)
N
Учитывая, что —1 = т0 и после ряда Ь0
преобразований уравнение (27) перепишем в ином виде:
й(є - х)
йі
+ (IVи + IV3,)(є- х) = Ш08
IV13 + V 3:
е ’ +1
[V13
- IV31 е~
е 3 + 1 Обозначим I =
(28>
и рассмотрим
V 1з - V з1 е
отдельно, учитывая уравнения (5) - (8), а также равенство из2 + из* = и*:
I =
= ^0(е
-и„ + у х -П -у х -и3
- е 3
... = У0 е-и“(е? х - е‘Г х ).
Ьв;[
Учитывая, что по определению shZ = гиперболический синус, то
) =... (29)
2
shy x =■
e - e
2
(30)
Тогда уравнение (29) можно переписать
в виде:
I = 2v0e 12sh(y x ) Учитывая выше сказанное: d (є - x)
(31)
dt
-U„ , у x U, -у x , ■.
+ V0 e 12 (e + e 3 )(£ - x) = ...
2m JS e"U + 1
-v e 12 sh(y x2).
( 32)
Пусть 0^ = —— - время релаксации, т.е.
среднее время прохождения системы через потенциальный барьер в состоянии при отсутствии нагрузки. С учетом того, что №^2 = г0 е~и" , получаем выражение для времени релаксации:
1
© =— eU
(33)
зависящее от температуры t
© '
р
Пусть выражение
2 т, 8
(34)
ет описать различные режимы деформирования, в частности, процесс термической усадки. Условия термической усадки:
ст = const = 0; T = f (t).
При условии х = 0 уравнение (35) переписывается в виде: ds
где
dr
eU* = A
- + (1 + A) є = 0,
(36)
Введем новое безразмерное время прохождения через барьер, которое будет выражаться в долях от времени релаксации © , и
(Т) - безразмерный параметр, являющийся функцией температуры.
Произведение то8 показывает максимальную деформацию, вызванную конформа-ционными переходами.
Параметр % - это максимальная де-
формация, вызванная конформационными переходами при определенной температуре. С учетом % уравнение (з2) можно переписать в виде:
— (е - х) + (е - х)(егх + еиз гх ) = xsh(y х2). (з5) йт
Уравнение (з5) работает при гипотезе о независимости кластеров друг от друга, т.е. при отсутствии корреляции. Таким образом, в рамках предложенной модели было получено определяющее уравнение (з5).
Обыкновенное нелинейное дифференциальное уравнение 1- го порядка (з5) позволя-
Уравнение (36) - это дифференциальное уравнение с разделяющимися переменными.
Здесь г = —, где і - время выдерживания об©
разца в печи. При этом выполняется условие: при г , е ^ £нас
Приведем решение уравнения (36): £=£йас[1 - е^(1+">'] (37)
Уравнение (37) с достаточной точностью соответствует экспериментальным данным по усадке мононити из модифицированного полиолефина, комплексной полиэтилентере-фталатной нити и текстильного полотна, состоящего из слабоусаживающейся основы в продольном направлении и указанных мононити и комплексной нити в поперечном направлении в диапазоне температур 115 - 200 °С [8]. Продольная усадка мононити и комплексной нити приводит к усадке полотна в поперечном направлении. В частности, можно привести зависимость деформации от времени для указанных систем при 150 0С. На рисунках 2 - 4 приведены экспериментальные и теоретические зависимости коэффициента продольной деформации £прод,% от времени выдерживания в
печи для мононити, комплексной нити и текстильного полотна соответственно, а также указаны соответствующие дифференциальные уравнения с их решениями.
Предложенная статистическая модель в достаточной степени точности описывает процесс усадки различных полимерных систем при воздействии температуры. Данная модель применима не только к “простым” моделям - мононити, комплексной нити, но и к более сложным конструкциям, в частности, к текстильному полотну. Возможность применения статистической теории обусловлена хорошей согла-сованностью с экспериментом в широком диапазоне температур 115 - 200 0С.
0
прод ’
eпpoд,
Рисунок 2. Зависимость коэффициента продольной деформации мононити от времени выдерживания в печи при температуре 150 0С: 1 - эксперимент; 2 -теория
— + 0,73s = 0; є = S7,5(1 - e“°,73t) dt
Рисунок з. Зависимость коэффициента продольной деформации комплексной нити от времени выдерживания в печи при температуре 150 0С: 1
- эксперимент; 2 - теория
ds • 0,S0t4
— + 0, Ss = 0; • • 11(1 • e ’ )
dt
Рисунок 4. Зависимость коэффициента поперечной деформации полотна от времени выдерживания в печи при температуре 150 °С: 1 - эксперимент; 2 - теория
d^ С\ С\ —0,76t\
----h 0,76s = 0; s = 50(1 - e ’ )
dt
Литература
1. Белошенко B.A. Эффект памяти формы в полиме-рах//Успехи химии:Обзорный журнал по химии.-2005. -Т.74. №3. -с.285 - 306.
2. Каюмов Р.А, Строганов И.В., Строганов В.Ф., Мухаметшин А. Т. Математические модели поведения полимерного материала с "памятью формы”// Известия КазГАСУ.-2009.- № 2 (12). -с.250 - 256.
3. Гросберг А.Ю., Хохлов А. Р. Полимеры и биополимеры с точки зрения физики. -М.: 2010.
4. Гросберг А.Ю., Хохлов А. Р. Статистическая фи-зикамакромолекул. М.:Наука, 1989.
5. Сандитов Д.С., Бартенев Г.М. Физические свойства неупорядоченных структур.-Новосибирск:
Наука, 1982.
6. Рымкевич П.П., Романова А.А., Сталевич А.М. Кинетическое описание релаксации механического напряжения в синтетических нитях//Изв. Вузов Тех-нол.текст.пром. 1999. №3 -с.56 -61.
7. Сталевич А.М., Макаров В.Г. Вариант прогнозирования процессов деформирования синтетических нитей // Хим.волокна. -2001.№4. - с.67 - 69.
8. Рымкевич О.В., Цобкалло Е.С. Влияние температурных режимов на усадку модифицированной по-лиолефиновой мононити термоусаживающейся трубки.//Изв. Вузов. Технол.лег.пром.// 2012. №2. с.13 - 16.
%
1 Рымкевич Ольга Васильевна - аспирант СПбГУТД, моб.: +7 921 316 92 26, e-mail:olga.rymkevich@ gmail.com;
2 Рымкевич Павел Павлович - кандидат физико-математических наук, доцент, профессор кафедры «Прикладная физика», СПбГУСЭ, ,моб.:+7(911) 224 59 13, e-mail: [email protected];
3 Романова Алла Александровна, кандидат технических наук, доцент, доцент кафедры «Прикладная физика», СПбГУСЭ, ,моб.:+7(9П)2ПЗ426, e-mail:[email protected].