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

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

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

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

( xˆ2 xˆ1) fi 0,

i 1, 2 ,

(3.46)

означающее параллельность обоих векторов сил линии, соединяющей материальные точки. Соотношения (3.44), (3.46) эквивалентны ньютоновскому принципу действия и противодействия для материальных точек, тогда как (3.44) справедливо и для тел, описываемых классическим континуумом. В более общем случае взаимодействия тел-точек или состоящих из них тел уравнения (3.44) ‒ (3.46) требуют дополнительного предположения

k2 k1 ,

вчем читателю предлагается убедиться.

Вмеханике континуума силы и моменты разделяют на массовые и поверхностные. Если на элементарный объем массой dm в текущий момент времени действует сила df , то массовая плотность силы определяется как

df

/ dm f .

(3.47)

Часто также требуется относить

df к объему в

текущей или отсчетной

конфигурации:

df

 

 

(3.48)

/ dV f ,

df

/ dV f .

(3.49)

Последние формулы следуют из (3.17) и (3.47). Силы (3.48), (3.49) называются объемными удельными силами. Массовыми (или объемными) силами описывается действие на частицы тела внешних по отношению к телу электромагнитного и гравитационного полей, а также сил инерции в неинерциальных системах отсчета. Действие же частиц тела на соседние частицы, хотя и имеет электромагнитную природу, в механике континуума описывается контактными, поверхностными силами. Такими же поверхностными силами осуществляется контактное внешнее воздействие на рассматриваемое тело через его поверхность. Если на элементарную

материальную площадку площади

ˆ

 

dS в текущей конфигурации действует

сила dt , то поверхностной силой считают

 

ˆ

ˆ

(3.50)

t

dt / dS .

Иногда силу dt относят к площади dS

прообраза рассматриваемой

площадки в отсчетной конфигурации:

 

t dt / dS .

(3.51)

 

51

 

 

 

 

 

B

 

 

 

 

 

Итак, силу, действующую

на тело

из его

внешности

B , можно

представить как

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

ˆ ˆ

 

 

 

 

f (B, B) df

 

 

,

(3.52)

fdV tdS ˆ fdV

tdS

 

 

B

V

S

ˆ

ˆ

 

 

 

 

 

 

V

S

 

 

 

 

где S, S — поверхность тела в отсчетной и текущей конфигурациях.

Так же разделяют собственные моменты тела, вызванные действием на

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

 

dk / dm k ,

 

(3.53)

 

 

,

(3.54)

dk / dV k

dk / dV k

 

(3.55)

имеют электромагнитную природу либо появляются при переходе в неинерциальных системы отсчета. Поверхностные моменты-пары

ˆ

(3.56)

mˆ dm / dS ,

m dm / dS

(3.57)

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

 

 

 

 

 

ˆ

ˆ

(3.58)

 

 

 

 

 

k(B, B) dk kdV mdS ˆ kdV

mˆ dS .

 

 

B

V

S

ˆ

ˆ

 

 

 

V

S

 

3.5. Принцип напряжений

Сосредоточим внимание на описании внутренних локальных контактных взаимодействиях в сплошном теле B . Пусть сначала это тело представляет собой классический континуум. В инерциальной системе отсчета для него справедливо уравнение баланса количества движения. Из

(3.35), (3.40)1, (3.9) и (3.52) следует

ˆ

ˆ

ˆ

(3.59)

tdS

ˆ (a f )dV .

ˆ

 

ˆ

 

S

 

V

 

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

принципе разрезания Эйлера ‒ Коши:

Действие любой части тела на оставшуюся часть равносильно действию поверхностной силы по поверхности раздела.

52

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

материальную площадку ˆ ˆ . Предположим, что ndS

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

Нормаль здесь направляется наружу к рассматриваемой части тела. Это — постулат Коши. Он утверждает, что

ˆ

ˆ

(3.60)

t

t (xˆ,nˆ) ,

при этом зависимость от локальной геометрии поверхности игнорируется.

Докажем основную лемму Коши:

 

 

ˆ

ˆ

(3.61)

t (xˆ

, nˆ) t (xˆ,nˆ) .

Для этого возьмем прямой круговой материальный цилиндр малого радиуса dr , высоты 2 и образующей, параллельной n , разрезаемый точно пополам касательной плоскостью к поверхности раздела в точке x . При 0 правая часть (3.59) исчезает, а левая сводится к интегрированию по круговым сторонам диска, имеющим противоположно ориентированные нормали. В итоге имеем tˆ(xˆ, nˆ) (dr)2 tˆ(xˆ,nˆ) (dr)2 0 и (3.61). Сейчас оказывается

возможным доказать фундаментальную теорему Коши:

 

Существует линейный оператор , такой что

 

t n .

(3.62)

