Научная статья на тему 'Анализ электрических цепей средствами компьютерной алгебры в машиноведении'

Анализ электрических цепей средствами компьютерной алгебры в машиноведении Текст научной статьи по специальности «Математика»

CC BY
140
18
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЛЕКТРИЧЕСКАЯ ЦЕПЬ / ГРАФИЧЕСКИЙ РЕДАКТОР / КОМПЬЮТЕРНАЯ АЛГЕБРА / МОДЕЛИРОВАНИЕ / УСТОЙЧИВОСТЬ / ПАРАМЕТРИЧЕСКИЙ АНАЛИЗ / ELECTRICAL CIRCUIT / GRAPHIC EDITOR / COMPUTER ALGEBRA / MODELING / STABILITY / PARAMETRIC ANALYSIS

Аннотация научной статьи по математике, автор научной работы — Банщиков А.В., Ветров А.А., Иртегов В.Д., Титоренко Т.Н.

Электрические цепи уже давно имеют разнообразные приложения в транспорте и других областях инже-нерных знаний. Создана оригинальная модель и на ее основе в интегрированной среде Embarcadero Delphi на языке Object Pascal разработан визуальный редактор графа электрической цепи. Особенность подхода заключается в том, что базовая информация о графе расширена данными о параметрах и функциях схемы цепи для автоматического формирования ее сим-вольного описания. Посредством графического пользовательского интерфейса для конкретной нелинейной электрической цепи сформированы граф и соответствующее ему описание в виде вложенного списка для дальнейшего исследования созданными ранее программными комплексами по моделированию и качественному анализу электрических цепей в символьном виде на РС. Для исследуемой цепи проведен динамический анализ, а именно, построена символьная модель (смешанный потенциал и диф-ференциальные уравнения состояния), найдены некоторые положения равновесия (не зависящие от времени решения уравне-ний состояния) и рассмотрен вопрос об их устойчивости по Ляпунову посредством выписанных уравнений возмущенного дви-жения в первом приближении. Условия асимптотической устойчивости в виде системы неравенств Льенара Шипара выра-жены через коэффициенты характеристического полинома шестой степени. При параметрическом анализе условий устойчи-вости использовались средства языка программирования и функции символьно-численного моделирования системы компью-терной алгебры Mathematica. Представлена графическая интерпретация результатов исследования устойчивости положений равновесия

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

Похожие темы научных работ по математике , автор научной работы — Банщиков А.В., Ветров А.А., Иртегов В.Д., Титоренко Т.Н.

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

ANALYSIS OF ELECTRICAL CIRCUITS BY MEANS OF COMPUTER ALGEBRA IN MACHINE SCIENCE

Electrical circuits have long been diversely applied in transport and other fields of engineering. The paper presents an original model and, based on it, a visual editor of the electrical circuit graph developed in the integrated environment of Embarcadero Delphi in Object Pascal language. The novelty of the approach is that the basic information about the graph is expanded by data on the parameters and functions of the circuit diagram for the automatic generation of its symbolic description. Through graphical user inter-face for a specific nonlinear electrical circuit, a graph and a corresponding description in the form of a nested list are formed for further investi-gation by previously created software complexes for modeling and qualitative analysis of electrical circuits in symbolic form on PC. A dynamic analysis was carried out for the investigated chain: namely, a symbolic model (mixed potential and differential equations of state) was constructed, certain equilibrium positions (time-independent solutions of the state equations) were found and the question of their Lyapunov stability by means of the outlined equations of perturbed motion in the first approximation... Electrical circuits have long been diversely applied in transport and other fields of engineering. The paper presents an original model and, based on it, a visual editor of the electrical circuit graph developed in the integrated environment of Embarcadero Delphi in Object Pascal language. The novelty of the approach is that the basic information about the graph is expanded by data on the parameters and functions of the circuit diagram for the automatic generation of its symbolic description. Through graphical user inter-face for a specific nonlinear electrical circuit, a graph and a corresponding description in the form of a nested list are formed for further investi-gation by previously created software complexes for modeling and qualitative analysis of electrical circuits in symbolic form on PC. A dynamic analysis was carried out for the investigated chain: namely, a symbolic model (mixed potential and differential equations of state) was constructed, certain equilibrium positions (time-independent solutions of the state equations) were found and the question of their Lyapunov stability by means of the outlined equations of perturbed motion in the first approximation. The conditions for the as-ymptotic stability in the form of a system of Lienard-Chipard inequalities are obtained in terms of the coefficients of the characteristic polynomial of sixth degree. During the parametric analysis of stability conditions, the programming language tools and the functions of symbolic-numerical modeling of the computer algebra system “Mathematica” were used. A graphic interpretation of the results of the study on the equilibrium positions stability is presented function show_eabstract() { $('#eabstract1').hide(); $('#eabstract2').show(); $('#eabstract_expand').hide(); } ▼Показать полностью

