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

Аннотация

Лапко, В.В. Чередникова, О.Ю. - Разработка методами частотного анализа аэродинамических моделей локальных ветвей шахтной вентиляционной сети. На основе методов частотного анализа для реального спектра производственных возмущений разработаны приближенные аэродинамические модели локальных ветвей шахтной сети. Показано, что длинные ветви шахтной сети с достаточной для практики точностью могут быть аппроксимированы цепочечной схемой из коротких ветвей с сосредоточенными параметрами. Установлена зависимость между числом коротких ветвей приближенной модели с сосредоточенными параметрами и протяженностью локальной ветви сети. Разработанные приближенные аэродинамические модели локальных ветвей могут быть использованы для синтеза математической модели шахтной сети как объекта управления.

Введение

В произвольных переходных режимах состояние горных выработок шахтной вентиляционной сети как аэродинамического объекта описывается телеграфными уравнениями волновых движений в частных производных. Однако, при относительно низком темпе переходных процессов эти уравнения могут быть сведены к более простой математической модели – к обыкновенным дифференциальным уравнениям. Для количественной характеристики темпа переходных процессов в выработках удобно использовать соотношение длины ветви l и минимальной длины бегущей волны сжатия воздуха вдоль выработки при действии производственных возмущающих воздействий λ_min=(λ_min= a.Tmin, где а ≅340м/с – скорость звука в воздухе, Tmin – период верхней гармоники существенного спектра производственных возмущающих воздействий в шахтной сети). Для реальных шахтных условий период верхней гармоники возмущающих воздействий составляет примерно 12,5 секунд. Поэтому минимальная длина бегущей волны сжатия воздуха вдоль выработки при реальных производственных возмущениях в шахтных выработках практически не превышает λ_min≅4250м. Для конкретной выработки сети этот темп возмущающих воздействий является либо низким, либо высоким, в зависимости от соотношения длины выработки (l) и минимальной длины бегущей волны сжатия воздуха в выработке (λ_min). В шахтных условиях большая часть выработок меньше 600м (90% всех выработок). В то же время, максимальная длина выработок может достигать 4000м. Следовательно, длина выработок реальной шахтной сети может быть как существенно меньше (в 7 и более раз), так и примерно равна длине бегущей волны сжатия воздуха вдоль выработки при реальных производственных возмущениях. Таким образом, в рабочей полосе частот возмущающих воздействий (примерно до 0.08 Гц) выработки шахтной сети по отношению к длине бегущей волны сжатия воздуха могут считаться как короткими (l<<λ_min), так и длинными (l≈λ_min). При этом аэродинамическая модель выработок должна воспроизводить тем более сложный характер бегущей волны сжатия воздуха при реальных производственных возмущающих воздействиях, чем больше геометрическая длина выработки по отношению к длине волны. В силу этого, для выработок шахтной сети различной протяженности должны использоваться различные (минимальные по сложности) математические модели.

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

Априорная разработка математической модели коротких выработок шахтной вентиляционной сети

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

Распространение аэродинамических воздействий вдоль любой выработки ограничено скоростью бегущей волны сжатия воздуха (скоростью звука) а≅340м/с. В силу этого время затрачиваемое на то, чтобы в выработке (ветви шахтной вентиляционной сети) произошло изменение потока воздуха на ее конце (выходе выработки), в общем случае должно быть не менее tд=l/a, где tД – запаздывание сигнала на границе (конце) ветви.

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

где λ_min=а Т_min=а (2π)/ω max; ωmax – круговая частота верхней гармоники спектра реальных гармонических сигналов.

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

где |S_q(jω)| – модуль спектральной плотности входного возмущающего воздействия в выработке по воздуху Q(t), возникающего при t = 0 и отсутствующего при t<0.

В соответствии с этим для наиболее типичных в шахтной вентиляционной сети апериодических аэродинамических возмущений по воздуху Q(t) полоса существенных частот составит

где α=1/τ_EXP; τ_EXP – постоянная времени экспоненциального возмущения по воздуху Q(t).

В силу этого, если в этом выражении не учитывать частоты, мощность которых составляет меньше 5%, для апериодических возмущений верхняя граница спектра существенных частот составит:

В общем случае при известной границе существенных частот ?max в соответствии с (1) выражение для критичной длины коротких выработок можно записать в виде

