Научная статья на тему 'МЕТОДИКА ИССЛЕДОВАНИЯ НЕСИММЕТРИЧНЫХ РЕЖИМОВ СИНХРОННЫХ МАШИН НА ОСНОВЕ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ВОЛЬТЕРРА ВТОРОГО РОДА'

МЕТОДИКА ИССЛЕДОВАНИЯ НЕСИММЕТРИЧНЫХ РЕЖИМОВ СИНХРОННЫХ МАШИН НА ОСНОВЕ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ВОЛЬТЕРРА ВТОРОГО РОДА Текст научной статьи по специальности «Математика»

CC BY
49
9
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СИНХРОННЫЕ МАШИНЫ / ПЕРЕХОДНЫЕ ПРОЦЕССЫ / НЕСИММЕТРИЧНЫЕ РЕЖИМЫ / УРАВНЕНИЕ ВОЛЬТЕРРА ВТОРОГО РОДА / SYNCHRONOUS MACHINES / TRANSIENTS / ASYMMETRIC MODES / VOLTERRA EQUATION OF THE SECOND KIND

Аннотация научной статьи по математике, автор научной работы — Данеев А.В., Данеев Р.А., Сизых В.Н.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A METHODOLOGY FOR RESEARCHING ASYMMETRIC MODES OF SYNCHRONOUS MACHINES BASED ON THE VOLTERRA INTEGRAL EQUATIONS OF THE SECOND KIND

In industry, the solution to many problems is directly related to the development and research of synchronous machines, which operate into a rectifier (asymmetric) load. Transients, which are described by a system of nonlinear differential equations in synchronous machines, play a fundamental role in it. In the study of synchronous machines, it suffices to consider electromagnetic transients due to the large inertial constant of the machine. The equations become linear, but with periodic coefficients, which also do not have a common solution, since they contain periodic coefficients. The possibilities for the practical application of such equations are generally limited by the difficulties associated with the determination of eigenvalues. In this work, we apply the transformation of differential equations with periodic coefficients to equations with constant coefficients based on the representation of a system with periodic coefficients by the Volterra matrix integral equation of the second kind. The consideration is carried out by the example of a three-phase magnetoelectric generator operating on an active-inductive load. Based on a comparison with the classical theory of systems with periodic coefficients, the article established that the constant matrix B is not strictly defined in the considered method. The proposed modeling method allows us to study both symmetric and asymmetric transients in synchronous machines. The method does not have a limitation associated with the sinusoidal spatial distribution of the magnetic driving force of the synchronous machine windings and can be applied taking into account the higher harmonics of the inductances of synchronous machines.

Текст научной работы на тему «МЕТОДИКА ИССЛЕДОВАНИЯ НЕСИММЕТРИЧНЫХ РЕЖИМОВ СИНХРОННЫХ МАШИН НА ОСНОВЕ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ ВОЛЬТЕРРА ВТОРОГО РОДА»

10. Mesarovi M., Mako D., Takahara Y., Theory of Hierarchical Multilevel Systems. Academic Press, New York, 1970, 294 p.

11. Saaty T. L. The Analytic Hierarchy Process: Planning, Priority Setting, Resource Allocation. New York, McGraw-Hill, 1980, 287 p.

12. Stanley R. P. Enumerative Combinatorics, Vol. 1. Cambridge University Press, Cambridge, 1997, 335 p.

13. Birkhoff G., Lattice Theory. American Mathematical Society, New York, 1967, 418 p.

14. Balagura A.A., Kuz'min O. V. Obobshchennaya piramida Paskalya i chastichno uporyadochennye mnozhestva [Generalized Pascal pyramid and partially ordered set]. Obozrenie prikladnoi i promyshlennoi matematiki [Surveys on Applied and Industrial Mathematics], Vol. 14, No. 1, pp. 88-91 (2007).

15. Kuz'min O. V. Obobshchennye piramidy Paskalya i ikh prilozheniya [Generalized Pascal Pyramids and their Applications]. Nauka Publ., Novosibirsk, 2000. 294 p.

16. Kuz'min O. V., Seregina M. V. Ploskie secheniya obobshchennoi piramidy Paskalya i ikh interpretatsii [Plane sections of the generalised Pascal pyramid and their interpretations]. [Discrete Mathematics and Application], 2010, Vol. 20, No. 4, pp. 377-389.

17. Kuz'min O. V., Khomenko A. P., Artyunin A. I. Discrete model of static loads distribution management on lattice structures. Advances and Applications in Discrete Mathematics, 2018, Vol. 19, Iss. 3, pp. 183-193.

