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

книги / Механика сплошных сред (теоретические основы обработки давлением композитных материалов с задачами и решениями, примерами и упражнениями)

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

1.4. ДИНАМИКА

Самостоятельно. Убедиться в правильности вычисления напряжений подстановкой их в уравнение равновесия (1.4.18).

Задача 1.4.2.3. По заданной функции напряжений Дж. Эри

Φ = a

x1x2

 

,

 

 

 

x2

+ x2

 

 

1

 

2

 

где a – константа с размерностью силы, определить компоненты тензора напряжений.

Решение. Распишем формулу (1.4.22) по цифровым индексам:

σ =

2Φ

= 2a

 

x x3

3x3 x

;

σ

 

=

2Φ

= 2a

x3 x 3x x3

;

 

 

1 2

 

1 2

 

 

 

1 2

1 2

x2

 

(x2

+ x2 )3

 

 

x2

(x2

+ x2 )3

11

 

 

 

 

 

22

 

 

 

 

 

 

2

 

 

 

1

 

2

 

 

 

 

 

 

 

 

1

 

1

2

 

 

σ = σ

 

= −

2

Φ

 

 

= −a

 

x4

6x2 x2 + x4

 

 

 

 

 

 

 

 

 

 

1

 

1 2

2

.

 

 

 

 

x x

 

 

 

 

 

(x2

 

 

 

 

 

 

12

 

21

 

 

 

 

 

 

 

 

 

+ x2 )3

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

1

2

 

 

 

 

Проверка. Вычислим слагаемые, входящие в уравнение равновесия (1.4.18):

∂σ

= 2a

x5

14x2 x3

+9x4 x

 

;

∂σ

 

= −a

2x5

28x2 x3 +18x4 x

;

11

 

2

 

 

1 2

 

1 2

 

 

 

 

21

 

2

 

 

1 2

1 2

 

x

 

 

 

 

(x2

+ x2 )4

 

 

 

x

 

 

 

 

 

(x2

+ x2 )4

 

1

 

 

 

 

1

 

2

 

 

 

 

 

 

2

 

 

 

 

 

1

2

 

 

 

∂σ

= −a

2x5

28x3 x2

+18x x4

 

∂σ

 

= 2a

x5

14x3 x2

+9x x4

 

12

 

1

 

1

2

 

1 2

;

 

 

22

 

1

 

1 2

1 2

.

x

 

 

 

 

(x2 + x

2 )4

 

 

 

 

 

x

 

 

 

 

(x2 + x2 )4

 

 

1

 

 

 

 

 

1

 

2

 

 

 

 

 

 

2

 

 

 

 

 

1

2

 

 

Нетрудно убедиться в том, что

∂σ11

+

∂σ21

0;

∂σ12

+

∂σ22

0.

x

x

x

 

 

x

 

 

1

 

2

 

1

 

2

 

1.4.3. Симметрия тензора напряжений

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

 

dV

 

 

σ× xdS = 0.

 

ρ F

 

× xdΩ +

(1.4.23)

dt

Ω

 

 

 

S

 

191

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

Здесь x – радиус-вектор материальной частицы. С помощью (П1.103) и (1.3.13) преобразуем поверхностный интеграл формулы (1.4.23) в декартовых координатах:

 

 

σ× xdS = σik

× xdS =

jmej

× x)

dΩ.

 

 

 

 

 

xm

 

 

 

 

 

S

 

 

 

Sα

 

 

 

 

Ω

 

 

 

 

 

 

После подстановки результата этого преобразования в (1.4.23) имеем

 

F dV ρ+ T

 

 

× xdΩ +

σ

jm

(e ×e

)dΩ = 0,

(1.4.24)

 

 

 

 

σ

 

 

 

 

 

m

 

 

Ω

 

dt

 

 

 

 

 

 

 

Ω

 

 

 

 

 

 

 

где учтено, что

 

x

=

xk

e

= δ

mp

e

p

= e .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xm

 

 

k

 

 

m

 

 

 

 

 

 

 

 

 

 

 

xm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первое слагаемое уравнения (1.4.24) при условии выполнения уравнения движения (1.4.16) равно нулю. Тогда, вследствие леммы об интегрировании по произвольному объему, из (1.4.24) получим

σjm (ej ×em ) = 0.

