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

книги / Моделирование систем управления

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
6.82 Mб
Скачать

Пуассона р. При данной конфигурации формы для определения всех ее раз­ меров достаточно задать один некоторый характерный размер /. Величина внешней нагрузки определяется силой F . Задача моделирования - выясне­ ние значения напряжения а, возникающего в некотором элементе конструк­ ции.

[/] = I,[F ] = IJWT-2, [е ] = Е~'МТ~2,

[у]=1Г2МТ-2, [а] = Г ,ЛУГ-2,[ц] = 1.

По известному алгоритму определяем критерии подобия:

EL2

у/3

о /2

* 1 - —

. * 2 = ^Г >

* Э = ( £ р *4=Ц-

Если модель и оригинал выполнены из одного материала, то значения р и Е у них одинаковы. В таком случае для постоянства TCI необходимо выполнение условия

/2

— = const.

(3.52)

Это означает, что при моделировании внешнюю нагрузку надо изменять пропорционально квадрату линейных размеров. При этом согласно выраже­ нию для я3 должно быть выполнено условие а = const. Это, в свою очередь, значит, что для подобных состояний модели и оригинала напряжения в сход­ ственных точках одинаковы. Однако согласно выражению для я2 при усло­ вии (3.52), должно быть выполнено равенство

у/ = const,

(3.53)

что невозможно (удельный вес у модели и оригинала один и тот же). Представим удельный вес в виде произведения у = рg, плотности р

и ускорения силы тяжести g . Тогда вместо (3.53) должно быть gl = const. При постоянном g = 9,81 м/с2 моделирование невозможно. Изменить

g можно только искусственным путем, установив модель небольших разме­ ров на специальную центробежную машину. При большом радиусе вращения модели центробежные силы инерции можно считать параллельными и анало­ гичными силам тяжести, но с другим ускорением. Выбором соответствую­ щей угловой скорости можно получить значение, большее и меньшее 9,81 м/с2.

Несколько слов о характерном размере. Конкретизация этого размера может быть весьма разнообразной. В зависимости от формы соответствую­ щего объекта в качестве «характерного размера» можно выбрать длину /, ширину b, высоту Л, диаметр d, радиус г и т.п. Допустимость такого про­ извола при физическом моделировании объясняется тем, что в любом случае отношение сходственных «характерных размеров» оригинала и модели опре­ деляет один и тот же масштаб геометрического подобия

Задача 5. Пользуясь системой относительных единиц, определить критерии подобия переходного процесса /(/) в цепи, образованной последо-

вательным соединением элементов с сопротивлением R и индуктивностью L, которая включается на постоянное напряжение.

Функциональные зависимости, характеризующие процессы в двух та­ ких системах, имеют вид:

Ah, h, ru Lu щ) = О, Ah* h, Гъ La, щ) = 0.

Из пяти параметров, характеризующих данный процесс, .три являются независимыми. Пусть независимые параметры /, г, t. Соответствующие этим параметрам базисные величины (i\bt R\b , ht) можно выбрать произвольно. Тогда базисные величины и[Ь, Lu, определяются из равенств:

Щь ~ htfiibl

L\b = R\bhb-

В этом случае

 

А_ = Л_,

А =А

*16 Rib

JL * А , hь ht

Щь «2b

A = A .

L\b ^26

Определим дополнительное условие, наложенное на выбор базисных величин i\b, Ru, hь, игы Ly» при выполнении которого процессы будут по­ добными.

Определим критерии подобия для обоих процессов:

.0)

~ . D ,

,(2)-Л а_.

 

п] -

. D »

 

 

 

i2R2

Проведем преобразования:

,о) _ _к_.

_(2) _

п . *

"2 “ » . »

п 2

Rih

Щ_

«2

 

 

 

 

 

я (1) _ U\ _ и \ь _

и 2Ь

 

_ u 2

 

h b R 2b _

д(2) h b R lb

i\R\

*2 Ri

i2Ri

 

v-2ь

u2b

hb R\k

ilb Rlb

 

 

 

 

