книги из ГПНТБ / Богатырев Ю.К. Импульсные устройства с нелинейными распределенными параметрами
.pdfсоответствует значение т = 1. Точки плоскости, заклю ченные внутри заштрихованных областей, соответст
вуют значениям |
параметров r/р |
и gp, |
при которых |
в отрезке линии возбуждаются |
колебания. Зона |
||
D (0)— область |
устойчивости. |
|
справедливы |
Проведенные |
выше исследования |
лишь для распределенных активных линий. Однако они качественно правильно описывают устойчивость активных линий задержки с дискретными параметрами (т. е. систем с пространственной дисперсией) для ти пов колебаний, длина волны которых много больше расстояния между ячейками. В случае их соразмер ности необходимо специальное рассмотрение, прове денное ниже.
Устойчивость волновых систем с дискретными па раметрами. Строгое исследование условий самовозбуждения волновых активных систем с периодической про странственной структурой необходимо проводить исхо дя из линеаризованных диф-
ференциально разностных уравнений.
Рассмотрим однородный отрезок активной линии пе редачи, состоящий из Л^оди-
наковых ячеек с произвольными импедансами Z(p) и Y (р) (рис. 3.5). В изображениях по Лапласу при ну левых начальных условиях линеаризованное уравне ние для п-го звена линии имеет вид
un+ i ( p ) — 2un(p) + u „ - i ( p ) = Z ( p ) Y ( p ) u ( p ) . |
(3.13) |
с граничными условиями для н= 1 |
|
Mi= [l + yZ+ Z/ (Zo-f-Zi)] i U z = k o ~ 1ti2 |
|
и для n = N—1 |
|
Ujv—1—[1 +YZ + Z/ (Zjy+ Zjy_i)]_1Ил'_2= |
2- |
140
Дисперсионное уравнение легко находится в ре зультате подстановки в (3.13) значения Ыгг~ехр(у/г). В итоге получим
е7 + е~т =:2 + Z y . |
(3.14) |
Из дисперсионного уравнения (3.14) видно, что в от резке линии возможно распространение двух волн одной и той же моды с постоянными распространения
±Y- Поэтому, отыскивая решение уравнения (3.13) в виде
ип(р) = А(р)е1п + В(р)е т",
с учетом граничных условий нетрудно составить ма трицу определителя А(р) (по аналогии с процедурой, использованной ранее для распределенной модели), а затем, раскрыв определитель А (р)=0, получить об щее характеристическое уравнение в виде
(е1 — К) (ет — kN) — (е-т — ka) (е-т — kN) exp X
Х [2(Л Г-2)у] = 0. |
(3.15) |
Как и ранее, положив в (3.15) Z0iiv = 0, оо, |
получим |
справедливое как для симметричных, так и несимме тричных нагрузок соотношение (для Z1; ZN_2 = = Z/2) 1—exp[2(iV—1)у] = 0, из которого и определяет ся постоянная распространения
у —dzjnm/ (N—1) =jk
(k = 2n/X, X = 2(N—l)/m ), m = 0, 1, 2, ..., N—1.
В качестве примера рассмотрим устойчивость ге нератора с эквивалентной схемой, представленной на рис. 3.2, при R = 0. Границы Д-разбиения найдем из дисперсионного уравнения, положив в нем р = /со и y = jnm/(N—1). После несложных преобразований, по-
141
лучим равенство
(Р^)3+ РЯ-^- + 2 -^ Г ( 1 - со5 j ^ T ) ==0>
структура которого подобна выражению (3.12). Из этого равенства определяем границы устойчивости для минимального значения т —О, pg = 0 и для максималь ного т= (N—1), (Х= 2 ячейкам):
г/р= —pgCo/lpg2+ 4C/Ci]Ct.
Таким образом, качественная картина разбиения плоскости параметров рg и r/р на области возбужде ния для генератора с дискретно-периодической струк турой осталась такой же, как и в случае распределен ного генератора (рис. 3.2). Отличие состоит лишь в том, что число различных областей в данном случае получилось конечным (равным N—1) и, кроме того, изменилась граница области, соответствующая т = = N— 1 (77].
На практике не всегда целесообразно использова ние волновых систем, полностью заполненных актив ной средой, или линией задержки, в каждое звено которой включены активные элементы. Действитель но, ллина отрезка линии, по существу, обеспечивает только необходимый период следования колебаний, а активные элементы — выполнение условий самовоз буждения генератора. Если потери в системе не вели ки, то нет смысла полностью заполнять ее активной спелой (включать активные элементы в каждое звено линии). Достаточно заполнить ее на такой длине, ка кая необходима (с небольшим запасом) для выполне ния условий самовозбуждения системы, а остальную часть линии оставить линейной. Такая система будет уже кусочно-однородной. Приведем (ввиду громоздко сти вычислений) лишь краткое изложение анализа устойчивости такой системы.
142
Устойчивость кусочно-однородных волновых систем. Пусть имеется система (рис. 3.6), состоящая из N отрезков линий пере дачи. В операторной форме уравнение для напряжения в п -м отрезке линии имеет вид
- - |
( Р) У« ( р ) Un ( Р ) - о. п = 1. 2......... |
N . |
На стыке отрезков должны выполняться условия непрерывности
« Л - 1 U n-I. />) = « n ( fn - l. PY>
U n - A l n - i P ) = « n ( P )
— Ря-I M 5 7 « n - i (/n -1. P ) '
а на границах системы краевые условия
«1 (0, р ) |
d |
(0, /?); |
=r a, ( p ) ^ J - U j |
||
|
d |
|
UN ^ N ' P ) = |
§N ( P ) ~dz |
UN V n - P Y |
где ctn — 1/ZnFn, P„ = l/ Z „F „+f.
Общее решение для каждого отрезка имеет вид
u „ ( z , р ) = А п ( р) |
ехр [уп ( p ) z ] + B n ( р ) |
ехр[— уп(/>)2], |
|
a Yп ( р ) — определяется |
дисперсионным |
уравнением \ п 2 ( р) = |
|
= Z n ( p ) Y n ( р ) . |
|
|
|
Для |
отыскания определителя системы |
А ( р ) необходимо со |
|
ставить |
блочную матрицу по рецепту, аналогичному использован- |
Рис. 3.6. Эквивалентная схема кусочно-однородной волновой системы.
ному выше для однородной системы. В общем виде матрица по лучается довольно громоздкой, однако, разбив ее на блоки и вос пользовавшись некоторыми свойствами блочных матриц [82], мож но в результате рекуррентных преобразований свернуть ее к сле-
143
дующему виду: |
|
|
|
|
|
|
|
|
^ ( р ) |
| |
* п - l M n - 1 , п - 1 |,Х |
||
|
|
|
\п=2 |
|
/ |
|
|
X (— 1)N |
|
N — 1 |
|
|
|
|
|
(I xn‘W„,n |- Л4п+1,„) + PiV^.v.v |
||||
где |
а,, |5д,— матрицы краевых |
условии |
О |
|||
|
|
|
-1 а |
P/V— |
||
|
|
|
О О |
W : |
||
хп, |
хп — матрицы граничных условий на стыках отрезков |
|||||
|
|
|
-1 |
0 |
1 |
0 |
|
|
|
-1 |
к ' хн — 0 |
ап |
|
и, |
наконец, М п,п и M n , n - i матрица |
вида |
|
|||
|
м |
= |
ехр(у„/и) |
ехр(— Y„/„) . |
||
|
|
|
Ynexp(Y„/„) |
—уп exp (—Yn/n) |
||
|
M |
^ |
exp(Y„/n-,) |
exp(—Yn/„_,) |
||
|
|
|
Yn exP (Y-n^n-i) |
— Yn exp (— '(n l n - i ) |
Здесь первый индекс матрицы указывает индекс постоянной распространения у (внутри матрицы совпадающий с номером отрезка системы), а второй — индекс, стоящий при длине отрезка.
Положив определитель этой матрицы равным нулю (Д (р)=0) и раскрыв его, получим общее характеристическое уравнение ку сочно-неоднородной системы. Оно определяется довольно слож ным выражением, однако на практике обычно используют систему из двух и гораздо реже из трех отрезков. Для этих случаев харак теристическое уравнение получается не слишком сложным и впол не исследуемым аналитически. Общее же выражение Д ( р) можно рассматривать как алгоритм для вычисления областей устойчиво сти произвольной кусочно-неоднородной системы на ЭВМ.
3.2. СТАЦИОНАРНЫЕ ИМПУЛЬСНЫЕ КОЛЕБАНИЯ
Спустя некоторое время после выполнения условий самовозбуждения в автоколебательной системе с рас пределенными параметрами устанавливаются стацио-
144
парные колебания в виде периодических волн, являюющихся решениями обыкновенных нелинейных диф ференциальных уравнений. В фазовом пространстве системы этих уравнений периодическим бегущим вол нам соответствуют предельные циклы. Число и пара метры предельных циклов, а следовательно, форма и скорость распространения стационарных волн опреде ляются видом нелинейности и параметрами линии.
В релаксационном распределенном генераторе установившиеся колебания представляют собой ре зультат взаимодействия многих типов волн. Возмож ность получения различных типов волн зависит от ве личины нелинейности (чем сильнее нелинейность, тем больше типов волн), а эффективность взаимодействия определяется близостью их скоростей распростране ния, т. е. зависит от дисперсионных свойств генерато ра (чем слабее дисперсия — пространственная или частотная-—тем выше эффективность взаимодей ствия).
Таким образом, уже в линейном приближении исходя из вида дисперсионной характеристики можно судить о возможности получения релаксационных ко лебаний в генераторе с распределенными парамет рами.
Так, например, для генератора с эквивалентной схемой, приведенной выше на рис. 3.2, скорость малых гармонических возмущений, найденная из дисперсион ных соотношений (3.11), равна
»«2=[1 +rg+ (Rg/k2)]/[LC+ (1 +rg)LCl + rRCCi].
Из этого выражения следует, что присутствие в гене
раторе низкочастотных потерь |
(/?=т^0) приводит к за |
|||
висимости |
скорости vK от волнового числа |
&= 2я/А |
||
(к = 1/т) |
собственных колебаний генератора. |
При |
||
этом скорости гармонических |
составляющих |
малых |
||
длинноволновых возмущений |
(соответствующие |
ма |
лым значениям т) отличаются существенно и практи-
10— 674 |
145 |
чески не участвуют во взаимодействии (дисперсия велика), следовательно, длинноволновые стационар ные волны не являются релаксационными. Скорости же коротковолновых возмущений (большие значения т) близки (дисперсия мала) и, эффективно взаимо действуя, образуют короткие релаксационные волны.
С другой стороны влияние дисперсии, т. е. зависи
мость vK= f(k) для |
волн, соответствующих |
одним и |
тем же значениям |
т, будет тем слабее, чем |
меньше |
величина Rg/k2 по сравнению с (1+ rg), т. е. чем боль ше длина линии (поскольку k = 2nm/l). Значит в рас пределенном генераторе с дисперсией возможны и длинные релаксационные волны, если протяженность генератора достаточно велика. В генераторе без низ кочастотных потерь (R = 0)vK=£-f(k), поэтому вид всех стационарных волн наиболее релаксационный незави симо от длины генератора.
Перейдем теперь к рассмотрению нелинейной за дачи об определении формы и основных параметров стационарных релаксационных волн генератора с экви валентной схемой рис. 3.2 (при /? = 0).
Строго говоря, вначале необходимо было бы убе диться, что стационарному решению (волне) на фазо вой плоскости соответствует предельный цикл, и дока зать его устойчивость **). Это особые задачи теории колебаний и останавливаться здесь на них нет особо го смысла. Тем не менее отметим, что для исследуе мой системы эта работа проделана в [77], где и пока зано, что для значений скоростей стационарных волн v2 - . \/L(Ci + C) соответствующий предельный цикл устойчив.
*) Вид дисперсии (временная или пространственная) при этом не имеет значения.
**> При этом, отнако, можно показать, что устойчивость пре дельного цикла является необходимым, но не достаточным усло вием устойчивости стационарного решения.
146
Уравнения, описывающие стационарные процессы в генераторе, имеют вид
duld\ — y,
* 4 г = |
{1- LC^ + - § ■ |
г («) ) у + ! (“)/№. |
|||||
|
|
|
|
|
|
|
(3.16) |
с краевыми условиями |
|
|
|
|
|
|
|
|
w(l) = w (| + /); |
|
|
|
|||
|
du(Q fd% = du(l + t) /dl, |
|
|
|
|||
где /( и ) — нелинейная |
зависимость тока |
активного |
|||||
элемента от напряжения, приложенного к нему; |
ц = |
||||||
= r(C/L) (1—LCiV2)/n0— малый |
безразмерный |
пара |
|||||
метр *). |
системы |
(3.16) |
должны |
удовлетворять |
|||
Решения |
|||||||
периодичности (с периодом |
Х = 1/т), |
что, |
очевидно, |
Рис. 3.7. Предельный цикл системы (3.16) на фазовой плоскости.
возможно лишь при дискретных значениях скоростей
v = v { m ) .
Найдем параметры стационарных бегущих волн
вслучае ji<C 1.
1.Если [х— Ю (г— й)) предельный цикл системы
(3.16) становится разрывным, состоящим из участков быстрых и медленных движений (рис. 3.7). Быстрые движения осуществляются по прямым у = const, а тра
*> Здесь нормирующий параметр имеет размерность (м/Ом) и определяется из выражения для конкретной зависимости 1(и).
10* 147
екторией медленных движений является кривая
у = —Ао1(и), |
(3.17) |
где
Ло« vL/[ (1—v2LC0)].
Длина волны стационарных импульсных колеба ний, бегущих в кольцевом генераторе со скоростью v, при достаточно малом ц близка к периоду движения изображающей точки по разрывному предельному циклу и определяется непосредственным интегрирова нием следующего выражения:
|
0__ |
1 |
С da |
I |
|
1 Г |
da |
|
|
|
|
“ATJ |
|
" A J |
"П ир |
|
|
||
|
|
|
Wi |
|
|
«3 |
|
|
|
В частности, для /( « ) = —g ( l—аи2)и длина волны |
|||||||||
|
|
|
2 |
“ 3,4 |
|
du |
|
|
|
|
Я = |
|
С |
|
|
|
|
||
|
|
gA0 |
J |
(1 — aus)u |
|
|
|||
|
|
|
|
|
|||||
|
|
|
|
« 1,2 |
|
|
|
|
|
Пределы интегрирования Mi,2 (или МзД легко нахо |
|||||||||
дятся с |
помощью |
рис. |
3.7. |
При Мгд, |
dy/du — |
||||
= —Ang(1—Зам2), |
следовательно |
М2,4= Т ( З а ) _1/2. Точ |
|||||||
ке и= и\ соответствуют значения |
у = А 01 (и2 л) ; |
отсюда |
|||||||
получается |
кубичное уравнение |
м3—м/га+ 2(За)_3/2=0, |
|||||||
единственный действительный |
корень которого |
равен |
|||||||
Mi,3 = + 2 (3 a)_1/2. |
|
|
|
|
|
|
|
про |
|
Подставляя пределы интегрирования, после |
|||||||||
стых преобразований находим длину волны |
|
|
|||||||
|
*,«2,1(1—v*LC0)/gvL. |
|
(3.18) |
Временной период колебаний равен Т = Ям-1. Скорость распространения различных типов волн, соответствую щих дискретным значениям длин волн K=l/m, опреде
148
ляется приближенным соотношением |
|
Vm ~ ~ ~ ^ с 7 + { ( l £ r ) + 1/LC»} ' • |
(ЗЛ9) |
Естественно, что при g = О, vm= (LC0)“1/2. |
Размах |
импульсных колебаний AU = u3—«i = 4/(3a) 1/2 и ампли туда импульсов U= u3—Пг= 3/(3а)1/2, как и следовало ожидать, оказались зависящими только от коэффици ента нелинейности а.
Проведенное рассмотрение справедливо, когда па раметр р при старшей производной системы (3.16) достаточно мал. При принятой аппроксимации I(и) параметр
И = г(1 —v2LCi)/gp* (ро= g ) ,
откуда следует, что значение p-Cl и, следовательно, импульсный режим работы распределенного генерато ра достигается в случае, когда произведение сопро тивления растекания г на сопротивление активного элемента, например туннельного диода, в рабочей точ ке вольт-амперной характеристики много меньше квадрата квазиволнового сопротивления р2=L/C.
При р<С1, но не равном нулю, время движения изображающей точки на участках и2—н3 и щ—щ пре дельного цикла равно, соответственно, длительности фронта и спада импульсных колебаний. Для опреде ления этих длительностей, строго говоря, необходимо получить точное решение системы (3.16). Поскольку сделать это аналитически затруднительно, можно вос пользоваться приближенными оценками их (как это было сделано в предыдущих главах) по корням ха рактеристического уравнения линеаризованной систе мы (3.16) или на основании решения укороченных уравнений.
Отметим, что в кольцевом генераторе структура фронта волны и его характерные временные и прост ранственные параметры будут очевидно такими же,
149