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

учеб реология Арет

.pdf
Скачиваний:
84
Добавлен:
15.03.2016
Размер:
2.21 Mб
Скачать

 

 

 

n x

n y

 

n y

 

vz2 (x, y)

 

sin

 

Dnch

 

Ensh

 

 

(34)

b

b

b

 

n 1,2,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем находить константы интегрирования по краевым условиям:

vz2 (x,0) 0; y 0;ch

n y

 

 

1;sh

n y

 

0

(35)

b

b

 

 

 

 

y 0

 

 

y 0

 

Запишем формулу (34) с учетом условия (35):

 

 

 

n x

 

 

 

 

 

 

 

 

 

0

sin

 

 

 

 

 

(36)

b

 

 

 

 

 

n 1,2,3

 

 

 

 

 

 

 

 

 

 

 

 

Используем далее теорию рядов Фурье :

 

 

2

b

 

 

n x

 

 

 

Dn

 

vz2

(x,0)sin

dx

 

(40)

b

 

 

 

 

0

 

 

b

 

 

Поскольку

vz2 (x,0) 0 , то Dn

0.

(41)

Краевое условие для скоростей на верхней стенке канала, предполагая как и ранее

условие прилипаемости среды к материалу насоса, имеет вид vz2 (x,h) Vz . (42)

Тогда выражение (34) можно записать так

Vz

Ensh n hsin n x

(43)

 

 

 

 

 

 

 

n 1,2,3

b

 

b

 

 

 

 

 

 

По теории рядов Фурье запишем

181

 

 

E

sh

n h

(44)

 

b

 

n

n

 

 

и

Ensh n h

2

Vz

sin n x dx

(45)

 

 

 

 

b

 

 

 

 

 

b

 

b

0

 

b

 

 

Проведя интегрирование в правой части уравнения (45), выразим

En

2Vz

1 cosn

(46)

n

 

sh

n h

 

 

 

 

b

Поскольку четные n=2,4,6… дают тривиальное решение En 0, следовательно будем учитывать только нечетные слагаемые n=1,3,5…, при которых числитель в правой части формулы (46) равен 2. Теперь распределение скоростей течения в ка-

нале определяется формулой вида:

 

 

 

n x

2V

1 cosn

 

n y

vz2 (x, y)

sin

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

sh

 

(47)

 

 

 

 

 

 

 

 

 

n h

 

 

 

n 1,3,5

 

 

b

n

 

 

sh

 

 

b

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4Vz

 

 

1

 

 

n x

 

sh

n y

 

 

 

 

 

 

 

 

 

b

 

 

 

vz2 (x, y)

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

.

(48)

 

n

 

b

 

 

 

 

 

n h

 

n 1,3,5

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

182

Двойным интегрирование получим расход среды, обусловленный движением верх-

ней стенки канала червяка:

 

 

 

 

 

4Vz

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

n x

 

 

sh

 

n y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

n h

 

 

 

 

 

 

 

 

0

 

0

 

 

n 1,3,5

n

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4Vz

 

 

 

 

 

 

 

 

1

 

 

 

 

 

n x

 

 

 

 

 

 

 

 

 

 

sh

n y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

dy dx

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n h

 

 

 

 

 

 

 

0

 

n 1,3,5

n

 

 

 

 

 

 

 

 

b

 

 

 

 

0

 

 

 

 

 

sh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4V

z

 

 

 

 

 

1

 

 

b

 

 

 

 

 

 

 

 

 

n x

 

b

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

n h

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1,3,5

n

 

 

n

 

 

 

 

 

 

 

 

 

b

 

 

0

 

 

 

n sh

n h

b

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4V

z

 

 

 

 

 

2b2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

n h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

n3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

n 1,3,5

2

 

 

sh

n h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

bh 16b

 

 

 

1

 

 

 

 

n h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(49)

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

th

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

3

 

 

 

 

n

3

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h n 1,3,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что сомножитель перед квадратными скобками в последнем выражении

(49) совпадает с первым слагаемым в формуле расхода среды в упрощенной ли-

нейной теории червячных нагнетателей :

Q

V

bh

 

 

bh3

P

 

