Вязкоупругая релаксация в полимерах
..pdfмов. Из данных рис. 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.