Текст научной работы на тему «Анализ электрических цепей средствами компьютерной алгебры в машиноведении»

ИРКУТСКИМ государственный университет путей сообщения

Современные технологии. Системный анализ. Моделирование № 4 (60) 2018

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

Милов А. В. Влияние флюса на точность измерений в процессе индукционной пайки алюминиевых волноводных трактов / А. В. Милов, В. С. Тынченко, А. В. Мурыгин // Современные технологии. Системный анализ. Моделирование. - 2018. - Т. 60, № 4. - С. 38-46. - Б01: 10.26731/1813-9108.2018.4(60).38-46

For citation

Milov A. V., Tynchenko V. S., Murygin A. V. Vliyanie flyusa na tochnost' izmerenii v protsesse induktsionnoi paiki alyuminievykh volnovodnykh traktov [Flux influence on the accuracy of measurements during induction soldering of aluminum waveguide ducts]. Sovremennye tekhnologii. Sistemnyi analiz. Modelirovanie [Modern Technologies. System Analysis. Modeling], 2018, Vol. 60, No. 4, pp. 38-46. DOI: 10.26731/1813-9108.2018.4(60).38-46

УДК 519.688:621.3.011.7 DOI: 10.26731/1813-9108.2018.4(60).46-54

А. В. Банщиков, А. А. Ветров, В. Д. Иртегов, Т. Н. Титоренко

Институт динамики систем и теории управления имени В. М. Матросова Сибирского отделения Российской академии наук, г. Иркутск, Российская Федерация Дата поступления: 3 октября 2018 г.

АНАЛИЗ ЭЛЕКТРИЧЕСКИХ ЦЕПЕЙ СРЕДСТВАМИ КОМПЬЮТЕРНОЙ АЛГЕБРЫ В МАШИНОВЕДЕНИИ

Аннотация. Электрические цепи уже давно имеют разнообразные приложения в транспорте и других областях инженерных знаний. Создана оригинальная модель и на ее основе в интегрированной среде Embarcadero Delphi на языке Object Pascal разработан визуальный редактор графа электрической цепи. Особенность подхода заключается в том, что базовая информация о графе расширена данными о параметрах и функциях схемы цепи для автоматического формирования ее символьного описания. Посредством графического пользовательского интерфейса для конкретной нелинейной электрической цепи сформированы граф и соответствующее ему описание в виде вложенного списка для дальнейшего исследования созданными ранее программными комплексами по моделированию и качественному анализу электрических цепей в символьном виде на РС. Для исследуемой цепи проведен динамический анализ, а именно, построена символьная модель (смешанный потенциал и дифференциальные уравнения состояния), найдены некоторые положения равновесия (не зависящие от времени решения уравнений состояния) и рассмотрен вопрос об их устойчивости по Ляпунову посредством выписанных уравнений возмущенного движения в первом приближении. Условия асимптотической устойчивости в виде системы неравенств Льенара - Шипара выражены через коэффициенты характеристического полинома шестой степени. При параметрическом анализе условий устойчивости использовались средства языка программирования и функции символьно-численного моделирования системы компьютерной алгебры Mathematica. Представлена графическая интерпретация результатов исследования устойчивости положений равновесия.

Ключевые слова: электрическая цепь, графический редактор, компьютерная алгебра, моделирование, устойчивость, параметрический анализ._

A. V. Banshchikov, A. A. Vetrov, V. D. Irtegov, T. N. Titorenko

Matrosov Institute for System Dynamics and Control Theory of Siberian Branch of Russian Academy of Sciences, Irkutsk, the Russian Federation Received: October 3, 2018

ANALYSIS OF ELECTRICAL CIRCUITS BY MEANS OF COMPUTER ALGEBRA IN MACHINE SCIENCE

Abstract. Electrical circuits have long been diversely applied in transport and other fields of engineering. The paper presents an original model and, based on it, a visual editor of the electrical circuit graph developed in the integrated environment of Embarcadero Delphi in Object Pascal language. The novelty of the approach is that the basic information about the graph is expanded by data on the parameters and functions of the circuit diagram for the automatic generation of its symbolic description. Through graphical user interface for a specific nonlinear electrical circuit, a graph and a corresponding description in the form of a nested list are formed for further investi-gation by previously created software complexes for modeling and qualitative analysis of electrical circuits in symbolic form on PC. A dynamic analysis was carried out for the investigated chain: namely, a symbolic model (mixed potential and differential equations of state) was constructed, certain equilibrium positions (time-independent solutions of the state equations) were found and the question of their Lyapunov stability by means of the outlined equations ofperturbed motion in the first approximation. The conditions for the asymptotic stability in the form of a system of Lienard-Chipard inequalities are obtained in terms of the coefficients of the characteristic polynomial of sixth degree. During the parametric analysis of stability conditions, the programming language tools and the functions of symbolic-numerical modeling of the computer algebra system "Mathematica " were used. A graphic interpretation of the results of the study on the equilibrium positions stability is presented