(1.4.25)

Здесь попарные векторные произведения ортов ej и em , учитывая (П1.14), запишем в виде ej ×em = jmp ep . Тогда из уравнения (1.4.25), представленного в

форме σ jm jmp ep = 0, следует закон парности касательных напряжений (боковых компонент тензора напряжений)

σik = σki.

(1.4.26)

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

 

Tσ = Tσт

(1.4.27)

в декартовых координатах справедливо для любого множества координат. Ясно, что тензор TΦ, обеспечивающий безусловное выполнение уравнения

равновесия (1.4.18), для получения с помощью (1.4.19) симметричного тензора Tσ должен быть симметричным.

192

1.4. ДИНАМИКА

1.4.4. Баланс мощности (работы)

Сначала рассмотрим мощность σαn VαdS, которую развивают поверхност-

Sα

ные силы, приложенные к телу Mα с объемом Ωα и поверхностью Sα. Используя (1.3.13) и (П1.103), запишем

σαn VαdS = σα Vα )dΩ.

(1.4.28)

Sα

Ωα

 

Упражнение 1.4.3. Доказать тождество

(Tσ V ) ≡ ( Tσ ) V + Tσ ( V ).

(1.4.29)

Упражнение 1.4.4. Используя симметрию тензора напряжений (1.4.26) или (1.4.27), доказать тождество

Tσ ( V ) ≡ Tσ Tξ ,

(1.4.30)

где Tξ определяется формулой Дж. Стокса (1.2.137).

Упражнение 1.4.5. Показать, что разложением тензоров напряжения Tσ и скоростей деформаций Tξ на девиаторную и сферическую части (П1.53)–(П1.56), формулу (1.4.30) при N = 3 можно представить в виде

Tσ Tξ = Dξ Dσ +3 σ0 ξ0,

(1.4.31)

где σ0 – среднее напряжение (1.3.20); ξ0 – средняя скорость деформации (1.2.148)

Подстановкой (1.4.29) в (1.4.28) при выполнении уравнения движения (1.4.14) получаем баланс мощности для α-среды:

σαn VαdS + ρα Fα VαdΩ =

Sα Ωα

= ρα

dVα

VαdΩ +

Τσα ( Vα )dΩ.

(1.4.32)

 

Ωα

dt

Ωα

 

 

 

 

В левой части уравнения (1.4.32) представлена мощность активных (поверхностных и массовых) по отношению к рассматриваемому телу Mα сил, а в правой – мощность реактивных (инерционных и внутренних) сил. Таким образом,

193

Vα , Vβ.

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

(1.4.32) следует рассматривать как баланс мощности всех действующих на тело Mα сил.

Теперь рассмотрим сплошную композитную среду M = Mα Mβ с объемом Ω и поверхностью S. Для объема Ωβ с поверхностью Sβ по аналогии с (1.4.32) запишем баланс мощности

σβn VβdS + ρβFβ VβdΩ =

Sβ Ωβ

= ρβ

dVβ

VβdΩ +

Τσβ ( Vβ )dΩ.

(1.4.33)

 

Ωβ

dt

Ωβ

 

 

 

 

Пусть в области Ωтензор напряжений Tσ, поле скоростей V , массовые силы F

иплотность ρтаковы, что в областях Ωα и Ωβ они принимают значения Tσα и Tσβ, Vα

иVβ , Fα и Fβ, ρα и ρβ соответственно. Кроме того, часть границ Sα и Sβ является общей границей Sαβ двух рассматриваемых тел, а оставшиеся части границ этих тел

S = Sα Sβ\Sαβ образуют замкнутую поверхность всего тела M. Баланс мощности всех действующих на тело M сил можно получить двумя способами.

Первый способ сводится к суммированию (1.4.32) и (1.4.33) с учетом обозначений параметров, действующих в объеме Ωи на его поверхности S. При суммировании приращений мощности, развиваемой поверхностными силами на соответствующих скоростях в точках s различных участков поверхностей, ограничивающих тела Mα и Mβ композитной среды M, необходимо учитывать действие на S

напряжения σn и скорости V , а на Sαβ– напряжений σαn , σβn и скоростей Прежде чем перейти к суммированию (1.4.32) и (1.4.33), рассмотрим зна-

