Систем дифференциальных уравнений в электрических цепях (2 видео)

№74 Расчет переходных процессов методом переменных состояния.

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

В более узком смысле в математике уравнениями состояния называют систему дифференциальных уравнений 1-го порядка, разрешенных относительно производных (форма Коши). Система уравнений состояния в обобщенной форме имеет вид:

Та же система уравнений в матричной форме:

или в обобщённой матричной форме:

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

Значения производных на к-ом шаге:

Значения переменных на к-ом шаге:

Для определения значений переменных и их производных на 1-м шаге ин¬тегрирова¬ния используются их значения на момент t=0, т.е. их начальные условия x1(0), x2(0). xn(0).

Уравнения состояния формы Коши для заданной схемы могут быть получены из системы уравнений Кирхгофа путем их преобразования. Для этой цели: а) из системы уравнений Кирхгофа методом подстановки исключаются »лишние» переменные, имеющие зависимые начальные условия, и оставляют переменные iL(t) и uC(t), которые не изменяются скачком и имеют независи-мые начальные условия iL(0) и uC(0); б) оставшиеся уравнения решаются относительно производных и приводятся их к форме Коши.

В случае сложных схем уравнения состояния формы Коши могут быть составлены топологическими методами с использованием матриц соединений [A] и [B].

Последовательность расчета переходного процесса методом переменных состояния выглядит так:

1. Производится расчет схемы в установившемся режиме до коммутации и определяются независимые начальные условия iL(0) и uC(0).

2. Составляется система дифференциальных уравнений по законам Кирхгофа для схемы после коммутации.

3. Методом исключения »лишних» переменных система уравнений Кирхгофа преобразуется в систему уравнений Коши, составляются матрицы коэффициентов.

4. Выбирается расчетное время (продолжительность переходного процесса) и число шагов интегрирования N.

5. Решение задачи выполняется на ЭВМ по стандартной программе. Выходную функцию получают в виде графической диаграммы x=f(t)или в виде таблицы координат функций для заданных моментов времени.

Пример. Для схемы рис. 74.1 с заданными параметрами элементов (e(t)=Emsin(ωt+ψE), R, R1, R2, R3, L1, L2, C) выполнить расчет переходного процесса и определить функцию uab(t).

1. Выполняется расчет схемы в установившемся режиме переменного тока до коммутации и определяются начальные условия i1(0), i2(0), uC(0).

2. Составляется система дифференциальных уравнений по законам Кирхгофа:

3. Система уравнений Кирхгофа преобразуется в систему уравнений Коши.

Для этой цели из (1) выражаем

и делаем подстановку в (1) и (2), а из (4) делаем подстановку в (1). Тогда получим:

Подсчитаем значения отднльных коэфициэнтов:

Составляем матрицы коэффициентов:

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

Число шагов интегрирования принимаем N = 1000,

Вводим исходные данные в ЭВМ и выполняем рассчет.

В качестве выходной функции принимаем:

Для выходной функции Uab(T) строим графическую диаграмму в интервале периода Т.

Видео:Переходные процессы | Классический метод расчета переходных процессов. Теория и задачаСкачать

Систем дифференциальных уравнений в электрических цепях

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

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

Основные методы анализа переходных процессов в линейных цепях:

Классический метод, заключающийся в непосредственном интегрировании дифференциальных уравнений, описывающих электромагнитное состояние цепи.
Операторный метод, заключающийся в решении системы алгебраических уравнений относительно изображений искомых переменных с последующим переходом от найденных изображений к оригиналам.
Частотный метод, основанный на преобразовании Фурье и находящий широкое применение при решении задач синтеза.
Метод расчета с помощью интеграла Дюамеля, используемый при сложной форме кривой возмущающего воздействия.
Метод переменных состояния, представляющий собой упорядоченный способ определения электромагнитного состояния цепи на основе решения системы дифференциальных уравнений первого прядка, записанных в нормальной форме (форме Коши).

Классический метод расчета

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

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

Таблица 1. Связь мгновенных значений напряжений и токов на элементах электрической цепи

;

при наличии магнитной связи с катушкой, обтекаемой током ,

;

Для последовательной цепи, содержащей линейные резистор R, катушку индуктивности L и конденсатор С, при ее подключении к источнику с напряжением u (см. рис. 1) можно записать

Резистор (идеальное активное сопротивление)

Катушка индуктивности (идеальная индуктивность)

Конденсатор (идеальная емкость)

Систем дифференциальных уравнений в электрических цепях

(1)

Подставив в (1) значение тока через конденсатор

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

В общем случае уравнение, описывающее переходный процесс в цепи с n независимыми накопителями энергии, имеет вид:

(2)

где х – искомая функция времени (напряжение, ток, потокосцепление и т.п.); — известное возмущающее воздействие (напряжение и (или) ток источника электрической энергии); — к-й постоянный коэффициент, определяемый параметрами цепи.

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

В общем случае порядок дифференциального уравнения определяется соотношением

(3)

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

Наличие индуктивных связей на порядок дифференциального уравнения не влияет.

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

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

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

В соответствии с вышесказанным, общее решение уравнения (2) имеет вид

(4)

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

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

Начальные условия. Законы коммутации

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

Таблица 2. Законы коммутации

Первый закон коммутации (закон сохранения потокосцепления)

Магнитный поток, сцепленный с катушками индуктивности контура, в момент коммутации сохраняет то значение, которое имел до коммутации, и начинает изменяться именно с этого значения: .

