
МММСС (От Кувыркина Г.Н.) / 33
.pdf
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
31 |
|
|
|
= 3p(σ1 |
− σ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
Непосредственным вычислением найти инварианты I1σ, I2σ, 31σ тензора напряжений
6−3 0
σij = −3 6 0 .
00 8
Найти главные напряжений для этого напряженного состояния и показать, что диагональная форма приводит к тем же самым значениям инвариантов.
Непосредственным вычислением находим:
I1σ = σii = 6 + 6 + 8 = 20.
I2σ = 1/2(σiiσjj − σij σij) =
=σ11 σ22 + σ22 σ33 + σ33 σ11 − σσ12 σ12 − σ23 σ23 − σ31 σ31 =
=36 + 48 + 48 − 9 = 123.

33
I3σ = det(σij = 6(48) + 3(−24) = 216.
Значения главных напряжений тензора напряжений σ1 = 9, σ2 = 8, σ3 = 3. Инврианты, подсчитанные через главные значения напряжений, равны
I1σ = σ1 + σ2 + σ3 = 9 + 8 + 3 = 20,
I2σ = σ1 σ2 + σ2 σ3 + σ3 σ1 = 72 + 24 + 27 = 123,
I3σ = σ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
|
|
dϕ |
|
ρϕ) |
|
|
||
ρϕ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 − bi◦ dV0 = 0, k, I, J = 1, 2, 3, |
||||||
|
|
− ∂aj |
|
||||||||
|
|
|
|
|
dvi |
|
|
∂ |
|
∂xi |
|
|
|
|
V0 |
|
|
|
|
|
|
|
|
где ρ0 и b◦i — плотность среды и проекции вектора 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 + λDVkk)δij + 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. При наличии моментов, распределенных по объему и/или по поверхности тела эта симметрия в общем случае отсутствует.