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

Руководство к решению задач

.pdf
Скачиваний:
268
Добавлен:
02.03.2016
Размер:
647.64 Кб
Скачать

41

ω

 

 

 

 

( ω x

 

− ω x ) − ω

 

 

 

 

 

 

 

( ω x − ω x

 

 

) − ω

 

 

 

( ω x − ω x

 

) +

2 x

 

 

2

 

 

x

 

 

 

3

 

x

3

 

3

 

1

 

2 1

 

3

 

3

 

3 1

 

1

 

 

1

 

3 1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

(ω x −ω x ) −ω

 

 

 

 

 

 

(ω x −ω x ) −ω

 

 

 

(ω x −ω x ) +

 

 

 

 

 

 

 

x

 

 

x

 

 

 

 

1 x

 

3 1 1 3

2

 

 

 

2 3

3 2

 

 

3

 

2 3

3 2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ω3ω1 − ω1ω3 = 0 ,

 

 

 

 

 

 

 

 

 

 

ω

 

 

 

( ω

x

 

 

− ω

x

 

) − ω

 

 

 

 

 

( ω x

 

− ω

x ) − ω

 

 

 

 

( ω x

 

− ω x )

 

 

 

 

 

3

2

 

 

 

 

 

 

 

2

2

x

 

 

2

3 x

 

 

 

2

 

3

 

 

1 x

 

 

 

 

1

 

2 1

 

 

 

2

1

2 1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− ω

 

 

 

 

 

 

 

( ω x − ω x

 

 

) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 1

1

3

 

 

 

 

 

 

 

 

 

 

 

= ω1ω2 −ω2ω1 = 0 .

Таким образом, rotaoc = 0 .

Если использовать символические формулы, то получим этот же результат. Покажем это

×(ω×r ) = ×(ωc ×V ) + ×(ω×V c ) =

=ωc ( V ) V( ωc ) + ω( V c ) V c ( ×ω) = = ωdivV V( ω ) + ω(V ) V c divω) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ω( vk

 

) V( ωk

 

) = ωk vk

 

ek vk ωk

 

ek = 0 ,

xk

xk

xk

xk

 

 

 

 

 

 

 

 

 

так как divV = 0 и divω = 0 , т.е. rotaoc = 0 , а вихрь поля ускорений будет равен

rota = rotaвр = 2ε.

3.13. Показать, что если поля a и b суть безвихревые поля, то поле a × b является соленоидальным.

Решение.

Согласно условию задачи rota = 0 и rotb = 0 , а согласно определению соленоидального поля его дивергенция должна быть равна нулю. Вычислим div( a ×b ). Это лучше сделать символически

div( a×b) = ( a×b ) = ( ac ×b) + ( a×bc ) = = −a ( ×b ) +b ( ×a ) = −a rotb +b rota = 0,

так как rota = 0 и rotb = 0.

Следовательно, векторное поле a × b будет соленоидальным.

42

3.14. Доказать, что поле b = ϕ ϕ, где ϕ − скалярная функция, является безвихревым, т.е. rotb = 0 .

Решение.

Доказательство проведем в символической форме

rotb = ×b = ×ϕ ϕ=ϕ( ×ϕ) rot grad ϕ=0.

Задачи для самостоятельного решения

 

3.15. Задан закон движения сплошной среды

x = A +

ebλ

sin λ( A + ωt ),

 

 

 

 

 

 

1

λ

 

 

 

eBλ

 

 

 

x2

= −B

cos λ( A + ωt ), x3 = X 3 . Показать,

что траектории частиц ок-

λ

 

 

 

 

 

 

ружности, а величина скорости постоянна. Определить также связь между X1 и

X 2

и константами A и B .

 

 

 

3.16. Доказать, что поле скоростей

V =( Ax3 Bx2 )e€1 +( Bx1 Cx3 )e€2 + (Bx1 Cx3)e2 + (Cx2 Ax1)e3

представляет вращение абсолютно твердого тела, так как для него D 0 .

3.17. Для вращения абсолютно твердого тела со скоростью

V = 3x3e€1 4x3e€2 +( 4x2 3x1 )e€3 определить вектор вихря скорости ω и пока-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зать, что V

= ω×

r

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.18. Дано стационарное поле скоростей

 

 

 

 

 

 

 

 

 

 

 

= ( x3 x x2

)e€

