Моделирование систем
..pdfСоставляем уравнения для определения неизвестных xi , yi , z i :
1. Для критерия 1 |
|
|
|
|
|
[L] [π1 ][i x1 U y1 t z1 ] . |
(3.51) |
||||
Уравнение размерности для (3.51) имеет вид |
|
||||
|
M 1 T |
2 |
I 2 |
|
|
L2 |
|
|
|
||
|
|
|
|
|
|
[1] L0 M 0 T 0 I1 x1 L2 M 1 T 3 I 1 y1 L0 M 0 T1 I 0 z1 . (3.52)
Составляем уравнения для каждой основной единицы измерения:
2y1 2; |
{L} |
|
|||
|
1; |
|
{M} |
|
|
y1 |
|
(3.53) |
|||
|
|
|
|
{T} |
|
3y1 z1 2; |
|
||||
x |
y |
1 |
2. |
{I} |
|
1 |
|
|
|
|
Так как 1-е и 2-е уравнения зависимые, то получаем следующую совместную систему:
y1 1; |
|
|
|||
|
|
|
|
|
(3.54) |
3y1 z1 2; |
|||||
x |
1 |
y |
1 |
2. |
|
|
|
|
|
||
Решение системы (3.54): |
|
|
|
||
|
|
x |
3, |
|
|
|
|
1 |
1, |
(3.55) |
|
|
|
y1 |
|||
|
|
|
1. |
|
|
|
|
z1 |
|
2. Для критерия 2
|
x2 |
U |
y2 |
t |
z2 |
(3.56) |
[Ñ] π2 i |
|
|
|
|||
91 |
|
|
|
|
|
|
Уравнение размерности для (3.56) имеет вид
|
|
|
|
|
|
|
|
2 |
M |
1 |
T |
4 |
I |
2 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.57) |
||||||
[1] |
0 |
0 |
T |
0 |
I |
1 |
x2 |
|
2 |
|
1 |
T |
3 |
I |
1 y2 |
0 |
0 |
T |
1 |
I |
0 |
z2 |
|||||
|
. |
||||||||||||||||||||||||||
L M |
|
|
|
|
L M |
|
|
|
|
|
|
L M |
|
|
|
|
Составляем уравнения для каждой основной единицы измерения:
2y2 2; |
{L} |
y2 1; |
{M} |
3y2 z2 4; {T} x2 y2 2. {I}
Совместная система имеет вид
y2 1;
3y2 z2 4;x2 y2 2.
Решение системы (3.59): |
|
x2 |
1; |
y2 |
1; |
z2 |
1. |
3. Для критерия 3
[R] 3 i x3U y3t z3 .
Уравнение размерности для (3.61) имеет вид
L2M 1T 3I 2
[1] L0M 0T 0 I1 x3 L2M 1T 3I 1 y3 L0M 0T1I 0 z3 .
(3.58)
(3.59)
(3.60)
(3.61)
(3.62)
92
Составляем уравнения для каждой основной единицы измерения:
|
|
|
|
2y3 2; |
|
|
|
|
|
|
|
|
|
|
{L} |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
y3 |
1; |
|
|
|
|
|
|
|
|
|
|
{M} |
|
|
|
|
|
|
|
|
|
(3.63) |
||||||
|
|
|
|
3y3 z3 3; |
|
|
{T} |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
x3 y3 2. |
|
|
|
|
{I} |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Совместная система имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
y3 |
1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
3y3 z3 3; |
|
|
|
|
|
|
|
|
|
|
|
(3.64) |
||||||||||||
|
|
|
|
|
|
|
|
x3 y3 2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Решение системы (3.64): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
x3 1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
y3 |
|
1; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.65) |
|||||
|
|
|
|
|
|
|
|
z3 0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
4. Для критерия 4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
x4 |
U |
y4 |
t |
z4 |
|
|
|
|
|
|
|
|
|
|
(3.66) |
||||
|
|
|
|
|
[ω] π2 i |
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|||||||||||||
Уравнение размерности для (3.66) имеет вид |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
0 |
T |
1 |
I |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
L M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.67) |
|||||||
0 |
0 |
|
0 |
|
1 |
|
x |
2 |
|
|
1 |
|
|
3 |
|
|
1 |
y |
4 |
0 |
0 |
|
1 |
|
0 |
|
z |
4 |
||||
T |
I |
|
4 |
|
|
T |
|
I |
|
T |
I |
|
. |
|||||||||||||||||||
[1] L M |
|
|
|
|
|
L M |
|
|
|
|
|
|
|
L M |
|
|
|
|
|
Составляем уравнения для каждой основной единицы измерения:
2y4 0; |
{L} |
|
||||
y4 |
0; |
|
{M} |
(3.68) |
||
3y4 z4 1; |
{T} |
|||||
|
||||||
x |
y |
4 |
0. |
{I} |
|
|
4 |
|
|
|
|||
|
|
|
93 |
|
|
Совместная система имеет вид |
|
|||||||||||
|
y |
4 0; |
|
|||||||||
|
|
|
z4 1; |
(3.69) |
||||||||
|
3y4 |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||
|
x4 y4 0. |
|
||||||||||
Решение системы (3.69): |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
x |
4 0; |
|
|||||
|
|
|
|
|
|
|
0; |
(3.70) |
||||
|
|
|
|
|
y |
4 |
||||||
|
|
|
|
|
|
|
1. |
|
||||
|
|
|
|
|
z4 |
|
||||||
Таким образом, в соответствии с (3.50) критерии име- |
||||||||||||
ют вид |
|
|
|
L |
|
|
|
|
|
|
|
|
π1 |
|
|
|
|
|
|
Li3U 1t 1; |
|
||||
|
|
|
|
|
|
|||||||
|
|
|
i 3U 1t1 |
|
||||||||
π2 |
|
C |
|
|
|
|
|
Ci 1Ut 1; |
|
|||
i1U 1t |
1 |
|
(3.71) |
|||||||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
R |
|
|
|
|
|
|
|
π3 |
|
|
|
|
|
RiU 1; |
|
|||||
i 1U1t 0 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|||||
π4 |
|
|
ω |
|
|
|
|
|
ωt. |
|
||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
i 0U 0 t 1 |
|
Рассмотренные положения относятся к случаю заведомо подобных процессов, т.е. определяют необходимые условия существования подобия. В связи с этим возникает вопрос об условиях, не только необходимых, но и достаточных для существования подобия. Такие обстоятельства, кроме равенства критериев, включают в себя требования подобия начальных и граничных условий сопоставляемых процессов.
Положения относительно необходимых и достаточных условий подобия изложены в основных теоремах теории подобия.
Теория подобия включает в себя три основные теоремы:
первая теорема подобия (теорема Ньютона– Бертранса) о небходимом условии подобия;
94
вторая теорема подобия ( -теорема); третья теорема подобия (теорема Кирпичева–
Гухмана) о необходимом и достаточном условии подобия.
3.7. Применение первой теоремы подобия для определения критериев подобия
Для того чтобы явления были подобны в том или ином смысле (полно, приближенно, физически, математически и т.д.), определенное сочетание параметров, называемых критериями подобия, должно иметь сходственную алгебраическую форму и одинаковые численные значения (без доказательства).
Первая теорема формулирует необходимые условия существования подобия (одинаковые критерии подобия), но не указывает способы установления подобия и способы его реализации при построении моделей.
Первая теорема применяется в следующих случаях:
1.Определение параметров для двух сравниваемых объектов.
2.Выделение группы независимых параметров и определение количество критериев.
3.Получение выражения для критериев подобия и сравнение их для двух определяемых объектов.
4.Вычисление численных значений алгебраических форм критериев подобия, если сравниваемые объекты совпадают для определенных точек пространства параметров. Одинаковые численные значения критериев позволяют судить о подобии этих объектов.
3.7.1.Определение критериев подобия по уравнениям процессов
Возможны два случая описания исследуемого процесса: 1. Все члены уравнения – однородные функции параметров, определяющих протекание этого процесса, и их про-
изводных.
95
С п р а в к а. В математическом анализе функция f(x1, ..., xn) от n аргументов, определенная в области D, называется однородной функцией k-й степени, если при умножении всех ее аргументов на множитель M она приобретает тот же множитель M в k-й степени, т.е. если тождественно вы-
полняется равенство |
|
f(Mx1, ..., Mxn) = M k f(x1, ..., xn). |
(3.72) |
При этом все члены уравнения имеют общий множитель, который может быть вынесен за знак функциональной зависимости.
2. Часть членов уравнения – неоднородные функции параметров, не допускающие вынос за знак функции общего
множителя [например, sin ( t + ), exp (–t/ )].
3.7.2. Определение критериев подобия процессов, описываемых уравнениями, содержащими однородные функции
Рассматриваются необходимые для определения критериев подобия преобразования, одновременно пример преобразования уравнения переходного процесса i(t) в последовательной цепи из активного сопротивления R и индуктивности L, которая включается на постоянное напряжение u.
Имеются уравнения двух подобных процессов 0 и 0, являющихся функциями параметров P1, P2, ..., Pj, ..., Pn и R1, R2, ..., Rj, ..., Rn соответственно, т.е. в общем случае
m
1 i m i 0 , (3.73)
i 1
m
Ô1 ... Ô i ... Ô m Ô i 0, (3.74)
i 1
или
|
f |
P , ..., P |
, |
|
|
i |
i |
1 |
n |
|
(3.75) |
|
|
|
|
|
|
Ôi fi P1, ..., Pn . |
|
||||
|
|
|
|
|
|
|
|
96 |
|
|
|
Для рассматриваемого примера уравнения двух переходных процессов это два дифференциальных уравнения:
i1R1 + L1di1 |
– u1 = 1 + 2 |
+ 3 = i = 0, |
(3.76) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
dt1 |
|
|
|
|
|
|
|
i 1 |
|
||||
i2R2 + L2di2 |
|
– u2 = 1 |
+ 2 |
+ 3 |
= Ôi , |
(3.77) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
dt2 |
|
|
|
|
|
|
|
|
|
i 1 |
|
|
где |
|
|
|
|
|
|
|
|
|
|
,u , |
|
|
|
||
i R f i ,t , R , L |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1 1 |
1 1 1 |
1 1 |
1 |
|
|
|
|||||||||
|
|
|
|
|
|
|
di1 |
|
|
|
|
|
|
|
|
|
|
|
L |
|
f |
|
|
|
|
|
|
|
|
||||
2 |
|
|
2 |
i ,t , R , L ,u |
, |
|
|
|||||||||
|
|
|
1 |
dt1 |
|
1 1 |
1 |
1 1 |
|
(3.78) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|||||||||||||
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 u2 f3 i2 ,t2 , R2 , L2 ,u2 |
. |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сопоставляемые процессы 0 |
и |
0 подобны, |
следова- |
тельно, между их сходственными параметрами должны существовать соотношения пропорциональности вида
P1 = m1R1, ..., Pj = mjRj, ..., Pn = mnRn |
(3.79) |
или |
|
R1 = mRR2; L1 = mLL2; u1 = muu2; i1 = mii2; t1 = mtt2, |
(3.80) |
где m1, ..., mj, ..., mn или mR, mL, mi и mt – масштабные коэффициенты.
В соответствии с первой теоремой подобия для подоб-
ных процессов 0 и 0, все члены уравнений которых однородные функции, должны существовать одинаковые критерии подобия. Их отыскивают приведением уравнения к безразмерному виду.
К безразмерному виду уравнения приводятся делением их на какой-либо, например на m-й, член ( m и m):
|
|
|
|
|
|
m |
|
|
|
1 |
|
i |
|
|
m 1 |
1 |
i |
0, |
(3.81) |
m |
m |
|
m |
||||||
|
|
m |
i 1 |
|
|
||||
|
|
|
97 |
|
|
|
|
|
Ô |
1 |
|
|
|
|
|
|
|
Ô |
i |
|
|
|
|
Ô |
|
|
|
|
|
|
|
m |
Ô |
i |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
m 1 |
1 |
|
|
|
|
|
0, |
(3.82) |
|||||||||||||||||
Ô m |
|
|
Ô m |
|
|
Ô m |
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
Ô m |
|
|
|
|
i 1 |
|
|
|
|||||||||||||||||
1 |
2 |
|
|
3 |
i1R1 |
|
|
L1di1 |
|
|
|
|
u1 |
|
|
0 , |
(3.83) |
||||||||||||||||||
|
i R dt |
|
i R |
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
i R |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
1 |
|
|
1 |
|
|
|
|
|
1 |
1 |
1 |
|
1 |
1 |
1 |
|
|
|
|
1 |
1 |
|
|
|
|
|
|
||||||
Ô1 |
|
|
Ô 2 |
|
|
Ô3 |
i2 R2 |
|
|
L2di2 |
|
|
|
|
u2 |
|
|
|
0 . |
(3.84) |
|||||||||||||||
|
|
|
i R dt |
|
|
i R |
|
|
|||||||||||||||||||||||||||
Ô |
1 |
|
|
|
Ô |
1 |
|
|
|
|
Ô |
1 |
i R |
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
2 |
2 |
|
|
|
|
|
|
2 |
|
2 |
|
|
|
Вследствие однородности (3.85) и (3.84) в выражениях
для i и i существуют некоторые общие множители, которые можно вынести за знак функции. Общий множитель для
i-го члена i исходного уравнения (3.83) – некоторая комбинация масштабных коэффициентов m1, ..., mj, ..., mn, т.е. в соот-
ветствии с (3.85) и (3.89):
i = fi (P1, ..., Pj, ..., Pn) =
|
|
|
|
|
= fi (m1R1, ..., mjRj, …, mnRn) = |
|
(3.85) |
||||||||||||||||||||||||
|
|
|
|
= Mifi (R1, ..., Rj, ..., Rn) = Mi . |
|
|
|||||||||||||||||||||||||
i R |
m i m R m m i R |
M Ô |
, |
|
|
||||||||||||||||||||||||||
|
|
1 1 |
|
|
|
|
i 2 R |
2 |
|
|
|
i R 2 2 |
1 1 |
|
|
|
|||||||||||||||
|
|
|
di |
|
|
|
|
dm i |
|
|
|
|
m m |
|
di |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
L1 |
|
|
1 |
mL L2 |
|
|
i |
2 |
|
|
L |
i L2 |
|
|
2 |
|
M |
2Ô2 , |
(3.86) |
||||||||||||
|
|
|
|
|
|
|
|
dt |
|
||||||||||||||||||||||
|
|
|
dt |
|
|
|
|
dm t |
2 |
|
|
|
m |
|
|
|
2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
1 |
|
|
|
|
|
|
t |
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
||||
u1 muu2 M3Ô3. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
Подстановка (3.85) в (3.81) дает |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
M1 |
|
Ô1 |
|
|
|
Mi |
|
Ôi |
|
Mm 1 Ôm 1 1 0 , |
(3.87) |
|||||||||||||||||||
|
Mm Ôm |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
Mm Ôm |
|
|
|
Mm Ôm |
|
|
|||||||||||||||||||
|
|
|
1 |
|
|
mLmi |
|
|
L2 |
|
di2 |
|
mu |
|
|
u2 |
|
0. |
|
(3.88) |
|||||||||||
|
|
|
|
m m m i R |
|
m m |
i R |
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
dt |
2 |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
i R t |
2 2 |
|
|
|
|
|
i R |
2 2 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
98 |
|
|
|
|
|
|
|
|
|
|
|
|
Поскольку уравнение (3.87) представляет собой сумму однородных функций, должен существовать общий для всех его членов множитель Mm:
|
M1 = ... = Mi = ... = Mm–1 = Mm |
(3.89) |
||||||||
или |
|
|
|
|
|
|
|
|
||
|
M1 |
|
Mi |
|
Mm 1 |
Mm 1 . |
(3.90) |
|||
|
|
Mm |
||||||||
|
Mm |
|
|
Mm |
Mm |
|
||||
Это означает, что |
|
|
|
|
|
|
|
|
||
1 |
|
|
mL |
|
mu |
. |
|
(3.91) |
||
|
|
|
|
|||||||
|
|
|
|
m m |
m m |
|
|
|||
|
|
|
|
|
R t |
i R |
|
|
Физический смысл результата: как исходное уравнение (3.76), так и преобразованное (3.88) описывают процесс0, т.е. переходный процесс i1(t1). Подобие процессов 0 и 0 означает, что они должны описываться одинаковыми уравнениями, так как имеют одинаковый качественный характер и различаются только масштабами.
Таким образом, формула (3.88) будет описывать 0 – переходный процесс i2(t2) – только в том случае, если будут равны единице комбинации масштабных коэффициентов при втором и третьем членах (3.88).
Тождественность (3.87) и (3.81) означает, что между соответственными членами (3.75) и (3.74) существуют отношения:
1 |
|
Ô1 |
,..., |
i |
|
Ôi |
,..., |
m 1 |
|
Ô m 1 . |
(3.92) |
|
m |
Ô m |
m |
Ô m |
|||||||||
|
|
|
|
m |
|
Ô m |
|
Обобщая (3.92) на произвольное число s подобных процессов, уравнения которых содержат только однородные функции, можно записать:
99
(1) |
|
(2) |
|
(s) |
|
|
|
1 |
|
1 |
|
1 |
1 idem, |
|
|
(1)m |
|
(2)m |
|
(ms) |
|
|
|
.................................................... |
|
|
|||||
|
|
||||||
(1) |
|
(2) |
|
(s) |
|
|
|
|
|
|
|
|
|||
i |
|
i |
|
i |
i idem, |
|
(3.93) |
(1)m |
|
(2)m |
|
(ms) |
|
|
|
.................................................... |
|
|
|||||
|
|
||||||
(1) |
|
(2) |
|
(s) |
|
||
|
|
|
|
||||
m 1 m 1 m 1 m 1 idem, |
|
||||||
(1) |
|
(2) |
(s) |
|
|
||
m |
|
m |
|
m |
|
||
|
|
|
|
где (1), (2), (i), ..., (s) – номера сопоставляемых процессов; idem означает «соответственно одинаково для всех рассматриваемых процессов».
Для рассматриваемого примера
2 |
|
L1 |
di1 |
|
Ô 2 |
|
L2 |
|
di2 |
idem, |
||||
i R |
i R |
|||||||||||||
|
|
dt |
|
Ô |
1 |
|
dt |
2 |
|
|||||
1 |
|
1 |
1 |
1 |
|
|
|
2 |
2 |
|
|
3 |
|
u1 |
|
|
Ô3 |
|
u2 |
|
idem. |
||
i R |
i R |
||||||||||
|
|
|
Ô |
1 |
|
|
|||||
1 |
|
1 |
1 |
|
|
|
2 |
2 |
|
В общем случае соотношения пропорциональности вида (3.79) справедливы на любых (и малых, и больших) интервалах изменения сопоставляемых функций. Поэтому символы дифференцирования и интегрирования при рассмотрении условий пропорциональности можно опустить, так как они не имеют размерности и не влияют на условия пропорцио-
нальности, заменив соответствующие члены уравнений j на их аналоги j*, которые называются интегральными, т.е. за-
менить d nx/dyn на x/yn и xdy на xy. Для рассмотренного примера
* |
= L1 |
i |
|
* |
= i1R1; |
|
*2 |
|
L1i1 |
|
L1 |
|
|
2 |
1 |
; |
1 |
1 = |
1* |
= |
|
= |
|
, |
|||
i1R1t1 |
R1t1 |
||||||||||||
t1 |
|||||||||||||
|
|
|
|
|
|
|
|
|
где 2* и 1* – аналоги 2 и 1. 100