Научная статья на тему 'Идентификация параметров динамических объектов посредством приложения Simulink пакета Matlab'

Идентификация параметров динамических объектов посредством приложения Simulink пакета Matlab Текст научной статьи по специальности «Механика и машиностроение»

CC BY
368
72
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по механике и машиностроению, автор научной работы — Бровцин В.Н.

В статье приведен пример идентификации параметров пахотного аг-регата по каналу управления положением линии тяги пахотного агрегата, составленного из колесного трактора МТЗ-82 и навесного плуга ПИН-4-35. Идентификация проведена с использование приложения Simulink® Parameter Estimation ППП MATLAB®.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Идентификация параметров динамических объектов посредством приложения Simulink пакета Matlab»

Технологии и технические средства механизированного производства продукции растениеводства и животноводства.

ЗАКЛЮЧЕНИЕ

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

Эффективность использования ИТ в с.-х. производстве будет достигаться за счет:

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

- оптимизации использования ресурсного обеспечения хозяйства, рационального распределения ресурсов по отраслям производства;

- принятия оптимальных решений при изменении условий производства и рыночной конъюнктуры цен на производимую продукцию;

- оптимального выбора технологий производства;

- оптимизации структуры машинотракторного парка.

Получено 17.03.2008.

УДК 631.531.17-52:633(470.31)

В.Н. БРОВЦИН, д-р техн. наук

ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ ДИНАМИЧЕСКИХ ОБЪЕКТОВ

ПОСРЕДСТВОМ ПРИЛОЖЕНИЯ SIMULINK ПАКЕТА MATLAB

В статье приведен пример идентификации параметров пахотного агрегата по каналу управления положением линии тяги пахотного агрегата, составленного из колесного трактора МТЗ-82 и навесного плуга ПИН-4-35. Идентификация проведена с использование приложения Simulink® Parameter Estimation ППП MATLAB®.

Наиболее часто используемые методы идентификации параметров моделей сельскохозяйственных процессов пригодны, как правило, для линейных по параметрам математических описаний, и, кроме того, распределенных по нормальному закону возмущающих воздействиях [1] (или имеющих специальный вид, например, единичный

196

ISSN 0131-5226. Сборник научных трудов.

ГНУ СЗНИИМЭСХРоссельхозакадемии. 2008. Вып. 80.

скачок). Тем не менее, для сельского хозяйства характерны существенно нелинейные процессы, экспериментальное исследование которых по существующим методикам организовать сложно, а часто и невозможно, в связи с чем целесообразно использовать метод, свободный от отмеченных ограничений - вариационный метод решения экстремальных задач. Метод позволяет провести идентификацию параметров моделей или модельных коэффициентов при любом удобном для исследователя наборе экспериментальных данных. Программная реализация метода представлена в приложения Simulink® Parameter Estimation ППП MATLAB®.

В качестве примера проведем идентификацию параметров пахотного агрегата по каналу управления положением линии тяги (рис. 1).

Рис. 1. - схема плуга:

е - смещение линии тяги, см; РТ - тяговое сопротивление, кН; Rex - усилие на бороздном колесе вдоль линии тяги, кН

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

197

Технологии и технические средства механизированного производства продукции растениеводства и животноводства.

Рис. 2. Зависимость тягового сопротивления плуга и усилия на бороздном колесе от смещения линии тяги на супесчаной почве [1]:

1 - тяговое сопротивление, кН; 2 - усилие на бороздном колесе

Из графика (см. рис. 2) видно, что с увеличением ширины захвата увеличивается тяговое сопротивление плуга. Также из графика видно, что при оптимальном положении линии тяги плуг имеет минимальное значение тягового сопротивления и максимальное значение усилия на бороздном колесе. Со смещением линии тяги в сторону от оптимального наблюдается увеличение тягового сопротивления плуга и уменьшение усилия на бороздном колесе. Это объясняется тем, что при смещении линии тяги в сторону от оптимального нарушается устойчивость плуга, тяги механизма навески трактора перекашиваются, и трактор начинает «уводить» в сторону вспаханного поля [1]. Если принять, что допустимое увеличение удельного сопротивления по потенциальной эксплуатационной характеристике должно быть не более 3-6%, то из графика (см. рис. 2) видно, что предел изменения линии тяги, при котором сохраняется устойчивое положение плуга находится в довольно узком диапазоне 50-100 мм, поэтому при работе необходимо добиваться более точной установки положения линии тяги. Для

198

ISSN 0131-5226. Сборник научных трудов.

ГНУ СЗНИИМЭСХРоссельхозакадемии. 2008. Вып. 80.

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

