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

Механика композитных материалов 5 1980

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

УДК 678.02:539.21

И. П. Бриедис, А. М. Савельев, В. С. Василенко

ВЛИЯНИЕ ИЗМЕНЧИВОСТИ МОЛЕКУЛЯРНОГО СТРОЕНИЯ СИНТЕЗИРУЕМОГО ПОЛИЭТИЛЕНА НА ЕГО ФИЗИЧЕСКИЕ СВОЙСТВА*

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

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

На рис. 1 приводятся значения показателя текучести расплава (ПТР) проб полиэтилена высокого давления (ПЭВД), взятых при нормальной работе реактора и при появлении отклонений в режиме синтеза от нор­ мального. В первом случае отклонение ПТР было в пределах от 12 до —14% от его номинального значения, во втором — от 35 до —70%. Ко­ лебания, вызванные отклонением режима синтеза ПЭВД от нормаль­ ного, продолжаются несколько часов, в течение которых установка син­ теза может произвести 40—50 т материала. В показанном на рис. 1 случае дисперсия ПТР для проб, взятых в течение 2,5 ч через каждые 15 мин, равна 4,62 при математическом ожидании 1,64. Для этих же проб значения коэффициента разбухания экструдата р колеблются от 1,38 до 1,65. Прочность на растяжение проб, определенная по стандарт-

Рис. 1. Изменения ПТР при синтезе ПЭВД в нормальном режиме синтеза (а) и при отклонении от него (б).

Рис. 2. Изменение предела прочности о, коэффициента разбухания экструдата Р и пленки п при отклонении режима синтеза от нормального.

Рис. 3. Изменения наибольшей ньютоновской вязкости т]о и характерного времени релак­ сации Тг при нормальном режиме синтеза (а) и отклонении от него (б).

* Доклад, представленный на IV Всесоюзную конференцию по механике полимер­ ных и композитных материалов (Рига, октябрь 1980 г.).

ной методике, характеризуется дисперсией 108,6 при математическом ожидании 133,3.

Изменчивость физических свойств полимера не только сказывается на качестве готовых изделий, но и ставит определенные ограничения на процесс переработки. Так, предельная скорость вытягивания пленки п менялась от 102 до 160 условных единиц (дисперсия 188 при математи­ ческом ожидании 130). От пробы к пробе изменялась стойкость к рас­ трескиванию, характеризующая долговременные свойства полимера, от которых зависит гарантированный срок эксплуатации полимерных из­ делий. Как следует из рис. 1 и 2, изменение всех отраженных на рисун­ ках характеристик происходит практически синхронно — ПТР, коэффи­ циент разбухания экструдата и предельная скорость вытягивания пленки достигают максимальных значений, а предел прочности при растяже­ нии — минимальных. Из этого следует, что изменения всех перечислен­ ных характеристик, несмотря на различные состояния материала, обус­ ловлены в основном одной причиной — изменчивостью молекулярного строения.

Для управления процессом синтеза полимера необходимо оперативно определять основные характеристики его молекулярного строения. Од­ нако из-за того что приготовление растворов для исследования требует длительного времени, прямые физико-химические методы для оператив­ ного контроля молекулярного строения полимера неприменимы. Возни­ кает необходимость установить набор легко измеряемых физических характеристик, неизменность которых гарантирует относительное посто­ янство основных характеристик молекулярного строения.

В [1] предлагается модель, согласно которой достаточно трех незави­ симых характеристик вязкоупругости для определения показателей молекулярного строения; в [2—4] разработаны основы способа определе­ ния этих показателей по наибольшей ньютоновской вязкости г|о, харак­ терному времени релаксации тг и зависимости между полидисперсностыо молекулярных масс и длинноцепной разветвленностью молекул (тг — величина, обратная частоте, соответствующей максимуму частот­ ного релаксационного спектра). Из перечисленных работ следует, что при поддержании постоянных значений ri0 и тг материал имеет практи­ чески неизменное молекулярное строение.

