Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Соломонцев Ю.М. Теория автоматического управления

.pdf
Скачиваний:
106
Добавлен:
02.05.2014
Размер:
5.48 Mб
Скачать

Рис. 2.9. Перенос сумматора:

а — исходная схема; б — перевес сумматора по

направлению передачи

информации;

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

 

 

против направления прохождения

информации

ввести элемент

с передаточным коэффициентом, равным В' = /С,.

 

Перенос сумматора. При переносе сумматора

по

направле-

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

через которое переносится сумматор (рис.

2.9, в).

2.2. ДИНАМИКА ЛИНЕЙНЫХ

СИСТЕМ

АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ

Сущность моделирования

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

 

Модель — некоторая другая

система,

сохраняющая

суще-

 

ственные черты оригинала и допускающая исследование физи-

 

ческими или

математическими

методами.

 

 

 

 

Моделирование — процесс проведения экспериментов

на мо-

 

дели

вместо

прямых экспериментов на самой системе. В настоя-

 

щее время широко применяют метод

моделирования как

способ

 

научного познания реальной действительности,

а в ряде случаев

 

он

оказывается единственным

средством

познания сложных

 

систем.

 

 

 

 

 

 

'

Чертеж детали, проект станка, система уравнений, описываю-

 

щих технологический процесс управления последним, и др. —

 

все это модели объекта проектирования, изготовления или управ-

 

ления.

 

 

 

 

 

 

 

Основой

моделирования является

теория

подобия,

которая

.

утверждает,

что абсолютное подобие

моделируемого объекта или

31

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

предсказать свойства проектируемого

объекта или про-

цесса.

 

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

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

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

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

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

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

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

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

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

Методы описания процессов в САУ

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

2 Теория автоматического

33

управления

 

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

нять

исследование.

 

 

 

 

Для получения математической модели в теории

автоматиче-

ского

управления используют один из двух путей:

 

1. Получение системы дифференциальных уравнений на основе

аналитического

анализа процессов

(физических, механических

и др.)

или экспериментальным

путем.

 

 

2.

Получение

косвенных

оценок

динамических

процессов,

к которым относятся

передаточные функции, временные характе-

ристики, частотные

характеристики.

 

 

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

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

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

чувствительного элемента—датчи-

У-, Уон,

Рис. 2.10. Модель механической обработки резанием на токарном станке

ка.

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

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

34

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

гиба резца с жесткостью /,. Жесткость детали несоизмеримо

боль-

ше, и поэтому

полагаем, что

(/г = со) — FBX

= /4хвх.

 

 

Но для любого момента времени,

вследствие динамики

явле-

ния, сила, вызывающая деформирование резца, оказывается

про-

порциональной

величине уг = хп — *ВЫх. являющейся

стрелой

прогиба

резца,

закрепленного

как

консоль.

Таким

образом,

для любого момента времени к резцу

приложена сила Ft

= /j#i =

/i (*вх — f/аых)- Силой, направленной в обратном направлении,

является

прежде

всего сила резания,

пропорциональная при про-

чих равных условиях производной отвеличины увых: Р = CdyBbIX/df, где С — коэффициент пропорциональности, значение которого

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

ленная наличием некоторой массы: / = тЛ/ВЫхА^8, где т — некоторая приведенная масса системы «резец—резцедержатель»;

d?y*va./dt2 — ускорение в

точке,

где

1/вых

— текущая

величина.

Уравнение сил

в общем

виде Рг

= Р

-f f

или /, (JCBX

1/вых) ==

= С dy^Jdt

+

nuPywJdt*.

 

 

 

 

 

 

 

При токарной обработке осевая сила (датчики контролируют

перемещения

в

осевом направлении)

 

 

 

 

 

 

 

 

P = CPt*p*Sl'p*vlp*KMo

-9,81.

 

 

 

 

 

*

 

 

 

рж

 

 

 

 

Если ур

= 1, а продольная подача

на оборот

детали

S =

= К dyvtaJdi, то при условии, что СРХ,

