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

Аннотация

Лапко, В.В. Чередникова, О.Ю. – Идентификация оптимальной по сложности структуры аэродинамической модели воздушного потока коротких горных выработок шахтной вентиляционной сети

Введение

Установлено, что аэродинамические процессы в каждом элементе шахтной вентиляционной сети – горной выработке – с достаточной для практики точностью отображаются моделью типа цепочечных акустических –C фильтров с сосредоточенными параметрами. Длина отдельных участков выработки с распределенными параметрами, которые замещаются одной ячейкой с сосредоточенными параметрами, зависит в основном от граничной частоты существенных по амплитуде аэродинамических возмущений. Экспериментальные исследования показывают, что для типовых шахтных условий граничная частота существенных частот не превышает 0.5 рад/сек. Установлено, что для этих условий правомерной является цепочечная схема типа –C фильтров, в которых одним –C фильтром замещается участок выработки протяженностью до 340 метров. Однако, до настоящего времени не исследованы вопросы упрощения структуры модели воздушного потока коротких выработок протяженностью менее 340 метров, количество которых на практике составляет не меньше 30% от общего количества шахтных выработок.

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

Исследование переходных аэродинамических процессов в коротких выработках

В общем случае, модель коротких горных выработок длиной до 340 метров (рис.1) для реального спектра производственных возмущений (от 0 до 0.5 рад/сек) может быть представлена одним акустическим –C фильтром (рис.2).

На этих рисунках обозначены:

x – расстояние от начала выработки, м; l – длина выработки, м;

Rн – сопротивление нагрузки рассматриваемой горной выработки, нс/м⁵;

H – напор вентилятора, н/м²; L – ρl/F – акустическая масса воздушного потока выработки, нс² нс/м⁵; ρ – плотность воздуха, кг/м³;

F – сечение выработки, м²; C=F/ρ/a² – акустическая гибкость воздушного потока, м⁵/н; а – скорость звука в выработке, м/с;

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

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

Для определения погрешности упрощенной модели (рис.3) будем сравнивать ее динамические свойства и характеристики оригинала (рис.2) в эквивалентных режимах проветривания. В качестве исходного режима проветривания используем равновесное состояние оригинала (рис.2) и упрощенной модели (рис.3) при определенном установившемся значении напора вентилятора Нв=Н0 и расхода воздуха во входном сечении выработки Qн(0). Из схем замещения (рис.2,3) следует, что в этом равновесном состоянии расходы воздуха в обеих моделях будут одинаковыми и равными: Q_0н=Н₀/ R_н ; Q_0L=Q_0н ; Q_0k=Q_0н.

Динамическую погрешность в динамике упрощенной модели оценим по собственному «движению» оригинала (рис.2) и упрощенной модели (рис.3) после выключения вентилятора и падении напора вентилятора до Нв=0.

Тогда в упрощенной модели (рис.3) собственное «движение» системы будет описываться дифференциальным уравнением

В соответствии с предыдущим в исходном состоянии Qн(t=0)=Q0н. С учетом последнего получим:

где T_L=L/Rн – постоянная времени переходного процесса.

Если сопротивление нагрузки выработки будет пропорционально ее волновому сопротивлению, т.е. Rн=кZв, где Zв – волновое сопротивление короткой выработки, то постоянная времени переходного процесса будет равна

Подставив (3) в (2) окончательно получим

Для расчета собственного «движения» оригинала (рис.2) при выключении вентилятора преобразуем исходную схему замещения оригинала в схему, показанную на рис.4. На схеме обозначены (эквивалентное операторное сопротивление в конце короткой выработки):

Zс = 1/рC – операторное сопротивление акустической гибкости воздушного потока выработки.

При Rн=кZв получим

При выключении напора вентилятора Н собственное движение полной системы (рис.4) будет, очевидно, описываться характеристическим уравнением:

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

Подставив параметры выработки в выражения для акустической массы L и гибкости C, получим окончательное выражение для характеристического уравнения оригинала:

Решая характеристическое уравнение (6) относительно коэффициента нагрузки k и длины выработки λ=a/l найдем

Следовательно, если входное волновое сопротивление нагрузки более чем в два раза превышает волновое сопротивление оригинала (рис.2), переходный процесс в L–C-оригинале качественно отличается от движения в упрощенной L–модели независимо от протяженности (l) короткой выработки. Следовательно, упрощенная L–модель при больших значениях сопротивления нагрузки может использоваться для замещения коротких выработок только в особых случаях, когда в определенных рамках можно пренебречь колебательным характером переходных процессов. В противном случае, для замещения короткой выработки необходимо использовать акустический L–C фильтр (рис.2), который адекватно отображает как динамику воздушного потока выработки, так и наличие гармонических составляющих в переходных режимах.

Заключение

Рассмотрена возможность упрощения математической модели коротких (длиною до 340 метров) горных выработок шахтных вентиляционных сетей, количество которых составляет около 30% от всех выработок. В общем случае для реального спектра аэродинамических возмущений в шахтной сети локальная короткая выработка замещается акустическим L–C фильтром. Аналитическими исследованиями и методами математического моделирования показано, что при определенных условиях структура модели этих выработок может быть упрощена до акустического –фильтра.

Установлено, что критерием перехода к упрощенной L–модели является соотношение между выходным волновым сопротивлением модели короткой выработки и входным волновым сопротивлением выработки-нагрузки. Замещение локальной короткой выработки простейшим L–фильтром адекватно отображает ее динамические свойства в том случае, если ее волновое сопротивление в два или более раза больше входного волнового сопротивления выработки-нагрузки.

Список использованной литературы

1. Матикашвили Т.И. Составление математической модели вентиляционной системы с распределенными параметрами и автоматическими регуляторами. – Сб. «Вопросы динамики шахтных турбомашин и их сетей». Тбилиси, «Мецниереба», 1967.
2. Фельдман Л. П., Святный В.А. Исследование динамики воздушного потока на аналоговой вычислительной мшине. – Известия ВУЗов. Горный журнал, 1966, №6, с.63–67.
3. Петров Н.Н., Ермолаев П.Н. Методы синтеза систем автоматического регулирования вентиляторов главного проветривания.- В кн. «Автоматическое управление в горном деле». Сб. научных трудов, - Новосибирск ИГД СО АН СССР, 1971, с.23-49.