книги / Моделирование систем управления
..pdfПуассона р. При данной конфигурации формы для определения всех ее раз меров достаточно задать один некоторый характерный размер /. Величина внешней нагрузки определяется силой 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- |
В этом случае |
|
А_ = Л_, |
|
hь |
hь ’ |
А =А
*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
где £ = -
Проведем преобразования: |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
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 , ) \