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

книги из ГПНТБ / Рогов И.А. Физические методы обработки пищевых продуктов

.pdf
Скачиваний:
15
Добавлен:
27.10.2023
Размер:
25.18 Mб
Скачать

 

ренные при

напряжениях,

 

превышающих

предельное

 

напряжение сдвига, опре­

 

делены

с ошибкой

до

±

 

(5—10%)

и пригодны

как

 

для

качественной

оценки

 

технологического

 

процес­

 

са, так и для расчета дви­

 

жения

продуктов

в рабо­

 

чих органах машин и аппа­

 

ратов. Поэтому таким из­

 

мерениям

следует

отдать

 

предпочтение.

 

 

 

 

 

 

Для исследования СМС

 

при

повышенных

давле­

 

ниях во МТИММПе на базе

 

вискозиметра

РВ-4

(рис.

 

6, позиции

1,

2,

3)

скон­

 

струирован новый измери­

 

тельный

узел

 

(позиции

 

410)\

стакан— статор

 

неподвижен,

 

вращается

 

рифленый

ротор,

который

 

выточен вместе

с

хвосто­

 

виком

из одной заготовки.

 

Отсутствие резьбового сое­

 

динения

ротор — хвосто­

 

вик

позволило

свести

ма­

 

ксимальную величину

би­

 

ения

на

конце ротора до

 

0,01 мм. Давление в стака­

 

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

 

вания

поршня,

который

 

снабжен упругой

мембра­

 

ной и манометром

для из­

 

мерения давления.

В рабо­

 

чем зазоре

давление изме­

 

ряется

тензометрическими

 

датчиками. Ротор

приво­

 

дится

во

вращение

от

 

падающих

 

грузов

 

через1

 

1 — станина;

2 — шкивы;

3 — бара­

 

бан; 4 — хвостовик

ротора;

5 — ро­

 

тор;

6 — тензодатчики

 

давления;

 

7 — неподвижный

стакан;

 

S — пор­

вышенном давлении:

шень;

9 — мембрана

с

манометром;

10 —■контрольный

краннк.

 

 

 

30

блоки и барабан. Влияние трения в уплотнениях учитывают при тарировке на разных давлениях и частотах вращения ротора.

Размеры измерительных роторов могут влиять на показания прибора 138]. В литературе пока нет обобщающих сведений по этому вопросу [161, однако некоторые рекомендации по выбору оптимальных размеров рабочих органов приборов можно дать [36]. Оптимальной считается на­ именьшая величина кольцевого рабочего зазора, при которой сдвиг бу­ дет распространяться на всю его толщину, но при равномерном распре­ делении давления по высоте зазора. Последнее уменьшает предпосылки проявления эффекта Вайссенберга [161. Толщина зазора должна быть на порядок выше размера частиц, чтобы исключить механическое «закли­ нивание» при измерениях и иметь равномерно изменяющийся градиент скорости. Для полусфероцилиндрических вискозиметров типа РВ-8 с уменьшением величины зазора будет уменьшаться ошибка от возмущения потока вязкой жидкости, создаваемых стыками цилиндрической и сфериче­ ской частей ротора прибора [16]. Например, для тонкоизмельченных видов мясного фарша оптимальное отношение наружного радиуса ротора RB к внутреннему радиусу стакана /?„ лежит в области 0,8—0,9 прирадиусё ротора 0,016—0,018 м [36].

Деформационное поведение реальных тел может описываться различ­ ными математическими моделями. Поэтому при разработке теории при­ бора исходная модель 1137], выбранная для вывода интегральных расчет­ ных формул, должна возможно более точно описывать существенные сдви­ говые свойства тела в данной области напряжений и градиентов скорости. Следует отметить, что в паспортах к приборам РВ-8 и РВ-4 для расчета сдвиговых свойств приведены формулы, основанные на модели Бингама. Эта модель, как показал 10. А. Мачихин 181], применима к пралиновым конфетным массам. Она была использована также довольно приближенно для обобщений экспериметальных данных по мясным фаршам [38,106].

