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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

31

 

 

= 3p1

− σ2)2 + (σ2 − σ3)2 + (σ3 − σ1)2 = r

 

 

 

 

(67)

 

 

 

3I2s,

 

 

 

1

 

 

 

 

 

 

 

 

2

b

 

 

где I2s — второй инвариант девиатора напряжений s.

 

 

 

-

ных осейb

,

 

 

 

 

,

 

 

 

b

 

 

 

 

 

Для площадки

 

произвольно ориентированной относительно глав

 

 

 

 

 

 

 

 

3

3

2

 

 

 

 

 

 

XX

τn2 =

σα2 nα2

σαnα2 ,

(68)

 

α=1

α=1

 

где nα — направляющие косинусы внешней нормали к этой площадке. В случае не равных между собой всех главных напряжений, когда σ1 > > σ3, наибольшее касательное напряжение τmax = (σ1 −σ3)/2 возникает в площадке, параллельной оси OX2 и равнонаклоненной к двум другим главным осям, т. е. при n2 = 0 и n21 = n23 = 1/2. Действительно, при n2 = 0, обозначая n21 = z и учитывая, что в этом случае n23 = 1 − z, из (68) получаем τn2 = σ12z + σ32(1 − z) − σ1z + σ3(1 − z) 2 и в соответствии с необходимым условием экстремума τn2 по z приходим к указанному результату.

Пример 5.2

Тензор напряжений в окрестности точки P в декартовых осях Ox1x2x3 имеет компоненты

σij =

1

0

2

.

 

3

1

1

 

1

2

0

Определить главные напряжения и главные оси тензора напряжений, с которыми будет связана система осей координат Ox1x2x3.

Главные напряжения σ определяют из уравнения

3 − σ 1 1

1 −σ 2 = 0,

1 2 −σ

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

(σ + 2)(σ − 4)(σ − 1) = 0.

Главные напряжения являются корнями этого уравнения: σ3 = −2, σ2 = 1, σ1 = 4. Пусть ось Ox3 совпадает с осью главного напряжения σ3 и пусть n(3)i — направляющие косинусы этой оси. Тогда,

(3 + 2)n(3)1 + n(3)2 + n(3)3 = 0,

32

n(3)1 + 2n(3)2 + 2n(3)3 = 0, n(3)1 + 2n(3)2 + 2n(3)3 = 0.

Отсюда n(3)j = 0, n(3)2 = −n(3)3 , а так как nini = 1, то (n(3)2 = 1/2. Поэтому

n)3)1 = 0, n(3)2 = 1/2, n(3)3 = −1/2.

Пусть точно так же ось Ox2 соответствует главному напряжению

σ2. Тогда

2n(2)1 + n(2)2 + n(2)3 = 0, n(2)1 − n(2)2 + 2n(2)3 = 0, n(2)1 + 2n(2)2 − n(2)3 = 0,

так что n(2)1 = 1/3, n(2)2 = −1/3, n(2)3 = −1/3.

Пусть, наконец, Ox1 соответствует σ(3). Тогда

−n(1)1 + n(1)2 + n(1)3 = 0,

n(1)1 − 4n(1)2 + 2n(1)3 = 0,

n(1)1 + 2n(1)2 − 4n(1)3 = 0,

так что n(1)1 = −2/6, n(1)2 = −1/6, n(1)3 = −1/6. ]

Пример 5.3

Непосредственным вычислением найти инварианты I, I, 3тензора напряжений

6−3 0

σij = −3 6 0 .

00 8

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

Непосредственным вычислением находим:

I= σii = 6 + 6 + 8 = 20.

I= 1/2(σiiσjj − σij σij) =

=σ11 σ22 + σ22 σ33 + σ33 σ11 σσ12 σ12 σ23 σ23 σ31 σ31 =

=36 + 48 + 48 − 9 = 123.

33

I= det(σij = 6(48) + 3(−24) = 216.

Значения главных напряжений тензора напряжений σ1 = 9, σ2 = 8, σ3 = 3. Инврианты, подсчитанные через главные значения напряжений, равны

I= σ1 + σ2 + σ3 = 9 + 8 + 3 = 20,

I= σ1 σ2 + σ2 σ3 + σ3 σ1 = 72 + 24 + 27 = 123,

I= σ1 σ2 σ3 = (72)3 = 216.

]

6. Законы сохранения количества движения и момента количества движения среды

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

