Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Конспект для итиутс.pdf
Скачиваний:
115
Добавлен:
11.05.2015
Размер:
766 Кб
Скачать
☆

yПЗ = -

СЗ

; yПИ = -

СПИ

-

 

1

; yПП

=

ССП +СЗП +СИП

+

1

+

1

;

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

UИП

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

UСП

UИП

 

 

 

ì

β

éU

 

- А(ϕ

С

-ϕ

И

)ù

при

 

ϕ

С

 

<

 

U

НАС

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï

 

ë НАС

 

 

 

 

 

 

 

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

= í

 

при

 

ϕС

 

>

 

UНАС

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ìβ é

А(ϕ

 

-ϕ

 

 

)ù

 

при

 

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

y2

ï

ë

 

С

 

 

 

 

И

û

 

 

 

 

 

 

 

 

 

 

С

 

 

= í

 

 

при

 

ϕС

 

>

 

UНАС

 

;

 

 

 

 

 

 

 

 

 

ïβU

 

 

 

 

 

 

 

 

 

 

 

 

 

 

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

é

 

 

 

АК (ϕС -ϕИ )

 

 

 

 

ù

 

 

ï-β êUНАС -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 А

ϕ

 

 

 

-ϕ

П )

+ 2ϕ

 

 

 

 

 

 

ï

ê

 

И

F

ú

 

y3

ï

ë

 

 

