- •1. Процесс проектирования электрических цепей
- •1.1. Уровни абстрагирования и аспекты описаний проектируемых объектов
- •1.2. Операции, процедуры и этапы проектирования
- •1.3. Классификация параметров проектируемых объектов
- •1.4. Классификация проектных процедур
- •1.5. Виды обеспечения в системах автоматизированного проектирования
- •2.3. Параметры и переменные электрических схем
- •2.4. Выбор базисных переменных
- •2. Основы автоматизированного проектирования электрических цепей
- •2.1. Основные характеристики программ
- •2.2. Возможности автоматизации схемотехнического проектирования электронных схем
- •2.5. Методы расчета электрических цепей на ЭВМ
- •2.5.1. Матрица сопротивлений
- •2.5.2. Матрица проводимостей
- •2.5.3. Матрица ЭДС источников
- •2.6. Матричные уравнения линейных электрических цепей
- •2.6.1. Матричные уравнения контурных токов
- •2.6.2. Матричные уравнения узловых потенциалов
- •2.6.3. Уравнения переменных состояния
- •3. Математические модели элементов схем
- •3.2. Модели двухполюсных элементов
- •3.3. Модели биполярных транзисторов
- •3.3.1. Обобщенная электрическая модель
- •3.3.2. Полная передаточная модель
- •3.3.3. Модель в матричной форме
- •3.4. Модель МДП - транзистора
- •3.5. Модели распределенных RC- и RLC-структур
- •3.5.1. Исходные уравнения элементов с распределенными параметрами
- •3.5.2. Разностная модель двухпроводной линии
- •3.5.3. Определение коэффициентов моделей
- •3.6. Модели элементов и функциональных узлов в виде эквивалентных схем
- •3.6.1. Модели зависимых источников
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