Измерительный ротор соприкасается с продуктом не только боковой поверхностью, но и торцом. Сопротивление на торце может быть учтено следующим образом. Измерения проводят при двух-трех высотах продук­ та в рабочем зазоре. Экстраполяция прямой на графике «крутящий мо­ мент-высота» до пересечения с осью абсцисс дает в отрицательной области отрезок эквивалентной длины Лэкв, который нужно прибавить к общей высоте продукта. Эквивалентная длина показывает высоту цилиндричес­ кой части ротора, на которой сопротивления будут такие же, как на тор­ це. Иногда, например, у вискозиметра «Реотест», плоский торец заменен вогнутой «воздушной подушкой», что оправдало себя при работе с вяз­ кими жидкостями. Для пластично-вязких систем такая форма торца не­ приемлема. При исследовании подобных систем в зоне плоского торца воз­ никают значительные нормальные напряжения, которые обусловливают выход продукта из зазора и наматывание его на внутренний цилиндр [16]. Сжатый воздух в подушке резко усиливает эти эффекты. Кроме наз­ ванных способов, используют сдвоенные роторы, роторы колокольного типа и др. Для пластично-вязких систем, на наш взгляд, целесообразно использовать вискозиметры с сочетанием измерительных поверхностей различной формы, например коницилиндрические или полусфероцилин­ дрические. Ротационные вискозиметры типа РВ системы Воларовича имеют сферический торец в сочетании с цилиндрической поверхностью

(см. рис. 5 и 6).

Величины крутящих моментов, возникающих на измери­ тельной поверхности ротора с радиусом R q [16, 106, 133], имеют следующие значения:

31

цилиндрическая поверхность

Mu = 2it/?*/i0;

(1—29)

полусферическая поверхность

**Я3В

Мпс= —

0;

(1-30)

плоская поверхность торца

= Y

(1-31)

где 0 — напряжение сдвига на поверхности.

Если форма измерительной поверхности составная,то общий крутящий момент М, измеряемый на валу ротора, равен сумме моментов на обеих поверхностях:

плоскоцилиндрическая

М = 2«Я? 0 (й +

;

(1-32)

полусфероцилиндрическая

М = 2r.R; 0 ^А +

.

(1-33)

По этим формулам эквивалентная высота для плоскоцилинд­ рического ротора составляет ^экв — Лв3 , для полусфероцилинд-

рического — /гэкв = -2—2-. Тогда комбинированные измеритель-

4

ные органы можно считать цилиндрическими, но высоту про­ дукта в^зазореопределять по приведенной высоте hnp /г + /гэк„. Общий крутящий момент определяется по углу закручивания торзиона (с предварительной тарировкой) либо, если привод ротора реализуется от падающих грузов (см. рис.5), по моменту на шкиве, т. е. M = mRg,

где т — общая масса грузов, за вычетом массы (точнее силы) на трение в подшипниках;

R — радиус шкива.

Формулы (I—32) и (I—33) можно использовать для

расчета

величины ПНС:

Мв

 

 

во

 

(1-34)

2*«!апр

 

 

где Мв— момент, при котором

начинается

вращение ротора;

обычно

его определяют на пересечении кривой течения с соответствую­

щей координатной осью методом экстраполяции до нулевого

значения градиента

скорости.

 

 

32

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

Толщина градиентного слоя определяет величину градиента скорости. Если обозначить толщину градиентного слоя Дг,

причем в общем случае Дг

RltRB, то, зная окружную ско­

рость

w (м/с)

боковой

поверхности ротора

 

 

 

w =

<и/?в == 2~NRB,

(1—35)

где со

— угловая

скорость;

ротора,

 

N — частота

вращения

 

градиент скорости ё(с-1) приближенно можно представить в виде линейной зависимости (по Д. С. Великовскому [66] — средний градиент):

W

(1-36)

е = — ■

Дг

Из формулы (I—36) видно, что толщина слоя определяет величину градиента. Поэтому в расчетных формулах под наруж­ ным радиусом Rn следовало бы понимать не геометрический ра­

диус стакана, а радиус градиентного слоя, т. е. RB — Дв+ Дг. Однако определить величину Дг чрезвычайно сложно, в связи с

чем названный метод расчета (по RH) не получил распростра­ нения.

Крутящий момент сопротивления, возникающий между двумя коаксиальными цилиндрами, определяется по формуле Маргулиса [14, 16, 29]:

Мц = т]-2тс/^ /г