Keywords: electrical circuit, graphic editor, computer algebra, modeling, stability, parametric analysis.

46

© А. В. Банщиков, А. А. Ветров, В. Д. Иртегов, Т. Н. Титоренко, 2018

Машиностроение и машиноведение

Modern technologies. System analysis. Modeling, 2018, vol. 60, no. 4

m

Введение

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

В последние годы существенно возросло значение систем компьютерной алгебры (СКА) [1, 2] и методов их использования [3] для научных вычислений в различных областях естественных наук и инженерных знаний. Многолетний опыт авторов показывает, что СКА оказались достаточно эффективными на этапах построения математических моделей и анализа их свойств. Сочетание методов аналитической механики и теории устойчивости движения с алгоритмами СКА (в данном случае системы МмИвтаИеа [4]) позволяет создавать достаточно мощные инструменты для исследования электрических цепей. На примере одной нелинейной электрической цепи в настоящей работе демонстрируется графический редактор [5] и применение комплексов программ [6, 7] для моделирования и качественного анализа электрических цепей в символьном виде на РС.

В качестве электрической рассматривается цепь, в которой произвольным образом соединены элементы сопротивления, катушки индуктивности, конденсаторы, источники тока или напряжения. Электромеханические аналогии [8, 9] позволяют использовать методы аналитической механики для описания электрических цепей и, тем самым, применять методы исследования механических систем к электрическим цепям. Построение функции Лагранжа и уравнений состояния линейной цепи в форме Лагранжа 2-го рода описано в [10].

В качестве нелинейной цепи рассматривается электрическая цепь, принадлежащая классу «полных» в смысле [11]. Дифференциальные уравнения такой цепи имеют вид:

дР(1, и)

4 (4)

dt

di.

du¡ dP(i, u)