ddtv = dt Z

ρvdV = Z pdS + Z bdV.

(69)

Q d

 

 

 

 

 

 

V

S

V

 

Левая часть (69) соответствует скорости изменения во времени t век-

тора количества движения Qv среды, занимающей в актуальной кон-

фигурации область объемом V , ограниченным поверхностью S (ρ и v — плотность среды и вектор скорости ее частиц соответственно). Правая часть (69) равна главному вектору системы сил, характеризуемых в данном случае векторами p и b плотности поверхностных и объемных сил соответственно. Необходимо отметить, что (69) является основным постулируемым динамическим соотношением механики сплошной среды, аналогичным закону Ньютона в механике материальной точки.

Для средней части (69) с учетом (38) и (40) получим

dt Z

ρvdV = dt Z

ρvJ dV0 ==

d

 

d

 

VV0

 

 

 

ρJ )

dv

 

 

dv

 

V0

 

 

 

 

 

 

 

 

 

Z

 

vd(dt

 

+ ρJ dt

 

dV0 = Z ρ dt dV.

(70)

Отметим, что проведенное преобразование с учетом (43) и выражения dϕ/dt = ∂ϕ/∂t+ v ·rxϕ полной производной (нижний индекс x у диффе-

ренциального оператора Гамильтона rx означает дифференцирование

34

по пространственным координатам) для любой непрерывно диффе-

ренцируемой в области V функции ϕ(x,t) (скалярной, векторной или тензорной) позволяет записать

d Z

dt

V

 

 

 

ρϕ)

 

 

ρϕdV = Z

ρ

 

dV = Z

 

∂(

 

+ rx · (ρϕv) dV.

(71)

dt

∂t

V

 

 

V

 

 

 

 

 

Если векторы в (69) и (70) представить в проекциях на оси прямо-

угольной системы координат Ox1x2x3, заменив, согласно (56), проек-

ции pi вектора p = σ(n) на σjinj, i, j = 1, 2, 3, где σji — компоненты тензора σb напряжений Коши, а nj — направляющие косинусы единичного вектора внешней нормали к поверхности S, то, объединив (69)

и (70), с учетом теоремы Остроградского — Гаусса получим

Z

 

ρ dti

− bi

dV − Z

σjinj dS = Z

 

ρ dti

∂xj

− bi dV = 0, (72)

 

 

 

dv

 

 

 

 

 

dv

 

∂σji

 

V

 

 

S

V

 

 

 

 

где vi и bi — проекции векторов v и b соответственно. Отсюда следуют

уравнения движения среды

ρ

dvi

=

∂σji

+ bi,

или ρ

dv

= rx · σb + b.

(73)

dt

∂xj

dt

Эти уравнения выражают локальную формулировку закона сохране-

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

этих уравнений

∂(ρvi)

=

∂(σji − ρvivj)

+ bi

(74)

∂t

 

 

∂xj

 

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

dvi

=

∂vi

+ vj

∂vi

. В частном случае неподвижной среды (v ≡ 0) из

dt

∂t

∂xj

(73) и (74) следуют ее уравнения равновесия.

Пример 5.4

Тензор напряжений также можно разложить на два тензора: шаровой тензор (или тензор гидростатических напряжений) с компонентами σkk δij/3 = σ0 δij и девиатор напряжений sij, т.е.

σij = σ0 δij + sij .

Разложить тензор напряжений

 

 

12

4

0

 

σij =

 

4

9

−2

 

 

0

−2

3

35

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

Имеем

σ0 = σkk/3 = (12 + 9 + 3)/3 = 8,

тогда

0

8

0

+

4

1

−2

,

σij = σ0 δij + sij =

 

8

0

0

 

4

4

0

 

 

0

0

8 0

−2

−5

 

причем sii = 4 + 1 − 5 = 0. ]

Представим уравнения движения в материальных координатах. Для этого преобразуем левую часть (72) с учетом (38), (40) и (61):

Z

ρ0 dti

− bi

dV0

Z

 

∂aI TJI NJ dS0 =

 

dv

 

 

 

 

 

 

∂xi

V0

 

 

 

 

 

S0

 

 

 

 

 

= Z ρ0 dt

 

 

 

∂ak Tjk − bidV0 = 0, k, I, J = 1, 2, 3,

 

 

∂aj

 

 

 

 

 

 

dvi

 

 

 

∂xi

 

 

 

V0

 

 

 

 

 

 

 

 