Аналогично можно получить

 

 

 

 

 

Œ0> -

J V

r

ь

 

 

n2

- n 2

 

 

L lb

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

 

Щь = infiib, L2b ~ Ribhb-

Базисные величины для

z2, R2i t2 могут быть выбраны произвольно.

Задача 6 . Обеспечить подобие систем. Пусть имеются три системы

(две электрических,

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

рис. 3.7, а - система 1 (цепь Z/Æ/C/), на рис. 3.7, б - система 2 (цепь £;С?)> на рис. 3.7, в - система 3 (изгиб нагруженной балки).

 

L2

 

 

 

RI

 

 

Cl

С2

 

 

а

б

в

 

 

Рис. 3.7

 

Уравнение для системы 1 (см. рис. 3.7, а) (цепь LJCJRJ)

 

 

dtf

Т' Л, Z.C,

(3.54)

 

 

для системы 2 (см. рис. 3.7, б) (цепь Z,2C2)

 

 

dt\ L2C2

(3.55)

 

 

для системы 5 (см. рис. 3.7, в)

 

 

 

 

Р

(3.56)

 

— г + ------v = 0.

 

л *

V o

 

Рассмотрим методику определения переменных масштабов на приме­

ре подобия систем (3.54) и (3.55).

 

 

Обозначим

 

 

 

 

^2 = ^ !>

Ог ~ FQQ\ »

 

где Ft и FQ неизвестны.

Подставим эти обозначения в исходное уравнение для 2-й системы (3.55).

p Q =о.

d (t,F f L2C2

Примем для времени постоянный масштаб, т.е. Ft = const = mt = — ,

тогда

2 1.2 mfdty

где £ = -

L\c \

Проведем преобразования:

 

 

 

 

 

 

 

 

 

 

 

 

d((FgQ,)\

 

 

d Q , d F e

 

 

 

 

 

 

 

 

dt

 

^ л Г + а ^ Г ’

 

 

 

d'Q,

dFe dQ,

dFq dQ,

 

d %

 

 

 

F° ^ T

+l

ï

l

ï

+l

ï l

 

ï

+e' W

 

1 eQ'

 

d*Q,

,

2

dFe dQ,

f

r d 2Ff

 

a= o.

 

 

 

Or+ m2к

 

 

di\

 

F2

 

dt}

 

 

 

 

 

 

 

Чтобы последнее уравнение стало почленно равно уравнению (3.54),

нужно обозначить

 

 

 

 

 

2dF,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.57)

 

 

 

 

 

 

 

 

 

— = Г

 

 

 

 

 

 

 

 

 

 

dt\FQ

 

 

 

 

 

 

 

 

 

 

 

d %

 

 

2 к

1

 

 

(3.58)

 

 

 

 

 

 

 

— + т;к =

------

 

 

 

 

 

 

 

 

dt'F{

 

 

 

L.Cj

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

Из (3.57) находится переменный масштаб

 

 

 

 

 

 

 

 

й

-

. U

^

,

 

 

 

 

 

 

 

 

 

а

г

 

1

 

1

 

 

 

 

 

 

 

 

 

\aFg = atu

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ’

 

 

 

 

 

 

 

 

 

 

F g = m e ee\

 

 

 

 

где mQ_ а

- постоянная величина, зависящая от начальных условий.

Ô2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнения (3.58) с учетом (3.57) можно найти критерий подобия %\:

 

 

 

 

 

 

 

 

 

т*к

 

 

 

 

 

 

 

 

 

 

я ,= -_1__г2

 

 

 

 

 

 

 

 

 

 

L\C\

4

 

 

 

 

Критерий %\ позволяет установить, что масштаб времени имеет

значение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т, =

 

1

-а Ч - =

 

1

г

L2 C2 .

 

 

 

 

— - - а

 

 

 

 

 

 

 

I.C ,

 

 

 

Особенностью критерия Я| является наличие коэффициента преобра­ зования, устанавливающего зависимость между временем t\ в одном уравне­ нии и временем t2 в другом.

Задача 7. Проделайте опыт с системой 1 и далее, перестроив получен­ ную кривую Qi = Д/|), с помощью масштабов т, и mg найдите зависимость Qi

