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

книги / Основы механики сплошной среды

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
13.75 Mб
Скачать

 

Итак, после введения отсчёт-

 

ной

и

актуальной конфигураций

 

можно определить движение тела

 

как

отображение* отсчётной кон­

 

фигурации в актуальную, т. е. за­

 

ниматься только

отображениями

 

“фотографий” тела, а не самим

 

телом

(рис. 3).

В силу гипоте­

 

зы

непроницаемости отображе­

 

ние (1.4) будет биективным.

 

 

Место, занимаемое частицей X в отсчётной конфигурации,

описывается радиусом-вектором TQ:

 

 

 

П) =

Го(£1,Ё2.£з»*о) = Я?(£ь£2.£з.*о)£г .

(1-5)

а в актуальной конфигурации — радиусом-вектором г:

 

 

 

?= ^ (£ь&.£з>*) = Яг(£ь£2.£з.*)Ач ,

(1.6)

где xi — эйлеровы координаты или координаты места в ящи­ ке Я, занимаемого частицей, которая в момент t = to занимала место с координатами xf. Из (1.5), ( 1.6) следует, что соотноше­ ние (1.4) можно записать в виде

Xi —Xj(xÇ, 3?2,Х3,t) ИЛИ f= f(fo ,t) .

(1.7)

В соотношениях (1.5) и (1.6), как и всюду далее в книге, используются следующие общепринятые правила суммирования [29, 55]:

а) индекс, изображающийся малой латинской буквой, изме­ няется от 1 до 3;

б) индекс, изображающийся большой латинской буквой, из­ меняется от 1 до 2;

в) латинские индексы могут встречаться в каждом одночлене либо один, либо два раза; если индекс встречается два раза, он называется немым и по нему производится суммирование от 1 до 3 (если он изображается малой буквой) или от 1 до 2 (если большой), причём для краткости знаки суммы опускаются; если индекс встречается один раз, он называется свободным (во всех одночленах данной формулы свободные индексы должны

совпадать);

 

 

 

 

 

г) индексы,

обозначающиеся греческими буквами,

могут

встречаться

в

каждом

одночлене произвольное

число

раз,

и по ним

суммирование

не производится (если,

разумеется,

специально не написан знак суммы); разным греческим буквам в индексах в данной формуле обязательно соответствуют разные числовые индексы.

Так, например, компактная запись

<Ч Л ьЬ акС аЬ = fc d J

эквивалентна системе восемнадцати уравнений (каждое из сле­ дующих шести уравнений надо взять при а = 1,2,3):

O llllb alC al

+ û n i2 ^ Q lco2 + a 112lba2Cal +

 

 

 

 

 

 

+ a il2 2 & a 2 C a 2 +

а П З |Ь а З с а1

+

« 1 1 3 2 ^ 3 ^ 2

f a ll '.

Û 1211& alca l

+

°1 21 2b a lC a 2 + ® 1221^a2C al +

 

 

 

 

 

 

+ O l222&a2Ca 2 +

а 1231^а 3с а1

+

ai2 32 & a 3 C a 2

=

/ a l 2 i

a 2l l l ^ a l ca l

+

Û2l 12^01 Ca 2 + 02121^

2 ^ 1 +

 

 

 

 

 

 

+ a 2122^a2ca2 +

a 2131^a3CQl + « 2 1 3 2 ^ 3 ^ 2

=

/a 2 b

“ 2211 ^ a l C al

+

0 2 2 12 ^ 1 Ca 2 + a 222lba 2Cal +

 

 

 

 

 

 

+ a2222^a2cce2 +

0 2 2 3 lb a 3C al

+

Û 2232ba3c a2

=

f a 22Î

a 311lbalCal +

Û3112balCa 2 + a3 1 2 lb a 2Cal +

 

 

 

 

 

 

+ Û3122ba 2ca2 +

®313lba3ca l

+

a 3132ba 3CQ2 =

f a 3b

Û321lbalCa l +

032i2balCa 2 + Û322lba 2Cal +

 

 

 

 

 

 

+ Û3222ba 2Ca 2 +

a 323lb a 3 ca l

+

Û3232ba 3ca 2

=

/a 3 2 -

Сравнение двух предыдущих формул, без сомнения, убедит чи­ тателя в преимуществах использования тензорной алгебры.

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

dxi

dxi д£к

дх)

dx°j

соотношения (1.7) можно обратить следующим образом:

X i = (Zi, Х 2 , я3, £) или fo = г0(г, <).

