Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Новая папка / Математ. модел-е РЦ с распределенными параметрами РЛ

.PDF
Скачиваний:
24
Добавлен:
10.04.2015
Размер:
1.46 Mб
Скачать

Шаг 20. Организуются циклы вычислений значений информативных признаков при варьировании проводимости изоляции, координаты места обрыва рельсовой нити и осуществляется формирование массивов m K контрольного режима.

Шаг 21. Печатаются результаты исследований и фиксируется окончание исследований.

Блок – схема сформулированного алгоритма формирования математических моделей на ЭВМ представлена в приложении 3.

Вопросы и упражнения для самопроверки

1. Определить матрицу передаточного сопротивления рельсовой цепи по рис. 3.1 в нормальном режиме

Z N = U 2 N .

no I1N

2. Определить матрицу передаточной проводимости рельсовой цепи по рис. 3.1 в нормальном режиме

Yno

=

 

 

 

 

I 2 N

.

 

 

 

 

 

 

 

 

 

 

 

E

3. Определить матрицу передаточного сопротивления рельсовой линии

по рис. 3.2 в шунтовом режиме

 

 

 

 

 

 

 

 

 

 

 

ZnoS

=

 

 

 

U 2S

.

 

 

 

 

 

 

 

 

 

I1S

4. Определить матрицу передаточной проводимости рельсовой линии

по рис. 3.2 в шунтовом режиме

 

 

 

 

 

 

 

 

 

 

 

Y S

=

 

 

I2S

.

 

 

no

U1S

 

5. Определить матрицу передаточной проводимости рельсовой цепи по

рис. 3.3 в контрольном режиме

 

 

 

 

 

 

 

 

 

 

 

ZnoK

=

U2 K

.

 

 

 

 

 

 

 

 

 

E

6. Определить матрицу передаточной проводимости рельсовой линии

по рис. 3.3 в контрольном режиме

 

 

 

 

 

 

 

 

 

 

 

Y K

=

I2 K

.

 

no

 

 

 

 

 

E

 

 

 

 

 

 

61

4. МАТЕМАТИЧЕСКИЕ МОДЕЛИ РЕЛЬСОВЫХ ЦЕПЕЙ ПРИ МНОГОПОЛЮСНОЙ СХЕМЕ ЗАМЕЩЕНИЯ РЕЛЬСОВЫХ ЛИНИЙ

4.1. Схемы замещения рельсовой цепи и ее компонент

При формировании математических моделей рельсовая цепь рассматривается как каскадное соединение устройства согласования и защиты аппаратуры в начале рельсовой линии (УСН), устройства согласования и защиты аппаратуры в конце рельсовой линии (УСК) и самой рельсовой линии. На рис. 4.1 представлена схема замещения рельсовой цепи.

 

 

l

 

 

 

 

 

İ1

 

İ3

 

İн

 

 

 

 

 

 

[А]ун

&

 

&

[А]ук

 

 

U1

[A]о

U3

&

Z н

Ė

&

&

 

Uн

УСН

U2 Рельсовая линия U4

УСК

 

 

İ2

 

İ4

 

 

Рис. 4.1. Функциональная схема рельсовой цепи с многополюсной схемой замещения рельсовой линии

Устройство согласования и защиты в начале рельсовой цепи состоит из последовательно соединенных: четырехполюсника [A]сн, замещающего аппаратуру регулирования и ограничения тока источника питания; (2×n) полюсника дроссель – трансформатора питающего конца рельсовой линии [A]дп; (2×n) полюсника [A]сп, замещающего сопротивления соединительных проводов и эквивалентное сопротивление между средним выводом основной обмотки дроссель – трансформатора и рельсовой линией. На рис. 4.2 представлена схема замещения устройства согласования и защиты аппаратуры в начале рельсовой цепи.

 

Io

İ1

Е

[А]сн

Uo [A]дп

 

 

İ2

 

İ1

 

&

[A]сп

U1

&

 

U 2

 

İ2

 

УСН

Рис. 4.2. Схема замещения устройства согласования аппаратуры в начале рельсовой линии

62