чения мощности поверхностных сил на Sαβ. На основании (1.3.55) и (1.3.16)

имеем

 

σαn = −σβn s Sαβ

(1.4.34)

или

 

pαn = −pβn ; ταn = −τβn s Sαβ.

(1.4.35)

Ранее (1.2.182) отмечалось, что разрыв (скачок) вектора скорости КВ-поля скоростей на какой-либо поверхности может осуществляться лишь за счет

194

1.4. ДИНАМИКА

тангенциальной к этой поверхности составляющей, а нормальная составляющая вектора скорости должна быть непрерывной (1.2.183):

V p V p s S

ΔΕ

.

(1.4.36)

Ε

 

 

Если теперь просуммировать

мощности на SΔΕ

 

 

 

³ n V ςΕn VΕ )dS,

(1.4.37)

SΔΕ

 

 

 

то с учетом (1.4.34) эту мощность Рис. 46. Статические и кинематические условия на

можно записать в двух эквивалент-

общей границе двух тел

 

 

 

ных формах (рис. 46)

 

 

³ ςn (V VΕ )dS

(1.4.38)

SΔΕ

 

 

или

 

 

³ ςΕn (VΕ V )dS.

(1.4.39)

SΔΕ

 

 

Более подробно первое подынтегральное выражение (1.4.38) на основании (1.3.16) имеет вид pn (V p VΕp ) + Ωn (V Ω VΕΩ ). К аналогичному виду приводится второе подынтегральное выражение (1.4.39): pΕn (VΕp V p ) + ΩΕn (VΕΩ V Ω ).

Тогда с учетом (1.4.36) в обеих суммах соответственно остаются Ωn (V Ω VΕΩ ) и

Ωn (V Ω V Ω ). Отсюда, принимая во внимание (1.4.35), следует, что оба подынтег-

Ε Ε

 

 

 

 

 

 

 

 

ральных выражения одинаковы и для всего тела M в (1.4.37) пути вычисления мощ-

ности на поверхности SΔΕ через (1.4.38) или (1.4.39) эквивалентны. Для определен-

ности здесь введем обозначения:

 

 

 

 

 

 

 

 

n

Ω

n

 

Ω

V

Ω

Ω

(1.4.40)

ΩΔΕ

 

; VΔΕ

 

VΕ

и перепишем (1.4.37):

 

 

 

 

 

 

 

 

 

³ ΩΔΕn

VΔΕΩ dS.

 

(1.4.41)

 

SΔΕ

 

 

 

 

 

 

195

ταβn

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

При решении конкретных задач о движении композитных материалов обозначения типа (1.4.40) назначают в зависимости от свойств контактирующих сред.

Окончательно сумма (1.4.32) и (1.4.33), или баланс мощности всех действую щих на композитное тело M сил, имеет вид

σn VdS +

ταβn

Vαβτ dS + ρF VdΩ =

 

S

 

 

Sαβ

Ω

 

=

ρ

dV

VdΩ + Τσ ( V )dΩ.

(1.4.42)

dt

Ω

 

 

Ω

 

Здесь мощность активных по отношению к телу M сил представлена первым и третьим слагаемыми левой части уравнения (1.4.42), а мощность реактивных

сил – слагаемыми правой части уравнения. Учитывая, что векторы и Vαβ

всегда направлены на Sαβ противоположно друг другу, второе слагаемое в балансе мощности (1.4.42), характеризующее ее рассеяние на межкомпонентной границе, относят к мощности реактивных сил.

Если на всей границе Sαβ отсутствует скачок вектора скорости (1.4.40), то для всего композитного тела M можно использовать единое векторное поле скоростей V (1.2.95), соответствующее балансу мощности всех действующих сил

σn VdS +

ρF VdΩ =

ρ dV

VdΩ + Τσ ( V )dΩ.

(1.4.43)

S

Ωα

Ω

dt

Ω

 

 

 

 

 

С учетом (1.4.30) уравнение (1.4.43) можно представить в виде

 

σn VdS +

ρF VdΩ =

ρ dV VdΩ +

Τσ TξdΩ.

(1.4.44)

S

 

Ωα

Ω

dt

Ω

 

 

 

 

Вывод уравнений (1.4.42)–(1.4.44) вторым способом предлагается выполнить самостоятельно в следующем упражнении.

