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

книги из ГПНТБ / Пак, В. В. Шахтные вентиляционные установки местного проветривания

.pdf
Скачиваний:
10
Добавлен:
21.10.2023
Размер:
10.03 Mб
Скачать

 

 

Т а б л и ц а 2

 

Условия эксплуатации

Трубы

хорошие

плохие (выработ­

 

(выработки

ки мокрые, вода

 

сухие)

агрессивная)

Металлические.......................................................

30-48

22-26

Капроновые и полнхлорвшшловые................. ...

14-16

10—12

Прорезиненные.......................................................

8-10

4 -6

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

Таблица 3

Сечение выработки, м2 ........................

о

8

10

12

14

16

Максимальный диаметр трубопровода,

0,5

0,6

0,8

0,8

0,9

0,9

м ........................................................

Коэффициент сопротивления

трубопровода

удобно

определять

с помощью формулы

 

 

 

 

 

 

а = (3,7-f-lg<fT)2’

где

0,0038 — для жестких труб с прямолинейной навеской;

0,0050 — то же, с волнистой навеской; а = 0,0030 — гибкие трубы с прямолинейной навеской;

0,0045 — то же, с волнистой навеской; 0,0120— то же, со складками.

§ 2. Влияние величины и характера утечек на поле скоростей в трубопроводе

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

10

динамического сопротивления а, так как при этом характер поля скоростей, в принципе, может меняться.

Для определения влияния величины и характера утечек на поле скоростей в трубопроводе рассмотрим задачу [2] о безвихревом течении идеальной жидкости внутри имитирующей трубопровод полуполосы, расположенной на плоскости комплексного перемен­ ного z = х -\-iij (рис. 2). Для этого вдоль действительной оси пло­ скости вспомогательного переменного t =• £ 4- й] на отрезке —1

1 разместим распределенный источник интенсивностью q ( £),

 

 

 

аа - а

- /

1

Ь.

со

 

 

 

С

A

D

F

В

В I

 

Рис. 2. К выводу формулы (9)

 

 

 

а на отрезках —а ^

^

—1 и 1 ^

^ Ъ— распределенные стоки

интенсивностью q1 (h,)-

В произвольных

точках

£ =

ck (к = 1,

2, . . ., п) разместим сосредоточенные стоки обильности Qk, имити­ рующие местные разрывы трубопровода.

Тогда комплексный потенциал течения на плоскости t, порожда­

емого этой системой особенностей,

будет иметь вид [3]:

 

 

-1

 

 

j d l —j

q1(l)ln(t — l)d l —

 

—a

 

 

b

n

 

 

— J ? i(£ )h (f —

- 2 Qkin(t — ck) .

(6)

1

fc=1

 

 

Для того чтобы найти выражение

комплексного

потенциала

для плоскости z, осуществим с помощью интеграла Кристоффеля — Шварца конформное отображение верхней полуплоскости t на вну­ тренность полуполосы CADFBC, которая представляет собой тре­

угольник с углами ас =

0, aD —aF =

. Вершинам треугольника

CDF

в

плоскости t

соответствуют

точки: tD =

= — 1;

tp =

=

1; tc — lc =

± oo . Тогда имеем [4]

 

 

 

z —M^-rM=r-Ar N — Mexcb.t-\-N.

(7)

 

 

 

 

 

il

Из условия соответствия точек D и F на плоскостях z и t находим

значения постоянных М и N ^М =

7V = 0^, подставляя которые

в уравнение (7), получим выражение для функции, осуществляющей конформное отображение верхней полуплоскости t на внутренность полуполосы плоскости z

z = -~arch£.

(8)

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

.

.

du>

dw dt

,

V(z)=Vx - l V y = —

 

откуда, дифференцируя по t выражения (6) и (8), получим

 

 

 

Г

1

 

 

- 1

 

 

 

 

 

v (z) = \

sh nz

Г

g(S)d$

Г

gi(SMS

 

 

 

2

 

J ch nz—|

J

ch nz —|

 

 

 

 

.-1

 

 

-o

 

 

 

 

 

f

91 (E)<%

V

 

Qk

 

(9)

 

 

J

ch nz—|

£

ch nz —ch nz£

 

 

 

 

 

 

ft-1

 

 

 

 

 

где

zk — координаты точек разрыва трубопровода

на плоскости z.

 

С помощью соотношения (8) имеем

 

 

 

 

 

 

 

 

a = h = ch яLt

 

 

 

 

где

L — длина трубопровода

в

калибрах, так как диаметр трубы

для простоты принят равным единице.

q (|) и

(:£) и размещая

 

Задавая различные формы функций

сосредоточенные утечки Qk в определенных

точках

трубопровода,

из уравнения (9) определяем величину модуля скорости потока

 

 

I v (х, у)|= У v% (х, у) + vI (х, у),

 

 

с помощью

которого определялся коэффициент поля

скоростей

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

J |о (х, у) |2 dy

 

 

(10)

 

 

 

ао(*) :

 

J И®.У) I dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

и изучалась его зависимость от величины и характера утечек.

 

При этом были рассмотрены следующие варианты функции q (£):

 

a)

