
книги из ГПНТБ / Чугаевский, Ю. В. Элементы теории нелинейных и быстропеременных волновых процессов
.pdfВид функции (40) при разных значениях |
ис , h |
показан на |
||||
р и с.18, а фазовая картина решений уравнения |
|
(37) - |
на рис. 19. |
|||
Особенностью кубической, уединенной волны является |
то , что |
она |
||||
может существовать как |
в виде |
горба, так и в |
виде |
впадины |
(со |
|
ли- и "антисолифоны"), |
Кроме |
того, в отличие |
от |
квадратичной |
||
нелинейности, здеоь стационарные волны любой |
заданной скорости |
|||||
распространения возможны и при амплитудах, |
больших |
солифонных |
||||
( - ис > и > ис ) . |
|
|
|
|
|
|
Параметр L , характеризующий сравнимость, или "степень конфликтности", дисперсии и нелинейности оценивается как отно шение нелинейного и дисперсионного членов в уравнении (351 и
равен для кубической цепочки
|
|
|
|
(43) |
Для кубичеокого |
солифона значение |
L , |
как нетрудно |
убе |
диться, пользуясь формулами (41, 4 2 ), |
равно |
Lc = 6 ,т „ е , |
в два |
|
раза меньше, чем для |
квадратичного (см. ниже). |
|
||
б) Квадратичная |
нелинейность |
|
|
|
Для решетки с квадратичной нелинейностью волновсе уравне ние принимает вид
|
|
.д 2 |
|
|
и |
ди |
___L. |
и„ |
=. о |
|
(4 4 ) |
||
d t 2Bx* |
|
* |
дх |
I |
дх |
|
,9 |
„. ? |
и > |
|
|
||
|
Л ‘ |
|
|
/ / |
dt |
|
|
|
|||||
что |
в случае |
стационарных процессов |
дает |
|
|
|
|||||||
|
|
|
|
|
|
|
|
+ |
|
|
|
|
( « I |
|
При' с |
г. ta а>0 |
последнее уравнение имеет отрицатель |
||||||||||
ный коэффициент при |
втором |
члене, |
что приводит к |
решениям в |
|||||||||
форме медленных кноидальных |
волн |
|
|
|
|
|
|
||||||
|
и - и л |
|
|
|
|
|
|
|
|
|
|
|
(46) |
где |
модуль |
£ |
и частота |
ш |
задаются |
соотношениями |
|
|
|||||
|
Г |
|
OJ 9=. ^Ф(и0- |
Ui) , |
Ug > и |
ut >UZ |
г |
(47) |
|||||
S |
\j *0 - |
|
|||||||||||
|
|
|
2V3C |
|
|
|
|
|
|
||||
а точки uf , |
u2 |
находятся |
по |
формуле |
|
|
|
|
|
и>.г |
|
|
1 U° ± i - |
\l(Uc+«o)Z~ 4 -uo (“c + Uo) ’ |
Uo = ~ ^ |
/ - - ' (48) |
||||||||||
|
При |
стремлении амплитуды |
u0 к |
кноидальный |
процесс |
|||||||||||
вырождается в медленный солифон |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
и = и„ |
|
|
|
|
|
|
|
|
(49) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
ширина которого |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
L |
|
_ |
^ |
С 10 |
|
|
|
|
|
|
(50) |
|
|
|
|
|
|
|
I/с% -сг |
|
|
|
|
|
|
|
||
неограниченно |
растет, |
а амплитуда |
падает |
до нуля при с — с0 . |
||||||||||||
|
В случае |
с г |
|
> |
с* |
уравнение описывает |
быстрые |
кноидаль |
||||||||
ны е волны, |
которые |
записываются |
теми же формулами, |
что |
и |
мед |
||||||||||
ленные, |
за |
исключением |
соотношения |
(4 8 ). |
Последнее |
следует |
за |
|||||||||
мен! |
„ . |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
// |
|
Uc -Ua |
|
|
|
|
|
|
|
|
3 (с г-Со ) |
(51) |
||||
— —£---£. ± |
Y / (“о - " с ) г - £ " о ( “ о - “ с ) >ис |
|
|
|||||||||||||
|
|
|
г |
|
|
|
|
|||||||||
|
В пределе |
и0 |
|
|
кс- |
формула |
(46) переходит в |
уединенное |
||||||||
решение - |
быстрый |
|
солифон |
|
|
|
|
|
|
|
|
|||||
|
и ш и0 |
$eeh2( A = c L ) , |
|
h « |
|
• |
|
|
(52) |
|||||||
|
|
|
|
|
|
|
I |
h |
j |
|
|
|
|
|
|
|
|
Оба |
типа |
волн - |
быстрые и медленные - вполне |
аналогичны |
|||||||||||
соответствующим процессам |
в жидкости, рассмотренным |
в главе I . |
В пределе очень быстрых возмущений ( сод — »- О ) уравнение
(45)принимает вид
^ |
м * ~э~2 |
— О |
|
(53) |
и быстрый солифон вырождается |
в короткий |
|
||
" ■ и° Sl c * ‘ |
( J!T L ) ’ |
- - J 7 ? |
’ |
<и > |
представляющий, |
как и короткий |
солифон |
в жидкости, |
элемеятар— |
61
ное возмущение, в котором в движении находится практически од
на частица (рис.20).
|
|
|
Параметр |
|
|
L |
для квацрат*ич- |
|||||
|
|
ной |
решетки равен |
|
|
|
|
|||||
|
|
|
|
|
L = Фи |
-Д |
|
|
(55$ |
|||
|
|
что |
для |
|
|
т О С* |
|
|
|
|||
-J i |
h |
солифона, |
как и следовало |
|||||||||
ожидать, дает |
1 С = 1 2 . |
|
|
|
||||||||
и(>г) |
|
|
|
|||||||||
|
Неустойчивость вибратора с ква- |
|||||||||||
Р и с . 20 |
|
|||||||||||
дратичным |
|
ангармонизмом |
(7) |
в |
об- |
|||||||
ластях больших амплитуд |
Зсод |
•с и |
и |
UJi |
|
\ |
|
|
|
|
||
|
|
I |
|
вляется™ ™ , |
|
|
||||||
|
|
2 у ~~ " ~~ |
|
— ) |
по-ви |
|||||||
|
|
|
|
) |
|
я5 6 |
* |
|||||
димому, одной из основных причин |
повышенной7 разрушаемости |
ма |
||||||||||
териалов с |
неоднородной |
структурой. |
|
|
|
|
|
|
|
|
|
в ) Простые волны, уравнение Кортевега - де Вриза
Поскольку уравнение (35) допускает стационарные волны с определенной скоростью распространения с } для возмущения,рас пространягощегооя в безграничной среде, например, вправо, зтс уравнение можно проинтегрировать один раз и представить в виде
ди |
|
|
д 3и |
-f-pu |
р |
Зс/ |
- |
|
(56) |
~dt |
+ т |
|
|||||||
<?,■> |
|
|
дх |
~ Сг |
|
|
|||
tn = с* -с* |
, |
|
|
Ct o J€ |
|
(57) |
|||
р -■ с г ~ с* |
|
|
|||||||
Умножая (56) |
на 2 р и , получим для функции |
|
|
||||||
|
|
|
«У= р и.-2 |
|
|
|
(58) |
||
в бездисиерсионном |
приближении уравнение |
|
|
||||||
|
|
|
|
|
|
|
|
|
(59) |
дающее простую |
волну s |
|
|
|
|
|
|
||
В случае квадратичного энгармонизма уравнение |
(44) |
для бе |
|||||||
гущих вправо волн |
можно |
записать |
в форме Кортевега |
- де |
Вриза |
||||
dD |
^ |
dD |
__ |
д 31 > |
_ |
|
(60) |
||
—. |
У" |
D |
£ |
+ т |
~з = |
О , |
|
||
d t |
|
|
дх |
|
дх< |
|
|
|
если положить
62
|
|
|
|
|
D = p и |
p i = . c i i P . |
(61) |
||
|
|
|
|
|
|
- 2 _ , |
|
||
|
|
|
|
|
|
|
|
|
|
|
Достаточно плавные процессы, |
следовательно, и в этом |
слу |
||||||
чае представляются |
в форме простых волн. |
|
|||||||
|
Поскольку |
уравне |
|
|
|
||||
ние |
ut |
+ригик = О |
от |
|
|
|
|||
личается |
от |
|
уравнения |
|
|
|
|||
й/, + иик = О |
инвариант |
|
|
|
|||||
ностью |
относительно |
за |
|
|
|
||||
мены |
и |
|
■ - |
и |
укруче- |
|
|
|
|
ние волн будет происхо |
|
|
|
||||||
дить |
по-разному: |
пер |
|
|
|
||||
вому |
уравнению соответ |
|
|
|
|||||
ствует |
р |
-пила, |
вто |
|
Р и с . 21 |
|
|||
рому |
- |
с^-пила |
(р и с.2 1 ). |
|
|
||||
|
|
§ |
4 . |
Распад на солифоны и квантование колебаний |
|
||||
|
|
|
|
цепочки [6 7 ] |
|
|
|
а) Спектры скоростей и энергий солифонов
Из факта дискретности цепочки вытекает, что длина распро страняющихся по ней волной в частности ширина солифонов, при нимает не сплошные, а дискретные значения Нп = 2 {с л , /? = = 1 , 2 , 3 , . . . Ид (41) и (42) вытекает, следовательно, что спектр скоростей распространения солифонов имеет вид
2 л
|
|
|
|
* * = |
со ^ 4 л 2-/ |
|
|
(62) |
|
|
|
|
|
|
|
|
|||
Максимальной |
скоростью обличает |
самый короткий |
солифон стд~ cf - |
||||||
= - — r.G |
ростом |
ширины солифона |
скорость |
его |
падает |
и в пре |
|||
деле |
п — - |
: о |
приближается к |
скорости |
линейных колебаний це- |
||||
п очки |
(с |
са )- |
|
|
|
|
|
|
|
|
Спектр |
амплитуд солифонов задается выражением |
|
||||||
|
|
|
и „Си')—■ |
d |
d = |
. П Г |
(63) |
||
|
|
|
|
\fd/7 |
|
|
|
|
|
что при больших |
л дает |
убывание по гиперболическому |
закону |
||||||
|
|
|
|
с/?у |
^ |
J L . |
|
|
(64) |
|
|
|
|
|
/7» j |
2 п |
|
|
|
63
Аналогично ннгладит спектр частот кноидальных волн
|
|
со, |
JEs. |
|
|
(6b) |
|
7 4 п г- 7 f t » ? |
2 ft |
|
|
|
|
|
|
|
|
|
||
Рассмотрим |
эволюцию |
распространяющейся |
по |
пепочке |
длинной |
|
( h » t 0 ) |
волны, считая длину |
цепочки |
для |
данного |
процесса |
практически безграничной. Для кристаллов "комнатных" размеров к
таким |
волнам относятся тепловые колебания решетки. Пока возму |
|||
щение |
остается медленнопеременным,дисперсионным членом д^и(<}1 гд ) |
|||
в уравнении |
(35) |
можно пренебречь. Нелинейная же тенденция, |
на |
|
капливаясь, |
рано |
или поздно приведет к укручению волны ,т.е.сде |
||
лает |
ее быстропеременной. Вступление в игру на этом этапе |
дис |
персии вызовет распад волны на солифоны и образование полисоди-
фояного состояния. Тепловые колебания цепочки будут, |
таким |
об |
|
разом, суперпозицией не гармонических, а солифонных |
или,точнее, |
||
кноидальных мод. |
|
|
|
Спектр возможных значений энергий осциллирующих узлов мож |
|||
но построить, приняв, что полная энергия каждого |
узла, попавше |
||
г о в поле действия солифона, равна потенциальной |
энергии |
узла, |
|
находящегося в вершине солифона: |
|
|
|
э(л ) 2-
|
Еп = |
|
Р и° |
у- |
|
___ , |
|
|
(66) |
|
|
|
|
|
2 |
|
2^ |
|
|
|
|
где f t - |
масса |
узла. |
При больших |
п |
вклад второго |
члена |
этого |
|||
выражения пренебрежимо |
мал, |
так что |
|
|
|
|
|
|||
у |
(о0 ft ttc |
* |
Ьи>о |
a; |
3tc?o{t |
|
|
( 6? ) |
||
|
2 |
|
~ |
/г* ’ |
° ~ 4-х ' |
|
|
|||
|
|
|
|
|
||||||
На "стелющемся" участке кривая |
|
1 |
быть |
с большой |
||||||
|
может |
|||||||||
точностью линеаризована. Действительно, выбирая |
некоторое |
до |
||||||||
статочно |
большое |
п - п д |
и 'вводя |
переменную |
А/- ^0 - ^ , |
|||||
энергетический спектр |
(67) можно представить как |
|
|
|||||||
|
Е |
|
|
ta >0 |
|
|
2А/ ) |
|
( 68) |
|
|
|
|
|
Ь U)q |
* 2 / |
|
||||
'ft ?>/ |
|
|
("о~Н )г |
|
|
|
|
|||
Если положить |
|
|
|
|
|
|
|
|
|
|
у = 26 _= |
|
|
|
|
|
Гы0 ??0 |
|
(69) |
||
* |
Л) |
|
|
|
|
|
2 |
’ |
|
|
|
|
|
|
|
|
|
64
то выражение (68) примет форму
£о +fb 'o M
совпадающую с |
выражением для |
спектра энергий |
квантового |
осцил |
|||
лятора с той, |
однако, |
разницей, |
что параметр |
у , |
играющий ь |
||
(70) роль "постоянной" |
Планка, |
как и "нулевая |
энергия" |
Ед о п |
|||
ределяются, как видно, |
с одной стороны, константами cvfffJa ,A г, |
||||||
а с другой - |
числом п07 т .е . |
выбираемым приближением. |
|
б) Фононы, дуализм, соотношение неопределенности
Поскольку для больших п сп я? са f можно условно ввести импульс полученного таким образом "классического кванта" как
2 Е* _ |
2 V |
С0 |
(71; |
* to |
т .е . трактовать его как возмущение с корпускулярно-волновым свойствами - "классический фонон". Взаимодействия этих класси ческих квазичастиц, как видно из (70) и (7 1 ), будет при этом подчиняться законам сохранения энергии и импульса
^' £м =
Поскольку |
каждый |
солифон представляет |
собой |
"облако" |
из |
|||||||
(2n0 - 2 N - 1 ) |
классических |
фононов |
(его |
энергия |
|
|
||||||
■= En ( 2 п - 2 Л/- 7 )), |
солифоны, |
очевидно, также обладают корпус |
||||||||||
кулярно-волновыми свойствами. |
Самый короткий солифон |
- |
это |
в о з |
||||||||
мущение, в котором в движении |
находится практически |
один |
узел |
|||||||||
цепочки, так что его распространение выглядит |
как |
движение пу |
||||||||||
стого |
узла |
- |
"дырки". |
|
|
|
|
|
|
|
|
|
|
Наконец, из (69) и (70) |
видно, что раз ища между нулевыми |
||||||||||
энергиями для двух разных приближений |
Ео (пд ) |
и |
|
Е0 (пд +1) |
||||||||
равна |
Ejp-0. Это та |
минимальная ошибка, |
с |
которой можно |
изме |
рить энергию рассматриваемых квазичастиц. С другой стороны,ми
нимальная |
ошибка |
в координате |
Л * mLn ~ |
>что соответствует |
|||
временам |
& |
= |
= |
. Возможные ошибки в |
определении энер |
||
гии А Е |
квазичастицы и времени ее прохождения |
через фиксиро |
|||||
ванную точку пространства |
A t |
связаны, |
..таким образом, клас |
||||
сическим |
аналогом |
соотношения |
неопределенности |
|
|||
|
|
|
A E A t |
* |
~ . |
|
(78) |
Зак.240 |
65 |
В рамках рассматриваемой |
теории линейный энергетический спектр |
||||||||||
вида (70) |
теряет законность при больших энергиях (N |
~ па ) , где |
|||||||||
уровни определяются |
формулой (6 6 ), |
дающей для модели |
кубическо |
||||||||
г о |
осциллятора (5) |
максимальное значение энергии |
Е = Етах~ |
||||||||
- |
|
. При колебаниях больших амплитуд в уравнение |
(35) еле,ду |
||||||||
ет, |
очевидно,вводить нелинейные члены более высокого |
порядка» |
|||||||||
|
Нелинейно-дисперсионная эволюция,описываемая уравнениями |
||||||||||
(35) и |
(4 4 ), в части |
взаимодействия |
встречных |
солифонов являет |
|||||||
ся |
задачей |
нерешенной |
(в |
литературе |
известны |
/?-солитонные |
ре |
||||
шения родственного уравнения u ft - ufX - 6 ( и г)хх - |
икххх - |
О |
|||||||||
[€ 8 ,6 4 ] |
) . |
Однако для моноориентированных возмущений |
процесс |
||||||||
становится кортевег-де-вризовским и на этом этапе картина |
рас |
||||||||||
пада ясна. |
|
|
|
|
|
|
|
|
|
§5. О нелинейной дифракции акустического пучка
втвердом теле
Ограничившись простейшим случаем изотропного упругого твердого тела, рассмотрим влияние на структуру пучка нелинейно сти тензора деформации (геометрической нелинейности).
Подставляя выражение тензора деформации
“ск |
= |
- |
дщ ^ |
дик |
. |
9и4 |
fa* |
\ |
|
|
дч |
|
dxL |
ч |
) |
|
|||
|
|
|
|
|
|
||||
в уравнение движения |
|
|
|
|
|||||
|
= |
4 k |
/ д \ |
> получим |
|||||
|
|
J>a r { |
|
дгч |
|
дл1/с |
|
||
|
|
3 / дхс дхг |
Xs |
|
|
||||
|
|
|
|
|
|
||||
|
|
|
|
Кк |
|
|
|||
|
|
|
|
|
|
|
|
дир |
диР |
|
|
|
|
|
|
|
|
dxL |
дхк |
Здесь |
к |
и |
{Л - модули |
сжатия и |
сдвига. |
|
|||
В случае отсутствия в волновом процессе сдвиговой |
|||||||||
ты может |
быть |
введен потенциал |
|
. |
после чего |
||||
ние (75) |
перепишется |
как |
|
L |
дхс • |
|
|||
|
|
|
|
|
(74)
(75)
компоненура вне-
dz<f |
д2(/ |
д2У (7 0 |
Я Г - ( к + У ) |
||
ч |
d*id*P |
дхк дкР |
Замечая, что
ik |
3z<f |
(77) |
|
3kl dxp dxk 3xp |
|||
# * t dxp.‘ |
|||
|
|||
|
|
и вводя оператор Даламбера, представим уравнение (76) в удобной для анализа форме
7 |
/ |
*■ |
1 |
/ |
11 ft ^ |
|
' V |
к |
\ |
/ |
д г(/ |
|
3 |
2 |
J\ |
dxt Зхр |
|
|
|
|
|
|
(78) |
1 |
3 г |
, |
|
г |
|
0 = 4 |
a t 2 |
|
с - у ! . - |
||
|
|
|
|
3j° |
Правая часть полученного уравнения даст в первом приближе нии члены, содержащие производные трех типов
Первая производная обусловит эффект самофокусировки |
- |
F ((/) , |
||||||||
где |
F (i/) |
- |
амплитуды пучка |
типа |
(162, г л .1 ), два |
|
остальных |
|||
типа |
слагаемых |
дадут дуплетное |
^ ( у ) |
и триплетное |
З'[у) |
рас |
||||
щепление |
пучка. |
|
|
|
|
|
|
|||
|
В заключение обратим внимание на |
один механизм, |
могущий |
|||||||
привести |
к |
выраженной самофокусировке |
акустического |
пучка |
- |
|||||
акустотермический эффект в сплавах |
типа |
нитинола, выделение |
т ы |
ла в которых под нагрузкой обусловлено структурными (мартенсит ными) превращениями.
ГЛАВА 111
НЕЛИНЕЙНЫЕ
ИДИСПЕРСИОННЫЕ ЯВЛЕНИЯ
ВГЕМОДИНАМИКЕ
Несмотря на длительный, почти трехвековой, период изучения [ульсовых осцилляций в кровеносных сосудах и огромное число как экспериментальных, так и теоретических работ в этой области, ограниченность числа точек соприкосновения между теорией и экс периментом остается основной трудностью гемодинамики.Эта труд ность в значительной мере обусловлена отсутствием сколько-нибудь
развитых нелинейных теорий |
пульсовых |
волн. |
Между |
тем |
уже |
из |
|||||||||||
предварительных |
оценок |
следует, |
что пульсации крови |
в |
сосудах |
||||||||||||
'осят существенно нелинейный характер. Действительно, |
согласно |
||||||||||||||||
опытным данным (например, С70 - |
72J) |
скорость |
пульсовой |
волны с |
|||||||||||||
7 |
человека |
колеблется |
в |
пределах |
|
от |
3 - |
4 м /с |
в аорте |
до |
|||||||
- |
14 |
м /с в |
артериях |
конечностей. Амплитуда |
же скорости пере- |
||||||||||||
'1ещения крови |
(форменных элементов) |
по сосудам |
v |
составляет |
|||||||||||||
жоло 60 см /с |
в |
аорте |
и около |
50 |
см /с |
в бедренных артериях.Чис |
|||||||||||
та |
Маха |
М—-jr, таким |
образом, |
лежат по |
этим данным |
в диапазоне |
|||||||||||
М |
I |
+ 10“ ^, |
что |
соответствует |
области |
существенно |
нелиней |
||||||||||
ных волн и делает линейное приближение малоприменимым. |
Послед |
||||||||||||||||
нее не означает, разумеется, что |
линейная и |
нелинейная |
теории |
||||||||||||||
должны расходиться по воем вопросам. |
|
|
|
|
|
|
|
|
|||||||||
|
В настоящей главе упругая трубка с жидкостью |
(кровеносный |
|||||||||||||||
сосуд , |
шланг и |
т ,п .) |
рассматривается |
как |
типичная кортевег - |
де- |
|||||||||||
вризовсквя нелинейная |
система |
со |
сложной, |
однако, |
инерционно |
геометрической |
дисперсией. Нелинейно-дисперсионная эволюция ге |
|
модинамического |
импульса, в частности, "солкптонный" |
распад, |
позволяют дать |
достаточно правдоподобную качественную интерпре |
|
тацию пульсовых осциллограмм и ряда эксперимента лысых |
э'№ектсв |
|
гемо,динамики. |
|
|
6Ь
§ I . Линейные одномерные волны
а) Инерционная дисперсия
Пусть сосуд представляет собой трубку с внутренним радиу
сом R0 , |
толщиной r0 ( R0 > > |
г0 ) |
и |
. плотностью материала р1 |
|
(р и с.22), |
плотность жидкости |
в сосуде |
p gj, Для |
длинных волн.' |
Р и с . 22
(л >> ) радиальной составляющей скорости течения можно прене бречь, считая осевую составляющую = v (х , t) однородной по се чению, и записать уравнение движения в линейном приближении
^ + ± М . |
— Р t |
( I ) |
d t р 0 дк |
|
|
где Д - вариация давления, усредненного по сечению трубки. Ра диальное смещение стенок трубки <$ задается уравнением непре рывности
|
|
3 . h + Ж = о , |
( 2) |
|||
|
|
2 |
дх |
d t |
|
W меж- |
которое вытекает из соотношения для |
изменения объема |
|||||
ду сечениями |
х |
и |
х + Ых: |
|
|
|
д\Р |
|
at |
dx |
3v_ |
|
|
3 t |
= 2 |
JlRn dt |
|
dx dx- |
|
|
Уравнение |
состояния |
|
запишем с учетом |
инерции |
трубки в виде уравнения радиальных колебаний элементарного объ
ема с |
сечением s |
= / |
и высотой |
г0 : |
|
|
|
$ + 4 |
$ - |
Р, Г0 Р |
сог |
я/ А - |
О ) |
|
о |
|||||
Здесь |
^ - упругая |
константа трубки, |
|
параметр, имеющий |
смысл частоты собственных колебаний стенки. Тангенциальными на пряжениями и смещения?.™ стенки в одномерном приближении следует,
очевидно, пренебречь.. Подстановка р