Упражнение 1.4.6. Исходя из записи мощности поверхностных сил

σn VdS +

ταβn Vαβτ dS,

(1.4.45)

S

Sαβ

 

действующих на композитное тело M = Mα Mβ, используя путь получения (1.4.32) из (1.4.28), вывести уравнение (1.4.42)

196

1.4. ДИНАМИКА

При выводе уравнения (1.4.42) предполагалось, что композитное тело M состоит из двух сред Mα и Mβ. Легко показать, что (1.4.42) обобщается для самого произвольного случая, когда M состоит из любого количества компонент. В этом случае во втором слагаемом баланса мощности (1.4.42) необходимо выполнить суммирование по греческим индексам α и β, как это делается во всем учебнике, полагая, что при α = β соответствующие слагаемые равны нулю. В частности, для гомогенных сред, когда α = β = 1, из (1.4.42) получаем (1.4.43) или (1.4.44).

Следует акцентировать внимание на том, что переход типа (1.4.28) с помощью (П1.103) возможен лишь тогда, когда поверхностный интеграл берется по замкнутой поверхности. Это требование вытекает из непременных условий теоремы М.В. Остроградского–К. Гаусса (П1.103).

Теперь приступим к выводу баланса работы всех действующих на гетерогенное тело M сил.

По аналогии с выводом баланса мощности (1.4.42), исходя из работы поверхностных сил на соответствующих перемещениях

σn UdS + ταβn Uαβτ dS,

(1.4.46)

S

Sαβ

 

получаем баланс работы всех действующих на рассматриваемое композитное тело M сил:

 

 

 

σn UdS + ταβn

Uαβτ dS + ρF UdΩ =

 

 

 

 

 

S

Sαβ

 

 

 

Ω

 

 

 

 

 

=

ρ d 2U2

UdΩ + Tσ ( U )dΩ,

 

(1.4.47)

 

 

 

Ω

dt

 

 

 

Ω

 

 

где скачок вектора перемещения на межкомпонентной границе

 

 

 

 

 

 

 

Uαβτ =Uατ Uβτ.

 

(1.4.48)

Если разрывы вектора перемещения внутри тела M отсутствуют, то баланс

работы (1.4.47) принимает вид

 

 

 

 

 

 

σ

n

UdS

+ ρF UdΩ =

ρ

d 2U

UdΩ + Tσ ( U )d

Ω .

(1.4.49)

 

 

2

S

 

 

Ω

 

Ω

dt

 

Ω

 

 

 

 

 

 

 

197

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

При описании движения гомогенной среды баланс работы имеет такой же вид. По аналогии с (1.4.30), используя симметрию тензора напряжений и кинематическую формулу О. Коши (1.2.70), можно показать, что

Tσ ( U ) Tσ Tε.

(1.4.50)

Поэтому баланс работы (1.4.49) обычно записывают, используя (1.4.50):

σn

UdS +

ρF UdΩ =

ρ d 2U

UdΩ +

T T dΩ . (1.4.51)

 

 

 

dt2

 

σ ε

 

 

 

 

 

 

 

 

 

S

 

 

Ωα

 

Ωβ

 

 

Ω

 

При решении многих задач ОМД, как отмечалось ранее, массовые и инерционные силы можно считать пренебрежимо малыми. В этом случае баланс мощности (1.4.42) применяется в упрощенном виде:

σn VdS + ταβn

Vαβτ dS = Dσ DξdΩ + 3σ0 ξ0dΩ , (1.4.52)

S

Sαβ

Ω

Ω

в котором использованы соотношения (1.4.30) и (1.4.31). При этом первое слагаемое в правой части уравнения представляет собой мощность, диссипируемую вследствие изменения формы композитного тела, а второе – вследствие изменения объема (1.2.151), (1.3.21), (1.3.22) (задача 1.4.4.2). Из (1.4.52) особо выделим тот частный случай, когда выполняется условие несжимаемости (1.2.98), и вследствие (1.2.146), (1.2.148) имеем

ξ0 = 0,

(1.4.53)

а на границе S тела M в поверхностном интеграле (1.4.52) рассматривается только мощность, развиваемая касательными поверхностными напряжениями. В этом случае из (1.4.52) получаем

τn V τdS + ταβn

