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

книги / Механика сплошной среды. Законы сохранения

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

Определим скорость изменения нормали к материальной площадке. Дифференцируя (1.57)

 

 

 

n F

 

n ,

 

 

 

 

 

 

 

T

 

 

 

 

где (nˆ V 2 nˆ) 1/2 0 , получаем

 

 

 

 

 

 

0 F

T

 

T

n

F

T

 

 

n F

 

 

n ,

 

 

 

 

 

 

 

 

 

 

откуда с использованием (2.15)

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

n L

n n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

0. Далее

Дифференцируя n n 1, получаем n

/ n LT n

и затем

 

n L nn n L.

(2.64)

n

2.9. Запись уравнений совместности в дивергентной форме

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

F v ,

(2.65)

записанное в лагранжевых переменных, следует из (1.11) и (2.2). Из (2.15) с учетом (2.18) получаем то же уравнение

ˆ ˆ

(2.66)

F F v ,

записанное в эйлеровых переменных.

Общие полевые уравнения континуума принято записывать в

дивергентной форме, в лагранжевых

 

 

 

 

(2.67)

либо эйлеровых переменных:

 

 

 

ˆ

ˆ

ˆ

 

t

(v ˆ

) ˆ ,

(2.68)

где ( x,t), ˆ( xˆ,t), ( x,t), ˆ ( xˆ,t), ( x,t), ˆ ( xˆ,t) ― плотности, потоки и источники тензорной полевой величины в лагранжевых и эйлеровых переменных. Данный формат записи изменения распределения полевой величины материального континуума в пространстве и времени оказывается удобным с математической точки зрения, в частности, для вывода условий на

41

границах и численного интегрирования. Больше примеров его применения будет дано в следующих главах. В качестве полевых переменных в (2.66), (2.67) могут выступать произвольные тензорные (полярные или аксиальные) величины. Данные уравнения называются уравнениями баланса (в локальной форме).

Для приведения уравнения (2.65) к форме (2.67) его необходимо транспонировать и применить тождество (I v) v , в результате чего получим

FT (Iv) 0.

(2.69)

Таким образом, в терминах (2.67) FT , Iv, .

Уравнение (2.66) приводится к форме (2.68) с помощью тождества Пиолы

ˆ

1

F ) 0

,

(2.70)

(J

 

которое доказывается с помощью (1.55) и теоремы Гаусса – Остроградского:

 

 

 

 

 

 

ˆ

 

 

1

 

 

ˆ

 

 

ˆ

 

 

 

 

1

 

 

 

ˆ

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(J F )dV

 

(J F )dS

 

IdS

IdV

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

в силу

 

произвольности

 

 

объема

 

 

следует

 

(2.70).

 

Уравнение

(2.66)

приводится к виду J

1

 

T

J

1

 

T

 

ˆ

 

 

 

 

 

 

 

Далее,

 

учитывая (2.6)

и (2.62),

F

 

 

 

 

F

 

 

v = 0 .

 

пишем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(J

1 T

) / t

(J

1 T

)

 

 

 

 

ˆ

 

 

 

1 T

) J

1 T

J

1 ˆ

 

 

T

 

 

 

ˆ

1 T

).

 

F

 

F

 

v (J

 

 

F

 

 

F

 

 

vF

 

v (J

 

F

Записываем

 

 

 

 

тождества

 

 

 

 

 

ˆ

 

 

 

1

 

 

T

) J

 

1 ˆ

 

 

 

T

 

 

 

ˆ

 

1 T

)

 

и

 

 

 

 

 

 

 

 