2to

(1-37)

Момент между двумя полусферами вычисляется по формуле Сокольского [16, 29]:

2—381

33

* 2к „

8<i>

^ПС —Ч'

(1-38)

 

1— E»\r

 

r J -

Сопоставим последние два уравнения с выражениями (I—29) и (I—30), в которых напряжение сдвига заменено произведе­

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

для коаксиальных цилиндров

4kN

(1-39)

для коаксиальных полусфер

16У

(1-40)

Rn_

где N — частота вращения ротора, об/с.

Представляет интерес совместное решение уравнений (I—39)

и (1—40):

4

(1-41)

1 — л-2 1 — х3

В результате решения получается один подходящий корень

О

х = ——= 0,675, показывающий, что при этом соотношении

Rh

радиусов величина градиента скорости в цилиндрическом и сферическом зазоре одинакова.

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

dw

d or)

г

йш

+ ш = г

da>

(1-42)

£

—-------- =

dr

dr

 

dr

 

dr

 

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

Интегрирование выражения (1—42) позволяет получить логарифмический закон для вычисления градиента скорости:

2т:

£

Rn

N.

(1-43)

2,3 lg

 

 

 

R0

 

 

34

Д. С. Великовский и др. [66] предложили строить кривую течения как функцию истинного градиента скорости на стенке вискозиметра. Этот градиент определяется по зависимости

Е

(1-44)

где

dlgO_

— величина,

пропорциональная индексу течения [см.

 

dlgw

уравнение

(1-24)].

 

 

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

Методика проведения экспериментов и их обобщения для различных ротационных вискозиметров принципиально не различаются между собой. Кратко рассмотрим методику работы и расчета СМС для ротационных вискозиметров системы М. П. Воларовича [106]. Эксперимент сводится к подготовке прибора, измерениям и разборке. Подготовка включает установку вискозиметра по уровню, определение трения в подшипниках, заполнение рабочего объема исследуемым материалом и его термостатирование. Во время измерений меняют вес грузов Р или их массу т и для каждого груза определяют частоту вращения ротора N. Минимальное время замера составляет 2—2,5 с, каждый опыт состоит из 30—40 замеров, про­ водимых несколько раз с постепенным увеличением и уменьшением гру­ зов. Одновременно строят реограммы N(m). Затем прибор разбирается для мойки и определения высоты продукта по оставленным им следам на роторе. После этого определяют константы опыта (рис. 7,а). У вискозимет­ ра «Реотест» высота продукта всегда постоянна, во время опыта меняют не крутящий момент, а частоту вращения ротора.

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

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

2* 35

сдвига, пользуясь реограммами. Обсчет проводят по равновес­

ной кривой,

проходящей около всех точек [28, 94].

 

 

 

 

 

 

 

 

 

 

.*

 

 

Вращение ротора вызы-

 

 

 

 

 

 

 

 

 

Но

вает появление внутренних

80 L._____

 

 

 

 

 

 

напряжений

в продукте,

 

 

 

 

 

// -3000

 

 

 

 

 

который

находится между

50 V

 

 

 

 

 

 

 

 

ротором

и стаканом.

Эти

40

 

v

И Л?

 

 

 

 

касательные

напряжения

 

 

1 / /

 

 

 

2000

пропорциональны

сдвига­

30

 

 

 

 

 

 

ющим усилиям,

поэтому

 

 

 

 

 

 

 

 

 

графическую

и математи­

го

 

 

 

 

 

 

 

1

то

 

ом

 

 

 

 

 

ческую обработку

опытов

0,04

 

0,08

Ь,м

 

 

можно проводить

в

кон­

к

 

 

 

 

 

 

 

т ГР,

систентных

переменных

8 \

 

 

 

 

 

 

 

кг

s (0) по уравнениям (I—29),

 

 

 

 

 

 

 

0,05

\

\ >.. О*

 

 

 

 

 

 

 

(1—30),

(1—39) и

(1—40),

 

 

 

 

 

0Д4