g(l) = g =

const;

 

б)

q (I) =

(1 — £) i;max;

12

Рис. 3. Влияние утечек (притечек) воздуха на характер его движения в трубопроводе при работе на:
а — всасывание; б — нагнетание

в сочетании со следующими видами функции qx (£):

a)

?i (Е) = <7i =

const;

 

б) q1(£) = qmax- k 1Z;

' в)

qx(£) = дтах^

х;

г)

?i(g) = ?n,aX — h i 2,

где постоянные q, vmax,

qv

qmax,

kv k2 и p выбирались таким обра­

зом, чтобы

суммарная

величина

утечек Q Q0 соответствовала

заданной величине коэффициента p Q = Q!Q0 = 3 (при L = 2000 м). Здесь Qо — расход воздуха на выходе из трубопровода.

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

(pQ -г; 3 при L — 2000 м)

является установившимся. Если в рассмотренной схеме заменить распреде­ ленный источник q (£) распределенным стоком той же интенсивности,

араспределенные стоки

qx ( |) — такими же рас­ пределенными источни­ ками, то нагнетательная схема проветривания за­ менится работой на вса­

сывание. При этом, как видно из выражения (9), изменится лишь направление скорости потока, все же остальное останется прежним. Таким образом, в рамках рассмотренной модели схемы всасывания - и нагнетания оказываются равноценными.

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

.с отрывом струй.

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

Из этого рисунка также видно, что условия подтекания потока к неплотности в трубопроводе при всасывании и нагнетании суще­ ственно различны. При всасывании поток подтекает к неплотности симметрично со всех сторон, а при нагнетании — сбоку, поворачи­

13

ваясь на 90°. Кроме того, в первом случае притечки эжектируются основным потоком в трубопроводе, а во втором случае утечки про­ никают в неподвижную атмосферу. Все это, как показывает анализ, 'увеличивает сопротивление той же неплотности в 1,3—1,5 раза при переходе от всасывания к нагнетанию. Поэтому можно ожидать, что такой переход будет сопровождаться не только уменьшением аэродинамического сопротивления, но и существенным снижением утечек.

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

Рпс. 4. К определешпо влияния величины притечек воздуха на аэродинами­ ческое сопротивление всасывающего трубопровода

трим задачу о подтекании притекающей струи расхода q через щель в трубопроводе и слиянии ее с основным потоком, имеющим расход Q. Схема рассматриваемого течения, происходящего на плоскости ком­ плексного переменного z = х + iy, показана на рис. 4, а. Характер­ ной особенностью рассматриваемого течения является отрыв при­ текающей струи в угловой точке С и образование свободной поверх­ ности СВ, вдоль которой модуль скорости, согласно постулатам

теории струй идеальной жидкости,

остается постоянным и равным v0.

Введем функцию Жуковского

 

 

со = In 4?- = In — -)- ей,

(И)

v

v

 

на плоскости которой (рис. 4, б) рассматриваемое течение происходит внутри полуполосы. Действительно, вдоль АВ аргумент скорости остается постоянным и равным нулю, а модуль увеличивается от vA — = QIH до у0. Вдоль AF также й = 0, а модуль скорости умень­ шается от vA до 0. Таким образом, отрезки АВ и A F на плоскости со располагаются вдоль действительной оси, причем точка В попадает в начало координат, а точка F уходит в бесконечность.

Вдоль ВС, как было сказано выше, модуль остается постоянным

и равным у0, аргумент изменяется от 0 до у , вследствие чего отре­

зок ВС на плоскости со располагается вдоль мнимой оси. Вдоль