(J

 

vF

 

 

 

 

vF

 

v (J

F

 

ˆ

 

1

 

ˆ

 

(J

1

 

 

 

 

 

1 T

 

 

ˆ

 

 

 

 

 

1 T

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(J

 

Fv)

F)v J

 

 

F

v J

 

F

 

v , где в последнем учтено (2.70),

и, собирая все вместе, приходим к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(J

1F T )

 

 

ˆ

 

J

1

(Fv

vF

T

)

0 .

 

 

 

 

 

 

 

 

 

(2.71)

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, в терминах (2.68) ˆ J 1F T , ˆ J 1Fv, ˆ .

42

3.Силы и напряжения в сплошной среде

3.1.Балансовые уравнения

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

Рассмотрим деформируемый материальный континуум, заполняющий

область пространства

ˆ

с границей

ˆ

V (t)

S (t) . Для произвольной тензорной

величины, определенной на материальном континууме и заданной объемной

плотностью ˆ

 

ˆ

 

ˆ

, записывается уравнение баланса

 

d ( xˆ,t) / dV

 

 

 

d

 

ˆ

ˆ

ˆ

ˆ

 

 

 

dt

 

ˆ dV nˆ dS

( ˆ i ˆ e )dV ,

(3.1)

 

 

ˆ

 

ˆ

 

ˆ

 

 

 

 

V (t )

 

S (t )

 

V (t )

 

 

 

 

 

 

 

ˆ

 

 

где nˆ ― внешняя нормаль к границе, а ― поток величины, втекающий

внутрь материала через границу, ˆ i (xˆ,t),

ˆ e ( xˆ,t) ― объемные плотности

внутренних и внешних источников рассматриваемой величины в текущей

 

 

ˆ

ˆ e

позволяют

учесть контактное

либо

конфигурации. Переменные и

полевое внешнее воздействие на материальный объем.

 

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

[4]

 

 

 

 

 

 

 

d

ˆ

 

ˆ

ˆ

 

 

dt

ˆ dV

 

( ˆ ˆ v)dV .

(3.2)

 

ˆ

ˆ

 

 

 

 

 

V (t )

V (t )

 

 

 

С помощью (2.6) выражение под интегралом в правой части можно преобразовать

ˆ ˆ ˆ v ˆ ˆ (v ˆ ) .t

Применяя к (3.1) с учетом (3.2) теорему Гаусса – Остроградского и

 

(

ˆ

ˆ

ˆ

ˆ

0 ,

t

(v ˆ

) ˆ i

ˆ e )dV

ˆ

 

 

 

 

 

V (t )

 

 

 

 

 

 

(3.3)

(3.4)

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

ˆ

ˆ

ˆ

ˆ e .

(3.5)

t

(v ˆ

) ˆ i

 

 

 

 

В лагранжевых переменных это уравнение будет выглядеть так

43

 

 

 

 

 

i

e ,

 

(3.6)

где ˆ J

1

ˆ

1

F , ˆ i

J

1

, ˆ e

J

1

 

 

,

J

 

i

e .

 

 

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

путем.

Рассмотрим произвольный

пространственный

ˆ

объем V

с

неподвижными границами

ˆ

 

в

пределах которого

располагается

S ,

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

 

d

ˆ

ˆ

ˆ

 

 

ˆ

 

 

 

dt

ˆ dV

nˆ ( v ˆ )dS

 

( ˆ i

ˆ e )dV

,

(3.7)

 

ˆ

ˆ

 

ˆ

 

 

 

 

 

 

V

S

 

V

 

 

 

 

в котором учтен конвективный поток

ˆ

 

 

ˆ

 

 

v величины вследствие движения

 

 

 

 

 

ˆ

 

 

ˆ

ˆ

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

неподвижную границу

S

.

Поток

 

 

величины ,

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

Ранее, в разделе 2.9, рассматривались уравнения совместности деформаций в дивергентной форме. Одно из этих уравнений, (2.71), следует

из (3.5) с учетом ˆ

J

1

F

T

ˆ

1

ˆ i ˆ e . Подставляя эти

 

 

, J

Fv,

выражения в (3.6) и используя (1.20), получаем второе уравнение ― (2.69).

3.2. Баланс массы в закрытой системе. Уравнение неразрывности

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

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

44

С помощью (3.5) и (3.6) можно записать уравнение баланса массы в локальной форме в эйлеровых и лагранжевых переменных. Для закрытой (не обменивающейся массой) системы

 

d

ˆ

 

 

 

dt

ˆ dV

0 ,

(3.8)

 

ˆ

 

 

 

 

V

 

 

где

 

 

 

 

 

 

ˆ

(3.9)

ˆ ( xˆ,t) dm / dV

плотность массы в актуальной конфигурации, поэтому в терминах (3.1)

будет ˆ ˆ,

ˆ

ˆ i ˆ e 0 ,

и мы

немедленно

получаем искомые

0,

уравнения

 

 

 

 

 

 

 

ˆ

ˆ

 

 

 

 

t

(v ˆ ) 0

(3.10)

и

 

 

 

 

 

 

 

 

ˆ

0.

(3.11)

 

 

ˆ ˆ v

Выведем балансовые уравнения для массовой плотности тензорной

ˆ

 

ˆ

 

 

 

величины ˆ ˆ d ( xˆ,t) / dm , которая требуется для записи балансовых

уравнений динамики и термодинамики континуума. Из (3.5) с учетом (3.10) следует

 

ˆ

ˆ ˆ

ˆ ˆ

 

 

ˆ (

 

ˆ e ,

(3.12)

t

v ) ˆ i

 

 

 

 

 

а из (3.6) с учетом (3.11) ―

 

 

 

 

 

 

ˆ ˆ

 

 

(3.13)

 

 

i e .

 

Уравнения (3.12) и (3.13) для массовых плотностей можно последовательно получить из (3.7) с учетом (3.10) или (3.11). Заметим, что из (3.12) и (3.13) невозможно вывести уравнения (3.10) и (3.11). Иногда для вывода (3.10) и (3.11) к локальной форме (3.7) (уравнению (3.5)), записанной для массовой

ˆ

 

 

 

 

 

 

плотности ,

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

( ˆ )

ˆ

ˆ

ˆ

 

 

 

t

(v ˆ ) ˆ i ˆ e ,

(3.14)

 

 

 

 

 

 

 

применяют тождество

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

( ˆ )

 

ˆ

ˆ

ˆ

 

 

 

 

t

(v ˆ ) ˆ

 

(3.15)

 

 

 

 

 

 

 

и последовательно получают (3.12) и (3.13).

 

 

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

 

 

 

 

( x) dm / dV .

 

(3.16)

 

 

 

 

 

45

 

 

 

Сохранение массы элементарного объема тела при его деформировании из отсчетной в текущую конфигурацию означает

 

 

(3.17)

dm dV dV .

Поскольку в силу (1.54)

 

 

dV JdV ,

 

 

то

 

 

 

 

(3.18)

J .

 

Уравнение неразрывности (3.18) — следствие закона сохранения массы. Дифференцируя (3.17)

 

 

 

 

 

 

 

 

 

0

dm dV

dV

 

 

 

 

 

 

с учетом (2.62)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV / dV

v ,

вновь получаем (3.11) ― уравнение неразрывности в дифференциальной форме в лагранжевых переменных, из которого нетрудно вывести его запись (3.10) в эйлеровой форме.

При помощи (3.10) можно выделить случай, когда плотность в точке пространства не зависит от времени:

 

 

(3.19)

( v ) 0 ,

а при помощи (3.11) ― случай, когда плотность в точке континуума не зависит от времени, т.е.

 

(3.20)

v 0 .

Последнее уравнение носит название условия несжимаемости.

 

3.3. Динамические структуры и балансовые уравнения системы материальных точек

Перед записью балансовых уравнений динамики континуума построим их для системы материальных точек [2] (на рис. 3.1 иллюстрируемой одной парой). Масса такой системы найдется как

m m ,

(3.21)

 

 

 

 

радиус-вектор ее центра масс

 

 

 

x

1

m x ,

(3.22)

 

 

m

 

скорость ее центра масс

46

 

 

v

 

1

 

m v .

 

(3.23)

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

С помощью x и v радиус-вектор и скорость любой точки представляется

 

 

 

 

 

 

 

,

v v v .

(3.24)

x

 

x x

 

 

 

 

 

 

 

 

 

 

 

 

 

Количество движения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q m v mv ,

 

(3.25)

главный вектор внешних сил

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

f

 

f

 

,

(3.26)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где f — сила действия на

-ю материальную точку других материальных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

l m ( xˆ y) v ( xˆ y) q ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.27)

 

 

 

 

 

 

m ( xˆ

y) v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и главный момент внешних сил

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

( xˆ

y)

f

 

 

( xˆ

y) f

 

( xˆ y) f + k,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(3.28)

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( xˆ y) f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

записываются

 

с

 

учетом

 

(3.24)

(3.26)

 

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

противодействия сил взаимодействия каждой пары точек системы f f

. В (3.27) и (3.28) выделяются моменты относительно центра масс системы —

собственный момент количества движения и собственный момент k .

Рис. 3.1. Система двух взаимодействующих материальных точек

47

В инерциальной системе отсчета для каждой материальной точки справедлив второй закон Ньютона

m a f ,

(3.29)

для всей системы материальных точек ― законы динамики Ньютона и Эйлера

q f , l m,

 

 

 

(3.30)

которые с использованием (3.25) ‒ (3.28) приобретают вид

 

 

 

 

q f , k .

 

 

 

(3.31)

 

 

 

Уравнение (3.31)2 получено вычитанием свертки ( x

 

y )

 

с (3.30)1 из

(3.30)2 и по существу является уравнением баланса собственного момента количества движения дискретной системы. Законы (3.31) независимы, однако следуют из законов (3.29). Внимание дискретной системе материальных точек здесь уделено по следующей причине. Нет оснований считать, что как угодно малая материальная частица не состоит в свою очередь из других частиц. С другой стороны, в нашем представлении среды, которую мы описываем континуумом, мы всегда ограничиваемся частицами определенного, хотя и малого по сравнению с размером тела масштаба, считая их неделимыми. Дискретная система материальных точек была рассмотрена выше как модель такой частицы. Для такой неделимой частицы (3.31)2 будет независимым от (3.31)1 фундаментальным законом. Такая модель частицы трактует ее саму малым деформируемым телом. Если, в частном случае считать моделью частицы абсолютно жесткое тело, называемое телом-точкой, то и в этом случае закон (3.31)2 будет независимым. В другом частном случае, когда континуумом описывается среда, состоящая из частиц — материальных точек, закон (3.31)2 пропадает. Такой континуум называют классическим, другие же, рассмотренные выше,

обобщенными: континуум, состоящий из тел-точек, называют континуумом Коссера, а континуум, состоящий из малых деформируемых тел, — микроморфным континуумом (Миндлина ‒ Эрингена).

3.4. Балансовые уравнения динамики континуума

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

48

Для каждого элементарного объема материала определена его масса dm, поэтому для всего тела B имеем

m dm .

(3.32)

B

 

Количество движения определяется как

 

q vdm ,

(3.33)

B

 

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

l

 

ˆ

 

y)

 

