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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 31.10.2024

Просмотров: 68

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Заменяя суммы в скобках соответствующими переменными <2, пре­ образуем уравнение (82) в систему из п уравнений канонической

формы (83) и дополнительное уравнение связи

(84):

k = n

pzn = Zn+1— «оУ + b0x -= - ~ d0y +

b 0x

при

PZn- 1 “ zn axy -f- bxx

при

k -

n -

1

 

 

при

 

(83)

pz-i = zz — an- 2 У + bn—2 X

 

k = 2

pz 1 = 4 dn—i у -\~ bn- 1*

 

при

k =

1

Zx = У - bn x.

 

 

 

(84)

Рассматривая уравнение (83), легко заметить, что характер их построения соответствует формуле (81). Уравнения (83) модели­ руются по структурной схеме, содержащей п интеграторов и не более трех инверторов, причем интеграторы имеют не более трех входов. Используем этот метод при исследовании следующего урав­

нения:

(85)

РгУ + aiРУ + «о У= bipx + b0x,

где аг=3, а0=4, bi—0,5 и Ь0~2. Начальные условия при исследова­ нии такие: г/(0) —2,5; рг/(0)=25; *(0) = 1; *тах=2. Входная перемен­ ная х поступает от внешнего источника. Исследование ведется в

реальном масштабе времени,

т. е. /С* =

1 -

 

 

 

(85),

при

Система уравнений (83),

эквивалентная уравнению

п= 2 и т —Г имеет вид

 

 

 

 

 

 

 

 

рг2 = гъ — а0у +

b 0x = а0у -f Ь0х

J

 

(86)

p z x= z2 — аху +

Ьхх

 

 

Г

 

 

 

где г3= 0, согласно условию (81).

(85) член

р2х

отсутствует,

Так как в правой части уравнения

т. е. &2=0, уравнение связи (84) имеет вид

 

 

 

 

 

Zi = У — Ь2х = у.

 

 

 

 

(87)

Учитывая это, преобразуем второе из уравнений

(86)

к виду .

 

ру = z2 а ху +

Ьгх,

 

 

 

 

(88)

где а2= 1 . Легко показать,

что исходному уравнению

(85) экви­

валентны два уравнения — первое из уравнений

(86)

и уравнение

(88). Структурная схема набора этих уравнений

показана

на

рис. 11. В схеме применяются два интегратора

(1

и 3)

и два

ин­

вертора {2 и 4).

 

 

 

 

 

 

 

 

34


Пользуясь указанной структурной схемой набора, можно на­ писать следующие соотношения:

 

 

 

и3 = — у ( M i +

k xv );

 

 

 

(89)

 

 

 

 

и, = —*3«S-.

 

 

 

 

 

(90)

 

 

щ = — - ( fe5«2 + M i

+

M );

 

 

 

(91)

 

 

 

 

и — knU.

 

 

 

 

 

 

(92)

После подстановок из уравнений

(90)

и

(92)

в уравнения

(89).

и (91)

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5Иг = — а0и +

%v\

 

 

 

 

(93)

 

 

 

=

а2ы2 — а,ы 4-

 

 

 

 

 

(94)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

« О

=

k б’'>

аа 12

==

^

7

-

f e l l Рs о'P. “i =

М

 

7 - 95)

 

 

«7

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя масштабные коэффициенты

Ку — ~,

Кх =

.

и

=

преобразуем уравнения

(93)

и (94):

 

 

 

 

 

г2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р * 2 =

-

«

Аг

о

£ .

 

У

+

Р о £

? *

;

98)

 

 

 

 

 

 

 

 

А*

 

 

 

 

 

 

 

Р У

=

“ 2

 

7

 

* 8

Yа х1 -У

+

P i

 

97)

 

 

 

Л у

 

 

 

 

 

 

А у

 

 

 

 

Уравнения (96) и

(97)

соответственно

эквивалентны

первому

из уравнений (86) и уравнению (88) при

наличии

следующих

соотношений:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(98)

На выходе интегратора 1 получается напряжение, соответст­ вующее новой переменной, но с обратным знаком, т. е. — z%. Поль­ зуясь уравнением (88), находим начальные условия по формуле

— 22 (0) = ру(0) + 3 i/(0) — 0,5х(0) = 25 + 3 • 2,5 — 0,5- 1 = 32.

35


Значения масштабных коэффициентов определяются с учетом следующих условий:

1) условие максимального значения входной величины

■^тах—2, Т. е.

Кх <

=

° , 5 0,25,' 1 [ м ' е ' ]

=

-'-max

^

 

 

где м. е- — машинная единица,

равная

100 В для

ABM МН-7

(с ламповыми ОУПТ); 2) начальное условие для машинной переменной и2, т. е.

“ а(0) = ^ z 2l 0 ) > 1 [м.е.];

отсюда

*, < ^ 0 , 0 3 ;

3)условие выравнивания машинных коэффициентов, т. е. cto=cc2- Подставляя соотношения (98), получаем:

отсюда

 

 

Кг

= 2;

= 2Кг = 2-0,03 = 0,06.

 

