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

книги из ГПНТБ / Богатырев Ю.К. Импульсные устройства с нелинейными распределенными параметрами

.pdf
Скачиваний:
5
Добавлен:
24.10.2023
Размер:
11.03 Mб
Скачать

Качественное исследование структуры ударных волн, проведенное выше, не позволяет непосредствен­ но получить необходимые для расчета формирующих систем количественные оценки характерных парамет­ ров волны: длительности фронта, амплитуды, периода и времени затухания осцилляций за ее фронтом и т. д. Такие оценки удается извлечь из аналитического или численного решения исходных или укороченных урав­ нений.

В качестве примера аналитического исследования структуры ударной волны найдем приближенное решение нелинейного диф­ ференциального уравнения (1.26), допуская аппроксимацию ха­ рактеристики Q(u) полиномом третьей степени: Q(u )=CoUD2u3 (Со, Ь 2— постоянные, положительные коэффициенты).

Укороченное уравнение (1.26) в этом случае запишем >в виде системы *>

i£dLJ2

du/dz = у; d y / d z д

(l/2— и2) и—

где согласно (1.24) А2= ( 1D2U2/Co).

При достаточно малом значении г (таком, что г=цг) для решения полученной системы можно применить метод усреднения

по негармоническим функциям {102].

В

нулевом

приближении

(ц = 0)

имеем

 

 

 

 

du/dx—y, dy/dx = 12D2(U2и2) ufAl Co,

 

откуда

уравнение фазовых траекторий

на

плоскости у, и будет

 

у 2= 6D (2U W —u'—UA),

 

где А — некоторая постоянная.

траекторий

описывает за­

При 0 < A < U 2 уравнение фазовых

мкнутые кривые типа эллипса, вложенные друг в друга. При за­ данном значении А представляющая точка совершает на плоско-

*> Поскольку качественное влияние сопротивлений R и г на структуру ударной волны одинаково, в (1.26) для простоты вы­ числений положим т н=0 и, кроме того, будем считать, что U| = 0

и Q{Ui) = 0, a U?=U .

40

сти у, и периодические движений с периодом

Т = § 1 Г = -2Л2/'

(6^>г/С0)1/2 ,

где F(k, я/2) — эллиптический интеграл

первого рода £=('1—

—Ы2212) V2, a ullzz= U 2± U ( U 2A) l/z.

 

При учете диссипации (р.=5^0) характер движения представ­

ляющей точки на фазовой плоскости существенно изменится: она

будет двигаться по спирали, соединяющей

неустойчивую особую

точку (и = у = 0 ) с устойчивой I(u=U, у = 0).

 

Для исследования этого движения будем считать параметр А

медленной функцией времени Л(т), которая

изменяется от 0 при

т— ►—оо до U1 при т—

-)- оо. Разрешая исходную систему отно­

сительно А(т), получаем уравнение

 

 

 

dA (х)

24г

f

3D,

\

 

 

Д2ро

T T “, ) ^ w +

24r

 

зд2

 

 

и*

+ Д2Ро ( -

Со

 

2и2

"О2

решение которого удается

получить,

если

заменить переменные

■их средними значениями и2

за период Т *>

т«,

- у - j" u2dx = - у - j* и2duly = u2E (k, я/2)/F (k, я/2);

 

0

 

u%

 

J

 

 

 

I

f ,

1

4

[

u\

«1 \ E ( k .

я/2)"1

T

J

3

“2

| ц2 ~

2y +

ul ) F

(ft.

« /2 ) |’

где £(&, я /2 ) — эллиптический

интеграл

второго

рода. В общем

виде это решение имеет вид

 

 

 

 

 

 

dA (т)

I- f И (х)]

где f[A(x)] — правая часть исходного уравнения, в которой заме­ нены переменные и* и4 на средние значения и* зависящие

только от А (т).

*> Отметим, что усреднение здесь ведется по периодическим функциям, отличным от гармонических, однако, так как p-C l, та­ кое усреднение вполне корректно [101].