которая определяет поведение системы 2.

Параметры системы /: С\ = 2-КГ3 Ф, L\ = 2,63 Г, R\ = 15 Ом;

Параметры системы 2: С2 = 5-10"3 Ф, L2 = 2 Г; e i( 0 ) = e 2(0).

Задана 8 . Определить масштабы на примере подобия систем (3.54) и (3.56). В системе (3.56) аналогом времени является изменение геометриче­ ской координаты x\= m tti, a аналогом заряда Q\ - прогиб балки у. При этом

y= QlmQeat.

3.9.Задачи для самостоятельного решения

Задача 1. Определить сходственные точки процессов /j(/|), i2(t2)

и/з(*з) (рис. 3.8), полагая соответственно

 

 

 

Ri = 10 Ом, L\ = 20

Г, U\ = 100 В;

 

 

 

 

R2 = 20 Ом, L2 =40 Г, U2 = 75 В;

 

 

 

 

 

R3 = 90 Ом, = 60

Г, С/3 = 500 В.

 

КО

С /.

--

Z,

 

+ 0

ял

0 - е

х ), т = ------постоянная времени.

»,— *

 

R

 

R

 

 

^

Л

г/

------------

Порядок расчета

Рис. 3.8

1.Определить критерии подобия.

2.Определить масштабные коэффициенты и степенные комплексы.

3.Определить сходственные точки и заполнить таблицу (табл. 3.2).

4.Построить зависимость %\ =Дп2).

Таблица 3.2

/

0,33

0,5

0,67

1,00

1,33

1,50

2,00

2,5

3,00

4,00

Я1(,)

-

-

-

2,00

1,50

1,33

1,00

0,80

0,67

0,50

*1*

 

 

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

 

 

0,17

*.<3>

-

-

-

0,67

0,50

0,44

0,33

0,27

0,22

 

-

-

-

2,54

2,06

1,89

1,58

1,40

1,29

1,16

 

 

 

 

 

 

 

 

 

 

-

*2<3)

-

-

-

1,29

1,16

1,12

1,05

1,02

1,01

1,00

 

Задача 2. Определить сходственные точки для процессов

Uiih)

для ЛС-цепи, Ri = 10 кОм, R2 = 50 кОм, С\ = 5 мкФ; С2 = 10 мкФ, U\ “ Ю В,

С/2=Ю В .

 

 

Задача 3.

Определить условия,

при которых сходственные функции

У\ х \ \ Уг ~ 8я22 подобны.

2

Задача 4.

Определить сходственные точки для функций у\ = xi + 4,

^2= &с22 + 6.

 

 

Задача 5. Определить критерии подобия для случая, когда имеются два тела массой М\ и М2) движения которых подобны и описываются одно­ родными уравнениями:

/ г £ к

А

. M É% 0.

(3.59)

dlf

 

Ml dû

 

Задача 6 . Определить критерии подобия переходного процесса рас­ пространения волны напряжения по длинной линии при включении ее на по­ стоянное напряжение.

Рассматриваемый процесс описывается уравнением

CL^

+CR ™ +W -d Ï +R G U ^

dt2

à

dt

dl2

где C, L, R к G - соответственно

емкость,

индуктивность, сопротивление

и проводимость на единицу длины; U- напряжение вдоль линии; / -время; /-длина»

Задача 7. Определить критерии подобия электрической цепи с взаи­ моиндукцией Ми, которая состоит из двух неподвижных относительно друг друга контуров, образованных последовательно соединенными элементами с активным сопротивлением R и индуктивностью I , причем второй контур короткозамкнут, а первый включен на постоянное напряжение U\ (рис. 3.9). Известно, что процесс в такой цепи описывается системой