В связи с этим была поставлена серия опытов для выявления чувст­ вительности т]о и тг к изменчивости молекулярного строения. Для 10 проб, взятых при нормальном режиме работы реактора, определяли час­ тотные зависимости составляющих комплексного модуля, по которым устанавливались г|о и тг. Пробы брались с интервалом в 15 мин, что позволило отразить изменение значений этих характеристик во время синтеза (рис. 3). Для приведенных данных диапазон изменения наиболь­ шей ньютоновской вязкости относительно среднего значения в нормаль­ ном режиме работы составлял от 66 до —15%, для характерного вре­ мени релаксации — от 172 до —66%. Для проб, взятых при нарушении нормальной работы реактора, отклонения в значениях г)о были равны от 6 до —46%, а для тг — от 82 до —49%. Диапазон относительного изме­ нения г|о и тг больше, чем ПТР и других физических характеристик по­ лиэтилена. Поэтому данные показатели наиболее пригодны для иденти­ фикации молекулярного строения полимера по его механическим свойствам. Для практического использования указанных характеристик необходимо разработать способ определения значений т)0 и тг, обладаю­ щий достаточной точностью. Наиболее целесообразно определять т]о и тг по результатам линейного периодического деформирования расплава по­ лимера. Используем для аппроксимации динамической вязкости рас­

плава уравнение [5]

Ло

Л =-

( 1 )

1 +

(о)Тг) а

В таком случае цо можно найти из результатов периодического деформи­ рования на трех частотах. Если частоты связаны соотношением

то

coi/сог= СО2/003 = т ,

(2)

2т1/1Т]/з- г]/1'П/2-'П/2Г1/з

 

1То = Г) 2'

 

/

/

/9

( 3 )

 

Ц

 

з—П т

 

а характерное время релаксации тг определяется уравнением

 

 

 

 

[ (Ло/Лз)-1

1

ОПо-тЫ-1 J

Шз ' ■

При соблюдении условия (2) необходимо оценить наиболее выгодное соотношение частот т , а также отношение одной из частот, например, частоты 02, к характерной частоте релаксации 0г= 1/тг, при которой до­ стигается наибольшая точность определения т)о и тг.

При измерении динамической вязкости г|г на трех частотах с ошиб­ кой бг вычисляется

Т1о(бь $2, бз) •

Расчет проводится по формуле (3), в которой

Vi = Л/1и(1 + 61); T|/2= TI'2H(1 + 62) I т],з==:'П/зи(1+ S3) •

(4)

г^'ш вычисляется по (1) при условии 6г = 0. Ошибки 6j варьируются в ин­ тервале |бтах —бт ш|- При всех расчетах было принято а = 0,65, что в среднем соответствует коэффициенту темпа снижения динамической вязкости полидисперсных полиэтиленов.

Вычисленная по (3) с учетом (4) величина т^о (бi, б2, б3) дает возмож­ ность определить отклонение наибольшей ньютоновской вязкости от фактического значения т)он (задаваемого при расчете r\iU) -.

АтДб!, 62, 63)

110(61, 52, б3) —Т)ои

(5)

 

TJou

Расчет проводится в интервале (бт т, бт ах), одинаковом для всех трех частот измерения с шагом б выбранным таким образом, что (бтах+ + |6min|)/6j — нечетное число. В выбранном интервале получаем для измерений на всех трех частотах [(бтах+ | бт т | )/бj]3 комбинаций возмож­ ных сочетаний ошибок.

Предполагается, что бгподчиняется нормальному распределению. Параметры распределения вычисляются на основании эксперименталь­ ных данных о точности измерения т|', которые выражают вероятность рл нахождения ошибки в интервале 16max—6min| • Построенное таким обра­ зом распределение делится на участки, длина которых соответствует шагу бj, и для каждого участка длиной 8j вычисляется вероятность на­ хождения в нем ошибки. Таким образом, для каждого значения относи­ тельной ошибки бг, определяемой с точностью до ± 6j/2, получаем соот­ ветствующую вероятность рг-.