(1.9)

Если в отсчётной конфигурации зафиксировать материаль­ ную частицу (например, имеющую лагранжевы координаты х? = = Cf), то из закона движения (1.7) получим уравнение линии:

Xi = Xi(Ci,C2,C3,t).

(1.10)

Эта линия называется траекторией частицы.

Если теперь положим в отсчётной конфигурации

 

х% = Сч = const, Х3 = С3 = const,

( 1.11)

т. е. зафиксируем прямую линию, составленную из материальных частиц, параллельную оси координат (Ох®), то в актуальной конфигурации в фиксированный момент времени t = T = const

получим

Л

(1.12)

 

xi = xi{x°l,C2,C3,T).

Уравнения (1.12) описывают некую кривую, составленную из ма­ териальных частиц, которые в отсчётной конфигурации лежали на прямой (1.11).

Если теперь положим в отсчётной конфигурации

х° = Сз = const,

(1.13)

т. е. зафиксируем плоскость, составленную из

материальных

частиц, параллельную координатной плоскости

(Ох®х®), то

в актуальной конфигурации в фиксированный

момент времени

t = T const получим

 

Xi = Хг(х°,Х2,С з,Т).

(1.14)

Уравнения (1.14) описывают поверхность, составленную из ма­ териальных частиц, в которую перешла при движении сплошной среды плоскость (1.13).

Из (1.11), (1.12) следует, что ортонормированный базис, составленный из материальных частиц в отсчётной конфигура­ ции, вообще говоря, превращается в актуальной конфигурации в “криволинейный” базис [36].

Заметим, что в силу (1.8) можно обратить соотношения (1.5) для мате­

риальных координат

 

& = &(х?,х5,хз,£0)-

(1.15)

Введём в рассмотрение вектор пе­

ремещения как разность

векторов г

и го (рис. 4):

 

U = f —fçj.

(1.16)

Компоненты щ в базисе к{ можно считать функциями лагранжевых координат и времени: щ = щ {х\,хjj.xjj.t).

Продифференцируем обе части закона движения (1.7) по времени и определим вектор скорости v с компонентами Vi

в базисе ïq, являющимися функциями лагранжевых координат

и времени: Й

ч й 4 4 ‘) = = ( 4 .4 4 *) = ^ ( 4 Л 4 Л

или v(fa, t) = ^ ( Р 0, i) = ^?(?о, () = ^ ( ?«. *) = ^ { ^ 0. t)-

В (1.17) полные производные по времени совпадают с частными, поскольку лагранжевы координаты х® от времени не зависят.

Продифференцируем обе части (1.17) по времени ещё раз и определим вектор ускорения w с компонентами го{ в базисе ki, являющимися функциями лагранжевых координат и времени:

^(а^ .хз.хз,*) = ^-(x?,xj),x§,i) = ^p(x?,x§,X 3,f)>

t

02-

Q2 -♦

(l- 18)

ИЛИ w (fo , t ) = - £ ( r 0 , *) =

0 =

g f i f ô ,

<)•

Компоненты тех же самых векторов перемещения, скорости и ускорения в силу (1.9) можно представить и как функции

эйлеровых координат х* и времени:

 

 

 

^(xj.xjbX jj.t) =Ui(xi,X2,x3,t);

(1.19)

^(xpXg.Æg.i) = Vi(xi,x2,xs,t);

( 1-20)

Wi(x°x,X*2, Æ3, t) = Wi(xi,x2,x 3,t).

(1.21)

Векторы w(r,t) и v(r, t) связаны между собой следующими

формулами:

dv

dv dxi

dv

dv

 

dv

( 1.22)

w = dt = m

+ dx{ dt

dt + dxi4 .

 

или покомпонентно:

 

 

 

 

 

dvi

dvi

, dvi dxj

dv:

dvi

(1.23)

S

= W

+ dxj dt

dt

+ d xjVj‘

 

В (1.22) и (1.23) видно различие между полной производной и частной производной по времени. Для любой физической величи­ ны F (x\,x2,x3,t), зависящей от эйлеровых координат и времени,

полная производная по времени (её также будем обозначать точкой справа вверху буквы) представима в виде

F (xi,x2, x3,t)

dF

 

= — (xi,x2,x3,t) =

 

 

dF

dF

= J ft(x i 'x*'x3't) + far(x i’xï ’x3’t) vi' (1-24)

т. е. в виде суммы частной производной по времени и конвек­ тивной производной по времени.

В зависимости от того, какие координаты — лагранжевы я® или эйлеровы Xi — выбраны в качестве независимых перемен­ ных, различаются два подхода к описанию движения сплошной среды, связанные с именами Ж.Л. Лагранжа и Л. Эйлера [45]. При лагранжевом подходе целью является нахождение закона движения (1.7), т.е. траекторий всех частиц тела. Лагранжево описание деформирования сплошной среды применяется чаще

вмеханике деформируемого твёрдого тела, где удобно следить за движением границы тела.

При эйлеровом подходе целью является нахождение поля вектора скорости v как функции г и t. Эйлерово описание применяется в случаях, когда необходима информация о том, как изменяется со временем та или иная физическая величина

вданной точке пространства (ящика Я). Такой подход чаще используется в гидроаэромеханике при исследовании жидкости [10].

Оба подхода эквивалентны. Действительно, зная закон дви­ жения (1.7), можно найти согласно (1.17) и (1.18) векторы скорости v(ro,t) и ускорения w(ro,t) и, воспользовавшись со­ отношениями (1.9), в силу (1.20), (1.21) получить векторные поля v(f, t) и w(f, t).

Пусть теперь имеются три функции Vi(xi,x2,X3,t). Выпи­

шем систему дифференциальных уравнений относительно функ­

ций Xj(£):

^

 

 

= Vi(xi,x2,xz,t)-

(1.25)

Решая систему (1.25) с начальными условиями я,(0) = я° (зада­ чу Коши), находим закон движения я^ЯрЯ^Яд, t).

Рассмотрим ниже три характерные задачи. Пусть имеется закон движения сплошной среды:

(1.26)

где а — постоянная. Надо найти поле скоростей ^(ац.яг.хз.Ё) в эйлеровом пространстве.

Заметим, что ii(0) = х\ и, кроме того,

 

 

1

at О

1 +

< А 2 > о (1.,27)

щ

—at 1

О

О

0

1

 

 

т.е. система (1.26) действительно представляет собой закон дви­ жения.

Дифференцируя (1.26) частным образом по f, получим

 

t)| = ai® . ^2 =

—a x j , из = 0 ,

(1.28)

а обращая (1.26), будем иметь

 

 

 

_0 _

i l

—a ti2

 

X l~

l + a 2t2 ’

 

< J ) _ x 2 + atxl

