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

2142

.pdf
Скачиваний:
2
Добавлен:
17.06.2024
Размер:
5.19 Mб
Скачать

ENGINEERING SYSTEMS

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

P4 x x4 a3 x3 a2 x2 a1x a0 .

Разложив по степеням (x ) , будем иметь:

P4 x p4 p4 x 12 p4 x 2 16 p4 x 3 x 4 .

Выберем из условия

p4 0 ,

получим

α

a3

.

 

 

 

 

 

 

 

 

 

 

Тогда:

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a3

 

 

 

 

 

 

 

 

 

 

 

P

z z4 pz2 qz r ; z x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

 

3

 

 

 

a a a3

 

a a

 

 

 

a a

2

 

3a4

p a

 

 

a2 ,

q a

 

, r a

 

 

 

 

8

2 3

3

1 3

2 3

3 .

 

 

2

 

3

1

 

2

8

0

4

 

 

 

 

16

 

 

256

 

Такое представление P4 z удобно при графическом определении действительных

корней (или установлении их отсутствия). При p 0 имеем:

 

 

 

 

 

 

q

2

 

q2

 

x

 

a

 

4

 

 

 

 

 

 

2

 

z4 pz2 qz r

= z4 p

z

 

 

 

 

r

 

 

 

=

 

3

p x

β

 

S ,

2 p

4 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

p a

 

 

3

a2 , q a

a a a3

, r

a

a a

a a2

 

3a4

 

 

 

 

 

 

2 3

3

 

1 3

 

2 3

 

3

 

;

 

 

 

8

 

 

256

 

 

 

2

 

3

1

 

2

 

8

 

 

 

0

 

 

4

 

 

 

16

 

 

 

 

 

 

 

 

 

 

a

 

 

q

 

 

 

 

 

q2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β 3

 

 

 

 

,

S r

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 p

4 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходя из (1), рассмотрим декомпозицию:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z4 pz2

qz r z2

bz z2

bz

r

 

.

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем полагать q 0 ; в противном случае легко определяются все четыре корня,

как и в случае r 0 . Сравнивая коэффициенты многочленов слева и справа в (2), получим систему уравнений для определения b и :

b q

 

 

 

 

 

, f

 

f

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

 

(3)

( f

r

 

, f

 

 

p q2

 

 

 

2

 

) .

 

 

 

 

 

 

 

 

 

r 2 2

 

1

 

 

 

 

2

 

 

 

 

 

 

Алгоритм вычисления

непосредственно вытекает из вида графиков функций

f1 и f2 γ (рис.1, 2). Функция

f1 имеет вертикальную асимптоту 0 и на-

клонную y γ. Функция

f2 γ

при r 0

имеет вертикальную асимптоту

r и

горизонтальную y p , а при r 0

есть только горизонтальная асимптота y p .

 

f1max f1

r 2

 

r , f1min

f1

r 2 r , r 0 ;

 

f2max f2

r p

q2

 

,

r 0 .

 

4r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Regional architecture and engineering 2013

 

 

№3

151

 

 

 

 

 

ИНЖЕНЕРНЫЕ СИСТЕМЫ

 

 

y

 

 

 

4

2

r

 

y

 

 

 

 

 

 

y z

pz

 

y=z4+pz2+r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = -qz

 

p

 

 

p

 

z

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

z

Рис. 1. Иллюстрация к графо-аналитическому

 

Рис. 2. Иллюстрация к графо-

определению корней многочлена

 

 

 

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

P4(z), (p, q < 0, r > 0)

 

 

 

 

многочлена Р4(z), (p, q, r > 0)

При графо-аналитическом определении корень уравнения f1 f2

(рис. 3, 4)

будем искать в одном из видов:

 

 

0 1 ;

 

 

 

r ,

если r 0

 

 

 

r , если r 0 и

 

p q2

 

0 1 ;

(4)

 

 

 

 

 

 

 

 

4r

 

 

 

 

r , если r 0 и

 

p q2

 

0 1 .

 

 

 

 

 

 

 

 

 

4r

 

 

При r 0

, p

q

2

непосредственной проверкой можно убедиться:

r .

4r

 

 

 

 

 

 

 

 

 

 

2 r

y

 

 

 

 

 

 

 

 

 

y=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = f1( )

 

 

 

 

 

 

 

 

 

 

 

 

 

