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