Метод узловых потенциалов

Материал из Википедии — свободной энциклопедии
Перейти к: навигация, поиск

Метод узловы́х потенциалов — метод расчета электрических цепей путём записи системы линейных алгебраических уравнений, в которой неизвестными являются потенциалы в узлах цепи. В результате применения метода определяются потенциалы во всех узлах цепи, а также, при необходимости, токи во всех ветвях.

Введение[править | править исходный текст]

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

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

Метод узловых потенциалов не привносит ничего нового к правилам Кирхгофа и закону Ома. Данный метод лишь формализует их использование настолько, чтобы их можно было применить к любой, сколь угодно сложной цепи и пригоден для расчёта посредством компьютеров. Иными словами, метод даёт ответ на вопрос «как использовать законы для расчета данной цепи?».

Теоретические основы[править | править исходный текст]

Если в цепи, состоящей из У узлов и Р рёбер, известны все характеристики звеньев (полные сопротивления R, величины источников ЭДС E и тока J), то возможно вычислить токи Ii во всех рёбрах и потенциалы φi во всех узлах. Поскольку электрический потенциал определён с точностью до произвольного постоянного слагаемого, то потенциал в одном из узлов (назовём его базовым узлом) можно принять равным нулю, а потенциалы в остальных узлах определять относительно базового узла. Таким образом, при расчёте цепи имеем У+Р–1 неизвестных переменных: У–1 узловых потенциалов и Р токов в рёбрах.

Не все из указанных переменных независимы. Например, исходя из закона Ома для участка цепи, токи в звеньях полностью определяются потенциалами в узлах:

\ I_i = \frac{\phi_A-\phi_B+E_i}{R_i} + J_i.

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

\ \phi_B = \phi_A + E_i + (J_i-I_i)R_i.

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

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

\ Y - P + K = 1

или

\ P = Y + K - 1,

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

Уравнение для потенциала в узлах[править | править исходный текст]

Рис. 1. Фрагмент цепи: узел с примыкающими звеньями

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

\sum_{i=1}^n I_i = 0.

Ток в звене определим, исходя из закона Ома для участка цепи:

 I_i = \frac{\phi_i-\phi+E_i}{R_i} + J_i

откуда

\sum_{i=1}^n \left( \frac{\phi_i-\phi+E_i}{R_i} + J_i \right) = 0;
\phi \sum_{i=1}^n \frac{1}{R_i} - \sum_{i=1}^n \frac{\phi_i}{R_i} = 
\sum_{i=1}^n \left( \frac{E_i}{R_i} + J_i \right).

Обозначив проводимости рёбер через

 Y_i = \frac{1}{R_i},

получим окончательное уравнение для узла

\phi \sum_{i=1}^n Y_i - \sum_{i=1}^n \phi_i Y_i = 
\sum_{i=1}^n (E_i Y_i + J_i).

Последнее уравнение получено, исходя из предположения, что все источники тока и ЭДС направлены в сторону рассматриваемого узла. Если какой-либо источник направлен в противоположную сторону, его ЭДС или ток необходимо взять с обратным знаком.

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

Практическое применение[править | править исходный текст]

Составление системы уравнений[править | править исходный текст]

Перед началом расчёта выбирается один из узлов (базовый узел), потенциал которого считается равным 0. Затем узлы нумеруются, после чего составляется система уравнений.

Уравнения составляются для каждого узла, кроме базового. Слева от знака равенства записывается:

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

Справа от знака равенства записывается:

  • сумма всех источников токов, примыкающих к данному узлу;
  • сумма произведений всех ЭДС, примыкающих к данному узлу, на проводимость соответствующего звена.

Если источник направлен в сторону рассматриваемого узла, то он записывается со знаком «+», в противном случае — со знаком «−».

Рис. 2. Пример электрической схемы

Пример системы уравнений[править | править исходный текст]

На схеме (рис. 2) четыре узла. Потенциал в узле 0 принят равным нулю (φ0 = 0). Записываем уравнения для узлов 1, 2 и 3:


\begin{cases}
  \phi_1 (Y_1+Y_4+Y_6) + \phi_2 (-Y_1) + \phi_3 (-Y_6) = E_6 Y_6 - E_4 Y_4 \\
  \phi_1 (-Y_1) + \phi_2 (Y_1+Y_2+Y_3) +\phi_3 (-Y_3) = 0 \\
  \phi_1 (-Y_6) + \phi_2 (-Y_3) + \phi_3 (Y_3+Y_5+Y_6) = J_5 - E_6 Y_6
\end{cases},

где проводимости рёбер равны

\ Y_1 = \frac{1}{R_1}; \quad Y_2 = \frac{1}{R_2}; \quad Y_3 = \frac{1}{R_3};
\ Y_4 = \frac{1}{R_4}; \quad Y_5 = \frac{1}{R_5}; \quad Y_6 = \frac{1}{R_6}.

Формальный подход[править | править исходный текст]

В матричном виде система уравнений для метода узловых потенциалов выглядит следующим образом[1]:

 \mathbf{A Y A^t U_0 = - A (J+YE)},

где

 \mathbf A — матрица соединений размера (q – 1) × p (q — количество узлов, р — количество рёбер) , в которой i–я строка соответствует узлу i, а j–й столбец соответствует ребру j, причём элемент Aij равен

  • 0, если ребро j не присоединено к узлу i;
  • 1, если ребро выходит из узла;
  • –1, если ребро входит в узел.

