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

Теоретические основы переработки полимеров (механика процессов)

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

Сопоставляя экспериментальные данные по температурной за­ висимости удельного объема с результатами расчета по уравнению (IV. 12) для Р = 0, можно определить значения V* и Г*

Величина Р* определяется путем сопоставления комплекса

вида PW/T с его теоретическим значением, рассчитываемым из уравнения (IV. 12) [при этом у определяется численным решением уравнения (IV. 13)]. Отношение экспериментального и теоретиче­ ского значений (Р/Р) дает величину Р*. На практике обычно вы­ числяют несколько значений Р* при разных температурах и берут среднее значение.

Для цепи бесконечной длины нормирующие параметры оказы­ ваются связаны соотношением:

P*V*/T* = R (C/ S ) / M Q

(IV. 14)

где Мо — молекулярная масса среднего сегмента.

Если принять, что для всех полимеров cjs = 1/3, то уравнение (IV. 14) принимает вид:

(P*V*/T*) М0 = 27,7

(IV. 15)

Значения характерных параметров для ряда полимеров при­ ведены ниже:

 

 

 

МПа

V* •10\ см3/кг

г, к

гг

к (P*V*/r*)Af0

 

 

 

Р \

 

 

 

 

Поли-о-метилсти-

745,8

0,9762

 

12 740

 

404

33,8

рол

[3] . .

 

 

Полистирол

[6]

763,8

0,9625

 

12 700

 

362

30,1

Полистирол

[3] .

745,3

0,9598

 

12 680

 

365

29,4

Полиметилметакри­

930,3

0,8350

 

И 890

376

32,7

лат

[6]

. . .

 

Поливинилхлорид [6]

1034,5

0,7105

 

11 320

349

20,3

Изотактический по­

 

 

 

 

 

 

 

лиметилметакри­

1008,8

0,8160

 

11 170

320

36,9

лат

[4]

 

 

Температурная

зависимость

В(Т)

может

быть вычислена из

теоретического выражения:

 

 

 

 

 

В = -

(0P/dV)T Voс -

Р

 

 

 

 

 

(IV. 16)

где

 

 

 

 

 

 

 

 

 

S/QC =

- (dP/dV)T/(d2PldV2)T

 

 

 

 

(IV. 17)

Легко показать, что с слабо зависит от Т. Однако можно при­ нять, что с = 0,0894 и рассчитать температурную зависимость В. Результаты такого расчета представлены на рис. IV. 1. Сопостав-

Рис. IV. 1. Зависимость нормированного
коэффициента В (Т) от нормированной температуры Т. Сплошная линия рассчи­ тана при c=0,0S94; ф — полистирол [6]; А — полиметилметакрилат 161; X — поли­ винилхлорид [61; О — полистирол [31; А поли-о-метилстирол [3].

ление теоретических (сплошная кривая) и экспериментальных данных (отдельные точки) сви­ детельствует об их удовлетвори­ тельном совпадении (ошибка ме­ нее 5%).

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

Термодинамическое уравнение состояния. Классическая термо­ динамика, используя представле-

ние о связи давления с измене! ем внутренней энергии и энтропии, дает следующее выражение д i термодинамического уравнения состояния:

 

(IV. 18)

Поскольку (dS/dW)T = (дР/дТ)у, уравнение (IV. 18)

можно запи­

сать в виде:

 

Р +

(IV. 19)

Используя полученное Франком [7] выражение для внутренней

энергии

 

U = - a / V n

(IV 20)

и известное из статической физики соотношение

 

PV = NkT = RT

(IV. 21)

(где N — число Авогадро, k — постоянная Больцмана, R — универ­ сальная газовая постоянная), получим:

Р +

па

(IV. 22)

V"+i

 

Если параметр п = 0, то уравнение (IV. 22) превращается в известное уравнение состояния для идеального газа (IV. 21). Если п = 1, то уравнение (IV. 22) превращается в уравнение Ван-дер- Ваальса:

а

RT

(IV. 23)

Р + уТ

V - Ь

Уравнение состояния Спенсера — Джилмора. Пытаясь обобщить извест­ ные экспериментальные данные по сжимаемости и температурному расши­ рению расплавов полимеров, Спенсер и Джилмор [9] установили, что это уда­ ется сделать с помощью модифициро­ ванного уравнения Ван-дер-Ваальса:

+ я) СV - Ь) = RT/M

(IV. 24)

Рис. IV. 2. Зависимость термиче­ ского коэффициента линейного рас­ ширения от температуры:

/ — полиэтилен

высокой плотности;

2—полистирол;