Cl (ul>~r = ~~

dt

du.

(k = 1,...,r),

(l = r + 1,...,r + m),

(1)

где Lk (ik) - индуктивности, C¡ (ul) - емкости, ik -токи, протекающие через катушки индуктивности, u¡ - напряжения на конденсаторах. Функция P(i,u) описывает часть цепи, содержащую элементы сопротивления и источники напряжения (тока), и называется смешанным потенциалом. Здесь i = (i1,...,ir), u = (ur+1,...,ur+m) - вектор-строки.

Метод построения смешанного потенциала по графу нелинейной цепи подробно изложен в [10].

Исследование конкретной электрической

цепи

1. Описание редактора для представления графа цепи. Создана модель электрической цепи, особенность которой заключается в расширении классических характеристик неориентированного графа (множество вершин и ребер) за счет параметров цепи. Такими параметрами могут являться постоянные емкости, индуктивности, токи, сопротивления и нелинейные функции. На основе созданной модели в интегрированной среде Embarcadero Delphi на языке Object Pascal разработан графический редактор [5] представления графа электрической цепи, который поддерживает стандартные функции, присущие такого рода системам (см., например, [12]). Создание программного обеспечения происходило в классической парадигме объектно-ориентированного программирования с использованием инкрементной модели разработки. Графический пользовательский интерфейс обеспечивает ввод схемы и данных графа электрической цепи путем интерактивного добавления связанных графических примитивов, предназначенных для построения линейной или специального класса нелинейной электрической цепи. Добавление параметров отдельных элементов цепи к соответствующим вершинам и ребрам графа реализовано через контекстное меню. Отличие от аналогичных разработок [13] заключается в том, что по построенному графу пользователь редактора может автоматически сформировать специальное символьное описание электрической цепи в виде вложенного списка для дальнейшего исследования другими программными комплексами по моделированию [6] и качественному анализу [7] электрических цепей в символьном виде.

Рис. 1. Электрическая схема

На рис. 1 приведена схема рассматриваемой нелинейной электрической цепи [14]. Здесь ¡1, ¡2 -токи, протекающие через катушки индуктивности; и3, и4, и5, ив - напряжения на конденсаторах; С1, С2, С3, С4, Ь1,Ь2, 10,101,102 - параметры цепи, соответствующие постоянным емкостям, ин-дуктивностям, токам. Нелинейными элементами цепи являются сопротивления, вольтамперные характеристики которых описываются функциями /[(ив), /2(и4), /3(и3), /4(и5) . Графики последних предполагаются качественно совпадающими: они находятся в первом и третьем квадрантах координатной плоскости (7,и) , проходят через начало координат и имеют по две экстремальные точки в первом квадранте.

Представим алгоритм использования редактора для создания графа рассматриваемой цепи.

Используя имеющиеся инструменты (например, «эллипс»), создают вершины графа. Следующим шагом является присвоение вершинам семантических данных (например, задать нумерацию).

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

На последнем шаге с помощью инструмента «экспорт символьного описания» автоматически формируется описание графа исследуемой электрической цепи (см. рис. 3) в виде списка пар

{{УьР]з) . Здесь К < У 2 - номера вершин графа, {]1,у2} - ребра графа, рУ - символьная характеристика электрического элемента, соответствующего данному ребру графа.

Машиностроение и машиноведение

Modern technologies. System analysis. Modeling, 2018, vol. 60, no. 4

Модифициро 100% (804, 306}

Рис. 2. Результат построения графа цепи

Рис. 3. Символьное описание графа цепи

ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ПУТЕЙ СООБЩЕНИЯ

Современные технологии. Системный анализ. Моделирование № 4 (60) 2018

2. Построение символьной модели цепи и ее положения равновесий

Исключив с помощью части уравнений (3) переменные и3, и5, ¡1,72 из остальных, получим

На входе пакета расширения ОголШ из [6, два полиномиальных алгебраических уравнения 10] сформированное редактором символьное опи- относительно и и и ■ сание графа исследуемой цепи выглядит следую-

щим образом:

{{{1,2}, Ll},{{2,3], fl},{{23'}, С}, {{3,4}, L2},{{4,5}, 102},{{4,5}, C4},{{4,5}, f4}, {{3,5}, C2},{{3,5}, f2},{{2,5}, I o}, {{1,5}, Сз},{{1,5}, f3 },{{1,5}, loi}} где fj (j = 1,...,4) - обозначения для функций

fj (uk ). Далее с помощью этого пакета строятся смешанный потенциал:

P(i,u) = i1(u4 + u6 -u3) + i2(u5 -u4) +10(u4 + u6) -

I01u3 I02u5

Cu6 Cu4

u5 -J f1(u)du -J f2(u)du

u3 u5

J 3 f3 (u)du - J 5 f4 (u)du

35

0 f3(u)du -J0

и уравнения состояния цепи вида (1):

__ di1 __ di2 L1 u4 + u^ u3 , L2 u5 u

dt

C3 dur = h + f3^)+1

dt

dt

01

C

du

4

2 ,, i2 i1 + f2 (u4) I0

(3)

101 -10 + g1u6(u6 -3s1) + + g3(u4 + u6)(u4 + u6 -3s3)2 = 0, g2u4(u4 -3s2)2 + g4u4(u4 -3s4)2 -- g1u6(u6 -3s1)2 +102 =

(4)

В [15] сделаны оценки на количество действительных решений уравнений (4) (т. е. положений равновесия) и выписаны некоторые из них в громоздком аналитическом виде. В настоящей работе для целей качественного анализа по понятным причинам выбраны достаточно компактные аналитические решения. Решения были получены в [16] при различных ограничениях, налагаемых на параметры задачи. А именно, для ограничений

101 10

102=0:

g 4 g 2

S3 S1 + 33 , S4 S2

имеется решение:

{ i1 I0 , i2 = 4g2S2

u3 3s + S2 , u4 S2 , u5 S2 , u6 3s } ,

(5)

(2)

dt

С4 7~ = 102 - 72 + 4(и5) , dt

С ^тт = -10 - 71 + Яив,)-dt

Подробнее о процессах построения функции Р(г,и) и уравнений (2) смотрите в [10].

Рассмотрим задачу выделения положений равновесия (т. е. не зависящих от времени t решений уравнений состояния цепи). Для нахождения указанных решений приравняем к нулю правые части уравнений (2): и4 + ив - и3 = 0 , и5 - и4=0 ,

для ограничений I01 = I0 + 4( g1s13 + g3 ^ + s2 )3),

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

I02 = 4(g2s3 - g^), s3 = s1 + s2, s4= — имеется

3

решение:

{ i1 = -I0 - 4g1S1, i2 = 4(g2S2 - g1S1) ,

u3 S1 + S^ , u4 S^ .

u5 S2 , u6 S1 } ;

(6)

(7)

201 + 71 + /з(из) = 0, 72 - 71 + /2 (л 4) - 10=0, 102 - 72 + /4(и5) = 0, 10 + 71 - /^в) = 0. Здесь /1 (ив) = §1ив (ив^1)2, /з (из) = gзUз (из-3яз)2,

/2(и4)= g 2и4(и4 - , /4(и5) = g4U5(U5 - 34)2 -

кубические интерполяционные многочлены Ла-гранжа, аппроксимирующие функции /1(ив ) ,

/2(и4), /3(и3), /4(и5). Яу - первая экстремальная

точка функции / (и,), gJ = )/^).

для ограничений 10 = 101 , 102 = 27g4я2(я2 -я4)2, s3 = s1 + s2 имеется решение:

{ 71= -101, 72=0, щ =3^ + ^2), и4 3^ , и5 3^ , ив 3sl } ; для ограничений 10 = 101 , s3 = 2s1,

102 = 272(я1 -я2)2 + g4(s1 -84)2) имеется решение:

{ 71= -10 , 72=21 g 2 ЯМ - 52^2 , и3=вя1, и4=3я1, и5=3я1, ив=3я1}.

3. Символьно-численный анализ положений равновесий на устойчивость. Возмущенное движение в окрестности положения равновесия запишем в общем виде:

_ о* -О • __•О _ >' О

^ 1 + , 7~2 ^ 2 + 7~2 , Лз ^ 3 + ЛЗ , и4= С4 + и°о, и5= С5 + и°, ив= Св + ив0

(9)

где i1, iO, u30, u°4, u50, uO - постоянные значения

Машиностроение и машиноведение

оо ео I

Modern technologies. System analysis. Modeling, 2018, vol. 60, no. 4

переменных в невозмущенном движении (причем и° = и° ), Zj (j = 1,6) - отклонения от невозмущенного движения.

Исходными данными для пакета Analysis из [7, 10] при исследовании устойчивости служат уравнения (2) и движение (9). По входным данным пакет формирует уравнения возмущенного движения в первом приближении:

L ^ + Z з-Z 4-Z 6=0, dt

l2 % + Z 4 -Z 5=0,

at

3 dt dZ± dt

Сз %-Z + Пз Z 3=0.

C2 %+Zi-Z 2 + П4 Z 4=0.

(10)

С4 ^ + с 2 + щ с 5=0,

т

С, ^ + С} + п6 С6=0, т

где

2 2 n3 = 3gз(35з2-4^3«3 + и3° ); «4 = 3^-+ ы°4 );

2 2 n5 = 3g4(354 - 454и° + и° );

«6 = 3g1 (3512-+ ы°° ) и соответствующие уравнениям (10) коэффициенты характеристического

полинома Л(Я) = ^ к Я6:

ко = С1С2С3С4 ¿1 ^2 ;

к1 = ЦЬ2 (С1 (С3С4 п4 + С2 (С4 п3 + С3п5)) + С2С3С4 п6); к2 = С (С4 (С3 (+ ¿2 ) + Ь1Ь2п3п4 ) +

+ С3 п4 п5 + С2 ((С3 + ¿2 п3п5) + С4 ¿2 )) + + Ь2 (Цп6 (С3С4 п4 + С2 (С4 п3 + С3п5)) + С2С3С4); к3 = п6 (С4 Ь1Ь2 п3п4 + С3 (С4 (Ь1 + ¿2) + Ь1Ь2 п4 п5) + + С2( ¿2(С4 + Ь1п3п5) + С3 ¿1)) + + С1(С3Ь1п4 + С4((11 +Z2)n3 + ¿2 п4) +

+ п5 (С3 (¿1 + ¿2 ) + ^¿2п3п4) + С2 (+ ¿¿2п5 )) +

+¿¿2 (С3С4п4 + С2 (С4п3 + С3п5)); к4 = (С2 + С4) С3 + (С3 п4 п5 + С4 п3 п4 + С2 п3п5) + + С1 (С2 + С3 + С4 + Ь2 (п3 + п4) п5 + Ь1п3 (п4 + п5)) + + п6(С 3Цп4 + С4((11 + Ь2)п3 + ¿¿2 п4) +

+ п5 (С3 (^1 + ¿2) + ¿^п3п4) + С2 (Дп + ¿¿2п5 )) ;

к5 = п6 (С2 + С3+С4 + Ь2 (п3 + п4 )п5 + Ь1п3 (п4 + п5)) + + п3 (С2 + С4 + Ь2 п4 п5) + С3(п4 + п5) + С1 (п3 + п4 + п5);

к6 = (п4 + п5 )п6 + п3 (п4 + п5 + п6). С помощью СКА МмквтаНеа в явном виде через параметры Д , С1-2, п (/ = 3,6) выписаны условия (см. критерий Льенара - Шипара [17]):

Г к

kl >0, (i = 0,6); det

k5 ^

^ k*2 0 k к

>0;

з j

det

Г к: к3 к5 0 01

к0 к2 к4 кб 0

0 к, к3 к5 0

0 к0 к2 к4 кб

v 0 0 к, к3 к5 j

(11)

>0

обеспечивающие существование корней с отрицательными вещественными частями полинома Л(Я) . Следовательно, на основании теоремы Ляпунова об устойчивости по первому приближению [18] исследуемое решение (9) будет асимптотически устойчиво.

Для демонстрации дальнейшего исследования предположим, что п4 = -п5, п3 = п6 . Численное моделирование проводилось с помощью функции Regi°nPl°t СКА МаМвтаНеа при следующих значениях емкости и индуктивности:

С1=1.1*10-3, С2 = 1.2*10-3 *

С3 =1.3* 10-3,

С4 = 1.4 *10-

^ 3 L = 2, L2 = 3 . Условия (11) выполняются в области, представленной на рис. 4.

£ оо!

Рис. 4. Область асимптотической устойчивости

Обратим внимание на решения (5)-(8). Заметим, что если постоянные значения и° (I = 3,6)

ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ПУТЕЙ СООБЩЕНИЯ

Современные технологии. Системный анализ. Моделирование № 4 (60) 2018

в невозмущенном движении принимают значения из решений (5) или (6), тогда все п1 будут равны нулю. Если из решения (7), тогда только п5 не равно нулю. Для значений из решения (8) будут

п3 = 0, п4 ф 0, п5 ф 0, п6 = 0 .

Анализ устойчивости этих частных решений приводит к критическим [19] по Ляпунову случаям. Так, для решений (7) и (8) у характеристического полинома Л(Я) = я(^ к1Я5~') имеется

один нулевой корень, а для решений (5) и (6) характеристический полином факторизуется и имеет

вид Л(Я) = Я2(к0Я4 + ~Я2 + к2) . Исследование устойчивости при п = 0 (I = 3,6) изложено в [16] с применением критерия [20], а при п1 = С1 -2п (I = 3,6) представлено в [15].

Подробнее рассмотрим анализ устойчивости решения (8), для которого, как замечено выше:

п3 = п6= 0, п4 ф 0, п5 ф 0 и Л(Я) = я(£5_0~Я5-) .

Если вещественные части корней характеристического полинома пятой степени отрицательны, то решение будет устойчиво, но не асимптотически. Моделирование с помощью функции Regi°nPl°t (при тех же вышеприведенных численных значениях емкости и индуктивности) соответствующих условий Льенара - Шипара определяет область значений параметров (см. рис. 5), обеспечивающую требуемые свойства корней. Отметим, что параметр п5 (или п4 ) может принимать близкое к нулю отрицательное значение, а для любых п4 > 0 и п5 > 0 будет устойчивость решения (8).

0.10 -

0.08

0.06

0.04

0.02

0.00

-0.02

0.10

Пл

Рис. 5. Область устойчивости решения (8)

Заключение

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, грант № 16-07-00201а.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Computer Algebra: Symbolic and Algebraic Computation / B. Buchberger et al. New York : Springer-Verlag, 1983. 284 p.

2. Davenport J.H., Siret Y., Tournier E. Computer Algebra: Systems and Algorithms for Algebraic Computations. New York : Academic Press, 1988. 313 p.

3. Cohen J.S. Computer Algebra and Symbolic Computation: mathematical methods. A K Peters, Ltd. 2003. 449 p.

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

4. Wolfram S. The Mathematica Book. Fourth Edition. New York : Cambridge University Press., 1999.

5. Графический редактор для формирования символьного описания электрических цепей : а. с. 2018614687 ; опубл. 13.04.2018.

6. Программный комплекс для моделирования в символьном виде механических систем и электрических цепей : а. с. 2016618253 ; опубл. 25.07.2016.

7. Программный комплекс для выделения и исследования устойчивости стационарных множеств механических систем : а. с. 2011612429 ; опубл. 05.07.2011.

8. Синицкий Л.А. Методы аналитической механики в теории электрических цепей. Львов : Вища школа, 1978. 139 с.

Машиностроение

Modern technologies. System analysis. Modeling, 2018, vol. 60, no. 4

9. Иртегов В.Д., Титоренко Т.Н. Об использовании электромеханических аналогий // Системы поддержки принятия решений для исследования и управления энергетикой. Новосибирск : Наука, 1997. С. 136-145.

10. Символьные вычисления в моделировании и качественном анализе динамических систем / А.В. Банщиков и др. // Вычислительные технологии. 2014. T. 19. № 6. С. 3-18.

11. Brayton R.K., Moser J.K. A theory of nonlinear networks-1 // Quarterly ofAppl. Math. 1964. Vol. 22. No. 1. P. 1-33.

12. Краткий обзор графических редакторов структурных моделей сложных систем / И.М. Якимов и др. // Вестник Казан. технолог. ун-та. 2014. Т. 17. № 17. С. 213-221.

13. Окунцев П.В., Соколов Е.А. Графический редактор для составления электрических схем // Материалы восемнадцатой Всерос. конф. Нижневарт. гос. ун-та. Нижневартовск, 2016. С. 1480-1484.

14. Галецкий Ф.П. Импульсный инвертор на туннельных диодах. М. : ИТМ и ВТ АН СССР. 1966.

15. Irtegov V., Titorenko T. On modelling and qualitative investigation of nonlinear systems with the aid of computer algebra // Proceedings of the Sixth International Workshop on Computer Algebra in Scientific Computing. Technische Universität München. 2003. P. 199-211.

16. Титоренко Т.Н. Исследование динамических систем с помощью компьютера // Аналитическая механика, устойчивость и управление движением : тр. IX Междунар. Четаевской конф. Иркутск, 2007. Т. 2. С. 271-283.

17. Гантмахер Ф.Р. Теория матриц. М. : Наука, 1967. 576 с.

18. Ляпунов А.М. Общая задача об устойчивости движения. М.-Л. : ГИТЛ, 1950. 472 с.

19. Четаев Н.Г. Устойчивость движения. Работы по аналитической механике. М. : Изд-во АН СССР, 1962. 535 с.

20. Кац А.М. К вопросу о критерии апериодической устойчивости // Прикладная математика и механика. 1951. Т. 15. Вып. 1. С. 120.

REFERENCES

1. Buchberger B., Collins G.E., Loos R. and Albrecht R. Computer Algebra: Symbolic and Algebraic Computation. Second Edition. Springer-Verlag, New York, 1983, 284 p.

2. Davenport J.H., Siret Y., Tournier E. Computer Algebra: Systems and Algorithms for Algebraic Computations. Academic Press, New York, 1988, 313 p.

3. Cohen J.S. Computer Algebra and Symbolic Computation: mathematical methods. AK Peters Ltd, 2003, 449 p.

4. Wolfram S. The Mathematica Book. Fourth Edition. Cambridge University Press, New York, 1999.

5. Banshchikov A.V., Vetrov A.A. Graficheskii redaktor dlya formirovaniya simvol'nogo opisaniya elektricheskikh tsepei [Graphic editor for the formation of a symbolic description of electrical circuits]. Certificate of state registration of a computer pro-gram. No. 2018614687. FGU FIPS, April 13, 2018.

6. Banshchikov A.V., Irtegov V.D., Titorenko T.N. Programmnyi kompleks dlya modelirovaniya v simvol'nom vide mekhanich-eskikh sistem i elektricheskikh tsepei [Software system for modeling in symbolic form of mechanical systems and electrical circuits]. Certificate of state registration of computer program No. 2016618253 FGU FIPS. July 25, 2016

7. Banshchikov A.V., Burlakova L.A., Irtegov V.D., Titorenko T.N. Programmnyi kompleks dlya vydeleniya i issledovaniya ustoichivosti statsionarnykh mnozhestv mekhanicheskikh sistem [The software package for the isolation and study of the stability of stationary sets of mechanical systems]. Certificate of state registration of computer program No.2011612429. FGU FIPS, July 5, 2011.

8. Sinitskii L.A. Metody analiticheskoi mekhaniki v teorii elektricheskikh tsepei [Methods of analytical mechanics in the theo-ry of electrical circuits]. L'vov: Vishcha shkola Publ., 1978, 139 p.

9. Irtegov V.D., Titorenko T.N. Ob ispol'zovanii elektromekhanicheskikh analogii [On the use of electromechanical analogies]. Sis-temy podderzhki prinyatiya reshenii dlya issledovaniya i upravleniya energetikoi [Decision Support Systems for Energy Research and Management]. Novosibirsk: Nauka Publ., 1997, pp. 136-145.

10. Banshchikov A.V., Burlakova L.A., Irtegov V.D., Titorenko T.N. Simvol'nye vychisleniya v modelirovanii i kachestvennom an-alize dinamicheskikh sistem [Symbolic calculations in modeling and qualitative analysis of dynamic systems]. Vychislitel'nye tekhnologii [Computational technologies], 2014, Vol. 19, No. 6, pp. 3-18.

11. Brayton R.K., Moser J.K. A theory of nonlinear networks-1. Quarterly of Appl. Math., 1964, Vol. 22, No.1, pp. 1-33.

12. Yakimov I.M., Abzalova L.R., Kirpichnikov A.P., Mokshin V.V. Kratkii obzor graficheskikh redaktorov strukturnykh mod-elei slozhnykh system [A brief review of graphical editors of structural models of complex systems]. Vestnik Kazanskogo tekhno-logicheskogo universiteta [The Herald of Kazan Technological University], 2014 Vol. 17, No. 17, pp. 213-221.

13. Okuntsev P.V., Sokolov E.A. Graficheskii redaktor dlya sostavleniya elektricheskikh skhem [Graphic editor for drawing up electrical circuits]. Materialy vosemnadtsatoi vserossiiskoi konferentsii Nizhnevartovskogo gosudarstvennogo universiteta [Proceedings of the Eighteenth All-Russian Conference of Nizhnevartovsk State University]. Nizhnevartovsk, 2016, pp. 1480-1484.

14. Galetskii F.P. Impul'snyi invertor na tunnel'nykh diodakh [Pulse inverter on tunnel diodes]. Moscow: IPM CE of AS USSR Publ., 1966.

15. Irtegov V., Titorenko T. On modelling and qualitative investigation of nonlinear systems with the aid of computer algebra. Proceedings of the Sixth International Workshop on Computer Algebra in Scientific Computing. Technische Universität München, 2003, pp. 199-211.

16. Titorenko T.N. Issledovanie dinamicheskikh sistem s pomoshch'yu komp'yutera [Study of dynamic systems using a comput-er]. Trudy IX mezhdunarodnoi Chetaevskoi konferentsii «Analiticheskaya mekhanika, ustoichivost' i upravlenie dvizheniem» [Proceed-ings of the IX International Chetaev Conference "Analytical Mechanics, Stability and Motion Control"]. Irkutsk, 2007, Vol. 2, pp. 271-283.

17. Gantmakher F.R. Teoriya matrits [Theory of matrices]. Moscow: Nauka Publ., 1967, 576 p.

18. Lyapunov A.M. Obshchaya zadacha ob ustoichivosti dvizheniya [The general problem of motion stability]. Moscow-Leningrad: GITL Publ., 1950, 472 p.

19. Chetaev N.G. Ustoichivost' dvizheniya. Raboty po analiticheskoi mekhanike [Stability of motion. Works on analytical mechanics]. Moscow: AS USSR Publ., 1962, 535 p.

ИРКУТСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ПУТЕЙ СООБЩЕНИЯ

Современные технологии. Системный анализ. Моделирование № 4 (60) 2018

20. Kats A.M. K voprosu o kriterii aperiodicheskoi ustoichivosti [On the criterion of aperiodic stability]. Prikladnaya matematika i mekhanika [Applied Mathematics and Mechanics], 1951, Vol. 15, Iss. 1, p. 120.

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

Банщиков Андрей Валентинович - к. ф.-м. н., доцент, старший научный сотрудник, Институт динамики систем и теории управления имени В. М. Матросова Сибирского отделения Российской академии наук, г. Иркутск, e-mail: bav@icc.ru

Ветров Александр Анатольевич - программист, Институт динамики систем и теории управления имени В. М. Матросова Сибирского отделения Российской академии наук, г. Иркутск, e-mail: vetrov@icc.ru

Иртегов Валентин Дмитриевич - д. ф.-м. н., доцент, старший научный сотрудник, Институт динамики систем и теории управления имени В. М. Матросова Сибирского отделения Российской академии наук, г. Иркутск, e-mail: irteg@icc.ru

Титоренко Татьяна Николаевна - к. т. н., старший научный сотрудник, Институт динамики систем и теории управления имени В. М. Матросова Сибирского отделения Российской академии наук, г. Иркутск, e-mail: titor@icc.ru

Authors

Banshchikov Andrei Valentinovich - Ph.D. in Physics and Mathematics, Assoc. Prof., senior research officer, Matrosov Institute for System Dynamics and Control Theory of Siberian Branch of the Russian Academy of Sciences, Irkutsk, the Russian Federation, Irkutsk, e-mail: bav@icc.ru

Vetrov Aleksandr Anatol'evich - programming specialist, Matrosov Institute for System Dynamics and Control Theory of Siberian Branch of the Russian Academy of Sciences, Irkutsk, the Russian Federation, e-mail: vetrov@icc.ru

Irtegov Valentin Dmitrievich - D. Sc. in Physics and Mathematics, Assoc. Prof., senior research officer, Matrosov Institute for System Dynamics and Control Theory of Siberian Branch of the Russian Academy of Sciences, Irkutsk, e-mail: irteg@icc.ru

Titorenko Tat'yana Nikolaevna - Ph.D. in Engineering Science, senior research officer, Matrosov Institute for System Dynamics and Control Theory of Siberian Branch of the Russian Academy of Sciences, Irkutsk, e-mail: titor@icc.ru

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

Банщиков А. В. Анализ электрических цепей средствами компьютерной алгебры в машиноведении / А. В. Банщиков, А. А. Ветров, В. Д. Иртегов, Т. Н. Титоренко // Современные технологии. Системный анализ. ние. - 2018. - Т. 60, № 4. - С. 46-54. - Б01: 10.26731/1813-9108.2018.4(60).46-54

For citation

Banshchikov A. V., Vetrov A. A., Irtegov V. D., Titorenko T. N. Analiz elektricheskikh tsepei sredstvami komp'yuternoi algebry v mashinovedenii [Analysis of electrical circuits by means of computer algebra in machine science]. Sovremennye tekhnologii. Sis-temnyi analiz. Modelirovanie [Modern Technologies. System Analysis. Modeling], 2018, Vol. 60, No. 4, pp. 46-54. DOI: 10.26731/1813-9108.2018.4(60).46-54

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