Понятие «входит» и «выходит» означает, что для каждого ребра задаётся направление, которое обычно ассоциируется с направлением тока в этом ребре;

 \mathbf Y — диагональная матрица проводимостей размера p × p, в которой диагональный элемент Yii равен проводимости i–го ребра, а недиагональные элементы равны нулю;

 \mathbf A^t — транспонированная матрица соединений;

 \mathbf U_0 — матрица-столбец узловых потенциалов размером (q – 1) × 1. Потенциалы измеряется относительно предварительно выбранного узла, потенциал которого считается равным нулю. Нулевой узел не входит ни в одну из перечисленных в данном разделе матриц;

 \mathbf J — матрица-столбец источников тока размером p × 1, где каждый элемент равен току соответствующего источника, причём эта величина нулевая, если в данном ребре источник тока отсутствует; положительная, если направление тока источника совпадает с направлением тока в ребре; и отрицательная в противном случае;

 \mathbf E — матрица-столбец источников ЭДС размером p × 1, где каждый элемент равен ЭДС соответствующего источника, причём эта величина нулевая, если в данном ребре источник ЭДС отсутствует; положительная, если направление ЭДС источника совпадает с направлением тока в ребре; и отрицательная в противном случае.

Пример системы уравнений[править | править исходный текст]

Для схемы рис. 2 матрицы имеют вид:


\mathbf A = \begin{pmatrix} 
  1 &  0 &  0 &  1 &  0 & -1 \\ 
 -1 &  1 &  1 &  0 &  0 &  0 \\
  0 &  0 & -1 &  0 & -1 &  1 
\end{pmatrix}; \quad
\mathbf U_0 = \begin{pmatrix} 
 \phi_1 \\ \phi_2 \\  \phi_3 
\end{pmatrix}


\mathbf A^t = \begin{pmatrix} 
 1&-1&0 \\ 
 0&1&0 \\
 0&1&-1 \\
 1&0&0 \\
 0&0&-1 \\
 -1&0&1 \\
\end{pmatrix}; \quad
\mathbf Y = \begin{pmatrix} 
 Y_1&0&0&0&0&0 \\ 
 0&Y_2&0&0&0&0 \\
 0&0&Y_3&0&0&0 \\
 0&0&0&Y_4&0&0 \\
 0&0&0&0&Y_5&0 \\
 0&0&0&0&0&Y_6 \\
\end{pmatrix}; \quad 
\mathbf J = \begin{pmatrix} 
 0 \\ 0 \\ 0 \\ 0 \\ J_5 \\ 0 
\end{pmatrix}; \quad
\mathbf E = \begin{pmatrix} 
 0 \\ 0 \\ 0 \\ E_4 \\ 0 \\ E_6 
\end{pmatrix}

Перемножаем матрицы в соответствии с матричным уравнением:


\mathbf {AY} = \begin{pmatrix} 
  Y_1 &  0 &  0 & Y_4 &  0 & -Y_6 \\ 
 -Y_1 &  Y_2 & Y_3 &  0 &  0 &  0 \\
  0 &  0 & -Y_3 &  0 & -Y_5 & Y_6 
\end{pmatrix};


\mathbf {AYA^t} = \begin{pmatrix} 
 Y_1+Y_4+Y_6 &  -Y_1 & -Y_6 \\ 
 -Y_1 & Y_1+Y_2+Y_3 &  -Y_3 \\
  -Y_6 &  -Y_3 & Y_3+Y_5+Y_6 
\end{pmatrix};


\mathbf {AYA^t U_0} = \begin{pmatrix} 
 (Y_1+Y_4+Y_6) \cdot \phi_1 -Y_1 \cdot \phi_2 -Y_6 \cdot \phi_3 \\ 
 -Y_1\cdot \phi_1 + (Y_1+Y_2+Y_3) \cdot \phi_2  -Y_3\cdot \phi_3 \\
  -Y_6\cdot \phi_1   -Y_3 \cdot \phi_2 + (Y_3+Y_5+Y_6) \cdot \phi_3 
\end{pmatrix};


\mathbf {J+YE} = \begin{pmatrix} 
 0 \\ 0 \\ 0 \\ Y_4E_4 \\ J_5 \\ Y_6E_6 
\end{pmatrix}; \quad 
\mathbf {-A(J+YE)} = \begin{pmatrix} 
 -Y_4E_4+Y_6E_6 \\ 0 \\ J_5-Y_6E_6 
\end{pmatrix}

Раскрывая матричную запись, получаем следующую систему уравнений:


\begin{cases}
  (Y_1+Y_4+Y_6) \cdot \phi_1 -Y_1 \cdot \phi_2 -Y_6 \cdot \phi_3 = - E_4 Y_4 + E_6 Y_6 \\
  -Y_1\cdot \phi_1 + (Y_1+Y_2+Y_3) \cdot \phi_2  -Y_3\cdot \phi_3 = 0 \\
  -Y_6\cdot \phi_1   -Y_3 \cdot \phi_2 + (Y_3+Y_5+Y_6) \cdot \phi_3 = J_5 - E_6 Y_6
\end{cases}

Ограничения[править | править исходный текст]

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

Примечания[править | править исходный текст]

  1. Нейман Л.Р., Демирчян К.С. Теоретические основы электротехники: в 2-х т. Учебник для вузов. Том I. — 3-е изд., перераб. и доп. — Л.: Энергоиздат. Ленингр. отд-ние, 1981. — 536 с., ил.

См. также[править | править исходный текст]