Случайные ошибки измерения на каждой из частот не зависят от ошибок на других частотах, и поэтому вероятность любой из комбина­ ций ошибок равна p = pi (61)^2(62)Рз(бз)- В силу этого для каждого зна­ чения Дп, определяемого формулой (5), можно вычислить соответствую­ щую вероятность р(бь 62, 63) . Если ошибке 62 задать ряд фиксированных значений бгь, то в трехмерной координатной системе бь 62, An для каж­ дого фиксированного значения бг/{получим поверхность Ал/«(бi; б3) . Если в этой плоскости отобразим условие р(бг, б3) ^ р Р при б2 = б2/{= const, то

оно образует замкнутый контур на плоскости Дт,ь, который при Pi(0)p2(0)/?3(62ft) =рР вырождается в точку (рис. 4). При постепенном изменении k замкнутые контуры образуют объем, внутри которого нахо­ дятся те комбинации ошибок измерения, вероятность появления которых больше рр.

При выполнении численного расчета была принята не очень высокая точность измерения динамической вязкости — предполагалось, что ошибки измерения б* с вероятностью ^ = 0,99 находятся в интервале [0,055-=----0,055]. Шаг по б принят равным 6j = 0,01. Рассматривалось II3 комбинаций ошибок б*, из которых учитывались только имеющие веро­ ятность более 10-4. В таком случае вероятность того, что результат из­ мерения будет отброшен, равна 8,576-10_3. Это значит, что при опреде­ лении т]о и тг в заводских условиях с 10-мин циклом, допустимым счита­ ется одно измерение пониженной точности за две с половиной смены.

Из проделанных расчетов следует, что при оз2=1/тг и т = 10 ошибка определения наибольшей ньютоновской вязкости Д^ по трем частотам с вероятностью 0,9914 находится в интервале 118,5-=- —12,2% |. С увеличе­ нием базы измерения, т. е. с возрастанием числа т , величина ошибок незначительно уменьшается. Если т становится меньше 10, происходит резкое изменение характера зависимости Дл от т, и ошибки катастро­ фически возрастают. С изменением частоты со2 в области 0,1/тг<:<о2< <С1/тг точность определения г)о меняется сравнительно мало, но если со2> 1 /тг, точность резко падает с увеличением со2 (рис. 5).

Аналогичные результаты были получены при анализе точности вы­ числения характерного времени релаксации. По сравнению с наиболь­ шей ньютоновской вязкостью точность определения характерного вре­ мени релаксации во всех случаях значительно ниже, и при некоторых комбинациях ошибки б* могут достигать 60%.

Для полиэтилена низкого давления полидисперсность вычисляется из соотношения [3]

Тг6/(1+6)

(6)

Ло

где 0 — показатель степени в урав­ нении зависимости т]о от средневе-

Рис. 4. Отклонения

бп в зависимости от ошиб.ок при измерении динамической вязкости

Рис. 5.

 

на трех частотах. 62=ai; 63= a 2.

 

 

Наибольшее

отклонение 6Л в зависимости

от отклонения между частотами из­

мерения

(а) при со = 1 /тг и от отношения сог/тг (б)

при т = 1 0 .

--------- наибольшее откло­

нение 6П при измерении динамической вязкости

с точностью

± 5 % ; -----------

наиболь­

 

 

шее отклонение с вероятностью р > 1 0 -4.

 

совой молекулярной массы. При Определении С с ошибками Дп и хг вы­ ражение (6) имеет вид

тге/0+9>t (1+Лт)0/(1+0)

Ло

1 + Дт]

При заданном способе расчета т]о и хг ошибки Дл и Дт всегда одинако­ вого направления, и их максимальные значения достигаются одновре­ менно. При 0= 3,4, Дл = 0,2 и Дх=0,6 характеристика С определяется с завышением на 19,8%. Для набора проб линейного полиэтилена, харак­ теристики которого приведены в [6], вычисление С с ошибкой 19,8% вызывает отклонение при расчете и на 11,2%, а вычисленная из

