автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование и выбор параметров механизмов в комплексе с приводными системами

10-7 2215

На правах рукописи

Мисюрин Сергей Юрьевич

Математическое моделирование и выбор параметров механизмов в комплексе с приводными системами

Специальность 05.13.18 - Математическое моделирование,

численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени док тора физико-математических наук

Автор

Москва-2010

Работа выполнена в Национальном исследовательском ядерном университете «МИФИ»

Научный консультант: - доктор физико-математических

наук, профессор Кудряшов H.A., НИЯУ «МИФИ»

Официальные оппоненты: - доктор физико-математических

наук, профессор Михайлов В.Б., ЦИТП РАН

- доктор технических наук, профессор Гогричиани Г.В., ОАО ВНИИЖТ

- доктор физико-математических наук, профессор Борисов A.B., УдГУ

Ведущая организация - Учреждение Российской академии

наук Институт радиотехники и электроники им. В.А.Котельникова РАН

Защита диссертации состоится «29» декабря 2010 года в 15 часов 00 мин на заседании диссертационного совета Д212.130.09 в Национальном исследовательском ядерном университете «МИФИ» по адресу: 115409, Москва, Каширское шоссе, д.31.

С диссертацией можно ознакомиться в библиотеке НИЯУ «МИФИ»