18. Kuz'min O. V., Khomenko A. P., Artyunin A. I. Development of special mathematical software using combinatorial numbers and lattice structure analysis. Advances and Applications in Discrete Mathematics, 2018, Vol. 19. Iss. 3, pp. 229-242.

19. Kuz'min O. V., Starkov B. A. Ierarkhicheskie struktury tipa treugol'nika Paskalya i postroenie navigatsionnykh marsh-rutov [Hierarchical structures of the Pascal triangle type and the construction of navigation routes]. [Current problems of science in the Baikal region]. Irkutsk State University Publ., Irkutsk, 2020, Vol. 3, pp. 119-123.

20. Lavlinskaya A. A., Fil' G. A., Kamnev M. D. Sozdanie modeli kvadrokoptera-ekologa [Creation of a quadcopter-ecolog model]. Prikladnye voprosy diskretnogo analiza: sb. nauch. tr. [Applied issues of discrete analysis: collection of research papers]. Irkutsk State University Publ., Irkutsk, 2020, Vol. 5, pp. 78-83

Информация об авторах

Кузьмин Олег Викторович - д. ф.-м. н., профессор, заведующий кафедрой теории вероятностей и дискретной математики, Иркутский государственный университет, г. Иркутск, e-mail: quzminov@mail.ru

Лавлинский Максим Викторович - соискатель, кафедра теории вероятностей и дискретной математики, Иркутский государственный университет, г. Иркутск, e-mail: LavlinskiMV@mail.ru

DOI 10.26731/1813-9108.2020.2(66). 143-150

Information about the authors

Oleg V. Kuz'min - Dr. Sci. in Physics and Mathematics, Prof., Head of the Subdepartment of Probability Theory and Discrete Mathematics, Irkutsk State University, Irkutsk, e-mail: quzmi-nov@mail.ru

Maksim V. Lavlinskii - external Ph.D. student, the Subdepartment of Probability Theory and Discrete Mathematics, Irkutsk State University, Irkutsk, e-mail: LavlinskiMV@mail.ru

УДК 62-97/98

Методика исследования несимметричных режимов синхронных машин на основе интегральных уравнений Вольтерра второго рода

А. В. Данеев 1, Р. А. Данеев 2, В. Н. Сизых 1М

1Иркутский государственный университет, г. Иркутск, Российская Федерация 2Восточно-Сибирский институт МВД России, г. Иркутск, Российская Федерация И sizykh_vn@mail.ru

Резюме

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

Ключевые слова

синхронные машины, переходные процессы, несимметричные режимы, уравнение Вольтерра второго рода

Данеев А.В. Методика исследования несимметричных режимов синхронных машин на основе интегральных уравнений Вольтерра второго рода / А.В. Данеев, Р.А. Данеев, В.Н. Сизых // Современные технологии. Системный анализ. Моделирование. - 2020. - Т. 66 № 2. - С. 143-150. - DOI: 10.26731/1813-9108.2020.2(66).143-150

Информация о статье

поступила в редакцию: 20.02.2020, поступила после рецензирования: 14.03.2020, принята к публикации: 22.03.2020

A methodology for researching asymmetric modes of synchronous machines based on the Volterra integral equations of the second kind

A. V. Daneev \ R. A. Daneev 2, V. N. Sizykh

1 Irkutsk State Transport University, Irkutsk, the Russian Federation

2 East Siberian Institute of the Ministry of Internal Affairs ofRussia, Irkutsk, the Russian Federation И sizykh_vn@mail.ru

Abstract

In industry, the solution to many problems is directly related to the development and research of synchronous machines, which operate into a rectifier (asymmetric) load. Transients, which are described by a system of nonlinear differential equations in synchronous machines, play a fundamental role in it. In the study of synchronous machines, it suffices to consider electromagnetic transients due to the large inertial constant of the machine. The equations become linear, but with periodic coefficients, which also do not have a common solution, since they contain periodic coefficients. The possibilities for the practical application of such equations are generally limited by the difficulties associated with the determination of eigenvalues. In this work, we apply the transformation of differential equations with periodic coefficients to equations with constant coefficients based on the repre-sentation of a system with periodic coefficients by the Volterra matrix integral equation of the second kind. The consideration is carried out by the example of a three-phase magnetoelectric generator operating on an active-inductive load. Based on a compari-son with the classical theory of systems with periodic coefficients, the article established that the constant matrix B is not strictly defined in the considered method. The proposed modeling method allows us to study both symmetric and asymmetric transients in synchronous machines. The method does not have a limitation associated with the sinusoidal spatial distribution of the magnet-ic driving force of the synchronous machine windings and can be applied taking into account the higher harmonics of the induct-ances of synchronous machines.