+( x2 x

 

+ x

 

)e€ .

Найти

скорости

относительно

точки

V

2

2

1

 

1

2

1

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

M (1,1,3 ) частиц в точках N1

(1,0,3 ), N2 (1,

3

,3 ), N3

(1,

7

,3 ), отнесенные к рас-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

8

 

 

 

 

 

стоянию от этих точек до точки M , и показать, что их величины стремятся к

 

 

 

 

 

 

 

относительной скорости, определенной формулой dV

=V

dr

.

 

3.19.

 

Доказать,

что

при

установившемся

движении среды,

когда

v

( ti = 0 ) линии тока и траектории совпадают.

43

3.20. Доказать, что для течения vi = 1x+i t линии тока и траектории сов-

падают.

3.21. Доказать, что поле скоростей vi = εijk bj xk + ci , где bi и ci – постоянные векторы, представляет вращение абсолютно твердого тела, и определить вектор вихря скорости ωi для этого движения.

Ответ: 2ωi = qi = 2bi .

3.22. Доказать, что для поля скоростей v1 = x12 x2 + x23 , v2 = −x13 x1 x22 , v3 = 0 линии тока будут окружностями.

44

Раздел 4

Динамика сплошной среды

Закон сохранения массы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

(δm)

=

 

 

d

 

 

(ρδτ)=

dρδτ

d

(δτ)= 0.

(4.1)

 

dt

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

dt

 

 

 

Уравнение сплошности (неразрывности)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dρ

 

divV

 

 

= 0 .

 

 

 

 

 

(4.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Другой вид уравнения неразрывности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

) = 0.

 

 

 

 

 

