Расчет переходного процесса численными методами.

2.3.1. Явный метод Рунге-Кутта 4-ого порядка.

Численные методы расчета переходных процессов предполагают представление системы дифференциальных уравнений Кирхгофа, записанных для схемы после коммутации в конечно-разностном виде. Представление временных производных в конечно-разностном виде в простейшем варианте соответствует так называемому методу ломаной Эйлера, но даже в этом случае возможны разные способы его реализации: явный и неявный метод Эйлера. При явном методе Эйлера приращение функции Δy на интервале времени Расчет переходного процесса численными методами. - student2.ru ,

Расчет переходного процесса численными методами. - student2.ru .

При этом значение Расчет переходного процесса численными методами. - student2.ru в начале интервала Расчет переходного процесса численными методами. - student2.ru , то есть в момент tk-1 считается известным. Например, для простейшего дифференциального уравнения

Расчет переходного процесса численными методами. - student2.ru .

Такое представление соответствует так называемой правой производной (рис. 10).

Расчет переходного процесса численными методами. - student2.ru
рис. 10

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

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

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

При неявном методе Эйлера правая часть дифференциального уравнения определяется по заранее неизвестному значению искомой функции y(tk) в конце интервала времени Δtk. В этом варианте используется так называемая левая производная при конечно-разностном ее представлении.

Расчет переходного процесса численными методами. - student2.ru

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

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

Наиболее распространенный из них является метод Рунге-Кутта. Современные программные системы имеют стандартные процедуры решения систем дифференциальных уравнений. В системе MathCAD процедура решения ДУ по методу Рунге-Кутта 4-ого порядка имеет вид:

S:=rkfixed(X,TO,Tmax,N,d),

где X – матрица-столбец неизвестных;

TO – начальное значение аргумента (начало переходного процесса, соответствующее моменту коммутации, );

Tmax – длительность переходного процесса;

N – число дискретных интервалов времени Δt на отрезке Tmax;

d – матрица-столбец правой части решаемой системы дифференциальных уравнений, записанной в нормальной форме Коши.

Таким образом, для численного расчета переходного процесса в рассмотренной ранее схеме необходимо представить систему уравнений Кирхгофа (1) в нормальной форме.

Для этого необходимо свести систему (1) к двум уравнениям относительно тока iL(t) и напряжения uC(t) с учетом уравнений связи Расчет переходного процесса численными методами. - student2.ru , Расчет переходного процесса численными методами. - student2.ru . Процедура сведения системы (1) к двум уравнениям, соответствует теоретическим положениям метода переменных состояния. Не вдаваясь в подробности этого метода проведем необходимые преобразования стандартными алгебраическими операциями.

Из первых трех уравнений системы (1) выразим i1(t), i5(t) и i6(t)

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru .

После подстановки полученных выражений в (4) и (5) уравнения системы (1) получим:

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru

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

Расчет переходного процесса численными методами. - student2.ru

или

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru , Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru , Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru , Расчет переходного процесса численными методами. - student2.ru

Результаты расчета переходного процесса с помощью процедуры rkfixedприведены ниже.

В качестве начальных условий взяты значения Расчет переходного процесса численными методами. - student2.ru и Расчет переходного процесса численными методами. - student2.ru , рассчитанные ранее (п.2.1).

Расчет переходного процесса численными методами. - student2.ru , , , , Расчет переходного процесса численными методами. - student2.ru

Графики iL(t) и uC(t) приведены на рис. 11

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Напряжение на емкости Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Ток в индуктивности Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru
рис. 11

2.3.2. Метод дискретных резистивных схем.

Численный метод расчета переходного процесса в ЭЦ на основе неявного метода Эйлера сводится к расчету дискретной резистивной схемы. Такая схема получается для дискретного момента времени tk с учетом конечно-разностного представления производных в уравнениях Кирхгофа в послекоммутационном режиме.

Действительно, напряжение на индуктивности в момент времени tk

Расчет переходного процесса численными методами. - student2.ru

Полученное выражение соответствует замене в ЭЦ элемента индуктивности L источником напряжения Расчет переходного процесса численными методами. - student2.ru и активным сопротивление RL (рис. 12)

Расчет переходного процесса численными методами. - student2.ru
рис. 12

Дискретный резистивный элемент сопротивление и источник ЭДС Расчет переходного процесса численными методами. - student2.ru заменяют в момент времени tk переходного процесса элемент индуктивность L.

Аналогичные выражения и схему замещения (рис. 13) можно получить для элемента емкость C.

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru
рис. 13

Дискретный резистивный элемент проводимость и источник тока Расчет переходного процесса численными методами. - student2.ru заменяют в момент времени tk переходного процесса элемент емкость C.

В результате для дискретного момента времени tk будем иметь дискретную схему ЭЦ (рис. 14)

Расчет переходного процесса численными методами. - student2.ru
рис. 14

Схема рис. 14 называется дискретной из-за наличия в ней дискретных источников Расчет переходного процесса численными методами. - student2.ru и Расчет переходного процесса численными методами. - student2.ru , значения которых изменяются на каждом последующем интервале времени Δt переходного процесса.

Расчет переходных токов в этом случае сводится к решению системы линейных алгебраических уравнений Кирхгофа для дискретной схемы (рис. 14) на каждом интервале Δt переходного процесса. Следовательно придется решать алгебраическую систему N раз, где N – количество дискретных интервалов времени Δt, на которые разбит весь отрезок длительности процесса Tmax. При этом на первом интервале Δt значения источников Расчет переходного процесса численными методами. - student2.ru , Расчет переходного процесса численными методами. - student2.ru определяются ННУ схемы. На каждом последующем интервале эти источники меняют свое значение, но матрица коэффициентов решаемой системы алгебраических уравнений не изменяется, меняется только матрица – столбец правой части системы.

Алгебраическая система уравнений Кирхгофа для дискретной схемы рис. 14 имеет вид:

1. Расчет переходного процесса численными методами. - student2.ru

2. Расчет переходного процесса численными методами. - student2.ru

3. Расчет переходного процесса численными методами. - student2.ru

4. Расчет переходного процесса численными методами. - student2.ru

5. Расчет переходного процесса численными методами. - student2.ru

6. Расчет переходного процесса численными методами. - student2.ru

В матричной форме система уравнений Кирхгофа примет вид

i1 iL iC iGC i5 i6

Расчет переходного процесса численными методами. - student2.ru ,

В матричной форме

Расчет переходного процесса численными методами. - student2.ru ,

Решение матричного уравнения

Расчет переходного процесса численными методами. - student2.ru

Численная реализация описанного алгоритма приведена в системе MathCAD ниже

Начальные значения токов схемы рис. 14 взяты из расчета ЗНУ классическим методом. Значение начального тока iGc в дискретной схеме определяется по 1 закону Кирхгофа для узла †.

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Расчет переходного процесса численными методами. - student2.ru Расчет переходного процесса численными методами. - student2.ru

Результаты расчета аналитическими методами (классическим и операторным) и численными методами (Рунге-Кутта и дискретной схемы) абсолютно совпадают.

Наши рекомендации