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

mehanika_spl.sr.tekst_lekciy_popov_savich

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

p

p

 

p

 

 

 

11

 

21

 

31

 

(4.10)

p12

p22

p32

.

p

p

23

p

33

 

 

13

 

 

 

 

Таблица (4.10) не удовлетворяет принятому порядку индексов: первый - номер строки, второй - столбец. Поэтому за тензор напряжений примем тензор Р с таблицей

 

p11

p12

p13

 

 

 

P

p

21

p

p

23

 

,

(4.11)

 

 

 

22

 

 

 

 

 

p

p

p

 

 

 

 

 

 

31

32

33

 

 

 

а тензору с таблицей (4.10) сопоставим тензор Pc , сопряженный с тензором Р. (отличающийся порядком индексов pik pki*C ).

В дальнейшем покажем, что различие между тензорами. P и P* отсутствует, поэтому операции умножения вектора на тензор справа или слева равны и формулы (4.9) можно переписать в тензорной форме так:

p

n nˆP Pcnˆ.

(4.12)

Эта форма равенств Коши не зависит от выбора системы координат (1827 г.)

Тензор второго ранга - это величина, определяемая в любой системе координат девятью числами (или функциями) pik , которые при изменении системы координат преобразуются в

pik' по закону

P'

'

'

m

P .

ik

il

k

lm

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

p11 1, p22 2, p33 3, p12 12, p23 23.

и т.д.

3. Теорема о взаимности касательных напряжений

Для равновесия тетраэдра MM1M2M3 необходимо равенство нулю и главного момента

приложенных сил. При составлении условий равновесия сил приложенных к тетраэдру с использованием главного вектора положение точек приложения поверхностных сил не имело значения, а для определения моментов поверхностных сил необходимо знать координаты точек приложения этих сил. Поскольку грани тетраэдра бесконечно малы, можно принять, что поверхностные силы по ним распределены равномерно и, следовательно, точки приложения равнодействующих сил находятся на пересечении медиан граней тетраэдра в точках N,N1,N2,N3 . Обозначим через r,r(1) ,r(2) ,r(3) радиусы - вектора этих точек. Чтобы не загро-

мождать чертеж, на рис. (4.2) показано положение точек N и N3 и их радиусы вектора. Со-

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

r pn Sn r(1) p1 S1 r(2) p2 S2 r(3) p3 S3.

Заменяя в правой части этого равенства pn его разложением (4.7), получим после сокраще-

ния на Sn

(

r

 

r

(1) )

p

n (

r

 

r

(2) )

p

2

n

2

(

r

 

r

(3) )

p

n

3

0.

(4.13)

1

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

Направленный отрезок

 

равен разности

r

 

и

r

(3)

 

 

 

 

 

 

N3N

 

(3) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N3N

MN

MN3

 

r

r

 

 

 

 

31

а так как точка N3 является проекцией точки N на плоскость MX1X3 , то направленный отре-

зок N3N параллелен вектору eˆ3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично, векторы

r

 

r

(1)

и

r

 

r

(2)

параллельны векторам eˆ

и

eˆ

2

, так что можно

положить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

(1)eˆ ,

 

 

 

 

 

(2)eˆ

 

 

 

 

 

 

(3)eˆ .

 

 

 

 

r

 

r

(1)

r

 

r

(2)

2

,

r

 

r

(3)

 

 

 

(4.14)

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

Подставляя значения этих разностей в (4.13), получим

 

 

 

 

 

 

 

 

(1)n (eˆ

 

p

) (2)n (eˆ

 

p

) (3)n (eˆ

 

p

) 0.

 

 

 

(4.15)

1

 

 

1

1

2

2

2

 

 

3

3

3

 

 

 

 

Остается определить скаляры (1) , (2), (3) . Для этого продлим отрезки MN1,MN2,MN3 до

их пересечения с плоскостью M1M2M3 в точках N1',N2' ,N3'