Здесь ( xˆ ) — определенное внутри тела поле тензора второго ранга, называемого тензором напряжений Коши, nˆ( xˆ ) — единичная нормаль к произвольно ориентированной малой материальной площадке, tˆ(xˆ,nˆ) — вектор поверхностной силы на этой площадке. Таким образом, определения девяти компонент тензора в точке среды достаточно, чтобы найти вектор поверхностной силы на любой площадке, проходящей через эту точку. Равенство (3.62) называют соотношением Коши. Если вектор-аргумент отображения (3.62) поставить справа (как привычнее), то в получаемых далее уравнениях локального баланса количества движения будет фигурировать транспонированный к тензор (а это неудобно).

О. Коши, который ввел тензор напряжений в 1822 г., доказал свою теорему так. В исследуемой точке среды выделяется элементарный материальный тетраэдр, три грани которого параллельны координатным плоскостям декартовой ортогональной системы координат, а четвертая — произвольна. Пусть n, t — внешняя нормаль четвертой грани и поверхностная сила на ней. Если ni — декартов орторепер, то внешние нормали первых трех граней и поверхностные силы на них суть ni и

53

tˆi ( nˆi ) tˆi (nˆi ) (использована основная лемма Коши). При стягивании тетраэдра в точку записанное для него уравнение (3.59) с точностью до малых более высокого порядка примет вид

 

 

 

 

tdS t dS 0,

 

 

 

 

 

 

 

 

i

i

 

 

 

где площади боковых граней тетраэдра

 

 

 

 

 

 

 

 

 

dS

/ dS

n n

,

 

 

 

 

 

 

i

 

 

 

i

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

n

ni ti .

 

 

 

Итак,

 

 

 

 

 

, то и

 

 

 

 

ni ti , а поскольку ti ti nj nj

 

 

 

 

 

 

 

 

ˆ

 

ˆ

ˆ ˆ

,

 

(3.63)

 

 

 

 

 

(t j

 

ni )ni nj

 

откуда

компонента ij

тензора

напряжений Коши

 

в произвольном

ортобазисе ni nj есть проекция поверхностной силы на площадке с нормалью ni на направление n j . Размерность компонент тензора напряжений Коши совпадает с размерностью компонент поверхностной силы — Па.

3.6. Локальная форма законов динамики

С помощью (3.62) и формулы Гаусса – Остроградского поверхностный интеграл в (3.59) можно свести к объемному:

tdS n dS dV ,

S

S

V

 

 

а затем переписать (3.59) в виде

 

 

 

 

ˆ

ˆ

0 .

(3.64)

 

( ˆ f ˆ a)dV

 

ˆ

 

 

 

 

V

 

 

 

Обсуждаемое уравнение справедливо для любой элементарной части тела, поэтому

 

 

 

(3.65)

f a .

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

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

 

 

 

 

f

 

 

 

 

a

 

 

(3.66)

x

tdS

x

 

 

dV

x

 

 

dV .

S

 

V

 

 

 

 

V

 

 

 

 

 

Из (3.65) получаем

54

xˆ

ˆ

ˆ

 

ˆ

 

ˆ

(3.67)

( )dV

xˆ f ˆ dV

xˆ aˆ dV .

ˆ

 

 

ˆ

 

 

ˆ

 

 

V

 

 

V

 

 

V

 

 

Вычитая (3.66) из (3.67), получим

 

 

 

 

 

 

 

ˆ

ˆ

xˆ

T

ˆ

0 .

 

xˆ ( )dV

 

nˆdS

 

ˆ

 

 

 

ˆ

 

 

 

 

V

 

 

 

S

 

 

 

 

Используя интегральное тождество

xˆ

T

ˆ

xˆ

 

nˆdS

ˆ

 

 

ˆ

S

 

 

V

ˆ

ˆ

,

(3.68)

( ) 2 dV

где

1

Є:

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

 

 

 

2

 

 

 

 

 

 

0,

(3.69)

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

Для классического континуума (3.69) есть локальное уравнение баланса собственного момента количества движения.

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

ˆ

ˆ

ˆ

ˆ

ˆ

 

f

 

k)

ˆ

ˆ

ˆ

 

a

 

)

ˆ

ˆ

(3.70)

( x

t

m)dS

