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

GOS for Iphone / mobile / ТАУ / Книги / ТАУ учебник

.pdf
Скачиваний:
228
Добавлен:
10.05.2015
Размер:
4.33 Mб
Скачать

Производная i-го порядка от обобщенной функции (5.58) определяется выражением

i

 

 

 

 

 

 

 

 

 

(

)

 

 

D {F x

0

(t ± τ), x

0

F±τ (nT ±

τ)δ

i

 

(t − nT ) +

 

 

 

(t ± τ) }=

 

 

 

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

r

 

 

 

 

( )

 

 

 

 

 

+∑ ∑

 

 

 

(t − nT ),

 

F+ j (nT ± τ) − F− j

(nT ± τ) δ

i

 

(5.59)

 

 

 

j=1 n=σ j

 

 

 

 

 

 

 

 

 

 

где δ (i)(t–nT) – производная δ -функции Дирака порядка i. Подставим формулу (5.59) в выражение (5.57), в результате

±

 

Bqiτ =

 

F±τ (nT ± τ)δ(t − nT ) +

 

 

Di

 

0

 

n=0

r ∞

 

 

 

 

 

 

−ρqt

 

+∑ ∑

 

 

 

 

dt,

 

+F+ j (nT ± τ) − F− j (nT ± τ)

δ(t − nT ) e

 

j=1n=σ j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

.

(5.60)

 

 

i = 0,1,

, u

q =1,2

,, m

 

С учетом (5.59), (3.52) – (3.54) выражение (5.60) принимает следующий вид:

τ±

τ±

i

,

 

Bqi

= Bq

ρq

i = 0,1,…,u ; q = 1,2,…,m.

где

 

B

±

 

F

(nT ±

 

τ

=

 

q

 

 

 

±τ

 

 

 

 

n=0

 

 

 

r

 

 

 

 

 

 

+∑ ∑

 

 

 

 

 

λ

 

F+ j (nT ±

T ) F

j=1 nj

λT )e−ρqnT +

j (nT ±

λ

 

e

T )

−ρqnT .

(5.61)

(5.62)

При рассмотрении линейных несинхронных САУ (см. п. 5.2) коэффициент запаздывания λ , определяющий значение τ через длительность периода прерывания T, может быть целым числом, правильной и неправильной дробью. Причем последний случай является наиболее общим. Поэтому вычисление аналитического выражения, определяющего интеграл (5.61), рассмотрим для случая, когда λ является неправильной дробью.

181

Предположим, что выполняется следующее условие (h – 1) < λ < h, а значение λ = h – m, причем h – целое число, m – правильная дробь. Тогда соотношение (5.62), определяющее интеграл Bqτ± , будет

τ± =

Bq

 

r

+ j

 

 

 

(n ±

+

 

 

 

F

 

j=1 n=σ j

 

 

(n

± h)T

mT e−ρqnT +

 

 

F

 

 

 

±τ

 

 

 

 

 

n=0

 

 

 

 

 

 

 

h)T mT

 

(n ± h)T mT

 

 

F

e−ρqnT .

(5.63)

 

 

 

 

 

 

 

 

 

 

 

 

− j

 

 

 

Если в формуле (5.62) ввести обозначение k = n ± h, то она принимает вид

 

B

±

 

F

[kT mT ]e−ρqkT

+

 

 

 

 

τ

= e±ρqhT

 

 

 

 

 

 

q

 

 

 

±τ

 

 

 

 

 

 

 

 

 

 

 

 

kh

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

−ρqkT

 

 

 

+

 

 

 

 

mT ]Fj [kT

 

 

e

 

.

(5.64)

 

 

 

 

 

F+ j [kT

 

 

mT ]

 

 

j=1 kj ±h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если в соотношении (5.63) положить m = 0, то получаем выражение, определяющее интеграл, когда коэффициент запаздывания является целым числом. При h = 1 соотношение (5.62) сводится к случаю, когда коэффициент запаздывания представляет собой правильную дробь. По аналогии с предыдущими соотношениями (5.37), (5.40), (5.45) для сигнала с запаздывающим аргументом

Bq

 

[kT + mT ]e−ρqkT +

τ

= e−ρqhT

Fτ

 

 

 

 

 

 

 

k =0

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.65)

+∑ ∑

 

 

 

 

 

 

 

 

−ρqkT

 

 

 

 

F+ j [kT + mT ]Fj [kT + mT ]

e

 

,

 

 

 

 

 

 

j=1 k j h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а для сигнала с опережающим аргументом

 

 

 

 

 

 

 

 

 

Bq

+

 

 

 

 

 

+

 

 

 

 

 

 

τ

 

= eρqhT

Fτ [kT mT ]eρqkT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k=h

 

 

 

 

 

 

 

 

r

 

 

F

[kT mT ]F

[kT

mT ]

 

 

 

 

 

+

 

 

eρqkT .