Keywords

synchronous machines, transients, asymmetric modes, the Volterra equation of the second kind

For citation

Daneev A.V., Daneev R.A., Sizykh V.N. Metodika issledovaniya nesimmetrichnykh rezhimov sinkhronnykh mashin na osnove integral'nykh uravnenii Vol'terra vtorogo roda [A methodology for researching asymmetric modes of synchronous machines based on the Volterra integral equations of the second kind]. Sovremennye tekhnologii. Sistemnyi analiz. Modelirovanie [Modern Technologies. System Analysis. Modeling], 2020, Vol. 66, No. 2, pp. 143-150. DOI: 10.26731/1813-9108.2020.2(66).143-150

Для цитирования

Article info

Received: 20.02.2020, Revised: 14.03.2020, Accepted: 22.03.2020

Введение <

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

В общем случае ПП в СМ описываются системой м нелинейных дифференциальных уравнений и не ч имеют общего решения: ю = var. Однако в большин- ] стве случаев исследования СМ достаточно рассмот- | реть электромагнитные ПП (ю = const) в силу боль- ] шой инерционной постоянной машины. Уравнения становятся линейными, но с периодическими коэффициентами, которые, однако, также не имеют об- ' щего решения, поскольку содержат периодические 1 коэффициенты. Возможности практического применения таких уравнений ограничиваются в общем д

случае трудностями, связанными с определением собственных чисел.

Для некоторых нестационарных систем оказывается возможным переход с помощью преобразования координат от уравнений с переменными коэффициентами к уравнениям с постоянными коэффициентами без решения задачи определения собственных чисел. К таким системам относятся СМ. Удобные для практики результаты получены в работе А.А Горева [1] на основе линейного ортогонального преобразования к вращающимся осям d, q, 0. В то же время такое преобразование имеет ряд ограничений:

- не позволяет избавиться от периодических коэффициентов в случае общей несимметрии фаз и при работе СМ на несимметричную нагрузку;

- СМ рассматривается как источник синусоидальной ЭДС без учета высших гармонических

составляющих.

Поэтому при исследовании несимметричных режимов СМ предпочтительно использовать систему в фазовых координатах. Однако, как показывают расчеты, обратная матрица индуктивности при приведении системы к нормальной форме оказывается слабо обусловленной, что ведет при моделировании на ЭВМ в физических единицах к потере точности вычислений и к большим затратам машинного времени.

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

К числу таких преобразователей относится преобразование, выполненное на основе представления системы с периодическими коэффициентами матричным интегральным уравнением Вольтерра второго рода [2-4].

Постановка задачи исследования многофазных синхронных машин

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

Влиянием демпферных контуров пренебрегаем. Остальные допущения те же, что и в работе А.Ш. Важнова [5].

Обозначим индексами а, Ь, с, т обмотки фаз а, Ь, с и «фиктивную» обмотку постоянного магнита т.

Тогда уравнения баланса напряжений СМ запишутся в виде

d у dt

= - Ri + u,

где

г> • т di

u = Rh i + Ln—, dt

(1)

(2)

U =k , Ub , Uc, Um Ï -

ца активных сопротивлений нагрузки.

Для линейной системы потокосцепления можно записать

(3)

где

L =

У = Li,

'4 Mab Mac Mam

Mab Lb Mbc Mbm

Mca Mcb Lc Mcm

Mma Mmb Mmc L m

матрица

индукгивностеи и взаимных индуктивностеи соответствующих обмоток -

La = ¡0 + Ls + ¡2 C0S 2У' Lb = l0 + Ls +

+/2 cos(2y +120), L = h + L + h cos(2y-120), Mba = Mab = ^0 + Ms + ¡2 cos (2y + 120) ,

Mb = Mbc = m0 + Ms + ¡2 cos 2y,

Mam = Mma = Mmd cos Y, Mmb = Mbm =

= Mmd cos(y- 120),

Mmc = Mm = Mmd cos (y + 120) ,

l = co^, m0 = -M = -—,

2' * 2

смысл коэффициентов l0, l2, L, Mmd тот же [5].

Уравнения (1), разрешенные относительно фазных токов, c учетом (2) и (3) имеют вид di

= A(t ) i, (4)

dt

где

dL

вектор мгновенных значений напряжений; у = , , , 1[ - вектор потокосцепления обмоток СМ;

Я = Я, Я, 0) - матрица активных сопро-

тивлений фазных обмоток; 1 = [^, ^, 1с, /т ] - вектор мгновенных значений токов; = diag (, , , 0) - матрица индуктивно-