(/,= /,* ,

^

+м 12^

,

 

О =i2R2 +L2 ^ - M n ^ .

 

2

2 ^

dt

12 dt

 

Задача 8 . Определить критерии подобия для процесса вынужденных

М12

 

механических колебаний в вязкой сре-

 

де:

груз

массой

М колеблется под

R\

 

действием

возбуждающей силы:

Ui

Ri

FB= Fsm(ù t на пружине жесткостью С;

при перемещении груза на расстояние

fi-

 

/ в вязкой среде появляется сила со*

Рис. 3.9

 

противления Fc,

пропорциональная

скорости V(FC = -kV ).

Задача 9. Рассмотрим гидродинамическую систему (рис. 3.10). Под­

водная лодка с

характеристическим размером d

движется с различной

 

 

скоростью в вязкой жидкости, испытывая

 

 

силу лобового сопротивления D.

 

I I >

0

Размерности параметров:

 

Скорость жидкости

У

19'*

V

 

Характеристический размер d

L

Р»М* ------ЧПВР \

Плотность жидкости р

ML'3

lilllMlil)(ПТЛШТТТГ

Вязкость жидкости р

м е - 'г ’

D-сила

 

Рис. ЗЛО

 

Сила лобового сопротивления D

M L9'3

 

 

 

 

Задача 10. Рассмотрим систему двух тел, вращающихся в открытом космосе друг относительно друга благодаря взаимному притяжению. Масса 1-го тела М\. Масса 2-го тела М2. Расстояние между телами R . Период об­ ращения Т.

Задача 11. Для оценки взаимных помех между различными радио­ электронными средствами, размещаемыми на судах, можно воспользоваться подобным физическим моделированием электромагнитной обстановки, кото­ рую определяют семь величин: частота/ [ Г -1], характерный размер / [L], на­ пряженность магнитного поля Я [ZT1/], напряженность электрического поля Е [ЬМТ~Ъ1~Х\, магнитная проницаемость е [L~3M~lT 4I2], удельная проводимость среды a [L~3M~lT 3I 2].

Задача 12. Обеспечить подобие зависимостей iii(#i), Ц2(Я2), р3(Я3) (рис. 3.11).

'

 

Рг'

 

 

k

 

Р * ‘

 

 

 

 

 

 

Рзо

 

 

 

Ц

Ч

"

 

1

 

1

 

 

 

1

 

1

 

 

1

 

 

 

1

 

1

 

 

 

 

 

1

 

1

 

 

1

 

ь

 

1

 

1

 

 

Ню

Hi

 

 

Нг

 

'

Н*

 

 

Яго

Язо Я 3

 

 

 

 

 

Рис. 3.11

 

 

 

Численные значения р и Яследующие:

 

 

 

Ц1

1

4

6

3

1

 

 

 

я,

2,5

5,5

7,0

8,5

10

 

 

 

я2

5

9

18

9,5

3,5

 

 

 

10

15

25

30

35

 

 

 

Цз

1,5

2,5

3

2,2

0,8

 

 

 

Я 3

6

8,5

11

14

17

 

 

 

Для кавдой из зависимостей вводятся соответствующие базис­ ные значения р/0, Щ 9 после чего характеристики строятся в новой системе

координат \i = ф(Я*), где р* = — и Н* = - ^ - . И/о Ню

4. ОСНОВНЫЕ ПОДХОДЫ К ПОСТРОЕНИЮ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ СИСТЕМ

Математическая схема - звено при переходе от содержательного к формальному описанию процесса функционирования системы с учетом воздействия внешней среды.

Имеет место цепочка «описательная модель - математическая схема - математическая (аналитическая или (и) имитационная) модель».

При построении математической модели системы необходимо решить вопрос о ее полноте. Также должна быть решена задача упрощения модели, которая помогает выделить основные свойства системы, отбросив второсте­ пенные, причем отнесение свойств системы к основным или второстепенным зависит от цели моделирования (анализ вероятностно-временных характери­ стик процесса функционирования системы, синтез структуры системы и т.п.).

Процесс функционирования системы можно представить в виде мно­ жества величин:

1. Совокупность входных воздействий на систему

х{е Х

2. Совокупность воздействий внешней среды vr е V.

3. Совокупность внутренних (собственных) параметров системы hk е Я.

4. Совокупность выходных характеристик системы

Я е Г .

Можно выделить управляемые и неуправляемые переменные.

Входные воздействия и внутренние параметры являются независимыми (экзогенными) переменными:

v ( 0 = (Vi(0> V 2(0 * ...V e ( 0 v ..) »

Выходные характеристики являются зависимыми (эндогенными) пере­ менными.

y(f) = O'!(0.^2(Оэ-эЗ'у(0.-).

y(t) = F (x,v, h, t) - закон функционирования системы. Это запись для динамической модели. Для статической модели закон функционирования вы­ глядит следующим образом: y(t) = F (xiv)h) . Эти законы могут быть заданы аналитически (формулами), графически, таблично и т.д.

Таким образом, математические модели систем описываются совокуп­ ностью входных воздействий х, е Х\ воздействий внешней среды vr е V; внутренних (собственных) параметров системы hk е Я; выходных характери­ стик системы yj е Y.

Модели делятся на детерминированные (не содержат элементы случай­ ности) и стохастические.

Приведенные математические соотношения представляют собой мате­ матические схемы общего вида и позволяют описать широкий класс систем.

4.1. Линейные регрессионные модели для одной переменной управления

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

4.1.1.Линейная регрессионная модель

Пусть имеется п выборок Xh У,.

Допустим, модель представляет собой прямую Ум = А0 +А\Х. Функцию

ошибки запишем в виде F = £ (Л 0 + А{Х ( - У))2 . Наилучшие оценки для ко­ эффициентов AQи А] получаются из системы

ал

^ = о

ЗА,

Решение системы уравнений дает значения для AQnA\:

i r i t x f - î x . î w

n i x f t - i x & Y i

л = — --------- 1— 4 -

Обычно мерой ошибки регрессионной модели служит стандартное (среднеквадратическое) отклонение S:

Для нормально

распределенных процессов » 68 % точек находится

в пределах ± S и 95 %

точек - в пределах ± 25. Большое значение S может

означать, что модель не соответствует процессу, который послужил источни­

ком экспериментальных данных, либо имеется большой разброс данных из­ мерений.

4.1.2.Преобразование линейной регрессионной модели путем

привязки к средним значениям

Дано п измерений Xh Yr

 

 

Определим средние значения X ,

Y и отклонения экспериментальных

данных от средних значений.

 

 

 

 

 

? ш Ш ,

 

х , - Х , - Х ,

п L

 

у,ш Г ,- Г .

 

Регрессионная модель в новых координатах имеет вид

Функция ошибки

Л . = 4 * -

 

 

 

 

 

F = I ( 4 * f -V i)2 -

 

Наилучшая оценка рдяА\ определяется из уравнения

 

9F

К л ,* , - у, )*. = °-

 

— = о,

 

CVij

 

 

 

Подставив А\ в уравнение, получим

 

Y = Y + Al( X - X ) = A0 +AlX i

 

 

AQ = Y - A }X .

 

В тех случаях, когда определитель матрицы »

Y.X, мал, преоб­

 

 

Ъ х ,

I x f

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

4.1.3.Достоверность регрессионных моделей

Требуется установить, является ли гипотеза, положенная в основу рег­ рессионной модели (о том, что коэффициент А\ = 0, или, что тоже самое, име­ ется наилучшая прямая) достаточно правильной.

Анализ точности и достоверности регрессионных моделей можно про­ водить с помощью следующих характеристик, учитывая, что число степеней свободы определяется соотношением v = w -& - 1 ,г д е и - число эксперимен­ тальных данных, к - число переменных управления:

- разброс данных вокруг их среднего значения (рис. 4.1, а)

S S , = Z ( V - Ÿ ) 2.

Число степеней свободы п - 1; - разброс данных вокруг линии регрессии (рис. 4.1, б)

S S 6 = £ ( F, - F , ) \

Соседние файлы в папке книги