z

 

 

 

 

 

 

(50)

2

 

 

 

 

12

z

 

Следовательно, выражение в квадратных скобках в последнем выражении формулы

(49) можно рассматривать как поправочный коэффициент, зависящий от отношения ширины канала b к ее глубине h и учитывающий тормозящее действие боковых не-

подвижных стенок глубоких червячного канала. По этому поправочному коэффи-

183

циенту можно расчетным путем оценить погрешность первого слагаемого в форму-

ле расхода упрощенной линейной теории червячных нагнетателей и определить применимость упрощенной теории в расчетах. Например, при h=b погрешность пер-

вого слагаемого в формуле (50) составляет около 50%.

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

ствием перепада давления была решена ранее, откуда запишем:

 

 

 

 

 

 

vz1

 

1

 

P y2

 

 

 

 

 

 

 

 

 

 

 

z 2

yh 4h2

2 3

 

n

 

 

 

ch

 

 

2x b

 

 

 

 

 

2h

 

 

 

 

 

 

 

 

 

 

 

 

n b

n 1,3,5

 

3

 

n ch

 

 

 

 

 

2h

 

 

 

 

 

n y

sin (51)

h

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

 

 

bh3

 

P

192h

 

1

n b

 

Q1

 

 

 

 

1

 

 

 

 

 

 

th

 

 

(52)

 

 

5

b

n

5

2h

 

12

z

 

 

n 1,3,5

 

 

 

 

Учитывая формулы (11), (49) и (52), получим окончательно для червячного нагнета-

теля с глубокими каналами выражение для расходно-напорной характеристики

 

 

V

bh 16b

 

 

 

1

 

n h

 

 

Q

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

th

 

 

 

 

 

 

 

2

 

 

 

3

 

 

 

 

n

3

 

 

 

b

 

 

 

 

 

 

 

h n 1,3,5

 

 

 

 

 

 

 

 

(53)

 

bh3

 

P

 

 

192h

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

n b

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

th

 

 

 

 

 

 

z

 

 

 

5

b

 

 

 

n

5

 

 

2h

12

 

 

 

 

 

 

 

 

 

n 1,3,5

 

 

 

 

 

Если представить выражения в квадратных скобках формулы (53) следующим об-

разом

184

 

 

16b

 

1

n h

 

kV

 

 

 

 

th

 

 

(54)

 

n3

b

 

 

3h n 1,3,5

 

 

 

 

 

192h

 

1

n b

 

 

kP

1

 

 

 

th

 

 

,

(55)

5b

n5

2h

 

 

n 1,3,5

 

 

 

 

то формулу (53) можно представить в виде:

 

V

bh

 

 

 

bh3

P

 

 

Q

z

 

kV

 

 

 

 

 

kP

,

(56)

2

 

 

 

 

 

12

z

 

 

а поправочные коэффициенты kV и kP , как функции отношения глубины канала к ее ширине, можно затабулировать.(табл.1)

Значения поправочных коэффициентов kV

и kP расходно-напорной

характе-

ристики червячного нагнетателя

 

 

 

 

 

 

 

 

 

 

 

 

 

Табл.1

 

 

 

 

 

 

 

 

 

 

 

h/b

0

0,4

0,8

1,2

 

1,6

2

 

 

 

 

 

 

 

 

 

 

 

 

kV

1

0,783

0,582

0,433

 

0,335

0,27

 

 

 

 

 

 

 

 

 

 

 

 

kP

1

0,748

0,515

0,346

 

0,239

0172

 

 

 

 

 

 

 

 

 

 

 

 

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

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

185

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

ным единице. В суммируемую функцию введен сомножитель, которого нет в формулах (54), (55) и который обеспечивает суммирование по нечетным индексам

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

дит в упрощенную, следовательно для таких червяков расчет можно проводить не по формуле (53), а по формуле (9).

Программа расчета поправочных коэффициентов для гидродинамической теории червячных нагнетателей.

N

 

50

 

 

 

k

 

1

 

cos(n )

 

n

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1(n)

1

 

1

 

tanh

f2(n)

1

 

1

 

 