(5.66)

 

 

 

 

+ j

 

 

j

 

 

 

 

 

 

 

 

j=1 k=σ j +h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

182

Далее, используя выкладки аналогичные вычислению интегралов Aqiτ± , Bqi, Bqi*, получены соотношения, определяющие интегралы Bqiτ± , для различных нелинейных характеристик при различных процессах на их входах.

Соотношения (5.65) используются для синтеза нелинейных непрерывных и импульсных САУ с запаздывающим аргументом методом ортогональных проекций. При этом наряду с аналитическими соотноше-

τ

τ

τ

τ

,

τ

ниями Aqi, Bqi, Cqi, Aqi

,

Bqi

, Cqi

, Aqi

Cqi при вычислении функ-

ционала требуется применять выражения, определяющие интегралы Bτ.

 

 

 

 

 

 

qi

Аналогично Aqiτ, Cqiτ, интегральные соотношения Bqiτсвязаны с Bqi

следующим образом:

 

 

 

 

 

 

 

 

Bτ

= eρqτ B .

 

 

 

 

qi

 

qi

 

 

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

183

Глава 6

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

С АЛГЕБРАИЧЕСКИМИ НЕЛИНЕЙНЫМИ ХАРАКТЕРИСТИКАМИ

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

6.1. Методы аппроксимации нелинейных характеристик

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

184

Рис. 6.1

тем управления, содержащих эле-

M

 

менты для которых характерны не-

 

 

линейные зависимости входных и

M 0

 

выходных величин, как правило,

 

стремятся по возможности упрос-

 

 

тить реальные характеристики.

0

 

Примером различных подходов к

ω = 0 ω

ω = ω хх

идеализации нелинейной характе-

 

 

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

единен с пружиной (т. е. вся система совершает колебания) [218], вид реальной характеристики которого показан на рис. 6.1.

В тех случаях, когда угловая скорость движения выходного вала невелика приближенно связь момента на валу M с угловой скоростью ω может быть установлена соотношением [218]

M = M 0 + αω + βω3.

(6.1)

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

M = M 0 + αω ,

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

M = M 0 + αω + βω3 + γω5 + υω7.

(6.2)

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

M = βω2.

(6.3)

Зависимость коэффициента сцепления µтормозящегося колеса транспортного средства с опорной поверхностью от величины относительного проскальзывания S имеет вид, показанный на рис. 6.2 (кривая 1).

185

 

 

 

 

э

 

 

 

б

3

1

 

2

 

0

S э

1

S

 

 

Рис. 6.2

 

 

Как видно из рис. 6.2, форма приведенной характеристики µ(S) аналогична рассмотренной выше зависимости момента на валу M асинхронного двигателя от угловой скорости ω (рис. 6.1). Следовательно, в основу аппроксимации характеристики µ = µ(S) могут быть положены те же принципы, что и для аппроксимации характеристики M(ω ).

При малых изменениях величины относительного проскальзывания колеса (в интервале от 0 до точки экстремума характеристики) справедливо соотношение аналогичное формуле (6.1)

 

 

 

µ = k

1 − kS 2

)

S ,

(6.4)

 

 

1