КмРх, ХРХ,

ПРХ

приняты

постоянными, определяемыми режимом обработки, материалами детали и резца и другими условиями, С = Ср t р*К.мр К, где К — некоторая постоянная величина. Из уравнения сил полу-

чим /2хвх = md2t/8blx/d^ + С dywJdt + /аг/Вых

и

после

деления

на / с введением обозначений т//а = 7\ и С//а

=

Т имеем

 

Полученное

дифференциальное уравнение

указывает

на то,

что рассматриваемая технологическая система

представляет собой

звено второго

порядка. Это теоретически

означает возможность

появления затухающего колебательного

переходного

процесса,

но на практике вследствие малости приведенной массы т и иногда возникающих незначительных ускорений величиной TWeux/^2

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

Уравнение движения поэтому может быть переписано так:

дгвх = Т dyvtn/dt +Увых1 это указывает на возможность представления рассматриваемой технологической системы как не-

2*

которого апериодического звена с передаточным коэффициентом К = 1 и постоянной времени Т = С//2.

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

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

чаемой детали и', следовательно, на точность станка.

 

В зависимости от связи колебательных систем узлов

станка,

механизма главного движения и механизма движения

подач

детали могут колебаться независимо друг от друга или

оказы-

вать влияние друг на друга.

 

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

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

На основании принципа Даламбера поведение системы описывает система обычных дифференциальных уравнений вида

mq (t) +nil (t) + jag (t) = P (t) + и (t);

где m — приведенная масса; п — коэффициент демпфирования; /п — жесткость подсистемы; Р (/) — изменение силы резания;

К — коэффициент пропорциональности.

Используя преобразование Лапласа при нулевых начальных условиях, это уравнение приводят к виду

где q (р),

Р

(р),

и (р), yt

(p)

преобразование Лапласа

перемен-

ных q (t),

Р

(t),

и

(t), yt

(t); p

— комплексная переменная.

Подставляя

в

полученную

систему выражения т//п

= Т*,

п/] = ДТ, где

£д

— коэффициент демпфирования, и

k =

1//п,

после преобразований находят

передаточную функцию q системы

W(p):

 

 

 

 

 

 

 

 

Я (Р)

=

 

 

 

 

 

 

 

Звенья системы станка представляют в виде передаточных функций. Строят частотные характеристики на основе передаточной функции с заменой параметра р на /со, где/ — мнимая единица; ш — частота. Для передаточной функции W (р)частотная характеристика имеет вид

После разделения на действительную и мнимую части передаточная функция примет вид

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

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

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

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

37

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

Для заданной функции

При у9 = 0 ряд Тейлора имеет вид

Рассмотрим систему, которая описывается дифференциальным уравнением второго порядка

F ( y , y , y , x ) = Q.

'

(2.2)

После линеаризации в окрестности заданной точки

(#„, уи,

Уо> *о) получим

 

 

floA£-h а1Д#-(-а1Д0 + &оДх = 0,

(2.3)

где &$ = у — у0; &у = у — у0; &у = у — у0;

Ьх = х —хв; а„ =

 

дР

 

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

зависимость, является гладкой.

Примем, что при постоянном входном воздействии х — х0 при / стремящемся к со выходная величина устанавливается и принимает постоянное значениех0. Тогда в установившемся режиме у — = 0, у — 0 и уравнение принимает вид F' (0, 0, $0, х0) = 0. Это уравнение называют уравнением статики в отличие от исходного, которое называют уравнением динамики. Коэффициентылинеаризованного уравнения являются постоянными, так как величины Уо,ние.ха не зависят от времени и время не входит в исходное уравне-

' Ограничения, накладываемые на уравнение (2.2)для получения

уравнения (2.3), следующие: отклонения Д»7, Д«/, Д#, Дх достаточно малы; функция F обладает непрерывнымичастными производными по всем аргументам в окрестности точки, соответствующей заданному режиму. При несоблюдении хотя бы одного условия линеаризацию проводить нельзя.

