Расчет следящей системы

1. Исходные данные

2. Выбор типоразмера двигателя и расчет редуктора

2.1 Предварительная оценка передаточного числа редуктора

2.2 Расчет редуктора с цилиндрическими колёсами

2.3 Проверка пригодности двигателя с учетом редуктора

3. Расчет усилителя мощности

4. Передаточная функция разомкнутой нестабилизированной системы

4.1 Функциональная схема разомкнутой нестабилизированной системы

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

4.3 Оценка общего коэффициента передачи

5. Коррекция следящей системы с использованием локальных обратных связей

5.1 Исходная ЛАЧХ

5.2 Коррекция с помощью ЖОС

6. Определение параметров предварительного усилителя

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

8. Принципиальная схема следящей системы

1. Исходные данные двигатель редуктор линейный синусоидальный

_нм=2.71 рад/c — наибольшая скорость вращения исполнительного вала при равномерной заводке и амплитуда скорости при синусной заводке.

_нм = 13.2 рад/c² — амплитуда ускорения исполнительного вала при синусной заводке.

М_нс = 1.25 Н^.м — статический момент сухого трения на исполнительном валу.

J_н = 26.9*10(-3) кг^.м² — момент инерции нагрузки исполнительного вала.

_ск=1⁰ — суммарная ошибка системы при равномерной заводке со скоростью _нм.

_д =2⁰ — динамическая погрешность при синусной заводке при амплитудных значениях скорости _нм и ускорения _нм.

_ip = 0.97 — КПД i-ой ступени редуктора.

_з > /4 — запас по фазе разомкнутой системы.

T_р = 0.1с — максимально допустимое время регулирования.

% - максимальная величина перерегулирования < 15.

Таблица 1

_д — без учета статической погрешности.

КПД одной пары зацепления редуктора _ip = 0,97.

серия двигателей ЭМ-М. f = 400 Гц.

УМ — двухтактный транзисторный усилитель в классе В.

ИУ — бесконтактные сельсины БД-160А и БС-155А в трансформаторном режиме. f = 400 Гц, К_и = 40 В/рад.

Тахогенератор: ТГ-5А. f = 400 Гц. К_тг = 1,1710^-2 Вс/рад, J_тг = 10^-7 кгм².

2. Выбор типоразмера двигателя и расчет редуктора.

Ptr1:=2*Mns*Wnm Ptr1 := 6.7750 Вт

Ptr2:=(2*Jn*Enm+Mns)*Wnm Ptr2 := 5.312 033 600 Вт

Ptr:=2*(Jn*Enm+Mns)*Wnm Ptr := 8.699 533 600 Вт Тогда самым оптимальным будет асинхронный двигатель ЭМ-8М.

Его технические характеристики приведены в табл. 2

Таблица 2

2.1 Предварительная оценка передаточного числа редуктора Диапазон допустимых значений передаточных чисел редуктора определяется корнями квадратного уравнения i_р1 и i_р2:

(1)

где

тогда решая уравнение (1) относительно ip получаем диапазон допустимых значений передаточных чисел редуктора [91.39 542 527, 212.3 663 658].

Примем ip=170.

Проверка выбранного передаточного числа редуктора:

— условие по скорости ;

170*2.71<1.2* 418.8 017 484 условие выполнено

— условие по моменту

(Jdv+Jn/(ip²*Kpd))*Enm*ip+Mn/(ip*Kpd)< 1.2*Mn т. е.737 784 1656e-2<.235 529 1026e-1 условие выполнено

— условие по перегреву

; (2)

Mn :=.196 274 2522e-1; sqrt ((Jdv+Jn/(ip²*Kpd))^2 * Enm*ip²*0.5+Mns²/(ip²*Kpd²))=929 837 0221e-2-> условие выполнено

2.2 Расчет редуктора с цилиндрическими колесами Передаточное число каждой пары зацепления i_n-1,n определяется как отношение числа зубьев ведомой шестерни (i_n) к числу зубьев ведущей (i_n-1):

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

Рис. 1

Для последней ступени редуктора выбираем i_n,n-1=10

