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

ОЦЕНКИ ВЕРОЯТНОСТНЫХ УРАВНЕНИЙ ВЗАИМОСВЯЗИ ПРОЦЕССОВ СУХОГО ДВУХСТАДИАЛЬНОГО ИЗМЕЛЬЧЕНИЯ

Адамбаев Мурат Джамантаевич

канд. техн. наук, профессор, КазНТУ им. К.И. Сатпаева, г. Алматы

Е-mail: adambaev т@тail. ru Шакенов Калижан Бахытжанович магистрант, КазНТУ им. К.И. Сатпаева, г. Алматы

E-mail: kalizhan-90@mail. ru Аденова Дана Бакытбековна магистрант, КазНТУ им. К.И. Сатпаева, г. Алматы

E-mail: dana. adenova@mail.ru

ESTIMATES OF PROBABLITY EQUATIONS INTERRELATION OF PROCESSES OF DRY TWO-STAGE SHREDDING

Murat Adambaev

Candidate of Technical Sciences, Professor, KazNTU after K.I. Satpaev, Almaty

Kalizhan Shakenov Graduate student, KazNTU after K.I. Satpaev, Almaty

Dana Adenova

Graduate student, KazNTU after K.I. Satpaev, Almaty

АННОТАЦИЯ

Экспериментально получены случайные реализации восьми основных переменных процесса сухого двухстадиального измельчения. Получен адекватный дискретный массив этих переменных с использованием предложенного метода определения времени реализации (Т) и квантования (А). Получены оценки вероятностных уравнений взаимосвязи переменных объекта исследования корреляционным и регрессионным методами.

ABSTRACT

Experimentally obtained random implementation of the eight main variables two-stage process of dry shredding. Obtained an adequate discrete array of variables, using the proposed method for determining execution time (T) and quantization (A). Obtained estimates of probability equations interrelation variables object of research with use correlation and regression methods.

корреляционный анализ; регрессионный анализ.

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

Это стало возможным после внедрения системы автоматического контроля покамерной и суммарной загрузок исследуемой мельницы [1], позволяющей контролировать такие основные переменные объекта, как уровни загрузки камер крупного и мелкого помолов, а также суммарную нагрузку мельницы по нагрузке ковшевого элеватора. Для данных исследований дополнительно производился контроль производительности тарельчатого питателя по частоте его вращения. Методом ручного пробоотбора и использования гр-нулометрических характеристик измельченного продукта на выходе мельницы определялись процентные содержания некондиционного класса — "крупки" и готового класса.

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

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

мельницы (определяется по ситовым характеристикам);

Случайные реализации переменных записаны на шеститочечном потенциометре (см. рис. 1). Для статистической обработки полученных

реализаций необходимо обоснованно выбрать время продолжительности реализации Т и интервал дискретизации или квантования Л1

Рисунок. 1 Схема регистрации переменных объекта

Для выбора Т получены автокорреляционные функции Яхх (т) переменных

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

Кж(т)< 0,05-Яхх(0). (1)

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

для определения Т, предложенной в [5]:

Т=10 иах

окончательно имеем:

Рисунок 2 Автокореляционные функции переменных

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

нестационарность процесса (рис. 3).

Рисунок 3 Случайные реализации Qn (—), II (—). Остальные реализации см. табл. П. 5 [2]

Интервал дискретизации или квантования осциллограмм выбран на основе применения теоремы отсчетов (теорема Котельника [3]) для наиболее высокочастотного параметра. По теореме Котельника непрерывная функция х(1:) с ограниченным частотным спектром Д^^х—1&П11П может быть дискретизирован последовательностью Хь Х2, .. .Хп с интервалом равным

Количество интервалов п на осциллограмме определяется по формуле

и = ^ = 2/-Д/. (4)

На кривой измеряется наибольший и наименьший периоды колебании и по НИМ определяется разность частот Дсй=Ютах-Ютт ПО формуле

А со = 2л

Т -Т ■

= 2л ^1ах ** . (5)

Т ■ Т . . ■

V тт тах / тах тт

Принимая Аа>х а/ и подставляя данные Лео в формулу (4), определяют величину п по формуле

. тах тт / /л\\

п =-----------------------------------. (6)