tic= k[Mwf(u)]°

(7)

величина Mw отличается от истинного значения на 5,5%. Функция поли­ дисперсности в (7) принята в виде [3] f(u) =ехр ( —knu) , где мера поли­ дисперсности молекулярных масс и= (Mw/Mn) —1. Таким образом, точ­ ность определения показателей молекулярного строения по характерис­ тикам вязкоупругости расплава полиэтилена не ниже точности при их определении физико-химическими методами.

Для разветвленных полимеров точность расчета показателей молеку­ лярного строения зависит от ширины распределения молекулярных масс. Для широких распределений, в самом неблагоприятном случае, Mw рас­ считывалось с ошибкой 12%, характеристика длинноцепной разветвлен­ ности <g)w — 5,7%, а индекс полидисперсности — 28%. Точность этих результатов также не уступает точности физико-химических методов при исследовании молекулярного строения разветвленных полимеров [7, 8].

Выводы. 1. Выпускаемые промышленностью полиэтилены обладают значительной изменчивостью механических свойств.

2.Характеристики линейной вязкоупругости расплава полиэтилена достаточно чувствительны к изменениям молекулярного строения при его синтезе.

3.Из результатов линейного периодического деформирования на трех частотах можно с достаточной для вычисления показателей моле­ кулярного строения точностью определить характеристики вязкоупру­

гости Г]о и Хт- 4. Лучшие результаты при определении характеристик вязкоупру­

гости по результатам периодического деформирования на трех частотах получаются, если средняя частота близка к обратному характерному времени релаксации расплава, низкая частота в 10—20 раз меньше, а высокая во столько же раз больше средней.

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

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

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

1.Бриедис И. П., Бурлий В. Н. Определение вязкоупругих характеристик расплава по молекулярно-массовому распределению. — В кн.: Структура и свойства полимерных материалов. Рига, 1979, с. 68—77.

2.Бриедис И. П., Файтельсон Л. А. Зависимость реологических свойств расплава

полиэтилена от его молекулярного строения. — В кн.: Реология полимеров и дисперс­ ных систем и реофизика. Т. 1. Минск, 1975, с. 78—84.

3. Бриедис И. П„ Файтельсон. Л. А. Реология и молекулярное строение расплавов полиэтилена. 2. Влияние молекулярного строения на вязкоупругие характеристики. — Механика полимеров, 1976, № 1, с. 120— 127.

4.Брауэр Э., Бриедис И. П., Бухгалтер В. И., Сульженко Л. Л., Файтельсон Л. А., Фидлер П. Реология и молекулярное строение полиэтилена. 4. Экспериментальное иссле­ дование. — Механика полимеров, 1977, № 2, с. 283—293.

5.Бриедис И. П., Файтельсон Л. А. Реология и молекулярное строение расплавов полиэтилена. 1. Идентификация расплавов по механическим свойствам. — Механика полимеров, 1975, № 3, с. 523—532.

6. Schroff R., Shida М. Effect

of long chain branching of the relation between

steady-flow and dynamic-viscosity of

polyethylene melts. — J. Polymer Sci., A-2, 1970

№ 8, p. 1917— 1925.

 

7. Strazielle C., Benoit H. Molecular characterisation of commercial polymers. — Pure and Applied Chem., 1971, vol. 26, p. 451—479.

8. Adams H. G., Ahad E., Chang H. S., Dabis В. B., French D. H., Hyer H. J., Lan R. D., Simkins R. J., Stuchbyry J. E., Tremblay M. A cooperative molecular weight

distribution test. — J. Appl. Polym. Sci., 1973, vol. 17,

№ 1, p. 269—283.

Институт механики полимеров

Поступило в редакцию 19.03.80

АН Латвийской ССР, Рига

 

Всесоюзный центральный научно-исследовательский

 

институт комплексной автоматизации, Москва

 

Охтинское научно-производственное объединение

 

«Пластполимер»

 

