Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Вязкоупругая релаксация в полимерах

..pdf
Скачиваний:
14
Добавлен:
15.11.2022
Размер:
11.51 Mб
Скачать

мов. Из данных рис. 11 очевидно, что использование смесей может служить эффективным средством создания механически инерт­ ных покрытий, поскольку для некоторых смесей характерно от-

V и с. 9. Вторичные динамические (релаксационные) переходы в сополиме­ рах метилметакрилата (ММА) и этилакрилата. Цифры у кривых указывают мольную долю ММА в сополимере. Сополимеры с высоким содержанием ММА растрескиваются при охлаждении, что не позволяет получить точные данные для стеклообразного состояния.

Температура, К

V и с. 10. Релаксационные максимумы потерь в сополимерах метилметакри­ лата (ММА) и этилакрилата. Величина tg6 = G"/G' проходит через максимум приблизительно при температуре стеклования. Цифры у кривых указывают мольную долю ММА в сополимере.

сутствие существенного влияния нагревания или охлаждения на их механические свойства. В едином материале можно соче­ тать домены, находящиеся в стеклообразном и высокоэластиче­ ском состоянии, так что в зависимости от соотношения ингре­

диентов можно создать наполненный каучук или пластифициро­

ванное

стекло.

Для

установления зависимости между G' и со0 следует

предположить, что нить в приборе имеет цилиндрическую форму. Тогда]

о-

(15)

Длина нити I и момент инерции I известны, радиус г нити с по­ крытием можно вычислить из веса единицы ее длины, но с не

Температура, К

Р и с. 11. Релаксационные максимумы потерь в смесях полиметилметакри­ лата (ПММА) и полиэтилакрилата. Положение температур максимумов по­ терь пс зависит от состава. Цифры у кривых указывают мольную долю ПММА в смеси.

очень высокой точностью, поскольку форма нити отклоняется

от строго цилиндрической и ее покрытие может

быть не вполне

равномерным. Выражение для модуля потерь имеет вид

Г>п 2// 2

( 16)

 

откуда

 

tg e = 4 - .