Т ■ Т ■

тах тт

Для выбранной реализации 7:(1:) имеем (см. рис. 3):

Тогда по формуле (6) вычисляем

= 4- 3,14- 240- (60- 20) 60-20

Увеличивая на 15 % полученное п для учета неравномерности частотной характеристики кривой имеем

п= 100,5 • 1,15=116,2.

Окончательно п принимаем равной 120. Тогда интервал дискретизации получается равным

Этот интервал дискретизации принят для квантования случайных реализаций переменных объекта. В табл. П.5 [2] приведены значения случайных реализаций параметров объекта с интервалом квантования 2 минуты, которые использованы для вычисления уравнений связи. Для решения поставленной задачи с применением ЭВМ использовалась стандартные программы из математического обеспечения ЭВМ-СПО10З и SPSS-16.

Результаты расчета экспериментальных данных на ЭВМ сведены в табл. 1. Анализ полученных результатов парной корреляции свидетельствует об одной их особенности — сравнительно невысоких значениях коэффициентов корреляции между параметрами двухстадиального сухого измельчения с объединенной сепарацией.

Причинами, обусловившими низкую коррелированность парных связей, являются:

этой связи.

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

Попутно отметим, что звукометрический сигнал камеры крупного помола 7^) зависит от уровня загрузки этой камеры рудой, которая, в свою очередь, определяется суммой двух материальных потоков — потока исходной руды Qn(t) и части потока циркулирующей нагрузки, поступающей в эту камеру kS. Следовательно, данный звукометрический сигнал 7^) зависит, по крайней мере, от двух входных переменных. Поэтому необходимо определить корреляцию трех переменных 7^), Qn(t) S1(t), т. е. искать зависимость вида 71=1^п^). Аналогично величина циркулирующей нагрузки S, процентные содержания вновь образованного готового класса у и некондиционного класса Р в единице веса разгрузки мельницы зависят а основном от уровней загрузки

обеих камер. Необходимо определить регрессионные зависимости вида:

гп); у=ґ(гь гп); р=ад, гп).

На производительность агрегата по готовому продукту Qг влияют уровни загрузки камер материалом, а также величина циркулирующей нагрузки в системе [6]. Регрессионная зависимость принимает вид: Qг=f(ZI, 2П, S).

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

Теснота связи между переменными оценивается корреляционным отношением г). Сравнивая по абсолютной величине полученные

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

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

том, что моменты измерения или квантования входных переменных сдвигаются на отрезки времени Д^, А{2,....от момента регистрации выходной координаты.

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

входными и выходными переменными с достаточной точностью не представляется возможным. В работе оценка величин сдвига Д^, А{2, ... произведена по взаимным корреляционным функциям (рис. 4). Каждая из величин А{ равна моменту времени, в котором взаимная корреляционная функция ЯУХ(т) достигает максимального значения.

Результаты расчета на ЭВМ корреляционных зависимостей с учетом разобранных причин, влияющих на коррелированность параметров объекта,

приведены в табл. 2. Анализ полученных данных позволяет сделать заключение, что статические характеристики двухкамерной шаровой мельницы хорошо согласуются по своим формам с известными исследованиями в этой области [4]. Некоторые регрессионные зависимости даны на рис. 5.

Оценка надежности коэффициентов корреляции производилась по формуле

Ч Чл]л &“1^Коэффициент корреляции значим при (1 — табличное значение). В программе СПО10З вычисляются критерии Стьюдента и Фишера, которые также сравниваются с табличными значениями этих критериев для Р=0,05 (Р — уровень значимости).

Рисунок 4. Взаимнокорреляционные функции переменных

Рисунок 5. Регрессионные зависимости параметров объекта

Таблица 1

Парная корреляция переменных объекта_____________________

№ Зависимост у=а + вх 2 у=а + вх + сх

пп ь а в ч а в с Г)

Таблица 2

Регрессионные уравнения взаимосвязи переменных объекта

№ пп Уравнение регрессии Ч12 Ч13 Ч23 Я

где: а, в,с — коэффициенты уравнения регрессии;

ч — коэффициент парной корреляции;

Ну — коэффициенты частной корреляции; Я — общий коэффициент корреляции;

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