cos(n

tanh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

3

 

2

2

 

 

 

 

 

n

5

 

 

 

 

2

 

 

 

2 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1(n)

0.969

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f2(n) 0.922

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

f1(n) 0.5

KV

 

 

1

 

 

 

 

192k

 

 

 

 

f2(n)

0.422

KP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

( )3 k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно также построить график вида

186

N

 

50

 

 

 

N

 

1

1

 

cos(n )

 

n

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tanh

 

 

 

 

 

 

 

 

 

 

 

 

 

n = 1 n

3

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )3 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1(x)

 

 

 

1

 

 

192x

 

 

 

 

 

 

 

1

 

1

 

cos(n

)

tanh

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

n = 1 n

5

 

 

 

2

 

 

 

 

 

2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

f(x)

 

 

f1(x)

0.5

 

 

1

2

 

 

x

Поскольку расходно-напорные характеристики различных каналов были рассмотре-

ны ранее, то можно перейти к расчета червячных экструдеров по методу совмещен-

ных напорно-расходных характеристик нагнетателя и формующей матрицы.

4_5_3

4.5.3 Расчет червячных экструдеров по методу совмещенных расходно-

напорных характеристик.

Пусть на экструдере формуется жгут круглого сечения. Тогда систему червячный нагнетатель-матрица можно рассматривать как систему генератор массы - потреби-

тель массы и найти рабочую объемную производительность и давление в предмат-

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

(50) и уравнения Пуазейля:

187

Q

V

bh

 

 

bh3

P

 

z

 

 

 

 

 

 

(57)

2

 

 

 

 

12

z

 

Q

R2 p

 

(58)

 

 

8L

Полагая падение давления в канале червяка линейным, запишем

 

P

 

Psin

 

 

 

 

 

,

(59)

 

 

 

z

 

Lч

 

где Lч –длина червяка;

- угол подъема винтовой линии червяка.

Совместное решение уравнений (57),(58) и (59) дает давление в предматричной ка-

мере:

P

12 VzbhLLч

(60)

2bh3Lsin 3 R2Lч

и объемную производительность экструдера

 

R2

 

3 VzbhLLч

 

 

Q

 

 

 

 

(61)

2L

2bh3Lsin 3 R2 L

 

 

 

 

 

 

 

ч

 

5_1

188

5.1.Проблема вероятности формосохранения пищевых изделий

Проблема вероятности формосохранения пищевых изделий, их товарного вида и

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

ческое состояние под действием собственной массы или массы расположенных вы-

ше изделий.

С помощью пластометра Ребиндера [1] можно найти предельное напряжение сдвига материала изделий по формуле :

 

0

k

P

,

(1)

h2

 

 

 

 

где 0 - предельное напряжение сдвига материала изделия, Н/м2; k - коэффициент прибора, функция угла при вершине конуса;

P -вертикальная сила , вдавливающая конус в материал, Н; h - глубина проникновения конуса в материал, м .

При нормальном законе распределения вероятностей напряжений нагружения и со-

противления можно рассчитать вероятность формосохранения [5]. Определяем квантиль нормального закона распределения по формуле суперпозиции

 

 

 

k

P

 

H

 

 

 

h2

UP

 

 

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

k

 

 

s

 

H n

h2

 

 

 

 

 

 

 

 

 

 

 

 

где UP

 

 

- квантиль нормального закона распределения;

 

 

 

 

-

 

удельный вес материала изделия, Н/м3;

H- высота штабелированных изделий, м;

S - вариация предельного напряжения сдвига материала изделия;

189

n - вариация предельного напряжения сдвига нагружения.

Теперь вероятность формосохранения изделия находим с использованием функ-

ции Лапласа по формуле:

 

1

 

UP

 

x2

(11)

 

 

 

e

 

 

 

P(UP) 0,5

 

2 dx

 

 

 

 

 

 

 

 

 

2

 

 

 

0

 

 

 

 

 

С помощью формул (2) –(11), и пакета программ Mathcad 7 Pofessional

были про-

ведены расчеты вероятности формосохранения при следующих данных

Z =500;

P=10 Н; h=1,2 м;

= 0,09 Н/м3; H=11 м; S

=0,03 ; n = 0,05.

 

190