но

можно

пользоваться и

 

 

 

 

 

 

 

 

003

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

М(т),

 

 

 

 

 

 

 

 

N(M), N(P).

 

 

 

/

Л

X

 

 

 

 

 

0,02

Модули

упругости

при

 

 

 

 

 

 

 

 

 

 

сдвиге определяют

по

за­

30

 

50

 

-70

 

2d

 

кону Гука

и графическим

 

 

 

 

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

относитель-

34

44

 

54

 

 

 

в-Ю3Па н°й

Деформации от време-

В-103Па

 

 

 

 

 

 

/

 

 

/

 

 

 

 

0

ни

действия

постоянного

 

 

 

 

 

?•

 

 

напряжения.

 

 

 

 

 

 

 

>

 

~ -

 

 

 

 

Наибольшую эффектив­

 

Г

 

 

 

4

 

 

 

ную вязкость, соответст­

 

/•

А > •

т

 

 

 

 

 

у

 

 

 

 

 

 

вующую

 

началу

пласти­

 

 

 

 

 

 

 

 

 

 

чно-вязкого течения, на­

 

10

20

 

30

40

 

ходят по уравнению Нью­

 

 

 

тона

 

 

 

 

 

 

Деления шкалы индикатора

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 7. Графики для определения кон­

71п= Д е /Д-с

 

 

 

стант приборов:

 

 

 

 

 

 

 

 

Дт

 

 

 

 

а — для

РВ-8;

1 — константа

конуса;

2 — со­

 

 

 

(1-45)

б — для

КП-3:

 

Д е /0

 

 

противление перемещению

штока

индикатора;

 

 

 

 

 

s — для мембраны

сдвигомера

Снмоняна.

 

 

 

 

 

 

 

где Дт — интервал времени, с; Де/Дт — скорость деформации для прямолинейного участка кривой млн аппроксимированного к пря­

мой криволинейного участка с небольшой выпуклостью.

Расчет эффективной вязкости для любого напряжения можно проводить по формуле (I—45), но удобнее при большом коли­

36

честве экспериментов пользоваться зависимостью, предложен­ ной М. П. Воларовичем [28, 29]:

 

т1э1)=

К -JJ- >

 

( 1 - 4 6 )

где К — постоянная прибора, м_1-с~2 (см.

рис. 7,

а);

N — частота вращения

ротора

вискозиметра,

с-1;

т —■масса вращающих

ротор

грузов,

кг.

 

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

описывается

моделью Шведова — Бингама. Ниже

приведены

основные расчетные зависимости для вычисления реологических

характеристик по вискозиметрическим данным,

полученным

на вискозиметре РВ-8.

 

 

 

 

 

 

 

При вычислении

реологических свойств возможны два слу­

чая [28]:

 

 

 

 

 

 

 

 

когда

1.

Сдвиг распространяется на всю толщину кольца,

масса сдвигающей нагрузки т0, вызывающей вращение ротора,

больше массы т' =

К\

0О * ПРИ которой сдвиг распространяется

на всю толщину кольца. Тогда пластическую вязкость у (Па-с)

вычисляют по формуле

 

 

 

 

 

 

 

 

 

 

Ti =~^~ (Km— К20о),

 

(I

47)

где 0и=

К0гпо— предельное

напряжение сдвига,

Па;

 

 

К, Ко,

Ки Ко— постоянные прибора, имеющие размерности соответствен­

 

 

но м_1'С_а; м-1-с-2; м-с2; 1;

кг.

 

 

 

т — масса

грузов, вращающих ротор,

 

 

II.

Сдвиг не распространяется на всю толщину кольца, ког­

да т0

т '.

Тогда

 

 

 

 

 

 

 

 

 

 

 

вп

1

(—

— 1 — In —

 

(1-48)

 

 

I = -^7 • —

 

или

 

 

N

4т:

\т й

т0

 

 

 

 

 

 

0п

 

 

 

 

 

 

 

 

 

7] =

F

 

(1—49)

 

 

 

 

 

 

 

 

 

 

N

 

U o

 

 

 

где величина

1 / т

,

т ,

 

обозначена

функциональной зави-

----

