CC BY
109
Применение общей термодинамической теории к решению проблем механики
УДК 539.18:541.64
Е.Я. Денисюк Институт механики сплошных сред УрО РАН
МЕХАНИКА И ТЕРМОДИНАМИКА ДЕФОРМИРОВАНИЯ ВЫСОКОЭЛАСТИЧНЫХ ПОЛИМЕРНЫХ СЕТЧАТЫХ МАТЕРИАЛОВ В ФИЗИЧЕСКИ АГРЕССИВНЫХ ЖИДКОСТЯХ 1
Abstract
The study is developed to mechanics and thermodynamics of deforming hyper elastic bounded swelling polymer materials containing dissolvent or interacting with dissolvent medium. The problem of description of stress and strain state in incompressible material being in thermodynamic equilibrium with dissolvent and subjected to an arbitrary static load is considered. To describe such case the general system of equations, constitutive relations and boundary conditions is stated. To model elastic and thermodynamic properties of polymer nettings the specific constitutive relations based on the Florey - Rener - Huggins theory are suggested. The theory is illustrated with an example of equilibrium swelling and elastic response of polymer netting which is subjected to a uniaxial tension in a solvent media.
Химически сшитые эластомеры и полимерные гели - это высокоэластичные полимерные материалы, которые могут испытывать гигантские упругие деформации, исчисляющиеся в сотнях и тысячах процентов. Другое их уникальное свойство состоит в том, что они способны поглощать низкомолекулярные растворители, многократно (в десятки и сотни раз) увеличиваясь в объеме. При этом они сохраняют свою форму и способность к обратимым упругим деформациям. Это явление называется ограниченным набуханием и объясняется молекулярной структурой данных материалов - они представляют собой пространственную полимерную сетку, состоящую из гибких макромолекулярных цепей, соединенных химическими связями [1].
Способность полимерных гелей поглощать, а затем высвобождать жидкость широко используется во многих современных технологиях, например, в биотехнологии (сепарация протеинов), в медицине и фармакологии (лекарственные гели), в сельском хозяйстве и ландшафтном дизайне (увлажнители почвы), в биохимии (гелевые мембраны, служащие для разделения и анализа биорастворов) и т.д. [2-5]. Что касается эластомеров, то эти материалы часто эксплуатируются в физически агрессивных средах
- органических растворителях и их парах. Поглощение растворителя вызывает набухание материала, приводит к изменению его физико-механических свойств, порождает в нем внутренние напряжения и в конечном итоге может вызвать разрушение изделия [6]. К этому следует добавить, что в физико-химии полимеров процессы набухания - это один из важнейших «инструментов» исследования физикохимических свойств и структуры полимерных сеток [1].
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Администрации Пермской обл. (проект РФФИ -Урал № 04-01-96039).
Применение общей термодинамической теории к решению проблем механики
Важной особенностью поведения эластомеров и полимерных гелей при их взаимодействии с физически агрессивными жидкостями является то, что деформации материала, вызванные внешней механической нагрузкой, могут существенно влиять на термодинамическое равновесие между полимерной сеткой и растворителем. Например, известно, что эластомер или гель в растянутом состоянии поглощает большее количество растворителя, чем в ненагруженном. Более того, неоднородное поле напряжений вызывает перераспределение растворителя, а неоднородное распределение растворителя порождает сложно-напряженное состояние материала. Все эти эффекты проявляются в полной мере, когда материал подвергается одновременному воздействию физически агрессивной жидкости и внешнему механическому нагружению.
В физико-химии полимеров термодинамика набухания и механическое поведение высокоэластичных сетчатых материалов в среде растворителя в основном изучены только при простых видах механического нагружения (одноосное растяжение или сжатие). Что касается нелинейной теории упругости, в которой изучению механических свойств эластомеров традиционно уделяется значительное внимание, то в этой дисциплине явление набухания материалов и их взаимодействие с растворителями практически не рассматриваются. Основная причина такого положения дел заключается в отсутствии общего теоретического формализма, позволяющего формулировать конкретные проблемы, касающиеся описания механического и термодинамического поведения высокоэластичного материала, содержащего растворитель или взаимодействующего со средой растворителя.
Общая теория, описывающая неравновесные процессы набухания и диффузии жидкостей в упругих материалах, сопровождающихся конечными деформациями упругой матрицы (теория механодифузии), развита в работах [7,8]. В качестве частного случая она содержит в себе механику и термодинамику деформирования упругих материалов, содержащих жидкость. Цели данной работы заключаются в следующем: 1) используя указанное обстоятельство, построить общую систему уравнений и определяющих соотношений, описывающих термодинамически равновесное напряженно-деформированное состояние упругого, набухающего материала, находящегося в контакте с жидкостью и подвергнутого произвольному статическому нагружению; 2) на достаточно простых моделях изучить основные закономерности механического и термодинамического поведения эластомеров и полимерных гелей при их деформировании в среде растворителя.
При описании физико-механических свойств эластомеров и полимерных гелей будем полагать, что полимерная матрица является однородным, изотропным, упругим несжимаемым материалом. Жидкость будем также рассматривать как несжимаемую среду. Приближение несжимаемости в данном случае означает, что единственной причиной объемных деформаций материала является поглощение жидкости. Оно вполне оправдано и объясняется тем, что под действием сравнительно небольших механических нагрузок эластомеры и полимерные гели могут испытывать значительные деформации сдвига, которые на много порядков превосходят объемные деформации, связанные с собственной сжимаемостью полимера и растворителя.
Рассмотрим упругий материал, который способен поглощать низкомолекулярную жидкость, называемую в дальнейшем растворителем. Если образец такого материала поместить в среду растворителя, то в результате поглощения растворителя, он перейдет в термодинамически равновесное набухшее состояние, в котором поглощенный растворитель равномерно распределен в объеме образца, а внутренние механические напряжения отсутствуют. Если на образец не действуют внешние механические силы, а параметры среды (температура, давление и т.д.) остаются неизменными, то в этом
Применение общей термодинамической теории к решению проблем механики
состоянии он может находиться сколь угодно долго. Выберем это состояние в качестве отсчетной конфигурации, которую зададим уравнением
Г = г(4, ц2,ц3), (1)
где (ц1, ц2, ц3) - система материальных координат, связанная с упругой матрицей.
Если образец подвергнуть статическому нагружению, то возникшее в первый момент напряженно-деформированное состояние материала в общем случае уже не будет термодинамически равновесным. Однако в результате дополнительного поглощения и перераспределения растворителя с течением времени образец перейдет в новое термодинамически равновесное напряженно-деформированное состояние, конфигурация упругой матрицы в котором описывается соотношением
Я = Я (ц1, ц2, ц3). (2)
Система уравнений, описывающая это конечное состояние, имеет следующий вид [7]:
V-Т = 0, |Д = 0, (3)
где Т - тензор напряжений Коши; д - химический потенциал растворенной в материале жидкости. Здесь и далее мы полагаем, что вся система находится в изотермических условиях. Первое уравнение выражает условие механического равновесия образца. Второе уравнение является условием термодинамического равновесия системы "упругий материал - растворитель". Оно означает, что в каждой точке образца химический потенциал растворителя одинаков и равен химическому потенциалу жидкой фазы растворителя, который, без ограничения общности, можно положить равным нулю.
Условия нагружения материала задаются граничными условиями: на части границы тела могут быть заданы перемещения w:
Я (ч\ ц2, ц3) = г (ц1, ц2, ц3) + w( ц1, ц2, ц3), (4)
а на другой части границы - поверхностное распределение сил Г:
N - Т = Г , (5)
где N - вектор внешней нормали.
Уравнения равновесия (3) справедливы как для несжимаемых материалов, так и для сжимаемых. В случае несжимаемого материала и растворителя к системе уравнений (3) и граничных условий (4), (5) необходимо добавить условие несжимаемости, которое с учетом выбора отсчетной конфигурации имеет следующий вид [7,8]:
з = Фе/Ф , (6)
где ф - объемная доля полимерной матрицы в актуальной (текущей) конфигурации; ф Е
- ее значение в отсчетной конфигурации, в качестве которой здесь принято равновесно набухшее состояние ненагруженного материала; 3 = ёй Г - третий инвариант
градиента деформаций, который характеризует объемные деформации материала,
вызванные поглощением или перераспределением растворителя. Градиент деформаций
Г определяется согласно [9].
Химический потенциал растворителя и тензор напряжений Коши для системы " несжимаемый материал - несжимаемая жидкость" имеют следующий вид [7]:
й = (д//дп)8 + V, р (7)
Т = 5,(-1/,Е)2(В2 -±/,(В2)Е)-рЕ (8)
5, = 2((, + /, д//дІ2) 52 = 2#/#2 (9)
/ = /(/,, Іг, и,, и,), 00)
Здесь / - плотность свободной энергии материала и растворенной в нем жидкости; р -гидростатическое давление, которое в данном случае является лагранжевым
Механика и термодинамика деформирования высокоэластичных полимерных сетчатых материалов
Применение общей термодинамической теории к решению проблем механики
параметром связи; Е - единичный тензор; Ух, У2 - мольный объем растворителя и вещества упругой матрицы соответственно. Мольная концентрация растворителя п1 и вещества упругой матрицы п2 с помощью условия несжимаемости (6) можно выразить через объемную долю упругой матрицы в отсчетной конфигурации и инвариант з :
П1 = КД1 -Фе/З), П2 = -Фе/З). (11)
Зависимость свободной энергии, химического потенциала и тензора напряжений от деформаций материала в соотношениях (7)-(10) выражена в терминах меры деформаций искажения формы
Б = Б/32/3 (12)
и ее инвариантов
/ = /,(В) = 1,/32'\ ¡2 = /2(В) = 1,13“\ 03)
где Б = ГТ - Г - мера деформации Фингера [9]; /, /2- первый и второй инварианты
меры Фингера соответственно.
Таким образом, определение термодинамически равновесного напряженно-деформированного состояния материала, находящегося в контакте с жидкой фазой растворителя, сводится к решению системы уравнений (3) с граничными условиями вида (4), (5). В случае несжимаемого материала и несжимаемой жидкости эту проблему можно представить в форме статической задачи нелинейной теории упругости для сжимаемого материала. Для этого с помощью второго уравнения (3) и соотношения (7) выразим гидростатическое давление и запиш)ем его в виде
р = п(ф (14)
где
п = -К1-1(д//д п1) (15)
- термодинамическая функция состояния системы. Подставив (14) в выражение для
тензора напряжений (8), получаем
Т = П, (16)
где тензор второго ранга
П = ^(Б - ^ ЛЕ) - 5 2 (В2 -1 /1(В2 )Е) - пЕ. 07)
С учетом (16) уравнение механического равновесия (3) и граничное условие (5) примут следующий вид:
V•П = 0, (18)
N-П = Г . (19)
Выясним физический смысл величины п. Для этого рассмотрим набухший образец материала, концентрация растворителя в котором меньше равновесной. Если образец заключить в жесткую оболочку, проницаемую для растворителя, и поместить его в растворитель, то образец будет стремиться поглотить некоторое дополнительное количество жидкости, недостающее для достижения термодинамического равновесия. Однако этому будет препятствовать жесткая оболочка, и образец окажется в условиях всестороннего сжатия. В такой системе возникнет давление, которое определяется соотношениями (14)—(15). В термодинамике растворов это явление называется осмосом, а давление п - осмотическим давлением [10]. Если несжимаемый материал подвергнут произвольному статическому нагружению и находится в
термодинамическом равновесии с жидкой фазой растворителя, то согласно (16) его
тензор напряжений совпадает с тензором (17), который естественно называть
осмотическим тензором напряжений.
Применение общей термодинамической теории к решению проблем механики
С помощью условия несжимаемости (6) зависимость осмотического тензора напряжений от концентрации растворителя можно выразить в терминах инварианта 3. В результате этого структура осмотического тензора напряжений полностью совпадет со структурой тензора напряжений упругого сжимаемого материала. Таким образом, проблема описания напряженно-деформированного состояния несжимаемого материала, находящегося в термодинамическом равновесии с растворителем, сведена к краевой задаче для уравнения (18) с граничными условиями (4) или (19), которая по своей форме полностью аналогична статической задаче нелинейной теории упругости для сжимаемого материала. При этом сжимаемость материала связана с поглощением или перераспределением растворителя, а роль тензора упругих напряжений в ней играет осмотический тензор напряжений, который содержит всю информацию о термодинамических и упругих свойствах системы "упругий материал - растворитель".
Для применения общей теории к описанию равновесных состояний и термодинамически равновесных режимов деформирования полимерных сеток в среде растворителя необходимо иметь конкретное выражение для осмотического тензора напряжений. Если свободная энергия системы "полимер - растовритель" известна и представлена в виде (10), то осмотический тензор напряжений находится по формулам (9), (15), (17).
Вывод выражений свободной энергии - это предмет статистической физики полимеров и полимерных растворов, в которой к настоящему времени разработано большое количество теоретических и эмпирических моделей высокоэластичности, призванных описывать различные аспекты механического и термодинамического поведения полимерных сеток. Для моделирования упругих и термодинамических свойств полимерных сеток можно обратиться к этим моделям.
Одной из первых и до сих пор наиболее «популярных» моделей полимерных сеток является теория Флори - Ренера - Хаггинса [11]. Она основана на классической теории высокоэластичности и теории полимерных растворов. Соответствующее выражение свободной энергии, представленное в терминах инвариантов меры деформации искажения формы, имеет следующий вид:
/ = КТ (п11п(1 - ф) + х(1 - фп + * п ф-2'3 /.), (20)
где Т - абсолютная температура; К - универсальная газовая постоянная; п1, п2 -
мольные концентрации растворителя и цепей полимерной сетки соответственно; х -параметр Флори - Хаггинса, характеризующий взаимодействие между полимерными цепями и растворителем.
Применив к выражению (20) формулы (9), (15), (17), получаем выражения для тензора напряжений Коши
Т = КТУ2"1ф1'3 (Б - 3 /1е)- рЕ, (21)
осмотического давления
п = -КТУ1"1(1п(1 -ф) + ф + хф 2 + £ 2 _1ф1'3 /1), (22)
и осмотического тензора напряжений
П = КТУ^2 Б - цЕ, (23)
ц = -кту1-1(1п(1 -ф) + ф + хФ 2). (24)
Величина ц имеет смысл осмотического давления раствора линейных (не сшитых) полимерных цепей; безразмерный параметр 2 = У2/Ух характеризует относительный объем цепи полимерной сетки по отношению к объему молекулы растворителя. Для высокоэластичных полимерных сеток и низкомолекулярных растворителей 1.
Применение общей термодинамической теории к решению проблем механики
Рассмотрим асимптотическое поведение осмотического тензора напряжений при больших степенях набухания полимерной сетки, то есть при ф~1. В силу малости ф можно воспользоваться разложением
1п(1 -ф) = -ф- 1 ф2 - 3 ф3 + 0(ф4). (25)
Подставив его в (24) и ограничившись квадратичными по ф членами, получим
ц = КТГ1~1Ьфа . (26)
Отсюда находим осмотический тензор напряжений
П = КТУ-1^-1ф1'3Б - ЬфаЕ), (27)
где а = 2 и Ь = 1'2 - х .
Теория Флори - Ренера - Хаггинса, а следовательно и асимптотические выражения (26) - (27) основаны на теории среднего поля и не учитывают наличие корреляций между полимерными цепями, которые имеют место в полуразбавленных растворах [12]. Более точное выражение для осмотического давления растворов полимеров предложено Де Клуазо в рамках концепции скейлинга [12]. Согласно его теории ц~фа , где скейлинговый показатель а = 9'4.
Соотношения (26) - (27) справедливы при Ь > 0, поскольку только в этом случае они описывают ограниченное набухание полимерных сеток. В рамках теории Флори -Ренера - Хаггинса константа Ь положительна в случае растворителей, для которых
параметр Флори - Хаггинса х < 1'2. Такие растворители называются хорошими [12].
Если х = 1'2, что имеет место для так называемых 9-растворителей [12], константа Ь = 0 . В этом случае в разложении (25) необходимо учитывать кубический член. Это приводит к тем же асимптотическим соотношениям (26)-(27), в которых Ь = 1/3 и а = 1'3. Дальнейшее снижение термодинамического качества растворителя, то есть рост параметра Флори - Хаггинса и превышение им порогового значения 1'2, приводит к резкому уменьшению степени набухания материала и стремлению ф ^ 1. Растворители с х > 1'2 являются плохими для данного полимера, и для описания систем «полимерная сетка - плохой растворитель» в рамках теории Флори - Ренера - Хаггинса необходимо обратиться к выражениям (22)-(24).
Продемонстрируем работу развитого выше формализма на простом примере одноосного растяжения образца эластомера или полимерного геля в среде
растворителя. Упругие и термодинамические свойства материала будем описывать в рамках модели Флори - Ренера - Хаггинса. Пусть в исходном состоянии образец не нагружен и находится в термодинамическом равновесии с растворителем. Тогда
Б = Е 11 = 12 = 3 ф = ф Е
и осмотический тензор напряжений (17) П = -пЕ . Условие термодинамического
равновесия сводится к уравнению п(фЕ ,3, 3) = 0. Подставив сюда выражение (22), получаем известное в физической химии полимеров уравнение Флори - Ренера [1,11]:
1п(1 -ф Е ) + ф Е +хф Е + 2 1фЕ/3 = ^ (28)
которое позволяет рассчитать равновесную объемную долю полимера фЕ в условиях изотропного набухания.
Не вынимая образец из растворителя, подвергнем его одноосному растяжению вдоль оси г. Будем полагать, что его боковая поверхность нагружению не
подвергается. При достаточно быстром растяжении количество содержащегося в нем растворителя практически не изменится. При этом его напряженно-деформированное состояние не будет термодинамически равновесным. В такой ситуации образец ведет себя как несжимаемый материал. Из соотношения (21) легко рассчитать его упругую
Применение общей термодинамической теории к решению проблем механики
реакцию. В результате приходим к хорошо известной формуле классической теории высокоэластичности
а3 = G(^2 -Г1), (29)
где X - степень растяжения образца; о3 - продольная компонента тензора напряжений Коши; G = ^ТУ2-1ф1/3 - модуль сдвига набухшего до равновесия материала.
После дополнительного поглощения растворителя и завершения диффузионных процессов образец перейдет в новое термодинамически равновесное состояние, в котором упругая реакция образца равна продольной компоненте осмотического тензора напряжений, при этом П1 = П1 = 0. Отсюда из соотношений (26), (27) получаем следующие уравнения:
а3 = П3 = G(k2 - J/А,), (30)
П1 =П1 = RTV1~1Z“Ч^Г1 - ^(фE / J) = 0. (31)
Последнее уравнение служит для определения равновесной объемной степени дополнительного набухания J растянутого образца.
В случае хорошего или 9-растворителя, когдаф = фE/J ~1, можно обратиться к соотношениям (26)-(27). Тогда уравнение (31) можно разрешить относительно J и найти асимптотическую зависимость равновесной степени набухания образца от степени его растяжения,
J=^1/а. (32)
Подставив ее в соотношение (30), находим равновесную упругую реакцию образца
а3 = П3 = G(X2 -Г1+1/а). (33)
Отсюда следует, что при деформировании полимерных сеток в хороших и 9-растворителях они могут испытывать весьма большие объемные деформации, вызванные дополнительным поглощением жидкости. В результате этого упругая реакция образца, деформируемого в растворителе, оказывается меньше упругой реакции образца, деформируемого в атмосфере. Это различие тем больше, чем лучше растворитель. Растяжение полимерной сетки в плохом растворителе не приводит к существенному поглощению растворителя, и ее равновесная упругая реакция практически такая же, как и при растяжении в атмосфере.
Таким образом, проблема описания напряженно-деформированного состояния высокоэластичного полимерного сетчатого материала, находящегося в термодинамическом равновесии с растворителем, может быть сформулирована как статическая задача нелинейной теории упругости для сжимаемого материала. Рассмотренные примеры показывают, что с точки зрения механики систему «полимерная сетка - хороший растворитель» можно рассматривать как физическую модель высокоэластичного материала, обладающего гигантской сжимаемостью. Более того, сжимаемость такого "материала" можно легко регулировать путем выбора растворителя.
Библиографический список
1. Тагер А.А. Физико-химия полимеров. - М.: Химия, 1978.- 484 с.
2. Полимерные системы для контролируемого выделения биологически активных соединений / Л.И. Валуев, Т.А. Валуева, И.Л. Валуев, Н.А. Платэ // Успехи биологической химии.- 2003.- Т. 43.- С. 307-328.
3. Kim J.J., Park K. Smart hydrogels for bioseparation // Bioseparation.- 1999.- Vol.7. -P.77-184.
Применение общей термодинамической теории к решению проблем механики
4. Nishino M., Gong J., Osada Y. Polymer gels as a chemical valve // Bioseparation.-1999.- Vol.7.- P. 269-280.
5. Osada Y., Gong J.P. Stimuli-responsive polymer gels and their application to chemomechanical systems // Progr. Polym. Sci.- 1993.- Vol. 18.- P. 187-226.
6. Моисеев Ю.В., Заиков Г.Е. Химическая стойкость полимеров в агрессивных средах.- М.: Химия, 1979.- 288 с.
7. Денисюк Е.Я., Терешатов В.В. // Прикл. механика и технич. физика.- 1997.- Т. 38, № 6.- С. 113-129.
8. Денисюк Е.Я., Терешатов В.В. // Высокомолек. соед. А.- 2000.- Т.42, № 1.-С.71-83.
9. Лурье А.И. Нелинейная теория упругости.- М.: Наука, 1980.- 512 с.
10. Ландау Л.Д., Лифшиц Е.М. Статистическая физика. Т.5.- М.: Наука, 1976.-584 с.
11. Flory P.J. Principles of polymer chemistry.- New York: Cornell univ. press., 1953.594 p.
12. Де Женн П. Идеи скейлинга в физике полимеров.- М.: Мир, 1982.- 368 с.
Получено 5.07.2005.