14

границ DF и DC аргумент скорости постоянен и равен ft =

а мо­

дуль изменяется от vD =

qlh до 0 и

v 0 соответственно, благодаря

чему на плоскости со эти отрезки располагаются на прямой CDF,

параллельной оси абсцисс и отстоящей

от нее

на величину,

.

Введем вспомогательную плоскость

t = |

+ ir| (рис. 4,

в) и ото­

бразим с помощью интеграла Кристоффеля — Шварца [4]

ее верх­

нюю часть (т] ^ 0) на полуполосу плоскости со таким образом, чтобы

границам полуполосы на

плоскости t

соответствовала ось

абсцисс

(q = 0). При этом точка В располагается в бесконечности, точка С

в нуле, точки A, F и D — соответственно в точках

= —a; %F =

= —1;

= — d.

 

Нетрудно видеть, что функция

 

 

“ - Т

<12>

осуществляет отображение верхней полуплоскости t на внутренность полуполосы плоскости со. Действительно, подставляя в формулу

(12) tB — оо; tF = —1 и tc = 0, получаем соБ = 0; a>F = сос =

. п

= i 2 ‘

Так как в точках А и D расположены источники с обильностями Q и q, а в точке В — сток с обильностью Q + q, то, продолжая течение из верхней полуплоскости t в нижнюю, напишем выражение для комплексного потенциала течения в этой плоскости

со —-2-In (г+ а) -(--1 ln(i-M ).

(13)

Заметим, что в формулу (13) не вошел потенциал стока, так как

он расположен в бесконечно удаленной точке.

 

Учитывая, что v = dwldz,

из формулы (И ), интегрируя,

имеем

2 = —

at

+

(44)

vq J

 

 

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

Подставляем в выражение (14) уравнения (12) и (13) и получаем функцию, отображающую верхнюю полуплоскость t на область искомого течения плоскости z,

 

z = _ 2 _ r 2 ( l + ^ i n i H l f - M l n "и A

 

 

 

 

nv0 L

\

Q J

 

И— l

и + Л

 

 

. 1

i Au- -1

1 " UQ

D In

и— Ш

 

 

<15>

+ T

l n'

A u -\-1

и-|-iD

]D

 

A

 

 

 

 

где для краткости обозначено

 

 

 

 

 

j, ,— I f Vt + i .

л _

1 f Va + 1 .

T) — \

f 1 + УГ^

(16)

 

 

V VT-i

 

V Va-1

D~V

1~vr

 

15

Постоянную интегрирования с найдем, сопоставляя координаты точки С на плоскостях t и z,

'= ^ [я ('1 + ?) - 2 А

,Mt*-А

~ т «rcte4"+

 

+ » ! ^ - l ( £

+ - L

) l n ! ± i ] . .

(17)

Сопоставляя далее координаты точки .F на плоскостях t и z, имеем

_ L _ o

Aarctg A+ -j arctg — -i

(18)

Q

л + ( д + ± ) ln -g ± I

 

 

h = S -D .

(19)

 

yo

 

На основании формул (11) и (12) имеем

/У"ti

Подставляя в последнее выражение t = —а, найдем соотношение между скоростями и_4 и у0

( 20)

В точке F на плоскости z модуль скорости обращается в нуль, и так как функция (15) в этой точке не теряет конформности, здесь должно быть

откуда, используя соотношение (13), получим связь между пара­ метрами а и d

d - l — -§-(« — !)•

(21)

Комплекс формул (15)—(21) решает поставленную задачу и поз­ воляет рассчитать искомое течение.

В дальнейшем наибольший интерес представляет выражение для

коэффициента сужения струи

е = Н*/Н (см. рис. 4, а), которое

на оснований формулы (20) имеет вид:

 

- <?+g

1

( 22)

vqH ~

А

 

где pq — коэффициент, учитывающий притечки воздуха.

Задавая различные значения коэффициента притечек, по форму­ лам (16), (18), (21) и (22) была рассчитана зависимость коэффициента

16

сужения струи от натурального логарифма коэффициента притечек (см. рис. 6, сплошная линия).

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

могут быть найдены с помощью теоремы Борда — Карно

[6].

Если

эту величину потерь

 

давления

отнести к динамиче­

 

скому давлению потока за вих­

 

ревой областью, то выражение

1,2

для коэффициента потерь будет

 