Устройство согласования рельсовой линии с нагрузкой (УСК) состоит из каскадно – соединенных: (2×n) полюсника, замещающего дроссель – трансформатор релейного конца рельсовой линии [A]др, четырехполюсника [A]ск, замещающего аппаратуру согласования и защиты нагрузки, в частном случае защитный блок фильтров (ЗБФ). На рис. 4.3 представлена схема замещения устройства согласования и защиты нагрузки в конце рельсовой линии.

İ3

 

İ3

 

İн

 

&

 

 

 

 

 

U 3

[A ]сп

[А]др

[А]ск

U н

Zн

&

 

 

 

&

 

U 4

 

 

 

 

 

İ4

 

İ4

 

 

 

 

 

 

УСК

 

 

Рис. 4.3. Схема замещения устройства согласования и защиты нагрузки в конце рельсовой линии

4.2. Схемы замещения и параметры рельсовых линий

Удельные значения первичных параметров рельсовых линий (r, L, C, g) и соответственно ее вторичных параметров: волновое сопротивление – Zв и коэффициент распространения γ – зависят от пространственной координаты РЛ, вследствие неоднородности качества и состояния шпал, высоты балластного слоя, наличия междупутных перемычек, применяемых для канализации тягового тока, отсасывающих фидеров тягового тока заземления опор контактной сети и других факторов [4]. При этом продольными параметрами являются удельное сопротивление (r) и индуктивность (L) рельса, а в качестве поперечных параметров выступают проводимость между рельсами (g) и емкость(C) между ними, соответственно.

Расчет электрических параметров РЛ основывается на ее представлении в виде двухпроводной или трехпроводной электрической линии с распределенными параметрами.

Аналитические методы расчета, используемые в работе [4], с двухпроводной схемой замещения рельсовой линии дают удовлетворительные результаты в нормальном и шунтовом режимах, но не могут обеспечить требуемой точности в контрольном режиме, поскольку сопротивление земли как провода принимается равным нулю.

63

Для учета сопротивления земляного тракта рельсовая линия представлена в виде трехпроводной схемы (рис. 4.4).

1R

U1

 

I1

0

 

 

R0

 

 

U12

 

 

 

 

 

 

 

 

 

U2

I2

2

R

L

1’

M12 gk1

C01

g01

C12

g12

 

 

 

0’

 

 

 

 

C02

g02

 

 

 

L

 

gk2

 

 

 

 

 

 

2’

Рис. 4.4. Шестиполюсная трехпроводная схема замещения рельсовой линии:

М12 – взаимная индуктивность между рельсовыми нитями; gк1, gк2 – проводимости слоя рельс-накладка; g01, g02 – удельные проводимости заземления рельсовых нитей; g12 – удельная проводимость верхнего слоя балласта и шпал; С01, С02 – емкости двойного слоя: рельс – накладная – рельс; R0 – сопротивление земляного тракта

Из всех видов физико – химических процессов, обусловливающих появление электрической емкости при протекании переменного тока промышленной частоты между рельсами по железобетонным шпалам, наибольшее значение имеют поляризационная емкость и емкость двойного слоя. На электрифицированных участках к одной из рельсовых нитей пути присоединяются опоры контактной сети, и эквивалентная проводимость изоляции может оказаться выше 1,0 См/км, вследствие утечки сигнального тока в землю через заземление опор.

В нормальном режиме возникает поперечная несимметрия, вследствие различия проводимости изоляции рельсовых нитей (Y1 ≠ Y2), которая оценивается коэффициентом

kg = gk1 / gk 2 1.

(4.1)

Продольная несимметрия рельсовой линии как многополюсника проявляется в шунтовом режиме, вследствие замыкания рельсовых нитей колесными парами подвижного состава. Следует отметить, что любая несимметрия приводит к увеличению числа независимых параметров эквивалентного многополюсника.

Матрица параметров линейного пассивного или активного неавтономного шестиполюсника (в данном случае 3×2 полюсника) имеет размерность равную четырем (рис. 4.5).

64

 

 

 

 

 

I1

 

 

 

 

 

l

 

 

 

 

 

 

I3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3×2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U1

 

 

 

 

 

