?????????????? ?????? ? ????????? ?????????? ??????? ?????

3.5. Алгоритм расчета электромагнитного поля во временной
области.
Приведем алгоритм расчета одномерного поля в системе SciLab.
Расчет мгновенных процессов производится в цикле по времени с шагом
приращения времени dt. Порядок расчета внутри цикла следующий.
• вычисление фазных токов Ia=sin(2π·f·t), Ib=sin(2π·f·t-2π/3), Iс=sin(2π·f·t+2π/
3),
• вычисление пазовых токов Ip=Cпф·Iф,
• вычисление источников тока, действующих на емкостных ветвях:
Jv=Uo·Cb + Сдиф·Ip
Источники тока в ветвях схемы – это производная от пазовых токов:
J v = ∂ I p
∂ x ≈ I p−I p− ,
2h
где индексы + и – соответствуют обозначениям в формуле (3.1.4). Вычисление
производной можно записать в матричной форме, связывающей
источники тока в ветвях Jv с пазовыми токами Ip: Jv = 1/2h ∙ Сдиф ∙
Ip, где Сдиф – «дифференцирующая» матрица размером v x p, имеющая
следующий вид:
пазы
ветви … p-1 p p+1 p+2 …
… … … … …
Сдиф = v -1 +1
v+1 -1 +1
… … … … …
• расчет источников тока Jr на узловых и контурных координатах подсхем
по формуле: Jr=At·Jv,
• из вектора Jr выделяются узловые источники тока Jо,
• расчет узловых напряжений отдельных подсхем Uo=Zor·Jo, Из этого
вектора узловых напряжений выделяются узловые напряжения на
внешних узлах Uoo. Эти узловые напряжения образуют контурные
э.д.с. цепи пересечений es,
• расчет контурных э.д.с. цепи пересечений es= - Cos·Uoo,
• расчет контурных токов цепи пересечений Is=Yss·es,
• контурные токи цепи пересечений – это результат соединения подсхем.
Контурные токи распределяются по подсхемам, вызывая приращения
узловых токов Jo, найденных ранее. Эти приращения dIo находятся
по формуле: dIo = Cos·Is,
• увеличение узловых потенциалов в подсхемах, вызванное их соединением
находится по формуле: dUo=Zor·dIo,
131