Русская Википедия:Метод узловых потенциалов

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

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

Введение

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

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

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

Теоретические основы

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

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

<math>I_i = \frac{\varphi_A-\varphi_B+E_i}{R_i} + J_i.</math>

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

<math>\varphi_B = \varphi_A + E_i + (J_i-I_i)R_i.</math>

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

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

<math>Y - P + K = 1,</math>

или

<math>P = Y + K - 1,</math>

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

Уравнение для потенциала в узлах

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

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

<math>\sum_{i=1}^n I_i = 0.</math>

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

<math>I_i = \frac{\varphi_i-\varphi+E_i}{R_i} + J_i</math>

откуда:

<math>\sum_{i=1}^n \left( \frac{\varphi_i-\varphi+E_i}{R_i} + J_i \right) = 0; </math>
<math>\varphi \sum_{i=1}^n \frac{1}{R_i} - \sum_{i=1}^n \frac{\varphi_i}{R_i} =

\sum_{i=1}^n \left( \frac{E_i}{R_i} + J_i \right).</math>

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

<math>Y_i = \frac{1}{R_i}</math>

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

<math>\varphi \sum_{i=1}^n Y_i - \sum_{i=1}^n \varphi_i Y_i =

\sum_{i=1}^n (E_i Y_i + J_i).</math>

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

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

Практическое применение

Составление системы уравнений

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

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

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

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

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

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

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

Пример системы уравнений

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

<math> \begin{cases}

 \varphi_1 (Y_1+Y_4+Y_6) + \varphi_2 (-Y_1) + \varphi_3 (-Y_6) = E_6 Y_6 - E_4 Y_4 \\
 \varphi_1 (-Y_1) + \varphi_2 (Y_1+Y_2+Y_3) +\varphi_3 (-Y_3) = 0 \\
 \varphi_1 (-Y_6) + \varphi_2 (-Y_3) + \varphi_3 (Y_3+Y_5+Y_6) = J_5 - E_6 Y_6

\end{cases}, </math>

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

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

Формальный подход

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

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

где

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

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

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

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

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

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

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

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

Пример системы уравнений

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

<math> \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}

\varphi_1 \\ \varphi_2 \\  \varphi_3 

\end{pmatrix} </math>

<math> \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} </math>

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

<math> \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}; </math>

<math> \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}; </math>

<math> \mathbf {AYA^t U_0} = \begin{pmatrix}

(Y_1+Y_4+Y_6) \cdot \varphi_1 -Y_1 \cdot \varphi_2 -Y_6 \cdot \varphi_3 \\ 
-Y_1\cdot \varphi_1 + (Y_1+Y_2+Y_3) \cdot \varphi_2  -Y_3\cdot \varphi_3 \\
 -Y_6\cdot \varphi_1   -Y_3 \cdot \varphi_2 + (Y_3+Y_5+Y_6) \cdot \varphi_3 

\end{pmatrix}; </math>

<math> \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} </math>

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

<math> \begin{cases}

 (Y_1+Y_4+Y_6) \cdot \varphi_1 -Y_1 \cdot \varphi_2 -Y_6 \cdot \varphi_3 = - E_4 Y_4 + E_6 Y_6 \\
 -Y_1\cdot \varphi_1 + (Y_1+Y_2+Y_3) \cdot \varphi_2  -Y_3\cdot \varphi_3 = 0 \\
 -Y_6\cdot \varphi_1   -Y_3 \cdot \varphi_2 + (Y_3+Y_5+Y_6) \cdot \varphi_3 = J_5 - E_6 Y_6

\end{cases} </math>

Ограничения

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

Примечания

Шаблон:Примечания

См. также

Шаблон:Расчёт электрических цепей

  1. Ошибка цитирования Неверный тег <ref>; для сносок Newman не указан текст