Автореферат разослан «Л £"""» ноября 2010 г.

Ученый секретарь диссертационного совета, доктор физико-математических наук, профессор

Леонов A.C.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

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

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

Научная новизна работы.

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

Определены базовые критерии подобия приводной системы.

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

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

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

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

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

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

• Практическими рекомендациями по применению цифровой техники для реализации сложной позиционной системы управления, в том числе на базе модели-наблюдателя (оценка

необходимого быстродействия, квантование управляющих сигналов по уровню и времени).

В ходе выполнения работы была разработана и внедрена «Методика расчета и анализа процессов движения и управления позиционированием объекта приводом с гидравлическим двигателем (программа gidroprivod)». Акт о внедрении представлен Институтом прикладной физики Российской академии наук (г. Н. Новгород), Нижегородский технический университет.

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

Использование модели-наблюдателя в системе управления пневмоприводом позволило обеспечить натяжение бумажной ленты в печатной машине NUR Expedio 3200 (HP) с относительной погрешностью 0,5%. Предложенным решением была упрощена конструкция и уменьшена ее стоимость.

Основные научные результаты, защищаемые автором.

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

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

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

• Установлена относительно малая чувствительность динамики привода к вариации кинематической модели механизма (вида его передаточной функции).

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

• Решена проблема анализа системы уравнений связей механизма в окрестности особых положений посредством построения его функции положения в явном виде в форме степенных рядов.

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

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

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

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

Апробация работы. Основные результаты диссертации докладывались на следующих симпозиумах, конференциях семинарах:

Международная конференция «SPATIAL MECHANISMSND HIGH CLASS MECHANISMS» (Theory and Practice), October 04-06, Almaty;

конференция «Новые технологии инновационному бизнесу» (первый российский форум). Москва, 09-20 февраля 2007 года;

конференция «Проблемы динамики и прочности исполнительных механизмов и машин». Астрахань 10-16 сентября 2007 года;

VIII Международная научно-техническая конференция по динамике технологических систем. Ростов-на-Дону, 09-13 октября 2007 года;

VI Международная научно-техническая конференция «Материалы и технологии XXI века». Пенза, 2008 год;

Международная научно-техническая конференция «Наука и образование - 2008». Мурманск, 2008 год.

Публикации. Основное содержание диссертации опубликовано в 40 статьях и докладах на конференциях, а также в более чем 10 научно-технических отчетах по исследовательским работам, выполненным по целевым программам и грантам РФФИ.

Структура и объем работы. Диссертация состоит из введения, пяти разделов, заключения и списка использованных источников. Основное содержание изложено на 235 страницах машинописного текста и содержит 62 рисунка и 8 таблиц. Список источников включает 206 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Введение

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

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

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

В литературе известно большое количество исследований, посвященных анализу, синтезу и математическому моделированию отдельно механической части привода, отдельно двигателей, отдельно систем управления. Большой вклад в решение проблем кинематики, динамики и синтеза плоских и пространственных механизмов был внесен известным ученым П.Л.Чебышевым, академиками И.И. Артоболевским, Н.Г Бруевичем, а также профессорами В.А. Зиновьевым, Ф.М. Диментбергом, А.Ф. Крайневым, В.А. Глазуновым, П.А. Лебедевым, G.J. Bennet'oM, С. Gosselin, J. Angeles, K.J. Waldron, V Parenti-Castelli, C. Innocenti и многими другими. Наряду с созданием графических графоаналитических методов исследования и синтеза механизмов они много внимания уделяли развитию аналитических и, что немало важно, численных методов исследования. Значение последних весьма велико, так как они дают возможность при помощи математического моделирования использовать современные вычислительные машины с программным управлением для решения практических задач.

В области теории двигателей (электрических, гидравлических, пневматических или других) мы пользуемся результатами, полученными многими авторами. Это направление представляют работы Е.В. Герц, Г.В. Крейнина, Е.А. Цухановой (Институт машиноведения РАН), Д.Н. Попова, В.М. Прокофьева (МВТУ им. Н.Э. Баумана), Б.И. Петрова, Н.С. Гамынина, A.B. Рабиновича, В.В. Саяпина (МАИ), Б.К. Чемоданова (МЭИ), а также зарубежных ученых W. Backe, В. Andersen, К. Araki, W. Frank, F. I. Takemura и многих других. Перед нами стояла задача выделить общие свойства двигателей различных типов, построить их обобщенную модель и сформировать критерии подобия, которыми можно было бы удобно пользоваться в решении задач синтеза приводной системы в целом.

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

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

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

Постановка задачи.

В первом разделе сформулирована основная задача диссертации разработка метода математического

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

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

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

синтеза состоит из нескольких этапов, рассматриваемых ниже во втором разделе.

| Механизм

^ Передаточная | функция:

| г = х!у ТТ./7777777 -

Рис. 1

Отметим, что в большинстве позиционных систем (дорожных машин, транспортирующих, погрузочно-разгрузочных и др. устройств, в том числе и в примере на рисунке 1, двигатель передает движение на объект через шарнирно-рычажный (или другой) механизм), имеющий переменное передаточное отношение / = х/у (где х скорость двигателя, у скорость перемещаемого объекта). Механизмы, у которых передаточное отношение есть величина постоянная, являются частным случаем.

В качестве примера механизма с переменным передаточным отношением приведем механизм, изображенный на рисунке 2.1а. Этот механизм предназначен для поднятия груза массой т на некоторую высоту Утах

В качестве входного звена можно выбрать одно из звеньев, соединенное с основанием (например, звено АС). В этом случае при движении механизма возникает особое положение (в момент, когда звено АР1 находится на одной прямой со звеном БЭ рис. 2.16). При параметрическом синтезе таких механизмов возникают дополнительные трудности, связанные с тем, что Якобиан системы уравнений связей обращается в ноль. Подробнее это проблема представлена в разделе 3.

Рис 2.1

Другим примером исследуемого класса механизмов может служить привод, изображенный на рисунке 2.2. Это схема действующего привода ПСДС 4, используемого в газопроводах для перекрытия подачи газа. Здесь 1 - блок управления привода, 2 -струйный двигатель, работающий за счет сжатого газа. Его вращение происходит под действием реактивной силы сопла, изображенного на рисунке 1-6. Выходным звеном механизма (привода) является вал 8. Передаточное отношение 1-х/у есть отношение скорости вращения двигателя 2 к скорости вращения вала 8. Оно зависит от параметров редуктора 6 и кулисы 9 и является величиной переменной. Таким образом, мы имеем механизм (6, 9) - двигатель (2) - управление (1).

12 3 4

Рис.2.2

На рис. 3 представлена общая расчетная схема позиционной системы, которая состоит из следующих подсистем: двигатель 1, распределитель 2, управляющее устройство 3 и передаточный механизм 4. С выходным звеном механизма связан перемещаемый (позиционируемый) объект 5. \

Простейшие математические модели позиционного механизма.

Во втором разделе представлены простейшие математические модели системы «механизм-привод-управление». Механизм заменяется его передаточной функцией. Вводится понятие идеальный двигатель двигатель, двиэ/сущая сила которого однозначно определяется только сигналом управления.

На первом этапе механизм 4 (рис. 3) представлен относительно простой типовой функцией передаточного отношения отношением скорости входного звена к скорости выходного i = х/у, а двигатель идеальной моделью, которая не учитывает зависимость движущей силы двигателя от скорости или других факторов. Предполагается, что «идеальный двигатель» мгновенно отрабатывает сигнал управления, подаваемый управляющим устройством. Получаемая таким образом простейшая математическая модель привода позволяет приступить к выбору в первом приближении закона управления с учетом взаимодействия всех подсистем.

Функция i в общем случае является переменной, зависимой от положения механизма; она может быть задана как функция г'(х) или i(y) В дальнейшем принимается, что передаточная функция задается в виде i(y). Если i = const (что характерно для зубчатых, винтовых и некоторых других механизмов, например, транспортера, изображенного на рис. 1-4), то имеет место соотношение i = s..lsv, где sY и sv взаимосвязанные

Л. ) Л у

перемещения двигателя и объекта. В позиционном приводе с ограниченными ходами sx полный ход двигателя, sy полный

ход объекта.

В случае i - var отношение полных ходов механизма на входе и выходе sx Is — ie называется эквивалентным

передаточным отношением и если функция i — i(y) известна, то

1 ^

\i(y)dy

^ л

где уК2 - координаты начального и конечного положений выходного звена механизма, Аяу — у2 — У\ - его полный ход.

Замена механизма передаточной функцией дает возможность выписать систему уравнений всей системы, взяв в качестве звена приведения перемещаемый объект ш:

ту = Ру + Рг

Р = Р

А у л X

(1)

где Рх движущая сила двигателя; Рг сила сопротивления, действующая на объект.

Если управляющее устройство следит непосредственно за движением объекта (по параметрам у^ и передаточное отношение задано в виде ¡(у), то при описании динамики системы к уравнениям (1) следует добавить уравнения изменения движущей силы двигателя под воздействием управляющего устройства

(2)

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

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

1) Г = К](уе-у)-К2 у-К3 у,

2) у = Кг 81§п(уе - уЦуе -у\-К2 у-К3 у (3)

где К], К2 и К3 - коэффициенты усиления в контурах обратной связи по положению, скорости и ускорению;

уе - координата конечного положения объекта.

Закон 1) это простейший закон пропорционального управления по трем параметрам состояния; однако вследствие ограничения \у\ < 1 здесь присутствует эффект насыщения, т.е. этот закон в общем случае является существенно нелинейным.

Что касается закона 2), то в нем, кроме эффекта насыщения, присутствует нелинейная обратная связь по положению.

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

У ~ ^ у ' V Ру-Ру Ху > Ру ~ Рх тах ' К >

Р.х ~ Рх тах ' У ' Р>- ~~ Ру Хг"> t — t -Т ,

/ = ¡е I, где /* = (шу /Ру)0'5 £ - максимальное перемещение входной координаты, Л-тах - максимальная движущая сила на входе в механизм.

Уравнение (1) в безразмерной форме переменных имеет вид:

(3)

\у = Ху + Хг

[Ху=Г 1

Законы управления (3) аналогично:

1) У (уе -у)-к2 \>-КуУ,

2) у = /с1 - У)л\\Уе-У\-К2'У-КЪ-У

где к2=(8у1?)- К2, къ = {зу/Р)- К3

Уравнения (1*) и (3 ) мы называем моделью первого уровня. В табл.1 представлены процессы, характеризующие зависимость динамики привода от характера изменения передаточного отношения и от структуры и параметров системы управления. Данные зависимости были получены путем численного решения системы 1* В качестве базовых выбрано четыре закона изменения передаточного отношения: 1) постоянное; 2) V— образное; 3) падающее; 4) возрастающее. Во всех четырех случаях было принято ^ =1,5, что эквивалентно переменной нагрузке в виде жидкостного трения х,- — (около половины максимальной величины движущей силы); минимальное значение безразмерного передаточного отношения = 0,5

Анализ представленных в таблице процессов позволяет сделать следующие выводы.

Время процесса относительно мало зависит от характера изменения I. Некоторое увеличение времени отмечено только при

ользовании механизма с передаточной функции вида 4, что зано с замедлением разгона. Время процесса заметно ращается при переходе от закона управления 1) к закону авления 2).

Таблица 1

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

Вид

передаточной Закон управления 1) Закон управления 2) функции I_

1 1 1 ■>

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

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

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

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

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

Модель привода второго уровня отличается от модели первого уровня включением в систему уравнений зависимости, которая учитывает влияние эффекта запаздывания в отработке сигнала управления распределительным устройством:

У = (Уте)-(Ус-У)

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

хе - постоянная времени распределительного устройства.

Модели следующих уровней строятся по нарастающей сложности в зависимости от типа двигателя. Более подробно это представлено в четвертом и пятом разделе.

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

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

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

Проблема изучения кинематики, кинетостатики и динамики механизмов тесно связана с проблемой функции положения механизмов, получаемой из системы уравнений связей механизма, Эта система имеет вид:

'/\(21,г2>->2п>У1>У2>->Ук>1\>12>->1т)=0

/2(2,, г2,...,г„, Ух,Уг,~>,Ук,1\>12>-~>1т) = 0

(4)

где координаты входных звеньев механизма;

У\>Уг>—>Ук ~ координаты выходных звеньев механизма; /1,/2,...,/т - параметры механизма.

Уравнение (4) можно также записать в векторном виде: Цг,у)=0 (4*)

Известны разнообразные методы представления уравнений связи. Например, метод В.А. Зиновьева, основанный на замкнутом контуре векторов, который выражает неизменность длин звеньев и соотношений углов. Метод Ф.М. Диментберга, основанный на

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

У\ =<Р](г]>22>->гп>1\>12>->1т)

I ^

Уг =

или у = Э(гьСистема уравнений (5) получается решением системы (4).

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

С проблемой изучения свойств функции положения механизма тесно связана проблема его особых положений. При движении плоские и пространственные механизмы как с одной, так и с несколькими степенями свободы могут оказаться в близи особого положения. С математической точки зрения особое положение характерно тем, что Якобиан системы уравнений связей механизма обращается в нуль. Автором данной работы впервые была разработана методика получения функции положения механизма в окрестности особого положения в явном виде путем представления функции в виде степенных рядов:

ук = ук° + Аук = ук° + ^ Д-Д^' где ук - координаты звеньев

механизма, определяющие его особое положение.

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

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

Здесь звено ОА (либо угол а) является входным звеном, звено ОВ (либо угол Р) выходным, пары О, А, В, И - вращательные, 00 - основание четырехзвенника. Уравнения связей механизма имеют вид:

х2+/-/,2 = О

- х)1 + (у} - у)2 -¡2=0 (хх-а)2+(ух-0)2-1% = 0

где 0(0,0), А(х,у), В(хьу!) Э(а,0) координаты шарниров механизма в декартовой системе координат; и, Ь> Ь длины звеньев ОА, АВ, ВЭ.

Приравняв к нулю Якобиан уравнений связей, находим особое положение 1-го и 2-го рода (рис. 4). Нас интересует положение второго рода, так как при этом положении система уравнений связей имеет неоднозначное решение. По разработанной методике получаем точное представление функции положения механизма в окрестности особого положения 2-го рода.

х

2

9-

У

Рис.4

Приведем их вид:

Аух = Ах Ах + А2Ах2 + А3 Ах3

Ах, = ВхАх1/2 + В2Ах + В3Ах

Ау] = С^Дх175 + С2Ах + А3Ах где

,1/2

.3/2

.3/2

Т2у°хТх-ВУ{х°+у°Ах) у°{у°Вх-СхТх)

_ У\Т2~У°С\С2

с _ т2вх+итх в =и-дс:

3 С{гх-у°в; 3 вх Т\ - х° - а, Т2=В1+А1С1,

и = 2у0х{В1 + А1у°х+С1Ах) в2 ^

3 где

2 2 .Г

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

Будем предполагать, что в механизме при подходе к особому положению скорость и ускорение малы, и ими можно пренебречь (0 = 0,иф = 0). Также мы не будем учитывать силы трения. Динамическое уравнение с учетом принятых допущений:

.. эп Л , м ш .

М--= 0 либо —--= 0 (7)

д (р у дх

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

у У + У\ У\

П =

где Ш1 т2 , т3 массы звеньев механизма; % постоянная величина. Продифференцировав выражение потенциальной энергии по обобщенной координате х, имеем.

дП дх

ду, Л ду\ , ч

— (/и, +т2) + —1-(т2 +т3) дх дх

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

дП дх

(/и, + т2 + 2 А2 Ах +...) + (т2 + тъ Сх Дх"1'2 + С2 +...)

После подстановки этого выражения в (7) получим уравнение следующего вида:

М + А(Дх) = 0 (8)

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

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

{тх +т2)А1 + {т2 + т^С.Ах^2 +С2)1|

дП

дх

Подставив это выражение в (6), получим следующее уравнение для определения «зон заклинивания» механизма:

д~= _(т2+т3 )С,_

|_2(2М /(Л)" (Щ + т2)Ах ~{т2 + тъ)С2) Другими словами, при ограниченном моменте М на входном звене ОА механизм не сможет подойти к особому положению ближе, чем на Дх. Все остальные координаты звеньев механизма вычисляются аналогично. Пользуясь полученным соотношением, можно вычислить, насколько следует увеличить момент, чтобы «зона ограничения» уменьшилась до требуемых размеров. Тот же результат можно получить уменьшением масс т2 и тъ. Если записать уравнения с учетом сил трения, то величина Дх будет несколько больше.

II. Решая поставленную задачу во втором приближении, т.е. отбросив все члены порядка выше второго, имеем:

— = [(#и, + т2 + (т2 + т3 )R2 )]£, дх 2

Щ =А1 +2Л2Ах, R3 = (1/2)С,ДЛГ1/2 + С2+(3/2)С3Дх1/2 + 2С4ДХ.

Подставив это выражение в уравнение (6) и умножив его

1/2

предварительно на Дх: , получим:

7j + Т2 Ах + Т3/ л[Ах + Г4 .л/Дх = 0, либо

Тх + Т2(VAX)3 + Т3+Т4 .(VAjc)2 - 0,

где 7] = 2М f(gy°) - (/77, + т2)Л, - (m2 + m3)C2,

T2 = -2A2(ml+m2)- 2C4(m2 + m3) + Mc° /(/ f2,

Г3 = ~(m2 +/я3)С,/2, Г4 = -(/w2 + m3)C3.

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

Заметим, что в обоих случаях решением для Дх будет положительное число. Как видно из рисунка, переменная х может отличаться от положения равновесия (х°) только на положительное значение.

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

Особое положение соответствует следующим значениям:

х° = 0,5, /=л/3/2, *,°=5/4, = у0 /2 = л/3/4

Пусть масса звеньев т{= тъ Приведём значения Ах для

различных величин т2(0,5 и 1 к^) и момента М (15, 20, 30,40 н).

^\М(н) 15 20 30 40

0,5 Дх =0,25 Дх =0,09 Дх =0,027 Дх =0,013

1 Ах =0,09 Дх =0,038 Дх =0,018 Дх =0,006 7

Как видно, размеры «зоны ограничения» вполне соизмеримы с размерами механизма.

Моделирование приводной системы с гидродвигателем.

В четвертом разделе на конкретных примерах показана полная процедура моделирования и последующего синтеза приводной системы «механизм-привод-управление» с использованием в качестве двю/сущей силы гидропривод. Моделирование строится по нарастающей сложности.

Рассматривается задача определения параметров технической системы, состоящей из гидродвигателя, механизма передачи движения от двигателя на перемещаемый объект и устройства управления. Механизм имеет одну степень свободы и переменное передаточное отношение. Система предназначена для быстрого перемещения объекта из одного положения в другое. Перемещаемый объект характеризуется массой и постоянной силовой нагрузкой, создаваемой, например, весом поднимаемого (опускаемого) объекта.

Постановка задачи: требуется определить основные параметры системы, перемещающей объект массой т, нагруженный силой Р>у, за время гс на величину хода с остановкой в конечной точке с заданной точностью, которая оценивается параметром отклонений уА0. К искомым параметрам относятся размеры двигателя и распределителя, давление питания

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

Моделирование передаточного механизма. На рис. 3 представлена общая расчетная схема позиционной системы (двигатель - 1 , распределитель - 2, управляющее устройство - 3, передаточный механизм - 4). С выходным звеном механизма связан перемещаемый (позиционируемый) объект 5.

Как было показано во втором разделе, на первом этапе передаточный механизм заменяется его простейшей кинематической моделью передаточной функцией механизма, выражающей закон изменения отношения скоростей на его входе и выходе в функции перемещения; например, выходного звена i = i0K6 I(у), где 1{у) - функция, пропорциональная i(y), которая равна единице при i = iDKe. Несколько функций 1{у) различного вида показано на рис. 6.

Определение движущей силы двигателя. При комбинированной нагрузке механизма (масса + сила сопротивления) уравнение движения его выходного звена

■Pry, (9)

Рис. 6

Функция 1(у) изменения передаточного отношения механизма в зависимости от перемещения объекта: 1-кусочно-линейная; 2- параболическая; 3-синусоидальная.

где т - масса перемещаемого объекта;

Р1у силовая нагрузка на объект, которая может быть как

положительной (например, при опускании груза), так и отрицательной (при подъеме груза);

2 - эффективные площади поршня гидроцилиндра (ГЦ) со

стороны напорной и сливной полости;

2 -давление в соответствующей полости.

my = iPi Р\~Рг Рг)'*экв 1

2 ,

3 - 1 •I

У1 У2

В общем случае (например, в ГЦ с односторонним штоком) Ф • Обозначим Т7 собственно площадь сечения ГЦ, а также введем относительные показатели эффективных площадей к\,2 ~ -^1,2 / ^ Уравнение статического баланса приложенных к объекту движущих сил и сопротивления получим из уравнения (9), положив у- 0. Одновременно представим силу сопротивления в виде Ру = 8,у рр Т7, где - сила сопротивления, выраженная в

долях от максимальной движущей силы ГЦ, которую определим как рр Т7, где р -давление питания. В результате имеем:

(Р\ Р\~Рг рг) ¿жв+дгу рр ^ = 0

или(6гк1-д2'к2)-1экв I + 51у = 0, (10)

где б] 2 =р12/рр.

При выборе величин <5 к:^,/жв следует учитывать, что силовая

нагрузка должна преодолеваться как в динамике, так и в статике. Поэтому должно выполняться условие

оке

(11)

Примем, что при 7 = 7т1п (наихудшие условия нагрузки на двигатель) указанное неравенство выполняется, если относительная движущая сила ГЦ <51 */с, —62 к2> 0,9; предельная величина этой силы равна единице. Тогда получаем следующее ограничение на выбор закона изменения 1(у) передаточного отношения механизма в зависимости от выбранных значений 5,у,1экв

> 1,11 • Д где Д = \5)у | / 1экв. В точке, где / = ¡экв, имеем 7 = 1. По выбранной с учетом других параметров величине находим Р = Ргу/51урр, предварительно задавшись давлением питания рр.

Как показало предварительное исследование, величину /т/Р\,2 или гДе /т~ эффективная площадь проходного

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

имеет передаточное отношение / = ¡экв (т.е. 1=1); скорость поршня постоянна и определяется при условии, когда в течение всего периода движения полностью открыты каналы распределителя входной в напорной полости и выходной - в сливной. Для описания движения этой упрощенной системы воспользуемся уравнением (11), которое для случая подъема груза (5^ < 0) можно

представить в виде

(Р\ Р\~Р2 р2)-1жв =\3,у\-РР Е (12)

К уравнению (12) следует добавить уравнения, описывающие процессы поступления жидкости в рабочую полость и отвода жидкости из сливной полости, которые определяют установившиеся давления в полостях ГЦ:

Р\ = Рр~2\\Р2 = Рр~ г г\

2 2 <13)

/] = (^1 *уст /Я) '>22 = 2~(Р2Хуст

где 2 = рр -ра; X — /т(2/р);

ра- давление в сливной линии; р - плотность рабочей жидкости. Решая систему уравнений (11) и (12) с учетом хуст — • ¡экд / (величины бу и 1С предполагаются заданными), будем иметь:

3 3

^ = (Зу *э* //с) ■ (-Р(К] +К2\ I-)0,5 (14)

р 2 Рр^х-Кг'да-^гуугэкв)

Далее, зная ^ = Р /\д1У\рр, определяем эффективную площадь /т каналов распределителя. Как уже указывалось выше, величина определяемая по соотношению (14), является первым приближением. Она может быть скорректирована по результатам моделирования динамики системы с учетом процессов изменения передаточного отношения механизма, открытия окон распределителя, влияния сжимаемости жидкости и других факторов. Для этих целей далее используется иерархия динамических моделей нарастающей сложности. Делая замену переменных

1 = у = у* V, х = у* р]2 = р*.3]2, (15)

принимая масштабы измерения безразмерных переменных р*=рр, Г =(тзу/РрР)°>5 (16)

а также учитывая, что масштабом измерения передаточного отношения механизма является 1экв (т.к. г = ¡экв /), а масштабом измерения у перемещения распределительного органа

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

у = {бх к^-52 К2) 1экв I + 81у;

¿>, =1-А;

с п * ч /А + -и2 =(1"^Н-^-}

a) при -и2-к142 §2 = \-32;

62(-^-)

b)при-и1 -/с2-<*2 <0 52=\ + 82-¿ГТ;--}

1 = тэкв=ну),£ = у- /,

у = -(1/(2УтУь)°'5) • (У + (21/,>|)0'5), (17)

где: За=ра1рр; 5,у = Р1у 1рр¥\уь = уь Iбу\ут =ут-(т/ррР) у0.(У) = fei /fe¡ - текущее значение эффективной площади / -го канала, /т- его максимальная эффективная площадь при полном открытии. Закон управления у оптимизирован применительно к задачам позиционирования объекта; он подробно рассмотрен в предыдущих двух параграфах. Величина = у^/

входящая в выражение у, характеризует относительную ширину

у1 области между линиями насыщения на фазовой плоскости «перемещение-скорость» выходного звена. Чем эта область уже, тем ближе подходит фазовая траектория к линии переключения, что способствует повышению точности отработки заданной позиции. С этой точки зрения желательно иметь величину уь в конце движения возможно меньшей. С другой стороны, если задать величину у1 малой на всей длине хода, то управляемость системы падает. В частности, затрудняется переход в режим торможения, который может оказаться достаточно резким. Кроме того, как будет показано ниже, необходимость выхода в зону активного управления (между линиями насыщения) часто возникает с самого начала движения. Поэтому ширина зоны активного управления была выполнена переменной (в форме «раструба») с сужением к концу хода:

^ "И» (18)

где у- ширина этой зоны соответственно в конце и в начале движения системы. Параметр

Р Р Sv

и = (Р//т) „ (19)

Р т (1-За)

является характеристикой инерционности системы. Подставив в (19) полученное ранее выражение (16) для /т /Т7 и

F = Ргу /\5)у |рр — а т- g /\5)у \рр, получим окончательно

и = (/с / 10Кв ){аё( 1 + / ¡жв - 5а) /(2зу |<?,7|(1 - ))°'5 (20)

где а - Р1у /т • g- сила сопротивления, выраженная в долях от

силы веса перемещаемого груза. Для случая подъема или опускания груза а = 1.

Если закон изменения передаточного отношения механизма имеет У-образную форму, то при выбранной величине максимальное значение тангенса угла наклона прямых /(г) к оси V будет Р — 4 • (1 — /т;п), (21)

где /т;п = 1,11 А. В результате получим законы изменения передаточного отношения механизма:

для - 1 < V < -0,5 1 = 1тт-р-(у + 0,5) (22)

для - 0,5 < V < 0 1 = 1т[п + р • (V + 0,5)

При использовании системы (16, 17) для моделирования динамического процесса вводятся исходные данные т,Рг„Аг,БЛ).

гу С )

Параметры ^Экв'^гу>ут»уЬ0'уы являются варьируемыми и выбираются из условия реализации поставленных требований к качеству процесса движения, к точности отработки и удержания заданной позиции, к структуре и параметрам передаточного механизма.

Процедура численного выбора параметров системы. Рассмотрим процедуру выбора параметров системы на конкретном примере. Заданы следующие значения исходных параметров: ¿С.=2с; 0,5 м; т- 500кг; а = 1 (подъем груза); задается

рр =4-106МПа, что соответствует 5а =0,025; д1у = -0,1 и -0,15 (два варианта); 1экв =0,2; ГЦ с ходом 5 = 0,2 • 0,5 = 0,1 м

имеет односторонний шток (при подъеме кг, = 1, кг2=0,8). Соответственно при опускании груза будем иметь 8,у - 0,1 и 0,15; кг, = 0,8, к2 — 1 •

Из соотношения (13) следует, что для реализации процессов подъема и опускания груза за одинаковое время необходимо иметь различные по величине /т /Т7 поскольку входящая в (13) относительная нагрузка <5 имеет знак "-" при подъеме груза и знак "+" при опускании. Различие в величинах /т / .Р при движении объекта вверх и вниз может возникнуть также по причине несовпадения требований к скорости подъема и опускания груза. Поэтому далее, учитывая, что при прочих равных условиях скорость движения груза вниз выше скорости подъема, отношение /т / Т7 выбирается по условиям движения груза вверх, а время опускания груза регулируется выбором значения входящего в (18) и определяющего ширину «раструба» в начальной позиции объекта. Расширение «раструба» означает увеличение зоны активного управления положением распределительного органа, что при достаточной ширине этой зоны вызывает частичное перекрытие каналов золотником с самого начала движения.

Пользуясь (13), (19), а также соотношением между величинами и Т7 получим для движений груза вверх

следующие значения /т / Г, Р, / и и

Таблица 2

5 гу ¿т Р,м2 Г 2 и

-од 0,000944 0,0123 13,4 80 30,2 2,14

-0,15 0,00136 0,0082 13,0 45 26,1 2,17

Далее пользуемся системой уравнений (16)-(18), а также (21)-(22) для моделирования процесса подъема груза, протекающего без учета сжимаемости жидкости. Кривые 1 и 2 на рис.7,а представляют характер изменения безразмерных скоростей объекта

и поршня ГЦ для 1^1 = 0,15 в функции их перемещений. При моделировании процесса было принято: у10 = 0,04 и уп = 0 (т.е. ширина «раструба» остается постоянной на всем участке движения); 1экв = 0,2 (ход поршня ГЦ в пять раз меньше хода объекта); ут = 0,25. В табл. 2 представлены также величины безразмерного времени движения тс и соответствующие им значения реального времени /с, которые достаточно мало отличаются от заданной величины /С=2с. При подсчете ¿с были использованы соотношения (14) и (15) между размерным и безразмерным временем.

Кривые 3 и 4 на рис.7,а представляют безразмерные скорости объекта и поршня ГЦ, полученные при моделировании процесса

опускания груза (также для относительной нагрузки = 0,15). Соответствующие значения параметров записаны ниже

Таблица 3

$гу /С,с

од 2 29,8 2,12

0,15 5 23,8 1,98

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

Для перехода от безразмерных скоростей и перемещений к размерным пользуемся приведенными выше зависимостями, из которых следует: х = 5,72 • £;у = 5,72 • V;х = 0,5 • £; у = 0,5 • у

Принятый за основу кусочио-линейный У-образный закон изменения передаточного отношения механизма существенно упрощает начальный процесс синтеза системы, позволяя установить допустимые пределы изменения передаточного отношения для конкретной постановки задачи. Однако резкое изменение угла наклона прямой в середине хода вызывает такое же резкое изменение в ускорении движения объекта, который при опускании груза испытывает воздействие, близкое к ударному. В реальных механизмах передаточное отношение изменяется более плавно, что способствует плавности движения системы. Это можно показать сравнением процессов движения систем: 1) с механизмом, имеющим У-образный закон изменения передаточного отношения (рис. 7, 1); 2) с механизмом, имеющим параболический закон изменения этой функции (рис. 7, 2). Из сравнения следует, что при одинаковых значениях 1экв и одинаковых пределах изменения 1{у) в том и другом случае с переходом к параболическому закону характер изменения скорости становится существенно более спокойным. При этом время процесса движения как вверх, так и вниз остается приблизительно одинаковым в обоих случаях. Поэтому использование У-образного закона в расчетах можно считать целесообразным с точки зрения упрощения процесса синтеза и обеспечения некоторого «запаса прочности» проектируемой системы.

Рис. 7.

Процессы подъема и опускания объекта, полученные моделированием с использованием упрощенной системы уравнений при 1{у) кусочно-линейного (а) и параболического (б) типа при и=45. Подъем (\'ц =0): 133

безразмерная скорость объекта; 2- безразмерная скорость поршня. Опускание =5): 3-безразмерная скорость объекта; 4- безразмерная скорость поршня.

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

Моделирование приводной системы с пневмодвигателем.

В пятом разделе показана полная процедура моделирования и последующего синтеза приводной системы «механизм-привод-управление» с использованием в качестве движущей силы пневмопривод. Решается задача использования информационной технологии (модель наблюдателя) в целях сведения к минимуму количества датчиков в динамической управляемой системе.

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

динамику привода в целом. Динамическая модель приводного устройства имеет вид:

<? = S} - 62 - vj - vyrsign(i) + 5r

Г = (}^е)<Гс-У) prptiy)

где т — t/t*]^ — х!х*\8 = р!р* 5)5 S2 изменения давлений в рабочих полостях пневмоцилиндра (имеющие достаточно сложную функциональную зависимость, учитывающую, в том числе, сжимаемость среды);

ус запаздывание в отработке сигнала распределительным устройством в виде изменения положения у органа распределения (золотника);

- законы открытия рабочих каналов.

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

Из всех видов такой технологии наиболее мощной и универсальной считается т.н. модель-наблюдатель (observer) математическая модель объекта управления, связанная обратными связями с физическим объектом на входе и выходе, функционирующая с ним совместно. Основным требованием к системе объект модель-наблюдатель (далее «наблюдатель»)

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

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

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

Наблюдатель и вся система управления могут быть реализованы в цифровом виде, поскольку к быстродействию цифровых блоков не предъявляется особо жестких требований. Наблюдатель может быть выполнен как универсальный цифровой модуль. Обобщенные (безразмерные) параметры универсального модуля выбираются по результатам исследования безразмерной модели «объект-наблюдатель» вместе с безразмерными параметрами собственно привода.

Основные результаты и выводы.

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

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

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

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

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

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

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

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

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

Основные результаты работы отражены в следующих публикациях:

1. Лунев В.В. Мисюрин С.Ю. Особые многообразия плоских и пространственных механизмов с несколькими степенями свободы. // Проблемы машиностроения и надежности машин. 1993. № 1, с. 102-109.

2. Лунев В.В. Мисюрин С.Ю. Применение метода многогранников Ньютона в задаче о положениях механизмов. // Проблемы машиностроения и надежности машин. 1994. № 1.

3. Лунев В.В. Мисюрин С.Ю. Решение задачи о положениях механизма методом многогранников Ньютона // Проблемы машиностроения и надежности машин. 1994. № 2.

4. Лунев В.В., Мисюрин С.Ю. К вопросу о функции положения плоских и пространственных рычажных механизмов с несколькими степенями свободы. // Проблемы машиностроения и надежности машин. 1994. № 6.

5. Bessonov А.Р., Lunev V V., Misurin S.Y. «THE ANALYSIS OF PECULIARITIES AND THE CONSTRUCTION OF EXACT SOLUTIONS IN POSITION PROBLEM FOR MECHANISM WITH A FEW DEGREES OF FREEDOM. // International Conference» SPATIAL MECHANISMS AND HIGH CLASS MECHANISMS" (Theory and Practice) October 4-6, 1994, Almaty.

6. Лунев В.В. Мисюрин С.Ю. Положение равновесия механизмов и способы их определения. // Проблемы машиностроения и надежности машин. 1996. №6, с.12-17.

7. Крейнин Г.В. Мисюрин С.Ю. Яшина М.А. Выбор параметров позиционной системы гидропривод-передаточный механизм. //Проблемы машиностроения и надежности машин. 2001. № 6, с. 8 -15.

8. Крейнин Г.В. Мисюрин С.Ю., Яшина М.А. К синтезу позиционной системы с гидроприводом и механизмом передачи движения с переменным передаточным отношением. // Проблемы машиностроения и надежности машин. 2002. № 5, с. 17-23.

9. Лунев В.В. Мисюрин С.Ю. Определение «зон заклинивания» в механизмах с замкнутыми кинематическими схемами. // Проблемы машиностроения и надежности машин. 2003. № 5, с. 95 -99.

10. Крейнин Г.В. Мисюрин С.Ю. Яшина М.А. К синтезу позиционной системы с гидроприводом, механизмом передачи движения с переменным передаточным отношением и комбинированной нагрузкой. //Проблемы машиностроения и надежности машин. 2004. № 5, с. 3-10.

11. Крейнин Г.В. Мисюрин С.Ю. Яшина М.А. Сквозной синтез технической системы с приводом. // Труды конференции: Новые технологии

инновационному бизнесу (первый российский форум). Москва, 09-20 февраля 2007, с. 60-62.

12. Мисюрин С.Ю. Многопараметрический синтез системы «механизм-привод» // Труды конференции: Новые технологии инновационному бизнесу (первый российский форум). Москва, 09-20 февраля 2007, с. 62-63.

13. Мисюрин С.Ю. Синтез многозвенного механизма с особым положением // Труды конференции: Проблемы динамики и прочности

исполнительных механизмов и машин. Астрахань 10-16 сентября 2007 г., 34-35.

14. Мисюрин С.Ю. Синтез сложной механической системы «механизм-привод» // Труды конференции: VIII Международная научно-техническая конференция по динамике технологических систем. Ростов-на-Дону, 9-13 октября 2007. Том II, с. 51.

15. Мисюрин С.Ю. Ивлев В.И., Бозров В.М. Тепловые приводы для клапанно-запорной арматуры // Проблемы машиностроения и автоматизации (engineering & automation problems). ISSN 0234-6206 2007. № 2, с. 91-96.

16. Крейнин Г.В. Кривц И.Л., Мисюрин С.Ю. Яшина М.А. Пневматический позиционный привод: оценка возможностей и перспектив применения. // Проблемы машиностроения и надежности машин. 2007. № 1, с. 27-35.

17. Крейнин Г.В., Кривц И.Л., Мисюрин С.Ю. Яшина М.А. Модель наблюдатель в системе управления пневматическим позиционным приводом. // Проблемы машиностроения и надежности машин. 2007. № 4, с. 81-90.

18. Мисюрин С.Ю., Котенко Д.С. Методы определения особых положений плоских пространственных механизмов // XIX Международная Интернет-ориентированная конференция молодых ученых и студентов по современным проблемам машиностроения. Москва, ИМАШ РАН, 5-7 декабря 2007 г.

19. Мисюрин С.Ю. Ивлев В.И. Односторонний газовый привод с термосорбционным источником питания // Проблемы машиностроения и автоматизации (engineering & automation problems). ISSN 0234-6206 2007. № 4, с. 84-87.

20. Мисюрин С.Ю. Сквозной синтез позиционной механической системы как составной части многих автоматизированных систем // Общероссийский симпозиум «Достижения современной механики, энергетики, машиностроения» «Информационные диагностические технологии - основа безопасности страны» 22 октября 2007 г., с. 51-52.

21. Мисюрин С.Ю. Котенко Д.С. Кинетостатика механизмов с особыми положениями // Общероссийский симпозиум «Достижения современной механики, энергетики, машиностроения» «Информационные диагностические технологии - основа безопасности страны» 22 октября 2007 г., с. 62-63.

22. Мисюрин С.Ю. Крейнин Г.В. Выбор схемы включения привода в структуру механизма при решении задач кинематического синтеза. // Проблемы машиностроения и надежности машин. 2008. № 1, с. 3-8.

23. Мисюрин С.Ю. К вопросу о многопараметрическом синтезе сложных механических систем // VI Международная научно-техническая конференция «Материалы и технологии XXI века». Сборник статей. Пенза, 2008 г., с. 131 -132.

24. Мисюрин С.Ю., Крейнин Г.В. Позиционный привод, базовая система управления в безразмерных переменных // Международная научно-техническая конференция «Наука и образование - 2008». Труды. Мурманск, 02-10 апреля 2008 г., с. 679-681.

25. Мисюрин С.Ю., Крейнин Г.В. О некоторых общих закономерностях в динамике приводов // Проблемы машиностроения и надежности машин. 2008 №6, с. 15-21.

26. Мисюрин С.Ю. Ивлев В.И. Косарев А.А., Костин А.В. Определение границ мертвых положений в механизмах с одной и несколькими степенями свободы. // Проблемы машиностроения и автоматизации (engineering & automation problems). ISSN 0234-6206 2008. № 3, c.50-54.

27 Мисюрин С.Ю. Оптимизация механической системы // Приоритетные направления развития науки и технологий. Доклады всероссийской научно-технической конференции. Тула, 2008 г., с.190-192.

28. Мисюрин С.Ю. Синтез рычажных механизмов // 8 краткий отчет об основных результатах научно-исследовательских работ за 2007-2008 г. ИМА1Н РАН, с. 64.

29. Мисюрин С.Ю. Крейнин Г.В. Позиционный привод с идеальным двигателем // 8 краткий отчет об основных результатах научно-исследовательских работ за 2007-2008 гг. ИМАШ РАН, с. 70.

30. Мисюрин С.Ю., Крейнин Г.В. Пневматический позиционный привод и цифровая система управления // Труды конференции, посвященной 70-летию ИМАШ РАН: Проблемы машиноведения, с. 287-291.

31. Мисюрин С.Ю. Синтез механической системы «Исполнительное устройство механизм управление». // Проблемы машиностроения и автоматизации (engineering & automation problems). ISSN 0234-6206 2008. № 4, с. 24-26.

32. Мисюрин С.Ю. Метод выбора параметров механического устройства с выделением критериев динамического подобия системы управления // Машиностроение и инженерное образование. 2008. №4, с. 39-43.

33. Мисюрин С.Ю. Крейнин Г.В. Динамика и синтез параметров позиционного привода // Проблемы машиностроения и надежности машин. 2009. №6, с. 3-11.

34. Мисюрин С.Ю., Крейнин Г.В., Косарев А.А., Овчинников А.А., Саяпин В.В. Привод с высокооборотным струйным пневмодвигателем. Динамика и управление. // Проблемы машиностроения и надежности машин. 2010. № 6.

Мисюрин С.Ю.

Подписано в печать «_»_2009 г.

Формат 60x84x1/16 Усл.печ. л. 2. Тираж 100 экз. Отпечатано на ризограф

Подписано в печать: 22.11.2010

Заказ № 4634 Тираж -100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

10 - 27 7 3 6

/

1 :Л

í /

I U

2010181611

2010181611

Пр^идяум БАК ЙЛичобрнэуки России

(решение ст« зЫЬ. ъ у/35'

решил выдать дишюм ДОКТОРА срч^икЬ' и<.ш%е ¿ид ташШ£л\>к

Начальник очде.п£ ЩьАии^'

71 11-1/214

Федеральное государственное Ьюджетное образовательное учреждение высшего профессионального образования «Национальный исследовательский ядерный университет «МИФИ»

На правах рукописи

Мисюрин Сергей Юрьевич

Математическое моделирование и выбор параметров механизмов в комплексе с приводными системами

05.13.18 - Математическое моделирование,

численные методы и комплексы программ

Диссертация на соискание ученой степени доктора физико-

математических наук.

Научный консультант - доктор физико-математических наук профессор

Кудряшов Николай Алексеевич

Москва-2010

ОГЛАВЛЕНИЕ

Введение ..................................................................... ^

1 Постановка задачи ...................................................... ^

2 Простейшие математические модели позиционного механизма

О 1

(модель первого и второго уровня) .................................

2.1 Передаточная функция механизма в качестве его кинематической

21

модели .....................................................................

Уравнения динамики идеального двигателя ........................... ^

? Ч 99 Законы функционирования системы управления .....................

2.4 Модель первого уровня с механизмом с постоянным

35

передаточным отношением ................................................

2.5 Модель привода первого уровня с механизмом с переменным

44

передаточным отношением ................................................

2-6 Математическая модель привода второго уровня сложности ...... ^7

1 58 Механизм........................................................................

о 1

-5-1 Кинематика механизма ......................................................

<5 /го

Особые положения и особые многообразия механизма ............

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

75

окрестности особого положения .......................................

3.4 Пример точного представления функции положения механизма в окрестности особого положения..........................................

3.5 Выбор схемы включения привода в структуру механизма при

решении задач кинематического синтеза.................................

4 Моделирование приводной системы с гидродвигателем......... ^^

Гидропривод................................................................... ^^

- Дроссельное регулирование ................................................ ^ ^

112

- Математические модели системы и ее подсистем ........................

4.2 Выбор параметров позиционной системы гидропривод -

119

передаточный механизм ...................................................

122

- Механизм с чисто инерционной нагрузкой.................................

127

- Механизм с комбинированной нагрузкой .................................

132

- Силовая нагрузка как движущая сила (опускание груза)..................

335

- Особенности синтеза параметров системы малой инерционности ......

4.3 Синтез позиционной системы с гидроприводом и механизмом

1 Л/

передачи движения с переменным передаточным отношением......

- Постановка задачи. Диапазоны варьирования параметров, ограничения

1 39

и оценки вариантов решении................................................

141

- Анализ результатов решений................................................

4.4 Синтез позиционной системы с гидроприводом, механизмом передачи движения с переменным передаточным отношением и

комбинированной нагрузкой.............................................. ^^

148

- Моделирование передаточного механизма.................................

152

- Механизм с эквивалентным передаточным отношением..................

152

- Определение движущей силы двигателя....................................

154

- Определение эффективной площади каналов распределителя............

5

5.1

Моделирование приводной системы с пневмодвигателем .

Пневматический позиционный привод: оценка возможности и перспективы применения ..............................................

167 167

- Структурная схема пневматического позиционного привода ............

173

- Законы управления .........................................................

177

- Коэффициенты динамического подобия .................................

- Постановка задачи синтеза. Безразмерные критерии качества............ ^^

5.2 Модель наблюдатель в системе управления пневматическим 188 позиционным приводом ...................................................

190

- Динамическая модель объекта .............................................

193

- Линеаризованная модель объекта управления ...........................

197

- Структура системы объект-наблюдатель .................................

- Исследование системы «объект-наблюдатель» ........................ ^^

211

Основные выводы и результаты .......................................

215

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

235

Приложение 1 ...............................................................

Приложение 2 ............................................................... ^^

Введение.

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

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

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

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

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

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

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

Большой вклад в решение проблем кинематики, динамики и синтеза плоских и пространственных механизмов был внесентакими учеными как П.Л.Чебышев, академиками И.И.Артоболевским, Н.Г.Бруевичем, а также профессорами В.А.Зиновьевым, Ф.М.Диментбергом, А.Ф.Крайневым, В.А.Глазуновым, П.А.Лебедевым, G.J.Bennet'oM, C.Gosselin, J.Angeles, K.J.Waldron, V.Parenti-Castelli, C.Innocenti и другими. Наряду с созданием графических и графоаналитических методов исследования и синтеза механизмов, они много внимания уделили развитию аналитических и, что не мало важно, численных методов исследования (в том числе на основе математического моделирования). Значение последних весьма велико, так как они дают возможность использовать современные вычислительные машины с программным управлением для решения практических задач. Так, при анализе структуры механизмов приходится иметь дело с проблемой определения функции положения [3, 15-18, 41-43, 46-49, 54, 55, 82, 83, 1211261, которая для большинства плоских и пространственных механизмов может быть задана только неявно в виде системы нелинейных алгебраических уравнений, выражающих условия связей [42, 43, 45-50, 121126]. В настоящее время не существует сколько-нибудь эффективных общих методов решения нелинейных систем алгебраических уравнений [27, 28]. Кроме того, по теореме Н.Абеля и Э.Галуа найти конечные формы для алгебраического уравнения степени п > 5 невозможно. А.Е. Кобринский в своей книге [60] при анализе возможностей получения замкнутого решения обратной задачи о положении манипуляционной системы подчеркивает, что

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

В области теории двигателей (электрических, гидравлических, пневматических или других) мы пользуемся результатами, полученными многими авторами. Это направление представляют работы Е.В. Герц, Г.В. Крейнина, Е.А. Цухановой (Институт машиноведения РАН), Д.Н. Попова, В.М. Прокофьева (МВТУ им. Баумана), Б.И. Петрова, Н.С. Гамынина, A.B. Рабиновича, В.В. Саяпина (МАИ), Чемоданова Б.К. (МЭИ), а также зарубежных ученых W. Backe, Andersen В, Araki К, W.Frank, Takemura F. I. и многих других. Перед нами стояла задача - выделить. общие свойства двигателей различных типов, построить их обобщенную математическую модель и сформировать критерии подобия, которыми можно было бы удобно пользоваться в решении задач синтеза приводной системы в целом.

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

наблюдателей (в целях сокращения количества датчиков в обратных связях) и других факторов.

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

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

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

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

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

• Формируется математическая модель (далее модель) привода первого уровня сложности, в которую входят: 1) упрощенная модель передаточного механизма (в виде функции передаточного отношения); 2) модель идеального двигателя и 3) одна или более моделей системы управления. Под идеальным двигателем имеется в виду двигатель, движущая сила которого не зависит от скорости и который отрабатывает команды системы управления без запаздывания. Как будет показано далее, модель первого уровня сложности, представленная в безразмерном виде, оказывается весьма удобной для исследования взаимодействия

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

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

• Модель привода третьего уровня сло