где ρ0 и bi — плотность среды и проекции вектора b= J b·Hbт плотности объемных сил на оси Oai системы материальных координат; Hb

пространственный градиент деформации; NJ — проекции на оси Oai

единичного вектора N внешней нормали к поверхности S0, ограничивающей область V0, занятую средой в начальной конфигурации; Tjk

компоненты тензора напряжений Пиолы — Кирхгофа Tb. Отсюда сле-

дует локальная формулировка закона сохранения количества движения в этих координатах:

 

dv

=

 

∂xi

Tjk + bi= 0,

 

dv

= ra · (F

т

· T) + b.

ρ0 dti

∂aj

∂ak

или ρ0 dt

 

 

 

 

 

 

 

 

 

 

b

 

b

Пример 5.5

Определяющие уравнения для некоторой сплошной среды имеют вид

σij = (−p + λDVkkij + 2µDVij ,

 

 

 

 

 

 

 

1

 

∂vi

 

∂vj

 

 

где p — давление; λD, µD

— постоянные; Vij =

 

 

 

+

 

, vi

2

∂xj

∂xi

компоненты вектора скорости.

 

 

 

 

 

 

 

 

 

 

Уравнения движения (73) в данном случае

 

 

 

 

 

 

 

dvi

 

∂p

 

∂Vkk

 

∂Vij

 

 

 

 

 

 

ρ

 

= −

 

δij + λD

 

δij + 2µD

 

+ bi .

 

 

 

dt

∂xj

∂xj

∂xj

 

 

 

36

Так как Vkk =

 

∂vk

, 2

∂Vij

=

2vi

 

+

2vj

, то

 

 

 

 

 

 

 

∂xj∂xi

 

 

 

 

∂xk

 

∂xj

∂xj∂xj

 

 

 

 

dvi

 

 

 

∂p

 

 

 

 

 

2vj

 

∂vi

 

ρ

 

=

 

+ (λD + µD)

 

+ µD

 

+ bi .

dt

∂xi

∂xi∂xj

∂xj∂xj

]

Обобщение закона сохранения момента количества движения ма-

териальной системы на случай сплошной среды с учетом второго равенства (30) дает

dt Z x × ρvdV =Z

(mV +x × b)dV +Z mS dS+Z x × pdS, (75)

d

 

 

 

 

V

V

S

S

где x — радиус-вектор частицы сплошной среды; mV и mS — моменты, распределенные по объему V и по поверхности S. Эти моменты учитывают при построении математических моделей сплошной среды, называемых микрополярными. Используя (71), интеграл в левой части (75) можно преобразовать к виду

d

Z

x

×

ρvdV = ρ

d(x × v)

dV =

 

 

 

 

dt

 

 

 

 

 

 

Z

dt

 

 

 

 

 

 

 

 

 

 

V

 

 

V

= Z ρ

 

dt

× v + x × dt

 

dV = Z ρx × dt dV,

 

 

 

 

 

 

dx

dv

 

dv

 

 

 

 

 

V

 

 

 

V

 

 

поскольку (dx/dt) × v = v × v = 0. Последний интеграл в правой части (75) с использованием (58), теоремы Остроградского — Гаусса и

символа Леви-Чивиты eijk (m = 1, 2, 3) примет вид

Z x × pdS = Z eijkxjσmknmei dS = Z eijk

∂xm

ei dV =

 

 

 

 

 

 

∂(xjσmk)

 

 

S

 

S

∂xm

V

eijk σjk

 

+ xj ∂xm

ei dV,

= Z eijk

 

∂xm σmk + xj

ei dV = Z

 

 

∂xj

∂σmk

 

 

 

∂σmk

 

V

 

 

 

V

 

 

 

где ei — орты репера системы координат Ox1x2x3. Подставив преобразованные интегралы в (75), с учетом (73) получим интегральную форму

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

Z (mV + eijkσjkei)dV +Z mS dS =Z x ×

ρ dt − rx · σ − b dV = 0.

 

dv

b

V S V

37

Отсюда при mV = 0 и mS = 0 следует eijkσjk = 0, т. е. σij −σji = 0, что соответствует симметрии тензора σb, а с учетом (66) — и симметрии тензора Tb. При наличии моментов, распределенных по объему и/или по поверхности тела эта симметрия в общем случае отсутствует.

Соседние файлы в папке МММСС (От Кувыркина Г.Н.)