|
Регуляторы скорости и тока на ПИ-регуляторах, Прохождение сигнала |
|
|
|
Oct 30 2013, 08:09
|
Местный
  
Группа: Свой
Сообщений: 447
Регистрация: 16-11-08
Из: Украина, Донецк
Пользователь №: 41 684

|
Здравствуйте, уважаемые форумчане! Требуется Ваша помощь. Специфика вопроса немного узкая, но все же спрошу. В университете преподаватель задал вопрос. Пока не могу разобраться в нем. Есть система регулирования шахтной подъемной машиной, показанная на следующей схеме.
Система двухконтурная, построенная на двух ПИ-регуляторах. Суть вопроса. На выходе задатчика интенсиности (на входе регулятора скорости) напряжение прямолинейно возрастающее. Что будет на выходе регулятора скорости, а что будет на выходе регулятора тока? Я нарисовал обычную передаточную характеристику ПИ-регулятора и сказал, что наклон будет изменяться в зависимости от постоянной времени. Преподаватель сказал, что это неправильно.
Спросил у людей, которые работают с этими регуляторами вживую - мне нарисовали следующие схемы: Слева при линейно возрастающем, а справа - при постоянном напряжени.
Спросил - где можно найти примеры этих характеристик, чтобы разобраться, как они получаются - они не знают. Регулятор скорости настроен на симмитричный оптимум, а регулятор тока - на модульный оптимум. Кто-нибудь сталкивался с этими характеристиками и не знают ли литературу или хоть какой-нибудь источник, где можно было бы найти эти графики и объяснить их построение
|
|
|
|
|
 |
Ответов
(120 - 134)
Guest_TSerg_*
|
Nov 10 2013, 18:06
|
Guests

|
Цитата(AlexandrY @ Nov 10 2013, 21:55)  В курсовом уже все спроектировано и проверено. Не надо топтаться на месте. Вопрос простой - где формулы для временной области? И с ясным физическим смыслом коэффициентов в них. И что? Я их обещал?  Не надо тень на плетень наводить.
|
|
|
|
|
Nov 10 2013, 19:24
|
Местный
  
Группа: Свой
Сообщений: 447
Регистрация: 16-11-08
Из: Украина, Донецк
Пользователь №: 41 684

|
Получил вот такую ПФ:
Может такое получиться? Как определить величину, на которую нужно разделить числитель и знаменатель? У Вас получалось: Цитата Окончательно выглядит так:
PC = K1*(T1*p+1)/p РТ = K2*(T2*p+1)/p ТП = K3/(T3*p+1) Э = K4/(T4*p+1) M = K5/p Ce - ОС по эдс
Итоговая ПФ по выходу скорости и входу задания:
Числитель: K1*K5*K6*(T1*p+1)*(T2**+1)
Знаменатель: K8*p^5 + K9*p^4 + K10*p^3 + (K11 + K12*T1*T2)*p^2 + K12*(T1+T2)*p+K12
K6 = K2*K3 K7 = Ce*K4*K5 K8 = T3*T4 K9 = T3+T4 K10 = T3*K7 + ДТ*K6*T2 + 1 K11 = K7 + ДТ*K6 K12 = ДС*K1*K5*K6 И еще вопрос - как Вы приводите подобные? Если я раскрою скобки, то у меня в знаменателе нужно очень много места, чтобы увидеть всю строку. У Вас также? Дробь получается вот такая и я боюсь приводить подобные из-за того, что чувствую, что у меня где-то ошибка:
|
|
|
|
|
Nov 10 2013, 21:15
|

Ally
     
Группа: Модераторы
Сообщений: 6 232
Регистрация: 19-01-05
Пользователь №: 2 050