y=f1( )

 

y = f2( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = f2( )

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

y=f2( )

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

y = f1( )

 

 

 

 

 

 

 

 

 

p q2

 

Рис.3. Зависимости

f1 ,

f2

 

 

 

Рис. 4. Зависимости f1 , f2

 

 

при r 0

 

 

 

 

 

при r 0

 

152 Региональная архитектура и строительство 2013

№3

ENGINEERING SYSTEMS

Введем в рассмотрение функции:

 

g1

1

 

 

f1 r

 

и g2

 

 

 

 

1

 

 

 

 

 

f2 r ,

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где знаки выбираются в соответствии с (4);

0,1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При r 0 получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g 1 2 ;

 

g p q

 

 

 

 

2

 

 

 

 

, p

 

p

,

 

 

q

q2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

r

 

 

 

r

 

 

 

2

2 r

 

 

 

 

 

 

r

 

 

 

r

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При r 0

и p

 

 

q

r

 

 

 

p

 

 

 

 

p

 

,

 

 

 

 

q

 

 

 

 

q2

,

 

 

 

q

 

 

0

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

2 1,

 

g p q

 

 

2

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

r

 

 

 

r

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При r 0

и p

 

 

q

r

 

 

p

 

 

 

p

,

 

 

 

 

q

 

 

 

 

q2

,

 

 

 

 

q

 

 

0

 

 

 

r

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

1 2 ,

 

 

 

 

 

g p

q

 

 

 

 

 

2

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

r

 

r

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

(6)

(7)

Соответствующие графики функций g1 и g2 представлены на рис. 5, 6, 7.

y

y=g1( )

 

 

 

 

 

2

p r

 

q r

 

 

4

 

y=g2( )

pr qr

 

 

 

 

 

 

1

 

 

pr

 

 

 

Рис. 5. Функции g1 и g2

Рис. 6. Функции

g1

и g2

при r 0

 

 

 

 

 

qr

 

 

при r 0

и

p

 

 

 

 

 

 

 

r

 

4

 

 

 

 

 

 

 

Regional architecture and engineering 2013

№3 153

ИНЖЕНЕРНЫЕ СИСТЕМЫ y

y=g1( )

pr q r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y=g2( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pr qr

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.7. Функции g1 и g2

при r 0 и

pr

qr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

В соответствии с предыдущим получим следующий алгоритм определения :

– построение графиков функций y g1 , y g2 ;

определение абсциссы 0 пересечения графиков этих функций;

определение 0 в соответствии с (4).

Для уточнения значения

 

 

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

аппроксимацией функции f

f2 f1 в окрестности точки 0 :

 

f

f 0

f 0 0

 

f 0

0

2

.

(8)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f f

 

f

p q2

 

2

 

 

 

r

=

 

(9)

2

 

 

 

2

 

 

 

 

1

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2 r

 

 

 

0 02 r

 

 

r 2

 

 

 

 

p q2

0

 

 

 

 

 

0

 

2q2

 

 

 

 

 

 

0

 

 

 

 

r 02

2

 

 

 

 

r 02 3

 

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 3 r 02 2 2r r 02

 

r

 

 

 

 

3

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

3

0

 

 

 

 

 

 

 

 

 

 

2

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

При уточнении

методом

частичной линейной

аппроксимации функцию

f1

 

r

оставим без изменения (в силу ее простоты), а

f2 заменим линейной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

p q2

2

2q2

0 02 r

 

 

 

 

r

.

 

 

 

0

 

 

 

 

 

 

 

 

r 02 2

r 02 3

 

0

 

 

 

 

 

 

 

 

 

 

 

154 Региональная архитектура и строительство 2013

 

 

№3

ENGINEERING SYSTEMS

С увеличением порядка системы интерпретация результатов значительно затрудняется. В связи с этим при идентификации систем большего порядка целесообразно использование подсистем порядка, не превышающего 4.

Вкачестве иллюстрации рассмотрим характеристический многочлен

x4 1, 48x3 6, 46x2 7,7x 0,087

для объекта с параметрами по таблице 2 в [2] для случая M 0,9 ,

h 12 км, или

z

4

7, 26z

2

 

2,56z 2,04 ,

z x

0,37,

 

x

a3

 

 

 

z

4

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значения g1 ,

g2

указаны в таблице.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g2

 

 

 

 

 

 

 

 

g1

 

 

 

 

 

 

 

 

 

 

0

 

 

-

 

 

 

 

3,8

 

 

 

 

 

 

0,2

 

5,2

 

 

 

 

3,88

 

 

 

 

 

 

0,3

 

3,63

 

 

 

 

4,01

 

 

 

 

 

 

0,27

 

3,97

 

 

 

 

3,95

 

 

 

Построив графики функций y 2,56z,

y z4

 

7, 26z2

2,04 , убеждаемся, что

все корни многочлена действительные (рис. 8).

 

 

 

 

 

 

 

 

z1 2, 41,

 

z2 0,781,

z3 0,381,

 

z4 2,81.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-2,41

 

0,381

 

 

 

 

 

 

 

 

 

 

 

 

 

-0,781

 

2,81

 

 

 

 

 

 

y=2,56z

 

y=z4-7,26z2+2,04

 

 

 

 

 

 

 

Рис. 8. Иллюстрация к декомпозиции динамической системы

Соответственно

x1 2,78, x2 1,151, x3 0,011, x4 2, 44.

Декомпозиция полинома представится в виде

z4 7, 26z2 2,56z 2,04 z2 3, 26z 1,088 z2 3, 26z 1,875 .

Решая квадратичные уравнения с точностью до 10-2, получим те же корни zi и соответственно xi .

Regional architecture and engineering 2013

№3 155

ИНЖЕНЕРНЫЕ СИСТЕМЫ

Рассмотренный подход успешно использовался при идентификации параметров кинетических процессов формирования физико-механических характеристик полидисперсных материалов [3].

Список литературы

1.Данилов, А.М. Математическое моделирование управляющих воздействий оператора в эргатической системе / А.М. Данилов, И.А. Гарькина, Э.Р. Домке //

Вестник МАДИ. – 2011. – №2. – С.18–23.

2.Красовский, А.А. Системы автоматического управления летательных аппаратов

/А.А. Красовский, Ю.А. Вавилов, А.И. Сучков. – М.: ВВИА им. Н.Е. Жуковского, 1986. – 479 с.

3.Системный анализ в строительном материаловедении: моногр. / Ю.М. Баженов, И.А. Гарькина, А.М. Данилов, Е.В. Королев. – М.: МГСУ: Библиотека научных разработок и проектов, 2012. – 432 с.

References

1.Danilov, A.M. Mathematical modeling of the control actions of the operator in ergatic system / A.M. Danilov, I.A. Garkina, E.R. Domke // Herald MADI. – 2011. – № 2. – P.18–23.

2.Krasovskii, A.A. The automatic control systems of the aircraft / A.A. Krasovskii, Y.A. Vavilov, A.I. Suchkov. – M.: The air force engineering academy named N.E.Zhukovsky, 1986. – 479 p.

3.Systems analysis in building materials: monograph / Yu.M. Bazhenov. I.A. Garkina, A.M. Danilov, E.V. Korolev. – M.: Moscow State University of Civil Engineering: Library of scientific developments and projects, 2012. – 432 p.

156 Региональная архитектура и строительство 2013

№3

ENGINEERING SYSTEMS

УДК 621.928.9

Волгоградский государственный архитектурно-строительный университет

Россия, 400074, г. Волгоград, ул. Академическая, д. 1,

òåë.: (8442) 96-99-78; ôàêñ: (8442) 96-99-78

Боровков Дмитрий Павлович, кандидат технических наук, доцент

кафедры «Безопасность жизнедеятельности в техносфере»

Бурба Иван Владимирович, аспирант кафедры «Безопасность жизнедеятельности в техносфере»

Пензенский государственный университет архитектуры и строительства

Россия, 440028, г. Пенза, ул. Германа Титова, д.28,

òåë.: (8412) 48-27-37; ôàêñ: (8421) 48-74-77

Чичиров Константин Олегович, кандидат технических наук, доцент кафедры «Теплогазоснабжение и вентиляция»

Volgograd state university of architecture and building

Russia, 400074, Volgograd, 1, Akademicheskaya St.,

tel.: (8442) 96-99-78; fax: (8442) 96-99-78

Borovkov Dmitry Pavlovich,

Candidate of Sciences, Associate Professor of the department «Safety in Technosphere»

Burba Ivan Vladimirovich,

Postgraduate of the department «Safety in Technosphere»

Penza State University of Architecture and Construction

Russia, 440028, Penza, 28, German Titov St., tel.: (8412) 48-27-37; fax: (8412) 48-74-77

Chichirov Konstantin Olegovich,

Candidate of Sciences, Associate Professor of the department «Heat and gas supply»

ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ПЫЛЕУЛОВИТЕЛЕЙ НА ВСТРЕЧНЫХ ЗАКРУЧЕННЫХ ПОТОКАХ

Д.П. Боровков, И.В. Бурба, К.О. Чичиров

Представлена кинематическая структура течения закрученного потока в сепарационной камере пылеуловителей ВЗП. Определены оптимальные значения характеристик первичного ввода пылеуловителя серии ВИП.

Ключевые слова: пылеуловитель, встречные закрученные потоки, эффективность, аэродинамическое сопротивление.

EFFICIENCY PRECIPITATORS AT COUNTER CORKSCREW FLOW

D.P. Borovkov, I.V.Burba, K.O. Chichirov

Kinematic structure of swirling stream flow in separation chamber of dust collectors is given. Optimal characteristic values of VIP series dust collector primary input are determined.

Keywords: dust collector, counter swirling flows, efficiency, pressure drop.

На предприятиях строительной индустрии в качестве пылеулавливающих средств сухой очистки получили широкое распространение вихревые пылеуловители на встречных закрученных потоках (ВЗП) [1]. Их внедрение обусловлено рядом преимуществ по сравнению с циклонными: более высокой степенью улавливания мелкодисперсной пыли; меньшей чувствительностью фракционной эффективности к колебаниям расхода газа и концентрации пыли в нем; меньшей степенью абразивного износа аппарата; большей удельной производительностью; возможностью более эффективной очистки горячих газов и регулирования процесса пылеулавливания изменением соотношений расходов газа через потоки; в некоторых случаях – меньшей затратностью. Вышеперечисленное делает задачу изучения и оптимизации пылеуловителей на встречных закрученных потоках актуальной.

В ходе анализа различных типов пылеуловителей ВЗП установлено, что за основу при их проектировании брались типовые пылеуловители циклонного типа [1, 2]. Однако ввиду наличия вторичного закрученного потока структура течения пылега-

Regional architecture and engineering 2013

№3 157

ИНЖЕНЕРНЫЕ СИСТЕМЫ

зового потока в сепарационной камере пылеуловителей ВЗП отличается от таковой в пылеуловителях циклонного типа.

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

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

Анализ экспериментальных результатов определения окружной и осевой составляющих скоростей газового потока в пристенной зоне сепарационной камеры позволяет сделать предположение об автомодельности кинематической структуры течения в пристенной зоне по отношению к начальному значению интенсивности закрутки потока Ф*0. Для получения данных о параметрах пристеночного течения в пылеуловителях на встречных закрученных потоках были проведены экспериментальные исследования. Все результаты, полученные при различных значениях Ф*0, ложатся на гладкую кривую (рис. 1). Аналогичная картина наблюдается с отношением начального значения окружной составляющей скорости к значениям, полученным в ходе измерений по высоте сепарационной камерыUτ/Uτ0.

1,00

 

 

 

 

 

 

 

0,98

 

 

 

 

 

 

 

0,96

 

 

 

 

 

 

 

0,94

 

 

 

 

 

 

 

0,92

 

 

 

 

 

 

 

0,90

 

 

 

1

 

 

 

0,88

 

 

 

 

 

 

0,86

 

 

 

 

 

 

 

0,84

 

 

 

 

 

 

 

0,82

 

 

 

 

 

 

 

0,80

 

 

 

 

 

 

 

0,78

 

 

 

 

 

 

 

0,76

 

 

 

 

 

 

 

0,74

 

2

 

 

 

 

 

0,72

 

 

 

 

 

 

 

0,70

 

 

 

 

 

 

 

0,68

 

 

 

 

 

 

h/D

0,5

1,0

1,5

2,0

2,5

3,0

3,5

4,0

Рис. 1. Зависимости, характеризующие изменение локального пристеночного параметра интенсивности закрутки tg(φ) / tg(φ)0 (1) и тангенциальной составляющей скорости воздушного потока в пристеночной зоне Uτ/Uτ0 (2) по высоте сепарационной камеры

Полученные кривые представляют собой монотонно убывающие зависимости вида

tg(φ) / tg(φ)0 = 0,902(h/D)-0,129;

(1)

Uτ/Uτ0 = 0,88(h/D)-0,183.

(2)

158 Региональная архитектура и строительство 2013

№3

ENGINEERING SYSTEMS

Основными силами, оказывающими значимое влияние на пылевую частицу, движущуюся в пристенной зоне сепарационной камеры пылеуловителя ВЗП, являются: сила сопротивления обтеканию газовым потоком; центробежная сила; кориолисова сила; сила тяжести.

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

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

m

d 2r

d 2

dr 2

ч dч2

,

(3)

dt

2

m

 

r

 

8

 

 

 

dt

dt

 

 

где r – радиальная координата; dч диаметр частицы; ρ – плотность воздуха; λч – коэффициент сопротивления обтеканию частицы газовым потоком; m – масса частицы; t – время.

Проекция сил на подвижную касательную ось τ складывается из проекции силы сопротивления обтеканию потоком Fμ и кориолисовой силы. Заменяя угловую скорость линейной на подвижной оси, получим:

m

d 2

 

d

2

ч dч2

m

2 dr d

.

(4)

dt

2

 

dt

U

8

r dt

dt

 

 

 

 

 

 

 

Проекция сил на вертикальную (продольную) ось складывается из силы тяжести Fт и вертикальнойпроекциисилысопротивленияобтеканиючастицыгазовымпотокомFμ:

m d 2

2y

F Fт dy

Ux 2

xч dч2

mg .

(5)

dt

 

dt

 

8

 

 

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

m

d 2r

d 2

dr 2

 

d 2

;

d

 

U

;

dx

Ux ,

dt

2

m

 

r

 

ч

ч

dt

 

dt

 

 

 

dt

dt

8

 

 

R

 

где R – радиус сепарационной камеры. Используя значения (1) и (2), имеем:

m

d 2r

 

0,88U

 

x 0,183 2

dr 2

 

d 2

 

dt

2

mr

 

 

0

 

 

 

 

ч

ч

;

 

 

 

 

R

 

 

dt

8

 

d 0,88U 0

x 0,183 ;

dx 0,794U

0

x 0,312 .

 

 

dt

 

R

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6)

(7)

Для оценки влияния интенсивности закрутки, создаваемой верхним вводом, на эффективность улавливания аппарата ВЗП проведены численные эксперименты. При

Regional architecture and engineering 2013

№3 159

ИНЖЕНЕРНЫЕ СИСТЕМЫ

численном интегрировании полученной системы приняты следующие начальные условия:

за начало отсчета принимается точка, расположенная на нижней границе тангенциального ввода, наиболее близкая к центру поперечного сечения сепарационной камеры х = 0, r = 0,6R (ширина тангенциального патрубка принята равной 0,4R);

тангенциальная составляющая скорости

частицы

d

U

 

, а начальное

 

0

 

 

 

 

 

dt

 

 

d

 

U 0

 

 

 

значение угловой скорости соответственно

 

;

 

 

 

 

dt

R

 

 

 

 

 

 

 

 

 

 

 

значение вертикальной (продольной) составляющей скорости частицы на входе

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

потока dxdt0 Ux0 ;

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

ного потока, т.е. drdt0 0 .

Условием, позволяющим сделать вывод о сепарации пылевидной частицы, является значение радиальной координаты, равное радиусу сепарационной камеры r = R. При этом значение угловой координаты φ не имеет практического значения. В качестве переменных приняты размер пылевой частицы dч и начальный параметр интенсивности закрутки потока Ф*0. Свойства пылевых частиц (плотность ρч и коэффициент аэродинамического сопротивления λч) приняты идентичными свойствам пыли кварцевого песка соответствующей фракции.

r/R

1,0

2

1

0,9

3

0,8

4

1

0,7

0,4

0,8

1,2

1,6

2,0

2,4

2,8

3,2

3,6

h/D4,0

Рис. 2. Результаты расчета зависимости радиальной координаты пылевой частицы

от вертикальной при движении в сепарационной камере пылеуловителя ВЗП: 1 – dч = 37, Ф*0 = 4; 2 – dч = 35, Ф*0 = 4,4; 3 – dч = 33, Ф*0 = 4,8; 4 – dч = 25, Ф*0 = 6

160 Региональная архитектура и строительство 2013

№3

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