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

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

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

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

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

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

30

ударные волны на сравнительно малой длине (при не­ большом числе звеньев) линии.

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

Для стационарных волн тока и напряжения, рас­ пространяющихся в линии в направлении возрастаю­

щих

п с постоянной скоростью

V = Vp, выполняются

соотношения: in±i(t) = in (/+ A /);

un+i = un(t4zM ), .

где

А ^=1/нр — постоянная, равная времени задержки

волны па одно звено формирующей линии. Эти равен­ ства позволяют от системы дифференциально-раз­ ностных уравнений ( 1.2) перейти к системе нелиней­ ных дифференциальных уравнений с запаздывающим аргументом.

Фазовое пространство системы дифференциальных уравнений с запаздывающим аргументом является, как известно, функциональным (бесконечномерным), что препятствует полному, хотя бы качественному исследованию их решений. Поэтому для исследования стационарных решений (стационарных ударных волн) применяются два приближенных метода [28, 22].

Первый метод предполагает переход к укорочен­ ным уравнениям системы (1.2) для стационарных ре­ шений в предположении малости параметра At- (по сравнению с некоторым характерным временным па­ раметром). Метод заключается в том, что в системе нелинейных дифференциальных уравнений с запазды­ вающим аргументом переменные в моменты t±At разлагают в ряд по степеням At. Затем, ограничи­ ваясь конечным числом членов разложения одинако­ вого порядка малости, приходят к системе обыкновен­

ных дифференциальных

уравнений

относительно

i(t)

ц u(t) [28]. Последняя

решается

аналитически

или

31

исследуется качественными методами теории колеба­ ний. Метод пригоден и для исследования нестацио­ нарных решений*) н, по-существу, является естест­ венным обобщением метода Р. В. Хохлова для вол­ новых систем, описываемых дифференциально-разно­ стными уравнениями.

Второй метод заключается в исследовании вида корней дисперсионного уравнения системы (1.1) — (1.2), линеаризованной в окрестности особых точек 1 и 2, характеризующих состояние равнозесия перед

фронтом

ударной волны (при i— »— 8) и далеко за

ним (при

t— *+оо). Известно [8], что

структура

ста­

ционарных ударных волн на фазовой

плоскости

(или

в фазовом пространстве) описывается фазовой траек­ торией, соединяющей особые точки 1 и 2. Эта траек­ тория определяется лишь теми корнями Xi^h) диспер­ сионного уравнения, у которых Im М(,{)>0; ImX2(fe)<0. При этом длительность начального участка фронта ударной волны определяется наименьшим по модулю корнем М- Второй метод, в отличие от первого, позво­ ляет получить количественные оценки некоторых па­ раметров структуры не только слабых, но и сильных ударных волн, когда М уже нельзя считать малым.

Рассмотрим теперь структуру стационарных удар­ ных волн для конкретных формирующих линий.

Линии с нелинейными элементами, описываемыми квазистатическими уравнениями связи. Поведение волн тока п напряжения, например, в линии, эквива­ лентная схема которой представлена на рис. 1.8, описывается системой нелинейных дифференциальных уравнений с запаздывающим аргументом для функ­ ций in (t) = i( t ) и un(t) = u (t) **).

*> В этом случае разложение ведется по п, т. е. осуществляет­

ся переход от дифференциально-разностных уравнений к укоро­ ченным уравнениям в частных производных [ 102].

**> Влияние сопротивления R и г качественно аналогично про­ явлению вязкости феррита. Оно обеспечивает диссипацию энергии в области фронта ударной волны.

32

и (t -j- Д^) ■—■и (t)

йф

dt

 

 

— г

