Файл: Севостьянов, А. Г. Основы математического моделирования механико-технологических процессов текстильной промышленности из цикла лекций заочного факультета по технической кибернетике в текстильной промышленности.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 31.10.2024
Просмотров: 66
Скачиваний: 0
устанавливают исходя из анализа экспериментальных данных ра нее проведенных исследований или из теоретического анализа мо делируемого процесса.
Если г/max и хт ах заранее неизвестны, то во избежание перена-„ пряжений следует масштаб выбирать более осторожно (например, Мтах=бО В), ,а затем его скорректировать. Для выбора масштаба времени надо знать длительность t реального процесса и задаться машинным временем т, в течение которого будет проводиться ис следование модели на АВМ.
Для однократного процесса решения задачи тШах обычно со ставляет десятки-сотни секунд. Очевидно, что при Mt> 1 моделиру ющий процесс на АВМ будет протекать быстрее, чем реальный мо делируемый процесс. При Mt< \ моделирующий процесс протекает медленнее, а при Mt= 1 — в реальном масштабе времени. При вы боре значения Kt учитывают необходимую степень изменения ско рости протекания процесса, частотные свойства модели, удобство наблюдения и регистрации выходных параметров.
Из соотношения (56) находим:
и= — — ;
му
II
M t ’
(60)
из соотношений (57) имеем:
и. = y / C v ; V = хКх,
II *1
(61)
Используя соотношения (60) и (61), можно установить следую
щие масштабные соотношения: |
|
|
|
||
для первой производной — |
|
|
|
||
М у |
и' |
_ _ Му . |
К у = м' _ |
Ку |
(62) |
У |
Mt ’ |
у' |
Kt |
|
для второй производной — |
|
|
|
|
|||
М у - |
v" = |
А ; |
к ,. |
и" - |
J |
l s l . |
|
у |
и" |
Мр |
|
У" |
|
Kt1 |
|
для п-й производной — |
|
|
|
|
|
||
|
v<"> |
_ |
Му |
/С<"> = |
и<п>_ |
/Су |
|
M v(n) =2L_ |
|
||||||
|
|
М,(п> ’ |
У |
yW |
|
A |
|
|
|
|
|
(63)
(64)
Взаимосвязь между операторами дифференцирования по ма- ‘ шинному и реальному времени следующая:
d .. |
(65) |
s = — - M t p |
(IT
и
(66)
Kt i t
30
Отсюда следует, что
s 2 ~~ М 2 р 2 и s n ~ M nt p n, |
(67) |
а также
р"
И S'
к? *7
Поэтому для производных'в операторной форме имеем:
M t |
9 |
А»? |
9,. |
8"и = МП^ > , |
ш = тг = /?v; |
S2U =т |
— |
Р 2У, |
|
Му |
|
Му |
|
Му |
8и '= %гРУ’ |
s 2u = |
|
р 2 у; |
s nu = 0 L p n у. ' |
At |
|
|
|
К? ■ |
(68)
(69)
(70)
Зная масштабные соотношения или масштабные коэффициенты, можно рассчитать значения коэффициентов передачи решающих устройств. Рассмотрим этот расчет на примере уравнений (42). и (51), а также уравнений в операторной форме (43) и (52); при этом будем пользоваться масштабными коэффициентами К.
С учетом соотношений (61), (62) и (63) уравнение (42) можно' преобразовать:
К* (Ра _ |
b0v |
а_^_ |
|
К t |
du |
|
|
||
K |
^ |
- |
A ^ - a °Ky |
|
' |
Р ’ |
|
|
|
ИЛИ |
|
|
|
|
|
|
|
|
|
— - = |
b 0 |
|
v — a 0 — |
а — |
— • |
— |
(71) |
||
(Н^ |
|
КхК2 |
Д2 |
|
1 /О |
dt |
|
||
Ура1внения (71) |
и (51) |
тождественны, если только |
|
||||||
?о—Ьо- |
|
|
Др |
и я. |
|
|
(72) |
||
Кх К? |
|
|
К( |
|
|||||
|
|
|
|
|
|
|
|||
Используя соотношения (70) |
и (61), можно преббразовать урав |
||||||||
нение (52): |
|
|
|
|
|
|
|
|
|
—Ъ ь К х Х - н К у У — Ч - Р У ,
К* Kt
или
|
к х к ] |
(73) |
Р2У = % |
х - * о К ] у — *x K t p y . |
|
Ку |
|
5* |
31 |
Уравнения (73) и (43) тождественны, если только
К |
к * |
74) |
60 = Р о — |
-К* 11а о «-=1 я о |
Л у
Из соотношений (74) легко можно получить соотношения (72)- При моделировании линейного дифференциального уравнения третьего порядка можно было бы получить вместо машинного урав
нения (52) уравнение вида
s3и — $0v — |
а0и — ajS« — |
а.^2и |
|
(75) |
||
и масштабные уравнения |
|
|
|
|
|
|
Ро |
ап |
«1 |
И as |
Д» |
(76) |
|
Kt |
Л7 |
|||||
Kx K*t |
|
|
|
Сопоставляя соотношения (72) и (76), нетрудно установить пра вила, по которым составляются масштабные уравнения при моде лировании линейного дифференциального уравнения с постоянными коэффициентами п-го порядка: 1) число масштабных уравнений равно я+ 1; 2) вид этих уравнений следующий:
Ро = |
„ „ |
а1 |
. „ |
~ Дя— |
Л1 |
~ |
Уап—1 |
~ L . (77) |
|
|
|
К"-‘ |
|
|
На основании уравнения (72) и соотношений к уравнению (51) устанавливаем, что для структурной схемы набора, изображенной на рис. 9, в, коэффициенты передачи решающих устройств опреде ляются следующими уравнениями:
M « = 6 o ~ V ’ |
(78) |
KxK2t |
|
В этих уравнениях величины Ьо, ао и ai заданы исходным уравнени ем и определяются по соотношениям к формуле (42). Масштабные коэффициенты Ку, Кх и Kt приблизительно определяются по форму лам (59). Задаваясь любыми двумя коэффициентами передачи, ос тальные три коэффициента рассчитывают по уравнениям (78). Коэффициенты передачи интегрирующих устройств мопут быть в пределах от 0 до 10 и для суммирующих и инвертирующих уст ройств от 0 до 20 ( в случае параметрической компенсации дрейфа нуля, например для УПТ в модели МН-7). При автоматической ста билизации нулевого уровня коэффициенты передачи могут быть со ответственно от 0 до 20 и от 0 до 100. Не следует выбирать очень малые значения коэффициентов передачи, поскольку установка их на АВМ затруднительна.
&
§ 4 МЕТОД ИССЛЕДОВАНИЯ НЕОДНОРОДНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ, ИМЕЮЩЕГО ПРОИЗВОДНЫЕ В ПРАВОЙ ЧАСТИ
Описанная выше методика не может быть непосредственно при менена для подготовки модели процесса, представляющей собой ли нейное дифференциальное уравнение с постоянными коэффициен тами, правая часть которого содержит производные от независимой
переменной х: |
. |
РпУ + ап- 1 р" - 1у -I-------\-а,ру -\-а0у='Ьтртх ~ |
|
-I- Ът- \ р т ~ 1 х ^ ... hbjpx -f Ь0х, |
(79) |
где р = — и т<^п. Такими дифференциальными уравнениями dt
описываются динамические системы с дробно-рациональными пере даточными функциями, работающие при любом заранее неизвест ном входном сигнале x(t), т. е.
bmPm+ b m_ l p ’n - i + ... |
+ b1p + b0 |
W{p) = |
(80) |
p n+ a n- i p n~ 1Н - . . . + а , / м а,
Бели уравнение (79) представить в виде системы машинных диф ференциальных уравнений первого порядка, т. е. использовать ме тод представления уравнения (42) и его машинного уравнения (51) двумя уравнениями (54) и (55), то в одном из них получим произ водные от независимой переменной х. Непосредственное моделиро вание этого уравнения посредством формирования правой части дифференцирующими устройствами практически невозможно ввиду сильного увеличения помех входного сигнала x(t). Именно поэтому уравнение (79) нужно представить такой системой уравнений, в структурных схемах набора которых использовались бы только ин теграторы, сумматоры и инверторы.
Рассмотрим наиболее эффективный из известных методов тако го представления — метод непосредственного интегрирования, или метод канонической формы [7]. Сущность этого метода заключается в том, что уравнение (79) при п — т преобразуется в систему урав нений канонической формы для новой системы переменных Zi, гъ t.., zn, получаемых по уравнению
|
pzk zk+i — an_ky+-b„-kX, |
(81) |
где k~ 1, 2, |
n, a Zn+1=0 |
|
Действительно, если порядок правой и левой частей уравнения (79) одинаков, то после деления уравнения на рп и решения его относительно у получим:
У = Ьпх + г[<-«л—1 |
У + bп—I х) + — ] ( — ия—гУ + |
|
-Ь Ь,1—2 х) -f- |
+ |
(— (iiy +■ Ьхх ) + |
Ч----(— апУ ~b Ь9х |
(82) |
|
Р |
|
° ] f -}]■ |
|
|
33