Текст
                    Государственный комитет СССР по народному обраао:
Московское ордена Ленина, ордена Октябрьской Революции к
ордена Трудового Красного Знамени
высшее техническое учнлнце им. Н.Э. Баумана
А.Г.КУЗНЕЦОВ
Утверадены
редсоветом МВТУ
СИНТЕЗ СИСТЕШ АВТОМАТИЧЕСКОГО РЕГШРОВАШ
ПНЕВМАТИЧЕСКОЙ КАМЕРЫ
Методические указания к до тине ну ааданив
по курсу "Оснош автоматики"
Под редакцией В.И.Крутова
Москва
1988


Дшные методические указан] издаются в соответствии с * ——• ~ 0-6 30.09.87 Г., х учвбно-иетодв- Дянмца методические указания издаются в сое тчебнш пданоы. Лсомотрены х одобрены кафедро! ыетодическо! ко«о<ш1 Факультета Э 25.04,88 г. ческим ленивы 20.06.88 г. Рецензент к.т.н^доц. Н.Б.Дешнков @ московское внаем техническое училище имени Н.Э.Баумана Александр Гаврилович Кузи»до» Редактор В.В.Оотваева Корректор *.И. Малютина ащ ЪМ . »*•« 1,76 н.я.(1,б уи.н»п.л.>. Тира* 300 зкз. Беопат». . Пояшоа* в веча» 29 04.06. План I960 г., i 1.3 доп. Типографа» МВТУ, 107006, Мооква, Б-6, 2не Бауияскм, 5.
Цель задания - спроектировать систему автоматического ре¬ гулирования (САР) давления в проточной пневматической камере, провести анализ статических и динамических свойств моментов САР (пневматической камеры как объекта регулирования и автома¬ тического регулятора) и всей систеш автоматического регулиро¬ вания в целом. ВЫБОР ВАРИАНТА ЗАДАНИЯ Номер варианта задания состоит из четырех цифр, показываю¬ щих номер столбца в соответствующем разделе таблица данных: первая цифра - номер столбца в разделе I, вторая - номер столб- на в разделе П и т.д. Часть исходных данных является общей для всех вариантов: площадь магистрали подвода J = 0,002 1^; температура воздуха Т- 293 К; атмосферное давление /з = К? Па; зона допустимой нестабильности давления Ар-^~ 0,005-105 На; максимальное допу¬ стимое отклонение давления' в переходном процесое лй = с m&jcQ = 0,03 -1Сг Па. Остальные исходные данные сведены в таблицу и определяются номером варианта. ПОРЯДОК ВЫПОЛНЕНИЯ ЗАДАНИЯ 1. По формулам (I) и (2) построить расходные характеристи- ки Gt(p) при рм , Gt(p)w р1гЫп и G3(pJ при /г . Опреде¬ лить значение Pm in * 2. Подсчитать значение фактора устойчивости /у, . * 3. Построить переходный процесс пневматической камеры apft) и определить показатели качества переходного процеоса £р>0?ая; ъ tp , 4. Сравнить показатели качества работы пневматической ка¬ меры ( pmiH , &Pmox я t^) й техническими требованиями ( Pmm0 > йРтол^ * *р0 > ® обооновать необходимость примене¬ ния регулятора. 5. Определить р2т1п , Лг'гкис • 6. По формуле (6) найти уравнение граничной прямой до тре¬ бованиям статики £(s) и провести ее на плоскости S . 7. Подсчитать величину степени устойчивости ex'. 3
Значения исходим данных для выполнения дошанего задания » | Я варианта раз¬ дела Параметр I г 3. 4 5 Номинальное давление в камере * 1Сг Па 1,4 1,5 1,6 1.7 1,8 Минимальное допустимое давление в камере /W ’ 1о5па 1,35 1,45 1,55 1,65 1,75 I Номинальное давление в магжотралж подвода ^ . I05 Па 2 2,2 2.2 2,4 2,4 Минимальное давление в магистрали подвода P/mi/i* Юг На 1,4 1.5 1,6 1.7 i.a ' Объем камеры V, i^ 2 3 4 5 6 П Максимальное допу¬ стимое время переход¬ ного процесса tpc , о I 1.5 2 2,5 3 ш Соотношение плеч ры¬ чага I 1,5 2 2,5 3 Масса подвижных ча¬ стей регуляторе т, кг 3 4 5 6 7 и Фактор торможения 1* , кг/с 100 120 140 160 Х80 8. ПО формул» ват-в штатам» Г£,Г бованжям динамики - SfS) ж прожияж ее да плоокоотн о, S • 9 Выбрать пять точек в области допустимых значений конст¬ руктивен Шфаштрож регуляторе А влоокоотж 6.S, Для выбран¬ ят точек провеет» аиалжв переходах продаооо* САР о опредале- нжем показателе! жх качеотве ДРтах* ^р • ^*,очв* пареходиих 4
процессов для всех пяти точек провести на ЭВМ. Переходный про¬ цесс для одной из точек построить по диаграше Выдаеградского и оопоставить с рассчитанным на ЭВМ. С учетом технических требо¬ ваний 49>хлг * Ш(3Рать параметры б' и 3 регуляторе. При необходимости провести анализ дополнительных переходных процес¬ сов САР. 10. В одной системе координат построить логарифмические амплитудные и фазовые частотные характеристики пневматической каюры и регулятора. 11. Путем сложения ооответотвувдих логарифмических харак¬ теристик пневматической камеры и регулятора построить логариф¬ мические амплитудную и фазовую частотные характеристики разомк¬ нутой систеш и определить запасы устойчивости САР по модулю и фазе. АНАЛИЗ СТАТИЧЕСКИХ И ДИНАМИЧЕСКИХ СВОЙСТВ ПНЕВМАТИЧЕСКОЙ КАМЕИ Пневматическая аИИШ «ьк у*гуяя,учш& объект Для шумильного протекания различных энергетических про¬ цессов часто требуется, чтобы значение какого-либо параметра, например давления, поддерживалось постоянным в пределах допу¬ стимых по техническим требованиям отклонений. Если изменение внешних условий работы агрег&та приводит к его неустойчивой ра¬ боте или недопустимым изменениям параметров рабочего процеооа, то агрегат оснащается регулятором, автоматически птягчарташтртм значение регулируемого параметра в заданных пределах. Необходи¬ мость установки регулятора следует из анализа характеристик рас¬ сматриваемого агрегата (объекта регулирования ОР). Объектом исследования в задании является проточная пневма¬ тическая каюра (рио. I). Воздух закачивается в каюру ко мире о- сором К и страниц веется из нее в атмосферу. Давление в магист¬ рали подвода р, , атюоферюе - ра. Давление р в пневматичес¬ кой каюре должно поддерживаться постоянным. Т РИо. I 5
Pi ПК РМо. 2 В теории автодатического регу¬ лирования для описания условий ра¬ боты рассматриваемого элемента поль¬ зуются функционально! схемой. Функ¬ циональная схема пневматической ка¬ меры приведена на рно. 2. Давления в магистрали подвода pf и атмосфер¬ ное ра являются дня пневматической камеры входнымк< или внешними, воз¬ действиями (координатами), а давле¬ ние р в хамере - выходной координа¬ той. Реши работы энергетического агрегата зависят от характера внеш¬ них воздействий. Ври неизменных внешних воздействиях параметры аг¬ регата оотаютоя постоянными во вре¬ мени, режимы работы называются установившимися, или равновеоными, а соответствующие характеристики - статическими. Изменение внешних воздействий приводит к изменению параметров во вреданв, агрегат работает на ^установившихся режимах. Процессы изменения пара¬ метров во времени называются переходными процессами и являются динамическими характеристиками энергетического агрегата. В тео¬ рии автоматического регулирования переходный процесо рассдатри- вам как реакцию на ступенчатое (мгновенное) изменение внешнего воздействия* Перехода! процесс пневматической камеры будем рассматри¬ вать как зависимость давления/jot времени t при ступенчатом изменении -»»«««■ р, в подводящем трубопроводе (рис. 3). Наи¬ более iff)™ показателями качества переходных процессов явля¬ ются время переходного процесса (время регулирования) £р и шксидальаое отклонение ординаты процесса pP^jc • Время пере¬ ходного процесса to определяется в тот момент, когда кривая шхшеооа (апериодического ми колебательного) вошла в зону до- ^стимсй нестабильности 2Лрд и в дальнейшем не «ходи. « 096 Большой интерес для практики представляет такие гармони- Рио. 3 6
чаский закон изменения внешнего воздействия, с помощью которого удобно описывать работу рассштриваемого элемента в условиях колебательных возмущений. При гармоническом изменении входного воздействия о круговой частотой колебаний cJ на выходе элемента устанавливаются колебания выходной координаты той же частоты о определенным фазовым сдвигом относительно входных колебаний. При различной частоте и) элемент ю-раэнсму пропускает колеба¬ ния, проявляя свои фильтрующие свойства. Для описания частотных свойств элементов используют частотные характеристики, показы¬ вающие зависимость параметров колебаний (в том числе относитель¬ ной амплитуда А и фазового сдвига f ) выходной координаты от ча- оготн колебаний и), Техдачеоксэ задание на гг,кект;гооад*ж9 САР пневматической камеры Давление р в пневматической камере определяется давления¬ ми ра и pf (рис. I и 2). Из условий эксплуатации известно, что атмосферное давление ра изменяется незначительно, а давле¬ ние Pf может уменьшаться от номинальной величины рш до мини¬ мальной Pfp-iin . В соответствии с техническими требованиями, предъявляемыми к пневматической камере, давление в ней должно поддерживаться постоянным и при номинальном давлении воздуха в магистрали подвода pf -pr/Y должно быть р . При уменьшении давления pf до р,тр) установившееся значение давления в камере не должно опуститься ниже р В наиболее опасном случае цри ступенчатом уменьшении давления в магистрали подвода от рм доPtmin вРемя установления давления в камере, т.е. время пере¬ ходного процесса, не должно превышать величину tpo , а ыахои- мальное отклонение изменения давления от рм не должно быть бодь- 1110 ДРтс/Хо • Одновременное удовлетворение всех технических требований представляет собой сложную, иногда противоречивую, инженерную задачу. Чтобы установить, соответствует ли работа пневматичео- кой камеры предъявляемым требованиям, необходима провести ана¬ лиз ее статических и динамических характеристик. Анализ статических саейотв даеаьтической камегы Из характеристик установившихся режимов пневматической ка¬ меры наибольший интерео представляют зависимости расходов в подводящей и отводацей магистралях 0, и от давления р в 7
в да3^®01 оаяаоти ощ»*»*®*0* мв~ дтнцжми зава о имост ям*: Х'Гл' .ЖТ д /О . Чтрт г, WW itfw **1 X пр*Л~ $*А/>> „м у _ / . площади подводящей и отводящей иагдстрадвй; W // и Jz а„ ртгияЛаты* А = 0.528 - критическое от- ^ s 1.4 - показатель адиаоатн, J3 в * ’ „„маний для воздуха. Если Отношение давлений в приве- - вшю формулах становится меньше критического, то проавоо денных ___ паоеходит в сверхзвуковую область, характеризую- "1ЯШМ ®nms"“ .-„та к воздуху уравнения состояния идеального газа yrsrzz -»& ' лг/л , а из уравнения адиабаты р,Ц /-></ УД объем воздуха в подводящей магистрали / Р If *Т V'm(o * = —7-7T3LГ ‘ * I Л/ / ft* Р * После подстановки выражений длк удельных объемов в ^рмулн расходов последние цринишют вид; ■№7r7F7"‘ “ 2* (2) * ’*J(j/-j}PT ! в щмвдь дашдаи «даю «кидал, от*>» Л™^'ж “ РИММ» [*““»• «Ив» “ *r'0It*““ ”>П0*‘ отвода на установившихся режимах. р ш По формула CI) определяем QtH HP* A Р’^СЧая/ Обходе» на ОДНОМ графике построить 3 характеристики. ,(р)
“Р* Р =Рн > Gt(f>) -р„у, и Ог(р) /г • При построении расходных характеристик до зна¬ чений Jb надо пользоваться формулами (I) в (2), при 0\р>&$кр расхода Q, и Gs остаются постоянными (рис. 4). Установившиеся режимы в ка¬ мере характеризуются равенством расходов воздуха в подводящей и отводящей магистралях, т.е. со¬ ответствуют точкам пересечения построенных характеристик. Точка пересечения характеристик Of(pJ при р,П1^ и Gp(p) показывает, какое давлениед,. установит¬ ся в пневматической камере оря снижении давления после коипрео- сора. С помощью расходных характеристик можно также оценить устойчивость работы ОР* т.е. его способность после внесения возмущений возвращаться в исходное или переходить в новое рав¬ новесное состояние. Вывод об устойчивости ОР делают по знаку фактора устойчивости на рассматриваемом режиме с дОг_ dQt п ~ ~др ср~ ‘ Если Fn ?Qi то ОР работает уотойчиво, если Рп<0% то не¬ устойчиво. , Определение производных, входящих в выражение для фактора устойчивости, удобнее всего сделать графически. Для этого «ядо провести касательные к характеристикам 0,(р) и G^fpJ и найти размерные тангенсы углов наклона касательных. Диапазон измене¬ ния давлений в пневматической камере составляет величину Рм~Pmln * Анализ построенных расходных характеристик свиде¬ тельствует об устойчивой работе пневматической камеры во воем диапазоне изменена давления Д, . так как на любых воз¬ можных режимах работы камеры углы наклона касательных к харак¬ теристике Gg(р) всегда положительные,'а углы наклона каса¬ тельных к характеристика G/pJ всегда отрицательные, поэтому рп Кроме оценки устойчивости работы пневматической камеры величина фактора устойчивости попользуется при построении пере- 9
ходшос провесов. Поскольку расходные характеристики нелиней¬ ные то касательные, проведенные к характеристикам Ot(pJ и Qjp) , например, при p-f>„ различный на- июн и величины факторов устойчивости различны. Возникает воп¬ рос, какое же значение Fn использовать при построении переход¬ ных'продаосов? Дм повышения точноотп раочета переходного про- цеоса F вычислить как разность частных производных, опреде¬ ленных по характеристикам: G,(р) , построенной для р,т:п . и GJP) . построенной дм }г . в точках, соответствующих значе- дав.»™ До */(А ^ ш рю-4)’ Для определения значений чаотных производных надо провести в точ¬ ках А и 3 касательные к расходным характеристикам и найти размерные тангенсы углов наклона касательных (РИ°* 5> Л 0/1 лГ dCh ■ -Jpz~ Tp‘ ^P 4° Рио. 5 дшда ддавдвю пга?*ита«№8 При реальной работе ОР характер изменения внешних возму- вродих воздействий может быть различным. При расчетном исследо- ниши принято рассматривать типовые возмудавдю воздействия, чаше всего ступенчатое и гармоническое. Теоретический анализ динамических свойств шявыатичеокой камеры проводится о помощь» дн^ренциальвого уравнения: Ж * °г°л > о применимости к воздуху £»«£>- стояния идэального газа, считая камере процессы изотермическими. В соответствии о методом ли теории автоматического регулирования функции расходов G, 7 Oz цредставхявтоя в виде линеаризованного ряда Тейлор®. Дифференциальное уравнение пневматической камеры в безраз¬ мерном виде [I» о.' 60 т dj' (3) +J> s«oF§ 10
где р - ДО /р^ - относительное изменение давления в камере, '’q*AQi/G,m - относительное изменение расхода в подводя¬ щемтрубопровода, Т0-Тп/Fn ~ поотояндая времени пневматической камеры, К г —Яж коэффициент усиления пненматичеокой камеры 0/ F"’P« по расходу Gi • При построении переходных процессов удобнее использовать шаференциалъное уравнение пневматической камеры в размерных координатах, дли частотных характеристик - в относительных (безразмерных) координатах. Рассмотрим переходный про песо изменения давления р в пнев¬ матической каш ре при ступенчатом изменении давления в подводя¬ щем трубопровода от pfM до p,min% что соответствует уменьше¬ нию расхода да величину dG/ . Дифференциальное уравнение пневматической камеры в размер¬ ных координатах daP ~аГ + Рп-йр (4) vjsji Лр - изменение давления в камере от значения рн, V- объ¬ ем камеры, Т - абсолютная температура воздуха в камере, R = 287 Дк/кг.град - газовая постоянная воздуха; aGf - изме¬ нение расхода воздуха в магистрали^додвода. при уменьшении дав¬ ления pf посла компрессор* (pic.4). Коэффициент пневматической камеры Тп = V/FT _ Решение дифференциального уравнения (4) при ступенчатом ' изменении давления р1 , т.е. расхода Qf , и нулевых начальных условиях t*0, &р*О да0^ выражение для построения переходно¬ го процесса изменения давления в пневматической камере: _ . — ' -• — По графику переходного процеооа определяется время переходного процесса tp (рис. 3). Обоснование необходимости шхименьния регулдтот Необходимость установки на OF регулятора может быть вы¬ звана либо неустойчивостью OF, либо несоответствием показате¬ лей качества работы OF предъявляемым к нему техническим требо¬ ваниям. Положительная величина фактора устойчивости пневмати¬ ческой камеры Fn свидетельствует об устойчивой работе OF. II
Используя полученные при анализе расходных и динамических характеристик числовые значения показателей качества работы пневматической камеры, необходимо показать, что величина статического отклонения давления до pmin и показатели качест¬ ва переходного процесса (Af>maX и tp > » соответствуют за- даичугм техническим УСЛОВИЯМ с р,Г)1л0 , * t/Ъ > • Несоответствие между полученными и заданными показателями говорит о необходимости оборудования пневматической кашры ре¬ гулятором. ПРОЕКТИРОВАНИЕ АВШЛАТИЧЗСГОГО РЕГУЛЯТОРА И АНАЛИЗ СТАТИЧЕСКИХ И ДИНАМИЧЕСКИХ СВОЙСТВ САР типа регулятора и ЗРоТРОешИ уаддд значений его ковотэтвдиад ЮРКШаМ Пяевматичеокая каме¬ ра и регулятор в совокуп¬ ности составляет систе¬ му автоматического регу¬ лирования (САРК Посколь¬ ку ОР (пневматическая ка¬ мера) задав и свойства его известны, то обеспечение заданных технических тре¬ бований достигается на¬ правленным выбором струк¬ туры и параметров регуля- I- Р„ Рис. 6 тора Для рассматриваемого ОР выбираем шмЗранный регулятор прямого действия, управлявший заслонкой, установленной в шги- SSL отвода воздуха. Схема получившейся САР приведена на оиГб. Основные конструктивные параметры регулятора - жест¬ кость пружины ё, эффективная площадь мембраны S , максималь¬ ный *одР= Мто*' Рассориваемого <ОР по юнструк- тивным соображениям удобно задать максишльный ход заслонки j/7 , который находится по разности максимальной и мини- шльво^влощадай магистрали отвода fz ~ fzmin • Минимальная площадь магистрали отвода^^ определяется из условия равенства ш установившихся режимах расходов вшд- ыодящей и отводящей магистралях О/min ~ P^min • Ч™/7опрвда“ ляется по формуле (X) приpf *р1ГП1Л и Р■ J2mis> И
определяется из формулы (2) при G2rrJin * Gf/77in и р ~pn7i0o . Для простоты примем, что в месте установки заслонки магистраль отвода имеет квадратное сечение поэтому площадь отводящей маги¬ страли до установки заслонки fs -а2 , где с/ - сторона квадра¬ та. Отсюда сторона квадрата а* Лр . После установки заслонки площадь шгистралн отвода определяется положением заслонки oh ~ \/^2 ‘h , а расход воздуха в магистрали отвода Л . _ /(7 2к (к ijRT р '7г р ХЧ х (5) Тогда максимальный ход заслонки а (л Л min) ~ ~7рГ (Л ~ fi'minj • При известном передаточном отношении <Jf рычага, связываю¬ щего шток с заслонкой, определяется величина АЕ„70^ ‘Zihmox/c//t Оставшиеся два параметра регулятора ( £ и S ) надо подобрать та¬ ким образом, чтобы САР удовлетворяла предъявляемым техническим требованиям. При выборе £ и S удобно использовать графическое изображение на плоскости с координатами в вида иокомых парамет¬ ров 6, S. Доя правильного выбора неизвестных параметров регулятора необходимо определить ограничения, которые накладывают на них требования по статическим и динамическим свойствам САР. Основ¬ ное требование по статике САР заключается в обеспечении стати¬ ческой ошибки (степени неравномерности), не превышающей заданную величину. Максимально допустимая по техническим требованиям ста¬ тическая ошибка ~pp,i/lo обеспечивается при равенстве соответ¬ ствующих изменений поддерживающей и восстанавливающей сил: /Э • <S ~ Prnln0 ' $ = & AZmoJc з откуда о о . с - f si г лил с с • «> В координатах 6.S зависимость (6) имеет вид прямой линии (I на рис, 7), вое точки которой соответствуют заданной статической ошибке. Для уменьшения статической ошибки необходимо увеличи¬ вать значение Pmir>» что приводит к уменьшению числителя выра¬ жения (6), т.е. к уменьшению утла наклона прямой I в координа¬ тах $ • Поэтому выбор любых значений £ и S , соответствующих точкам, ледащим на прямой I или под ней, обеспечивает приемле- 13
мую для данной САР статическую ошибку регулирования. Определение ограничений, на¬ кладываемых на конструктивные пара¬ метры регуляторами) качеству дина¬ мических процессов может проводить¬ ся либо с использованием частотных характеристик [2, 3], либо непо¬ средственно по переходным процес¬ сам. В данном задании будем прово¬ дить расчет по переходным дроцеосам с последущей проверкой частотных показателей качества. Область на плоскости двух параметров, обеспечивающих в пер¬ вом приближении затухание переходного процесса за время tp, найти методом -разбиения по заданной степени устойчиво¬ сти [il Математическая модель рассматриваемой САР включает диф¬ ференциальные уравнения ОР . регулятора и уравнение связи. Дифференциальное уравнение пневматической камеры о учетом заслонки в магистрали отвода: где &h - изменение положения заолонхи, (7) QG> ah <2- dh - определяется дифференцированием выражения (5) по h цри р *рм . Дифференциальное уравнение мембранного регулятора в отно¬ сительных координатах при неизменной настройке [I, с. 87 3 имеет над где Р Г2 Р dt2 лг/Zo + т — * 'к di + 3 «р</>> (8) - относительное изменение положения итога мембраны регулятора, Z0 -h0 /и* - положение итога мвиЗрагацга уотавовашавмоя режиме (примем. h0 ='hrrKfX~Cb^^ I • 7j^ • /~g ” постоянная эрамвн. paw, Z - постоянная времени S'- Ау катарагаа.Л^ * > - коэффициент усиления регулятора, т- Р О- 2г> 14
(леса движущихся чаете! регулятора, if - коэффициент гидравли¬ ческого трения. Дифференциальное уравнение регулятора в размерных коорди¬ натах (при неизменной настройке) т^Ш~,‘>'дм~^6йг Дифференциальное уравнение САР в размерных координатах, полученное из оистеш дифференциальных уравнений ОР , регулято¬ ра и уравнения связи Др - и, А 2 , имеет вид и Р5АР л с/гйР . л с/ЛР л _ _ c/*aG< Аъ -atrtA2-dP А'~ж~ tAo^Ps/n ~ррр~ ' где * ^ » Ai = Аг - Fnrr) At * Fn& + Тпё, ^~n & * ^os ' ‘G , (9) (10) ■ n v ' ''OS Общий интеграл дифференциального уравнения CAP (9) пред* ставляет собой алгебраическую о ушу трех составляющих, чячуч» из которых соответствует одному из корней характеристического уравнения. Время переходного процесса определяется наиболее мед- ленно затухающей составляющей &р~Др0 € , соответствующей корню характераотичеокого уравнения с наименьшей но модулю дей¬ ствительной частью< , где Др0=р„ -P/r,inQ ~ начальное откло¬ нение в шраходном проце.оое. Ори £ - Q Ар£„0 * лро . . Запиши условие скончания переходного процесса за время Л, - alt а _ **> др^ =дрд = Ар0 • & Рйзделжв начальное отклонение Aptx0 наРA ptn • подучим соотношение =e~°itgt откуда / . Ж, 4% Сп ■=£-&. , называем называемую степенью устойчи находим величину ос' = - —г- юи«. % . Обычно метод 0-разбиения применяется дхя построения границы облаете! устойчивости CAP (od-O) [I, о. 281]. В в том случае «пользуется уравнен» (9). Дкя построения границ обла- 15
отей заданного по времени tp качества переходного процеоса не¬ обходимо преобразовать уравнение (9) в смещенное уравнение. В характеристическом уравнении, соответствующем уравнению (9) Ад г 5+ А2 X 2 i-A, I +А0 - О производится переход к новой координате f путем подстановки ex'. Характеристическое уравнение в новой системе коор¬ динат получает вид -о, где ^2 * ^2. ^ ^ 4 - ^ «?о<Ai+.3*(aAz , зп (II) (12) ,0 - А0 7* Af * ъ(*Аг + ЫЛА& . На границе устойчивости WK0) переходный 'процесс стано¬ вится колебательным незатухающим, а хотя бы один из корней ха¬ рактеристического уравнения чисто мнимый и равен с и) , где г - мнимая единица, й) - круговая частота незатухающею переходного процесса. Бели рассматривается смещенное (<*<<?) характеристи¬ ческое уравнение (II), го наличие чисто мнимого корня гаранти¬ рует не только устойчивость, но и заданную степень устойчиво¬ сти U . Таким образом, в соответствии с методом 0-разбиения в уравнении (II) делается подстановка jr = iu) , после которой с учетом зависимостей (10) и (12) смещенное характеристическое уравнение принимает вид комплексного числа Q~32a)s + QfS+хог с/,6 * Ы(/?-А5а)2 + 7h$) ‘С, Если комплексное число равно нулю, то равны нулю его действи¬ тельная и мнимая части, т.е. получаем систему двух уравнений \ф -%,*)** Q, 6 * </,J = О, \R-A3u)z т* Т„6 = О о двумя неизвестными & ъ s , ^ Q - о(гА2 + oC3/4j + ы. Fn г? Qfs£/r+°(7/?t # = 2о(Аг * 3°С% * Fn г?. Решение системы о помощью определителей имеет вид 16
где V 4sc 6 ’ я, Xq2 Uf Гп 6 -Q+32oJ2 -A+AscJ2 s* = -к ог *оги1 'utr„, Q* -Q+d2a)£ V~Q,(fi-Abuj*)+Tn(Q-&»<*>*)* :-(M-Z,0)+(0,Aj-b$№ Tn -R+A50)2 Тогда. t-~4 Tn Tn и )* £ л -TnQ+Qtg _ QfA6 -Tt&z ^j2 XQS i/f Tn Ka2 uf Tp 'os Приравняв полученные из последних уравнений щ лсения для (А2 , находим связь между неизвестными конструктивнши па¬ раметрам регулятора / -R Aj (Tp Q - Я,Я) Aj Хог i/j ^jgj TJOA-WОА-Ъв; ’ имеющую в координатах £ S вид прямой линии (2 на рис. 7). Штриховка обращена к области значений £ и S » обеспечивающих заданную степень устойчивости U , т.е, время переходного про¬ цесса tn . Зои, образованная линиями I и 2, на рис. 7 показывает об¬ ласть возможных значений £ и s i ЧРД которых САР в первом при¬ ближении удовлетворяет заданным требованиям по статической ошибке и времени переходного процесса. 17
Анализ переходных процессов САР. Эыяь кснстьуктаыил дацеметшв -«сгуллтоюа Область возможных значений конструктивных параметров регу¬ лятора на рис. 7 расположена между осью S , отрезком £?/) пря¬ мой I и прямой 2 (вправо от точки А ). Если прямая I является точной границей допустимой статической ошибки процесоа регули¬ рования, то прямая 2 дает лишь приближенное представление о возможных значениях £ и S, обеспечивающих заданное качество работы САР в динамике. Окончательный выбор значений параметров регулятора £ и 5 делается после анализа переходных процессов рассматриваемой САР при различных значениях £ и S , т.е. дал различных точек области возможных значений параметров регулято¬ ра. Переходные процессы САР рассматриваем при ступенчатом изме¬ нении расхода в подводящей магистрали на величину AGt (рис. 4). Анализ переходных процессов САР проводится с помощью ЭВМ для нескольких (пяти-шести) вариантов значений , соответствую¬ щих выбранным точкам. Инструкция по использованию программы ана¬ лиза переходных процессов приведена в приложении. При анализе переходных процессов определяются показатели их качества tp и ДРтаХ' с учетом конструктивных соображений выбираются значения £ л S, при которых переходный процеоо удовлетворяет заданным требованиям tp0 и является наилучшим да одному из показателей А, или й^э^цри условии, что другой показатель не выходит за ограничения или/J/^j^). Если показатели качества переходных процессов не удовлетворяют заданным ограничениям, то выбираются новые значения & и S и про¬ водится анализ дополнительных переходных процессов. Окончатель¬ но выбираются такие значения £ и S , при которых показатели ка¬ чества tp и Afmox переходного процесса САР не превышают задан¬ ных значений tpQ и йртол. • Для выбранных значений § и 3 построить переходные процеооы да данным распечатки ЭВМ и по диаграмме Вышнеградского. Срав¬ нить построенные переходные процессы. В результате направленного выбора неизвестных конструктив¬ ных параметров регулятора £ и S статические и динамические свойства САР пневштической камеры соответствуют предъявляемым техническим требованиям. 18
Анализ лагагдДшчеоквх частотных ИМЮВШШ звеньев и разомкнутой САР Для построения частотных характеристик используют дифферен¬ циальные уравнения ОР я регуляторе в относительных координатах /формулы (3) и (8)/. Изменения параметров колебаний давления в пневматической камере при гармонической изменении расхода воздуха во впускной магистрали отражаются частотными характеристиками: амплитудной (АЧХ) ^ / / An(j) , ГАе Квг^ k° и фазовой (ФЧХ) /7, u>‘V Fn Гн ‘ -Oxcty(u)To) , которые строятся при изменении и) от 0 до со , АЧХ регулятора строится по формуле Ар (и)} = ФЧХ - по формуле /(/-Г0гсЛг/+ &(<А) = 'arcto 7-rfuJ* -Я - cnctq oJTA V /- ТргоОг при 0<u)i А/Тр.,, при //fc< cJ< оо при изменении <А от 0 до 00 . Анализ частотных свойств пневматической камеры, регулятора и САР в целом проведем по логарифмическим частотным характери¬ стикам: амплитудной (ДАХ) и фазовой (ЛФХ). Единицами измерения логарифмических координат являются де¬ када (дек) - по оси логарифмов частот и децибел (дБ) - по оси логарифмов амплитуд £(рис. 8). Декадой называется интервал ча¬ стот, соответствующий изменению чаототы в 10 раз. На логарифми¬ ческой шкале частот декада изображается отрезком единичной дли¬ ны, так как СдЮои-фи)-/, Поэтому относительна величины фи) логар^мичеокая шкала является равномерной, а относительно ча¬ ототы и) - неравномерной (рио. 8), Децибел используется при вве¬ дении логарифмической шкалы до оси ординат и! ЛАХ (рио. 8). 19
Чтобы получить ЛАХ из АЧХ, необходимо взять функции Г LJ3 б и построит* ее график, ио- пользуя логарифмическую шкалу частот (рис. 8). При построении Ш>Х лога¬ рифмическая шкала применя¬ ется только по оси частот, а по оси ординат использу¬ ется натуральный масштаб (в градусах или радианах). Для практических расчетов £((.(') А(и)) 40 20 L Рис. 8 оказывается удобным при изображении ЛАХ и ЛФХ использовать одну и ту же ось частот, совместив точку [Я оси ординат Ш о точ¬ кой 0 дБ оси ординат ЛАХ (рис. 8). Для удобства последующего анализа частотных свойотв САР ло¬ гарифмические частотные характеристики пневматической камеры и регулятора надо построить в одной системе координат. Логарифмические амплитудные характеристики некоторых звень¬ ев мало отличаются от своих асимптотических прямых, поэтому вме¬ сто точных ЛАХ с достаточной степенью точности можно строить со¬ ответствующие асимптоты. Асимптоты ЛАХ пневматической камеры строятся по следующим формулам: - низкочастотная ип„(и>)= 2ССу Л; / при и> < и>а = - высокочастотная C„a(u')r?CCj£%W и) > и>с -//Г*. ИРИ U)0 - J/to низкочастотная и высокочастотная асимптоты сопрягаются Спн ('x>cj * /иу, поэтому и >0* //rQ называется сопрягающей частотой. При построении асимптот пневматической камеры удобно пользоваться таким приемом: на осиL откладывает¬ ся значение кс, и проводится горизонтальная прямая (низ¬ кочастотная асимптота) до значения &и'с =Ц(//т0)> после ко¬ торого делается излом и при Clj Cl' ?ши ^ прш&я идет под углом - 20 дБ/дек (высокочастотная асимптота). Ввд амплитудной частотной характеристики регулятора зави¬ сят от соотношения его конструктивных параметров. Если I?** Ьтб\ то регулятор является апериодическим звеном 2-го порядка. ЛАХ 20
такого эвена весьма точно может быть заменена тремя асишгготамж, для построения которых надо определить две сопрягающие частоты Левее первой сопрягающей частоты ( Cd<u), ) раоподожева горизонтальная низкочастотная асимптота cPrt *206} *р . На пер¬ вой сопрягающей частоте U.J, горизонтальная прямая делает первый излом,и в интервале частот О.), <и> < и)г среднечастотная асимп¬ тота имеет наклон 20 дБ/дек. На второй сопрягающей частоте с<Ог прямая делает второй излом,и при cJ >al, высокочастотная асимп¬ тота имеет наклон 40 дБ/дек. Если 4/п£, то регулятор является колебательным эве¬ ном, В атом случае ЛАХ лучше строить не по асимптотам, а непо¬ средственно по выражению 6р(и.) *2000Ар (и>) где Ар(и}) под¬ считывается ш формуле (14). Ддя САР в целом наиболее важным чаототным показателем ка¬ чества является зато устойчивости. Один из удобных алоообов оценки samca устойчивости основан на частотном критерии Найкви¬ ста [ I, о. 2&0 1 и заключается в определении го логарифмическим частотный характеристикам разомкнутой САР запасов устойчивости го амплитуде (модулю) и фазе. Разомкнутая САР пнев¬ матической камеры (рис. 9) отличается от замкнутой отсутствием главной обрат¬ но* связи (пунктирная ли¬ ния), т.е. представляет . собой последовательное соединение двух звеньев САР; ОР и регулятора. Одно из преиму¬ ществ логарифмических частотных характеристик заложено в свой¬ стве логарифма произведения: логарифм произведения неокодьхих величин равен сумме логарифмов этих величин. Это свойство в значительной степени упрощает построение результирующих частот¬ ных характеристик последовательно соединенных звеньев. Известно [4, о. 122], что при последовательном соединении звеньев их АЧХ перемножаются, а ФЧХ складываются: ОР(ПК) !_ П7 . j Рио. 9 21
А (и>) * Ап (и>) Ар (и>), гм - /п (и>; +/0(и>). Прологарифмировав выражение дйя А(иО . получим ^ (а//а (и>) - ^^\Ап (ojJ Ар( u)j\p6jAn (и>) г ГдАр (и>). Последнее выражение объясняет смысл введения логарифмической шкалы по оси ординат АЯХ. После построения логарифмических ха- 'ракгеристик пневматической камеры и регулятора производят их графическое сложение, в результате которого получается ЛИ и jffiX разомкнутой САР. В соответствии с критерием Найквиста опре¬ деляют запас устойчивости по амплитуде Нт как расстояние (ве¬ личию в дБ) в точке пересечения ЛФХ с горизонталью, проходящей через I8CP, от оои логарифмов частот до ЛАХ разомкнутой САР, и запао устойчивости по фазе у4 как расстояние (величина в гра¬ дусах или раднаюх) в точке пересечения оси логарифмов частот о ЛИ от JBX разомкнутой САР до горизонтали, проходящей через I8CP (рис. 8). ЛИТЕРАТУРА 1. Основы теории автоматического регулирования /В.И.Кру¬ тов, Ф.М.Данилов, П.К.Кузьмик и др.; Под ред. В.И.Крутова. - 2-е кзд., перераб. и доп. - М.: Машиностроение, 1984. - 368 с. 2. Теория автоматического регулирования: Кн. I: Математи¬ ческое описание, анализ устойчивости и качества систем автома¬ тического регулирования /Под общ. ред. Солодовникова В.В. - М.: Машиностроение, 1967. - 768 с. (Техническая кибернетика). 3. Теория автоматического регулирования: Кн. 2: Анализ и синтез линейных непрерывных и дискретных систем автоматического регулирования /Под общ. ред, Солодовникова В.В. - М.: Машино¬ строение, 1967. - 680 с. (Техническая кибернетика). 4. Основы автоматического регулирования и управления /Под ред. Пономарева В.М. и Литвинова А.П. - М.: Высшая школа, 1974. - 439 с. 5. Крутов В.И. Переходные процессы систем автоматического регулирования. - Ы.: Машиностроение, 1965. - 252 с. 22
Приложение I гаочет Декхолнвд дродвссов на 3iM Переходный процесс САР определяется как решение оистеш дифференциальных уравнений, ошоыващих звенья САР о учетом связи между ними: c/i г Fn3p - - л'6; т сРс\2 dt2 & ctAZ с/Г + баг = sap Решение данной системы уравнений проводится на ЖМ числен¬ ным методом Рунге-Кутта. Для этого приведенные уравнения запи¬ сываются в вике систем! дифференциальных уравнений первого по¬ рядка: dpP_ dT dt f ~ ''vs d £ аг - - Я ; <&L*- йС, Тл d с/дг __ _6 dt ’ dt- ’ лч c/i * rri Список параметров и их идентификаторов, ме, приведен в таблица. иошльзуешх в црограм- Параметр £ т„ кол d <уС, 2>’ /77 ё S t Идентификатор FK TV Р и С О м S SX т Про грамм, построена в диалоговом режиме - общение о сШ црожо- ходит путем ввода числовых значенжй и ответа на вопросы. Поолв входа в вычислительную систему и запуска программ на экран* терминала появляется сообщение: щщгЕ исходна дтш: fp, tv, р, а4 с$,м. После этого с клавиатуры дисплея вводятся перечисленные пара¬ метры. В программе ясгользуетоя 1C? Па в качестве в дятлу дав¬ ления, поэтому цри вводе в 341 Fn я Тп надо увеличить в ICr роз (например, волн Fn « 2-1СГ5 м*о, то о клавиатуры надо вводить &)• Числе иные значения иоходных дя”™* вводятся черва запятую. После сообщения на экране 23
ВВЕДИТЕ ШАГ СЧЕТА "Н" И ПЕЧАТИ ”МР"; - о клавиатуры вводятся 0.001,50. Далее ва вкране появляется сообщения: ВВОД 5 ЗНАЧЕНИЙ В: - вводятся 5 выбранных значений 6(черев запятую); ВВОД 5 ЗНАЧЕНИЙ SK: - вводятся 5 выбранных значений S (через запятую). * После втого программа осуществляет расчет ординат переход¬ ных процессов и определяет для каждого варианта время переход¬ ного цроцесоа Гр и максимальное отклонение AW. На а кран вы¬ даются значения ГР и РМ ял 5 введенных вариантов. Затем на экране появляется сообщение: БУДЕТЕ РАСПЕЧАТЫВАТЬ? f У/N J Выбирается вариант (или несколько вариантов) значений 6 и S , удовлетворяющих заданным требованиям tpc и бр/ПО:Хо.Но оконча¬ тельный выбор значений ё и S делается после построения переход¬ ного процесса. Для зтого на вопрос "Будете распечатывать?" надо нажать клавишу " У ". Далее появляются сообщения: СКОЛЬКО ВАРИАНТОВ? - вводится количество выбранных вариантов (от одного до пяти); КАКИЕ ВАРИАНТЫ? - вводится номер выбранного варианта; если их несколько, то номера вариантов вводятся через запятую. Если ни один из пяти введенных вариантов значений £ и S не удовлетворяет заданным требованиям tpQ и Др/похд. то на вопрос "Будете распечатывать?" надо нажать клавишу "/У". Далее появляется сообщение: БУДЯ ОДЕ ВАРИАНТЫ? С У/Vj Если значения C*S выбрана и расчет закончен, то налиж¬ ется клавиша "/У", что завершает работу программы. Если есть необходимость в дальнейшем анализе переходных процессов САР, то берутся 5 юных вариантов 6 я S . На вопрос "Будут еще вариан¬ ты?" нажимается клавиша " у " и программа возвращается в свое исходное состояние - на экране появляется первое сообщение: ВВОД 5 ЗНАЧЕНИЙ В: ж повторяется.весь процесс расчета о новыми пятью вариантами значений 6 я S . Итогом работы на ЗЗМ являются выбранные конструктивные па¬ раметре регулятора ёж S ж график соответствующего им переход¬ ного процесса САР, построенный по распечатке. 24
Приложение 2 I&ffisi. пзрвхддных процессов пр ^ияграмме Вышнвграяскотд Дм построения переходного продаоса о аомощьв диаграмм! Вшие градского необходимо рассмотреть дийферендаадьное тж» ние САР в нормированном виде Пооле определения число на в печений коэффициентов X х у точка с координата»* (X,f ) наносятоя на дваграшу Вшавградского. Оолоиенмв точи на wafftua* позволяет определить тип переход¬ ного щюпеооа. Если получввнаж точка попадает в область I, то переходный продаос апериодический и строится как сум* трех вксповент: Если же точка располагаетоя в области П, то переходный пропеоо колебательный: По известным ча&ювыи значениям козффвциентов X и ^ можно без дополнительных расчетов определить вое необходимы» & ет~ ры переходных цродаооов. Значения втих парадетров в зависимости от значений X и £ либо сведены в таблиц»С 53, либо ш пола диаграи* Выше градского в шда *р<гжгг щщ постояв тх зна¬ чениях выбранного параметра Cl, о. 260-272J. Шпальные отдвонз- j)* С,ег,т+ Сг . 25
ния составляющих переходного процесса определяются по формуле: 'Г* к где начальное отклонение л <-'с 6'-4Z Y7>OJC начальная скорость начальное ускорение V %■ К — > AG, (*)/< £ 1‘ 'Ас > Ар/ ?/( \г Gt» м А0 / Значения и <;*берутся непосредственно с диаграммы Вышкеградского или из таблиц. Построение экспоненциальных составляющих переходного процеоса проводится с гомощью безразмерного отрезка времени , в те¬ чение которого ордината составляющей уменьшается вдвое (рис, П.1). Числовые значения такие берутся о диаграмм или из таблиц. Построение колебательных составляющих проиллюстрировано на рио, П.2. Сначала строятся огибающие экопоненты. На оси ординат откладываются значения ± С> . По оси Z откладываются периода ко¬ лебаний Тт. Ординаты огибающих экспонент определяются в соот¬ ветствии со степенью затухания: ордината после первого перио¬ да 7^ равна J?Cf , ордината после второго периода L ?гС, и т.д. После построения огибающих экспонент строят колебательные составляющие - синус и косинус. Для определения значений орди¬ нат огибающих экспонент после каждой четверти периода использу- 26
от коэффициенты ^ {с/с/ - 1/4, 1/2, 3/4, I). Эхотрашлыше ординаты составляющих смвщеш от точек иаоания ооотааяящкх о огибающими зкспонентами нв интервал дТс влево. Экстремальны* ординаты колебательных составляющих определяются путем увеличения ординат касания и <£у . рав. А д5 "косинусоида 'инусоиЗа АЪ Рис. П.2 Переходный процвсо находят графически прелфо соотавляпщх. Переход к размерному времена осуществляет ал в ооответот трех с ооотноиеннем t * Г Подробное описание методик востровнм переходных продао- сов в оиетешх 3-го порядка дяш в [I. с. 260-2723.
Длгапмяии» Выбор варианта а ада raw 3 Порядок выполнены задания 3 Анализ статических и джнамичеоких овойств пневматической и“иР“_- л* • • м: • • 5 Проектировании автоматического, регулятора и анализ стати- чеоких и динамических свойств САР 12 Литература 22 Приложены 23 I