{Q [и + ДО] — Q [«(011;

Х - ^ { ф И О ] - ф И*-ДО]}-

(1-23)

Нетрудно видеть, что любые постоянные u=U = const и i = I = const являются решениями этой системы. Поэтому естественно искать в рассматриваемом слу­

чае

стационарные решения (стационарные волны)

Рис. 1.8. Эквивалентная схема

звена

нелинейной искусствен­

 

ной линии.

в виде перепада между двумя состояниями равнове­

сия: / ь

U1 при t— »— оо и

h, U2 при t— v+oo (I2,

U2> h,

Ui) с соблюдением условия dhuldtk= dhi/dth= 0

при t— *± оо для любых k = 1, 2, 3, ...

Интегрируя полученную

систему по достаточно

большому интервалу времени {/ь t2}, охватывающему весь процесс перехода из первого состояния равно­ весия во второе, можно определить граничные усло­ вия и, следовательно, время задержки At=\Atp, равное

А/Р(£/а—С/*) = Ф (/г) —Ф (Л ),

(1-24)

А*р (/2- / 1) = Q( U2) —Q(Ul).

(1,24а)

Откуда

Д/2= [Ф(/2) —<D(/i)]X

iXtQ(И2 ) - Q i U f i m h - h ) {Uz-Ui)}-

3-G74

33

Выражение (1.24) можно рассматривать как гра­ ничное условие для ударной волны, распространяю­ щейся в формирующей линии с нелинейными дискрет­ ными параметрами *>. Оно полностью совпадает с гра­ ничными условиями для волн в линии с распределен­ ными параметрами. Это значит, что пространственная

Рис. 1.9. Графическое определение состояний равновесия:

а) в линиях с нелинейной емкостью или индуктивностью; б) в линии с не­

линейной емкостью и индуктивностью.

дисперсия (так же, как и частотная [8]) не изменяет вид граничных условий для фронта ударной волны.

Выражение (1.246) допускает простую геометри­ ческую интерпретацию: скорость распространения фронта ударной волны обратно пропорциональна квадратному корню из произведения тангенсов углов наклона отрезков, секущих кривые Ф(/) и Q(u),

*> Граничные условия могут быть получены и непосредствен­ ным интегрированием системы дифференциально-разностных урав­ нений ( 1.2), т. е. для нестационарных ударных волн, в предполо­

жении, что напряжение и ток изменяются медленно вне области фронта ударной волны [102]. Задача определения величины At, входящей в коэффициенты уравнений (1.23), может быть также сведена к нахождению собственных значении параметра At, для которых система (1.23) имеет стационарные решения конечной амплитуды [28].

34

проведенные через точки, соответствующие положе­ ниям равновесия До и UU2 (рис. 1.9,а).

Достаточно очевидно, что во втором положении равновесия /2 или U2 не могут быть заданы незави­ симо в силу вытекающего из (1.24) условия

(/2- / 1) [Ф (/2) - Ф (/1)] = ( U2- £/i)[Q(U2) - Q ( U 0].

Из граничных условий нетрудно найти уравнение для «волнового сопротивления» липни г,, для удар­ ной волны, которое можно определить как отношение перепада напряжения Up к перепаду тока / р на фрон­ те ударной волны:

 

7 __

Щ _ Ut - 1/. _

Ф (/,)-Ф (Л )

'

 

р

/ , - Л

Q (и.г) — Q (I/,)

В

частности,

если положить U1= I l= 0,

а магнит­

ный

поток

и электрический

заряд определить как

Ф= L{h)Iz, Q = C(U2)U2,

Ф(0) =Q(0) =0, то Zp= [I(/2)/C (t/2) p .

Исследуем структуру стационарных волн сначала первым методом. Предположим для простоты, что с])—линейная функция тока <l> = L0i и сведем систе­ му (1.23) к одному уравнению относительно и. Затем, разложив в степенной ряд по малому параметру At переменные с аргументами t±At, получим укорочен­ ные нелинейные уравнения вида

2Д2ь

d2* f

, г

d Q \

(2k)\

dz2*

x0

dz J

(1.25)

где безразмерное время x = t/x« (xr = Lo/Rx0, Л=А//т0). Если в правой части (1.25) ограничиться только членом, соответствующим к= 1, то получим выраже­ ние, совпадающее с уравнением липни с распределен­

3*

35

ными параметрами. Поэтому простейшей модели, учитывающей пространственную дисперсию, соответ­ ствует учет в уравнении (1.25) по крайней мере двух членов суммы по k. Проинтегрировав дважды полу­ ченное выражение (с учетом условий для т— > ± о о ) , приходим к нелинейному дифференциальному урав­ нению четвертого порядка.

Вреальных формирующих линиях сопротивление

гмало (r = pr) и, кроме того, всегда выполняется неравенство r<^R. Поэтому проведем второй этап укорочения полученного уравнения, отбросив произ­ водные 3-го и 4-го порядка, коэффициенты при ко­

торых— величины

2-го

порядка

малости ( ~ ц 2).

В итоге получим укороченное уравнение*)

d2u I

12 /

. rQ' (м)л ' da

dx2

Д2

то )

Лх

§-[Q («) -Q (£/,)]-% -

-

п,)д2=о.

(1.26)

 

 

 

 

хо

 

 

 

 

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

(1.26) па

фазо­

вой плоскости определяются из условий:

 

 

■ ж = у = о;

 

 

 

 

 

= о,

из которых также можно найти

величину A (At).

Линеаризуя (1.26) в окрестности особых

точек,

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

 

 

Я2 +

-4 т

+

хг) я -

 

(xj — Д2) =

О,

(1.27)

*> Отметим, что уравнение линии с распределенными пара­

метрами получается

из

(1.26),

если

отбросить

члены

~ А 4.

Его решение:

1 Г [Q(u)-Q(t/,)I£oA§-(«-t/.)A*

 

а (х) =

 

Д2 J

 

xR +

rQ'

(и)/х0

 

d’*1-

 

36

корни которого определяют характер состояний рав­ новесия U0. Здесь xr=rQ'(U0)/xo, xq2=L0Q'(U)/x«2.

Вид и знаки корней (1.27) в основном зависят от со­ отношения между производными Q'(U0) и постоян­ ной Д. Более полное представление о характере со­ стояний равновесия, определяемых этими корнями

Рис. 1.10. Плоскость параметров Л0, В00= 6(тя -|-Тг)/Д2;

.во= 12(t2q—Д2)/Д4].

при всех возможных значениях параметров тг, Xr, хq и Д дает диаграмма, иллюстрирующая разбиение плоскости параметров (рис. 1.10).

Не конкретизируя вида функции Q(u), рассмот­ рим качественно структуру стационарной волны для двух наиболее интересных случаев.

1. Характеристика Q{u) имеет вид, изображенный на рис. 1.9,а. Состояниями равновесия будут точки

пересечения Uo=Uyt2 кривой

Q(u)Q(Ui) с прямой

Q (u)~Q (U i)= A 2x02{u— Ui)IL0

(tg а = (Дт0)2/10). Со­

стояние равновесия в особой точке 1, т. е. перед

фронтом ударной волны

(t— ъ—о о ), всегда неустой­

чиво (седло), так как

при U0=U i, корни

(1.27) дей­

ствительные,

разных

знаков.

Состояние

равновесия

в особой точке 2, т.

е.

за фронтом ударной волны

(t— >-+

о о ),

устойчиво

(устойчивый узел, если Во2>

>4Л 0),

или

устойчивый

фокус

(при обратном нера­

венстве). Корни (1.27) соответственно действитель­ ные, отрицательные и комплексные с отрицательными

37

реальными частями. Особой является ситуация, когда в одной из точек равновесия (например, Ui) секущая совпадает с касательной. При этом отклонение от положения равновесия не определяется линеаризован­ ным уравнением (1.26) и надо учитывать квадратич­ ные члены в разложении Q(u).

Структура фазовых плоскостей для рассмотренных выше случаев изображена па рис. 1.11. Видно, что

Рис. 1.11. Фазовые портреты и форма стационарных ударных волн.

существует только одна интегральная кривая — се­ паратриса седла, — описывающая переход из состоя­ ния 1 в 2. Этой фазовой траектории соответствует стационарная ударная волна, экспоненциально нара­ стающая на начальном участке фронта и апериоди­ чески (рис. 1.11,а) или колебательно (рис. 1.11,6) приближающаяся к постоянному значению за фрон­ том волны.

2. Характеристика Q{u) имеет три точки пересе­ чения (для ВОЛНЫ ОДНОЙ полярности) Uq— U1,2,4* по­ явилось новое состояние равновесия t/4. В этом слу-

38

чае состояние равновесия в особых точках U1)2 устой­ чиво (устойчивый узел или фокус), а в особой точке Щ неустойчиво (седло). Следовательно, в линии мо­ гут существовать стационарные волны в виде пере­ пада напряжения с амплитудой U>U!l (от u = f/4 при t— >— оо до u = U 2 при t— >-+ оо) или U<Ui (от w= J74

при t— v+oo до u —U1 при t— *—оо).

Если нелинейны оба реактивных элемента, а вид зависимостей Q(u) и Ф(/) произволен, число и харак­ тер состояний равновесия можно определить, лишь

построив

обобщенную

характеристику

Q(u)=f(i)

(или

Ф (1)= ф («)). Из

граничных условий имеем

 

(1—h )

[ф (I) - Ф (/i)] = (и—C/i)[Q (и) - Q

(t/i)].

Задавшись значениями величин i = l i ,2

и u —Ui, а сле­

довательно,

Ф (/1,2) и

Q(U1) и определив из этого

равенства или графиков функций Ф(/)

и Q(u) значе­

ния и = ии

 

и.2 , и3, ...

и

Q(ui),

Q { u2),

Q ( u3), ... при

изменении

тока =

i2,

t3, ...

строят

кривую Q(w) =

=Д(г)

(рис.

1.9,6). Через точки Д, Q(h)

и /2, Q(/2)

проводят прямую, тангенс угла наклона которой ра­

вен Ai/2=[Q(£/2)—Q(C/i).]/(/2—/ 1),

Число точек

пере­

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

равновесия

(числу

особых точек на фазовой плоскости).

Стационарные волны особого вида. Интересно от­ метить, что при отсутствии потерь в обоих рассмот­ ренных выше случаях уравнение (1.26) допускает другие стационарные решения. Действительно, при тн= т,.=0 особые точки с состояниями равновесия типа устойчивый узел или фокус вырождаются в центры. Фазовые траектории теперь соответствуют не удар­ ным (в принятом нами определении), а стационар­

ным волнам особого

типа: сепаратриса — уединенной

волне

(солитону) *1,

замкнутые

траектории внутри

сепаратрисы — периодическим

квазигармоническим

волнам [28, 102].

 

 

*)

Уединенные волны будут рассмотрены в § 5.1.

39

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