'

Подставляя

найденные

значения масштабных коэффициентов

в соотношения

(98), найдем:

0,03

2;

•э; «г

0,06

 

 

. 3„ .= 2 •. 0,03

—0,24;

?! = 0,5-

0,25

 

 

0,06

0,03

0,06

0,12.

0,25

Используя соотношения (95) и принимая Кз=Кт=\, найдем значения коэффициентов передачи по каждому входу решающих устройств:

fei = Зо = 0,24; /г2 = а0 = 2; Л* = ?, = 0,12;

fe6 = ч ■= 2; fee = а1 = 3.

Подставляя найденные значения коэффициентов в уравнения (93) и (94), получим систему машинных уравнений в явном виде:

su2= — 2и + 0,24 v

|

(99)

su *= 2м2 — Зи + 0,12 v /

-36


Начальные условия при 1

м. е.= 100 В следующие:

 

и,(0) =

- Ку у(0) =

-

0,06-2,5-'= — 0,15 [м. е] =

- 15 Вг

«з(0)

К 2г2ф) г

-0 ,0 3 -3

2 = - 0,96_ [м. е] =

— 96 В;;

г'(0) == К х *(0)

= L0,25

-1 = к0,25

[м.*Ге] = 25 В;

 

*W “ Кх л'тах = 0,25-2 = 0,5

|м. е| = 50

В.

§ 5 МЕТОД ИССЛЕДОВАНИЯ НЕЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

При исследовании нелинейных обыкновенных дифференциаль­ ных уравнений используется такая же методика, как и при иссле­ довании линейных уравнений. Введением вспомогательных пере­ менных система нелинейных дифференциальных уравнений преоб­ разуется в систему линейных дифференциальных уравнений и до­ полнительных нелинейных уравнений связи между вспомо­ гательными переменными и исходными. Нелинейные уравнения связи воспроизводятся на АВМ нелинейными блоками, при ,этом: масштабные коэффициенты вспомогательных переменных зависят не только от масштабных коэффициентов исходных переменных, но и от коэффициентов нелинейных преобразований нелинейных блоков.

Если

уравнение связи

представляет операцию

умножения

у = ах1*2 ,

то, переходя к

машинным переменным и = Куу, Vi-KXi

и v2

= Кх,х2>получим машинное уравнение

 

 

 

 

аКу

ТЧТ/а.

(100)-

 

 

И = ----- ----

 

 

 

КХ1 Кх„

 

 

Для

блока перемножения,

имеющего

коэффициент

передачи ku,.

связь между входными щ и v2 и выходным и переменными опре­ деляется уравнением

u ^ - k „ v t v 2.

(101>

Уравнения (100) и (101) эквивалентны при условии

 

К у = —1КХ\ КХ2.

(102)-

а

 

Величина kn определяется конструкцией БП; в машинах с лам­ повыми усилителями со шкалой 100В £п = 0,01, в полупроводни­ ковых АВМ со шкалой 50В &п=0,02 и со шкалой ЗОВ Ап=0,03.

При необходимости изменения величины kn на АВМ, имеющих диодные блоки перемножения, меняют величину сопротивления в цепи обратной связи усилителя, суммирующего выходные токи квадраторов (см. рис. 6).

37


Ранее указывалось, что для осуществления операций деления,

извлечения корня и возведения

в

степень можно

использовать

блок перемножения в цепи обратной связи усилителя

(см. табл. 2,

строка 3).

 

 

 

Если уравнение связи

 

 

(103)

у = а — ,

 

Х

9

 

то машинное уравнение имеет вид

 

 

Ку«х2

 

(104)

и. = а

 

 

Для решающего устройства, осуществляющего операцию деле­

ния (см. табл. 2, строка 3), можно написать:

 

и =

k Ay_t }

(105)

где /гд — коэффициент передачи делительного устройства.

Уравнения (104) и (105) эквивалентны при условии,

что мас­

штабный коэффициент

 

а Кх2 .

0 об)

где &д=10. Для схемы, избраженной в табл. 2 (строка

3, схема

б), можно написать:

 

Переходя к исходным переменным и при сопоставлении с уравне­ нием (103), получаем формулу для определения соотношения соп­ ротивлений в схеме делительного решающего устройства:

 

Л'уА'

(Ю7)

 

ak

 

 

Ri

\хг

 

Масштабные коэффициенты для схем извлечения корня и воз­

ведения в степень определяют аналогичным образом.

функции

Если уравнение связи

задается

в виде нелинейной

то масштабный коэффициент для нее находят по формуле

7S ____wmax_

1,10

(108)

/l

/i(-T)ma

ZlWn

 

Для иллюстрации подготовительных операций к исследованию системы нелинейных дифференциальных уравнений на АВМ рас­ смотрим следующие два уравнения:

РУг = а2уа —

а,у!у2 -

aty\

|

РУг = сцУг +

«з У1У2 +

а\ У\

X109)

(здесь а2=Ю ; а3=2; а4=0,5; а1= 0 г 1; ^3=0,05; а^=0,1).

38