vdm

 

 

 

dm

 

ˆ

 

y)

 

vdm

 

B .

(3.34)

 

( x

 

 

 

 

 

( x

 

 

 

 

 

B

 

 

 

 

 

B

 

 

 

B

 

 

 

 

 

 

 

Здесь ( x y ) vdm — момент количества движения элементарного объема тела B, сосредоточенного в точке x , относительно точки y , dm — момент количества движения того же элементарного объема относительно точки y x , т.е. собственный момент количества движения элементарного объема. Вектор B следует считать частью собственного момента количества движения тела B, образуемой собственными моментами количества движения составляющих это тело элементарных объемов. Заметим, что первое слагаемое в (3.34) не разделяется на моменты количества движения тела относительно y и относительно его центра масс, так как впоследствии нас будет интересовать собственное вращение элементарного объема, а не тела. Для нахождения материальной производной q ,l нужно принять во внимание, что dV dm — масса элементарного объема материала и

.

 

 

 

 

 

 

 

 

 

поэтому dm 0 , тогда

 

 

 

 

 

 

 

 

 

 

 

q adm,

 

 

 

(3.35)

 

 

 

B

 

 

 

 

 

 

l

 

ˆ

y)

 

adm

 

B .

(3.36)

 

( x

 

 

 

 

B

 

 

 

 

 

 

 

Для каждого элементарного объема материала, сосредоточенного в x ,

определена действующая на него сила df

со стороны тела А,

поэтому для

всего тела B имеем

 

 

 

 

 

 

 

 

 

fB f ( B, A ) df df ,

(3.37)

 

 

 

 

 

B

 

B

 

где df — действующая на элементарный объем сила со стороны тела А и оставшейся части тела B. Для момента же имеем

mB m( B, A ) ( x

y ) df dk ( x y ) df kB .

(3.38)

B

B

B

 

 

49

 

 

Здесь ( x y ) df — момент силы df , действующей со стороны тела А на элементарный объем тела B, сосредоточенный в x , относительно точки y , dk — момент, действующий со стороны тела А на тот же элементарный объем относительно точки y x , т.е. собственный момент, действующий на элементарный объем, называемый еще моментом-парой. Вектор kB следует считать частью собственного момента тела, вызванной действием моментовпар. Заметим, что первое слагаемое в (3.38) не разделяется на моменты сил относительно y и относительно центра масс тела, так как впоследствии нас будет интересовать собственный момент, действующий на элементарный объем, а не на тело. Для получения выражений (3.37), (3.38) через внешние силы использован принцип действия и противодействия (3.44), (3.46) для телточек, доказываемый ниже.

В инерциальной системе отсчета для элементарного объема тела B

adm df

 

(3.39)

а для тела B

 

 

 

 

 

q f

B

,

l m

B

(3.40)

 

 

 

 

Из (3.40) с использованием (3.35) ‒ (3.38) получаем соотношения

 

q fB ,

 

B kB ,

(3.41)

аналогичные соотношениям (3.31) для системы материальных точек. Уравнение (3.41)2 является балансовым уравнением той части собственного момента количества движения континуального тела, которая вызвана действием моментов-пар. Это уравнение не следует из (3.39). Если же частица тела сводится к материальной точке, то для нее

0, k 0 ,

(3.42)

а уравнение (3.41)2 для тела пропадает.

Принцип действия и противодействия для материальных точек и тел,

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

действию внешних сил. В инерциальной системе отсчета

 

q

q

2

0

,

l

l

0 .

(3.43)

1

 

 

 

1

2

 

 

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

f2 f1.

(3.44)

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

m2 m1 .

(3.45)

50