Как известно, для условий реальных шахт верхняя граница полосы существенных частот ω_МАХ≈0.5 рад/с. Следовательно, согласно в условиях реальных производственных возмущений в соответствии с (5) в качестве предельной длины ветви, в которой можно пренебречь распределенностью параметров, ориентировочно можно принимать выработку длиною

Отметим также, что в соответствии с (4), эти выработки (600м и меньше) могут рассматриваться как схемы с сосредоточенными параметрами, т.е короткими, также для непериодических экспоненциальных возмущающих воздействий с постоянной времени τ_EXP≥40c.

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

В качестве исследуемой модели короткой выработки выберем простой четырехполюсник, например, симметричный Т – четырехполюсник с сосредоточенными акустическими и аэродинамическими параметрами (рис.1):

Изменения расхода и давления воздуха относительно установившихся значений соответственно на входе (Qi, Pi) и выходе (Q0, P0) Т – модели короткой выработки могут быть описаны (в матричной форме) в виде:

Z1 = (рL+R)l – продольное операторное аэродинамическое сопротивление короткой выработки;

Y2 = (g + рC) – поперечная операторная проводимость короткой выработки;

Z2 = 1 / Y2 – поперечное операторное аэродинамическое сопротивление схемы замещения короткой выработки;

р – оператор Лапласа.

Кроме того, в силу обратимости построенной модели (рис. 1) для (7) и (8) справедливо соотношение

A*D–B*E = 1 (9)

Если ввести для рассматриваемого Т – четырехполюсника вектор входных (Xi) и выходных (Y0) переменных, уравнения (7), (8) и (9) окончательно могут быть записаны в форме:

Заключение

На основании методов частотного анализа в работе предложен оригинальный подход к разработке приближенных аэродинамических моделей локальных ветвей шахтной вентиляционной сети. Для оценки точности приближенных аэродинамических моделей локальных ветвей сети использовались частотные характеристики согласованной по входу и выходу выработки без потерь, вид которых хорошо известен и теоретически доказан для любого спектра аэродинамических возмущений. Приближенные модели выработки строились из конечной цепи коротких участков, каждый из которых аппроксимировал достаточно малый конечный элемент выработки. Адекватность приближенной цепочечной модели выработки из определенного количества коротких участков определялась сравнением реальной и теоретически известной частотных характеристик согласованной по входу и выходу выработки. Для реального спектра производственных возмущений шахтной сети показано, что 90% всех ветвей шахтной сети (длиною до 600м) с достаточной для практики точностью может быть аппроксимировано одним коротким участком с сосредоточенными акустическими и аэродинамическими параметрами. Установлена зависимость порядка приближенной цепочечной модели ветви от ее протяженности для нетипичных длинных ветвей шахтной сети. Показано, что длина коротких участков цепочечной модели длинных выработок не должна превышать 300м. Построенные модели выработок могут быть использованы при разработке аэродинамической модели шахтной сети как объекта управления.

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

1. П.Н. Ермолаев и др. Частотные свойства шахтной вентиляционной сети как объекта автоматического регулирования /– В кн. : Автоматическое управление в горном деле. – Новосибирск: Б. И., 1971, с.16 – 23.
2. Петров Н.Н., Ермолаев П.Н., Пономарев П.Т. Электронная модель системы автоматического управления проветривания шахт. В кн.: Автоматическое управление в горном деле. Новосибирск: Б.И.,1971, с.89 – 93. .
3. Угрюмов Е.П. Цифровая схемотехника. Санкт-Петербург, 2000. – 517с. .
4. Уэйкерли Дж. Проектирование цифровых устройств, том 2. – М.: Постмаркет, 2002.- 1064с. .
5. Berstein H. Handbuch der praktischen Elektronik. Band 2. – Fransis’Verlag – 2003. – 567c. .
6. Гоноровский И.С. Радиотехнические цепи и сигналы. – М.: Советское радио, 1963. – 684с. .
7. Воронов А.А. Основы теории автоматического управления. Часть 2. Специальные линейные и нелинейные системы автоматического регулирования одной величины. – М.: Энергия, 1966. – 325с. .
8. Атабеков Г.И. Теоретические основы электротехники, часть 1. – М.: Энергия, 1964. – 31c.