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

3.2. Элементарная схема.
Схема не содержит источников напряжения, поэтому уравнение, описывающее
отдельные ветви можно привести в узловой форме:
ib + Jb = Ybb ∙ Ub, (3.2.1)
где матрица проводимостей Ybb при расчете в частотной области:
в
в 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
1 G0
2 G0
3 G0
4 G0
5 G0
6 G0
Ybb= 7 G0
8 G0
9 G0
10 G0
11 Ki Ki YC
12 Ki Ki YC
13 Ki Ki YC
14 Ki Ki YC
15 Ki Ki YC
здесь использованы следующие обозначения
Для расчета в частотной области Y c = jω 1
G0= 2δ'
μ 0
4ah 2
ρ
2 ah
, K i =
ρ ⋅V .
(3.2.2)
, где ω1 − угловая ча-
стота питания обмоток индуктора.
Для расчета во временнóй области методом Эйлера 1−го порядка емкость
заменяется разностным аналогом с уравнением в узловой форме с
источником тока:
i n1 C b
Δt U n = C b
Δt U n1 .
(3.2.3)
Таким образом, в матрице Ybb емкостная проводимость вычисляется
по формуле:
Y c = 4ah2
Δt⋅ρ
Дополнительно к независимому источнику тока Jb, действующему на
емкостной ветви надо добавить величину 4ah2
Δt⋅ρ ⋅U n , где Un − напряжение,
вычисленное на предыдущем расчетном шаге.
3.3. Уравнения подсхемы.
Сформируем уравнение в контурном базисе. Для этого потребуется матрица
сопротивлений ветвей Zbb:
Zbb = (Ybb) -1 .
122