книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 2.6. Модель линейно-упругой среды |
121 |
Исключение составляют процессы распространения упругих волн боль шой амплитуды в твердых телах, приводящие к значительному разогреву и даже плавлению и испарению твердых сред при ударе. К таким процессам относятся, например, удары частиц космического мусора, движущегося со скоростями 8 ... 16 км/с вдоль элементов космических аппаратов.
Если принять допущения (2.6.75), то из (2.4.3) получаем следующее выражение для плотности энтропии:
V= V0 + ^7Г+ (1Гр)а • • 4с • • (в - а # ) . |
(2.6.76) |
Подставляя (2.6.76) в (2.6.74), находим выражение для $ в адиабатиче
ском процессе: |
|
д = - х р - - е , х = — — °Л-------. |
(2.6.77) |
С ур — OL • • (Зво |
|
В частности, для изотропных сред, для которых тензор (3 шаровой (см. (2.6.69)),
$ = -x/fe, |
----- , |
(2.6.78) |
c v p |
— а [36о |
|
где s = е - -Е.
Из (2.6.78) следует, что изменение температуры $ в адиабатическом про цессе определяется только объемным расширением-сжатием упругой среды, причем при расширении, когда г > 0, упругая среда охлаждается: $ < 0, а при сжатии, когда г < 0, среда нагревается: $ > 0. Отметим, что это свойство справедливо только для идеальных сред, для неидеальных наблюдается иной эффект (см. § 5.3).
Приведем характерные значения $ для стали. Если |
|
|
||
cv = 0,3 • Ю3 ■Щ;. |
Р = 7 , 8 1 0 3 ^ . а = |
11 • 1СГ6 |
Е = 2 ■1011 |
Па, |
кг • К |
м3 |
|
|
|
и = 0,3, Ai = 1,15- |
1011 Па, А2 = 1,54 • 1011 |
Па, 90 = 293 К, (3 = 7,1 |
• 106, |
|
|
|
|
(2.6.79) |
то при характерном значении |
объемной деформации г = 0,01 (1 %) имеем |
|
х = 0,125 |
• 1СГ3 К2/Па, д = -8,9 К. |
(2.6.80) |
Следовательно, при заданных условиях сталь охладится всего на 8,9 К, т. е. $/$о = 0,03 <С 1, для сравнения температура плавления стали 6т рз 1900 К.
Таким образом, допущения (2.6.75) оказываются вполне приемлемыми для многих адиабатических волновых процессов в упругих средах.
Подставляя (2.6.77) в (2.6.2), находим тензор напряжений в адиабатиче ском процессе:
сг = 4С • • (е — Oi'ff) = 4С • • (е + ксх ® (3 - - е) = 4Са • • в, |
(2.6.81) |
где 4Са — тензор адиабатических модулей упругости,
122 |
Глава 2. Упругие среды с малыми деформациями |
|
|
4с а = 4С + х а • • 4С (X) /3 = 4С + я/3 (X) /3. |
(2.6.82) |
|
Используя (2.6.28) и (2.6.65), для изотропной среды получаем |
|
|
4Са = А!аЕ ® Е + 2Л2А, |
(2.6.83) |
|
А1в = А! + я р 2, |
(2.6.84) |
где Aia — адиабатическая константа упругости Ламе.
Для сталей, используя значения (2.6.79) и (2.6.80), находим характерные
величины Ai и Aia: |
|
Ai = 1,15 • 109 Па, х/32 = 6,3 • 109 Па, А1а = 121,3 • 109 Па. |
(2.6.85) |
Иначе говоря, для упругих сред отличие адиабатических 4Са и изотермиче ских 4С модулей упругости обычно невелико, и им часто можно пренебречь.
Учет различия между 4С и 4Са не составляет труда: для этого необходимо подставить (2.6.78) в уравнение движения (2.6.60) или, что то же самое, определяющие соотношения (2.6.81) в (2.2.5а), в результате приходим к постановке динамической задачи линейной теории упругости для адиабатиче
ских процессов:
(°ри = V - ( 4C a - - V c x ) u ) + pf,
< п • 4Ca • • V (g) u L = t ne, |
u L = ue, |
(2.6.86) |
IZja |
IZJU |
|
t = 0: u = u0, u = v0.
Сравнивая соотношения (2.6.60) для изотермических процессов и (2.6.86) для адиабатических процессов, заключаем, что эти постановки задач отлича ются только формальной заменой тензора изотермических модулей упругости 4С на 4Са, поэтому для решения задачи (2.6.86) можно использовать решение изотермической задачи (2.6.60) при д = 0 с заменой в полученном решении 4С —►4Са. В силу этого, далее не будем делать различия между 4С и 4Са.
2.6.9. Единственность решения для динамической задачи теории упругости
В п. 2.4.9 доказана теорема единственности для квазистатических задач теории упругости. Для динамической задачи линейной упругой среды имеет место аналогичная теорема.
Теорема 2.6.1. Для линейно-упругой среды, обладающей положительно определенным тензором модулей упругости, т. е. удовлетворяющей усло виям (2.4.75), решение динамической задачи (2.6.86) является единствен ным.
▼Пусть имеется два различных решения задачи (2.6.86).
Выберем произвольное непрерывно-дифференцируемое поле w (х,£), удо
влетворяющее условию w L |
= 0, и домножим на w уравнения движения |
'Xи |
|
§ 2.6. Модель линейно-упругой среды |
123 |
в (2.6.86), записанные для иа, а затем проинтегрируем их по У:
p(Ua - f) • W dV |
w • (V • 4C • • V <X>ua) dV = 0, a = 1, 2. (2.6.87) |
v |
v |
Преобразуем второй интеграл стандартным образом по формулам (2.4.3)- (2.4.4), тогда с учетом граничных условий задачи (2.6.86) получим
р(иа —f) • w dV + |
s(w) • • 4С • • s(ua) dV |
w • t ne dT, = 0, a = 1, 2. |
у |
у |
(2.6.88) |
|
|
Вычитая уравнение (2.6.88) при се = 1 из уравнения (2.6.88) при а = 2 и выбирая w = U2 —йь получаем
рй • w dV + е(й) • • 4С • • е(и) dV = 0, |
(2.6.89) |
уу
где и = U2 —ui и w = й. Из первого и второго интегралов в (2.6.89) можно выделить полные производные по t:
d |
о |u|2 |
^ e(u) • • 4C • • e(u) йП = 0. |
(2.6.90) |
|
dt |
^ |
|||
d V + i |
|
уу
Тогда, интегрируя по t, имеем |
|
|
о |й|2 |
е(и) • • 4С • • е(и) dV = const. |
(2.6.91) |
|
||
у |
у |
|
Поскольку и — разность двух решений задачи (2.6.86), то при t = 0: и = 0 и й = 0 в V, а значит оба интеграла в правых частях (2.6.91) обращаются в нуль, и константа в этом уравнении равна нулю. Тогда, поскольку первое слагаемое в (2.6.91) всегда неотрицательно (это кинетическая энергия К тела), то второе слагаемое должно быть неположительным, что противоречит условию (2.4.75) положительной определенности тензора 4С. Следовательно, принятое допущение о двух различных решениях задачи (2.6.86) является неверным, и решение этой задачи единственно. А
2.6.10. Следствие из условия устойчивости Адамара
Как отмечалось в п. 2.4.9, чтобы обеспечить единственность решения задачи МДТТ, для упругого потенциала Ф принимают условие устойчивости Адамара (2.4.73), эквивалентное условию существования положительного касательного модуля (2.4.76).
Пусть линейно-упругая среда обладает положительным касательным мо дулем. Тогда, подставляя выражение (2.6.1) для ф в (2.4.76), получаем, что
124 Глава 2. Упругие среды с малыми деформациями
тензор модулей упругости 4 С должен удовлетворять неравенству
h • • 4С • • h > 0 |
Vh ф 0, |
(2.6.92) |
ИЛИ |
|
|
Cijk% hkl > 0 |
Vhij, |
|
т. е. тензор 4С должен быть положительно определенным.
Если линейно-упругая среда обладает положительной касательной подат ливостью, то, подставляя выражение (2.6.96) в (2.4.77), получаем, что тензор упругих податливостей 4П должен удовлетворять следующему неравенству:
|
|
|
h • • 4П • • h > 0 |
Vh ф 0, |
|
|
(2.6.93) |
|||
ИЛИ |
|
U ^ h ifik i > 0 |
Vhij- |
|
|
|
||||
|
|
|
|
|
|
|||||
|
Выбирая в неравенстве (2.6.93) сначала тензоры h |
с |
отличными от |
|||||||
нуля |
только диагональными компонентами h{j = |
Sai5ajh a, |
а затем только |
|||||||
с |
ненулевыми диагональными |
компонентами: h{j |
= (l/2)(5ai5pj + |
|||||||
( а |
ф |
/3 ф у ф а ) , получаем шесть соотношений для упругих податливостей: |
||||||||
|
|
n QaaQ> 0 , |
n Q/3a/3> 0 , |
а ф |
ф , |
а , 1 3 = |
1,2,3. |
|
(2.6.94) |
|
|
Аналогичные соотношения можно получить из (2.6.92) |
и для компонент |
||||||||
pjijkl. |
|
|
|
|
|
|
|
|
|
|
|
|
> 0, |
Са^ |
> 0 , |
а ф |
ф |
а , ( 3 = |
1, 2, 3. |
|
(2.6.95) |
Однако для некоторых случаев анизотропии (в частности, для изотроп ных сред) границы допустимых значений упругих констант, вытекающие из (2.6.95), можно сузить, если предварительно перегруппировать в (2.6.92) слагаемые следующим образом. Подставим в (2.6.92) спектральное представ ление (2.6.41) для тензора 4С:
h . . 4C . . h = £ |
Г1 |
|
|
о |
i ^ ( a e .. h)(a/3. . h ) + У] |
c aah • • 4Та • • h |
|||
а,/3=1 |
а |
^ |
a = m + 1 |
|
|
= |
£ |
Cal3Ya(h)Yf3(h) + |
СааГ^Ъ) > 0. (2.6.96) |
|
|
о;,/3=1 |
а = т + 1 |
Здесь учтено определение линейных инвариантов Y^(h) симметричного тен зора второго ранга (см. (2.5.60), а также т. 1, (4.10.11)), определение ор
тогонального проектора |
(2.5.56): |
= 4Та • • h и свойство |
ортогональных |
|
проекторов: h • • |
= |
Y^(h) (см. т. 1, |
(4.10.15)). |
|
В силу взаимной ортогональности всех ортопроекторов Р ^ |
(а=1, ..., п), |
всегда можно выбрать такие значения тензора h, чтобы отличным от нуля был только один спектральный инвариант Y^(h) (для инвариантов I a \ h ) так можно сделать не всегда; например, так нельзя сделать для инвариантов
§ 2.6. Модель линейно-упругой среды |
125 |
Ji(h), / 2(h) в группе изотропии). Тогда из (2.6.96) следует, что должны выполняться неравенства
Саа > 0 , а = 1, ..., п. |
(2.6.97) |
Рассмотрим неравенства (2.6.94) и (2.6.97) для основных групп симмет рии.
Изотропные среды. Используя матричное представление (2.6.32) для ком понент тензора упругих податливостей, находим
п «««« = ! , й а/3а/3 = 2 ^ = ! ± Д а фр. а> (3 = 1,2,3.
Подставляя эти выражения в (2.6.94), приходим к следующим неравенствам:
Е > 0, G > 0, ту > —1. |
(2.6.98) |
Используя соотношения (2.6.45) и (2.6.48), для изотропной среды имеем
а „ = л , + | л 2 = А - = з ^ . |
а 22= 2л2 = ± . |
Подставляя эти выражения в (2.6.97), получаем еще одно дополнительное неравенство: v < 1/2 (второе неравенство снова приводит к G > 0). Объ единяя это неравенство с (2.6.98), окончательно получаем два независимых неравенства:
Е > 0, |
(2.6.99) |
ограничивающие допустимые интервалы упругих констант для изотропной линейно-упругой среды (неравенство G > 0 следует из (2.6.99)).
Ортотропные среды. Подставляя выражения упругих податливостей IPJ^ (2.6.17) через технические константы в (2.6.94), получаем
Еа > 0, Ga/3> 0, осф(3\ се,/3=1,2, 3. |
(2.6.100) |
Используя представление компонент Саааа через технические константы Еа и на(3 (см. упр. 1 к §2.6), из условий (2.6.95) получаем еще три неравен ства:
--------- V°pViot> 0, |
= |
1 — ^ 23^32 — ^ 13^31 — ^ 12^21 Н- |
2^ 12^ 23^ 31- |
( 2.6. 101) |
Вводя обозначение |
z1 |
= туа(3^(3а и учитывая, что |
все z1 > 0 |
(поскольку |
z1 = iy^p(Ep/Еа) > 0), из (2.6.101) получаем следующую систему неравенств:
1 > z1 > 0, |
1 + 2 y/ziz2 z3 > z\ + Z2 + Z3 . |
(2.6.102) |
|
Несложно проверить, |
что еще одна система неравенств, |
вытекающая из |
|
(2.6. 101): |
1, |
1 + 2y/z\z2zs < z\ + z2 + Z3 , |
(2.6.103) |
z7 > |
является несовместной.
126 |
Глава 2. Упругие среды с малыми деформациями |
Таким образом, из (2.6.100) и (2.6.102) имеем следующую систему нера венств для технических констант:
Еа 0, Gа(3 ^ 0 , 0 < |
^ 1? |
|
^12^21 + ^13^31 + ^32^23 < |
1 + 2^12^23^31 * |
(2.6.104) |
Трансверсально-изотропные среды. Подставляя выражения упругих по датливостей W ikl (2.6.26) через технические константы в (2.6.94), получаем
Е = Е х >0, E l = E 3 >0, G = G\2 > 0, Gl = G l 3>0. |
(2.6.105) |
Если подставить в неравенства (2.6.95) выражения упругих модулей Саааа через технические константы (см. упр. 3 к §2.6), то из (2.6.95) полу чим еще два дополнительных неравенства:
1 -VUV21 |
Q |
1 -^23^32 Q |
(2.6.106) |
|
Д„ |
|
Д„ |
v |
; |
A i/ = 1 - ^12^21 |
- |
2^23^32 + 21/12^23^32- |
|
|
Используя те же рассуждения, что |
и для ортотропной среды, из (2.6.105) |
и (2.6.106) имеем следующую систему |
неравенств для технических констант: |
Е\ > 0 , Е3 > 0, G\2 > 0, G\3 > 0 , |
0 < Щ2^21 < 1, |
0 < z/23^32 < 1, |
^12^21 + 2^23^32 < |
1 + 2^12^23^32* |
(2.6.107) |
Упражнения к § 2.6
Упражнение 1. Используя матричное представление (2.6.17) для тензора 4П и тензора 4С (см. т. 1, (4.6.64)), путем обращения матриц 6 x 6 показать, что компо
ненты тензора модулей упругости Cl^kl для ортотропной среды выражаются через технические константы (2.6.15) следующим образом:
С 1111 = Л, = Си = ^ (П 22Пзз - |
П |3) = f y i |
- |
^23^32), |
с 2222 = Л2 = С22 = ^(ПпПзз - |
П2з) = |Щ |
- |
^з^зО- |
С3333 = Л3 = С33 = 3-(П„П22 - |
П22) = |
- |
v{2v2{), |
|
С2233 = Л4 = С23 |
= ^(П 13П12 - |
П„П23) = | Щ |
2з + щмО . |
|
||||
|
с шз = Л5 = С\з |
= Т (п 12П2з - |
П22П13) = ^-(г/31 + г^г'зг), |
|
|||||
|
С1122 = Лб = |
С12 |
= -^(П 1 зП2з - |
ПззП12) = ^-{v\2 + ^1 3 ^3 2 ). |
|
||||
Ху1313 |
\ |
,55, |
Ху2323 |
\ |
Ху1212 |
|
\ |
, |
|
О |
= Ag = Ст1з = 0 |
L |
= А7 = |
O2 3 = O44, О |
|
= Ag = & \ 2 = Og6 |
|||
|
А„ = 1 - У'гз^зг - |
^13^31 |
- vn v2\ + 2 ^12^23^31. |
А = |
Д . |
|
|||
|
|
|
|
|
|
|
|
O1O2O3 |
|
128 |
Глава 2. Упругие среды с малыми деформациями |
Упражнение 8. Используя (2.6.28) и (2.6.35), показать, что тензор (3 = 4С • • а. для изотропной среды является шаровым, т. е. справедливы формулы (2.6.69).
Упражнение 9. Используя спектральное представление тензора модулей упругости
(2.6.41) и выражения для тензоров аа и 4Га (см. т. 1, (4.10.44) и (4.10.56)), а также
о ^
представления (2.6.19) и (2.6.20), показать, что соотношения между Сар и Сгэк1 для трансверсально-изотропных сред (группа Тз) имеют вид (см. также т. 1, (4.11.46)):
Си = О3333 = А2, |
С22 + С33 = 2Спп = 2 (A i |
+ 2 А 5), |
С33 = 2Ст2 = 2А5, |
||
С12 = С 2 д тг = у/2 С 2233 = V2 Аз, С 44 |
= |
2 С 1313 = 2С 2323 = 2А4, |
|||
для ортотропных сред (см. формулы т. 1, (4.11.56) |
и упр. 1 |
к §2.6): |
|||
|
|
о |
|
|
|
для изотропных сред: |
|
С'а(3 Оад, |
|
|
|
|
|
|
|
|
|
С\\ = Ai + |
- А 2 |
1111 + 2С1122) |
С22 = 2А2 = 2 С 1212. |
||
|
ГХ |
С |
|
|
|
§ 2.7. Фундаментальные решения в теории линейной упругости
Рассмотрим некоторые общие представления решения квазистатической задачи (2.6.63) для линейно-упругого тела.
2.7.1. Фундаментальное решение Кельвина
Проанализируем частный случай, когда линейно-упругое тело — все трех мерное пространство Eg, и в некоторой точке х = XQ этого пространства действует сосредоточенная массовая сила, для которой поле вектора f(x) имеет вид
|
f(x) = fo |
^(х _ х о)’ fo = /оёь |
(2-7.1) |
где fo — константа; |
— один |
из векторов декартова базиса; |
<5(х —XQ) — |
трехмерная дельта-функция Дирака (см. т. 1, формула (3.7.12)), удовлетво ряющая соотношению
А(ХЖХ - |
х0) dV = I |
(2.7.2) |
J |
10, |
х0 ^ V, |
V |
4 |
|
для произвольного тензорного поля А(х).
Рассмотрим решение задачи линейной теории упругости (2.6.63) для та
кого тела с нагружением (2.7.1) в изотермическом случае: |
|
|||||||
J V |
• (4С • • V |
0 11*4 ) + |
°р |
|
<5(х —х0) = о в q , |
|
||
l u ( fc)| |
= 0 , |
’ |
V x u W I |
la |
= 0 . |
1 |
j |
|
V. |
loo |
|
|
|
|
|
Здесь, поскольку рассматриваемая область Eg не имеет границы, для един ственности решения в качестве «граничного» условия вводим условие обра щения в нуль перемещения и его ротора V х на бесконечности. От-
130 Глава 2. Упругие среды с малыми деформациями
Согласно определению 3.7.3 из т. 1, производная 5-функции удовлетворяет
следующему интегральному соотношению: |
|
|
—V С) пА(х0), |
х Е У , |
(2.7.12) |
пА(х) ® V5(x —XQ) dV |
х £ V, |
|
n0, |
|
|
v |
|
|
V<5(x - х0) = -£J <5(X - x 0)el = Vj5(x - |
x0)r\ |
(2.7.13) |
Выражение (2.7.11) называют двойной сосредоточенной силой, поскольку вектор можно представить как предел отношения суммы двух сосре
доточенных сил f^ 5 (x —XQ —hei) —fg^5(x —XQ), приложенных к точкам XQ + hei и X Q , имеющих одинаковую величину, но противоположно направ ленных, к расстоянию h между точками приложения сил (рис. 2.7.2):
|
|
,ИЧ0/( Й ’ |
(2.7.14) |
|
|
|
|
fS i = |
(5)*)<5(х - хо - |
hei) - fok)S(x - х0)). |
(2.7.15) |
Если индексы к и I совпадают, |
то силу с плотностью |
называют |
|
двойной сосредоточенной силой без момента (рис. 2.7.2, а). Если же к и I |
|||
не совпадают, то силу с плотностью |
называют двойной сосредоточенной |
||
силой с моментом (рис. 2.7.2, б). |
|
|
Массовую силу в виде суммы трех двойных сосредоточенных сил без момента
(2.7.16)
а =\
называют центром расширения-сжатия или центром дилатации (рис. 2.7.2, в). Массовые силы, образованные комбинацией трех двойных сосредоточенных сил с моментом
f(t) = 4„. |
(2.7.17) |
называют центрами вращения (рис. 2.7.2, г) |
относительно координатных |
осей. |
|
Если подставить какое-либо из выражений для f(x) (2.7.11), (2.7.15)- (2.7.17) в уравнения равновесия системы (2.6.63), то получим аналог задачи (2.7.3) для всего пространства 8 %, например, для плотности силы (2.7.11):
V • (С • • V ® и<$) + p f f }-£J 8 {X - х 0) = 0 в q , |
(2.7.18) |
||
|
= 0, |
V X ujj? I = 0. |
|
ujj? I |
|
||
(/) loo |
’ |
( l ) loo |
|