|
Цитата(TSerg @ Nov 10 2013, 20:06)  И что? Я их обещал?  А кто тут обещал всему научить?  Цитата(Andbiz @ Nov 10 2013, 21:24)  Дробь получается вот такая и я боюсь приводить подобные из-за того, что чувствую, что у меня где-то ошибка:
Возможно и ошибка. Но это не важно. Вот простой скрипт в Matlab который считает переходную характеристику вашей модели и показывает отклик на ступень: Код SC1 = 0.199; SC2 = 0.08; SC3 = 0.00008158;
CC1 = 0.0559; CC2 = 0.0664; CC3 = 0.036;
DRV1 = 76.986; DRV2 = 0.01;
ROT1 = 17.889; ROT2 = 0.0664;
MOT1 = 0.0559; MOT2 = 5.099;
CFB1 = 0.023; SFB1 = 0.306;
s = tf('s'); SC = (SC1*SC2*s+SC1)/(SC3*s); CC = (CC1*CC2*s+CC1)/(CC3*s); DRV = DRV1/(DRV2*s+1); ROT = ROT1/(ROT2*s+1); MOT = MOT1/(MOT2*s); CFB = tf(CFB1); SFB = tf(SFB1);
Sum1 = sumblk('ver = vin - vfb'); Sum2 = sumblk('cer = vout - cfb');
% Define block I/Os ("u" and "y" are shorthand for "InputName" and "OutputName") SC.u = 'ver'; SC.y = 'vout'; CC.u = 'cer'; CC.y = 'cout'; DRV.u = 'cout'; DRV.y = 'u'; ROT.u = 'u'; ROT.y = 'cur'; MOT.u = 'cur'; MOT.y = 'v'; SFB.u = 'v'; SFB.y = 'vfb'; CFB.u = 'cur'; CFB.y = 'cfb';
T = connect(SFB,CFB,Sum1,Sum2,SC,CC,DRV,ROT,MOT,'vin','cout'); step(T), grid Строился скрипт по этой модели:
А вот собственно и сама переходная функция для тока: Код >> T
T = From input "vin" to output "cout": 20.12 s^5 + 2870 s^4 + 9.789e04 s^3 + 1.271e06 s^2 + 5.704e06 s ----------------------------------------------------------------- s^5 + 115.1 s^4 + 6425 s^3 + 2.141e05 s^2 + 3.858e06 s + 2.635e07 Continuous-time transfer function. А вот переходная функция для скорости: Код >> T
T = From input "vin" to output "vout": 195.1 s^5 + 2.489e04 s^4 + 1.534e06 s^3 + 3.013e07 s^2 + 1.807e08 s ------------------------------------------------------------------- s^5 + 115.1 s^4 + 6425 s^3 + 2.141e05 s^2 + 3.858e06 s + 2.635e07 Continuous-time transfer function. Как видно все функции 5-го порядка. У кого получилось иначе значит допустил ошибку либо вообще неправильно считал.  Можете эти функции подставить в блок Matlab и они дадут точно такой же график как ваша исходная модель. Да , и присмотревшись к документу на который ссылался TSerg я таки увидел, что они решали численно. Т.е. аналитического решения или нет или оно не практично из-за необозримости результата. Ну и наконец вот преобразование отклика скорости на ступенчатое воздействие в функцию от времени:
MathCAD легко решает численно такое преобразование. Можно сравнить с графиком по модели в Matlab, будет идеальное совпадение. Численное инверсное преобразование Лапласа можно сделать и в Matlab-e, кстати, только формулы будут не презентабельные.
|
|
|
|
Guest_TSerg_*
|
Nov 11 2013, 06:27
|
Guests

|
Цитата(AlexandrY @ Nov 11 2013, 01:15)  А кто тут обещал всему научить?  1. Не надо мне приписывать то, чего я не говорил. 2. Не научить, а подсказать путь решения, а именно - операторным методом через теорему разложения, поскольку препод отказывается признавать решения через математическое ПО или намекает на простой путь. Готовые же решения я приводить не буду, принципиально - это задача студента. Путь, предлагаемый Таней я вовсе не отвергаю, о чем уже говорил, но с учетом некоторых оговорок. В условиях установившихся процессов ( разгон, стаблизация скорости, после наброса нагрузки) возможно получение простых выражений для любой точки схемы. В переходных процессах это возможно только путем решения полных диф.ур-ов.
|
|
|
|
|
Nov 11 2013, 08:05
|