Vαβτ dS = Dσ DξdΩ.

(1.4.54)

S

Sαβ

Ω

 

Такой случай для несжимаемых сред возможен, если на участках границы S с ненулевыми поверхностными напряжениями нормальная к этим участкам составляющая вектора скорости равна нулю. Важность этого частного случая состоит в том, что он довольно часто встречается при решении задач ОМД и баланс мощности с учетом (1.3.27) или (1.3.28) полностью определяется девиаторами напряжений Dσ и скоростей деформаций Dξ, касательными напряжениями

198

1.4. ДИНАМИКА

и вектором скорости. Аналогичные частные виды баланса работы можно получить из (1.4.47).

Задачи к пп. 1.4.4

Задача 1.4.4.1. По заданному в прямолинейной полосе 2 E1 2 ;

h E2 h полю скоростей V с компонентами

 

 

 

2

 

 

 

 

V1

= a 1

E2

 

; V

 

= 0,

h2

2

 

 

 

 

 

 

где a – константа с размерностью скорости, и полю тензора напряжений Tσ с компонентами

σ11 = σ22 = −p1 +

p p

2

 

2E

 

; σ12

= −

( p p )E

2

 

1

1

1

 

1 2

,

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где p1, p2 – давление, приложенное к границам E1 = 2 и E1 = 2 соответствен-

но, определить составляющие баланса мощности.

Решение. Начнем с вычисления мощности внутренних сил в (1.4.44):

Int = Tσ TξdΩ = σik ξkidΩ.

ΩΩ

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

ξ11 =

V1

= 0; ξ22

=

V2

= 0; ξ12

= ξ21 =

1

 

V1

+

V2

 

= −a

E2

.

 

 

E1

E2

2

E2

 

 

 

 

 

 

 

 

E1

 

h2

Подстановкой этих значений и заданных напряжений в формулу мощности внутренних сил получим

Int = (σ11ξ11 + 2σ12ξ12 + σ22ξ22 ) dΩ =

 

 

Ω

 

 

 

 

 

 

h

 

2a( p1

p2 )h

 

 

2

 

 

= 2 ∫ ∫

σ12ξ12dE1dE2 =

.

 

3

 

 

h

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

199

1. МЕХАНИКА СПЛОШНЫХ КОМПОЗИТНЫХ СРЕД

Теперь вычислим мощность внешних, в данном случае на поверхности S, сил в (1.4.44):

Ext = σn VdS = ( pnV p + τnV τ ) dS,

S S

где pn, τn – модули нормального и касательного напряжений на поверхности S с

единичной внешней нормалью n ; V p , V τ – модули нормальной и касательной составляющих вектора скорости на поверхности S.

На границе E1 = 2 нормаль n имеет компоненты n1 = –1; n2 = 0. Поэтому по формуле (1.3.14) с учетом заданных напряжений находим модуль нормального напряжения pn = – p1, а по формуле V p =V n с учетом заданного поля скоростей имеем модуль нормальной составляющей вектора скорости

V p = −V1

 

 

E2

 

 

= −a 1

 

2

 

. На этой же границе по формуле (1.3.29), учитывая, что

h

2

 

 

 

 

 

 

диагональные компоненты девиатора для заданных напряжений равны нулю,

находим модуль касательного напряжения τn = −σ12 = ( p1 p2 )E2 , а с помо-

щью формулы (1.2.175) находим V τ = 0.

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

 

 

 

 

 

2

 

 

 

ложной границе E1

=

 

имеют вид pn = – p2; V p = V1

= a 1

E2

 

,

а на грани-

2

2

 

 

 

 

 

h

 

 

 

 

= ± h p

n

= − p

+

p p

 

 

1

2E

 

; V p = 0; V τ = 0;

τ

n

= ±

( p1 p2 )h

цах E

 

1

2

 

1

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соответственно.

Теперь, используя найденные условия на границе, вычисляем мощность внешних

сил:

 

 

 

 

 

 

 

 

 

h

 

 

E2

 

 

2a( p

p )

 

Ext = a( p1 p2 )

1

 

2

dE2

=

1

2

.

h

2

3

h

 

 

 

 

 

Подставим в формулу баланса мощности (1.4.44) без учета массовых и инерционных сил

JБ = Int – Ext = 0

200