Файл: Севостьянов, А. Г. Основы математического моделирования механико-технологических процессов текстильной промышленности из цикла лекций заочного факультета по технической кибернетике в текстильной промышленности.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