(1.29)

*2 _

1 + 02t2 '

 

^ ^3

*з-

 

Подставим обратный закон движения (1.29) в (1.28) и выпишем решение задачи:

v\ = аХ2 + ОХХ J

 

I + аЧ2

 

atxo —xi

(1.30)

щ — а 1 + a2t2

 

из = 0.

 

Как следует из (1.28)iи (1.30), вид функций

üi(i®,i®,i®,t)

и V{(x\,X2,xz,t) существенно различен (хотя эти функции и обо­ значаются одной буквой).

Пусть теперь дано другое поле скоростей Vi{x\,X2,x^,t)

в эй­

леровом пространстве:

 

V \= U )X 2 , V2 = U)X[ , v3 = v 0 ,

(1.31)

где ÜJ и vo постоянны. Требуется определить закон движения частиц.

Для этого необходимо решить задачу Коши для трёх уравне-

dxi

 

d i2

= tUI1,

dx3

(1.32)

dt = UX2

,

dt

dt — v0

с начальными условиями хДО) = i®. Следствием (1.32) является система

(fxi

о

$ X 2

о

dx3

 

dt2

= - a r i i ,

dt2 =

 

~ d t= v°’

(1.33)

Выпишем её общее решение:

 

xi =

С\ cosuit + С2 sinuit,

 

 

х2 Сз cos u)t + С4 sin ut,

(1.34)

 

a* = Щt + C5.

 

Подставим теперь

решение (1.34) в уравнения (1.32) и начальные