Для минимизации приведённого к валу двигателя момента инерции редуктора будем рассчитывать передаточные числа ступеней пар зацеплений с использованием соотношения:

С учетом того, что i_1,2 < i_3,4 <.< i_{n-1,n получаем передаточные числа:}

i1,2=1.107 265 672, i3,4=1.854 818 960, i5,6=2.327 663 134, i7,8=3.765 295 059, i9,10=10 и так мы получили пяти ступенчатый редуктор.

Рассчитываем модуль зуба m= (sqrt ((0.2 635*Mns)/(Ksi*Rn)))=> m=.2 869 560 071 из стандартной шкалы модулей (по ГОСТу) m=0.3 мм Тогда число зубьев на ведущих шестернях равно 26, т.к. D= 8.0 mm, а диаметр и число зубьев ведомых шестерен равно соответственно: D2 := 8, Z2 := 27; D4 := 15, Z4 := 50; D6 := 19, Z6 := 63; D8 := 30, Z8 := 100; D10 := 80, Z10 := 267.

Расчет моментов инерции колес редуктора тогда момент инерции колес равен:

J1 :=.928 906 1159e-8 кг/м² J6 :=.295 546 8114e-6 кг/м²

J2 :=.928 906 1159e-8 кг/м² J7 :=.928 906 1159e-8 кг/м²

J3 :=.928 906 1159e-8 кг/м² J8 :=.183 694 8130e-5 кг/м²

J4 :=.114 809 2581e-6 кг/м² J9:=.928 906 1159e-8 кг/м²

J5 :=.928 906 1159e-8 кг/м² J10 :=.928 906 1159e-4 кг/м²

Где b — толщина колеса равна 3 мм, а = 7.7 * 10³ кг/м³

Приведённый момент к валу двигателя редуктора:

Jред = 839 351 1764e-7 кг/м²

Kpdred = 8 587 340 257

2.3 Проверка пригодности двигателя с рассчитанным редуктором

— условие по скорости ;

ip2*Wnm<1.2*Wnm т. е. 452.5 488 281<502.5 620 981 условие выполнено.

— условие по моменту

((Jdv+Jn+Jред+Jtg)/(ip2²*Kpdred))*Enm*ip2+Mn/(ip2*Kpdred)<=1.2*Mn

т.е. 261 327 7749e-2<.235 529 1026e-1 условие выполнено.

sqrt ((Jdv+Jред+Jtg+Jn/(ip2²*Kpdred))^2*Enm*ip2²*0.5+Mns²/(ip2²*Kpdred²))<=Mn

т.е. 890 692 9114e-2<.196 274 2522e-1-> условие выполнено.

3. Расчет усилителя мощности Схема двухтактного транзисторного УМ, работающего в классе B, представлена на рис. 2.

УМ нагружен на обмотку управления двигателя. Нагрузка замещается последовательно включёнными активным R_у и реактивным X_у сопротивлениями. Входной сигнал на переменном токе поступает от предварительного усилителя, в качестве которого используется операционный усилитель К140УД7. Напряжение насыщения и допустимое минимальное сопротивление нагрузки операционного усилителя составляют U_s=11,5B и R_{к min}=2kOм.

Рис. 2

Комплексный коэффициент усиления УМ на несущей частоте ₀=2f₀ (f₀=400 Гц):

Upumax = 5 В, тогда Kum := 16, Uunom:=80 В.

KpdTr1:=0.9:KpdTr2:=0.9 КПД трансформаторов Тр1 и Тр2 соответственно.

КПД коллекторной цепи Kpdk:=1/(0.1+(1/Kpd2))=8 256 880 735

Общая мощность нагрузки коллекторной цепи

Py = 63.99 497 659

Наибольшая мощность, выделяемая на переходе к-э транзистора Где Pk:=Py/Kpdk P := 24.2 655 843Вт Выбираем транзистор КТ855 В, его технические характеристики приведены в табл. 3.

Таблица 3

Напряжение питания U₀<0,5U_КЭ.ДОП

Примем напряжение питания равным U0 := 72.00 В Коэффициент трансформации Тр₂:

n2 = 1.903 077 509

Rэ:=0.1*Ry/n2², Rэ = 1.604 217 239 Ом. Примем Rэ равным 2Ом.

Амплитуда токов коллектора и базы:

Ikmax := 2.152 917 421 А, Ibmax =861 166 9684e-1 А.

Амплитуда входного напряжения транзистора:

Uvxtr := 4.305 834 842 В.

Коэффициент трансформации Тр₁:

n1 :=.9 568 521 871

Ivxmax:=Ibmax*n1 Ivxmax :=.824 009 4972e-1А.

Входное сопротивление усилителя

Rvxus = 60.67 891 228 Ом, т.к. входное сопротивление получилось меньше 2кОм, то необходимо на входе предварительного усилителя поставить эмитторный повторитель.

4. Передаточная функция разомкнутой нестабилизированной системы

4.1 Функциональная схема разомкнутой нестабилизированной системы Для системы с асинхронным двигателем эта схема представлена на рис. 3. В нее входят измерительное устройство ИУ, предварительный усилитель ПУ, усилитель мощности УМ, двигатель ДВ, редуктор Р и нагрузка Н.

Рис. 3.

В расчетах принимаем ИУ, ПУ, УМ и Р безинерционными элементами с соответствующими коэффициентами передачи: К_иу; К_пу; К_дм; К_р=1/i_p.

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

0 <_дв < (1,1…1,2) _ном =(1.1…1.2)* 418.8 017 484

аппроксимируется системой уравнений где b:=Mp/Uyn=.392 500 0000e-3 (н*м/в);

a:=(Mp-Mn)/Wn=.281 101 3761e-4 (н*м/(р/с));

Mном:=г*Mп=.196 274 2522e-1(н*м) щном:=г*щx=418.801 7484(р/с)

Tdv — постоянная времени двигателя Будем искать связь, для этого приравняем:

.

Перенесем все слагаемые, содержащие в левую часть уравнения, а остальные в правую часть, затем разделим уравнение на :

.

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

Tdv:=((Jdv+Jn+Jp+Jtg)/(ip²*Kpdred))/a=.399 658 5280e-1(с/р)

Kdv:=b/a= 13.9 629 3414(р/с)/в — коэффициент передачи двигателя по напряжению.

Km:=1/a= 35 574.35449(р/с)/(н*м) — коэффициент передачи двигателя по моменту.

Для уравнения применим преобразование Лапласа.

На рис. 4 показана структурная схема асинхронного двигателя с учётом редуктора и нагрузки:

Рис. 4. Структурная схема двигателя с учетом редуктора и нагрузки

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

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

где _.

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

В режиме равномерной заводки ошибка системы _сист складывается из трех составляющих:

· ошибки измерителя рассогласования _изм,

· статической ошибки _ст,

· кинетической ошибки _кин.

Люфтом редуктора в расчетах пренебрегаем.

Ошибка измерителя (у нас измерителем является измеритель на сельсинах) складывается из:

· технологической погрешности, которая рассчитывается по формуле:

_сд — ошибка сельсина-датчика,

_сп — ошибка сельсина-приемника;

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

· динамической погрешности (связанной со скоростями изменения угла) — не будем учитывать.

Тогда для ошибки системы можем записать:

Для нахождения ошибки системы изобразим структурную схему (рис. 5) замкнутой системы с двигателем, учетом редуктора и нагрузки (см. рис. 4). При этом введем обозначение .

Рис. 5

Рассмотрим статическую ошибку данной системы:

Будем считать все входы постоянными, то есть

— это режим включения.

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

Значит,

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

Будем считать все входы постоянными, то есть

— это режим, когда вал сельсин-датчика вращаем с постоянной скоростью.

Но сигнал по другому входу остается постоянным:

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

Тогда коэффициент усиления оценивается по следующему соотношению:

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

Таблица 3. Классы точности сельсинов

Вычислив по формуле (11) для 3 разных классов, подставим значения в соотношение (12) и получим три оценки К_р:

Сельсин 3-го класса не подходит. Тогда получаем два значения Kp которые равны соответственно:

Kp1 := 242.7 766 786 и Kp2 := 535.8 340 540

Динамическая ошибка _дин(t) системы определяется при отработке синусоидального воздействия с параметрами эквивалентной синусоиды

.

Рассчитаем параметры эквивалентной синусоиды:

Отсюда видно, что, значит, мы можем определить параметры эквивалентной синусоиды:

Тогда Щ= 4.87 084 8708(1/с) и A =.556 371 2121(р) Динамическая ошибка меняется во времени по гармоническому закону с амплитудой, Поэтому в расчете амплитудное значение этой ошибки связано с параметрами эквивалентной синусоиды в виде:

.

Это соотношение определяет положение модуля частотной характеристики на частоте :

=24.4 914 574

Полученные соотношения (13−14) задают на логарифмической амплитудно-частотной характеристике (ЛАЧХ) разомкнутой системы три контрольных точки (см. рис. 6). Причём на частоте = 1 [1/c] откладываются 2 точки, определяемые (13) и одна точка на частоте = .

ЛАЧХ разомкнутой нескорректированной системы строится согласно выражению:

где

Так как

При этом значение выбирается так, чтобы ЛАЧХ разомкнутой нескорректированной системы проходила выше контрольных точек.

Тогда из графика (Рис.6) видно что лучше всего принять Kp2 := 535.834 0540Kp:=535.8 340 540

Рис. 6.

После оценки требуемого значения определим частоту среза _ср путем решения уравнения:

Частота среза определяется в точке пересечения ЛАЧХ оси частот.

ср=114.446 1294(р/с)

Оценим запас по фазе на частоте среза разомкнутой системы по выражению: .21 524 3202(р), что меньше р/4,то необходимо решить вопрос о синтезе корректирующей обратной связи.

5. Коррекция следящих систем с использованием локальных обратных связей

5.1 Исходная ЛАЧХ

Исходная ЛАЧХ разомкнутой системы, проходящая выше контрольных точек, характеризуется требуемым минимально допустимым коэффициентом передачи и частотой среза ср, которая определяет полосу пропускания системы и её быстродействие. Время переходного процесса (ПП) системы оценивается по формуле

В техническом задании (ТЗ) указано время ПП —. Таким образом, видно, что время ПП системы, вычисленное с учетом выше полученного значения частоты среза удовлетворяет времени ПП, заданному в ТЗ.

Это означает, что — это время, за которое переходной процесс (ПП), должен входить в 5% зону установившегося значения. То есть, чтобы соотношение (16) выполнялось необходимо иметь:

.

5.2 Коррекция с помощью ЖОС Коррекция в системе с асинхронным двигателем осуществляется с помощью жесткой обратной связи (ЖОС), передаточная функция которой имеет вид. При коррекции вводится дополнительный предварительный усилитель К_ПУ1.

_K3:=._{588 235 2941e-2}

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

Kpu2:=Kp/(Ki*Kum*Kdv*(1/ip))=> Kpu2 := 10.19 348 220

Значит, K2:= 2277.29 4730(р/с)/в В скорректированной системе для обеспечения требуемого запаса по фазе на частоте среза реализуют

Введение

жёсткой обратной связи уменьшает постоянную времени исходной системы:

где К_ПУ2 = К_ПУ исходной системы.

Исходя из выражения, можно определить :

Kos:=(Tdv*Kp-1)/(Kpu2*Kum*Kdv*Ktg)=.766 206 2013(1/р) Найдем передаточную функцию охваченной части.

То есть, охватывая при помощи ЖОС, которое является инерционным звеном, мы уменьшаем его коэффициент усиления и постоянную времени.

Kpck:=K1*K2*K3/(1+K2*Kos*Ktg)= 25.213 6025(1/р)

Tdvck:=Tdv/(1+K2*Kos*Ktg)=.186 624 9435e-2(с/р) Из соотношения видно, что постоянная времени скорректированной системы уменьшилась в — раз по сравнению с исходной нескорректированной системы. Во столько же раз уменьшается и общий коэффициент усиления разомкнутой системы:

Отсюда находится

Kpu1:=1+Kpu2*Kum*Kdv*Ktg*Kos= 21.4 150 6493(1/р)

Lck (w) = 54.58 060 622−20 * ln (w)/ln (10)-20 * ln ((.348 288 6954e-5*w²+1)^(½))/ln (10)

Рис. 7. ЛАЧХ скорректированной системы Частота среза скорректированной системы.

421.246 6800(р/с) Запас по фазе на частоте среза разомкнутой скорректированной системы равен

Значит, условие, которое требует, чтобы запас по фазе разомкнутой системы был больше, выполнено.

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

На рис. 9 показаны ФЧХ исходной нескорректированной системы и ФЧХ системы (рис 9. а) после введения жёсткой обратной связи, где учтено значение .

Рис. 8.

Рис. 8а.

6. Определение параметров предварительного усилителя Сигнал с сельсина приёмника (СП) проходит через ограничитель, показанный на рис. 9. Он состоит из балластного резистора R_б и двух встречно включенных стабилитронов VD1 и VD2.

Рис. 9

Выбор R и типа стабилитрона должен обеспечивать ограничение U_пу по модулю до 10 В. По этой причине выбираем стабилитрон, напряжение стабилизации которого не превышает 10 В. Выберем стабилитрон Д814А, данные для которого приведены в таблице 4.

Таблица 4.

Оценим величину R_б, исходя из условия ограничения тока через стабилитрона при максимальном напряжении, снимаемом с выхода СП:

Выберем R_б из номинального ряда Rb := 750.0(Ом).

Здесь U_ст — напряжение стабилизации; I_ст — ток стабилизации; U_{ст мax} = 36 В.

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

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

Предварительный усилитель удобно строить на двух операционных усилителях типа К140УД7. На первом усилителе (D1) имеется два входа. На первый вход поступает сигнал с измерительного устройства. При этом реализуется коэффициент передачи По второму входу реализуется коэффициент передачи по тракту обратной связи:

Зададимся величиной R₀ = 10 кОм.

Тогда можем вычислить R₁ и R₂

R1:=R0/Kpu1= 466.9 609 937 выберем R1 из номинального ряда R1=560(Ом).

R2:=R0/Kos=13 051.31697 выберем R2 из номинального ряда R2=1500(Ом).

Рис. 10

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

Зададимся величиной R₃ = 10 кОм.

Предположим, что коэффициент передачи упругого звена, тогда R₄ = 10 кОм.

Из выражения

можем найти

выберем из номинального ряда: .

Из выражения

можем получить :

Выберем из номинального ряда: .

Постоянная времени С₁R₃ вводится для ограничения усиления высокочастотных помех усилителя.

Если входное сопротивление усилителя мощности менее 2 кОм, то на выходе предварительного усилителя ставится эмиттерный повторитель на транзисторе VT1.

7. Моделирование следящей системы и построение графиков переходного процесса и ошибок при линейной и синусоидальной заводке На рис. 12 представлена схемы моделирования следящей системы: схема для снятия переходного процесса (на вход системы подается скачок с единичной амплитудой):

для снятия ошибки при равномерной заводке (на вход системы подается линейнонаростающий сигнал):

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

0.556*sin (4.87*t)

График переходного процесса в системе изображен на рис. 13. Необходимо проверить перерегулирование, значение которого по ТЗ не должно превышать .

График ошибки при равномерной заводке изображен рис. 15

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

Из полученного графика видно, что ошибка не выходит за коридор

5,5*10^(-3)[рад]=

то есть требуемое условие выполнено.

График ошибки при синусоидальной заводке изображен рис. 14.

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

Из полученного графика видно, что ошибка не выходит за коридор

5*10^(-3)[рад]=

то есть требуемое условие выполнено.

8. Принципиальная схема следящей системы Принципиальная схема следящей системы изображена на рис. 11.

Рис. 11

Рисунок 12. Переходной процесс

Рисунок 13. Динамическая ошибка Рис. 14 Кинетическая + статическая ошибки.