Второй закон коммутации (закон сохранения заряда)

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

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

На практике, за исключением особых случаев (некорректные коммутации), допустимо использование указанных законов в другой формулировке, а именно:

первый закон коммутации – в ветви с катушкой индуктивности ток в момент

коммутации сохраняет свое докоммутационное значение и в дальнейшем начинает изменяться с него: .

второй закон коммутации – напряжение на конденсаторе в момент

коммутации сохраняет свое докоммутационное значение и в дальнейшем начинает изменяться с него: .

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

Действительно, при переводе в схеме на рис. 2,а ключа из положения 1 в положение 2 трактование второго закона коммутации как невозможность скачкообразного изменения напряжения на конденсаторе приводит к невыполнению второго закона Кирхгофа . Аналогично при размыкании ключа в схеме на рис. 2,б трактование первого закона коммутации как невозможность скачкообразного изменения тока через катушку индуктивности приводит к невыполнению первого закона Кирхгофа . Для данных схем, исходя из сохранения заряда и соответственно потокосцепления, можно записать:

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

Пример. Определить токи и производные и в момент коммутации в схеме на рис. 3, если до коммутации конденсатор был не заряжен.

В соответствии с законами коммутации

и .

На основании второго закона Кирхгофа для момента коммутации имеет место

и .

Для известных значений и из уравнения

определяется .

Значение производной от напряжения на конденсаторе в момент коммутации (см. табл. 1)

Корни характеристического уравнения. Постоянная времени

Выражение свободной составляющей общего решения х дифференциального уравнения (2) определяется видом корней характеристического уравнения (см. табл. 3).

Таблица 3. Выражения свободных составляющих общего решения

Вид корней характеристического уравнения

Выражение свободной составляющей

Корни вещественные и различные

Корни вещественные и

Пары комплексно-сопряженных корней

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

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

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

которое называется декрементом колебания, или натуральным логарифмом этого отношения

называемым логарифмическим декрементом колебания, где .

Важной характеристикой при исследовании переходных процессов является постоянная времени t , определяемая для цепей первого порядка, как:

где р – корень характеристического уравнения.

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

Основы теории цепей: Учеб. для вузов /Г.В.Зевеке, П.А.Ионкин, А.В.Нетушил, С.В.Страхов. –5-е изд., перераб. –М.: Энергоатомиздат, 1989. -528с.
Бессонов Л.А. Теоретические основы электротехники: Электрические цепи. Учеб. для студентов электротехнических, энергетических и приборостроительных специальностей вузов. –7-е изд., перераб. и доп. –М.: Высш. шк., 1978. –528с.
Теоретические основы электротехники. Учеб. для вузов. В трех т. Под общ. ред. К.М.Поливанова. Т.1. К.М.Поливанов. Линейные электрические цепи с сосредоточенными постоянными. –М.: Энергия- 1972. –240с.

Видео:2020 г. Дифференциальные уравнения для электрических цепей. Лекция и практикаСкачать

Расчет электрической цепи постоянного и переменного тока

Расчет переходных процессов методом численного интегрирования дифференциальных уравнений на ЭВМ

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

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

В дифференциальных уравнениях дифференциалы функций заменяются их конечными приращениями, а производные функций — отношениями приращений:

На каждом шаге интегрирования решается система дифференциальных уравнений, в результате решения определяются численные значения производных и самих функций. В качестве исходных данных для их определения используются значения этих же функций на предыдущем шаге, а на начальном 1-ом шаге – их значения в момент коммутации при t =0 , т.е. начальные условия. В результате расчета для функций и их производных составляются массивы их значений в исследуемом интервале времени Т, которые после завершения цикла подвергаются соответствующей математической обработке, а именно: строятся графические диаграммы функций, составляются необходимые таблицы, исследуются функции на наличие максимумов и минимумов, устанавливается продолжительность переходного процесса и его характер, и т.д.

Пример. Рассчитать переходный процесс в схеме рис. 153 с заданными параметрами элементов: , R1, R2, R3, L1, L2, С.

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

По законам Кирхгофа для схемы после коммутации составляется система дифференциальных уравнений:

Выбирается шаг интегрирования h (например, из расчета N=1000 шагов на период Т=0,02 с переменного тока, тогда h=Т/ N =2·10-5 с).

Составляется алгоритм вычислений для произвольного к-го шага:

Далее следуют вычисления по тому же алгоритму для (к+1)-го шага и т. д.

В соответствии с составленным алгоритмом на любом языке составляется программа вычислений на ЭВМ, что представляет собой несложную инженерную задачу.

В настоящее время метод численного интегрирования является наиболее универсальным и наиболее простым методом расчета переходных процессов в электрических цепях. Достоинствами метода являются:

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

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

К недостаткам метода следует отнести необходимость составления индивидуальной расчетной программы для каждой конкретной задачи и решение ее на ЭВМ, что сегодня уже посильно каждому инженеру.

Метод эквивалентных преобразований. Сущность метода заключается в том, чтобы сложную разветвленную цепь с помощью эквивалентных преобразований привести к простейшей одноконтурной цепи, включающей ветвь с искомым током, значение которого определяется затем по закону Ома. К эквивалентным преобразованиям относятся: а) преобразование представления источников электрической энергии; б) замена последовательных и параллельных соединений однотипных элементов эквивалентными одиночными элементами; в) преобразование соединений «звезда»–«треугольник» и «треугольник»–«звезда».