{

 

 

условия, тем самым связывая константы C i, . . . ,65

и х°. В ре­

зультате получим закон движения

 

 

' xi =

Х[ cosut + х\ sin art,

 

 

< X2 =

—x [ sin ait + x® cos ut,

( 1.35)

, X 3 = x° + VQL

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

скоростью Vo и постоянной угловой ско­

ростью и

(рис. 5).

Ось

же хз движет­

ся вдоль

самой себя со

скоростью vo.

В третьей задаче имеется поле скоростей

Vi(x1,x2, 13, t) в эйлеровом пространстве:

V\ = СХ2 —6x3 ,

V2 = ОХз —СХ1,

 

V3 = 6x 1

—а х г,

(1.36)

 

 

где а, Ь, с — некоторые постоянные. Необхо­ димо найти траектории частиц.

Исследуем систему уравнений задачи Коши, возникающей в данном случае:

dxi

СХ2

,

 

 

^

0X3»

 

 

dx2

= 0x3 —сх 1,

- Л

(1.37)

dt

Xi(0) =

,

 

 

 

dxz

ах2,

 

 

= bxi -

 

 

at

 

 

 

 

Умножим первое уравнение (1.37) на а, второе на Ь, а третье на с и просуммируем:

dx\ ( и dx2 , dx3

Интегрируя (1.38) по времени с учётом начальных условий:

ах1+ Ьх2 + сжз = ах® + bx® + ex®,

( 1.39)

придём к выводу, что каждая частица, имевшая в начальный момент координаты Xj, х®, х®, будет оставаться в плоскости, определяемой (1.39), т.е. траектория каждой частицы является плоской.

Умножим далее первое уравнение (1.37) на х ь второе на хг, а третье на хз и просуммируем:

dx1

dx2

dx3

_

1

dx\

dx\

= 0. (1.40)

Xl~ d t+X2 dt

+ x3 dt

~

2

dt

dt

Интегрируя (1.40) по времени и опять же учитывая начальные условия, будем иметь

х\ + х\ + х | = (х?)2 + (х§)2 + (х§)2 ,

(1.41)

т.е. каждая частица, имевшая в начальный момент координа­ ты х®, х®, Хд, будет оставаться на сфере, определяемой уравне­ нием (1.41). Соотношения (1.39)

и (1.41) называются первыми инте­ гралами системы (1.37).

Итак, траекториями будут окруж­ ности (для каждой частицы своя), являющиеся пересечением плоскос­ ти (1.39) и сферы (1.41) (рис. 6). Ис­ ключение составляют лишь частицы, находившиеся при t = to на прямой I, которая проходит через начало коор­ динат и параллельна вектору {а;Ь\ с} (эти частицы, очевидно, будут нахо­ диться в покое).

Отметим некоторые часто встречающиеся законы движения сплошной среды.

а) Трёхмерное растяжение-сжатие:

Ха = (1 + аа (4))Хц, « = 1 , 2 , 3 ,

(1-42)

причём а*(0) = 0 и aj(t) ф 1. Соответствующее поле скоростей v(f, t) имеет вид

0>аХа

а = 1,2,3.

(1.43)

Va = 1 + Ûa

Если аз(£) = 0, то говорят о двумерном растяжении-сжатии

в плоскости (Ох\Х2). Если аз(£) = аг(£) = 0, то имеет место

одномерное растяжение-сжатие вдоль оси (Oxi).

б)

Одномерный сдвиг:

 

 

 

 

х\ = xÇ + û(£)xj>,

Х2 = х \ , хз = х®,

(1-44)

причём а(0) = 0. Поле v(f, t) имеет вид

 

 

vi = àx2,

v2 = V3 = 0.

(1-45)

в)

Однородное состояние:

 

 

 

 

Xi = A ij(t)xj,

(1.46)

причём

-4v,'(0) = Sij] det A(£) ф 0; Sÿ

символы

Кронекера

 

Г 1,

если

г = j,

(1.47)

 

\ 0 ,

если

г ф j.

 

 

Для скоростей v(f, t) справедливы соотношения

 

 

Vi = Bikxk ,

В = À • А - 1 ,

(1-48)

где А -1 — обратная к А матрица, т. е. А • А -1 = А -1 • А = I; I — единичная матрица, компонентами которой в декартовой системе координат являются Sij.

Л Е К Ц И Я 2

ЭЛЕМЕНТЫ ВЕКТОРНОГО АНАЛИЗА

Векторной линией данного векторного поля называется та­ кая кривая, в каждой точке (х\,Х2,х$) которой её касается вектор а(хьХ 2,хз), принадлежащий этому векторному полю

Рис. 7

(рис 7). Если a(x],X2,x$,t) явно зависит от времени, то картина векторных линий будет со временем меняться. Для векторного поля скоростей векторные линии называются линиями тока.

Если через точку пространства проходит более одной линии тока, то эта точка называется особой. Примеры таких особых точек изображены на рис. 8. Проведём в момент времени t через некоторую неособую точку М с радиусом-вектором г линию тока (рис. 9). Вдоль этой линии выберем естественный параметр

Рис. 8

Рис. 9