книги / Основы механики сплошной среды
..pdfРассмотрим линейное упругое тело [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, v° = 0) |
|
||||
|
|
|
dû |
|
|
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), в литературе их часто называют уравнения ми совместности в напряжениях.