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

книги из ГПНТБ / Чугаевский, Ю. В. Элементы теории нелинейных и быстропеременных волновых процессов

.pdf
Скачиваний:
12
Добавлен:
19.10.2023
Размер:
7.46 Mб
Скачать

Вид функции (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“ ^,

что

соответствует

области

существенно

нелиней­

ных волн и делает линейное приближение малоприменимым.

Послед­

нее не означает, разумеется, что

линейная и

нелинейная

теории

должны расходиться по воем вопросам.

 

 

 

 

 

 

 

 

 

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

(кровеносный

сосуд ,

шланг и

т ,п .)

рассматривается

как

типичная кортевег -

де-

вризовсквя нелинейная

система

со

сложной,

однако,

инерционно­

геометрической

дисперсией. Нелинейно-дисперсионная эволюция ге­

модинамического

импульса, в частности, "солкптонный"

распад,

позволяют дать

достаточно правдоподобную качественную интерпре­

тацию пульсовых осциллограмм и ряда эксперимента лысых

э'№ектсв

гемо,динамики.

 

 

из (2) в (I ) дзет

§ 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 Р

сог

я/ А -

О )

 

о

Здесь

^ - упругая

константа трубки,

 

параметр, имеющий

смысл частоты собственных колебаний стенки. Тангенциальными на­ пряжениями и смещения?.стенки в одномерном приближении следует,

очевидно, пренебречь.. Подстановка р

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