В связи с тем, что зависимости РТ и Лбх как функции смещения линии тяги имеют ярко выраженный унимодальный характер, целесообразно в качестве управляющего устройства использовать экстремальный регулятор, для синтеза которого необходимо разработать модель по каналу управления положением линии тяги и определить численные значения ее параметров. С этой целью были проведены экспериментальные исследования на супесчаной почве при скорости движения агрегата 2,5 м/с с дискретностью записи 0,05 с. В качестве входного воздействия записывали величину напряжения (и,в) на обмотке электромагнитного клапана (см. рис. 4, а) , а выходными - движение штока (Ь,см) гидроцилиндра (см. график 1 на рис. 4, б) и усилие (Лбх, кН) на бороздном колесе (см. график 1 на рис. 4, в).

При определении структуры примем во внимание, что порядки передаточных функций по каналам управления для плугов, культиваторов и лущильников my, пу<2 [2].

Предположим, что динамические свойства объекта Wo6 (усилие

на бороздном колесе - положение штока гидроцилиндра) описывается выражением

Wоб (s)

1

г2 s2 + Ts + г

(1)

где Т1 и Т2 - параметры, с,

а экстремальная характеристика процесса (нелинейное звено) -простой параболой.

Рб = аЬ2 + bL + c, или Рб = Рб0 + А(Ь-Ь0)2, (2)

где А = а, Рб = -Ь/2а, Ь0 = (4ас - b2)/4a.

Передаточную функцию электромагнитного клапана запишем

в виде:

WK (s) = —^, (3)

к T +1

где Тк - обобщенная постоянная времени электромагнитного клапана и маслопровода, с.

199

Технологии и технические средства механизированного производства продукции растениеводства и животноводства.

Движение штока гидроцилиндра соответствует интегрирующему звену:

Wg (s) =-k (4)

Tg

где Tg - постоянная интегрирования гидроцилиндра, с; Tg = Tgs при движении штока в прямом направлении и Tg = Tgi - в обратном.

В выражениях (1), (2) и (3) коэффициенты Т1, Т2, Тк, A, L0, Tgs, Tgi и R0 являются параметрами объекта управления и подлежат определению по результатам обработки экспериментальных данных.

Модель по каналу управления S, в соответствии с уравнениями (1) - (4), представлена на рис. 3.

Рис. 3. Модель по каналу управления положением линии тяги

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

200

ISSN 0131-5226. Сборник научных трудов.

ГНУ СЗНИИМЭСХРоссельхозакадемии. 2008. Вып. 80.

Рис. 4. Результаты идентификации параметров пахотного агрегата:

а - процесс изменения напряжения на обмотке электромагнитного магнитного клапана; б - смещение линии тяги;; в - усилие на бороздном колесе

1 - экспериментальные данные; 2 - идентифицированные данные

Итак, задача заключается в том, чтобы на решениях S-модели получить численные значения идентифицируемых параметров, обеспе-

201

Технологии и технические средства механизированного производства продукции растениеводства и животноводства.

чивающих при заданных на них ограничениях минимальное значение квадратичного критерия качества.

Для решения задачи идентификации использовали нелинейный метод наименьших квадратов (Nonlinear least squares), реализованный в алгоритме Гаусса-Ньютона (Gauss-Newton), программная реализация которого представлена в приложении Simulink® Parameter Estimation ППП MATLAB®.

Результаты идентификации параметров технических средств канала управления положением линии тяги представлены на табло S-модели (см. рис. 3).

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

Малые возмущения начальных условий приводят к тому же результату, что подтверждает корректность (адекватность) модели экспериментальным данным [3].

СПИСОК ЛИТЕРАТУРЫ

1. Степанов А.Н. Повышение эффективности вспашки путем пользования плугов с изменяемыми параметрами /А.Н. Степанов: Дис. ... канд. техн. наук. - СПб. - Пушкин, 1999.

2. Лурье А.Б., Еникеев В.Г., Теплинский И.З., Смелик В.А. Сельскохозяйственные машины (машины для обработки почвы, посева, посадки, внесения удобрений и хим. защиты растений) /А.Б. Лурье и др. - СПб.: изд-во СПбГАУ, 1988.

3 Бровцин В.Н. Исследование и оптимизация динамических объектов сельскохозяйственного назначения средствами вычислительного эксперимента /В.Н. Бровцин. -СПб.: СЗНИИМЭСХ, 2004. -ISBN 5-88890-027-3.

Получено 24.03.2008.

202

i Надоели баннеры? Вы всегда можете отключить рекламу.