—1—In——

 

 

 

4л \ "'о

 

то

 

 

 

 

 

симостыо F

т

 

 

приведена в табл.

1.

 

 

; величина ее

 

 

37

 

 

 

 

 

Т а б л и ц а 1

т

 

т

 

т

F ( - )

гп0

\'Н0 )

то

\ М*о }

w0

1,1

0,00037

2,1

0,02849

3,1

0,07706

1,2

0,00140

2,2

0,03276

3,2

0,08249

1,3

0,00299

2,3

0,03716

3,3

0,08801

1,4

0,00505

2,4

0,04174

3,4

0,09358

1,5

0,00752

2,5

0,04644

3,5

0,09923

1.6

0,01034

2,6

0,05139

3,6

0,10494

1,7

0,01348

2,7

0,05623

3,7

0,11072

1,8

0,01688

2,8

0,06131

3,8

0,11656

1,9

0,02054

2,9

0,06645

3,9

0,12245

2,0

0,02441

3,0

0,07173

4,0

0,12839

При вычислениях по формулам (I—47), (I—49) значения величин удобнее брать непосредственно по кривой N(m), ко­ торую в процессе расчета корректируют. В окончательном виде кривая должна плавно проходить около всех опытных точек. При этом количество точек с каждой стороны кривой должно быть примерно одинаковым или, более точно, сумма положи­ тельных и отрицательных отклонений должна быть равна нулю (по среднеарифметическому методу).

Постоянные в приведенных уравнениях определяют по за­ висимостям, в которые входит поправка на сферический торец ротора

К =

 

Rg

 

 

RhК

RhRj

 

 

8л2

 

 

 

R l - R l

h +

 

 

 

R l - R t

 

/с0 =

-

Rg

 

 

2

-r:

 

 

 

 

 

 

 

 

2izRBh -J-

 

Ki =

2kR~

 

 

1

/?„

— — A;

K, = 2.303 — lg - r 2-

 

Rg

 

 

 

Ra

где R — радиус шкива (R =

0,02235

м)*;

 

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

м);

Ra— внутренний

радиус

стакана

(/?„= 0,019

(1-50)

(I—50а)

(1—506)

* В скобках даны геометрические размеры РВ-8; для приборов различ­ ных выпусков возможны некоторые колебания размеров.

38

Rи— радиус

ротора (RB — 0,01605

м);

h — глубина

погружения ротора

в продукт, м.

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

К =

 

Rg

 

(1-51)

RlRl

 

 

‘пр

 

 

R l - R l

 

 

 

 

где приведенная длина /tnp=

/1 +

Лэ,ш

и

 

Rв + Rh

Rb + Rh

(1-52)

^экв — R

 

+1

 

— + —

 

 

Rn

Ru

 

 

 

При соотношении радиусов 0,7 и более общая ошибка при определении константы по приближенному значению эквива­ лентной длины не превышает 0,5%.

Для удобства расчета на рис. 7 приведены зависимости констант от глубины погружения ротора в продукт.

Предельное напряжение сдвига 0о(Па) определяют по выра­ жению

 

Од -- КдШд ,

 

(1-53)

где гпд— масса нагрузки, при которой ротор

начинает вращаться, кг;

для

расчета эту величину определяют

графически

как отрезок

от

нуля до точки пересечения кривой N(m) с

осью абс­

цисс.

При обработке экспериментальных данных, полученных на вискози­ метре «Реотест», отличия от описанной выше методики будут следующие. Поправка на сферический торец не вводится. Вращающий момент М = = k*a( где k* — константа торзиона определяется тарировкой; а — угол закручивания по шкале прибора). ПНС определяется по формуле (I—34), градиент скорости — по уравнению (I—39), эффективная, или ньютонов­ ская, вязкость — по выражению (I—37). Высота продукта в приборе по­ стоянна и определяется высотой ротора.

Для измерения предельного напряжения сдвига широкое распространение в реологических исследованиях получили вследствие простоты устройства и надежности в работе кони­ ческие пластометры (см. рис. 5, г), сконструированные по мето­ ду акад. П. А. Ребиндера [102].

39

Соседние файлы в папке книги из ГПНТБ