Геометрически линеаризация нелиней-

у

 

ной зависимости

между

двумя перемен-

 

 

ными

(рис. 2.11)

означает

замену

исход-

 

 

ной кривой АВ

отрезком касательной CD

*"

 

к А В

в точке

С,

соответствующей

задан-

 

 

ному

режиму

и

переносу

начала

коор-

 

 

динат

в

эту

точку.

 

 

о

 

Таким

образом, полученное алгебраи-

Yo

Y

ческое уравнение является аналогом урав-

нения (2.2), описывающего поведение объекта в малой области изменения ко-

ординат относительно установившегося

Рис 2.11. Графическая

интерпретация процеду- Ры линеаризации

режима

(например, равновесного,

характеризуемого

значе-

ниями

*0, у0). Полученное уравнение

динамики, связывающее

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

преимущество, которое дает линеаризация.

Описание процессов через передаточные функции. Дифференциальное уравнение звена САУ в общем виде запишем так:

.+ to, (2.4)

где у — выходная величина звена (в отклонениях от состояния равновесия); х — входная величина звена (в отклонениях от состояния равновесия); а„, ап_ь ..., аь а0, bm, bm-lt ..., Ьь Ь0 — постоянные коэффициенты, определяемые конструктивными особенностями и параметрами звена.

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

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

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

39

Если в уравнение (2.4), содержащее функции времени у (t) и х (t), ввести функции х (р) и у (р) комплексного переменного р, поставив условием, что эти функции связаны зависимостями

СО

00

 

У (Р) = [У (t) е-" dt, х (р) = J x (t) ег* dt,

(2.5)

то оказывается, что дифференциальное уравнение, содержащее

функции у (t) и х (0,

равносильно линейному алгебраическому

уравнению, содержащему функции у (р) и х (р):

 

а*-*Р"~1У (РН ----- 1- fliW (Р)+ о,» (Р) = *mPm* (Р) +

 

«х (р)+ • • • + MX (р)+ box (р).

(2.6)

Такой переход от дифференциального уравнения к однозначно соответствующему ему алгебраическому уравнению называют преобразованием Лапласа, интеграл (2.5) — интегралом Лапласа, комплексное переменное р — оператором. Сообразно с этим алгебраическое уравнение (2.6) является записью исходного дифференциального уравнения (2.4) в операторной форме.

Функцию у (р) называют изображением функции у (t), а функцию у (t) — оригиналом функции у (р). Операция перехода от исходной функцииу (t) к ее изображению у (р) (нахождение изображения по оригиналу) называют прямым преобразованием Лапласа. Математически прямое преобразование Лапласа записывают условно с помощью символа L [у (t) ] = у (р). Операцию перехода от изображения у (р) к искомой функции у (t) (нахождение оригинала по изображению) называют обратным преобразованием Лапласа. Математически обратное преобразование Лапласа записывают с помощью символа L"1 [у (р)] — х (t). Практически переход от дифференциального уравнения к алгебраическому происходит без каких-либо вычислений.

Если сравнить уравнение (2.4) с уравнением (2.6), то. нетрудно заметить, что формально переход от дифференциального уравнения к алгебраическому операторному уравнению при нулевых началь-

ных условиях 1 получают путем замены символов

дифференци-

рования оригиналов функций dn/dtn, dn~l/dfl~l, ...,

d/dt соответ-

ственно символами р", рп~1, ..., р и функций у (t) — их изображениями у (р). С оператором р можно, как и с другими членами алгебраического уравнения, производить различные действия (умножение, деление, вынесение за скобки и т. д.). Возможность записи дифференциального уравнения в операторной алгебраической форме значительно упрощает все расчеты.

Каждое звено САУ в общем случае описывается дифференциальным уравнением вида (2.4). Следовательно, при выводе диффе-

1 Нулевые начальные условия для дифференциального уравнения п-го порядка характеризуются тем, что для 1 = 0 значения самой функции у (t) и всех ее производных равны нулю.