Материал: 5_Лекции_по_ТОЭ_20-и тд

Таблица 1. Кусочно-линейная аппроксимация ВАХ нелинейного резистора

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

В качестве примера определим напряжение  в цепи на рис. 2, в которой     . ВАХ нелинейного резистора  приведена на рис. 3, где .

Решение

            1. В соответствии с заданной ВАХ нелинейный резистор на участке 1-2 заменяем линейным резистором с сопротивлением

,

на участке 2-3-источником тока с током  и на участке 4-1-источником тока с током .

            2. На основании данной эквивалентной замены для тока на участке 1-2 ВАХ можно записать:

откуда

При движении изображающей точки по участку 2-3 ВАХ имеем

,

при движении по участку 1-4 ВАХ-

.

            3. Определяем интервалы движения изображающей точки по отдельным участкам ВАХ. Для точки излома 1 на основании (1) справедливо уравнение

или

.

            Отсюда получаем два значения мгновенной фазы питающего напряжения на одном периоде, соответствующих точке 1: . Первое значение определяет переход изображающей точки с участка 4-1 на участок 1-2, второе – с участка 2-1 на участок 1-4.

            Аналогично записываем для точки 2 излома ВАХ

или

откуда  (значение, соответствующее переходу с участка 1-2 на участок 2-3) и  (значение, соответствующее переходу с участка 3-2 на участок 2-1).

            Таким образом, получаем для одного периода питающего напряжения

;

;

;

;

               .

            В соответствии с периодичностью синусоидальной функции данные решения повторяются через 360°n.

            На рис. 4 представлен график зависимости искомой величины.

Метод гармонического баланса

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

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

            Методика расчета нелинейной цепи данным способом включает в себя в общем случае следующие основные этапы:

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

            2. Выбирается выражение аналитической аппроксимации заданной нелинейности.

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

            4. Осуществляется подстановка функций, определенных в пунктах 2 и 3, в уравнения состояния с последующей реализацией необходимых тригонометрических преобразований для выделения синусных и косинусных составляющих гармоник.

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

            6. Осуществляется решение (в общем случае численными методами на ЭВМ) полученной системы уравнений относительно  и .

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

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

Для решения определяем аналитическое выражение характеристики  для первых гармоник:

откуда

После подстановки выражения тока и соотношения (2)  в уравнение состояния цепи

получаем

или

            На основании последнего получаем систему уравнений

из которых находим искомые параметры и .