(

 

 

 

 

 

 

û

 

= í

é

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ù

 

ï

 

 

 

 

АКUНАС

 

 

 

 

 

 

ï-β êUНАС -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 А(ϕИ -ϕП ) + 2ϕF

 

 

ï

ê

 

ú

 

î

ë

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

û

y4 = − y1 − y2 − y3.

< UНАС ;

при ϕС < UНАС ;

при ϕС ³ UНАС ;

3.5. Модели распределенных RC- и RLC-структур

Математическим описанием элементов, параметры которых имеют распреде- ленный характер, служат дифференциальные уравнения в частных производных.

Это обстоятельство создает определенные трудности при построении моделей распределенных элементов, так как используемые в программах АПЭЦ матема- тические модели базируются на методах решения обыкновенных дифферен- циальных уравнений. В связи с этим опишем алгоритмы представления моделей распределенных элементов в виде РИ - моделей многополюсников, в которых локализуется решение уравнений в частных производных.

Прежде чем перейти к способам представления моделей, рассмотрим математи- ческое описание неоднородных структур RC-и RLC-типа.

3.5.1. Исходные уравнения элементов с распределенными параметрами

Различают два основных типа распределенных RС-элементов: R-C-NR- и C- R-NC- структуры (рис.3.5). Учитывая свойство дуальности этих структур, в даль- нейшем будем рассматривать только структуры типа R-C-NR. Уравнения этой структуры записываются следующим образом:

∂u

= -r(x)éi

(x,t) - Ni

 

(x,t)ù;

(3.14)

x

 

ë R

 

NR

û

 

∂iR = -

∂iNR

= -C(x)

u

,

 

(3.15)

 

t

 

x

 

x

 

 

 

 

где r(x) и С(х) - соответственно погонные сопротивления и емкость, приходящиеся на единицу длины структуры l.

53

r(x) x

r(x+ x) x

Ia

R

Ic

 

a

С

с

Ib

 

Id

b

nR

d

 

 

uC(x,t)

С R

NС

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ir(x,t)

 

 

 

 

 

 

 

 

 

ir(x+

x,t)

С(x)

 

x

 

 

 

 

 

 

С(x+

x)

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i(x,t)

 

 

 

u(x,t)

 

 

 

 

 

 

 

 

 

 

 

u(x+ x,t)

 

 

 

 

nr(x)

 

x

 

 

 

 

 

 

 

 

nr(x+

x) x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

 

 

 

inr(x,t)

 

 

 

 

 

 

 

 

inr(x+

x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iC(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iC(x+

 

x,t)

 

 

С(x)

 

x

 

 

 

 

 

 

uC(x+

x,t)

 

 

 

 

 

 

С(x+

x) x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ir(x,t)

 

 

 

ir(x+ x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r(x)

x

 

 

 

r(x+

x) x

 

 

 

 

 

 

 

 

 

 

 

 

 

NС(x)

 

x

 

 

 

 

 

 

uNC(x+

x,t)

 

 

 

 

 

 

NС(x+ x) x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iNC(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iNC(x+ x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Рис. 3.5. Распределенные RC-структуры: а - условное обозначение и эквивалентная схема R-C-NR-структуры; б - условное обозначение и эквивалентная схема С-R-NC-структуры

Дифференцируя (3.14) по х и подставляя результат в чить основное уравнение R-C-NR- структуры

∂2u −σ (x)

∂u − μ(x)

∂u = 0

∂x2

∂x

∂t

с граничными и начальными условиями

(3.15) нетрудно полу-

(3.16 )

u(0,t) = Ф(i),u(l,i) = Ψ(t),u(x,0) = u0 (x).

(3.17 )

Коэффициенты уравнения σ(х) и μ(x), характеризующие неоднородность струк-

туры, определяются через погонные параметры:

 

σ (x) =

1

∂r(x); μ(x) = (N +1)r(x)C(x).

 

 

 

 

r(x) ∂x

 

Уравнение (3.16) является уравнением параболического типа, решение кото- рого совместно с заданными граничными и начальными условиями представляет со- бой краевую задачу первого рода. Дуальное математическое описание структуры С- R-NC–типа рассмотрено ранее.

Уравнение, описывающее RLC – структуру – неоднородную двухпровод-

54

ную линию с потерями (рис. 3.6, а) через первичные погонные параметры L(x), С(х), r(х) и g(x), можно получить из эквивалентной схемы элементарной ячейки

(рис. 3.6, б):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i1(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r(x)

 

 

 

 

 

 

 

 

 

 

Ia

i1(x,t)

 

Ic

 

 

 

 

 

 

 

 

 

 

 

x

L(x)

 

x

 

 

 

 

 

 

 

 

 

С(x) x

 

 

 

 

g(x) x

 

 

2

 

2

 

 

 

 

 

 

 

 

 

a

 

 

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r(x)

 

L(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ib b

i2(x,t)

d Id

 

 

 

 

 

u(x,t)

 

 

x

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

 

i2(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.6. Распределенная RLC–структура: а– условное обозначение; б–

 

 

 

 

 

 

 

эквивалентная схема

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

∂u = r(x) (i

−i ) +

L(x)

 

∂(i1

−i2 )

;

(3.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x

2

 

1

2

 

2

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−

∂i1

 

= ∂i2

= g(x)u +C(x) ∂u .

 

 

 

 

(3.19)

 

 

 

∂x

 

 

 

 

 

 

 

 

 

∂x

 

 

 

 

 

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

 

 

 

 

Объединяя (3.18) и (3.19), нетрудно выразить основное уравнение распре- делённой линии:

∂2u −σ (x)

∂u − μ(x)

∂2u −ε(x)

∂u −η(x)u = 0.

(3.20)

∂x2

∂x

∂t2

∂t

 

с граничными и начальными условиями (3.17). Коэффициенты уравнения (3.20) для однородной линии с потерями (R не равно 0, G не равно 0) определяются через полные параметры R, G, L, С:

σ = 0, μ = LC, ε = LC + LG,η = RG. Для неоднородной линии без потерь (r(x)=0, g(x)=0 ),

σ (x) = (1/ L(x))(dL(x)/ dx); μ(x) = L(x)C(x), ε(x) = 0,η(x) = 0.

Уравнение (3.20) представляет собой уравнение в частных производных гиперболического типа, решение которого так же, как и для RCL–структур, сводится к краевой задаче первого рода.

Учитывая, что рассматриваемые элементы с распределёнными параметрами являются линейными и инерционными, будем строить разностные (Р) модели этих элементов в форме (3.3).

Чтобы получить матричное уравнение (3.3) Р – модели элемента, необхо- димо явно выразить связь между полюсными токами и напряжениями элемента,

причем при построении модели в базисе узловых потенциалов аргументами должны быть напряжения. Полюсные токи RC– и RLC–структур, как следует из (3.14) и (3.18), можно выразить через значения производных на концах структур дu/дх|x=0

и дu/дx|x=1.

Значения этих производных в свою очередь, могут быть вычислены в резу-

55

льтате численного решения уравнений (3.16) и (3.20) и являются функциями заданных в форме (3.17) граничных условий, т.е. напряжений между полюсами. Ниже показано, что производные дu/дх|x=0 и дu/дx|x=1 на каждом временном шаге численного решения уравнений (3.16) и (3.20) приближённо выражаются в виде линейной комбинации граничных условий Ф(t) и ψ(t):

∂u

 

 

≈ u

 

 

=Tϕ(t

n+1

) − PФ(t

n+1

) +

 

 

 

 

∂x

 

x=0

x

 

x=0

 

 

 

 

 

 

 

 

 

 

∂u

 

 

≈ u

 

 

=Tϕ(t

n+1

) − KФ(t

) +

 

 

 

 

∂x

 

x=0

x

 

x=0

 

 

n+1

 

 

 

 

 

 

 

 

E(tn ,tn−1);

(3.21)

F(tn ,tn−1);

С помощью (3.21) можно получить из исходного описания распределенных структур их Р–модели.

3.5.2. Разностная модель двухпроводной линии

Рассмотрим методику составления Р – модели для однородной линии с потерями. Пусть полюсные токи линии имеют направления, показанные на рис. 3.6, а. Для фор- мирования искомой модели нам необходимо связать полюсные токи Ia, Ib, Ic, Id с напряжениями между полюсами uab, ubd, uac, ucd. Из (3.19) очевидно, что сумма токов i1(x,t) и i2(x,t) не зависит от координаты x и определяется внешними напря- жениями uac и ubd , откуда следует справедливость выражения

R

(i

−i ) + L

∂(i1 −i2 )

= u

 

+u ,

(3.24)

2

∂t

ac

1

2

2

 

bd

 

где R и L - полные сопротивления и индуктивность линии. Полагая длину линии l равной единице, получаем, что значения погонных параметров R, L, G, С равны значениям полных параметров линии.

Подставим в (3.18) и (3.24) токи Ia = i1(0,t), Ib = i2 (0,t), Ic = −i1(l,t), Id

= −i2 (l,t),

и их производные по времени:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

(I

 

− I

 

) +

L

(I

 

− I

 

) = − ∂u

 

 

; R

(I

 

− I

 

) +

L

(I

 

− I

 

) = −

∂u

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

a

 

b

2

 

a

 

b

∂x

 

x=0

2

 

 

d

 

 

c

2

 

 

d

 

 

c

 

∂x

 

x=l

(3.25 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

(I

 

− I

 

) +

L

(I

 

− I

 

) = u +u ; R (I

 

− I

 

) + L (I

 

− I

 

) = −u −u .

 

2

 

 

 

 

 

 

 

 

 

 

 

a

 

b

2

 

a

 

b

ab

 

bd

2

 

d

 

 

c

 

2

 

d

 

 

c

 

 

ac

 

bd

 

Здесь I означает производную по времени dI/dt. Предположим, что в качес-тве метода интегрирования в уравнении ММС выбран неявный метод первого порядка, тогда производные в (3.25) можно заменить следующими разностными эквива- лентами: I ≈∆I/∆t, ∂ u/∂ x ≈∆u/∆x. При этих условиях, учитывая (3.21), получаем систему уравнений, связывающую полюсные токи и напряжения линии, решение которой относительно вектора полюсных токов дает искомую Ρ-модель однородной

двухпроводной линии в виде (3.26), где Ia*, Ib*, Ic*, Id* полюсные токи в момент

времени t-∆t . Данная Ρ–модель характеризует состояние двухпроводной линии на некотором (n+1)-м шаге численного интегрирования уравнения (3.20).

56