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

Механика сплошной среды

..pdf
Скачиваний:
33
Добавлен:
15.11.2022
Размер:
13.71 Mб
Скачать

h = Y {l2' ~ VliVi,)’

(5l35)

I3=\Vi/\-

Объему Vo в (5.28) теперь соот­ ветствует объем частицы в момент t (р0 — плотность в этот момент), объему V — объем в момент t+ dt (р—плотность). Следовательно, dQ= = —dp/p, и потому из (5.28) имеем

J6_= ___1_ _Ф_=

dt

р dt

(5.36)

= div v = - ^ - + —L+ —

дхх дхг дх3

причем dp/dt — субстанциональная производная:

dt

dt

dxi

Из (5.36),находим

 

 

 

0 = In

(5.36')

 

 

Р

Следовательно, мы получили условие сохранения массы части­ цы в эйлеровом пространстве, выражающееся одним из уравнений

dt

dxi

(Pyi)= 0,

dt -f-divv=0.

(5.37)

Условие объемной несжимаемости среды имеет

вид dpjdt = Qy

или

 

divv = 0.

 

(5.38)

 

 

 

Условия совместности компонент

тензора

скоростей

деформаций

Vij получаются из

(5.33)

заменой e^j на v

и х на х.

потенциаль­

Движение среды называется

безвихревым, или

ным,, если J2=rotv = 0. При этом из условий

 

 

^ L - ^ L = 0

(г, / = 1, 2, 3)

 

 

dxj

dxi

 

 

 

следует существование потенциала скоростей ф(лт, t)

Vi = -^~, v=grad<p. dxt

Если к тому же среда несжимаема, то из (5.38) следует

divgrad ф=Дср = 0,

где

 

 

 

л

а2

, а2

, в2

Д = ------- 1---------1-------

 

дх\

дх\

дх\

и называется оператором Лапласа.

Потенциал скоростей ср — гар­

моническая функция. Этим свойством обладает движение идеаль­ ной несжимаемой жидкости в потенциальном поле сил, если в ка­ кой-нибудь момент rotv = 0.

Движение любой сплошной среды, рассматриваемое в эйлеро­ вом пространстве, обладает некоторыми свойствами, вытекающи­ ми из определений линий токов, вихрей и закона сохранения мас­ сы. Как уже отмечалось в § 3, линией тока в момент t называется траектория вектора скорости v(x, t), проходящая через какую-ни­ будь точку хо, т. е. линия, определяемая дифференциальным урав­

нением

 

 

 

 

 

 

d * = v (* , 0

(Vi

5Щ

.V3

3

9

)

где Я— параметр. Решение этого уравнения

 

 

 

x = f ( x 0i X, t);

X i = f i ( x о Ху t)

(i= 1, 2, 3)

 

(5.40)

является параметрическим уравнением линии тока. Если в момент /=const рассмотреть какую-нибудь линию ж0=ф(р), то из (5.40) получим х=\ (|х, Х)у т. е. уравнение поверхности. Если линия х0= =ф(р) замкнута, то поверхность *=f(p, X) называется трубкой тока.

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

^ pvvd2=^ piyfS,

&&

где v — нормаль, w = vv— нормальная составляющая скорости. V,

— элемент площади поверхности.

Вектор скорости v лежит на поверхности трубки тока, и пото­ му поток массы среды внутрь трубки равен нулю. Рассмотрим в фиксированный момент t внутри области движения среды произ­ вольный объем Vy ограниченный замкнутой поверхностью 2. Про­ интегрируем по V умноженное на dV= dxxdx2dxz уравнение (5.37); получим

d V = - [ ^ - d V .

J dt

Но по формуле Грина — Остроградского

Jd iv A d F = j - | ^ dK =jA vd2.

(5.41)

Следовательно,

(5.42)

т. е. поток массы через замкнутую неподвижную поверхность Е равен секундному ее изменению в объеме V с обратным знаком.

Из определения вектора вихря Я (5.34) следует

div Я=

dQi

:0.

(5.43)

дх{

 

 

 

 

Внося Q вместо А в (5.41), получим теорему Стокса

 

J Q vdS=0,

 

(5.44)

т. е. поток вихря через

любую

замкнутую

поверхность равен

нулю.

поверхность, ограниченная замкнутым

Если ЯГ— незамкнутая

контуром L, то поток вихря через ЯГ равен

интегралу по 6Г от

QvdE. Но для любого вектора В по формуле Стокса

J rotBvdE = j BdS,

(5.45)

if

 

L

 

 

где dS обозначает вектор-элемент длины дуги контура L; выра­ жение, стоящее в правой части (5.45), называется циркуляцией вектора В по L. Полагая B = V получим

J Q vd2=^ vdS =

Г,

(5.46)

т. е. еще одну теорему Стокса: поток

вихря через

незамкнутую-

поверхность ZT равен циркуляции Г вектора скорости по ограни­ чивающему ЯГ контуру L.

Отметим еще важное для дальнейшего свойство интегралов от различных функций ф(ж, t) по объему Vg фиксированной массы среды, т. е. ограниченному замкнутой поверхностью 2*, состояще­ му из неизменных физических частиц. Найдем полную производ­ ную по времени от интеграла

1= |ф(лг, 0 р (*. 0 dV

(5.47).

где

л

а2 ,

а2 ,

д2

Д = ------- 1---------1-------

 

дх\

дх\

дх\

и называется оператором Лапласа. Потенциал скоростей ср — гар­ моническая функция. Этим свойством обладает движение идеаль­ ной несжимаемой жидкости в потенциальном поле сил, если в ка­ кой-нибудь момент rotv = 0.

Движение любой сплошной среды, рассматриваемое в эйлеро­ вом пространстве, обладает некоторыми свойствами, вытекающи­ ми из определений линий токов, вихрей и закона сохранения мас­ сы. Как уже отмечалось в § 3, линией тока в момент t называется траектория вектора скорости v(x, t), проходящая через какую-ни­ будь точку *о, т. е. линия, определяемая дифференциальным урав­

нением

 

 

 

dAT=v(j»r, OdV —

Щ

v3

(5.39)

t>i

 

где К— параметр. Решение этого уравнения

 

х = Ц х 0, к t)\ x£= f i ( x 0i к

t) (i= 1, 2, 3)

(5.40)

является параметрическим уравнением линии тока. Если в момент /=const рассмотреть какую-нибудь линию ж0=ф(р), то из (5.40) получим *=f(|x, К) у т. е. уравнение поверхности. Если линия х0= =<р (р) замкнута, то поверхность x=f(p, X) называется трубкой тока.

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

^ pvvd2=^ рuvd2,

&&

где v — нормаль, vv = t\, — нормальная состазляющая скорости. V, d2 — элемент площади поверхности.

Вектор скорости v лежит на поверхности трубки тока, и пото­ му поток массы среды внутрь трубки равен нулю. Рассмотрим в фиксированный момент t внутри области движения среды произ­ вольный объем V ограниченный замкнутой поверхностью 2. Про­ интегрируем по V умноженное на dV=dx\dx2dx^ уравнение (5.37); получим

Jp -dV.

J dt

Но по формуле Грина — Остроградского

 

j" div AdV s= J

1

dV = ^AvdZ.

(5.41)

V

v

z

 

Следовательно,

 

 

 

 

 

j pvvd2= —J -|£- dK,

(5.42)

2V

т.e. поток массы через замкнутую неподвижную поверхность 2 равен секундному ее изменению в объеме V с обратным знаком.

Из определения вектора вихря Я (5.34) следует

 

d i v Q = - ^

= 0.

(5.43)

 

 

dxt

 

 

Внося Я

вместо А в (5.41), получим теорему Стокса

 

 

Q vd2= 0,

(5.44)

 

 

2

 

 

т. е. поток вихря через

любую замкнутую

поверхность равен

нулю.

— незамкнутая

поверхность, ограниченная замкнутым

Если

контуром L, то поток вихря через

равен

интегралу по ST от

Qvdl,. Но для любого вектора В по формуле Стокса

 

f rotBvdS =

j BdS,

(5.45).

 

if

 

L

 

где dS обозначает вектор-элемент длины дуги контура L; выра­ жение,стоящее в правой части (5.45), называется циркуляцией вектора В по L. Полагая B = V получим

Q vd2=^ v d S ^ r,

(5.46)

*L

т.е. еще одну теорему Стокса: поток вихря через незамкнутую» поверхность вГ равен циркуляции Г вектора скорости по ограни­ чивающему ф0контуру L.

Отметим еще важное для дальнейшего свойство интегралов от

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

/ = J ф (ЛГ О Р (•*. О d.V

(5.47),

Так как область интегрирования

Vg меняется

со временем, то

на основании закона движения

частиц х=х(х,

t) преобразуем

(5.47) к лагранжевым координатам. Из условия сохранения мас­

сы физической частицы имеем

 

 

рdV=p(x,t) dxldx2dx3=

р0dV0— р0 (х) dxtdx2dx3;

(5.48)

на основании закона движения

 

 

ф(*. 0 =

x(x, t).

(5.49')

■Область VQ (для фиксированной массы) не изменяется со време­ нем. Тогда получим

I J P<pdV= J х (х, t) Ро (х) dV0

Vg Vo

и, следовательно,

dx(s. t) Ро (*) dV0. dt

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

 

дХ(х, 0 _

4ф(*, 0

/5 49)

 

dt

dt

К’ ’

получим

 

 

 

J ф (х, О Р (х, t) dV= J dq> (*’-Q PXX>t) dV

ve

'

vg

 

Если объем V и граница области 2 фиксированы и неизменны

в эйлеровом пространстве, то, обозначая их Vc, 2 С

получим

 

■ f

 

 

 

Ус

vc

 

§ 6. ТЕОРИЯ НАПРЯЖЕНИИ

Основная идея теории внутренних сил в МСС уже изложена в

§3: любой объем, мысленно выделенный внутри среды, находится

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

Выделим в некоторой точке М (х) среды в деформированном состоянии бесконечно малый тетраэдр МАВС (рис. 6.1), основные

ребра которого МА, МВ и МС направлены по векторам репера

Эг, Эз.

Обозначим удвоенные площади треугольника МВС через d2*,

треугольника MAC — через d22

треугольника МАВ — через d23 а

треугольника АВС — через d2v.

Единичные векторы нормалей

площадок d2“, очевидно,

коллинеарны векторам эа (а=1, 2, 3).

Величины этих площадок определяются через начальные

d2“= d x pdxv

 

 

соотношениями (4.27)

 

 

d2“= ] /§ £ ““ d 2“

(6.1)

Напомним, что по греческим индексам суммирование не проис­ ходит. Вектор-площадь треуголь­ ника АВС подсчитывается по фор­ муле (4.31):

d2v • х= У ]гс1 Ь э‘,

(6.2)

 

Рис. 6.1

где v — единичный вектор внешней нормали к площадке АВС:

 

x

= v ia ‘ = v ‘ 9 i ,

(6.3)

причем

Vj — ковариантные,

v* — контравариантные

компоненты v.

Из

формул (6.1) и (6.2)

находятся выражения d2“ через d2v

и v:

 

 

 

 

d2a= V ^ v ad2v

(6.4)

Как установлено в § 3, взаимодействие рассматриваемой нами частицы-тетраэдра с окружающей средой реализуется за счет векторов внутренних сил, действующих по граням, с точностью до малых порядка |dx|, равномерно по ним распределенных. На каждую из площадок d2a действует поверхностная сила плотно­ сти Р<а> а на площадку d2v— поверхностная сила плотности P<v). Эти векторы называются векторами истинных внутренних напря­ жений. С точностью до малых высшего порядка силы, действую­ щие по граням тетраэдра, равны

—P(a)d2a -f-0(|dx|3) ( a = l, 2, 3).

(6.5)

Мы здесь учли, что векторы внешней нормали к площадкам d2“ направлены в обратную сторону по отношению к векторам эа. Очевидно, элементарный объем тетраэдра dV'=d2v-A/3, где h — величина перпендикуляра, опущенного из точки М на d2v. Учиты-

вая еще равенства (6.4), (6.5), получим на основании принципа Даламбера

 

з

J+ О (| dx|3)= 0 , (6.6')

dSv (pF—pw) h + P(v)—

Р(аЧ V I™

 

a—\

 

где w — ускорение; слагаемое

порядка h учитывает силу инерции

и массовую силу F, действующие на частицу. Сокращая на d2v и

устремляя Л->-0 вместе с \dx\, получим условие

равновесия

бес­

конечно малого тетраэдра

 

 

 

 

__

__

__

з

__

 

p (v)= ] / g u р ° Ч + у

g™ P (2)V 2 + у

g 33 p (3)v3=

V

У g aa p <a)v a .

(6 .6)

a=l

Эта очень важная формула называется выражением вектора ис­ тинного напряжения на косой площадке с нормалью v через ос­ новные координатные векторы напряжений. Она доказывает, что

вектор внутреннего напряжения P(v) на площадке с нормалью v —

линейная функция v. Формула (6.6) справедлива как для внут­ ренних площадок d2v так и для площадок, расположенных на границе 2 области движения среды; из закона равенства дейст­

вия и противодействия заданная

на границе тела

сила

Рвнеше-

равна P<v> имеющей выражение

(6.6). Естественно,

что

Рвнешя.

может быть любой функцией координат и нормали v на 2.

 

Обозначим векторы, называемые

основными контравариант-

ными векторами напряжения:

S1=Vg™ Р(1)

 

 

 

 

 

sa= V g ™ v{a),

s 2=Vg™ р (2)

 

(6.7)

 

83=

]/^зз р(3)

 

 

Как видим, векторы Sa не зависят от v и отличаются от векторов

истинных напряжений, хоть

и коллинеарны

 

им, так как

gaa~-ф-\

при всевозможных значениях t.

 

 

 

 

Из (6.6) с помощью новых обозначений получим

 

P(v) =

S'v,= S‘<vt9,,

 

 

(6.8)

при этом векторы S’ и Р(а) мы представили в репере эj

 

 

sa

Sal3j

(6.9)

 

У ^™ '

 

V i

 

 

Нормальная составляющая

вектора

P(v> на площадке

dSv яв-

.ляется скалярной величиной и с учетом

(6.3)

имеет значение

Поскольку

v — вектор и ЛЛУ>— скаляр,

то из (6.10)

следует (по

обратному

признаку тензора), что S ij

представляют

контравари-

антные компоненты тензора напряжений S в лагранжевом репе­ ре 9i\ в таком представлении S называется тензором напряжений Коши Лагранжа. Касательная составляющая Pv равна

Tm = V (V (v))2—(N(v))2.

 

(6.11)

Так как в (6.7) множитель Y %,аа при Р<а>различен для

разных а,

то выражение Р<а) в репере 9j

 

 

 

р(“)= р(“/>Э/

 

 

(6.12)

приводит к объекту P(lj), не являющемуся тензором.

к

частице

Тензор S'i симметричен, поскольку

приложенные

внешние (массовые) силы дают момент относительно

центра ча­

стицы, являющийся величиной, малой

более высокого

порядка

сравнительно с бесконечно малым порядка dV\

Обычно рассмат­

риваемые в механике силы,

как указывалось,

обладают этим

свойством.

компонент

следует из уравне­

Доказательство симметрии

ния моментов для параллелепипеда, если учесть, что момент век­ торов сил Р(а>d2a, действующих по всем шести граням, должен равняться нулю с точностью, включающей малые порядка dV. При этом напряжения Р(а) равномерно распределены по соответ­

ствующим граням

и, значит,

векторы Р(a)d2a

приложены

в цен­

трах граней параллелограммов. Например, вектор

Р(3)d23

прило­

жен в точке с координатой

(относительно М)

 

г3=

1/2 (э^х1 +

+ 32dx2) на нижней площадке

d23 и равный

ему (с точностью до

малых высшего порядка) — в точке с координатой

г '=

1/2 (э^хг +

+ 3 2dx2) + a3dx3—

на верхней d23. С точностью до

малых выс­

шего порядка момент сил, действующих по нижней и верхней гра­ ням dS3 равен

m 3 = P ( 3 ) d 2 3 X

г з — Р (3) d l a х r 3 = P

(3) X ( г 3 — г 3) d 2 3 = P (3)

х

9 3dx3d S 3.

 

 

 

 

(6.13)

Но

_

_

 

 

 

dS3= y g V g33 dx1dx2,

 

 

и потому, учитывая (6.7) и (6.9), получаем

 

 

m3= V g'dx1dx2(ix3 S3 х 93= dVS3’9j X Э3.

 

(6.14)

Моменты nil и ш2 по граням dS1 и dE2 аналогичны,

и потому

условие т = т ( + m2 + m3 = 0 принимает вид

 

 

m =dV {Slj9 / х

+ S2'3j X 9 2 + S3i3 j x 93}— dVS'iaf X

9t= 0. (6.15)

Учитывая, что

X Э у= —Э/ X Эи

из (6.15) имеем

 

 

 

( S 1 2 —

S 2 1 J Э3 + (S23— S32) э 1 +

(S31—513) э2= 0,

 

(6.16)

откуда и следуют условия симметрии S^ = S^.

взаимности

Из свойства симметрии тензора

5

следует закон

напряжений на

основных площадках

d2a; умножая

Sa

(6.7)

на

эр и учитывая

(6.7), получим выражение компонент SaP

через

ос­

новные векторы истинных напряжений и единичный базис кр:

S“p= l / g “agPPP(a)kp, кр= —

 

1*р1

 

Из симметрии SaP = SPa находим закон взаимности

 

p(a)kp= P (P)k°:

(6.17)

проекция Р(а) на направление

равна проекции

Р(Р) на направ-

ление эа.

 

 

С помощью метрического тензора ^ и аффинора А тензор на­ пряжений S можно представить в репере э1 другими компонента­ ми, например ковариантными и смешанными

S u = S ”» g tmg nl= S ? 4 g in = S \ ? g mh

(6.18)

причем вектор напряжений P<v) согласно (6.8) при замене э,= =9mgmj и Vi = vngin получит выражения

 

 

Piv)= Snmvn9m= S ‘ivldi.

 

(6.19)

Тензор S можно представить также и в других реперах. Напряже­

ние можно представить и другими матрицами,

не являющимися

собственно тензорами. Приведем две из них.

 

площадкам.

Напряжения, отнесенные к

недеформированным

Пусть

начальная

элементарная

вектор-площадка

О

0 0

dS= d2n с еди-

ничным

вектором

о

о

 

 

нормали п=/гг е2 в результате деформации

превращается в вектор-площадку d2 = d2v

с единичным

вектором

нормали v = Vi3\ Тогда, очевидно,

0

 

 

0

0

d2г в

репере ег равна dXnh

и потому, полагая в формуле

(6.2)

d2v = d2

и умножая

обе части

(6.2) на 3j (с учетом 3i9j = 6ji),

получим

 

 

 

 

г - 0

о

dZ

г-

0

(6.20)

dI,vi= Y g n id'E,

\ j = Y g n h

<i2 Вектор t(v), определенный равенством