Ally
     
Группа: Модераторы
Сообщений: 6 232
Регистрация: 19-01-05
Пользователь №: 2 050

|
Цитата(TSerg @ Nov 11 2013, 08:27)  В переходных процессах это возможно только путем решения полных диф.ур-ов. Ну кто в наше время решает диф.ур-ы? Привожу порядок действий в Matlab для того чтобы узнать отклик скорости вращения двигателя на линейно нарастающее напряжение на задатчике: (в предположении, что выше опубликованный скрипт был уже запущен) Код >> T = connect(SFB,CFB,Sum1,Sum2,SC,CC,DRV,ROT,MOT,'vin','v')
T = From input "vin" to output "v": 4.575e05 s^2 + 1.261e07 s + 8.613e07 ----------------------------------------------------------------- s^5 + 115.1 s^4 + 6425 s^3 + 2.141e05 s^2 + 3.858e06 s + 2.635e07 Continuous-time transfer function.
>> RAMP = 1/(s^2)
RAMP = 1 --- s^2 Continuous-time transfer function.
>> series(RAMP, T)
ans = From input to output "v": 4.575e05 s^2 + 1.261e07 s + 8.613e07 ----------------------------------------------------------------------- s^7 + 115.1 s^6 + 6425 s^5 + 2.141e05 s^4 + 3.858e06 s^3 + 2.635e07 s^2 Continuous-time transfer function.
>> impulse(ans) В результате имеем такой график:
Что вполне правдоподобно.
|
|
|
|
Guest_TSerg_*
|
Nov 11 2013, 10:58
|
Guests

|
Цитата(AlexandrY @ Nov 11 2013, 12:26)  Они и сейчас без matlab-а работают. Эт я тут просто развлекаюсь. В жизни никто моделей объектов в готовом виде не дает и в лучшем случае остается применять fuzzy logic.  ТС объяснил, что Матлаб и Маткад ему не в помощь и для чего тогда Вы тут развлекаетесь? В жизни, инженерной, модели разрабатываются, тестируются, верифицируются, а всякие фази-лоджик - это имеет место быть, но не принципиально заметно на общем фоне. Лучше бы Таня привела пару формул для установившихся процессов. Тоже помощь студенту. Я "хлеб" отнимать не буду  Цитата(Andbiz @ Nov 10 2013, 23:24)  И еще вопрос - как Вы приводите подобные? Если я раскрою скобки, то у меня в знаменателе нужно очень много места, чтобы увидеть всю строку. У Вас также? По мере свертывания структурной схемы (СС) из отдельных звеньев происходит последовательное преобразование в более простые выражения с учетом ввода новых объединяющих коэффициентов Ki. Из-за этого при каждой очередной трансформации СС, выражение не "расползается" в бесконечность. Можно было бы еще и их упрощать для красоты, но в этом нет смысла - так быстрее и меньше ошибок. Ну, разве, что чуть увеличиться погрешность численных вычислений. В конечном итоге все равно отдадим на расчет компьютеру для вычисления корней. P.S. Кстати хотел взять за образец расчета САУ вышеприведенный мной документ BDC-SAC.doc, принять за основу сделанные там расчеты и показать наэтом примере преобразования, но оказалось, что и там - полная лажа. Итоговые цифры не соответствуют исходным данным, в итоге все расчеты - подгонка. Достаточно привести вычисление номинального тока Iн = Pн/(Uн * КПД) = 56.18 А На самом деле расчет дает 66.39 А Ну и дальше. Хотя, методически, документ выглядит гораздо лучше курсового студента.
|
|
|
|
|
Nov 11 2013, 11:16
|
Гуру
     
Группа: Модераторы
Сообщений: 8 752
Регистрация: 6-01-06
Пользователь №: 12 883

