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

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
19
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 2.6. Модель линейно-упругой среды

121

Исключение составляют процессы распространения упругих волн боль­ шой амплитуды в твердых телах, приводящие к значительному разогреву и даже плавлению и испарению твердых сред при ударе. К таким процессам относятся, например, удары частиц космического мусора, движущегося со скоростями 8 ... 16 км/с вдоль элементов космических аппаратов.

Если принять допущения (2.6.75), то из (2.4.3) получаем следующее выражение для плотности энтропии:

V= V0 + ^7Г+ (1Гр)а • • • • (в - а # ) .

(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, для сравнения температура плавления стали рз 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С (см. т. 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) = | Щ

+ щмО .

 

 

с шз = Л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 = • • а. для изотропной среды является шаровым, т. е. справедливы формулы (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