иметь вид:

 

 

(23)

о,a

о,ч

Рис. 5.

Схема рассасывания вихре­

Рис. 6. Зависимость коэффициента су­

 

вой области

жения струи е и коэффициента потерь £

 

 

от величины притечек воздуха:

 

 

--------- г = 1 (lnpq);-----------£=Ф(1прд)

На

рис. 6 пунктиром показана зависимость £ = ф (In р 0).

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

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

2 Заказ 902

1,7

такой переход вполне допустим. Так, например, в работе [7] решена осесимметричная задача об истечении струи из конической воронки и приведены результаты численных расчетов для угла раскрытия конуса 2р = 90°. Полученный при этом коэффициент сжатия струи е0 = 0,750 мало отличается от коэффициента сжатия в соответству­ ющей плоской задаче епл = 0,745 [8].

В работе Роуза и Абуль-Фету с помощью электрогидродинамической аналогии было изучено осесимметричное струйное течение из отверстия диаметром d в дне бесконечно длинного круглого сосуда диаметром D. Значения коэффициента сжатия струи е0 при различных отношениях (dlD)2 приведены в табл. 4. Здесь же даны

коэффициенты

епл эквивалентной

плоской

задачи, рассчитанные

Р. Мизесом для тех же отношений

(ЫВ)2, где b — ширина

отвер­

стия, В — ширина сосуда [8].

 

 

Та б л и ц а 4

 

 

 

 

Коэффициент сжатия

Отношение (d/J5)2 = (b/B)3

0

0,25

0,50

0,75

1

 

 

Плоский случай.................................................

случай

0,611

0,644

0,689

0,757

1

Осесимметричный

0,612

0,644

0,691

0,757

1

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

§ 3. Расчет шахтных неплотных трубопроводов

При абсолютно плотном (герметичном) воздуховоде, когда произ­ водительность вентилятора Q равна расходу воздуха (?„, подава­ емому в забой, развиваемое вентилятором давление подсчитывается просто:

H = R Q l

Однако практически все шахтные воздуховоды негерметичны. Наличие утечек (притечек) существенно усложняет картину про­

ветривания и приводит (для обеспечения

заданного расхода

Q0)

к необходимости увеличения расхода (Q

Q0) и давления (h

Н)

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

(24)

и давления

(25)

18;

Как установлено рядом исследований

[9—11], утечки

воздуха

из

жестких воздуховодов

пропорциональны квадратному

корню

из статического давления.

По мнению Г. М. Леви

[12], для эластич­

ных

воздуховодов закон

утечек имеет,

якобы,

другой характер,

а именно: с увеличением давления относительная величина утечек снижается в результате самоуплотнения стыков отдельных звеньев трубопровода. Однако такие результаты получены в лабораторных условиях на новых трубах с проклейкой мест их сшивки. Испытания же серийно выпускаемых эластичных труб в шахтных условиях [10— 14] показали, что все-таки имеет место квадратичный закон утечек, так как основная их часть определяется местами сшивки полотнищ, множеством трещин, пор и разрывов, образующихся при эксплу­ атации.

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

циальной форме принято [9] представлять как

 

d Q = k d TVhdl,

(26)

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

Если произвольно принять, что pq = ph, то из формул (24)—(26) можно получить известную формулу проф. В. Н. Воронина [9]

PQ = ( i + i l d TV R y . .

(27)

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

Кроме формулы (27), где коэффициент p Q представляет собой параболу третьей степени от длины трубопровода, существует целый ряд других существенно отличающихся друг от друга формул и реко­

мендаций [12, ,15—18]. По Г. М. Леви [12], например,

зависимость

P q = /

(l) представляет собой степенную функцию с

показателем

меньше единицы. По данным Н. А. Богомолова [15],

 

 

pQ = 1 + aZli8B.

(28)

По

рекомендациям А. И. Ксенофонтовой [17] величину утечек

следует принимать из расчета Ъ % на каждые 100 м трубопровода

' P Q = 1 + 10-Hb%. _ (29)

Сравнение указанных формул и рекомендаций между собой пока­ зывает, что для воздуховодов относительно небольшой длины (Z < < 300 м) разница между значениями коэффициентов p Qлежит в пре­ делах допустимой погрешности и не является принципиальной. Что

же касается трубопроводов большой протяженности (Z

1000 м),

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

2*

19

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