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

книги / Основы механики сплошной среды

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

Рассмотрим линейное упругое тело [1, 22, 28, 33,43, 54, 57, 58], т. е. среду, в которой тензоры о и £ связаны линейной, вообще говоря, анизотропной тензор-функцией. Общий вид та­ кой функции следующий:

= Сумей, или о = Q

(10.3)

Шесть независимых соотношений (10.3) носят название закона Гука для анизотропного упругого тела [25], а тензор четвёртого ранга С = Cijkiki ® kj ® /с* ® называется тензором модулей упругости. В силу симметрии <ту = ац он, очевидно, симметри­ чен по первым двум индексам, а в силу симметрии £ÿ = ец — по последним двум.

Пусть тензор-функция (10.3) является потенциальной, т. е.

существует скалярный потенциал деформаций W(Ê),

 

 

 

 

 

(10.4)

такой что

dW \

 

dW

 

_ 1,( 8W

, или a =

(10.5)

&тп — 2 1

+

.

\двтп

9бпт/

 

Ч&

 

Подставим выражение (10.4) в (10.5. Тогда имеем

деtj

& ТТ1П 2 Cijkl ^ dSmn dSvnn Ekl'j =

77X71

= T C 'ijk l ( 2 A f c i m n £ j j "Ь 2 Д ijm n ^ k l) ~ п

К ^-'m nkl^kh

где Д = Aijkiki ® kj <2>kk ® k[ — единичный тензор четвёртого ранга с компонентами

A{jkl — g (fiikfijl 4" fiüiïjk)-

( 10.6)

Легко проверить, что

о д

dajj

dajj

(10.7)

2Л« ы = ^

+ 5 ^

 

для любого тензора второго ранга & Из определения (10.6) кроме того следует, что для любого симметричного по первым и пос­ ледним двум индексам тензора четвёртого ранга А справедливы

соотношения

Aijkl&klmn — Д ijklAklmn Aijmn ИЛИ ^ • Д — Д А А,

( 10.8)

т. е. тензор Д в самом деле является единичным.

Из (10.4) следует, что тензор модулей упругости перестано­ вочен ещё и по парам индексов. Эта перестановочность эквива­

лентна равенству смешанных производных от W по

и а д

„ _ д а ц

d2W

d2W = д<гы =

(10.9)

^ijkl г\

л л

 

деы

деыдец

 

 

 

Cijkl — Сjikl Cijlk CkUj

(10.10)

Подсчитаем, сколько независимых компонент в Ж3 имеет тен­ зор модулей упругости, обладающий симметрией (10.10). В об­ щем случае в Ж3 тензор четвёртого ранга имеет З4 = 81 неза­ висимую компоненту. Поэтому для наглядности на плоскости набор этих компонент можно формально представить в виде матрицы 9 x 9 :

^С пн

Q122

С2211

Q222

С ззп

С3322

С1211

Q222

Q311

Q 322

Q l l l

C3122

Q111

Q l22

Q211

Q222

\C m i

Q 322

С и 33

C1112

Q233

Q212

Q333

C3312

Cl 233

C1212

Q333

Q312

Q133

Q112

Q133

Q l I 2

CO CO CM

Q212

C1333

C t312

Q l 23

C1131

Cl 121

Cl 132

С ш з \

Q223

C2231

C2221

Q232

С221З

C3323

C3331

C3321

C3332

C3313

C1223

Cl 231

Cl 221

Ci 232

C i213

Q323

Q331

C2321

C 2332

C2313

Q l 23

C3131

C3121

C3132

C3113

Q l 23

C2131

C2121

C2132

Q113

Q223

C3231

C322I

С3232

C3213

Q 323

Ci 33i

C l321

C1332

C1313/

В силу того что Сф 1 = Cjiki, последние три строки матри­ цы совпадают соответственно с её четвёртой, пятой и шестой строками, а в силу того что Cijki = Cijik, последние три столбца матрицы совпадают с её четвёртым, пятым и шестым столбцами. Следовательно, независимыми будут элементы, составляющие лишь верхний левый минор 6 x 6 . Наконец, условие Cijki = Сыц означает, что этот минор симметричен относительно главной диагонали. Симметричная же матрица шестого порядка имеет 6 • (6 + 1)/2 = 21 независимую компоненту.

Рассуждая аналогично, нетрудно установить, что в Ж2 тензор модулей упругости Q имеет шесть независимых компонент.

Обратный

закон Гука,

т.е.

тензор-функция,

обратная

к (10.3), записывается в виде

 

 

 

 

 

 

 

£ij = Jijkim,

или

£ = j / : 2 ,

(10.12)

где

■—*

—*

 

тензор

упругих податливос­

J = Jijkih ® kj ® kk® h —

тей.

Тензоры

четвёртого ранга

Q

и J

взаимообратны, т.е.

CijklJklmn JijklCklmn —Души> ИЛИ

Аналогично

(10.4)

введём

скалярный

ний w(a):

 

 

 

 

 

 

 

 

w ~ 2 kjkPblV ij.

такой что

 

 

 

 

 

Efnn

2

dut

+

dw

или

да.mn

 

 

 

danm

Q j/ —J : Q —Д.

(10.13)

потенциал напряже­

(10.14)

dw

e = -7—. (10.15)

да

Тогда типы симметрии тензоров J и (7 совпадают:

Jijkl Jjikl Jijlk Jklij-

(10.16)

Отметим, что объекты С и J являются материальными функ­ циями (в данном случае материальными тензор-константами)

определяющих соотношений (10.3), (10.12).

 

Если упругая среда изотропна, то компоненты

пред­

ставимы в виде линейной комбинации всевозможных свёрток символов Кронекера [36]:

C{jkl = ÀÔijSkl P&ik^jl P fiil&jk'

(10.17)

С учётом

симметрии (10.10) необходимо в (10.17) положить р, =

= р!, так

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

тропной среды имеют вид

r^ijkl ~ ^ij^kl "b p{&ik$jl

 

= Xôijôfcl "b 2//ДijkP

(10.18)

Подставим выражения (10.18) в

(10.3) и получим закон Гука

для изотропного материала:

 

 

 

a i j = X d ô ij

+

2ju .£ ij,

( 10.19)

где в — дилатация (5.23). Коэффициенты Л и д , называемые по­ стоянными Ламе, представляют собой независимые материаль­ ные константы определяющих соотношений (10.19) изотропного упругого тела.

Разобьём тензоры деформаций и напряжений на шаровые части и девиаторы (см. (5.25)):

E ij —

j ,

( Т у — & Ô ij +

3 ÿ ,

6

j

 

(10.20)

0 = €kb

o — g^fcjk,

et*; = 0,

Skk = 0.

Умножая обе части (10.19) на 5у, получим связь шаровых частей:

а = ( \ + ^ в = Кв,

(10.21)

где постоянная К — модуль объёмного растяжения-сжатия.

Умножая далее обе части равенства (10.21) на Æy и вычитая из (10.19), получим связь девиаторов g и е:

ау — 2peij.

( 10.22)

Постоянную Ламе р называют также модулем сдвига. Подставляя в (10.20) выражения в из (10.21) и еу из (10.22),

получим

= ш { у + h Sii= h ( " з л т ^ ; a k i+ "«)• (10-23)

Введём вместо Л и р так называемые технические постоян­ ные: Е — модуль Юнга H I / — коэффициент Пуассона, из соотношений

Ev

_

Е

А = (l + i / ) ( l - 2 и)'

 

2(1 + 1/)-

В конце этой лекции будет выявлен механический смысл упругих постоянных Е и и. Пока же отметим, что для встречающихся в природе изотропных упругих сред коэффициент Пуассона ме­ няется в интервале от 0 до 1/2.

Подставим (10.24) в (10.23) и после преобразований придём к обратному закону Гука для изотропного материала:

Eij = [—3i/crÆy + (1 + 1/)сгу]. (10.25)

Заметим, что из пяти уже введённых упругих постоянных — Л, р, К, Е, и — в качестве независимых можно выбрать любые две. Выражения для трёх остальных следуют из формул (10.21) и (10.24).

Найдём вид скалярных потенциалов W{E) (10.4) и w(a) (10.14) д л я изотропной среды. Пользуясь соотношениями (10.19)

и (10.25), запишем

 

 

 

 

 

W — g

+ 2fj£ij)eij — g Л и

- г t i t

i j c i j

 

 

 

 

 

 

 

= ^К в2 + fieijetj,

(10.26)

= 2 1 °‘ÿ [ - 3l/<T% + 0

+ 1У)£7у1=

9"

2

 

 

0-2 +

 

 

+

1 + м

3(1 —2M)

2 1 + v

 

 

2E

2E

’ a‘ + ~2E~Sii Sii'

(Ю-27)

Из (10.26) следует, что условия положительной определён­

ности квадратичной формы W(e) следующие:

 

 

 

 

К > 0,

ц> 0,

 

(10.28)

а из (10.27) видно, что квадратичная форма w(cr) положительно определена в каждом из двух следующих случаев:

 

1

 

(10.29)

 

Е > 0, —1 < м < -•

 

либо

У,

Ч

 

Е < 0,

и 6 (—оо; —1) U f g’ +°°J

(Ю.ЗО)

С учётом ранее сделанного замечания о значениях коэффициента

Пуассона

для

реальных

материалов ( 0 < м < 1 / 2 )

можно ут­

верждать, что

из систем (10.29), (10.30) физически реализуется

только первая.

ц = £/[2(1 + м)], К = Е/[3(1—2м)],

 

 

Так

как

то

усло­

вия (10.28)

совпадают

с объединением условий

(10.29),

(10.30). Это естественно, поскольку, как следует из (10.3), (10.4) и (10.12), (10.14), W(g) =w(g) для любого напряжённодеформированного состояния в упругом теле.

Итак, шесть определяющих соотношений линейного упругого материала — (10.3) либо (10.12), а в случае изотропии — (10.19) либо (10.25), замыкают систему (10.2), (5.5) пятнадцати уравне­ ний в области V

О начальных и граничных условиях уже говорилось в преды­ дущей лекции. В отличие от (9.27), теперь надо задавать и пе­

ремещения, и скорости

частиц в

начальный момент

времени:

4 = 0

и = й°(х),

— = v°(x).

(10.31)

Это вызвано тем, что уравнения (10.2) имеют второй порядок по времени. Граничные же условия на кинематической и

статической Es частях поверхности Е = dV

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

аналогии с (9.55), (9.56):

 

 

х б Е ц

й — üo(x,t).

(10.32)

х е Е s

S W = S0(x,t)-

(Ю.ЗЗ)

В компонентах векторное условие (10.33) имеет вид

х € Es

crijNj = Soi(x,t).

(10.34)

Как и ранее, если Eu = Е, то задача называется первой краевой, если Es = Е, то второй краевой, если Ev ф 0 и Es ф 0 — то смешанной.

Подставим соотношения Коши (5.5) в (10.19):

Oij — + pi^-ij W'j,г)I (10.35)

а затем (10.35) — в уравнения движения (10.2). В результате по­ лучим уравнения Ламе движения изотропного упругого тела [5]:

 

р

*= (Л + p)uj,ji + рАиi + pF{

(10.36)

Их можно написать и в векторной форме

 

 

Q2

 

 

 

р

- ^ = (\ + p ) g r z â d \ v û + р А и + рР.

(10.37)

Для записи статических граничных условий (10.34) в терми­

нах перемещений подставим в них выражения (10.35):

 

х € Es

\u kikNi + р

+ UjjNj'j - Soi{x, t).

(10.38)

Таким образом, начально-краевая задача теории упругости в перемещениях состоит в решении трёх уравнений Ламе (10.36) при выполнении начальных условий (10.31) и граничных усло­ вий (10.32), (10.38).

Если изучается не движение, а состояние покоя тела в поле внешних сил, т. е. (dujdt){x, t) = 0, то вместо уравнений движе­ ния записываются уравнения равновесия

+ pFi = 0,

(10.39)

а вместо динамических уравнений Ламе (10.36) — статические уравнения

Начальных условий (10.31) уже задавать не требуется. Соот­ ветствующие краевые задачи (10.39), (5.5), (10.19), (10.32), (10.34) и (10.40), (10.32), (10.38) называются статическими задачами теории упругости.

Если массовые и поверхностные силы зависят от времени столь незначительно, что силами инерции (правыми частями уравнений (10.2) и (10.36)) можно пренебречь, то говорят о ква­ зистатике и квазистатических задачах теории упругости.

Обратимся теперь к интегральному равенству (7.20), выра­ жающему теорему живых сил, или теорему об энергии в МСС. Для изотропной упругой среды изменение работы внутренних сил (7.19) имеет вид

MW = - ацецJ

dV = - J | ^E ij dV =

V

v v

= - 1 dW dV = - d WdV=-d<p (10.41)

vv

иявляется полным дифференциалом. Назовём величину

А®

WdV

(10.42)

 

v

 

работой внутренних сил, а интегральный оператор ц>— потен­ циальной энергией деформаций. Для изотропного тела

l:K92 + peijeij j dV.

(10.43)

Пусть массовые и поверхностные силы не зависят от переме­ щений. Тогда изменение работы внешних сил (7.21) также есть полный дифференциал и согласно (7.17), (7.18)

Л (е) = pF ■ÜdV + J

udZ.

(10.44)

E

 

 

Таким образом, каждое слагаемое в (7.20) является полным дифференциалом, поэтому можно записать первый интеграл тео­

ремы об энергии:

.,

(10.45)

 

К + <р-А& = С.

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

интегрирования соотношения (7.20) и, следовательно, является постоянной по времени величиной.

Пользуясь этим фактом, можно доказать теорему един­ ственности решения динамической задачи теории упругости. Действительно, предположим, что существует два решения й' и и" динамической задачи в перемещениях. Для разности этих решений й = й" — и' имеем однородную задачу, а именно:

однородные уравнения Ламе (F = 0)

д2щ

 

 

 

 

' dt2

(Л + p ) u j j i

+ р А щ ,

(10.46)

однородные граничные условия (гГо = 0, SQ 0)

 

х £ Eu

и = 0,

 

 

f £ Es ЛUk,kNi + р

+ UjjNj'j = 0

(10.47)

 

и однородные начальные условия (и 0 = 0, = 0)

 

 

 

 

 

t =

0 :

û — 0,

dt = 0.

(10.48)

Для рассматриваемой системы А ^

= 0 и из (10.45) и посто­

янства лагранжиана следует, что

 

 

 

Щ

+ ip(t) =

К(0) + <р(0).

(10.49)

Но кинетическая энергия /С (7.16), зависящая лишь от скоростей, в силу (10.48) равна нулю при t = 0, а при t > 0 неотрицательна. Потенциальная энергия tp (10.43), зависящая от деформаций, также в силу (10.48) равна нулю при t 0, а при t > 0 неотри­ цательна. Правая часть в (10.49) равна нулю, а следовательно, и левая часть этого равенства в любой момент t равна нулю, т. е. JC(t) = ip(t) = 0. Если упругие постоянные таковы, что квадра­ тичная форма W(s) положительно определена, то все скорости и деформации в теле равны нулю, или

du' _ dû"

(10.50)

dt dt

В случае второй краевой задачи равенства (10.50) и заклю­ чают в себе теорему единственности. В случае первой либо смешанной краевой задачи из (10.50) дополнительно следует, что

 

Рассмотрим одну из простей­

,лг(в“ 1

лгсч

ших краевых задач статической

гя-*-

теории упругости:

растяжение-

*1

сжатие стержня (рис. 38). Пусть

■ki—

Рис. 38

 

к

торцам однородного стержня

 

с

постоянным по

длине круго-

 

 

вым сечением площади Е приложена продольная сила X — Х к{ на Е^) и —X = —Хк\ на Е(2). Внешние нормали N М и N&) к обоим торцам и нормаль № 6оК) к свободной боковой поверхно­

сти имеют компоненты N^

= —N ^ = 6ц, /vj60^ = 0 ,

поэтому

граничные условия (10.34)

записываются следующим

образом:

я е Е * 1*, х €Е<2> а ц = ^ 6 ц

(10.51)

 

 

х € Е (бсж)

ai2N2 + cr&N3 = 0.

 

Поскольку на всей поверхности стержня задаётся вектор напря­ жений, данная задача является второй краевой задачей.

Решением уравнений равновесия (10.39) при отсутствии массовых сил, удовлетворяющим граничным условиям (10.51),

будет, очевидно, поле напряжений

X

<Тц= — = <То, <722 = ^33 = <712 = 023 = СГз| = 0.

(10.52)

Этому полю напряжений согласно обратному закону Гука (10.25) соответствует поле деформаций

<7о

иао

Л

£ц = -£>

£22 = £зз = — Ё~ = ~ Uen’

е *2 = £23 = £3i =

 

 

(10.53)

Из (10.53) выведем механический смысл технических посто­ янных Е и гл Модуль Юнга Е есть коэффициент пропорцио­ нальности между растягивающим (сжимающим) напряжением <то и продольной деформацией eu при одномерном растяжениисжатии стержня. Коэффициент Пуассона и, взятый со зна­ ком “минус”, есть коэффициент пропорциональности между по­

перечной £22 и продольной £ц

деформациями при одномер­

ном растяжении-сжатии стержня. Если сила X растягивающая

(X > 0), то

из физических соображений ясно, что стержень

удлиняется,

т. е. ец

> 0 ,

а если

сжимающая (X < 0), то стер­

жень укорачивается,

т. е.

ец < 0 .

Таким образом, модуль Юн­

га — величина положительная, неравенства (10.30) для реальных

упругих материалов невыполнимы и условиями положительной определённости W(g) и w(a) будут неравенства (10.29).

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

Рассмотрим, как и ранее, область V, занимаемую линейно упругим телом, с замкнутой границей Е. Трёх уравнений рав­ новесия в V относительно шести компонент тензора напряже­ ний а недостаточно для замыкания системы. Попытаемся её замкнуть. Для этого воспользуемся уравнениями совместности в виде (5.51).

Подставим в уравнения (5.50) выражения обратного закона Гука (10.25) для изотропной среды и после некоторых преобра­ зований придём к шести уравнениям совместности уже в компо­

нентах тензора напряжений:

 

 

(1 + 1/)Дсг^ + 3(7,ÿ = 31/Дa5ij + (1 +

+ Vjk.kù-

(10.54)

Из уравнений равновесия следует, что

aik,kj = —p^ï,j\

°~jkM =

= —pFjj. Таким образом, правую часть (10.54) можно выразить через F :

AtJij +

= “ Y ~ div F Sij - P(Fi.j + F3,i)-

(10.55)

Получены

недостающие для замыкания системы

уравне­

ния в области. Они носят название уравнения БельтрамиМичелла 0 .

Классическая постановка квазистатической задачи теории упругости в напряжениях [3,61] состоит в решении трёх уравне­ ний равновесия и шести уравнений Бельтрами-Мичелла (10.55) в области V относительно шести компонент <гу при удовлетво­

рении трёх граничных условий:

 

«i(2)ls, = « ?, <W,-|E2 = S°.

(10.56)

Уравнения Бельтрами-Мичелла (10.55) были получены для замыкания трёх уравнений равновесия в области V относительно шести компонент Оц. Однако самих уравнений (10.55) всего шесть, и, таким образом, в V имеются теперь девять уравнений относительно тех же шести неизвестных функций. На границе области Е выполняются всего три условия (10.56).

') Так же как и (10.54), в литературе их часто называют уравнения­ ми совместности в напряжениях.