(4.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ div( ρV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В проекциях на координатные оси x1 ,x2 ,x3

 

 

 

 

 

 

ρ

+

 

 

(ρv

 

 

)+

 

 

 

 

 

 

 

(∂v

 

 

 

)+

 

 

 

v

 

)= 0.

(4.4)

t

x

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

3

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для несжимаемой жидкости (ρ = const) уравнения (4.2) и (4.4) примут

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divV

 

 

 

 

= 0, или

 

 

 

vi

 

 

= 0.

 

 

(4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема об изменении

 

 

 

главного

 

 

 

вектора

 

количества

движения

''жидкого'' объема

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

dV

 

δτ =

 

 

 

 

 

 

δτ+

pn€δs .

 

 

(4.6)

 

 

 

 

 

ρf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

Дифференциальные уравнения динамики сплошной среды ''в

напряжениях'':

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

dV

 

 

 

 

 

 

 

 

+ DivP,

 

 

 

 

 

(4.7)

 

 

 

 

 

 

 

 

 

 

 

= ρf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или с учетом индивидуальной производной

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+(V

 

 

 

 

 

 

 

 

 

 

 

 

= ρf + DivP ,

(4.8)

 

 

ρ

 

 

 

 

 

V )V

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в аналитической форме

45

 

vi

 

 

 

 

pki

 

 

 

 

vi

= ρ fi +

( i =1,2,3 ).

 

ρ

 

+ υk

 

 

 

(4.9)

t

 

xk

 

 

xk

 

 

 

Теорема об изменении момента количества движения ''жидкого'' объема

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

×ρ

dV

δτ =

r

 

 

 

 

 

 

 

δτ+

 

 

 

 

 

×ρ f

r

× pn€δs .

(4.10)

dt

τ

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема об изменении кинетической энергии ''жидкого'' объема

 

 

 

 

 

 

d

 

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

= ρ f V

 

 

+V DivP .

(4.11)

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Другая форма этой теоремы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ρ f V + div( PV ) P D .

 

ρ

 

 

 

 

 

 

(4.12)

 

 

 

 

 

 

 

 

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Объемная плотность распределения мощности внутренних сил по среде

 

Nвн

 

 

 

 

 

 

 

) = −P D.

(4.13)

 

=V

DivP div( PV

Уравнение баланса полной энергии в дифференциальной форме

 

 

d

 

 

v2

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

= ρ f V

+ div( PV ) q .

(4.14)

 

U +

 

 

 

 

dt

 

 

2

 

 

 

 

 

 

 

 

 

 

 

Уравнение баланса внутренней энергии

 

ρ

dU

 

= ρq Nвн = ρq + P D .

(4.15)

dt

 

 

 

 

 

Теорема Бернулли

 

 

 

 

 

 

 

v2

+ ρ + Π = const.

(4.16)

 

 

2

 

 

 

 

Примеры решения задач

 

 

 

 

 

 

 

 

 

 

A( x2

x2

)

 

 

 

 

4.1. Дано плоское течение несжимаемой жидкости

v =

1

2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

r 4

 

 

 

 

 

2Ax1x2

 

 

 

 

 

 

 

 

v

2

=

, v

3

= 0, где r 2 = x2

+ x2

. Доказать, что поле скоростей удовле-

 

 

 

r 4

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

творяет уравнению неразрывности и является безвихревым.

46

Решение.

Для несжимаемой жидкости уравнение неразрывности имеет вид

 

 

 

 

 

 

v

 

 

v

2

 

 

v

3

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

x

2

x2

 

2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divV =

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

2

 

 

1

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

= 0 .

 

В

нашем

случае

 

=

A

4x

 

 

 

 

+

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

4

 

 

 

 

 

 

x1

 

 

x2

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

1

 

 

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

2x

 

 

 

 

 

 

x2

 

 

 

v

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1

8x

 

2

 

 

и

 

= 0. Складывая, получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= A

 

4

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

r

 

 

 

 

1

r

 

 

 

x3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4x3

 

 

 

8x x2 + 4x r 2 )A

 

 

 

 

 

- 4x1

x2 x2

+

 

2x

 

2x

8x x2

=

+ 4x x2

=

 

A

1

6

 

 

 

 

4

 

 

 

+ A

 

4

 

6

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1 2

 

1 2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

= rA6 (4x13 4x1 x22 + 4x13 + 4x1 x22 )= 0. Таким образом, уравнение неразрывно-

сти удовлетворяется.

Для безвихревого течения необходимо выполнение условия rotV = 0. Проверим его вычислением определителя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e€1

 

 

 

 

 

 

 

 

 

e€2

 

 

 

 

 

e€3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(x2 x2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Ax x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rotV = ×

r

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

1

2

 

e€

=

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

x

 

 

x

 

 

r 4

 

 

x

 

 

r 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(x2

 

1

 

 

 

2 )

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

3

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

2Ax x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 4

 

 

 

 

 

 

 

 

 

 

 

 

 

r4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2x2

 

 

 

 

 

 

2

x2

 

 

 

2x2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

4x2 r

2

 

 

 

2

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

8x1

 

 

 

 

 

 

 

 

 

4x2 (x1

 

x2 )

 

 

 

 

 

 

 

 

 

 

8x1

x2 + 4x1

x2 4x2

 

 

 

 

 

 

= A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e€3 =

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e€3

 

=

 

r

4

 

 

 

 

 

r

6

 

 

 

 

 

 

r

4

 

 

 

 

 

 

 

 

 

 

 

r

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

4x2 (x1

+ x2 )4x1 x2

 

4x2

e3

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следовательно, движение безвихревое.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.2.Определяющие

уравнения

для

некоторой

среды

имеют

 

 

 

вид

p

i j

= (p + λ* D

kk

)δ

ij

+ 2µ* D , где λ* и µ*

 

– постоянные. Написать уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

движения, выраженные через вектор скорости V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения движения ρv&i

 

= ρfi +

 

 

в этом случае будут иметь форму

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρv&

 

= ρ f

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

*

 

 

 

 

 

 

D

kk

 

 

 

 

2µ

*

 

Dij

.

По определению 2D

 

v

i

 

 

 

vij

 

,

 

 

 

−δ

 

 

 

 

 

 

 

 

+ λ δ

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

=

 

 

 

+

 

 

 

 

 

 

 

 

ij xij

 

 

 

 

 

 

x j

 

 

 

 

x j

xij

 

xi

 

 

 

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

47

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2D )

 

 

2

 

 

2v

 

 

 

 

так

что

 

Dkk

=

vk

=

vi

 

 

 

и

 

ij

=

 

vi

+

 

 

j

.

Поэтому

xk

xi

 

 

 

 

 

 

 

x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2j

xi x j

 

 

ρv&i

= ρfi

p

+ (λ*

* )

 

2v j

*

2vi

 

. Это уравнение можно записать в

 

 

 

x j x j

 

 

 

xi

 

 

xi x j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p + (λ*

+ µ* ) ( V

)* 2V

.

символических обозначениях ρV&

= ρf

4.3. Вычислить материальную скорость изменения кинетической энергии среды, занимающей объем τ, и объяснить физический смысл полученных интегралов.

Решение.

T

 

Согласно формуле

= ρvi v&i dτ,мощность внешних поверхностных

 

dt

τ

сил равна vi pni ds, что можно записать иначе: vi pij n j ds.Применим к этому

s s

интегралу теорему Гаусса – Остроградского и воспользуемся уравнением

движения.

Тогда

vi pij ni ds = pij

vi

dτ+ ρ(vi v&i fi vi )dr. Отсюда получим

xi

 

dT

 

 

s

vi

τ

τ

 

= ρfi vi dτ− pij

dτ+ vi pni ds .

Интегралы этой суммы представляют

 

dt

x j

 

τ

τ

 

s

 

 

мощности массовых сил и внешних поверхностных сил соответственно.

4.4. Доказать, что для движения абсолютно твердого тела с полем скоростей vi = εijk ωj xk интеграл кинетической энергии сводится к аналогичному

выражению из динамики твердого тела.

Решение.

Воспользуемся формулой интеграла кинетической энергии, в которую поставим значения скоростей данного движения

T = ρ

vi vi

dτ =

1

ρεijk ωj xk εipq ωp xq dτ =

1

ρωp ωj (δip δkq −δjq δkp )xk xq dτ =

2

2

 

τ

 

τ

2

τ

 

 

 

 

=

 

ωj ωp

ρ(δ jp xq xq x p x j )dτ =

ωj ωp I jp

,

 

 

 

 

 

 

 

 

 

 

 

 

 

2 t

2

 

где I jp

= ρ(δip xq xq

x p x j )dτ – тензор моментов инерции.

 

τ

 

 

 

 

 

 

 

 

 

 

4.5. В некоторой точке сплошной среды даны тензоры скоростей деформации и напряжения соответственно

 

 

 

48

 

 

 

1 4 5

 

 

 

4 0 1

 

4 2 6

 

и

 

0 5 7

 

D =

 

P =

.

 

5 6 3

 

 

 

1 7 6

 

 

 

 

 

 

Определить в этой точке величину λмощности напряжений.

Решение.

Умножая каждую компоненту тензора D на соответствующую компоненту тензора P , а затем складывая полученные произведения, получаем

λ = D P = 4 + 0 5 + 0 +10 + 42 5 + 42 +18 =106.

4.6. Найти функцию давления P(p) при баротропном течении жидкости с уравнением состояния p = λρk , где λ и k – постоянные величины.

Решение.

По определению функции давления имеем

 

 

 

1

 

1

 

 

 

k 1

 

 

 

 

 

 

 

 

 

p dp

p

p

 

 

kλ

k

 

 

 

p

k

 

p

 

p

 

 

k

 

 

0

 

 

 

 

 

 

 

P(p)=

 

=

 

 

 

 

 

 

 

 

 

 

k

=

 

 

 

 

 

 

 

 

dp =

 

 

 

 

p

 

 

 

 

 

ρ

k 1

 

 

k 1

 

ρ

ρ

 

.

p0

p0

λ

 

 

 

p

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

4.7. Газ с уравнением состояния p = λρk , где λ и k – постоянные ве-

личины, вытекает из большого закрытого резервуара через гладкую тонкую трубку; течение баротропное. Давление в резервуаре равно N атмосферам. Определить скорость истечения газа, считая давление в струе на выходе из резервуара равным атмосферному.

Решение.

Применим интеграл Бернулли к установившемуся течению в двух точках линии тока газа: A – в покоящейся внутри резервуара и B – на свободной поверхности вытекающей струи

 

v2

 

v2

 

 

ПA + pA +

A

= ПB + pB +

B

= const.

(1)

2

2

 

 

 

 

Но vA = 0, и если пренебречь силой тяжести и учесть, что течение баротропное, т.е.

 

p

 

 

p

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp

 

 

p k

 

kλk

 

k 1

p

 

k

 

p

 

p

 

 

 

p =

 

=

 

 

 

dp =

 

 

p

 

 

=

 

 

 

 

0

 

,

(2)

ρ(p)

 

k 1

k

 

k 1

ρ

 

 

 

 

λ

 

 

p

 

 

 

ρ0

 

 

 

p

 

 

p

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где p0 и ρ0 – значения величин, относящиеся к какой-нибудь точке покоя-

щегося газа, а k = c p – показатель адиабаты, то уравнение (1) дает cV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

49

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

pA

 

 

 

 

 

 

 

 

1

2

 

 

2

 

2k

 

 

 

pA

 

 

pB

 

 

 

 

 

 

 

 

 

 

 

 

pB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

vB , или vB =

 

 

 

 

 

 

 

 

 

.

 

(3)

 

 

 

 

k 1

ρA

 

 

2

k

 

 

ρA

 

 

 

 

 

 

 

 

 

 

 

 

ρB

 

 

 

 

1

 

 

ρB

 

 

Вынося отношение

 

pB

за скобки, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρB

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2k pB pAρB

 

 

 

pB NρB

 

 

 

 

 

 

 

 

 

 

vB

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

=

 

 

 

 

 

 

 

 

 

1 .

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

ρA

 

 

 

 

 

 

 

 

 

 

 

 

1 k ρB pBρA

 

 

 

ρB

 

 

 

 

 

 

 

ρB

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

 

pB

 

 

 

= N

 

 

 

K

, то окончательно vB =

 

 

 

 

2k pB

 

k

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

1 .

 

ρA

pA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1 ρB

 

 

 

4.8. Адиабатическое течение идеального (без трения) совершенного (подчиняющегося уравнению Клапейрона) газа представляет собой баро-

тропное течение и описывается уравнением p = cρk , где c и k – постоянные,

причем k = c p – отношение теплоемкостей при постоянном давлении и по- cv

стоянном объеме. Найти зависимость между температурой и плотностью, а также между температурой и давлением для данного процесса.

Решение.

 

 

 

 

 

 

 

 

Подставляя p = cρk в уравнение состояния p RT ,

находим соотно-

шение между температурой и плотностью:

 

 

 

 

 

ρk 1

=

R

= const.

 

 

 

 

 

T

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

p

 

 

 

 

 

 

k

.

 

Кроме того, из уравнения баротропного течения ρ =

 

 

Подставляя её в

 

 

 

 

c

 

 

уравнение состояния, получим зависимость между температурой и давлением

 

k 1

 

 

 

 

 

p

k

 

=

R

 

= const.

T

 

1

 

 

 

 

 

c k

 

 

 

 

 

 

 

 

4.9. При движении твердого тела с неподвижной точкой поле скоростей имеет вид vi = εijk ωj xk . Доказать, что для такого движения формула теоремы

момента количества движения для континуума сводится к известному уравнению моментов динамики твердого тела вокруг неподвижной точки.

50

Решение.

Формула теоремы моментов для континуума имеет вид

S εijk x j pnk dS + τ εijk x j ρfk dτ = dtd τ εijk x j ρvk dτ.

В этом выражении слева стоит сумма моментов всех поверхностных и объемных сил относительно начала координат. При vk = εipq ωp xq это уравнение

преобразуется следующим образом:

M i = dtd τ εijk x jρεkpqωp xq dτ = dtd τ ρωp ( δipδjp − δiqδjp )x j xq dτ =

 

d

 

 

 

d

(ωp Iip ),

=

ωp ρ( δip xq xq xp xi )dτ

 

 

 

dt

τ

dt

 

где Iip = ρ(δip xq xq xp xi )dτ – тензор моментов инерции.

τ

Задачи для самостоятельного решения

4.10. Найти вид уравнения неразрывности для безвихревого движения, для которого вектор вихря скорости тождественно равен нулю rotV 0.

Ответ: ddtρ 2Y = 0.

4.11. Показать, что поле скоростей V = rA2 (x1e€1 + x2e€2 + x3e€3 ), где

r 2 = x12 + x22 + x32 и A – производная константа, удовлетворяет уравнению неразрывности несжимаемой жидкости.

4.12. Пусть функция Pij ...(r ,t) представляет любую скалярную, векторную или тензорную величину, отнесенную к единице массы сплошной среды,

так что Pij ...(

 

,t)= ρPij* ...(

 

,t).

 

d

ρPij* ...(

 

,t)δτ = ρ

dP* ...(

r

,t)

r

r

Показать, что

r

ij

 

 

 

δτ.

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

τ

 

 

 

 

4.13. Если

pij

= −pδij ,

то мощность напряжений можно представить

выражением D P

=

p

 

dρ

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij ij

 

ρ dt