Дифференциальные уравнения неустановившегося движения воздуха по рудничным воздуховодам (стр. 2 из 2)

(18)

где 2а – линеаризованный коэффициент аэродинамического сопротивления [формула (16)], но его численное значение пока не известно.

Значение этого коэффициента определим из условия равенства площадей, заключенных между вертикальными линиями граничных скоростей, между осью абсцисс и параболой (площадь, ограниченная контуром aABba) и между осью абсцисс и прямой линией (площадь, ограниченная контуром aαβba). Площадь, ограниченная параболой (площадь, ограниченная контуром аАВba , рис.2), равна:

(19)

Площадь, ограниченная прямой линией (площадь, ограниченная контуром aαβba , рис.2), равна:

(20)

По условию площади, описываемые формулами (19) и (20), равны. Приравняем их. Разность квадратов величин и разность кубов величин разложим на множители. Затем одинаковые множители слева и справа от равенства сократим. Проделаем эти выкладки.

(21)

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

(22)

Формулы (16) и (22) являются эквивалентными. Из сравнения этих формул следует среднее интегральное значение скорости, определяющей линейный режим аэродинамического сопротивления, при изменении этой скорости от v1 до v2, то есть

(23)

Подставим значение 2а из формул (16) и (22) в систему уравнений (15). В результате получим линеаризированную систему дифференциальных уравнений, описывающую переходные процессы воздуха в воздуховодах и выраженную через давление и скорость движения воздуха, т. е.

(24)

Для анализа переходных процессов при движении воздуха в воздуховодах систему уравнений (24) следует выразить через расход воздуха в данном поперечном сечении S. Для этого каждое слагаемое этой системы уравнений, содержащее скорость движения воздуха, умножим и разделим на площадь поперечного воздуховода. Тогда произведение скорости движения воздуха на площадь поперечного сечения воздуховода будет представлять собой расход воздуха (м3/с) в данном поперечном сечении, т. е.

(25)

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

(26)

Эта система уравнений была получена И.А. Чарным в его книге [6], но только другим путем.

Первое уравнение системы (26) учитывает как аэродинамическое сопротивление движению воздуха, так и его инерционные свойства. Каждый член этого уравнения представляет собой силу, приходящуюся на один кубометр движущегося воздуха (Н/м3).

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

Каждый коэффициент системы (26), стоящий перед зависимой переменной, имеет свой физический смысл и свое название в технической литературе, а именно:

2аρ/S – линеаризованное аэродинамическое сопротивление, приходящееся на единицу длины воздуховода, численно равное давлению, необходимому для создания единицы скорости одному кубометру воздуха в стационарном режиме,

;

ρ/S – коэффициент, учитывающий инерционность воздуха, численно равный давлению, необходимому для создания единицы ускорения одному кубометру воздуха,

, иногда этот коэффициент называют акустической массой или инерционностью [7];

S/(ρc2) – коэффициент, учитывающий сжимаемость воздуха, численно равный количеству воздуха, которое необходимо сжать для создания единицы давления на одном метре длины воздуховода,

, иногда этот коэффициент называют акустической гибкостью или податливостью [7].

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

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

Пример. Рассчитаем переходный процесс движения воздуха в коротком трубопроводе. Для этого используем только первое уравнение системы (26), так как второе уравнение для коротких воздуховодов обращается в нуль. При этом расход воздуха не зависит от длины потому, что из-за несжимаемости воздуха в коротких воздуховодах он по всей длине трубопровода не меняется. Поэтому частные производные в первом уравнении системы (26) можно заменить на обычные. Расчет проведем в операторной форме в виде интегрального преобразования Лапласа – Карсона [8]. После разделения переменных это уравнение будет иметь следующий вид:

(27)

где

– оператор Лапласа – Карсона ( – мнимая единица), Q(q) – расход воздуха в операторной форме (он не зависит от расстояния по оси воздуховода), р(x,q) – давление в операторной форме (оно является также и функцией расстояния по оси воздуховода). Остальные условные обозначения были приведены ранее.

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

(28)

Постоянную интегрирования А определим из граничных условий на конце воздуховода: при x = L давление p(L,q) = 0. Тогда для этого случая уравнение (28) будет иметь вид:

откуда постоянная интегрирования

(29)

Подставив значение А в формулу (28), получим выражение распределения операторного давления воздуха в трубопроводе в функции расстояния от начала воздуховода, т. е.

(30)

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

Перепишем выражение (30) для начального сечения воздуховода (х = 0). При этом учтем, что на вход воздуховода скачком подается постоянное давление ро = constant. В операторной форме Лапласа – Карсона это давление имеет выражение: p(0,q) = po. После подстановки в выражение (30) указанных значений величин и после преобразований получим формулу операторного расхода воздуха в начальном сечении воздуховода, т. е.

(31)

Применим к операторному выражению (31) обратное преобразование Лапласа – Карсона [8, формула 21.3]. Тогда получим выражение расхода воздуха в функции времени (t) в начальном сечении короткого воздуховода, а значит и расхода воздуха на протяжении всего трубопровода, т. е.

(32)

где Rаэр = 2аρL/S – линеаризованное аэродинамическое сопротивление короткого воздуховода длиной L [(1/c)∙(кг/м3)∙(м/м2) = =(кг∙м/с2)∙(с/м5) = Н∙с/м5].

Формулу (32) можно записать в относительных единицах (в безразмерном виде), а именно:

(33)

График расхода воздуха в относительных единицах, формула (33), показан на рис.3.

Рис.3. График изменения во времени расхода воздуха по длине

короткого трубопровода. Q*(t) – расход воздуха в относительных

единицах, 2at – безразмерное время.

Расход воздуха, обозначенный графиком на рис. 3, является одним и тем же по всей длине короткого воздуховода. График распределения давления по длине короткого воздуховода в приведенном примере во времени не меняется. Давление по длине короткого воздуховода от его начала до конца распределяется по линейному закону. Наибольшее постоянное во времени значение давления находится в начале воздуховода; а давление, равное нулю, – на конце воздуховода.

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