|
Цитата(TSerg @ Nov 11 2013, 14:58)  Лучше бы Таня привела пару формул для установившихся процессов. Тоже помощь студенту. Я "хлеб" отнимать не буду  Уже же ж... И сам студент приводил формулу. Там только зачем-то есть член, учитывающий выпрыгивание пассажиров на ходу... Я обычно словами объясняю... Вот еще и про формулы. Арнольд часто приводил формулу для причаливания корабля (уравнение) с весьма полезными для данной дискуссии комментариями. Это я к тому, что невозможно, глядя на многоэтажные формулы, понять смысл происходящего. Без упрощения. Еще... вот даже расчет классической механической задачи (бильярд) очень быстро приводит к стохастизации. Вот Ваши любимые корни... А будем исследовать, как эти корни будут меняться при небольших отклонениях (точность параметров и пр... ), и к чему это приведет? Это я к вопросу о робастности...
|
|
|
|
|
Nov 11 2013, 13:09
|

Ally
     
Группа: Модераторы
Сообщений: 6 232
Регистрация: 19-01-05
Пользователь №: 2 050

|
Цитата(TSerg @ Nov 11 2013, 12:58)  ТС объяснил, что Матлаб и Маткад ему не в помощь и для чего тогда Вы тут развлекаетесь?
В жизни, инженерной, модели разрабатываются, тестируются, верифицируются, а всякие фази-лоджик - это имеет место быть, но не принципиально заметно на общем фоне. Ну вы еще поучите инженерной жизни  Наши инженерные системы работают по всему миру от Австралии до Перу. На интуиции (она же fuzzy logic) все и держится. Но судя по всему преподаватель имел в виду вывод вот такой трансцендентной функции:
Ну выведете ее без MatCAD-а на бумаге (или в Word-е), и без ошибок, я буду смеяться.
|
|
|
|
Guest_TSerg_*
|
Nov 11 2013, 17:51
|
Guests

|
Цитата(AlexandrY @ Nov 11 2013, 17:09)  Ну вы еще поучите инженерной жизни  Наши инженерные системы работают по всему миру от Австралии до Перу. Как же чел выпучивает свое офисное планктонство на основе Маткада и пр.  К примеру: "Термостат для проращивания яиц крокодилов с целью последующего пожирания в респектабельный ресторанах". Цитата На интуиции (она же fuzzy logic) все и держится. Бред для домохозяек. Цитата Ну выведете ее без MatCAD-а на бумаге (или в Word-е), и без ошибок, я буду смеяться.  Легко, если бы мне это было надо. То, что Вы сейчас вымучиваете с софтом, злясь на невозможность символического решения после 5-й степени, мы спокойно решали с логарифмической линейкой и свыше 5-го порядка. P.S. Да, я владею современным софтом и пользуюсь им, но здесь учебный случай, оговоренные условия и все ваши потуги напрячь студента на софт или дискретные системы, тем более использовать теорию Лютфи Заде (между прочим, моего земляка) - идут лесом. Извините, но от Вас в этой ветке только белый шум, не более.
|
|
|
|
|
Nov 11 2013, 18:23
|

Ally
     
Группа: Модераторы
Сообщений: 6 232
Регистрация: 19-01-05
Пользователь №: 2 050

|
Цитата(Tanya @ Nov 11 2013, 17:07)  Инженеры учитывают члены в таких задачах какого порядка и с какой точностью? Таня, вы серьезно хотите поговорить о членах? Я бы не фокусировался слишком на этом. Современные инструменты позволяют автоматически сократить количество членов без ущерба качеству. Так что это вопрос второстепенный для численных методов. Такая фича в частности есть в Control System Toolbox. Цитата(TSerg @ Nov 11 2013, 19:51)  Извините, но от Вас в этой ветке только белый шум, не более. Ок, полегчало? Но все таки не отвлекайтесь от существа дела (если есть что сказать), а то ветку снесут
|
|
|
|
|
  |
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|