(на рис. 4.2 показана только точ-

ка N3' ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

MN|

 

3

MN

 

 

3

r(1),MN|

 

3

MN

 

 

3

r(2),MN|

 

3

MN

 

 

3

r(3) .

 

 

 

 

 

 

 

 

 

1

2

 

1

2

2

2

 

2

2

3

2

 

3

2

 

'' '

Поскольку концы векторов MN1,MN2,MN3 так же как и MN , лежат в наклонной плоскости M1M2M3 , проекции их на нормаль n к плоскости M1M2M3 будут равны одной и той же величине кратчайшему расстоянию h точки М до плоскости M1M2M3 , так что

3

r

(1)

n

 

3

r(2)

n

 

3

r(3)

n

h.

(4.16)

2

2

2

 

 

 

 

 

 

 

 

 

 

 

Умножая обе части каждого из равенств (4.14) скалярно на n , получим с использованием

(4.16)

 

 

 

 

r

 

r

(1) (1) eˆ ,(

r

 

 

r

(1))

n

 

 

(1)eˆ

 

n

(

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

r

 

r

(k) )

n

r

 

n

 

r

(k)

n

h

2

h

1

h,(k 1,2,3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

 

(1)eˆ

n

(1)eˆ

(eˆ

n

eˆ

n

2

eˆ

n

) (1)n

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

1

 

 

2

 

 

 

 

 

3

3

 

 

1

Отсюда следует, что (1)n

(2)n

 

 

(3)n

 

1

h.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя полученные значения (k)

 

в равенство (4.15) и сокращая обе части на h, при-

дем к следующему векторному соотношению:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eˆ1

p1 eˆ2

p

2

eˆ3

 

p

3

0.

(4.17)

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

p12 p21, p23

p32 , p31

p13.

(4.18)

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

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

32

X3

M 3

 

eˆ3

 

 

N

 

 

 

 

eˆ1M

 

 

 

eˆ2

X 2

 

 

 

 

M2

N3

N3|

 

 

 

 

 

 

 

M1

X1

33

Лекция 5

Уравнения статики сплошной среды “в напряжениях”

Принцип затвердевания, примененный к деформируемому телу, является необходимым, но недостаточным условием равновесия.

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

Выделим в сплошной среде произвольный объем (рис.5.1) ограниченный

поверхностью S . К элементарному объему приложена сила F , а к элементарной площади S , ориентированной согласно единичному вектору нормали nˆ к поверхности Sx

сила n S , где n вектор напряжения в точке N поверхности S .

f

nˆ M

N

S

pn S

S

Рис.5.1

Запишем необходимое, согласно принципу затвердевания, условие равновесия равенство нулю главного вектора приложенной к объему и k ограничивающей его поверхности S совокупности сил. Будем иметь

f

p

n S 0,

(5.1)

 

S

 

где первый интеграл, взятый по объему , представляет собой главный вектор внешних,

объемных сил f , приложенных к элементам объема , а второй, взятый по поверхности

S , главный вектор внешних поверхностных сил pn S , приложенных к элементам поверхно-

сти S .

 

 

 

Применяя равенство Коши в форме

p

n

 

nP, приведем уравнение (5.1) к виду

 

 

 

nˆP S 0.

 

f

(5.2)

 

 

S

 

Используя формулу Гаусса – Остроградского, перейдем от поверхностного интеграла к объемному его эквиваленту. Будем иметь

nˆP S DivP ,

S

и тогда уравнение равновесия (5.2) приведется к виду

 

f

DivP 0.

(5.3)

 

 

Пользуясь произволом в выборе объема , применим равенство (5.3) к элементарному объему . Тогда интеграл, стоящий в левой части, сведется к одному слагаемому, получим

( f DivP) 0.

Замечая, что представляет собой произвольную бесконечно малую величину, но не равную нулю, разделим обе части предыдущего равенства на и получим искомое уравнение статики «в напряжениях»:

34

DivP f

0,

(5.4)

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

x3 , приведется к системе трех дифферен-

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

 

 

p11

 

 

 

p12

 

 

 

 

p13

 

 

f 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

x2

 

 

 

 

 

x3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p21

 

 

p22

 

 

 

p23

f2 0.

(5.5)

 

 

x1

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3

 

 

 

p31

 

p32

 

p33

f3 0.

 

 

 

 

 

 

 

 

x1

 

 

 

 

x2

 

 

 

 

 

x3

 

 

Согласно свойству симметрии тензора напряжений,

касательные напряжения,

 

отличающиеся порядком индексов, равны между собой:

 

 

p12 p21,

 

p23 p32 ,

p31 p13 .

 

И система трех уравнений равновесия заключает в себе не девять, а только шесть неизвестных компонент тензора напряжений, что делает ее недостаточной для определения этих шести неизвестных.

Эту статическую неопределенность можно избежать, вводя дополнительные допущения о физических, свойствах сплошной среды, например о ее упругости, подчиняющейся закону Гука о пропорциональности тензора напряжений тензору деформаций.

Наиболее просто решается вопрос о равновесии идеально текучей среды, в которой касательные напряжения отсутствуют

 

p12 p21 p23 p32 p31 p13

0,

(5.6)

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

 

 

 

 

 

 

p

n pnnˆ,

p

1 p11eˆ1,

p

2 p22eˆ2,

p

3

p33eˆ3.

(5.7)

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

 

 

 

 

 

 

 

 

 

p11 p22

p33 pn ,

 

 

 

(5.8)

т.е. для идеальной жидкости три нормальных напряжения, приложенные к трем взаимно перпендикулярным площадкам, как угодно ориентированным в пространстве, равны между собой. Этот закон изотропии нормальных напряжений в точках сплошной среды, находящейся в равновесии, был открыт в середине 17 века Паскалем.

Общее значение нормальных напряжений в данной точке среды, взятое со знаком минус, назовем давлением в этой точке и обозначим р, так что

p11 p22 p33 p. (5.9)

Чтобы подчеркнуть, что принятое определение соответствует лишь случаю равновесия среды, это давление называют гидростатическим давлением.

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

p

n pnˆ ,

(5.10)

приложенное в точках поверхности мысленно выделенного объема, направлено в сторону, противоположную орту внешней нормали к этой поверхности, т.е. внутрь объема.

Давление p представляет собой физический скаляр, такой же как и плотность, температура, концентрация и др. скалярные характеристики.

При равновесии среды тензор напряжений в ней, согласно (5.6) и (5.9) имеет таблицу

p 0

0

1

0

0

 

 

p

 

 

 

1

 

 

0

0 p 0

0 ,

 

0

0

 

 

0

0

1

 

 

p

 

 

т.е. обладает сферической симметрией, что соответствует свойству изотропии нормальных напряжений в покоящейся среде и можно положить (Е - тензорная единица)

35

 

 

 

 

 

P pE .

 

 

 

 

(5.11)

В этом случае, легко видеть, что DivP gradp,

и уравнение равновесия (5.4) в вектор-

ной форме примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.12)|

 

 

 

 

 

gradp ,

 

 

 

 

 

 

f

 

 

 

 

а в проекциях на оси координат,

 

p

 

 

 

 

p

 

 

 

p

 

 

f

 

, f

 

 

, f

 

 

,

(5.12)

 

 

 

 

 

1

 

x

 

2

 

x

2

 

3

 

x

 

 

1

 

 

 

 

 

 

 

3

 

 

эти уравнения носят наименование уравнений Эйлера гидростатики.

Уравнения (5.12 ) легко интегрируются для случая объемной силы тяжести равной f g . Жидкость подверженную силе тяжести называют «тяжелой» жидкостью. В этом случае уравнение (5.12) примет вид

 

g

gradp.

 

(5.13)

Направляя ось Ox3 вертикально вниз, будем иметьg1 g2

0,

g3 g,и уравнения (5.12)

сведутся к одному такому

 

 

g

p

.

 

(5.14)

 

 

 

 

 

x3

 

 

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

p gx3 const

 

(5.15)

Погрузим в тяжелую жидкость с плотностью твердое тело объема

и с поверхностью

S . Главный вектор

 

пов сил давления жидкости на поверхность тела, согласно равенству Га-

R

усса - Остроградского, будет равен

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïîâ

np S gradp ,

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

Но по (5.13)

 

 

 

gm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

,

 

 

 

 

 

 

gradp

g

g

 

G

 

 

 

 

 

 

gm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

где

 

- вектор силы тяжести жидкости в объеме , так что

 

G

 

 

 

 

 

 

 

 

 

ïîâ

 

G

 

g

 

.

(5.16)

 

 

 

 

 

 

 

R

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

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

всвоем весе столько, сколько весит жидкость в объеме погруженного в нее тела.

Главный вектор Rпов называют архимедовой силой или гидростатической подъемной силой в знак того, что она стремится вытолкнуть тело из окружающей его жидкости.

Легко вычислить также момент сил давления тяжелой жидкости на поверхность погруженного в нее тела. Имеем

L

r

 

np S

n

 

rp S

rp rot(pr

) .

(5.17)

S

 

 

S

 

 

 

 

 

Замечая, что по известной формуле векторного исчисления и очевидному тождеству rotr 0

rot(pr) protr gradp r r gradp r g ,

перепишем (5.17) в форме

 

L r gradp r g .

(5 18)

Вводя выражение вектор - радиуса rc центра тяжести вытесненного телом объема жидкости

rc

1

 

r

 

1

 

r

g ,

 

 

 

 

 

G

 

 

 

 

36

убедимся, что момент L можно представить как

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

r

g

 

 

r

 

 

,

 

 

 

 

 

 

L

 

 

 

G

G

 

 

 

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пов по (5.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или, т.к. G

R

 

 

 

 

 

 

 

 

пов .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.18)

 

 

 

 

 

 

 

 

 

L

r

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

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

Из уравнения (5.12) можно исключить плотность и давление. Для этого возьмем сначала от обеих его частей операцию вихря rot; тогда р исключится, так как

rotgradp 0,

будем иметь rot( f ) 0. Отсюда, раскрывая скобки по известному правилу, получим

rot( f ) x ( c f ) x ( fc) ( f ) grad f rotf grad f 0.

Умножая обе части этого равенства скалярно на f , и, замечая, что второе слагаемое, как

векторное произведение, перпендикулярно к своему сомножителю f будет равно нулю, найдем следующее общее ограничение, накладываемое на класс сил, под действием которых возможно равновесие жидкостей и газов:

f rotf 0.

(5.19)

К числу объемных сил, удовлетворяющих условию (5.19), относятся силы имеющие потенциал П, так как для них

 

 

f

gradП, rotf 0.

(5.20)

В этом случае

 

 

 

rotf

 

grad

 

 

0 grad ( gradП) 0

или

 

f

 

grad gradП 0,

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

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

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

Если в движущемся или покоящемся газе плотность является функцией только давления, то движение или равновесие газа называют баротропным Из предыдущего следует, что баротропное равновесие среды возможно при наличии только потенциальных сил, так как при условии (p)изобары и изостеры совпадут, а, следовательно, силовое поле должно быть потенциальным. В этом случае уравнения Эйлера (5,12) представляют собой так называемую систему в полных дифференциалах и имеют общее решение вида

p dp

(p)

p0

или при наличии потенциала П (м)

p dp p0 (p)

M

f dr,

M0

П(M0) П(М),

(5.21)

(5.22)

где p0 p(M0 ). Чтобы задача была определенной, надо использовать дополнительные до-

пущения о форме и зависимости плотности от давления.

37

Q V ,
V
M

Лекция 6

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

Напомним формулу этой теоремы в дифференциальном виде

dQ

 

 

e,где

 

 

e F

e

,

R

R

 

dt

 

 

k

 

 

 

 

 

 

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

Выделим в движущейся сплошной среде произвольный жидкий объем , ограниченный поверхностью S и применим к нему теорему (рис. 6.1).

Напомним, что в сплошной среде вместо обычных сил вводятся плотности их распределения, соответственно в объемах и на поверхностях: f - для объемных и pn - для поверхно-

стных сил. Если элемент объема, содержащий точку М со скоростью V , то элементарное количество движения этого элемента равно

Q V ,

где - плотность распределения массы в объеме .

Главный вектор количества движения сплошной среды Q , равный векторной сумме элементарных количеств движения будет определяться вычисленным по объему интегралом

(6.1)

F

n

N

S

Pn

 

 

 

Рис.6.1

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

запишется так:

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

f

p

n S.

(6.2)

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

Таким образом, применительно к механике сплошных сред, формулировка претерпевает небольшие изменения и становится следующей: индивидуальная производная по времени от главного вектора количества движения «жидкого» объема равна главному вектору внешних объемных и поверхностных сил, приложенных к частицам расположенным в объеме и на ограничивающей его поверхности.

Производная по времени, стоящая слева, следит за всеми изменениями со временем – локальными и конвективными, главного вектора количества движения среды в движущемся вместе со средой объеме . Эту производную можно вычислить по общим правилам дифференцирования интеграла

 

d

 

 

 

 

dV

 

 

 

d

( ).

(6.3)

 

V

 

V

 

dt

dt

 

 

 

 

 

 

 

dt

 

Но в силу общего закона сохранения массы

движущегося элементарного объема

будет

38

 

d

( ) 0.

(6.4)

 

 

 

 

 

dt

 

и формула теоремы в окончательном виде запишется так:

 

 

 

 

 

 

 

 

 

 

dV

F

p

S.

(6.5)

 

dt

 

 

n

 

 

S

 

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

Воспользовавшись теоремой об изменении главного вектора количества движения, выведем общее уравнение динамики сплошной среды - так называемое «уравнение в напряжениях». Для этого преобразуем поверхностный интеграл правой части формулы (6.5) в объемный, произведя замену pn n1 p1 n2 p2 n3 p3 и используя формулу Гаусса-Остроградского,

будем иметь эквивалентные друг другу результаты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p1

 

p2

 

p3

 

 

 

 

p

 

 

(n p n

 

p

 

n p

 

 

 

 

 

 

 

DivP .

(6.6)

 

 

 

 

 

)

 

 

 

 

 

 

 

n

 

 

1 1

2

 

2

3 3

 

 

x

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

3

 

 

 

 

Подставим в (6.5) приведенные к объемным интегралам члены и объединим их, перенеся в левую часть уравнения

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

F DivP 0,

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

откуда в силу произвольности объема придем к уравнению динамики сплошной среды «в напряжениях»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

DivP,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или учитывая, что индивидуальная производная равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(V

 

)V

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

окончательно получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(V )V f DivP,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Перепишем его в аналитической форме,

взяв проекции левой и правой частей

равенства

(6.8) на оси прямоугольной декартовой системы координат, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

ki

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vk

 

 

 

 

 

 

fi

 

 

 

 

 

,(i 1,2,3),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в развернутом виде,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V1

 

 

 

 

V1

 

 

 

 

 

 

 

 

 

 

V1

 

 

 

 

 

 

 

 

 

V1

 

 

 

 

 

 

 

 

p11

 

p21

 

p31

 

 

 

 

 

 

V1

 

 

V2

 

V3

 

 

 

 

f1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

t

 

 

 

 

x

x

 

 

x

 

 

x

 

 

x

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

3

 

 

 

 

V

 

 

 

 

V

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

p

 

 

 

p

22

 

 

 

p

23

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V1

 

 

 

V2

 

 

 

V3

 

 

 

 

 

 

f2

 

 

 

 

 

 

 

 

 

 

 

 

,

(6.9)

 

 

t

 

 

 

x

 

 

x

 

 

x

 

 

 

 

x

 

 

x

 

x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

V

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

p

 

 

 

p

23

 

 

 

p

 

 

 

 

 

 

3

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

3

3

 

 

 

 

 

 

 

 

13

 

 

 

 

 

 

33

 

 

 

 

 

 

 

 

V1

 

 

 

V2

 

 

V3

 

 

 

 

 

f3

 

 

 

 

 

 

 

 

 

 

 

.

 

 

t

 

 

 

 

x

x

 

x

 

 

 

 

 

x

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

К этой системе уравнений присоединяется уравнение неразрывности, имеющее в указанных координатах вид

 

 

 

( V )

 

 

( V

)

 

 

( V ) 0.

(6.10)

 

x

x

 

x

 

t

1

2

2

 

3

 

 

 

1

 

 

 

 

 

3

 

 

Уравнения (6.9) и (6.10) не образуют замкнутую систему уравнений для определения тринадцати неизвестных: ,V1,V2,V3, p11, p21, p31, p21, p22, p23, p13, p23, p33. Для замыкания этой

39

системы уравнений необходимо применять те или иные допущения о математических моделях среды: идеально упругая среда, подчиняющаяся линейному закону Гука, идеальная жидкость, лишенная внутреннего трения, вязкая жидкость, движение которой описывается законом Ньютона, и т.д. Так, например, если принять модель идеальной жидкости как жидкости, в которой нет внутреннего трения, то все касательные составляющие тензора напряжений будут равны нулю

p12 p21 p23 p32 p31 p13 0,

а тогда, как это будет следовать из равенств Коши, все нормальные составляющие тензора напряжений будут равны между собой и можно положить

p11 p22 p33 p,

где р - скалярная величина, выражающая гидродинамическое давление. Согласно принятой модели идеальной жидкости, совокупность уравнений (6.9) и (6.10) перейдет в следующую замкнутую систему уравнений:

 

 

V1

V

V1

V

 

 

 

 

V1

V

V1

 

f

 

 

 

1

 

 

 

 

p

,

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x

 

2 x

2

3

x

 

 

 

 

1

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

1

 

 

 

 

V2

 

V

V2

 

V

 

 

 

V2

V

V2

f

 

 

1

 

 

 

p

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

1 x

2 x

2

3

x

3

 

 

 

 

2

 

 

 

x

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V3

V1

 

V3

V2

V3

V3

V3

f3

 

 

 

1 p

 

(6.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

t

 

 

x

x

 

 

x

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

3

 

 

 

 

 

 

 

3

 

 

 

 

 

 

( V )

 

 

 

 

 

( V )

 

 

( V ) 0,

 

 

 

 

 

 

 

 

 

x

 

 

 

t

 

x

1

 

x

 

 

2

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(p)

здесь зависимость плотности от давления задается последним уравнением этой системы, представляющим собой уравнение баротропности движения (если =const, то жидкость не-

сжимаема; если p k const, где к - показатель адиабаты, то происходит адиабатическое движение).

Возвращаясь к векторной форме уравнений, имеем

dV

 

 

 

 

 

 

 

 

 

 

 

1

gradp

 

(V

)V f

 

 

 

 

dt

 

 

 

 

 

 

d

divV

 

 

0.

(6.12)

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

(p)

3.Теорема об изменении момента количества движения и ее применение для доказательства равенства касательных напряжений

Напомним теорему: производная по времени от главного момента количества движения системы относительно какого-либо центра равна главному моменту внешних сил, приложенных к системе относительно того же центра

ddtK0 m0 (Fke ).

В динамике сплошной среды этот закон формулируется так: индивидуальная производная от главного момента количества движения конечного «жидкого» объема равна главному моменту внешних объемных и поверхностных сил, приложенных к этому объему и ограничивающей его поверхности. Обозначая через r текущий вектор-радиус точек объема и ограничивающей его поверхности S , выразим сформулированную теорему следующим равенством:

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]