стей нагрузки; Ян = diag (Ян, Ян, Яд, 0 ) - матри-

A(t) = (L - LH )-1 " R " dLJ = A(t + T) (5)

- периодическая матрица параметров СМ и нагрузки c периодом T=2n и размерности 4*4. Решение (4) имеет вид

i = Ф(/ ,т)0,

где ф(/,т) - фундаментальная матрица, имеющая ту же размерность, что и матрица А, она является решением матричной системы

Ф( t, т) = A(t ) Ф( t, т), Ф(т, т) = e;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где i0 - вектор начальных токов в момент времени т. Обозначим

Ф(г) = Ф(г, о) = Ф(г,т\_о, Ф(о, о) = Ф(о) = e,

где Е - единичная матрица размерности 4*4. Тогда

ФФ (t ) = А(ОФ(0. (6)

При выполнении условия (5) фундаментальная матрица системы (4) удовлетворяет тождеству

Ф(г+T )=Ф(г )ф(т ), (7)

где ф(Г ) = ф(г ) - матрица монодромии (матри-

цант системы (4)).

Справедливость тождества (7) проверяется подстановкой (7) в (6) при выполнении условия (5).

Аналитическое решение интегральных уравнений Вольтерра второго рода на основе метода последовательных приближений

Основное положение классической теории систем с периодическими коэффициентами базируется на теореме Флоке - Ляпунова: для уравнения (4) с периодической непрерывной периода Т матрицей А существует непрерывная, неособая периодическая периода Т, имеющая кусочно-непрерывно интегрируемую производную матрицу преобразования

¥(г ) = Ф(г )ехр (- В г), (8)

такая, что замена

i = F ( t ) y

(9)

приводит матричное уравнение (4) к уравнению с постоянными коэффициентами

dy

— = By. dt

(10)

В формуле (11) заменим t на t+T: ф(— + T ) = V (t + T ) exp ( Bt + BT ) =

= Ф( t ) exp ( BT ). Из сравнения (12) и (7) имеем

B = lin Ф^ ).

(12)

(13)

Следует отметить, что постоянная матрица не является строго определенной, так как значение in Ф^ ) многозначно [6].

Покажем, что последнее утверждение справедливо для метода [2, 3].

Продифференцируем уравнение (8) с учетом (11) и (13)

V(t) = A(t)v(t) - V(t ) exp (B t) B exp (- B t). (14)

С другой стороны подстановка (9) и (10) в (4) дает

V (t ) = A(t )V (t )- V (t )B. (15)

Нетрудно заметить, что формула (15) является частным случаем (14) при

exp(Bt ) = E, exp(-Bt) = E.

Разложим матричные экспоненты (экспоненциалы) в ряд Тейлора и подставим в выражение (14):

B

2

Bn

exp(Bt ) = E + Bt + — t2 +... + — t

2!

n!

B2

Bn

exp(- Bt ) = E - Bt +—t2 +... +(- 1)П— tn.

2!

n!

Тогда

B2

Г (г) = А(г У (г)- V (г )|в + (-1)" ^ г2п | (16)

Уравнение (16) эквивалентно интегральному матричному уравнению Вольтерра второго рода

V ( t ) = V ( o)+j{ A (e).F (e)

+

+V (e)

B + (-!)'

B

2n+1

\2n

( n!)2

>d e, (17)

где V(0) = ф(o)exp(- B • 0) = E, n = 1, 2, ...

Сделаем замену 9 = вф ( y = — )и перепи-

s

шем (17) в виде

Теорема сформулирована в работе «Линейные дифференциальные уравнения с периодическими коэффициентами и их приложения» [6], а первый пример ее практической реализации в приведен в статье В.М. Лупкина [7].

Согласно этой теореме фундаментальная матрица может быть выражена через матрицу преобразования

ф(г ) = V (г )ехр (вг). (11)

у

V (у, s) = V ( 0, s) + sJ{ A (ф)-Г (ф, s)

Б1-

-V (ф, s)

B(s) + (-l)n s2пЩ ^n (n!)

2n 2 ф2

^ф, (18)

где V (0, в) = Е, в =--«малый» параметр.

ю

Заметим, что в уравнении (18) являются неизвестными матрица преобразования V(у,в) и постоянная матрица В (в) .Согласно [3] решение (18)

может быть получено в виде абсолютно сходящегося относительно параметра в ряда методом последовательных приближений Пикара:

V (у, s) = V ( 0, s) + ^skVk (у)

(19)

k=1

где

i

Vk (y) = J{A(ф). Vk-1

-V

k-1

\2n+1

B (s) + (-1)n s2 ф2 n

( n!)

dф ,(20)

к = 1,2,... п = 1,2,.... V = Е(при к = 1).

В [8] показана сходимость ряда (19) через существование постоянной матрицы В, записанной аналитически относительно параметра в

n

0

0

0

B(e) = £ek-% .

(21)

k=1

Подставим (21) в матричное уравнение (20) и сравним в нем коэффициенты при одинаковых степенях к в исходном (19) и преобразованном с учетом (20) и (21) уравнениях

ад ад у

(у)^{{А(фК- -

k=1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-V,

k-1

k=1 0

f да

ч 2n+1

Zek-1Bk + (-i)n e2n

Z Bk

\k=1 y

k=1

( n!)2

dф .(22)

Bk+1 =

B1 =1 j А(Ф^ Ф,

T 0

y

V (y) = j(А(ф) - B )dф,

0

1 T

- j{ А(ф)^+1 (ф)-Vk+1 (ф) B1 -...-V (ф) Bk+1$ dф,

I

Vk+1 (y) = j{A(ф)Vk(ф)-Vk (ф)B1 -...-V(ф)Bk -Bk+ }dф.

2. Решается система уравнений с постоянными коэффициентами

£ = В (8) у

аг

с начальными условиями у0 = V(0, б) /0.

3. Определяются переменные состояния системы (4) в фазных координатах

г = V(у, е)у .

Линейное приближение экспоненциалов (п = 1)

Тождество (22) принимает вид

ад ад у

(У)^{{А(ф).V- -

В тождестве (22) при п = 1, 2 определяем число членов разложения экспоненциалов в ряд Тейлора в формуле (16) без учета нулевой составляющей этого разложения.

Используя условие периодичности системы (4)

¥к(Т) = ¥к(0) = 1 I 1.2.....

получим формулы, определяющие матричные коэффициенты в разложениях рядов (20) и (21).

При определении матриц V(у, в) и

В (в) ограничимся двумя важными для практики случаями [4].

Нулевое приближение экспоненциалов

Случай нулевого приближения экспоненциалов

ехр (В г) = е, ехр (-В г) = Е

рассмотрен в работе [3]. Действительно, из формулы

(8) следует: V(г) = Ф^), т. е. предложенный метод [3] является по существу алгоритмом определения фундаментальной (переходной) матрицы системы (4). Алгоритм вычислений сводится к следующему:

1. Матрица приведения V (У, 8) и постоянная матрица В(в) представляются в виде абсолютно сходящихся рядов по формулам (20) и (21), в которых

k=1

k=1 0

-Vk-,

Ze*-1Bk -0,25e21 £e*-1Bk

к=1

к=1

ф

>d ф.

Для первого приближения метода ( к= 1) имеем

у

V (у) = |{ А(ф) - В + 0,25в2В3ф2}аф.

0

Из условия периодичности матрицы V (у)

т

V (т) = |{А(ф) - В + 0,25в2В3ф2} аф = 0

определим матрицу В^:

2 2 I T

B -— B13 = - jA(ф)dф.

л

T

0

Таким образом, матрица В (в) не является строго определенной, что согласуется с результатами, полученным из классической теории (формула (13)). Однако ввиду малости параметра в можно считать

1 T

B1 =1 j A(ф)d ф.

Для последующих приближений алгоритма можно поступать следующим образом: матрицу В1 определять приближенно, используя условие малости параметра в, а матрицу преобразования V(г) вычислять с учетом величин, вносимых в алгоритм вычислений от последующих членов разложения экспоненциалов в ряд Тейлора.

Недостатком метода является ограничение в области низких частот. С уменьшением частоты ю параметр в увеличивается, сходимость рядов (20) и (21) ухудшается и, следовательно, возникает необходимость рассмотрения большего числа членов ряда.

0

0

Пример расчета переходных процессов синхронных машин на основе решения интегральных уравнений Вольтерра второго рода

Применим к системе (4) изложенный метод. При этом ограничимся первым приближением алгоритма

V (у, s) = E + s V (У),

где

у

v, (у) = {{ а(ф) - в,} d ф,

0

b(s) = b,, 1 T

b, = - j a(ф)d ф, T = 2 л. T 0

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

А *0 Q 0

а0 0

Q а0 0

0 0 0 1

* («0 =

(y, 6) v|2 (y, 6) 6) ^fa б)

v31 (y, 6) v32 (y, 6)

/41^ 6) v42^ 6)

где элементы матрицы преобразования определяются выражениями:

V(y,6) =

uy, 6) v4(y, 6)"

6) v24(y, 6)

vö(y, 6) v34 (y, 6)

v43(y, 6) v44(y, 6)

v! (y, 6) = 1 + — (A + a sin 2y - A cos 2y),

^ 6) = 2(-a+7+

+A sin(2y -120) - A cos(2y -120)),

^s(Y, 6) = "

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

-ДА A

+a sin(2Y +120) - a cos(2Y +120));

V22(y, 6) = 1 + -

A_ V3A2 2 2

+A sin(2Y +120) - A cos(2Y +120) ) , 6

V23 (y, 6) = — (A + A sin2y- A cos2y),

V33(y, 6) = 1 + -

A y^A

2

+A sin(2y -120) - A cos(2y -120)), vm( y, 6) = ^(y, 6),v31(y, 6) = ^(y, 6),

V23(y 6) = V32(y, 6),

V14(Y, 6) = -Де(1 + cos y), v2 4 (y, 6) = - 6(1 + cos(Y-120)), V34 (y, 6) = -D16 (1 + cos( y +120) ), v41 (y, 6) = 6 (d + d3 sin y - D cos y),

V42^ 6) = 6

V3

D — D + 2 3 2 2

+D3 sin(Y -120) - D2 cos(Y -120)),

(

V43(Y, 6) = 6

Л D +1D + 2 3 2 2

+Д яп(у +120) - Д С08(у +120)), Г44 (у, в) = 1.

Коэффициенты А0, В0, С0, Аь А2, А, В2, Б3 зависят от параметров СМ и нагрузки [9].

Ряд близких и смежных вопросов моделирования объектов такой физической природы рассмотрен в работах [10-19].

Заключение

Таким образом, на основе сопоставления с классической теорией систем с периодическими коэффициентами установлено, что постоянная матрица В в методе, изложенном в работе «Линейные дифференциальные уравнения с периодическими коэффициентами и их приложения» [6], не является строго определенной. Предложенная методика моделирования позволяет:

1. Исследовать как симметричные, так и несимметричные переходные режимы в СМ.

2. Снять ограничения, связанные с синусоидальным пространственным распределением МДС обмоток при учете высших гармонических составляющих в матрице индуктивностей и взаимных индук-тивностей СМ.

Список литературы

1. Горев А.А. Переходные процессы синхронной машины. М.: СЭИ, 1950. 551 с.

2. Еругин Н.П. Приводимые системы // Труды МИАН им. М.А. Стеклова, 1946. Т. 12. С. 3-96.

3. Бреус К.А. Об одном классе линейных дифференциальных уравнений с периодическими коэффициентами // Укр. матем. журн. 1960, Т. 12. № 4. С. 25-32.

2

2

2

4. Лупкин В.М. Обобщение методов приведений и аналитического решения уравнений несимметричных электрических машин // Электричество, 1985. № 2. С. 22-29.

5. Важнов А.Ш. Переходные процессы в машинах переменного тока. Л.; Энергия, 1980. 320 с.

6. Якубович В.А., Старжинский В.М. Линейные дифференциальные уравнения с периодическими коэффициентами и их приложения. М.; Наука, 1972. 720 с.

7. Лупкин В.М. Аналитическое решение линейных дифференциальных уравнений вентильного двигателя // Электричество. 1981. № 6. C. 22-31.

8. Бреус К.А. О приводимости канонической системы дифференциальных уравнений с периодическими коэффициентами // ДАН СССР, 1958. Т. 123. № 1. С. 21-25.

9. Александров А.А., Данеев Р.А., Сизых В.Н. К вопросу моделирования вентильных синхронных машин на основе квазианалитического метода // Известия Самарского научного центра РАН, 2019. Т. 21. № 4. С. 63-69.

10. Данеев А.В., Русанов В.А. Об одном классе сильных дифференциальных моделей над счетным множеством динамических процессов конечного характера // Известия высших учебных заведений. Математика. 2000. № 2. С. 32-40.

11. Rusanov V.A., Antonova L.V., Daneev A.V. Inverse problem of nonlinear systems analysis: a behavioral approach // Advances in Differential Equations and Control Processes, 2012. Т. 10. № 2. Pp. 69-88.

12. Данеев А.В., Лакеев А.В., Русанов В.А., Русанов М.В. К теории реализации сильных дифференциальных моделей. I // Сибирский журнал индустриальной математики, 2005. Т. 8. № 1 (21). С. 53-63.

13. Данеев А.В., Русанов В.А., Шарпинский Д.Ю. Нестационарная реализация Калмана-Месаровича в конструкциях оператора Релея-Ритца // Кибернетика и системный анализ, 2007. № 1. С. 82-91.

14. Сизых В.Н., Мухопад А.Ю. Ассоциативный автомат адаптивного управления технологическими процессами на основе нейронных сетей // Научный вестник Новосибирского государственного технического университета, 2014. № 1 (54). С. 34-45.

15. Сизых В.Н. Итерационно-релаксационный метод нелинейного синтеза регуляторов // Автоматика и телемеханика, 2005. № 6. С. 47-58.

16. Мухопад Ю.Ф., Пашков Н.Н., Сизых В.Н. Адаптивный подход к нейронному управлению одним классом абсолютно устойчивых систем // Фундаментальные исследования, 2011. № 8-1. С. 139-147.

17. Данеев А.В., Русанов В.А., Русанов М.В., Сизых В.Н. К апостериорному моделированию нестационарных гиперболических систем // Известия Самарского научного центра РАН, 2018. Т. 20. № 1 (81). С. 106-113.

18. Сизых В.Н. Итерационно-релаксационный метод приближенно-оптимального синтеза регуляторов // Доклады Академии наук, 2000. Т. 371. № 5. С. 574-576.

19. Агеев А.М., Сизых В.Н. Синтез оптимальных регуляторов системы управления самолетом через решение обратной задачи АКОР // Научный вестник Новосибирского государственного технического университета, 2014. № 3 (56). С. 7-22.

References

1. Gorev A.A. Perekhodnye protsessy sinkhronnoi mashiny [Transients of a synchronous machine]. Moscow: SEI Publ., 1950,551 p.

2. Erugin N.P. Privodimye sistemy [Reducible systems]. [Proceedings of Steklov Mathematical Institute of the RAS], 1946, Vol. 12, pp. 3-96.

3. Breus K.A. Ob odnom klasse lineinykh differentsial'nykh uravnenii s periodicheskimi koeffitsientami [On a class of linear differential equations with periodic coefficients]. Ukr. matem. zhurn. [Ukrain. math. journal], 1960, Vol. 12, No. 4, pp. 25-32.

4. Lupkin V.M. Obobshchenie metodov privedenii i analiticheskogo resheniya uravnenii nesimmetrichnykh elektricheskikh mashin [A generalization of the reduction methods and the analytical solution of the equations of asymmetric electric machines]. Elektrichestvo [Electricity], 1985, No. 2, pp. 22-29.

5. Vazhnov A.Sh. Perekhodnye protsessy v mashinakh peremennogo toka [Transients in AC machines]. Leningrad: Energiya Publ., 1980, 320 p.

6. Yakubovich V.A., Starzhinskii V.M. Lineinye differentsial'nye uravneniya s periodicheskimi koeffitsientami i ikh prilozheniya [Linear differential equations with periodic coefficients and their applications]. Moscow, Nauka Publ., 1972, 720 p.

7. Lupkin V.M. Analiticheskoe reshenie lineinykh differentsial'nykh uravnenii ventil'nogo dvigatelya [Analytical solution of linear differential equations of a valve motor]. Elektrichestvo [Electricity], 1981, No. 6, pp. 22-31.

8. Breus K.A. O privodimosti kanonicheskoi sistemy differentsial'nykh uravnenii s periodicheskimi koeffitsientami [On the reduc-ibility of the canonical system of differential equations with periodic coefficients]. Doklady Akademii Nauk SSSR [Papers of the Academy of Sciences of the USSR], 1958, Vol. 123, No. 1, pp. 21-25.

9. Aleksandrov A.A., Daneev R.A., Sizykh V.N. [On the issue of modeling synchronous valve machines based on the quasianalyt-ic method]. [Bulletin of the Samara Scientific Center of the Russian Academy of Sciences], 2019, Vol. 21, No. 4, pp. 63-69.

10. Daneev A.V., Rusanov V.A. Ob odnom klasse sil'nykh differentsial'nykh modelei nad schetnym mnozhestvom dinamicheskikh protsessov konechnogo kharaktera [On a class of strong differential models over a countable set of dynamical processes of a finite nature]. Izvestiya vysshikh uchebnykh zavedenii. Matematika [Bulletin of Higher Education. Mathematics], 2000, No. 2, pp. 32-40.

11. Rusanov V.A., Antonova L.V., Daneev A.V. Inverse problem of nonlinear systems analysis: a behavioral approach. Advances in Differential Equations and Control Processes, 2012, Vol. 10, No. 2, pp. 69-88.

12. Daneev A.V., Lakeev A.V., Rusanov V.A., Rusanov M.V. K teorii realizatsii sil'nykh differentsial'nykh modelei [On the theory of the implementation of strong differential models]. Sibirskii zhurnal industrial'noi matematiki [Siberian Journal of Industrial Mathematics], 2005, Vol. 8, No. 1 (21), pp. 53-63.

13. Daneev A.V., Rusanov V.A., Sharpinskii D.Yu. Nestatsionamaya realizatsiya Kalmana-Mesarovicha v konstruktsiyakh operato-ra Releya-Rittsa [Unsteady Kalman-Mesarovich implementation in the constructions of the Rayleigh-Ritz operator]. Kibernetika i sistemnyi analiz [Cybernetics and System Analysis], 2007, No. 1, pp. 82-91.

14. Sizykh V.N., Mukhopad A.Yu. Assotsiativnyi avtomat adaptivnogo upravleniya tekhnologicheskimi protsessami na osnove nei-ronnykh setei [The associative automaton of adaptive control of technological processes on the basis of neural networks]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta [Scientific Herald of Novosibirsk State Technical University], 2014, No. 1 (54), pp. 34-45.

15. Sizykh V.N. Iteratsionno-relaksatsionnyi metod nelineinogo sinteza regulyatorov [The iterative-relaxation method of nonlinear synthesis of controllers]. Avtomatika i telemekhanika [Automation and Telemechanics], 2005, No. 6, pp. 47-58.

16. Mukhopad Yu.F., Pashkov NN, Sizykh V.N. Adaptivnyi podkhod k neironnomu upravleniyu odnim klassom absolyutno ustoichivykh sistem [An adaptive approach to neural control of one class of absolutely stable systems]. Fundamental'nye issledovaniya [Fundamental Research], 2011, No. 8-1, pp. 139-147.

17. Daneev A.V., Rusanov V.A., Rusanov M.V., Sizykh V.N. K aposteriornomu modelirovaniyu nestatsionarnykh giperbolich-eskikh sistem [On a posteriori modeling of non-stationary hyperbolic systems]. Izvestiya Samarskogo nauchnogo tsentra RAN [Bulletin ofthe Samara Scientific Center ofthe Russian Academy of Sciences], 2018, Vol. 20, No. 1 (81), pp. 106-113.

18. Sizykh V.N. Iteratsionno-relaksatsionnyi metod priblizhenno-optimal'nogo sinteza regulyatorov [The iterative-relaxation method of approximate optimal synthesis of regulators]. Doklady Akademii Nauk [Papers of the Academy of Sciences], 2000, Vol. 371, No. 5, pp. 574-576.

19. Ageev A.M., Sizykh V.N. Sintez optimal'nykh regulyatorov sistemy upravleniya samoletom cherez reshenie obratnoi zadachi AKOR [Synthesis of optimal controllers for an airplane control system by solving the inverse problem of analytical design of optimal controllers]. Nauchnyi vestnik Novosibirskogo gosudarstvennogo tekhnicheskogo universiteta [Scientific Herald of Novosibirsk State Technical University], 2014, No. 3 (56), pp. 7-22.

Информация об авторах Information about the authors

Данеев Алексей Васильевич - д. т. н., профессор кафедры Aleksei V. Daneev - Doctor of Engineering Science, Professor

информационных систем и защиты информации, Иркут- of the Subdepartment of Information Systems and Information

ский государственный университет путей сообщения, г. Protection, Irkutsk State Transport University, Irkutsk, e-mail:

Иркутск, e-mail: daneev@mail.ru daneev@mail.ru

Данеев Роман Алексеевич - к. т. н., доцент кафедры инфор- Roman A. Daneev - Associate Professor of the Subdepart-

мационно-правовых дисциплин, Восточно-Сибирский ин- ment of Information and Legal Disciplines, East Siberian In-

ститут МВД России, г. Иркутск, e-mail: romasun@mail.ru stitute of the Ministry of Internal Affairs of Russia, Irkutsk, email: romasun@mail.ru

Сизых Виктор Николаевич - д. т. н., профессор кафедры Viktor N. Sizykh - Doctor of Engineering Science, Professor

автоматизации производственных процессов, Иркутский of the Subdepartment of Automation of Production Processes,

государственный университет путей сообщения, г. Ир- Irkutsk State Transport University, Irkutsk, e-mail:

кутск, e-mail: sizykh_vn@mail.ru sizykh_vn@mail.ru

Б01 10.26731/1813-9108.2020.2(66).150-157 УДК 656.2

Модель оптимизации работы городских логистических систем с учетом стратегии краткосрочного планирования

О. А. Лебедева И

Ангарский государственный технический университет, г. Ангарск, Российская Федерация И kravhome@mail.ru

Резюме

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

i Надоели баннеры? Вы всегда можете отключить рекламу.