(|7>

где Д — логарифмический декремент затухания свободных коле­ баний. На значения tg6 из указанных выше неточностей влияет лишь равномерность покрытия на нити.

Соотношение между компонентами динамического модуля и статистическим модулем* G имеет вид [201

(1)2Т2

G'=G 1 02 2 ДО, (18)

+ С т

G

сот

(19)

со-т

 

 

где ДG — разность между статическим модулем G и асимптоти­ ческим значением G' при высоких частотах. G' возрастает с по-

Р м с. 12. Оценка

времени релаксации по точке переп

Га зависимости

G' от

частоты (о и по максимуму зависимости G"

от со.

вышением со от величины G (которая может быть равна пулю) до некоторого предельного значения. Предполагая наличие лишь одного релаксационного механизма (одного значения т), получаем на зависимости G' от w точку перегиба при со —- 1/т, а на за висимости О" от о — максимум (рис. 12).

Метод измерения отраженных импульсов

Когда сдвиговая волна создается

в плоской подложке, величи­

на деформации в граничном

слое

пропорциональна тому, на­

сколько согласован сдвиговый

механический импеданс Z с им­

педансом подложки. Если чувствительный датчик установлен так, чтобы измерить энергию отраженного импульса, то по отно­ сительному ослаблению интенсивности колебаний можно судить о свойствах исследуемого материала, поскольку ослабление

* Здесь под статическим модулем имеется в виду равновесный модуль, характеризующий величину напряжений, сохраняющихся в материале после завершения релаксационных процессов.— Прим. ред.

пропорционально Z, а твердые подложки характеризуются высо­

кими значениями

импеданса.

 

 

Было показано

[11], что оптимальное значение угла падения

составляет

11°

В этом случае использование

подложки трапе­

цеидальной

формы, изготовленной из

материала с известным

импедансом,

обеспечивает выгодные

условия

распространения

сдвиговых волн в приклеенной пленке (рис. 13).

Используемую подложку изготовляли из плавленого кварца

(Zq = 8,29* 10б

г/с-см2). Частоту измерений

варьировали от 2

до 100 МГц. Импульсы длительностью 4 мкс, преобразованные и прошедшие через подложку к плоскости покрытия, собирали и

Р и с. 13. Схема, иллюстрирующая

метод измерения

характеристик пленки

по отражению сдвиговых волн на

границе раздела

покрытие

подложка.

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

Величина А связана с коэффициентом отражения г формулой

lg r= —А/40.

(20)

Изменение значения Д как функции продолжительности от­ верждения давало характерные S-образные зависимости [21], которые при желании можно линеаризовать [22] и использовать для исследования кинетики процесса удаления летучих ве­ ществ из пленки. Абсолютные значения импеданса рассчитывали по упрощенной формуле, следующей из более общего выраже­ ния [23], в предположении неизменности фазы импульса на границе раздела

Z ^ l,58-105[j ^ f ,

(21)

где численный коэффициент находили как Z?cos0 (0 — угол па­ дения). Величину G' вычисляли по формуле

G' ^ pZ2 « 2,5 • 1010

)2

(22)

Здесь принято, что р ^ 1 . Следует помнить, что точность исполь­ зуемой формулы резко снижается с уменьшением г.

Без независимых

измерений сдвига фазы на

поверхнос­

ти раздела невозможно определить G" покрытия, так как оно

связано с величиной

coscp, где <р — фазовый угол

отраженной

волны.

 

 

Применение метода измерений А для контроля за переходом

из жидкого состояния в твердое вполне приемлемо,

если не за-

Поглощение А, дБ

Р и с. 14. Зависимость модуля упругости от поглощения сдвиговых импуль­ сов. Цифры у кривых указывают величину отношения фазового угла (в гра­

дусах) к поглощению (в децибелах).

бывать об ограничениях, возникающих при использовании еди­ ничного измерения для характеристики свойства, зависящего от двух материальных параметров. До тех пор пока значение <р невелико, расчетное значение G' относительно нечувствительно ни к чему, кроме интенсивности поглощения [24], как это еле-

дует из рис. 14. Когда А = 23 дБ, модуль G' становится пол­ ностью независимым от <р. Непрерывный контроль за структур­ ными превращениями можно осуществлять, наблюдая за изме­ нением А во времени [25].

Иногда S-образная форма зависимости А от времени при вы­ сушивании пленки искажается из-за накоплений в пленке на­ пряжений, которые могут исчезнуть только лишь при образова­ нии разрывов, происходящих вследствие нарушения адгезионных

икогезионных связей [26]. Наблюдаемая картина явлений в дей­ ствительности связана с неполной релаксацией напряжений в пленке при приближении к отвержденному состоянию [27], поскольку пленка не может одновременно оставаться и сплошной,

иприлипшей к подложке.

ЗАКЛЮЧЕНИЕ

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

Правильное сочетание различных механических свойств яв­ ляется необходимым условием образования пленок. Величиной, определяющей напряжения в материале, является модуль упру­ гости.

Регулирование механических характеристик пленки в про­

цессе ее нанесения и формования

обеспечивается

изменением

геометрической формы изделия, скорости нанесения

покрытия,

а также варьированием его состава.

Цель такого регулирования

состоит в создании покрытия с требуемым соотношением демпфи­ рующих и упругих характеристик.

Модуль потерь как жидкой, так и отвержденной пленок мень­ ше, чем их модуль упругости. В процессе формования пленки модуль упругости возрастает. Максимум потерь соответствует времени наблюдения по порядку величины, близкому ко време­ ни релаксации покрытия. Продолжительность наблюдения варьи­ руется в пределах от 10“8 до 10® с, причем покрытие может прояв­ лять вязкоупругие свойства, если за время экспериментов успе­ вает произойти изменение конформации цепи. В противном случае может происходить либо разрушение материала, либо осуществ­ ляется процесс течения.

Оптимальное отношение между временем наблюдения и вре­ менем релаксации близко к единице. Это условие выполняется в момент, когда пленка при отверждении проходит через кау­ чукоподобное (или липкое) состояние.

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

СПИСОК ЛИТЕРАТУРЫ

 

 

 

 

 

 

1.

Reiner М ., Phys. Today, 17, 62 (1964).

Math-Naturw.

K1.,

Abt. II, 69,

2.

Stefan М. У., Sitzber. Akad. Wiss. Wien,

3.

713 (1874).

 

 

 

 

 

 

Griffith A. A., Phil. Trans. Roy. Soc. London, Ser. A, 221, 163 (1921).

4.

Banks U". Я ., Mill С. C., J. Colloid Sci., 8, 137 (1953).

 

 

5.

Gaskell R. £ ., Trans. ASME, 72, 334 (1950).

 

 

 

 

6. Miller J. C. , Myers R. R., Trans. Soc. Rheology, 3, 77 (1959).

(1961).

7.

Myers R. R., Hoffman R. D., Trans. Soc. Rheology, 5, 317

8. Myers R. R., Tappi, 46, 745 (1963).

 

18,

101 (1951).

 

 

9.

Bergen J. T ., Scott G. Yf., J. Appl. Mech.,

1963.

10.

Ферри Дж.y Вязкоупругие свойства полимеров

ИЛ, М.,

11.

Mason tt;'. P .t Backer U;? О., McSkimin Н. У., Heiss У. Я ., Phys. Rev., 75,

 

936 (1949).

 

 

 

 

 

 

12. Kutz D. Я ., Senior Thesis, Lehigh University, Bethlehem, Pa. (1959).

13.

Nielson L. E. , Rev. Sci. Instr., 22, 690 (1951).

 

 

 

 

14.

Lewis A. F., Gillham J. K. , J. Appl. Polymer Sci., 6, 422 (1962).

15.

Gillham У. K ., Lewis A. F., Nature, 195, 1199 (1962).

 

 

16.

Lewis A. F., Tobin M. C., Trans. Soc. Rheology,

6, 27 (1962).

A, 2, 5017

17.

Stratta J. J.y Reding F. P.y Faucher

J. A.,

J.

Polymer Sci.,

18.

(1964),

 

 

 

 

 

 

Gillham J. /(., Science, 139, 494 (1963).

 

 

 

 

 

19.

Myers R. R., Knauss С. У., Shroff R. N .y Proc. 5th Intern. Congr. Rheology,

20.

1, 473 (1969).

 

 

 

 

 

 

Passaglia E., SPE Trans., 4, 1 (1964).

 

Technol.,

33,

940 (1961).

21.

Myers R. R .y Offic. Dig. Federation

Soc. Paint

22.Myers R. R., Schultz R. Д\, Offic. Dig. Federation Soc. Paint Technol., 34, 801 (1962).

23.O'Neill H. T .t Phys. Rev., 75, 928 (1949).

24.

Myers R. R.,

Knauss

С. У., Polymer Colloids, Pienum,

New York

25.

(1971).

KUmek J.,

Knauss C. J . t J. Paint Technol., 38,

479 (1966).

Myers R. R.,

26.Myers R. R., Schultz R. /(., J. Appl. Polymer Sci., 8, 755 (1964).

27.Myers R. R., Knauss С. У., J. Paint Technol., 40, 315 (1968).

ЧИСЛЕННОЕ ПРЕОБРАЗОВАНИЕ ФУНКЦИОНАЛЬНЫХ ЗАВИСИМОСТЕЙ ДЛЯ ВЯЗКОУПРУГОГО ТЕЛА МЕТОДОМ ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ

Шен, Джемисон, Дрэнделл

S h е п

М., J a m i e s o n R. Т., Department

of Chemical

Engineering, Univer­

sity of

California,

Berkeley, California 94720;

D г a n d e 11

M., University Sys­

tems, Los Angeles,

California 90025

 

 

ВВЕДЕНИЕ

Математическая теория линейной вязкоупругости развита весьма детально, в частности получены точные уравнения, свя­ зывающие друг с другом различные функции, которые характе­ ризуют вязкоупругое поведение материала [1—3]. Известны также два приближенных метода расчета этих функций: прямой пересчет известной функции в другую, представляющую инте­ рес, и пересчет с использованием на промежуточной стадии вы­ числений спектров времен релаксации и запаздывания.

Результаты экспериментов по релаксации напряжения обыч­ но описываются аналитически с помощью обобщенной модели Максвелла, которая приводит к следующим формулам [4]:

 

( 1)

 

( 2 )

—оо

 

Е" (ш) = f [Я(х)о>т/(1 + co2x2)]d In т.

(3)

Здесь £(/), £'(а>) и Е''(ы) — зависящий от времени релаксацион­ ный модуль и зависящие от частоты динамические модуль упру­ гости и модуль потерь; Ее — равновесный модуль, равный для несшитых полимеров нулю; Н(т) — функция распределения вре­

мен

релаксации

(т). Частота и

время связаны соотношением

со =

1/т.

 

 

При описании экспериментов по ползучести обобщения модель

Фойхта приводит

к следующим

соотношениям [4]:

СО

D (i) = Dg + f L (т)[ 1- pxp (— //x)]d In x + t/тц,

(4)

D '(a)= D g+

j [L (т)/( 1 + co2T2)]d In т,

(5)

оо

— -О

 

 

 

D"(o>)= j>[L(T)coT/(l + a> V )]dlnT + -^-.

(6)

— СО

 

 

Здесь D(t), D'(a>), Z)"(со) и

соответственно податливость

при

ползучести, действительная и мнимая компоненты динамической податливости и податливость, отвечающая стеклообразному со­ стоянию; L(T) — функция распределения времен запаздывания; т], — вязкость при установившемся течении (для сшитых поли­ меров т^-^оо).

Из соотношений (1)—(6) ясно видно, что по известным функ­ циям распределения времен релаксации или запаздывания тео­ ретически можно рассчитать податливость или модули. Однако L(т) и Н(т) не могут быть измерены непосредственно, а рассчи­ тываются из экспериментальных наблюдений таких зависимостей, как E(t), D(t), £'(со) и т. д. Более того, не существует удобных и точных соотношений, по которым из этих измеренных зависи­ мостей можно было бы рассчитать функции #(т) и L(T). Было предложено огромное число приближенных соотношений, позво­ ляющих определять Н(т) по любому модулю и Ь(т) по любой подат­ ливости. Обобщение этого материала с указанием погрешностей каждого метода было выполнено Ферри [4].

Более надежны соотношения, позволяющие пересчитывать одни зависимости модулей или податливостей в другие, минуя использование спектров распределения времен релаксации или запаздывания. Применимость и надежность этих соотношений была недавно рассмотрена* в работе Шварцля и сотр. [5]. Одна­ ко возможность использования прямых пересчетных формул в значительной мере зависит от интервала времен, для которых были получены экспериментальные данные, и вида измеренных функций. Кроме того, исходные экспериментальные данные должны быть представлены в определенной форме. Во многих случаях точность пересчета остается неизвестной и требуется проведение дополнительных перекрестных расчетов. Поэтому необходимо сформулировать общую для всех расчетов такого рода методику, позволяющую использовать технику программи­ рования для счетнорешающей машины, не требующую пубо'

* См. также статьи, посвященные анализу численных методов взаимны^ преобразований переходных и динамических функций, характеризующих дк

ханические свойства вязкоупругих материалов: Schwarzl F.

8

, Rheol. Acta

8, 6 (1969) и Struik L. C. £ ., Schucirzl F. R., Rheol. Acta,

, 134 (1969).—

Прим. ped.

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

МЕТОД ЛИНЕЙНОГО ПРОГРАММИРОВАНИЯ

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

min

 

 

 

 

 

или

z= c 1x1 + c2x2-{------- Vcnxn

(7)

max

 

 

 

 

 

при соблюдении

связующих

ограничений

 

 

an xi + а12хо + • • •

+ а1пхп ^

Ьг

 

 

 

 

 

( 8 )

 

^mlXl ~Т“ 6lm2X2 " f “ * *

’ “ Ь

&тггХп ^

Ьт

 

X j^ 0,

/= 1 ,

 

/7.

 

Предполагается, что все ср аи и

bt — величины известные.

Все переменные ограничены областью положительных значе­ ний от нуля включительно. Это условие известно под названием «ограничение неотрицательности». Согласно терминологии ли­

нейного программирования,

z — искомая функция, Cj — со­

ставляющие результирующего

вектора с и bt — составляющие

правостороннего вектора Ь. Решение задачи ЛП основывается на стандартной программе «ALPHAC», разработанной в вычисли­ тельном центре Калифорнийского университета для машины СDC 6400.

Исходным моментом для решения была запись соотношений

(1)—(6) в виде дискретных сумм, что непосредственно вытекает из обобщенной модели Максвелла и Фойхта [4].

Исходным методом для решения была запись соотношений

Е (0 = Ее+ 2 Ej exp (— t/Tj),

/-i.