ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 2000, том 42, № 1, с. 71-83
ТЕОРИЯ
УДК 541.64:532.73
НЕЛИНЕЙНАЯ ТЕОРИЯ ПРОЦЕССОВ НАБУХАНИЯ ЭЛАСТОМЕРОВ В НИЗКОМОЛЕКУЛЯРНЫХ ЖИДКОСТЯХ1
©2000 г. Е. Я. Денисюк, В. В. Терешатов
Институт технической химии Уральского отделения Российской академии наук
614000 Пермь, ул. Ленина, 13
Поступила в редакцию 30.12.98 г. Принята в печать 26.08.99 г.
Предложена трехмерная теория процессов набухания эластомеров в низкомолекулярных жидкостях, которая строится как теория массопереноса в нелинейных упругих материалах. Общая теория применена для исследования кинетики набухания плоского слоя эластомера при больших степенях набухания. Для этого случая получены явные асимптотические выражения, описывающие началь-' ную и конечную стадии процесса. Обнаружена возможность существования качественно различных режимов набухания. Установлена связь кинетики набухания с физическими свойствами материала. Эти результаты позволяют объяснить механизм аномалий кинетики набухания и предложить новые методы исследования диффузионных и термодинамических свойств эластомеров.
Как известно, кинетика набухания многих полимерных материалов в низкомолекулярных жидкостях при больших степенях набухания заметно отличается от законов нормальной сорбции, вытекающих из классической теории диффузии Фика [1]. Эти отклонения в литературе называются аномалиями [1,2]. По современным представлениям одной из основных причин наблюдаемых аномалий является влияние на диффузионный процесс напряженно-деформированного состояния полимерной матрицы, возникающего из-за набухания материала [3-7].
В настоящей работе представлена геометрически и физически нелинейная теория, позволяющая рассматривать набухание сшитых эластомеров как совокупность взаимосвязанных процессов диффузии жидкости и деформирования материала. Ее основой служит ранее предложенная авторами теория механо-диффузионных процессов в упруго деформируемых средах [8].
ТРЕХМЕРНАЯ ТЕОРИЯ НАБУХАНИЯ ЭЛАСТОМЕРОВ В НИЗКОМОЛЕКУЛЯРНЫХ ЖИДКОСТЯХ
Эластомер и растворенную в нем жидкость будем для краткости называть смесью. Выберем некоторую фиксированную конфигурацию смеси, которую назовем отсчетной. В этой конфигура-
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 98-ОЗ-ЗЗЗЗЗа).
ции введем систему материальных координат ф (к - 1,2,3), которую свяжем с полимерной матрицей. Пространственное положение частиц матри-цы в отсчетной конфигурации определим радиусом-вектором
г = г(д\ V),
(1)
а их движение в пространстве будем описывать радиусом-вектором
= Щд\д2,д3,1)
(2)
Этой функцией определяется актуальная конфигурация матрицы.
В отсчетной и актуальной конфигурациях определим векторные базисы гк = дг/ддк, К* = ЭИ/Э^ и операторы градиентов
V = г*Э/Э/, V = И*Э/Ъц
Деформации материала будем характеризовать с помощью тензоров второго ранга .
о = УК, Е = <2Г <2,
которые называются градиентом места и мерой деформации Фингера соответственно [9].
Тензор Г симметричен и имеет три независимых инварианта (у = 1, 2, 3), которые выража-
ются через главные относительные удлинения материала Хк '
11 — А,1 + Х^ + Я-2,
(3)
Инвариант градиента места
] = ск* <2 = Х,12Х3 (4)
характеризует объемную деформацию материала.
Определим новую меру деформации и скалярные параметры
£ = ¥Гт, Г, = 1х/]т, Г2 = /2//'3,
Шл
(5)
где со = -отношение объемов У0 и ^"сухого" и набухшего полимера в ненагруженном состоянии; (Оо - значение со в отсчетной конфигурации.
Из определений (5) следует, что тензор £ и параметры Г,, Г2 нечувствительны к объемным деформациям смеси и характеризуют искажения ее формы. Сжимаемость смеси под действием механической нагрузки характеризует параметр Г3.
Набухание эластомеров в низкомолекулярных жидкостях - это частный случай механо-диффу-зионных процессов, которые в упруго деформируемых материалах (в изотермических условиях) описываются следующей системой уравнений [8]:
йпх/(¡1 + п,У ■ у2 - -V • ^ с1п2/с1г + п2 V • \2 = О V • <г = О j = -п,К • У|Х
(6),
(7)
(8) (9)
Здесь па и уа - мольные концентрации и скорости движения компонентов смеси соответственно (значение индекса а = 1. соответствует жидкости, а значение а = 2 - полимерной матрице); а - тензор упругих напряжений Коши (тензор "истинных" напряжений)^ = - у2) - плотность диффузионного потока жидкости относительно полимерной матрицы; |х - химический потенциал жидкости; К - неотрицательный, симметричный тензор второго ранга, характеризующий проницаемость полимерной сетки; материальная производная ¿/ей = Э/Эг + уг • V описывает изменение зависящих от времени величин в системе отсчета, связанной с деформируемой полимерной матрицей.
К уравнениям (6)—(9) необходимо добавить физические (определяющие) соотношения для тензора упругих напряжений, химического потенциала и тензора проницаемости. Для.изотропного упругого материала они имеют вид [8]'
<т = ^-^Г.Е^ + Ц^/.^Е^-рЕ (10)
^ = 2(Э//ЭГ, + Г,Э//ЭГ2), ^2 = -2(Э//ЭГ2)
р = -г3э//аг3 (12)
II = Э//Эл, + У1Р, V, = ю^Эю/Э«! (13)
К^коЕ + к^ + к^2 (14)
к, = МГ1)Г2,Гз,п1), 5 = 0,1,2 (15)
/ = /(Г\, г2, Г3, пъ п2), (16)
где Ё - единичный тензор; р - давление; Ух - парциальный мольный объем жидкости; / - плотность свободной энергии смеси. Соотношения (14)-<<16). являются наиболее общими представлениями тензорной и скалярной функций, характеризующих физические свойства изотропной упругой смеси.
Определяющие соотношения (10)-( 16) записаны в терминах меры деформации £ и параметров Гг Такая параметризация делает их инвариантными относительно выбора отсчетной конфигурации. Благодаря этому заметно упрощает описание термодинамических свойств смеси в сложно-напряженном состоянии.
Если смесь считать несжимаемой, то следует положить
Г,= 1
(17)
При этом все уравнения и определяющие соотношения, за исключением выражения (12), остаются в силе, но в них пропадает зависимость от Г3. Давление в этом случае является свободным параметром, который определяется из решения полной системы уравнений с условием несжимаемости (17).
Уравнения (6)-(9) описывают процесс набухания в актуальной конфигурации, которая, как правило, неизвестна и подлежит определению. В некоторых случаях математическое описание процесса заметно упрощается, если все уравнения сформулировать в отсчетной конфигурации. Это достигается путем применения к исходной систе-
ме уравнений известных преобразований [9]. Приведем результат этих преобразований [8]
Эл°/Э/ = -V •
Эп^/Эг = О, V • <ти = 0, 1° = -п^К" • Уц,
ОжгО
где
на = Зпа, а = 1,2
<т° = /(СГ1)Т-о
(18)
(19)
(20)
? = т~1)Т-За, К0 = (О-'/'КО"1 (21)
Здесь полагается, что все величины, входящие в уравнения, являются функциями материальных координат, а частные производные по времени берутся при фиксированных ф, поэтому они эквивалентны материальной производной. Величина Па (а =1, 2) численно равна количеству молей й-го компонента, приходящемуся на единичный объем смеси в отсчетной конфигурации. Тензор <т° называется тензором напряжений Пиола [9], а - диффузионным потоком Пиола. Они определяют механические напряжения и диффузионный поток, отнесенные к единичной площадке отсчетной конфигурации.
Из второго уравнения (18) и формулы (19) следует, что
п2 = п2/У,
(22)
где п2 не зависит от времени. Эта величина не зависит и от координат, так как отсчетная конфигурация является неискаженной'.
Приложение общей теории к описанию процессов набухания эластомеров предполагает задание конкретных определяющих соотношений для тензора упругих напряжений, химического потенциала жидкости и тензора проницаемости. Воспользуемся теорией полимерных сеток Флори [10] и запишем ПЛОТНОСТЬ свободной энергии набухшего в жидкое ти эластомера в инвариантном относительно выбора отсчетной конфигурации виде [8]
/ ЯТ^Пу 1пф, + п2ЬЩ2 + %ф2"1 + ^Ф^Г,/^ (23)
Здесь ф„ ф2 - объемные доли жидкости и полимера соответственно; % - параметр взаимодействия Флори-Хаггинса; Ь - безразмерный параметр, связанный с функциональностью полимерной сетки. Последний член описывает упругую реакцию полимерной сетки.
В теории Флори полагается, что смесь несжимаема, а объемы компонентов смеси аддитивны, поэтому (0 = ф2. С учетом этого из условия несжимаемости (17) следует, что 7= ф2/ф2, где ф2 - объемная доля полимера в отсчетной конфигурации. Если отсчетной конфигурацией служит состоя^
ние "сухого" эластомера, то ф2 = 1 и
/=1/ф2 (24)
Применив к соотношению (23) формулы (10), (11) и (13), получим следующие выражения для тензора напряжений и химического потенциала:
а = ДГ^'фГ^-^Г.Е^-рЕ (25) |£ = *г(ьф1 + ф2 + хф2 + ^Ф210Г1] + рУ1, (26)
где X - У21Уу, V,, У2 - мольные объемы жидкости и субцепей полимерной сетки соответственно. В формуле (26) опущен малый член -Ь2~1 ф2, который в дальнейших выкладках не играет существенной роли.
Выражение для тензора напряжений Пиола следует из соотношений (20) и (25). Его физические компоненты в главных осях имеют вид
а°к .= Я7^'(ф°) -1/,А*-')-к = 1,2,3
(27)
Что касается тензора проницаемости, то его общее выражение дается формулой (14), причем скалярные функции (15) в силу несжимаемости смеси не зависят от Г3. Чтобы установить конкретный вид этих функций, необходимо обратиться к эксперименту.
ПОСТАНОВКА ЗАДАЧИ О НАБУХАНИИ БЕСКОНЕЧНОГО ПЛОСКОГО СЛОЯ
Рассмотрим две задачи о набухании бесконечного плоского слоя эластомера в хорошем растворителе. В первой из них будем полагать, что слой эластомера свободен, т.е. не подвергается внешним механическим нагрузкам. Во второй задаче будем считать, что эластомер скреплен с не-деформируемой подложкой. Пусть в начальный момент времени эластомер не содержит жидкости и находится в неискаженном состоянии. Это состояние выберем в качестве отсчетной конфигурации, а декартовы координаты, связанные с полимерной матрицей в этом состоянии, - в качестве материальной системы координат. Ось х на-
правим поперек слоя, остальные оси - вдоль слоя. Обе задачи рассмотрим в одномерной постановке, т.е. будем считать, что все величины, характеризующие состояние эластомера, зависят только от одной поперечной координаты х.
В случае свободного набухания слоя начало координат поместим в центр слоя и движение полимерной матрицы представим в виде
X = Х(х, 0, У = \^)у, г = (28)
где (х, у, г) - материальные, а (X, У, 7) - пространственные координаты, определяющие актуальную конфигурацию полимерной матрицы. Тогда относительные удлинения эластомера
(29)
Хх = % = дХ/дх, Ху = Xz = дУ/ду = dZ/dz = v(f)
Отсюда и из формулы (4) имеем
J = Xv2 (30)
Уравнения и краевые условия, описывающие свободное набухание плоского слоя в отсчетной конфигурации, имеют вид
(31)
Ъгц/Ы = -д]°х/дх, Л = -пХхЪ\1/дх, да°х/дх = 6
л?(х,0) = 0,' Х(л:,0) = х (32)
Эц(0, г)/Ьх = 0, Х(0,0 = 0 (33)
ц(й, 0 = 0, <т"(А, /) = 0 (34)
<о> = '<<# =0 (35)
здесь А - полутолщина слоя, а угловыми скобками обозначена операция интегрирования по координате х:
<...> = й ..
dx
Начальные условия (32) очевидны. Условия (33) означают отсутствие диффузионного потока и перемещения полимерной матрицы вдоль оси х в центре слоя. Условия (34) выражают наличие термодинамического равновесия на границе эластомер - жидкость и отсутствие механического на-гружения эластомера в поперечном направлении, а выражение (35) есть условие отсутствия внеш-
него нагружения слоя в продольном направлении. К этим уравнениям необходимо также добавить условие несжимаемости смеси (24). '
Из уравнения механического равновесия (31) и второго условия (34) следует, что а°х = 0. Данное уравнение и соотношение (27) дают выражение для давления, подстановка которого в формулы для химического потенциала (26) и продольных напряжений (27) с учетом (30) приводит к
О0 = o°z = RTV2\\-J2/V5)
(36)
ц = ■ДГ(ЬФ1+Ф2+хФ2+;Г,Ф;1У-4) (37)
Химический потенциал, как следует из формулы (37), зависит от двух переменных - концентрации жидкости и относительного продольного удлинения слоя. Это справедливо и для коэффициента проницаемости. Продольное удлинение слоя V не зависит от координаты х, поэтому диффузионный поток можно представить в виде
fx = -Dbn\/dx,
(38)
где коэффициент диффузии О определяется соотношением
Щп1 V) = пМ"?, v)Э^l(n?, V)/Эп? (39)
и является функцией концентрации жидкости и относительного продольного удлинения слоя. Уравнение переноса жидкости (31) с учетом формулы (39) запишем в виде
эГ = З^"" у)йГ
(40)
Подставив уравнение (36) в выражение (35), для продольного удлинения слоя получим
V = {J\x,t))m
(41)
Краевые условия с учетом формул (37) и (38) приобретают следующий вид:
дпЧ/дх = 0 при х = 0 (42)
1пф, + ф2 + хф2 + Тх§2\~4 = 0 при х = А (43)
Таким образом, исходная задача о набухании бесконечного плоского слоя эластомера сведена к краевой задаче (40)-(43) с начальным условием (32)- Ее решение позволяет полностью описать процесс набухания свободного слоя, т.е. определить текущее распределение жидкости в слое и с помощью соотношений (28)-(30), (36) рассчитать
напряженно-деформированное состояние материала.
Уравнение (40) не является обычным уравнением диффузии. Коэффициент диффузии в нем зависит от концентрации жидкости и продольной деформации слоя. Граничное условие (43) также не имеет аналогов среди тех граничных условий, которые используются при традиционном описании процессов сорбции. Интересно отметить, что это условие задает в системе положительную обратную связь, которая в процессе набухания материала порождает нестационарный граничный режим. Нелинейные распределенные системы с положительной обратной связью, как известно, называются активными средами и характеризуются сложностью и многовариантностью поведения [11]. При свободном набухании эластомеры в определенном смысле ведут себя как активные среды. Наиболее ярко это проявляется при больших степенях набухания; рассмотрим данный случай подробнее.
Обозначим объемную долю полимера в равновесно набухшем состоянии эластомера через е. Она определяется из уравнения
1п( 1 - е) + е + хе2 + Т1гш = 0, (44)
вытекающего из уравнения (43), если в нем положить V = £"ш. При больших степенях набухания величина е мала и может рассматриваться как малый параметр. Разложив 1п(1 - е) в ряд по степеням е до второго порядка, из формулы (44) получим известную асимптотическую формулу
е = [2(0.5-х)Г3'5
Действуя аналогичным образом, из уравнения (43) получим выражение для объемной степени набухания слоя на границе эластомера с жидкостью. С учетом выражения (24) и последнего соотношения запишем его в виде
У = г-\гтм)6", (45)
где (1 = 2/9. Отметим, что это приближенное соотношение определяет деформационную зависимость равновесной степени набухания эластомера в жидкости в условиях однородного двухосного растяжения и сжатия при больших степенях набухания материала.
ч
Подставив выражение (41) в формулу (45), получим условие на границе эластомер-жидкость
У(А, /) = е2"~1 и2(х, о/ (46)
Окончательную формулировку задачи представим в виде
и, = (к(и,1)их)х, (0,1), i>0 (47) m(jc, 0) = 0, ux(l,t) = 0 (48)
(1 - е)м(0, о + е = < [(1 - еМх, 0 + е]У (49)
Здесь осуществлен переход к безразмерным переменным. Единицами измерения расстояния и времени служат соответственно величины h и h2/D0, где D0 - значение коэффициента диффузии для предельно набухшего состояния эластомера. Для удобства здесь также произведено преобразование пространственной переменной / —- 1 - jc. Угловыми скобками в выражении (49) обозначено интегрирование по jc в пределах от 0 до 1. Функция и(х, 0 принимает значения в интервале (0, 1) и является безразмерной концентрацией диффундирующей жидкости. С концентрацией жидкости и локальным объемным расширением материала она связана соотношениями
п\ = v^V-iM*,*),
(50)
J(x,t) = е 1 [(1-£)«(*, 0+ е]
Величина / = el/3v представляет собой нормированное на единицу продольное удлинение слоя. В силу выражений (41) и (50) .
/(/) = <[(l-e)«(jc,i) + e]2>1/6 (51)
Безразмерный коэффициент диффузии определяется формулой
к(и, I) = D(M, l)/D0
Здесь аргументы правой части равенства выражены через ми/.
Продольные напряжения слоя (36) выразим через безразмерные напряжения q(x, t) по формуле
а0у = = RTV~2q{x, t), (52)
где в силу выражений (36) и (50)
(53)
= e~ml(t){ 1 - [(1 - е)и(х, t) + e]2/f(t)} Рассмотрим две функции
= <u(x,t)>, g2(t) = (u2(x,t)), (54)
которые являются интегральными характеристиками кинетики набухания свободного плоского слоя и могут быть определены экспериментально. Первая функция характеризует относительное количество жидкости, поглощаемое эласто-
и(х, О
X
Рис. 1. Распределение диффундирующей жидкости (а) и продольных напряжений (б) в свободном плоском слое при его набухании с постоянным коэффициентом диффузии к(и, /)= 1 ие = 0.1. г=0.05 (7), 0.2 (2), 0.4 (5), 0.6 (4), 1.0 (5) и 1.8 (6).
мером за время г, вторая функция, как следует из формулы (51), связана с продольной деформацией слоя. Для сильно набухающих эластомеров
82«) = /6(0 (55)
Результаты численного решения задачи (47)-(49) при постоянном коэффициенте диффузии представлены на рис. 1. Здесь показана эволюция концентрационного поля диффундирующей жидкости и продольных напряжений. Видно, что граничная концентрация жидкости в процессе набухания слоя монотонно растет.
При набухании эластомера, скрепленного с не-деформируемой подложкой, продольная деформация отсутствует. Постановка соответствующей задачи практически полностью повторяет вывод
предыдущей и приводит к уравнению (40) и граничным условиям (42), (43), в которых V = 1. Все характеристики скрепленного слоя, включая химический потенциал, оказываются зависящими только от концентрации жидкости. Отсюда следует, что граничная концентрация жидкости не зависит от времени и соответствует равновесно набухшему состоянию скрепленного эластомера. Эта величина определяется из уравнения (43), его нужно решить при V = 1. Коэффициент диффузии жидкости в скрепленном слое определяется формулой (39), в которой следует положить V = 1, поэтому количественно он отличается от коэффициента диффузии жидкости в свободном слое.
Перейдя к безразмерным параметрам по тем же правилам, что и в случае свободного слоя, краевую задачу, описывающую процесс набухания скрепленного слоя, запишем в виде
ы, = (к(и)их)х, и(х,0) - 0,
«^(1,0 = 0, и(0,0.= 1
Здесь единицами измерения расстояния и времени служат соответственнЪ величины й и к2Юь где Л - толщина слоя, £>, - значение коэффициента диффузии для предельно набухшего состояния скрепленного слоя.
Полученная краевая задача и ее решение хорошо известны. Они описывают кинетику сорбции нормального типа [1,2]. Продольные напряжения определяются формулой (52), в которой
Я{х,ь) = 1- - еТа[( 1 - е^иСх, О +
(в! - объемная доля полимера в предельно набухшем состоянии скрепленного слоя).
ДИФФУЗИОННАЯ КИНЕТИКА НАБУХАНЙЯ СВОБОДНОГО ПЛОСКОГО СЛОЯ
Рассмотрим две стадии набухания свободного слоя - начальную и конечную. На начальной стадии процессы у противоположных границ слоя конечной толщины незначительно влияют друг на друга, поэтому диффузию жидкости в слой можно рассматривать как диффузию в полупространство.
В самом начале набухания количество поглощенной слоем жидкости мало. Это позволяет в формула (49), (51) положить и(х, 0 « О. В результате имеем
м(0,/) = Уо = (ем-е)/(1-е), /(/) - еш,
и уравнение (47) приобретает вид обычного параболического уравнения, описывающего диффузию на полупрямой с постоянной граничной концентрацией \|/0. Оно имеет автомодельное решение вида
и(х, 0 = у060с/7О- Функция 9(£) удовлетворяет уравнению (к(в, е1/3)0')' + £в'/2 = 0 и краевым условиям
6(0) = 1, 0(+°°) = 0 (56)
Отсюда следуют выражения для интегральных характеристик процесса:
.1/2
gl(t) = щ Mxt
Ш = \flM2tm,
где
Мр = р =1,2
(57)
Полученные соотношения характеризуют асимптотические свойства процесса набухания при ¿ —0.
По мере поглощения жидкости продольные деформации слоя растут. Это вызывает повышение граничной концентрации жидкости. Для сильно набухающих материалов при достаточно больших временах в выражении (49) можно пренебречь всеми членами, содержащими е в качестве множителя, в результате чего оно приобретает вид
м(0, t) = (и\х, t))"
(58)
Здесь угловыми скобками обозначено интегрирование по х в пределах от 0 до +«>.
Анализ краевой задачи для уравнения (47) на полупрямой с граничным условием (58) при произвольной зависимости к(и, I) представляет собой самостоятельную математическую проблему, выходящую за рамки настоящей работы. По этой причине здесь мы рассмотрим модельную задачу с коэффициентом диффузии вида
к(и, I) = и, s > 0
(59)
Анализ данной упрощенной задачи позволяет качественно понять многие закономерности диффузионной кинетики свободного набухания эластомеров.
Уравнение (47) на полупрямой с коэффициентом диффузии (59) имеет богатый набор автомодельных решений вида
и(х, I) = \|/(0в(*/ф(0),
где функция 0(^) удовлетворяет условиям (56) и определяет профиль диффузионной волны,
описывает граничный режим, а <р(/) - глубину проникновения диффузионной волны [12]. Среди этих решений есть такие, которые удовлетворяют граничному условию (58). Именно они описывают диффузионную кинетику набухания материала.
При набухании свободного плоского слоя эластомера могут существовать качественно различные автомодельные диффузионные режимы: степенной, экспоненциальный и режим с обострением. Степенной режим имеет место при достаточно малых значениях 5. В данном случае количество поглощаемой жидкости, граничная концентрация, глубина проникновения диффузионной волны и продольная деформация слоя являются степенными функциями времени. Если параметр 5 достигает критического значения
= г/й-4, (60)
то возникает режим, при котором процесс набухания протекает по экспоненциальному закону. И наконец, если 5 > то набухание развивается в
режиме с обострением. Ниже приведены решения, описывающие эти режимы.
I. Степенной режим (ж < я*):
,2т.
М(*,0 = % = x/[M?s(t + t0)n]
gl(t) = AM^-^+fo)«; g2(t) = M22r(t + t0y
m =
4 =
S*и s
l/d-l
n =
r =
l/d-2
Sic — s
1
(61)
(62)
(63)
S*. — s
d(s*-s)
П. Экспоненциальный режим (s =
u(x,t) = ехр[лф-*0)]е(4), £ = xM2exp[(-sM22(t-10))/2]
gi(t) = (AVM2)exp[(i/2+l)A^(i-i0)],
g2(t) = exp[(s/2 + 2 )M2(i - r0>] III. Режим с обострением (s > st): u(x,t) = Mf(/0-i)m6(4),
% = x/[M?s(t0-t)n] gl(t) = M.M^-'Co-')*,
g2(t) = м22Ъо-*У,
где степенные показатели определяются формулами (63), но поскольку s > s*, величины т, п, q,
г< 0.
Для степенного режима и режима с обострением функция 8(4) определяется из уравнения
(е50')'+ |п|^в'-|т|е = о
Для экспоненциального режима в этом уравнении следует положить п = sJ2 и т = 1. Константы
Л/, и М2 определяются по формуле (57), а константы t0 можно оценить по порядку величины из условия и(0, f0) ~ Yo - Для экспоненциального
режима ц = -2dM~2 Ine, а для степенного и режима с обострением t0 - М^' г2Щт. Как следует из последней формулы, при степенном режиме с увеличением равновесной степени набухания материала Y0 —► 0, поэтому , при больших степенях набухания можно считать, что t0 ~ 0.
Интересно отметить, что при т = l/s (соотношение выполняется при s = l/d- 2) решение (61) выражается через элементарные функции и описывает диффузионную волну, распространяющуюся с постоянной скоростью, равной d:
u(x,t) = (l-2d)lls[d(t + t0)-x)]l's
Решение, описывающее режим с обострением, за конечное время обращается в бесконечность. Краевые эволюционные задачи, допускающие такие решения, являются глобально неразрешимыми по времени и в последнее время находят применение для моделирования различных физи-ко-химических процессов [12]. Очевидно, что при набухании слоя конечной толщины все решения ограничены, а представленные здесь автомодельные, решения описывают асимптотические свойства диффузионных режимов на начальной, достаточно развитой стадии набухания.
На конечной стадии набухания и(х, t) —► 1. Приближенное решение задачи (47)-(49) в данном случае можно получить, линеаризовав ее в окрестности точки равновесия и = 1. Для этого введем новую переменную vOc, 0=1— и(х, t). После преобразований имеем
= vx(l,t) = 0, v(0,f) = 2d(v(x,t))
Применив метод разделения переменных, получим решение этой задачи в виде
о©
v(x,t) = £а*ехр(-аМсо8[(а*(1-х)], (64)
к- 1
где ак - положительные корни характеристического уравнения
ак = 2dщ(<ak) (65)
Ограничившись первым членом ряда (64), для конечной стадии набухания запишем уравнение кинетической кривой
.#,(*) = 1 -(й1/а[)ехр(-а1*)8та, (66)
Как показано выше, по теорий Флори й = 2/9. В этом случае первый корень уравнения (65) а! = = 1.2220. Сравнив его с соответствующим значением в модели нормальной сорбции (а! = тс/2 ~ = 1/571), нетрудно видеть, что конечная стадия свободного набухания эластомера оказывается более затянутой во времени по сравнению с нормальной сорбцией.
Полученные выше выражения позволяют описать форму кинетических кривых ¿^(О в целом. В частности, из соотношений (62)-(63) следует, что при q< 1 кинетические кривые в координатах (г, #,) являются выпуклыми вверх и имеют вид, характерный для псевдонормальной сорбции. Это имеет место при .$ < ¿*/2 - 1. Если 5 = л*/2 - 1, то # = 1, что соответствует линейному режиму поглощения жидкости. При 5 > 5^/2 - 1 начальный участок кинетической кривой становится выпуклым вниз, а кривая в целом приобретает Б-образую форму. Отметим, что в координатах (г|/2, #,) при 5 > 0 все кинетические кривые имеют Б-образный вид.
На рис. 2 представлены результаты численного решения задачи (47Н49) с коэффициентом диффузии вида (59). Кинетические кривые набухания при различных значениях ^ изображены на рис. 3. Они имеют разную форму и описывают различные типы аномалий.
Как показано выше, коэффициент диффузии может зависеть от продольной деформации слоя, поэтому возникает вопрос о влиянии этой зависимости на диффузионную кинетику набухания эластомера. Оказывается, что наличие деформационной зависимости коэффициента диффузии вида к(и, /) = и?!?? ~ и*(и(х, гУ не порождает новых диффузионных режимов. Все три типа автомо-
0.5
(В)
—Л—
\25 \
20 | 1
0.5
Рис. 2. Диффузионная кинетика набухания свободного плоского слоя с коэффициентом диффузии к(и) = и5 и е = 0.1: а, б - степенной режим набухания при 5 = 1 и 2.5 соответственно; в - экспоненциальный режим (я = 5); г - режим с обострением = 5.5). Числа у кривых соответствуют моментам времени /.
дельных решении остаются в силе, но в этом случае критическое значение
** = (2-р)/</-4
Решения, описывающие диффузионную кинетику набухания плоского слоя, могут служить для экспериментального исследования диффузионных и термодинамических свойств эластомеров. В частности, по формуле
5 = 2Л/-4-(1Л/~ 1)/<7,
(67)
следующей из выражений (60) и (63), можно оценить концентрационную зависимость коэффициента диффузии жидкости в эластомере. Параметр д в силу (62) определяется по начальному участку кинетической кривой набухания. Оценив параметр г из деформационной кривой набухания 1(1) с помощью соотношения (55), по формуле
¿/=1 -ф, (68)
вытекающей из выражения (63), можно найти параметр й, который характеризует определяемую
соотношением (45) деформационную зависимость равновесной степени набухания эластомера в условиях однородного двухосного растяжения. Стоит отметить, что непосредственное определение й в опытах по равновесному набуханию деформированных образцов представляет* собой трудоемкую экспериментальную задачу.
Значение коэффициента диффузии в равновесно набухшем состоянии эластомера можно определить по конечному участку кинетической кривой набухания с помощью уравнения (66), которое р размерных переменных имеет вид
*,(*) = 1-Сехр(-<х?£у//г2),
(69)
где а, определяется из уравнения (65). При Л = 2/9 а, = 1.2220.
Подчеркнем, что все эти соотношения справедливы только цри достаточно больших степенях набухания эластомеров.
$1
Рис. 3. Кинетические кривые набухания свободного плоского слоя при различных значениях показателя концентрационной зависимости коэффициента диффузии: 5 = 0 (7), 1 (2), 1.5 (5), 2.5 (4), 5 (5), 5.5.(6). 1-4 - степенной режим; 5 - экспоненциальный режим, 6 — режим с обострением. е = 0.1.
РЕЗУЛЬТАТЫ ЭКСПЕРИМЕНТАЛЬНОГО ИССЛЕДОВАНИЯ КИНЕТИКИ НАБУХАНИЯ ЭЛАСТОМЕРОВ
Нелинейная теория набухания, как показано выше, предсказывает ряд эффектов, которые должны проявляться при больших степенях набухания эластомеров и могут быть обнаружены экспериментально.
В качестве объектов исследования использовали сшитые аморфные полибутадиенуретаны (ПБУ) с различной плотностью химической сетки: 0.3 кмоль/м3 (ПБУ-1); 0.05 кмоль/м3 (ПБУ-2); 0.2 кмоль/м3 (ПБУ-3); 0.1 кмоль/м3 (ПБУ-4). В качестве сшивающего агента при отверждении форпо-лимера олигобутадиендиола применяли олигоокси-
81
Рис. 4. Кинетические кривые набухания скрепленного (I) и свободных (2,5) образцов эластомеров в ДБС (а) и ДОС (б): 1,2 - ПБУ-1; 3 -ПБУ-2.
пропилентриол - Лапрол 373. Методика изготовления эластомеров изложена в работе [13].
Образцы эластомеров были изготовлены в виде дисков диаметром 35-60 мм и толщиной 2 мм. Часть образцов наклеивали на металлические пластины. Они служили моделью скрепленного слоя. Кинетика набухания эластомеров определялась в следующих низкомолекулярных жидкостях: толуоле, дибутилсебацинате (ДБС) и диок-тилсебацинате (ДОС).
Определить кинетические кривые набухания скрепленных образцов удалось только для умеренно набухающих эластомеров ПБУ-1 и ПБУ-4 в ДБС и ДОС. Скрепленные образцы других эластомеров разрушались в самом начале процесса. На рис. 4 показаны кинетические кривые набухания
НЕЛИНЕЙНАЯ ТЕОРИЯ ПРОЦЕССОВ НАБУХАНИЯ ЭЛАСТОМЕРОВ Таблица 1. Характеристики кинетики набухания свободных и скрепленных образцов эластомеров
Эластомер/жидкость £ Ч £>0х 108, см2/с £1 Я £>! х 108, см2/с
свободные образцы скрепленные образцы
. ПБУ-1/ДБС 0.345 0.67 6.2 0.460 0.50 2.6
ПБУ-1/ДОС 0.461 0.67 1.8 0.541 0.50 1.0
ПБУ-4/ДОС 0.318 0.63 3.4 0.446 0.51 1.7
скрепленных и свободных образцов эластомеров. Б-образная форма кинетических кривых набухания свободных образцов в координатах. (г1/2, ¿»^ свидетельствует об аномальной сорбции. Наиболее существенно она проявляется в сильно набухающей системе ПБУ-2-ДБС. Напротив, кинетические кривые набухания скрепленных образцов на начальной стадии имеют ярко выраженный линейный участок, что является признаком нормальной сорбции.
Количественные характеристики кинетики набухания свободных и скрепленных образцов приведены в табл. 1. Коэффициенты диффузии жидкостей для свободных образцов определяли по уравнению (69), а для скрепленных - стандартным методом [1, 2], основанным на модели нормальной сорбции. Показатель ц оценивали методом наименьших квадратов по начальному участку кинетических кривых с помощью уравнения gl = с^. Величины коэффициента корреляции при этом составили 0.998-0.999.
Из табл. 1 видно, что коэффициенты диффузии жидкостей в свободных и скрепленных образцах имеют разные значения, а показатель для свободных образцов заметно превосходит величину 0.5, которая соответствует нормальной сорбции. Таким образом, экспериментально подтверждается полученный ранее вывод о качественном и количественном различии кинетики набухания скрепленных и свободных образцов.
Типичные кинетические и деформационные кривые свободного набухания эластомеров приведены на рис. 5. Как видно, они имеют 5-образный вид. Количественные оценки (табл. 2) показывают, что процесс свободного набухания исследованных эластомеров соответствует диффузионному режиму степенного типа. При этом концентрационная зависимость коэффициента диффузии жидкостей, оцслшосшя параметром я, является весьма слабой. При определении параметров дигпо кинетическим и деформационным кривым коэффициенты корреляции составили 0.997-0.999 и 0.994-0.998 соответственно. Экспоненциальный режим и режим с обострением при набухании данных эластомеров не обнаружены.
Интересно отметить, что экспериментальные значения параметра с1, характеризующего деформационную зависимость равновесной степени набухания эластомеров при однородном двухосном растяжении, оказались близкими к теоретическому значению (1 = 2/9 = 0.22. Это свидетельствует о правомерности использования теорий Флори для описания процессов набухания эластомеров.
/1/2, мин"2
Рис. 5. Кинетические (а) и деформационные (б) кривые свободного набухания эластомеров в толуоле: 1 - ПБУ-3, 2 - ПБУ-4, 3 - ПБУ-1, 4 -ПБУ-2.
Таблица 2. Характеристики кинетики набухания свободных образцов эластомеров
Эласгомер/жидкосгь i е Я г s d £>0 х 108, см2/с
ПБУ-1/Толуол 0.278 0.68 0.84 0 0.19 140
ПБУ-2/Толуол 0.093 0.83 1.09 0.8 0.24 49
ПБУ-З/Толуол, 0.230 0.78 1.02 0.5 0.24 130
ПБУ-4/Толуол 0.179 0.71 0.97 0 0.27 140
ПБУ-1/ДБС 0.345 0.67 0.92 0 0.27 6.2
ПБУ-2/ДБС 0.128 0.78 1.04 0.5 0.25 2.5
ПБУ-З/ДБС 0.316 0.67 0.83 0 0.19 6.6
ПБУ-4/ДБС 0.267 0.69 0.88 0 0.21 7.4
ПБУ-1/ДОС 0.461 0.67 0.81 0 0.17 1.8
ПБУ-4/ДОС 0.318 0.63 0.81 0 0.22 3.4
ЗАКЛЮЧЕНИЕ
Таким образом, в данной работе предложена геометрически и физически нелинейная теория набухания эластомеров и приведены примеры ее приложения. Следует подчеркнуть, что используемая для описания термодинамических и механических свойств эластомеров модель полимерных сеток Флори не является неотъемлемой частью этой теории. Последняя сформулирована в общем виде и позволяет использовать любые термодинамические модели эластомеров.
Исследованы процессы набухания эластомеров в низкомолекулярных жидкостях, сопровождающиеся большими деформациями полимерной матрицы. Описано напряженно-деформированное состояние материала, порождаемое диффузионным процессом и изучено его влияние на диффузионную кинетику набухания.
Теоретически и экспериментально установлена зависимость диффузионной кинетики набухания эластомеров от условий их внешнего механического нагружения. Показано, что набухание скрепленного слоя эластомера протекает в режиме нормальной сорбции, а кинетика набухания свободного слоя является аномальной.
Изучен механизм аномальной сорбции при набухании эластомеров в низкомолекулярных жидкостях. Установлено, что главная причина аномалий - нестационарный граничный режим. Этот режим порождается деформациями, которые возникают в процессе набухания эластомера и повышают его термодинамическую совместимость с жидкостью. Аномалии кинетики сорбции - ха-
рактерное явление, присущее в той или иной степени все упругим Набухающим материалам.
Из теории следует возможность существования качественно различных диффузионных режимов свободного набухания эластомеров. Конкретный вид процесса определяется комплексом механических, термодинамических и диффузионных свойств материала. Для свободного плоского слоя теория предсказывает возможность существования трех режимов набухания — степенного, экспоненциального и режима с обострением. Экспериментально установлено, что в исследованных эластомерах реализуется степенной режим свободного набухания. Вопрос о существовании двух других режимов набухания в реальных материалах пока остается открытым.
Предложены методы, позволяющие оценить концентрационную зависимость коэффициента диффузии жидкости в эластомере и деформационную зависимость равновесной степени набухания эластомера в условиях однородного двухосного растяжения по кинетическим и деформационным кривым свободного набухания.
СПИСОК ЛИТЕРАТУРЫ .
1. Роджерс К. // Проблемы физики и химии твердого состояния органических соединений. М.: Мир, 1968. С. 229.
2. МалкинАЯ., Чалых А.Е. Диффузия и вязкость полимеров. Методы измерении. М.. Химия, 1979.
3. Petropoulos J.H. // J. Polym. Sei., Polym. Phys. Ed. 1984. V. 22. № 2. Р. 183.
4. Korsmaer R.W., Lustig SR., Peppas NA. // J. Polym. Sei., Polym. Phys. 1986. V. 24. № 2. P. 395.
5. Neogi P., Kim M., Yang Y. // AIChE J. 1986. V. 32. № 7. P. 1146.
6. Sanapoulos M., Roussis P., Petropoulos J.H. // J. Polym. Sei., Polym. Phys. 1995. V. 33. № 7. P. 993.
7. Edwards DA., Cohen D.S. // SIAM J. Appl. Math. 1995. V. 55. № 3. P. 662.
8. Денисюк Е.Я., Тереиштов B.B. II Прикл. механика и техн. физика. 1997. Т. 38. № 6. С. ИЗ.
9. Лурье А.И. Нелинейная теория упругости. М.: Наука, 1980.
10. Flory PJ. Principles of Polymer Chemistry. New York: Cornell Univ. Press, 1953.
11. Васильев В.А., Романовский Ю.М., Яхно В.Г. Автоволновые процессы. М.: Наука, 1987.
12. Самарский А.А., Галактионов В.А., Курдю-мов С.П., Михайлов А.П. Режимы с обострением в задачах для квазилинейных параболических уравнений. М.: Наука, 1987.
13. Тереиштов В.В. // Высокомолек. соед. А. 1995. Т. 37. № 9. С. 1529.
A Nonlinear Theory of the Process of Elastomer Swelling in Low-Molecular-Mass Liquids
E. Ya. Denisyuk and V. V. Tereshatov
Institute of Technical Chemistry, Ural Division, Russian Academy of Sciences, ul. Lenina 13, Perm, 614000 Russia
Abstract—A three-dimensional theory of the elastomer swelling in a low-molecular-mass liquid, constructed as a theory of mass transfer in nonlinear elastic materials, is proposed. The general theory is applied to study the swelling kinetics of a flat elastomer layer at large degrees of swelling, in which case explicit asymptotic expressions describing the initial and final stages of the process are obtained. The possibility of existence of the qualitatively different regimes of swelling is established. Relationships between the kinetics of swelling and the physical properties of the material were studied. The obtained results allow the mechanism of the anomalous swelling kinetics to be explained and new methods to be proposed for the study of diffusion and thermodynamic properties of elastomers.