УДК 539.376.001:678

А. Ф. Крегерс

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

Для решения краевых задач часто необходимо иметь связь между напряжениями и деформациями одновременно в двух разных формах: е= е(сг ) и сг = а(е), что в области физической нелинейности иногда соз­ дает определенные трудности. Рассмотрим некоторые специальные усло­ вия, при которых можно получить взаимно обращаемые нелинейные со­ отношения материала.

Вначале остановимся на зависимостях нелинейного упругого мате­ риала, так как это является частным случаем связи вязкоупругого материала, например, при о = const и ^->0 или t-*-оо. В подавляющем большинстве случаев эта связь представлена в виде тензорного поли­

нома [1]

 

8ij QijlilGhl4"^ijftimnCT/i/CTmn&ijhlmnopGhlOmnOop

где i, /,..., р= 1,2, 3. Положительной стороной такого

полинома явля­

ется простота математической структуры и возможность двумя-тремя его членами (не обязательно взятыми подряд) обеспечить хорошую аппроксимацию экспериментальных данных. Особенно отметим качест­ венное соответствие характера кривой е ~ а уравнению (1), т. е. рост производной \de/do\ с увеличением численных значений напряжений. Постепенное уменьшение касательных и секущих модулей упругости при росте растягивающих напряжений свойственно в большей или меньшей мере всем жестким конструктивным материалам (для биоматериалов наблюдается обратная закономерность [2]). Отрицательной стороной со­ отношения (1) является катастрофически быстрый рост количества не­ зависимых параметров материала [1] с ростом степени нелинейности и невозможность выразить аналитически из (1) напряжения через дефор­ мации в конечной форме.

Связь а ~ е также можно представить в виде полинома [3, 4]

( 2)

Структура выражения (2) такова, что им удобно описывать деформативные свойства биоматериалов, но не жестких конструктивных материа­ лов. Чтобы в случае растяжения получить кривую о ~ е с убывающим касательным модулем, необходимо сохранить много членов в (2), при­ чем кривая будет иметь несколько точек перегиба. Кроме того, кривая сг~е может иметь экстремумы и даже менять знак при е > 0, например, при а = ае —6е2. В результате можно говорить лишь о применимости уравнения (2) в каких-то пределах по е.

Изложенные выше недостатки уравнений (1) и (2) привели к по­ пытке найти такие нелинейные формы связи в классе простых функций, которые были бы обратимы и имели бы сравнительно небольшое коли­ чество параметров (рис. 1). Понятно, что поставленная цель приведет к

Доклад, представленный на IV Всесоюзную конференцию по механике полимер­ ных н композитных материалов (Рига, октябрь 1980 г.).

1.5

Рис.

1. Некоторые элементарные

функции, используемые

для

ап­

проксимации

кривых

а ~ е :

1

j/=Arsh(x);

2 —

г/= 1п(1+х);

3

i/= arctg(x); 4

у = th(x);

5

у = 1 —ехр ( —х );

6

х =

 

= у + 4 у 3; 7 у =

2

 

 

 

—arctg(x).

менее общим зависимостям, чем уравне­ ния (1) и (2), но создаст определенные преимущества при решении таких задач, в которых одновременно необходимо поль­ зоваться аналитическими зависимостями

типа сг = сг(е) и е= е(а).

Вопросам применения экспоненциаль­ ной функции для описания физической нелинейности посвящены работы [2, 5], а функции гиперболического синуса — в [6, 7]. Вопросы возможности обращения физически нелинейных уравнений состоя­

ния в

области упругости обсуждаются

В [5].

базе экспоненциальной функции

На

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

ст = 6[1 —ехр( —c|e|)]sign(e),

(3)

из которой можно выразить деформацию

In ( 1---- ^ L )s ig n (a ); o < b .

(4)

Изотропный материал. Для обобщения зависимостей (3) и (4) на случай произвольного объемного напряженного состояния нелинейно­ упругого изотропного материала воспользуемся следующими обозначе­ ниями:

Ро—llWал —~jЬJ \а2

Ь212а\

I\e= Eij&ijt

 

Ре У ^ е л —

 

В 212г‘,

l i e =

GijEjil

(5)

Lij(o)

dWan '

2Ь{ I

2 b 2Gij>

I —Gijbijt

 

 

j

 

 

 

 

 

Lij(e)

dWzл

2 S ] / i eS ij ”f- 2 B 2E{ji

I2a OijOji,

 

Оец

 

где /,/= 1,2,3; Wo11, WEn — квадратичный потенциал линейно-упругого материала, выраженный через инварианты тензора напряжений или де­ формаций соответственно.

Представим выражение напряжений в следующей форме:

* i = - ^ [ 1- e x p ( - B ,A ) ] .

(6)

В случае сдвига

а12 = sign (е12)1В2/2[1 - exp ( - В ^2 В 21ei21)],

(7)

что по структуре соответствует уравнению (3). Из (7) получаем, что де­ формация сдвига равна

 

BQ]/2B2

[ 14-Ы)

 

£12 =

sign(ai2) In

 

( 8)

С учетом (4), (5) и (8) для произвольного напряженного состояния компоненты тензора деформаций имеем

е « = -

1п(1-Р„); 0<Р„<1.

(9)

 

ZP0

 

 

Связь между параметрами В0, В\,

В2 уравнения (6)

и параметрами

Ьо, Ь\, Ь2 уравнения (9)

установим из рассмотрения двух видов напря­

женных состояний. Согласно (9) в случае сдвига

 

ei2= - Ь 0У62/21п(1—У262| 0121) sign (ai2).

(10)

Из сопоставления параметров уравнений (10) и (8) находим, что

 

Ь2=\/В2 и

Ьо=\/В0.

(11)

При нагружении гидростатическим давлением ап = 022 = сгзз = сгср бла­

годаря

изотропии

материала ец = е22 = езз= Еср. Из (6) выпишем

выра­

жения для ап, 022 и а3з, откУДа

 

 

 

или

оСр = sign (еср)1В{+ В2/3[1 —exp (—В^9В\ + ЗВ21еСр | ) ]

(12)

 

 

 

 

 

 

Вер —

s ig n (a c p )

 

3

(13)

 

Bd]/9Bi + ЗВ2

 

3Bi + В2

 

 

 

 

Из (9)

после определения ец> егг, езз при гидростатическом давлении

находим

 

 

 

 

 

еср= —Ь0sign(ocp)ybi + 62/З ln[l-Y9fci —3&2|а ср|].

(14)

Сопоставляя (13) и (14),сучетом

(11) имеем

 

 

 

 

или

Ь\

(15)

 

bi = -

(ЗВ\ + В2)В2

В\ = — (ЗЬ{ + Ь2) ь2

Изучение уравнений (6) и (9) показывает, что коэффициент Пуас­ сона при одноосном растяжении не зависит от величины напряжения и равняется

 

v= v0 =

Б22 (CTi 1)

Ь\

= const.

(16)

 

fin (on)

b\ + Ь2

Модули упругости при Oij->-0 выражаются как

 

 

 

Е0=-

1

G0=- 2b b

(17)

 

bo(b\-\-b2)

 

 

 

 

0^2

 

Несжимаемый изотропный материал получим, если вместо (5) примем

------

д

 

-------

д

 

Ра — УЬ12аО\

dsn ш(Ро2) —2bSij\

Pt = ~jBI2eD\

^<j(e)= ~Tоец

{Pt2) = 25e>j">

 

 

?2<JD = SijSji',

l 2eD = £ij&ji'i

 

 

i вместо компонент e,j и а,-,- будем пользоваться девиаторными компонентами е,-, и s , j .

1ри этом Ь=\/В и Ь0= \/В 0.

С учетом (5) уравнение (6) можно привести к виду

 

дР

(18)

G ij= [\ - e x p { - B 0P e ) ] ~ -

UBij

 

■+■ BQP)E—1]■