книги из ГПНТБ / Богатырев Ю.К. Импульсные устройства с нелинейными распределенными параметрами
.pdfКачественное исследование структуры ударных волн, проведенное выше, не позволяет непосредствен но получить необходимые для расчета формирующих систем количественные оценки характерных парамет ров волны: длительности фронта, амплитуды, периода и времени затухания осцилляций за ее фронтом и т. д. Такие оценки удается извлечь из аналитического или численного решения исходных или укороченных урав нений.
В качестве примера аналитического исследования структуры ударной волны найдем приближенное решение нелинейного диф ференциального уравнения (1.26), допуская аппроксимацию ха рактеристики Q(u) полиномом третьей степени: Q(u )=CoU—D2u3 (Со, Ь 2— постоянные, положительные коэффициенты).
Укороченное уравнение (1.26) в этом случае запишем >в виде системы *>
i£dLJ2
du/dz = у; d y / d z д |
(l/2— и2) и— |
где согласно (1.24) А2= ( 1—D2U2/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— |
—Ы22/Ы12) V2, a ullzz= U 2± U ( U 2—A) 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.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 |