41

Нетрудно показать, что. начиная с некоторого момента вре­ мени т>т* значения Mi«0,14£/, средние значения а2 и ГА с до­ статочной степенью точности приближаются к величине U2 и /У4. В этом случае решение для А (т) выглядит значительно -проще:

Г (

1 —

и*

24г

А (т) ^ U2 1 —

=2

ехр

 

' и 2

Д2Ро

(U0= A (т = т *). Аналогично может быть найдено приближенное

выражение А (т) для начального участка колебаний, где иi2«

« 2 t/2, и22~0.

Изменение амплитуды колебаний за фронтом ударной волны определяется как разность значений щ—«2, в частности, для

 

 

 

ui—и2« (Ul2—Ы22) /2£/= (t/2—А*(т)) V2,

 

а

время

затухания колебаний т3« Л 2р0/24л. Частота их,

начиная

с

момента -времени, где

т>т*,

слабо зависит от параметра А.

В первом приближении ее можно считать постоянной:

 

 

 

 

о)3« 2я£/ (6D2/C0) ‘/2/Д2.

 

 

Рассмотрение не годится, очевидно, в области начального

отклонения

изображающей точки от состояния равновесия и—

= (/=0.

Это

отклонение

можно

оценить промежутком

времени

Тнач. характеризующим изменение напряжения на фронте удар­ ной волны от 0,1Ы| до 0,9«ь

 

0,9Wi

 

 

da

ча 0,7(1 — 3D2U y C 0)/U (D2/C0)I/2 .

-

У (А =.: 0)

 

O.lKi

 

Кроме того, метод усреднения справедлив лишь при выпол­ нении неравенства Тнач'СТз, что всегда выполняется в реальных линиях. С другой стороны, как уже было отмечено, время запаз­ дывания фронта стационарной ударной волны на одно звено ли­ нии А должно быть много меньше некоторого характерного вре­ мени т, равного по порядку величины тНачЭто выполняется тем лучше, чем меньше параметр D2, т, е. чем меньше нелинейность характеристики Q(u). Следовательно, вышеизложенное справедли­ во для слабых ударных волн.

В большинстве практических случаев основной ин­ терес представляют сильные ударные волны. Однако при,исследовании их структуры почти всегда необхо­ димо учитывать времена релаксации нелинейных ве­

42

Рис. 1.12. Схема тороидаль­ ного ферритового сердеч­ ника.
вид

ществ, поэтому метод исследования стационарных ударных волн рассмотрим на примере формирующих линий с ферромагнитным и сегнетоэлектрическим за­ полнением.

Структура стационарных волн в линии с нелинейны­ ми элементами, описываемыми динамическими уравне­ ниями связи. Рассмотрим систему нелинейных уравне­ ний (1.5) и (1.7) предпола­ гая, что емкость звена ли­

нейна (Q(u)=C0u), а ин­ дуктивность — катушка, на­ мотанная на тороидальном ферритовом сердечнике, ди­ намика перемагничивания которого полем ударной волны Hi,=PoihZo (рис. 1.12)

описывается уравнением (1.5). В этом случае систе­ ма укороченных уравнений

1+ '

d

 

2Д2н d*4_

4пт]M

d^m

I

d2i

(1.28)

dz

i

(2k ) \

di2h

p 0

dx2

dx2 ’

 

 

k=\

 

 

 

 

 

 

 

dm

 

VT8xo(l — tr?) i,

k = l , 2,

3, .

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

где тг= гС 0/то,

m(x) =Mh/M

(Mh — проекция

вектора

намагничивания M на направление действующего по­

ля Н,

М — модуль

вектора

М).

Запаздывание А, со­

гласно

(1.24), равно

 

 

 

 

 

 

 

 

 

Ф2 — Ф, __- , 4ЩМ (т2 т ,)

(1.29)

 

 

 

 

1\

Ръ (Iг — 7,)

 

 

 

 

 

 

 

где mi, 2 — значения

т до фронта и за

фронтом удар­

ной волны

(т.

е. до и после перемагничивания ферри­

та). Рассматривая укороченные уравнения для разных

43

k, исследуем структуру стационарных ударных волн в формирующих линиях различного вида.

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

нии (1.28)

Т г= 0, найдем

его решение для k — \, для

случая, не

учитывающего

влияние пространственной

дисперсии. Это решение представляет интерес в том отношении, что позволяет исследовать и сравнить раз­ дельно влияние дисперсии линии и диссипации ферри­ та на структуру фронта ударной волны.

Укороченные уравнения для k= \ имеют вид

Сведя эту систему к одному уравнению и проинтегри­ ровав его дважды по времени с учетом начального условия i — I i= 0, m = mi при t-*— оо, получим

Подставив выражение для Д2 и положив m%= 1 (что

соответствует полному насыщению феррита при (=

44

= /2), найдем в неявном виде |(т):

In (1

m,) [1 -(-/и, -f- (1 — /п,)/// 2]<т‘ 1)12if \1— т1—

— (1 — тО ///21(т,+1)/2 /, =

1/2 8х0 (1 - т.) / 2х.

В частности,

при mi = 0 это

выражение принимает

более простой вид

 

i =

/ 2

{1+ е х р [5х0/ 2т]}-1/2

и совпадает с формулой, полученной в работе [5]. Видно, что форма стационарной ударной волны то­

ка t(т), распространяющейся в нелинейной линии с распределенными параметрами, представляет моно­ тонно изменяющуюся функцию от / = 0 при t^>оо до

/ = /2 при t - > + оо. Длительность фронта

ударной

вол­

ны конечна и определяется параметрами

феррита

(а,

у), его начальной намагниченностью яд и установив­ шимся значением амплитуды действующего в феррите магнитного поля Н = р012. Безразмерное время нара­ стания начального участка фронта волны при этом равно

Тн=[ауроТо(1—nii)I2l (1 + а 2)]-1.

Отмстим, что для линии передачи с нелинейными сегнетоэлек-

трическими

конденсаторами (<T» = Loi)

система уравнений (1.28)

при £=1

и То= 0 имеет вид

 

 

 

йги

йгр

 

 

-^=[4*Р /М Л 2- 1 ) ] ^ г ;

 

 

^ - = Ро*о (1 —Рг)!(и),

где Д2 =

1 +

4пР (1 — Pi)fb0Ui (t/, = 0,

р г = 0).

Для

сравнительно слабых полей

f (u ) — exр (—ha/u) решение

системы относительно напряжения и(т), выраженного в неявном виде, с учетом граничных условий, находится сравнительно просто:

expI—2/i0/ ( l —P i)t/2]Ei{(l/«— (l—/Ji)/(l+Pi)t/2]} —

 

—Ei {/to{l/u—1 /£/2]} —2РоТо exp (—/loT/L^),

 

где Ei — интегральная показательная функция.

л

45

д нелинейной линии передачи с сосредоточенными па­ раметрами, соответствует учет в (1.28), по крайней мере, двух членов суммы по k * \ Система (1.28), про­ интегрированная дважды по времени с учетом гра­ ничных условий /i = 0, di/dx = 0 и т = т 1 при — оо; приобретает в этом случае вид

 

 

 

dH

 

di

-i — m- m.i

 

 

 

 

 

dz2

 

dz

 

 

 

 

 

 

dm

:5, (1 — rrf) i,

 

 

(1.30)

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

где

обозначено

тс= т г[12/(Д2—1)]1/2;

6i = <xypo£ W (1 +

+ a )2(l—nii)\

Q = [12/(A2—1)]№/Д2

и

введено новое

безразмерное

время гНов= Йтстар и безразмерный ток

бгов= (1

fTli) icr&p/I

 

 

 

 

описываемой

Состояния

равновесия системы,

i,

(1.30), в фазовом пространстве di/dx,

т находятся

из условий d2i/dx2=di/dx = dm/dx=0,

из

которых сле­

дует:

1)

i = 0,

di/dx=0,

m — mi,

2,

3) / = ± 1—ти

di/dx = 0, т —± 1

(знак ±

соответствует волнам с раз­

личным направлением тока). Эти состояния равнове­ сия являются «грубыми», поэтому топологическую структуру их исследуем по линеаризованным уравне­ ниям.

Характер состояний равновесия /(0, 0, пц),

2(1 —

т, 0, 1) и 3(—1—nii, 0, —

1) определяется соответ­

ственно характеристическими

уравнениями

 

А.3+ТсЯ + Х—б1(1—tni2) =0;

 

(Л2+ тД +1)[М -2б1(1+т1)] = 0.

(1.31)

Здесь и в дальнейшем мы будем рассматривать линию, нелинейным элементом которой является только индуктивность с ферритом.

**> Отметим, что постоянная тг предполагается малой' по

сравнению с А. В противном случае в (1.30) надо учитывать производные более высокого порядка малости.

47

Нетрудно видеть, что состоишь равновесия 1 яв­ ляется седлом-фокусом (при ъу<2) или седлом (при тс>2), а состояния 2 и 5 ^ устойчивыми фокусами (при тс<2) или устойчивыми узлами (тс>2).

Структура фазового пространства изображена на рис. 1.14. Видно, что в рассматриваемом случае суще-

Рис. 1.14. Структура фазового пространства (а, б) и зависимость

намагничивания и тока от времени в стационарной ударной вол­

не (в).

ствует одна интегральная кривая — сепаратриса седла (седла-фокуса), описывающая переход из состояния 1 в состояние 2. Соответствующая этой кривой функция

i(t) — зависимость

тока от

времени

в стационарной

ударной

волне при

условии

(г/ро) <

яг| (1—m^M/polz

(тс< 2),

приведена на рис. 1.14,6. Заметим, что намаг­

ничивание изменяется монотонно [это видно и непосред­ ственно из (1.30) при t> 0], а закон приближения т к единице не зависит от параметра тс[^= —2Si(l—т)].

48

Величина тока колеблется около стационарного зна­ чения /2, причем коэффициент затухания колебаний и их частота полностью определяются параметрами ли­ нии н не зависят от релаксационных свойств феррита.

При (г/р0) >лт| (1—>щ)М1ра1г (т0> 2 )

изменение i(x)

за фронтом ударной волны апериодическое.

Как уже отмечено, величина т (т )

в обоих случаях

изменяется монотонно от значения т = т1 при %—>— оо до m = m2 = 1 при т->-+оо (рис. 1.14,в). При этом, как нетрудно видеть из второго уравнения системы (1.30), максимальная скорость изменения т(х) во времени определяется значением параметра q. Большим вели­ чинам q соответствует более быстрое перемагничивание феррита за время, малое по сравнению с пери­ одом колебаний тока. В фазовом пространстве (рис. 1.14,в) это соответствует быстрому движению изображающей точки по сепаратрисе седла-фокуса на плоскость т = 1, где в дальнейшем она по спирали приближается к состоянию устойчивого равновесия

(т = 1, / = /2).

Величины, характеризующие структуру стационар­ ной ударной волны, можно оценить вторым методом, рассмотрев корни характеристических уравнений (1.31). Нарастание тока на начальном участке фронта ударной волны определится корнем первого характе­ ристического уравнения, имеющим положительную действительную часть. Для безразмерного времени на­ растания фронта ударной волны тн~;А,-1 при тс= 0

получим *>

 

тн = 2/(3) ‘/*sh{arsh[3 (3)

(1- m i2)/2]/3}.

*> Для случая тсф 0 выражение для тн более громоздкое и поэтому не приводится. Отметим, однако, что исходная система уравнений (1.30) справедлива лишь при тс <С1; следовательно, учет потерь в линии дает в первом приближении малую поправку к длительности фронта волны, но существен для описания про­ цессов за ее фронтом.

4— 674

49

Соседние файлы в папке книги из ГПНТБ