(

 

где S = 1

 

ωк

, здесь ωк (0) = ωкпос = ωспос – угловая скорость тормо-

 

 

 

ωс

− ωст

 

 

 

 

зящегося колеса в начале процесса торможения; ωс (0) = ωспос

– угловая

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

ωст = lст + rк dϕ ст ,

rк dt

где lст – длина стойки; rст – радиус колеса; ϕ ст – угол колебаний стойки. Зависимость, построенная в соответствии с формулой (6.4), показана на рис. 6.2 (кривая 2).

186

Аппроксимация вида (6.4) может использоваться в тех случаях, ког-

да необходимо изучить динамику выхода системы торможения в экстре-

мум характеристики µ(S).

 

При необходимости рассмотрения режима устойчивых автоколеба-

ний рабочей точки в окрестности экстремума характеристики µ(S), про-

исходящих с малой амплитудой, целесообразно использовать параболи-

ческую зависимость вида (6.3)

 

µ = µэ −k1 (S −Sэ ),

 

в соответствии с которой построена кривая 3 (рис. 6.2).

Данный режим соответствует торможению транспортного средства

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

са ωс = const и достаточно часто используется при исследованиях дина-

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

этапе их разработки.

 

При моделировании процесса торможения, когда ω c изменяется от

максимального значения ωс (0) = ωспос до минимального, величина уг-

ловой скорости свободно катящегося колеса ω

к также может меняться в

широких пределах. В каждый момент времени ω к может принимать

любое значение от ωк = ωс при полном растормаживании колеса, до

ωк = 0 – при возникновении состояния глубокого юза. Следовательно,

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

пределах от 0 (при полном растормаживаниии) до 1 (при юзе). В дан-

ном случае для аппроксимации характеристики µ(S) можно использо-

вать соотношение аналогичное формуле (6.2).

Еще одним примером различных аппроксимаций одной и той же ре-

альной нелинейной зависимости может служить сила вязкого трения,

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

F ( x )

жидкостях и газах можно считать пропор-

a 1

циональной скорости и направленной на-

 

встречу последней. Однако при больших

b 1

скоростях абсолютную величину силы вяз-

x

кого трения можно принять пропорцио-

b

нальной квадрату скорости.

 

Наконец, что касается силы сухого тре-

a

ния, то его зависимость от скорости дви-

 

жения v имеет вид, показанный на рис. 6.3.

Рис. 6.3

 

187

F

В тех случаях, когда спад характеристи-

a1

ки (участок a – b) не играет существенного

 

b1

значения при рассмотрении динамики кон-

кретной системы управления, то обычно

x

 

зависимость упрощают, полагая силу тре-

b

ния постоянной при v ≠ 0 (рис. 6.4). Такая

 

аппроксимация справедлива, когда в про-

a

цессе движения не имеется конечных ин-

Рис. 6.4

тервалов времени, в течение которых ско-

 

рость v = 0.

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

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

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

управления.

 

Достаточно часто качество аппроксимации оценивается по точности

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

ловно, малое отклонение кривой, построенной по формуле от аппрок-

симируемой характеристики, является достоинством. Однако, как от-

мечается в [218], нередко не соображения точности воспроизведения

реальной нелинейной характеристики являются определяющими при

выборе способа аппроксимации. Часто требуется невысокая точность

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

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

нейных черт явления. Поэтому основным критерием достоинства апп-

роксимации при решении той или иной задачи является достижение

достаточно точного и в то же время простого по форме и доступного

простому анализу результата. Удачная аппроксимация позволяет про-

интегрировать дифференциальное уравнение и получить уравнение

движения системы в конечном виде, т. е. осуществляется алгебраиза-

ция задачи, что достаточно часто ее облегчает.

А. А. Фельдбаум [218] делает по вопросу аппроксимации характери-

стик нелинейных элементов следующий основной вывод – универсаль-

ной аппроксимации не существует и вместе с тем определяющими кри-

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

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

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

относительно простой пригодный для исследования вид решения.

188

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

Рациональная функция

Примером такой аппроксимации служит математическое выражение следующего вида:

 

x + F (x) =

 

aF (x)

,

 

 

 

 

b

+ cF (x)

 

 

 

 

 

 

 

 

 

 

 

 

 

либо относительно выхода нелинейного элемента

 

 

F (x) = −

a + b cx

±

 

 

a + b cx

2

 

 

2c

 

 

 

 

 

2c

 

+ 4bcx ,

 

 

 

 

 

 

 

 

 

 

 

введя обозначение, d = − a + b ,

получаем

 

 

 

 

 

2c

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x2

 

d

 

b

 

2

 

F (x) = −d +

 

±

 

 

 

 

 

x + d

 

.

2

 

4

2

 

 

 

 

 

 

 

c

 

 

 

Аналитическое выражение характеристики скольжения S от момента асинхронного двигателя µд, выраженного в относительных единицах, является частным случаем данного вида аппроксимации

µд =

2

 

,

S +

1

 

S

 

 

 

 

где S – скольжение; µд – момент в относительных единицах.

Алгебраическая (степенная) функция

Примером может служить аппроксимация

F (x) = kx2 , –∞ < x < +∞ .

либо

F (x) = kx3, –∞ < x < +∞ ,

189

показанная на рис. 6.5.

которые обобщаются соотношением

F (x) = kxn , –∞ < x < +∞ ,

где n = 2, 3, 4, ..., – любое целое положительное число.

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

F (x) = k (x a )2 c , –∞ < x < +∞ ,

где k, a, c – любые как положительные, так и отрицательные числа.

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

 

 

F ( x )

плоскости.

 

 

 

 

 

 

 

Кроме того, к аппроксимации данно-

F ( x ) =

k x 2 ,

го вида относятся разнообразные мате-

при

x > 0 ;

матические соотношения, содержащие

 

 

 

 

F ( x ) = 0,

квадратичную и кубичную формы. На-

при

x

≤ 0

пример, нелинейная характеристика

 

 

 

x

 

2

, x > 0;

 

 

 

0

 

 

 

kx

 

 

 

 

 

 

F (x) =

x

0

 

 

 

 

0,

 

 

Рис. 6.5

 

 

 

 

Аналитическая аппроксимация

В некоторых случаях целесообразно представление характеристики нелинейного звена в виде конечной комбинации аналитических функций [218]

 

 

x

F (x) = a ln 1

+

 

, x >

 

 

 

b

F( x) = a arctg bx + cx, − ∞ ≤ ≤x+∞

;

 

x

 

F (x) = e

a+bx

1, x >

 

F(x) = a tgx3.

(6.5)

190