3— полиэтилен низ­

кой плотности;

4—полиамид [6, 10].

Здесь V — удельный объем полимера; М — молекулярная масса структурной единицы, об­ условливающей межмолекулярное взаимодей­ ствие; л и b — константы (по своему физиче­ скому смыслу л — это значение внутреннего давления, а b — объем, занимаемый собственно молекулами полимера); л, b и М опреде­ ляются экспериментально.

Если

продифференцировать уравнение (IV. 24)

по переменной

V, то после несложных преобразований получим:

 

л Ум +

Р Ум = д / - RT (дР/дУ)т

(IV. 25)

Очевидно, что величины я и М можно рассчитать, если изве­ стны значения (дР/дУ)т и Р. Константу Ь нетрудно вычислить по значению удельного объема при атмосферном давлении и извест­ ной температуре. Константы уравнения состояния для ряда наи­ более распространенных полимеров приведены ниже.

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

Р = - J (дР/дР)т

(IV. 26)

К = - ( д р /д Т ) р

(IV. 27)

Р

 

Значение р для всех материалов положительно. Для расплавов полимеров значение К также всегда положительно.

Дифференцируя уравнение (IV. 10), можно получить следую­ щие выражения для р и К:

р = (Р + я Г 1[ 1 + ^ - ( Я + я)]

(IV 28)

Я = [г + ^ - ( Р + я)]

(IV 29)

Интересно сопоставить значения коэффициентов сжимаемости и температурного расширения, определенные непосредственно из эк­

спериментальных данных, с результатами, полученными при ис­ пользовании уравнений (IV. 28) и (IV. 29).

Ниже приведены константы уравнения состояния:

 

 

 

 

 

м

я -К )-8, Па

6-103, м3/кг

Полиэтилен

низкой

плот­

28,1

3,20

0,875

ности [9 ].......................

высокой

плот­

Полиэтилен

281,0

6,77

1,110

ности [10] . .

 

Полистирол

[ 1 1 ] ...............

[11]

104,0

1,80

0,822

Полиметилметакрилат

100,0

2,10

0,734

Полиамид

6,

10

[11]

. .

2,7

1,08

0,738

Этилцеллюлоза

[9]

60,5

2,37

0,720

Ацетобутират

целлюлозы

54,4

2,80

0,688

[9]

 

 

 

 

По данным [И], коэффициент температурного расширения для полиэтилена низкой плотности при атмосферном давлении в интер­ вале температур 493—473 К равен 6,75-10_4К_1. Расчет по урав­ нению (IV. 28) при Т = 393 К дает значение X = 7,35- 10-4К-1, т. е. примерно на 10% больше экспериментального значения.

На рис. IV. 2 приведена температурная зависимость термиче­ ского коэффициента линейного расширения при атмосферном дав­ лении для четырех различных полимеров, рассчитанная по фор­ муле (IV. 29). Из этого рисунка следует, что применение линейных формул, описывающих зависимость удельного объема от темпера­ туры, является в ряде случаев довольно грубым приближением.

IV.3. КИНЕТИКА ПРОЦЕССОВ УСАДКИ И РАСШИРЕНИЯ

Уравнение (IV. 10) описывает равновесные изменения удельного объема, соответствующие достаточно медленному разогреву — ох­ лаждению или сжатию — разгрузке расплава. При быстром разо­ греве или охлаждении расплава, так же как и при быстром сжа­ тии или разгрузке, изменения удельного объема не будут успевать за изменениями температуры (или давления). Это явление отра­ жено на рис. IV. 3, на котором представлена кинетика изменения удельного объема расплава при резком повышении его темпе­ ратуры.

Результаты экспериментальных исследований кинетики темпе­ ратурного расширения расплавов позволяют предложить следую­ щее уравнение, описывающее временную зависимость удельного объема:

V =

а _ а + ЬТ - -

а + ЬТг - V*) ехр ( - р АТ)

(IV. 30)

где

Ух — равновесный

объем, соответствующий температуре

7Y, АТ = Т — 7\.

Константа а, определяющая до­ лю приращения удельного объема, зависящую от скорости изменения температуры, находится из выра­ жения:

a = r ( b - b ') / K v

(IV. 31)

Рис. IV. 3. Кинетика изменения удель­ ного объема расплава полимера при резком повышении температуры.

где

г — скорость изменения температуры

(г =

dTJdt); Ъ' — коэффициент, характе­

ризующий мгновенное, не зависящее от вре­ мени изменение удельного объема ДУ< (см.

рис. IV. 3), прямо

пропорциональное изме­

нению температуры;

 

или

(dWldt)T = bf

(IV. 32)

Kv — константа, характеризующая

скорость

изменения удельного объема;

 

(dV/dt)T = Kv (V0- V )

(IV. 33)

Для определения константы Kv используют данные, получен­ ные при наблюдении за кинетикой расширения полимера при скач­ кообразном изменении температуры. Легко показать, что при этом величина К» связана с полупериодом времени расширения %зави­ симостью [11, с. 140]:

Kv = 1п2/х

(IV. 34)

Константа f$, характеризующая влияние температуры на изме­ нение удельного объема, определяется выражением:

Р = Kv/r

(IV. 35)

В случае бесконечно медленного изменения температуры урав­

нение (IV. 30)

принимает вид:

Ve = a + bT

(IV. 35а)

Видно, что константы а и b характеризуют равновесное значе­

ние удельного объема Vc и легко рассчитываются из уравнения со­

стояния

(IV. 24).

 

 

В том случае, если имеет место изотермическое расширение,

уравнение (IV. 16) сводится к виду:

 

Ve- v

= ехР ( - V )

(IV. 36)

y

~ v

е

' I

 

 

 

Константа k v зависит от температуры. Поэтому при

построе­

нии кривых изменения объема необходимо вести расчет методом последовательных приближений, принимая каждый температурный интервал равным примерно 20 К. Типичные значения константы скорости kv для изотермического расширения полистирола приве­ дены ниже [12]:

Т, к •

310

330 350 370 390 410 430 450 470 490

ftv, ч"1

0,068

0,26 0,8 2,3 6,0 13,2 27,5 57 105 190

Рис. IV. 4. Неравновесные кривые из­ менения удельного объема полистирола при разогреве (или охлаждении) при

различных

скоростях

охлаждения

(/, 2, 3) и нагревания (4—7):

1

м и н у с оо;

2 — м и н у с 1,0 ° С / с ; 3 — м и ­

н у с

0,1 ; 4 оо; 5 — 1,0; 5 — 0,1;

7 - 0 , 0 1 , С / с .

Значения

остальных

констант

для

полистирола:

а = 0,845 X

X Ю3

см3/кг,

b =

0,375

см3/(кг. К),

Ь' = 0,25 см3/(кг-К).

 

 

На

рис. IV. 4 приведены кривые

изменения

удельного объема при

равной

скорости разогрева и охлаж­

дения,

рассчитанные

по

уравнению

(IV. 30). Заметим,

что

увеличение

скорости

нагрева

 

(или

охлажде­

ния) приводит к повышению темпе­ ратуры, при которой соответствую­ щая кривая отклоняется от равно­ весной кривой (жирная сплошная линия на рис. IV. 4).

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

Допустим, что лист из полистирола получают методом экструзии в водяную ванну. Затем из листа вырубают на прессе изделия, причем требования к точ­ ности размеров этих изделий таковы, что максимально допустимое изменение линейных размеров за счет усадки не должно превышать 0,33%.

Пусть температура расплава на выходе из экструдера составляет 500 К, а температура охлаждающей воды 300 К. Определим, сколько времени надо вы­ держать лист при 300 К, чтобы остаточная усадка не превышала 0,33%.

Когда горячий лист попадает в холодную воду, скорость охлаждения бы­

вает максимальной; без большой ошибки можно

принять,

что

г — —с». При

300 К первоначальный удельный объем составляет 0,9825

10~3

м3/кг.

Посте­

пенно удельный объем уменьшается и достигает

равновесного

значения

Ve =

= 0,9575

Ю-3 м3/кг.

 

 

 

 

При допустимой усадке 0,33% лист необходимо выдерживать в воде до тех

пор, пока

удельный объем не станет равен 0,9672

10_3 м3/кг. Это вытекает из

условия, что если максимальное допустимое изменение линейного размера равно а, то изменение объема равно

Vi/Ve =

(1 + а)3 »

1 + За

 

(IV. 37)

При

изотермической усадке следует использовать уравнение

(IV. 36);

V* -

V

0,9575

— 0,9672

= 0,388 =■ exp (— &Y/)

 

Ve-

V,

0,9575

— 0,9822

 

При 300 К значения k v , рассчитанное методом экстраполяции по приведен­ ным выше данным, составляет 0,035 ч-1. При этом время, необходимое для усад­ ки до удельного объема V = 0,9672 составит 30 ч. Если повысить темпера­ туру выдержки листа до 350 К, то изменение удельного объема будет описы­

ваться прямой аб (см. рис. IV. 4) В точке

б удельный объем становится рав­

ным 0,98 м3/кг. При понижении температуры

материала до

300 К его

удельный

объем определится точкой с.

ky увеличится

до 0,84-1,

время от­

Поскольку при этом значение константы

жига уменьшится до 2 ч.

 

 

 

IV.4. ЭНЕРГЕТИЧЕСКИЙ БАЛАНС ПРОЦЕССА УСТАНОВИВШЕГОСЯ ТЕЧЕНИЯ

Рассмотрим систему, схема которой приведена на рис. IV. 5. По­ скольку течение установившееся, то внутри системы не происходит накопления массы или энергии.

Первое условие математически формулируется так:

-щ- (pAU) = 0

(IV. 38)

где / — направление течения; А — площадь поперечного сечения

потока; 0

средняя скорость течения жидкости.

 

Второе условие выражается уравнением:

 

bE + bEb + H E p ^ Q - ^ W

(IV. 39)

где Ек и ЕР представляют собой соответственно кинетическую и потенциальную энергию жидкости.

Изменение потенциальной энергии единицы массы жидкости при подъеме на высоту Z составляет:

A£p = gZ

(IV 39а)

где g— ускорение свободного падения.

 

Приращение кинетической энергии при возрастании средней скорости от Ua до Оъ определяется выражением:

(IV. 40)

где k — коэффициент, зависящий от распределения скоростей течения в потоке.

Если рассматривается течение в круглой трубе радиусом /?, то средняя кинетическая энергия равна:

Ek = - ^ - \ n 3rdr

 

(IV. 41)

о

 

 

где G — Мас£°&ый расход.

Следовательно,

коэффициент k определяется из соотношения:

--

-

(IV. 42)

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

ТЕПЛОПРОВОДНОСТЬ

Основной задачей теории теплопроводности является установле­ ние распределения температур внутри тела. Если распределение температур не зависит от времени, то задача теплопроводности является стационарной; если распределение температур зависит от времени, то задача становится нестационарной.

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

q = kVT

(IV. 48)

где q — количество тепла, проходящего через единичную поверхность, перпен­ дикулярную направлению теплового потока; k — коэффициент теплопроводности.

Полагая в уравнении энергетического баланса (III. 4) V = О, получим:

9cv ( ^ f ) ------

(Vf)

(IV. 49)

Уравнение (IV. 49) представляет собой уравнение теплопро­ водности для изотропного твердого тела.

Если внутри изотропного тела имеется источник тепла, то урав­ нение (IV. 49) необходимо дополнить членом, учитывающим теп­ ловыделение:

dT/dt = aV2T + (G/pcy)

 

(IV. 50)

где а = k/pcp — коэффициент

температуропроводности [замена

с\ на сР в

уравнении

(IV. 50)

возможна

для несжимаемых твердых тел];

V2 — оператор

Лапласа в прямоугольной системе координат

 

у2Т - & Т

д2Т

д2Т

 

(IV. 50а)

дх2

ду2

дг2

 

 

 

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

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

Следовательно, распределение температур носит установившийся характер и определяется интегрированием дифференциального уравнения

aV2T + G/pcp = О

(IV. 51)

Нестационарная теплопроводность. В большинстве случаев в реальных процессах переработки -приходится иметь дело с неста­ ционарным режимом теплопроводности, когда полимер подвер­ гают нагреву или охлаждению (например, охлаждение в форме отлитого изделия). Теоретические исследования процесса неста­ ционарной теплопроводности представляют собой обширный раз­ дел математической физики. Решения, получаемые в результате интегрирования уравнения (IV. 50), представляют собой функции времени и пространственных координат, удовлетворяющие на­ чальным и граничным условиям. Различают четыре рода гра­ ничных условий [13, с. 24].

Условия первого рода: задано распределение температур на поверхности, которое может либо быть постоянным, либо зави­ сеть от времени; в простейшем случае, если положение границ определяется одним числом (например, расстоянием L), такие гра­

ничные

условия математически определяются выражением

вида:

Т (/ L) =

Т (L)

(IV 52)

Условия второго рода: задана плотность теплового потока для

каждой точки поверхности тела как функция времени:

 

Яп (0 = / (О

(IV. 53)

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

а (/, L) - Tf] = — X ~ (/. L)

(IV 54)

Условия четвертого рода соответствуют теплообмену тела с окружающей средой по закону теплопроводности или теплооб­ мену системы тел, находящихся в тепловом контакте (темпера­ тура соприкасающихся поверхностей одинакова):

Tn (t) = Tc (t)

(IV. 55)

- Л (дТ/дп)п = — X {дТ/дп)с

(IV. 55а)

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]