U3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

четырехполюсника

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U2

 

 

U4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I2

 

 

 

 

 

 

 

 

 

 

 

I4

 

 

Рис. 4.5. Шестиполюсная схема замещения участка рельсовой линии

 

Это

означает,

 

 

что

из

восьми

 

 

активных

параметров

&

&

&

&

&

&

 

&

 

&

, ) четыре любых параметра (вектор воздействия)

(U1

, I1

, I2

,U2

,U3

, I3

, I4 ,U4

можно задавать произвольно, а остальные четыре (вектор реакции) будут функционально зависеть от заданных. Общее число возможных видов систем параметров для рассматриваемого проходного шестиполюсника равно

C84 = 70 . Каждую из этих систем можно записать в нескольких формах, в зависимости от выбранного порядка внутри групп входных и выходных величин многополюсника. Кроме того, возможны различные способы выбора направлений токов. В рассматриваемом случае удобно выбрать направления токов I&1 и I&2 как втекающих в 2×4-полюсник, а токов I&3 и I&4 как вытекающих из него. С целью получения симметричной структуры матрицы выбран следующий вид записи системы линейных уравнений, определяющих коэффициенты матрицы А – параметров:

 

&

&

&

 

&

 

&

,

 

U

= a U

+ a I

 

+ a I

 

+ a U

 

 

& 1

11& 3

12 &3

13 &4

 

14 & 4

 

 

I1 = a21U3 + a22 I3 + a23I4 + a24U4

,

(4.2)

 

&

&

&

 

&

 

&

,

I2

= a31U3

+ a32 I3

+ a33 I4

+ a34U4

 

 

&

&

&

&

&

 

 

U2

= a41U3

+ a42 I

3

+ a43 I

4

+ a44U4.

 

Система А – параметров по сравнению с другими системами имеет ряд преимуществ. Во – первых, при каскадном соединении нескольких многополюсников расчет параметров эквивалентного многополюсника сводится к выполнению операции матричного умножения. Во – вторых, формулы для определения входных сопротивлений нагруженного многополюсника имеют простую структуру и трудоемкость вычислений по этим формулам относительно невелика.

Эквивалентным шестиполюсником может быть представлен как дискретный участок небольшой длины, схема замещения которого содержит элементы со сосредоточенными параметрами, так и участок рельсовой линии

65

произвольной длины, который должен рассматриваться как длинная электрическая линия с распределенными параметрами.

Для определения А – параметров шестиполюсника рельсовая линия представлена в виде большого числа элементарных участков длиной dx каждый, в пределах которого схема замещения может быть представлена набором типовых дискретных полиномов (рис. 4.6).

dx

i1+di1

Z1 dx i1

eM1=-ZM dx i2

i1

 

 

 

+

 

 

 

 

 

 

 

 

 

u1+du1

 

 

g1dx

u1

 

 

g12dx

 

 

u2+du2

 

 

g2dx

u2

 

 

+

 

 

i2+di2

Z2 dx i2

e =-ZM dx i1

i2

 

 

M

 

 

 

 

2

 

 

x+dx

x

Рис. 4.6. Трехпроводная схема замещения элементарного участка рельсовой линии

Для удобочитаемости формул в этом разделе условно будем обозначать комплексы пассивных и активных параметров в упрощенном виде.

На основании 2-го закона Кирхгофа (с учетом того, что ось x

направлена влево) для участка пути длиной dx получаем уравнение

 

z1dx i1 + u1 (u1 + du1 )= eM1 ,

 

 

 

 

учитывая, что eM1 = −zM dx i2 , получаем

 

 

du1

 

 

= z i

+ z

M

i

2

.

 

 

 

 

 

 

 

(4.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании 1 – го закона Кирхгофа получаем уравнение

 

 

di1

 

 

= (g

1

 

+ g

12

)u

 

 

 

g

12

u

2

.

(4.4)

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично получаем еще два уравнения

 

 

du2

 

 

= z

 

i

+ z

M

i ,

 

 

 

 

 

 

(4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

2 2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

di2

 

= (g

2

+ g

 

)u

2

g

12

u .

(4.6)

 

 

 

 

dx

 

 

 

 

 

12

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

66

Запишем соотношения (4.3) – (4.6) в виде системы обыкновенных дифференциальных уравнений

du

 

 

 

 

 

 

 

 

 

 

1

 

 

 

= z1i1 + zM i2 ,

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

di1

 

 

= (g + g )u g u ,

 

 

 

 

 

 

1

12

1

12

2

 

dx

 

(4.7)

 

 

 

 

 

 

 

 

 

di2

 

 

= −g12u1 + (g2 + g12 )u2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

du2

 

 

 

= zM i1 + z2i2.

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получение характеристического уравнения упрощается, если

выполнить формальную

 

замену вида

du

 

= p

u и

di

= p i , тогда

 

dx

dx

 

 

 

 

 

 

 

 

 

 

характеристическая матрица системы имеет вид

 

 

 

 

 

 

 

 

p

z1

zM

 

 

0

 

 

 

 

 

g

 

+ g

p

0

g

 

 

 

 

(4.8)

 

1

12

 

 

 

 

 

12

 

 

 

 

g12

0 p

g1 + g12

 

 

 

 

 

 

 

 

 

0

zM

z2

 

 

p

 

 

 

 

Собственные значения этой матрицы являются корнями характеристического уравнения. Их можно найти, приравняв определитель матрицы (4.8) к нулю. Характеристическое уравнение имеет вид

p4 [z (g + g

12

)+ z

2

(g

2

+ g

12

)2z

M

g

12

]p2

+

+ (z z

 

1 1

 

 

 

 

 

 

 

 

 

 

 

2

z2

)(g g

2

+ g g

12

+ g

 

g

)= 0.

 

 

 

1

M

 

1

 

 

 

1

 

 

2 12

 

 

 

 

 

 

Решив это биквадратное уравнение, получим выражение для определения коэффициентов распространения сигнала в рельсовой линии

γ k = ± a2 ± a22 a0 ; k=1…4;

причем

a0 = (z1z2 zM2 )(g1g2 + g1g12 + g2 g12 );

a2 = 12 [g12 (z1 + z2 2zM )+ g1z1 + g2 z2 ].

67

4.3. Матрица [A]ON – параметров несимметричной трехпроводной

рельсовой линии в нормальном режиме

Класс образов нормального режима характеризуется свободным и исправным состоянием рельсовых линий. Схема замещения рельсовой линии в виде трехпроводной линии соответствует рис. 4.5.

Общее решение для напряжения U&1 (x) (рис. 4.6) для классов образов нормального режима описывается выражением

 

 

 

 

 

 

 

u1(x) = A1chγ1x + A2 shγ1x + A3chγ2 x + A4 shγ2 x ,

 

(4.9)

где γ

1

=

a

a

2

a

;

γ

2

=

 

 

 

a

 

+ a2

a

,

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

0

 

 

 

 

 

 

 

 

 

2

 

2

 

0

 

 

 

 

 

 

 

 

 

 

A1, A2 , A3, A4

произвольные постоянные,

зависящие от граничных условий

задачи.

 

Для

 

 

каждой

 

 

корректно

 

поставленной

 

задачи

величины

&

 

&

(x),

&

 

 

 

 

 

 

 

 

 

 

выразить через значения этих постоянных и

U 2 (x), I1

I2 (x) можно

 

некоторые

дополнительные

 

 

 

коэффициенты

 

h1, h2 и

квадратную

матрицу

[yij ]второго порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2 (x)= h1(A1chγ1x + A2shγ 1x)+ h2 (A3chγ 2 x + A4 shγ 2 x),

(4.10)

 

 

h =

γ

2

z

(g + g

12

)

+ z

M

g

12 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h =

 

γ

2

z

 

(g + g

 

 

)

+ z

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

1

 

 

12

 

 

 

 

M

12 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i1(x)= y11(A1shγ 1x + A2chγ 1x)+ y12 (A3shγ 2 x + A4chγ 2 x),

 

(4.11)

 

 

i2 (x) = y21(A1shγ 1 x + A2chγ 1x)+ y22 (A3shγ 2 x + A4chγ 2 x),

(4.12)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[yij ]=

 

 

 

 

1

 

 

 

 

 

γ1(z2 h1zM )

γ2

(z2 h2 zM )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

(h z z

 

)

 

 

(h z

 

 

 

)

.

 

 

 

 

 

 

 

 

 

z z

2

 

z

2

 

 

 

1

M

γ

2

1

z

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

M

 

 

1 1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

&

&

 

&

 

 

&

 

 

 

 

 

 

 

 

 

Определим значение U1 , I1 , U2 ,

 

I2 , через параметры рельсовой линии.

Для этого продифференцируем по (x) обе части выражения (4.3)

 

 

d 2u

 

 

di

 

 

 

 

di

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

= z

1

 

+ z

 

 

 

2

.

 

 

 

 

 

 

 

 

 

 

 

 

(4.13)

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 dx

 

M dx

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь подставим в правую часть (4.13) соотношения (4.4) и (4.6)

 

 

d 2u

 

 

[(g

 

 

 

)u g

 

 

 

]+ z

 

[(g

 

 

)u

 

 

u ]=

 

 

1

 

= z

+ g

2

12

u

2

M

2

+ g

2

g

 

 

 

 

 

dx1

1

1

 

 

1

 

 

 

 

12

 

 

12 1

(4.14)

 

 

u1[z1(g1 + g12 )z M g12 ]+ u2 [zM (g2 + g12 )z1 g12 ].

 

 

 

 

68

Выразим из (4.22) переменную u2

 

 

 

 

 

 

u

 

=

 

1

 

 

d 2u1

u [z (g

+ g

 

)z

g

 

] .

2

 

(g2 + g12 )z

 

 

12

12

 

 

zM

 

 

1 1 1

 

M

 

 

 

 

 

1 g12 dx1

 

 

 

 

 

 

 

Дважды

продифференцируем

выражение

(4.9) по

x

и

(4.15)

подставим в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

правую часть выражения для

 

 

 

1

 

 

и

u1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

2

=

 

 

 

 

1

 

 

 

 

 

{γ2 A chγ x

+ γ2 A shγ

x + γ2 A chγ

2

x

+ γ2 A shγ

2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

zM

(g2 + g12 )z1g12

 

 

1

1

 

1

 

 

 

 

1

2

 

1

 

2

 

3

 

 

 

 

 

 

2

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[A1chγ1x + A2shγ1x + A3chγ2x + A4shγ2 x][z1(g1 + g12 )zM g12 ]}=

 

 

 

 

 

 

 

 

 

γ2

 

+ z

g

 

z (g

 

+ g

)

A chγ x

 

 

 

γ2

+ z

 

g

z

(g

2

+ g

 

)

 

A shγ

 

x +

=

1

 

 

M 12

1 1

 

 

12

 

 

+

1

 

M 12

1

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

1

 

1

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ2

 

+ z

g

 

z (g

 

+ g

)

A chγ

 

x +

 

γ2

+ z

g

z

 

(g

+ g

 

)

A shγ

 

 

 

x =

+

2

 

 

M 12

1 1

 

 

12

 

 

2

 

2

 

 

M 12

1

 

1

 

 

12

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

3

 

 

 

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= F(γ1 )(A1chγ1x + A2shγ1x)+ F(γ2 )(A3chγ2x + A4shγ2x),

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.16)

где F(γ) =

γ

2 + z

M

g

12

z

(g

1

+ g

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

12

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь найдем формулы

 

для

токов

i1

 

и

 

i2.

 

Продифференцируем

выражение (4.4) по x и подставим в левую часть соотношения (4.3). Для удобства обозначим

в =

du1

 

= γ

1

A shγ

1

x + γ

1

A chγ

1

x + γ

2

A shγ

2

x + γ

2

A chγ

2

 

 

 

1

 

 

dx

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

3

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в

2

=

du2

 

= h γ

1

A shγ

1

x + h γ

1

A chγ

1

x + h γ

2

A shγ

2

x + h γ

 

 

 

 

 

dx

 

 

 

 

1

 

1

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

2

 

3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h = F(γ

 

 

) =

γ2

+ z

M

g

12

z (g + g

12

)

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

1 1

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

zM (g2 + g12 )z1g12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h = F(γ

 

 

) =

 

γ2

+ z

M

g

12

z

 

(g + g

12

)

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

1

 

1

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

zM

(g2 + g12 )z1g12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x,

2 A4chγ2 x,

В результате получим систему линейных алгебраических уравнений второго порядка относительно искомых переменных i1 и i2

z1i1 + zM i2 = в1,zM i1 + z2i2 = в2 ,

Решив эту систему, получим

69

i1 = в1z2 в2 zM z1z2 zM2

= z2 γ1 zM h1γ1 (A1shγ1x + A2chγ1x)+ z1z2 zM2

 

 

 

 

 

 

 

 

 

 

 

 

 

+

z1γ2 zM h2γ2

(A shγ

2

x + A chγ

2

x)

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z1z2 zM

 

 

 

3

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= y11(A1shγ1x + A2chγ1x)+ y12 (A3shγ1x + A4chγ2 x),

где y

 

= γ

 

z2 h1zM

;

y

 

 

= γ

 

 

z2 h2zM

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

1 z z

2

z2

 

 

 

 

12

 

 

2 z z

2

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

M

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

M

 

 

 

 

 

 

 

i

 

=

 

z1в2 zM в1

 

= γ

 

 

z1h1 zM

(A shγ

 

x + A chγ

 

x)+

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

z z

2

 

z2

 

 

1 z z

2

z2

 

 

 

1

 

 

1

2

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

M

 

 

 

 

 

1

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

+ γ

 

 

h2 z1 zM

(A shγ

 

 

x + A chγ

 

x) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

z z

2

z2

3

 

 

2

 

 

 

4

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= y21(A1shγ1x + A2chγ1x)+ y22 (A3shγ2 x + A4chγ2 x),

где

y

 

= γ

 

z1h1 zM

 

;

y

 

 

= γ

 

 

h2z1 zM .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

 

1 z z

2

z2

 

 

 

 

 

22

 

 

2 z z

2

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

M

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

M

 

 

 

 

 

Возвращаясь к обозначениям токов и напряжений в схеме, изображенной на рис. 4.5, а также учитывая направление оси x, указанное на рис. 4.6, полагая в уравнениях (4.9), (4.10), (4.11) и (4.12) x=0, получаем

&

 

= A

+ A ,

U

 

&

3

1

3

U 4 = h1 A1 + h2 A3 ,

&

 

= y11 A2 + y12 A4 ,

I3

 

&

 

= y21 A2 + y22 A4 .

I 4

Выразим параметры A1, A2 , A3, A4 через U& 3 ,U& 4 , I&3 , I&4

 

&

&

 

&

&

 

&

&

 

&

&

 

A

=

U 3h2

U 4

; A =

I3 y22

I4 y12

; A

=

U4

U3h1

; A =

I4 y11

I3 y21

. (4.17)

1

 

h2 h1

2

 

d

3

 

h2 h1

4

 

d

 

 

 

 

 

 

 

 

 

 

d = y11 y22 y21 y12 .

 

 

 

 

 

 

 

 

 

 

 

 

Полагая в уравнениях (4.9), (4.10), (4.11) и (4.12) x=l (l – длина линии),

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I&1 = y11S1 A1 + y11C1 A2 + y12 S2 A3 + y12C2 A4 ,

 

 

 

 

 

 

 

&

= y21S1 A1 + y21C1 A2

+ y22 S2 A3 + y22C2 A4

,

 

 

 

 

 

 

I2

 

(4.18)

 

 

 

&

= C1 A1 + S1 A2 + C2 A3 + S2 A4 ,

 

 

 

 

 

 

 

U1

 

 

 

 

 

 

 

 

 

&

= h1C1 A1 + h1S1 A2 + h2C2 A3 + h2C2 A4 ,

 

 

 

 

 

 

 

U 2

 

 

 

 

где S1 = sh(γ 1l );

С1 = сh(γ 1l );

S2 = sh(γ 2l); С2 = сh(γ 2l);

 

 

 

 

70