
- •Применение mathcad в теории электрических цепей
- •Введение
- •1 Расчёт сложных разветвленных электрических цепей
- •2 Расчет линейных электрических цепей синусоидального тока
- •3 Построение графиков функций в Mathcad
- •3.1 Построение двумерных графиков функций, заданных аналитическими выражениями
- •3.2 Построение в одном шаблоне графиков нескольких функций
- •3.3 Построение графиков кусочно-непрерывных функций
- •3.4 Форматирование графиков
- •3.5 Построение двумерных графиков функций, заданных таблично
- •3.6 Примеры построения графиков
- •3.7 Построение потенциальных и топографических диаграмм
- •4 Расчет переходных процессов классическим методом
- •5 Расчет переходных процессов операторным методом
- •6 Расчет переходных процессов с помощью интеграла Дюамеля
- •7 Расчет переходных процессов спектральным методом
- •Список литературы
- •050013, Алматы, ул.Байтурсынова 126 алматинский университет энергетики и связи Кафедра теоретических основ электротехники
- •Применение MathCad в теории электрических цепей Учебное пособие
5 Расчет переходных процессов операторным методом
Исходные данные:
Е = 125 В
R1= 50 Ом
R2= 200 Ом
R4= 250 Ом
L = 0,01 Гн
С = 5 мкФ
Определить ток iL в переходном режиме.
Для составления эквивалентной операторной схемы определим независимые начальные условия iL(0) и uC(0) по законам коммутации:
iL(0)
=
iL(0-)
=
= 0,25 А;
uC(0)
= uC(0-)
=
=
50 В.
Составим эквивалентную операторную схему:
Рисунок 5.2
Составим уравнения методом контурных токов относительно изображений контурных токов по схеме рисунок 5.2.
(5.1)
В среде MathCad решить систему уравнений (5.1) можно различными способами.
Способ 1. Набрать систему уравнений в шаблоне матрицы, с числом строк - 2 и числом столбцов - 1. Между левой и правой частями уравнений поставить логический знак равенства. Применить к системе уравнений, записанных в виде матрицы, операторы «solve», «explicit», «collect» которые находятся на панели «Символика» и используются для решения системы уравнений в символьном виде относительно заданной переменной.
Способ 2. Решить систему уравнений в матричной форме путем обращения матрицы коэффициентов. Где Z – квадратная матрица контурных сопротивлений; Ek – матрица-столбец контурных ЭДС; Iкк (p) – матрица-столбец контурных токов.
Таким образом:
. (5.2)
Для определения корней характеристического уравнения приравняем знаменатель выражения (5.2) к нулю.
Следовательно один корень р=0, остальные корни определяются из выражения
(5.3)
Решение квадратного уравнения (5.3) дает корни p1 = - 3000 +j4000 c-1;
p2 = - 3000 - j4000 c-1
Для определения оригинала используем следующий вид теоремы разложения:
, (5.4)
,
,
.
Подставив в полученные выражения p1 = - 3000 +j4000 c-1, получим:
.
Подставив рассчитанные значения в теорему разложения (5.4), получим:
(5.5)
Определение корней характеристического уравнения (5.3) и функций входящих в теорему разложения (5.4), (5.5) можно произвести в среде MathСad изложенными выше способами.
При
построении графика в среде MathCad
рассчитаем время окончания переходного
процесса t
= 5τ = 5
=
5
=0,0017
c
для определения наибольшего значения
по оси х.
Для определения порядка значений по
оси у
рассчитаем значения iL(0)
и iLсв(0)
с помощью оператора «substitute»
при t=0.
Подробно методика построения графиков
приведена в разделе 3 данного пособия.
Рисунок 5.3
6 Расчет переходных процессов с помощью интеграла Дюамеля
На вход цепи (см. рисунок 6.1) с заданными параметрами R1= R3=40 Oм, R2= 20 Oм, L = 10 мГн, U = 20 В, tи = 4 мс подается импульс напряжения u(t) (см. рисунок 6.2). Определить ток i2(t) в переходном режиме.
Рисунок 6.1 Рисунок 6.2
Этап1.
Для определения переходной проводимости
рассчитываем
переходной процесс при подключении
цепи при нулевых начальных условиях
i2(0)=0
к источнику единичного постоянного
напряжения (см. рисунок 6.3) и определяем
ток i2(t).
Рисунок 6.3
Ток i2 в переходном режиме можно представить в виде i2= i2пр+ i2св.
Определим принужденное значение тока i2:
,
А.
,
А.
Для записи свободной составляющей найдем корни характеристического уравнения, составленного для послекоммутационной схемы относительно ветви с индуктивностью:
Таким
образом,
При t=0: 0=0,0125Е+А, откуда, А=-0,0125Е, а значит
.
Следовательно,
переходная проводимость
равна:
Этап 2. Определение тока i2(t) в переходном режиме с помощью интеграла Дюамеля при подключении цепи на импульс напряжения u(t):
Рассчитаем
i2(t)
в интервале времени
,
не включая скачка.
,
,
Проверить правильность расчета интеграла Дюамеля в среде MathCad в первом интервале можно с помощью операторов expand и collect.
Рассчитаем
i2(t)
в интервале времени
.
Проверить правильность расчета интеграла Дюамеля в среде MathСad во втором интервале можно с помощью операторов expand и collect.
График i2(t) представляет собой кусочно-непрерывную функцию, заданную на первом интервале выражением (6.1) и на втором интервале выражением (6.2). График строится в интервале времени от 0 до 2tu. Для определения порядка значений по оси Y рассчитаем значения i2(tu) для обоих интервалов с помощью оператора «substitute» при t= tu= 0.004c, где i21(t)- значение тока i2(t) в первом интервале, i22(t)- значение тока i2(t) во втором интервале. Подробно методика построения графиков изложена в разделе 3 данного пособия.
Рисунок 6.4
Рисунок 6.5