( x

 

 

 

dV

( x

 

 

 

dV .

ˆ

 

 

 

ˆ

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

S

 

 

 

V

 

 

 

 

 

 

V

 

 

 

 

 

 

 

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

m n ,

(3.71)

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

произвольном ортобазисе ni nj есть проекция поверхностного момента-пары на площадке с нормалью ni на направление n j . Размерность компонент тензора моментов-пар совпадает с размерностью компонент вектора поверхностного момента-пары — Па м.

55

Уравнение (3.66), справедливое для любого континуума, вычитаем из (3.70) с учетом (3.71) и используем формулу Гаусса – Остроградского и следующее из нее тождество (3.68), в результате чего получаем

ˆ

0

( 2 ˆ k ˆ )dV

ˆ

V

или для любого элементарного подобъема тела

2

 

ˆk ˆ .

(3.72)

 

 

 

Заметим, что сейчас нет оснований говорить о симметрии тензора напряжений Коши. Для рассматриваемого здесь континуума в уравнении (3.65) фигурирует несимметричный . Более подробно уравнение (3.72) мы обсудим, когда приступим к отдельному изучению континуума Коссера. А до той поры будем заниматься исключительно классическим континуумом.

Для многокомпонентного классического континуума, согласно теории Трусделла, баланс количества движения (3.59) записывается для каждой компоненты в форме (3.7)

 

 

ˆ

σi

 

 

 

ivi dV

t

ˆ

 

 

ˆ

 

ˆ

 

V

 

S

ˆ

ˆ i

fi

ˆ

,

ˆ ivivi nˆdS

χˆ i dV

 

ˆ

 

 

 

 

V

 

 

 

где суммирование по двойным индексам i не производится, vi — скорость компоненты i, σi — напряжения, действующие на i компоненту, fi — вектор массовых сил, действующий на i компоненту. Локализация к произвольному пространственному объему дает

ˆ ivi

ˆ

ˆ

 

t

ˆ ivivi σi

ˆ i fi χˆ i ,

 

 

 

что с учетом уравнений баланса парциальных плотностей в отсутствии химических реакций

ˆ i

ˆ

 

 

ˆ ivi 0

,

t

 

 

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

 

divi

ˆ

 

 

divi

 

vi

ˆ

ˆ i

dt

σi

ˆ i fi

χˆ i ,

dt

 

t

vi vi .

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

ˆ dv m ˆ σ ˆ f dt

56

при условии уравновешенности внутренних сил i χˆ i 0 , если ввести

следующие обозначения

ˆ v m i ˆ ivi , σ = i σi ˆ i vi v vi v , ˆ f i ˆ i fi , ˆ i ˆ i ,

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

3.7. Уравнение баланса количества движения в отсчетной конфигурации

Иногда уравнение баланса количества движения элементарного объема в текущий момент времени относят к отсчетной конфигурации. Чтобы получить требуемое уравнение вместо (3.59) следует использовать

tdS (a f )dV .

(3.73)

S

V

 

Для того чтобы в левой части (3.73) перейти от интеграла по поверхности к интегралу по объему в K0 , сначала найдем тензор, связывающий t с n:

t n P .

(3.74)

Для этого примем во внимание (3.50), (3.51), откуда следует

ˆ ˆ

tdS tdS ,

соотношение Коши (3.62)

t n

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

(1.55)

ˆ

ˆ

ˆ T

ndS .

 

ndS JF

 

В результате получим

 

 

 

 

 

 

ˆ

 

 

t n (JF ) ,

 

т.е. тензор в (3.74) имеет вид

 

 

 

 

 

 

 

.

(3.75)

 

P JF

Этот тензор называется тензором напряжений Пиолы. Смысл его компонент позволяет прояснить соотношение (3.75), где он выступает как линейный оператор. Сравнивая (3.75) с (3.62) и вспоминая смысл векторов t и n, можно заключить, что компонента Pij тензора напряжений Пиолы P в

произвольном ортобазисе ni nj есть проекция силы dt , действующей на

элементарную площадку ˆ ˆ в текущей конфигурации (прообраз которой в ni dS

отсчетной конфигурации есть ni dS ) и отнесенной к dS , на направление n j .

57

Подставляя (3.74) в (3.73) и используя формулу Гаусса – Остроградского, приходим сначала к

( P f a)dV 0 ,

V

а затем к

 

P f a .

(3.76)

3.8. Мощность напряжений и закон сохранения механической энергии

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

U

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

 

 

(3.77)

K U P ,

в котором

K 21 v vdV ,

V

P f vdV t vdS ,

V

S

и будем полагать

U u dV ,

V

где u — массовая плотность внутренней энергии. Собирая все вместе,

d

 

1

ˆ

ˆ

ˆ

ˆ

 

dt ( 2

ˆ v v u)dV ˆ f vdV t

vdS .

(3.78)

 

ˆ

ˆ

 

ˆ

 

 

 

V

V

 

S

 

 

Перед взятием материальной производной слева (3.78) нужно перейти от интеграла по объему в текущей конфигурации (которая, конечно, зависит от времени) к интегралу по массе тела, а после него вернуться обратно:

dt (

2 v v u) ˆ dV dt (

2 v v u)dm (v v u)dm

 

d

 

1

ˆ

d

 

1

 

 

 

ˆ

 

 

 

B

 

B

 

 

V

 

 

 

 

(3.79)

ˆ ˆ

(v v u) dV .

ˆ

V

Мы будем полагать, что (3.78) записано в инерциальной системе отсчета, и потому справедливо (3.65):

f v ,

которое сворачиваем с v и интегрируем по объему текущей конфигурации, в результате получая равенство

ˆ

ˆ

ˆ

ˆ v vdV

[ ˆ f v ( ) v]dV .

ˆ

ˆ

 

V

V

 

 

58

 

Вычитая его из (3.78) с учетом (3.79), получим

 

ˆ

ˆ

ˆ

ˆ

ˆ

 

 

vdS ( ) vdV .

u ˆdV t

ˆ

 

ˆ

 

ˆ

 

V

 

S

 

V

 

Для преобразования интеграла по поверхности к интегралу по объему сначала используем соотношение Коши (3.62)

t n ,

а затем формулу Гаусса – Остроградского:

ˆ

ˆ

 

ˆ

 

ˆ

ˆ

vdS

nˆ vdV

( v)dV .

t

ˆ

 

ˆ

 

ˆ

 

 

S

 

V

 

V

 

 

Далее воспользуемся тождеством

 

 

 

 

 

ˆ

 

ˆ

 

ˆ

 

 

( v) ( ) v : v ,

 

в результате чего получаем

 

 

 

 

 

 

 

 

ˆ

ˆ

0

 

 

 

( ˆ u : v )dV

 

 

 

ˆ

 

 

 

 

 

 

V

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

(3.80)

 

 

u

:v .

 

Используя симметрию тензора напряжений Коши (для изучаемого сейчас классического континуума) (3.69) и определение тензора деформации скорости (2.18), (2.19), можно (3.80) переписать в виде

 

:D .

(3.81)

u

Левая часть (3.81) есть скорость изменения внутренней энергии, отнесенная к элементарному объему в текущей конфигурации. Мы видим в этом выражении в качестве силовой переменной тензор напряжений, благодаря чему u называют локальной мощностью напряжений. Выражение (3.81) подчеркивает, что локальный объем континуума взаимодействует со своим ближайшим окружением контактными силами, которые описывает тензор напряжений. Справедливость тождества сохранения механической энергии можно доказать на основе закона динамики (3.65) и определения скорости изменения внутренней энергии (3.81).

Аналогичным образом, опираясь на (3.74) и (3.78), можно показать, что

локальную мощность напряжений, отнесенную к прообразу dV

в отсчетной

конфигурации, можно выразить и как

 

u P:v .

(3.82)

Из (3.81) и (3.82) получаем записи массовой мощности напряжений:

u

1

: D,

(3.83)

 

 

 

 

 

u

1

P:v .

(3.84)

 

 

 

 

59

 

3.9. Меры напряжений, энергетически сопряженные мерам деформаций

Ранее нами вводились меры деформации и их скорости. Следуя похожей формальной методике [12], можно было бы сейчас ввести аналогичный набор мер напряжений. Любой мере деформации соответствует единственная мера напряжений, вместе с которой они должны давать

мощность напряжений. Например, для меры деформации e

правого

семейства находится мера напряжений t , такая, что

 

t:e u ,

(3.85)

где мощность напряжений единицы объема отсчетной конфигурации согласно (3.81) и (3.18) записывается еще и как

 

 

 

 

 

 

u : D ,

 

 

(3.86)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

(3.87)

тензор напряжений Кирхгоффа. Итак, вводимая мера

напряжений t

должна быть связана с мерой деформаций e условием

 

 

 

 

 

 

 

t e : D .

 

 

(3.88)

В этом случае меры деформаций

 

 

e

и

напряжений

t называются

энергетически сопряженными.

 

 

 

 

 

 

 

 

 

 

 

 

 

Записывая условие (3.88) в компонентах с использованием (2.40), (2.41)

t : e i f ( i )tii Dii

 

2 i j

(ei

ej )tij Dij = ij Dij

 

2

2

 

 

 

 

 

 

 

i j

i

j

 

 

 

и сравнивая коэффициенты при произвольных компонентах

Dij , получим,

что нормальные компоненты меры t в базисе правого трехгранника равны

 

 

 

 

ii / i

f ( i ) ,

 

(3.89)

а сдвиговые —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

2i

2j

1

 

 

 

при

i j

(3.90)

e

e

 

 

 

2

 

 

 

 

 

i

 

j

 

 

 

i

 

 

j

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

/ i f ( i )

при i j .

(3.91)

В (3.89), (3.90) компоненты

ij

 

записаны в базисе левого трехгранника. В

частных случаях получаем меры напряжений, сопряженные некоторым мерам деформаций семейства (1.43), имеющие компоненты:

m 1

ij / i j ,

(3.92)

m 1

ij i j ,

(3.93)

 

60