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

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

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

§ 2.5. Нелинейно-упругие среды с малыми деформациями

101

 

р ( сг) _

р ( £)

 

(2.5.71)

 

О)

 

у(е) («)

 

 

 

 

а

 

 

укороченное спектральное представление, или

 

гп

 

 

п

<Ра р(е)

 

= Е

аа

^

^

(2.5.72)

уС

а=1

 

 

а=т+1

а

 

полное спектральное представление, где

 

 

YY = ра(х5\

..., уп<е>) = дъ/дг^)

(2.5.73)

скалярные функции от спектральных инвариантов; Ф — потенциал, представленный как функция этих инвариантов.

▼В т. 1, гл. 4 было установлено, что для всех групп Gs полный функцио­ нальный набор скалярных инвариантов тензора второго ранга можно предста­ вить в виде алгебраических инвариантов 1^ не выше третьей степени, т. е. линейных, квадратичных и кубических инвариантов. Следовательно, потен­ циал ф (е) можно представить как функцию (2.5.21) от этих инвариантов и определяющие соотношения идеально-упругой среды можно записать в виде (2.5.22). Поскольку, по предположению, среда является квазилинейной, то Ф не может зависеть от кубических инвариантов l!f\e ), так как их производные dl!f^/де дают тензорно-квадратичные слагаемые в (2.5.22), что запрещено в силу теоремы 2.5.2. Тогда если Ф зависит только от линейных и квад­ ратичных инвариантов тензора в, то в качестве таковых можно выбрать не l\ s\e ), а спектральные инварианты Y1(e) в соответствующей группе Gs, т. е. рассматривать Ф как функцию вида

Ф = Ф Ц /Д

...,

Ц(£)).

(2.5.74)

Тогда определяющие соотношения (2.5.71) примут вид

 

_

5Ф _ ^

дф

dY/e)

(2.5.75)

а

де

 

де

а

 

 

а = \

 

 

Определим производные d Y ^ /де. Очевидно, что имеют место следующие

соотношения:

1

dY/£) 2 _

1

p (g)

дГ (а)

 

 

 

 

 

 

( е )

 

де

2

де

(е)

(

де

(2.5.76)

у(£)

1а)

 

 

а

 

± а

 

 

 

Вычислим тензорную производную от левой и правой частей (2.5.55) по е:

п др ( е )

 

Е

(<*)

(2.5.77)

де

а = \

и умножим скалярно получившуюся формулу на Р (*).

102

 

Глава 2. Упругие среды с малыми деформациями

 

р (е) _

TD(£)

П

 

4-р

 

4-р

 

 

= Е

 

 

 

(Р )

Е рW)

Г (/9) ' '

Г («) =

 

 

а =\

а =\

 

 

 

 

 

 

 

 

= £ ■ 4Г ( й| ■ • 4Г (й)(а)

= Р $

• ■(№<*>/*)■

(2.5.78)

 

 

 

 

 

 

(а) ’ ’ V— («)

 

Здесь использовано свойство (2.5.57) ортогональности тензоров 4Г(а).

Подставляя формулу (2.5.78) в (2.5.76), находим

 

 

 

 

д ¥ & /д е = *>%/¥&,

а =

1 ,

П.

(2.5.79)

 

 

(а)

 

 

 

 

 

 

С учетом (2.5.79) определяющие соотношения (2.5.75) примут вид

 

 

 

1

д/Ф_р(е)

 

(2.5.80)

 

 

у(£) Яу(е)

(а)’

 

 

 

 

 

 

 

а=1

а

а

 

 

 

 

Если выделить в (2.5.80) линейные инварианты при а = 1, ..., ш и учесть формулу (2.5.59), то получим в точности представление (2.5.72).

Поскольку для тензора сг имеет место свое спектральное представление (2.5.54), то, сопоставляя (2.5.54) и (2.5.80), в силу свойства (2.5.57а) взаим­ ной ортогональности, получаем, что имеет место и укороченное спектральное представление (2.5.71). Действительно, пусть каждый ортопроектор p|^j вы­

ражается через все ортопроекторы р |^ : p|^j =

^ /сад Р |^ . Тогда, умножая

это соотношение на р |^

(а ф у), находим

 

Г) __ ■р(СГ)

т э ( сг) __ и т э ( £ )

т э ( £ )

и -

*(а)

' ' ^(7) ” ^7-^(7)

• • -^(7)’

т. е. ка1 = 0, следовательно, отличными от нуля являются только каа. Умножая соотношение (2.5.71) само на себя, в силу (2.5.59) и (2.5.61),

получаем

(2.5.81)

 

УУ)2 = УУ)2Га1У У 2-

Отсюда

следует, что ipa = уУ \ т. е. действительно

имеет место формула

(2.5.73).

 

С помощью теоремы 2.5.3 определяющие соотношения (2.5.21) или (2.5.46) для изотропной квазилинейной среды можно представить в виде

(2.5.72):

- ^

,

^(е.еи)

 

„ _ 1

(2.5.82)

ст— —=

еи)Уь-\-------

V3

 

 

ei

 

(в (2.5.46) и (2.5.82) функции у7 различны).

 

Укороченное спектральное представление (2.5.71) имеет вид

 

Р $ = ( 1 / л / З Щ Ц

е и ) В ,

P g

= Щ ( е , £и) / Щ Р ( 2) 4

(2.5.83)

а соотношения (2.5.72) между инвариантами —

 

а = л/3 Lp\{e, еи), еи = (р2(е,£и)-

(2.5.84)

§ 2.5. Нелинейно-упругие среды с малыми деформациями

103

Здесь выполнена замена обозначений:

 

<р{2Зе, еи) -►<р(е, £и).

 

Применяя теорему 2.5.3 для трансверсально-изотропной

квазилинейной

среды, получаем, что определяющие соотношения (2.5.21) или (2.5.49) можно представить в виде

 

ст=

-2 ,

1 _ г™ - 2ч

+

Ч>Ъ р(£) I ^ 4 „ ( е )

(2.5.85)

 

+

- ^ ^ ( Е -

с|)

Ч

(е)

(3)

^л(е) (4)’

 

 

 

V 2

 

 

 

 

м

 

где

скалярные функции инвариантов тензора е:

 

 

 

 

 

 

О)

 

У,(£)>

(2.5.86)

 

Для ортотропной квазилинейной среды определяющие соотношения

(2.5.21) или (2.5.52) можно представить в виде

 

 

 

 

 

з

-2

 

 

б

VOL 0 (e)

 

 

 

 

Е

 

|

 

(2.5.87)

 

 

 

<Р*Са + 2^

у(£)

(°0’

где

 

 

а =\

 

 

а=4

 

 

 

 

 

 

 

 

 

 

 

 

(2.5.88)

 

 

 

^

 

 

 

. . . . n

w )

(подчеркнем еще раз, что функции уу* в (2.5.52) и (2.5.88) различны). Варианты конкретного задания скалярных функций (2.5.84), (2.5.86) и

(2.5.88) приведены в т. 1.

2.5.9. Следствие из условия устойчивости Адамара

Вспомним, что для обеспечения единственности решения задачи МДТТ в п. 2.4.8 было введено понятие сверхустойчивости среды по Адамару (эл­ липтичности). Это понятие для идеальных сред эквивалентно требованию (2.4.75) или (2.4.76) существования положительного касательного модуля.

Учитывая, что потенциал Ф(е,0) является функцией (2.5.21) инвариантов

l\ s\e ), требование (2.4.76) можно записать в виде

 

 

УЗ

^ а ^ / З

+ УЗ (p^h2 ^

т Ь2 > 0,

(2.5.89)

а,(3=1

 

 

7=1

 

 

 

 

где

 

 

 

 

 

 

 

дгЧ!

 

dlis)

—о

д 2 Т^

о

(2.5.90)

фа(3= д1(а8)д ф У ka

Ь

де

h2=

h • • —

• • h,

h2= h • • h,

1

ds

 

 

h — произвольный ненулевой симметричный тензор второго ранга. Неравенство (2.5.89) накладывает ограничения на вид потенциала Ф. Более детально эти ограничения будут рассмотрены для линейно-упругих

сред в § 2.6.

104

Глава 2. Упругие среды с малыми деформациями

Упражнения к § 2.5

Упражнение 1. Показать, что для изотропной несжимаемой квазилинейной упругой среды определяющие соотношения можно представить в виде

а = Е - (р2е, ip2 = д

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

а = —рЕ + ^2сз + 40 • • е,

д

r(3)/

V>7 = ^ * ( 4 3)00.

1Г(е),в),

а для ортотропной несжимаемой квазилинейной упругой среды —

<7 = -рЕ + 5 >

7г Щ 4о г . . е> ^ =

.... у ° \ /<°>),

7=1

 

Т

к которым следует присоединить соотношение (2.5.42).

Упражнение 2. Доказать, что спектральные инварианты (2.5.61)

действительно

являются инвариантами относительно группы Gs.

 

Упражнение 3. Показать, что в базисе сi спектральные представления (2.5.62) для группы I (изотропная среда) имеют следующий вид:

 

<?ij —(l/3)(T(^j ~\~Sij,

&ij —(1/3)s5jjj ~\~Cij,

где s~ij = P(2)ij>

GJ = P(2)ij ~

компоненты девиаторов, а инварианты имеют вид

сг = 1\{ст) = ап

+ & 2 2 + ^зз,

£ = I i(s) = €и +?22 + £33,

= SijS13 = |((стц

- ст22)2 + (ст22 -

Стц)2 + (стп - ст33)2 + 6(ст?2 + ст|3 + of3)),

£ и = e i j г 3

= 2 ^ £ п

— £ 2 г)2 + ( £ 2 2 —

£ ц ) 2 + ( £ ц — е з з ) 2 + 6 ( £ 22 + £ 2з + £ 2з ));

для группы Тз (трансверсально-изотропная среда):

Ртц = ^ З зМ ^ З ,

Р(2)Ц(2)ii

=- 2 ^

 

 

- Si36j3),

P & =

& i j -

^(<711

+ ^ 2 2 ) № j - S i s S j s )

+

( f s s f i i s d j S - ( ^ i 3 ^ j 3 + 7 '3 ^ г з ) .

 

■P(4)^ = ^ ] 2№17> + M il) + a23(^i2^j3 + ^i.3^j2)>

 

 

F,(CT)

= <T33,

У2(ст)

= (1Д /2)(а„

+ a 22),

У3(ст) =

(1 ^ ) ( ( а ц - а 22)2 + 4 а 22)1/2>

Y4(<7)

= \/2 (of3 + 5f3)1/2;

для группы О (ортотропная среда):

 

 

 

 

э(а) О

ja.i ^ —1?2, 3,

 

—cr23(^2^j'3 д~ $i2>$j2)i

тэ(а

 

 

 

P(Q)ij = ^\2^i\^j2 + ^ 2^l),

*(5)ij = ^13(^l^j3 + SisSjl),

Y ^ = d aa,a =

1,2,3;

Г У

=C2\a23\,

Y p ] = ^2 |a13|, Y6(ff) = V2 \an \

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

105

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

2.6.1. Определение линейно-упругой среды

Рассмотрим наиболее простую, но вместе с тем самую широко используе­ мую в приложениях модель твердой среды с малыми деформациями, называ­ емую линейно-упругой средой.

Определение 2.6.1. Твердую среду с малыми деформациями, для которой оператор (2.1.39) свободной энергии Гельмгольца ф имеет вид

ф(е, в)

=

+

\ ( е - е)

■■4с ■■(е - е),

(2.6.1)

 

 

 

2Р

 

 

 

 

 

в

 

 

сЛв')

 

Фв

Фо +

Cv{e')de' -

в

 

в'

 

 

 

 

 

 

 

называют л и н е й н о - у п р у г о й

средой.

 

 

 

Здесь фо начальное значение свободной энергии, фо = е о —вщ\ cv(0) — теп­ лоемкость; 4С(0) — тензор четвертого ранга, называемый тензором модулей упругости; е{в) — тензор второго ранга, называемый тензором тепловой деформации. Все эти величины являются заданными.

Подставляя выражение (2.6.1) в (2.1.396) и (2.1.39в), находим определя­

ющие соотношения для линейно-термоупругого тела:

 

 

 

<х = 4С-- (е- 1),

(2.6.2)

v = т + ijidO' +

- а

1

(е-

(2.6.3)

• • (е

#0 в

Р

¥Р

 

 

 

 

 

е = ф + Ор

 

(4С —в Q) • • (в —в)+

 

Здесь

 

+

(1/ P ) OLQ • • 4 С • • (в —в).

(2.6.4)

 

в

 

 

 

 

 

 

 

о

а(в') d0',

e{9')d9'.

(2.6.5)

 

е =

Тензор ос называют тензором теплового расширения.

 

Если относительное

изменение температуры (9 — 9Q)/9Q не слишком ве­

лико, то для большинства твердых сред тензоры 4С и а

можно считать не

зависящими от температуры. В этом случае

 

а = «о = const, 4С<? = О,

(2.6.6)

(2.6. 11)

106

Глава 2. Упругие среды с малыми деформациями

 

 

е = а ( 0 - 0 о).

(2.6.7)

Модель линейно-термоупругой среды, для которой справедливы соотноше­ ния (2.6.6) и (2.6.7), называют линейно-термоупругой средой с не зависящи­ ми от температуры свойствами. Учет тепловой деформации в соотношении (2.6.7) называют законом Дюгамеля — Неймана.

Если рассматривают изотермические процессы с0 = $о>тов = Ои модель (2.6.1) называют моделью линейно-упругой среды.

Модель линейно-термоупругой среды предполагают всегда обратимой, т. е. всегда существует обратное соотношение к (2.6.2):

е = е + • •сг,

(2.6.8)

где тензор 4П, называемый тензором упругих податливостей, является

обратным к 4С:

(2.6.8а)

4С • • 4П = А .

С учетом соотношений (2.6.8) и (2.6.8а) свободную энергию ф (2.6.1)

можно выразить через тензор напряжений:

 

Ф = фв{0) + Л а • • 4П • •<Т’

(2.6.9а)

2Р

 

а с помощью (2.1.42) можно получить выражение для свободной энергии Гиббса в модели линейной термоупругой среды:

с = Фв(0) — (

сг • • 4П • • сг —\ сг ■■£.

(2.6.96)

2Р

Р

 

Очевидно, что если подставить (2.6.9.б) в (2.1.47), то после вычисления тензора производной е = —р(д(/дст) снова получим соотношение (2.6.8).

Таким образом, для линейно-упругой среды ^-модель (2.6.1) и ("-модель (2.6.96) являются эквивалентными и представляют собой различные формы одной и той же модели.

2.6.2. Число независимых компонент тензора модулей упругости

Тензор модулей упругости 4С, как всякий тензор четвертого ранга, имеет З4 = 81 компоненту Сг^ы в базисе сц. Однако вследствие наличия симметрии по первым двум индексам, а также по третьему и четвертому индексам, т. е.

(2 6 10)

число независимых компонент сокращается до 36, а условие того, что тензор 4С образует квадратичную форму в потенциале (2.6.1), означает наличие еще одного вида симметрии — по паре индексов:

Qljkl _ Qkllj

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

107

В результате максимально возможное число независимых компонент сокра­ щается до 21 (см. подробно т. 1, п. 4.6.9).

Для различных групп симметрии Gs число независимых компонент может быть еще меньше (см. т. 1, §4.6, а также п. 2.6.3). Число независимых компонент W ikl тензора упругих податливостей 4П равно числу независимых компонент Сг^ы соответствующего тензора модулей упругости 4С. Для тен­ зора 4П справедливы те же соотношения симметрии (см. т. 1, упр. 4 к § 4.6).

Тензор теплового расширения а является симметричным и имеет не более шести независимых компонент.

2.6.3. Тензорно-инвариантное и матричное представления модели линейно-термоупругой среды

Для линейно-термоупругой среды (2.6.1) можно использовать эквивалент­ ное тензорно-инвариантное представление (2.5.21). В этом случае потенци­ ал Ф задается в виде квадратичной функции от инвариантов I^s\ s ) (см. т. 1, п. 4.9.4).

Для ортотропной линейно-упругой среды

ф = р ф = р^фв +

2(е) +

+ 2Аб+а/^ .з(е)),

 

а=1

(2.6. 12)

 

 

где а ф /3 ф у ф а; а, /3, 7 = 1,2, 3; Аь ..., Ад — независимые упругие константы, число которых для ортотропной среды равно девяти (см. т. 1, п. 4.6.9).

Соответствующий тензор модулей упругости 4С для ортотропной среды можно представить в следующем виде (см. т. 1, п. 4.6.9):

з

4С = У^(Аас^ 0 с2а + А3+а(сг| 0 с Щ 0 е |) + А6+аОа 0 Оа) =

а=\

_

(g,

(g) cfc 0 СI- (2.6.13)

Для компонент Сг^к1 тензора модулей упругости 4С часто используют удобное матричное представление в виде матрицы 6 x 6 (см. т. 1, п. 4.6.9):

с 12

Схг

0

0

0

\

^22

С*23

0

0

0

 

(4С) = {Соф)

Сзз

0

0

0

 

 

2С-23

0

0

 

СИМ.

 

 

 

 

 

2С13

0

 

 

 

 

 

2С п )

108

Глава 2. Упругие среды с малыми деформациями

 

 

 

£1122

£1133

0

0

 

0

 

 

 

£2222

£2233

 

 

 

 

0

0

 

0

 

 

 

 

£3333

 

 

 

 

 

0

0

 

0

 

 

 

СИМ.

 

2£2323

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

2С 1313

0

 

 

 

 

 

 

 

 

2С 1212/

 

 

 

 

/А‘

А5

0

0

0

\

 

 

 

^2

а4

0

0

0

 

 

 

 

СИМ.

Аз

0

0

0

(2.6.13а)

 

 

 

 

2Xj

0

0

 

 

 

 

 

 

 

 

 

 

 

2А8

0

 

 

J

Тензор упругих податливостей 4П

ортотропной среды имеет точно такую

же структуру, как и тензор модулей упругости (см. т. 1, упр. 4 к §4.6):

з

 

4П = ^ (A 'QC2 ® с2а + \'3+а(Щ eg)

с х )с2р) + Х'6+аОа eg) Оа) =

а = \

= H^klCi <g)Cj <S)Cp <g)Ci, (2.6.14)

причем константы Xfa, естественно, отличаются от Аа и выражаются од­ нозначно через них по правилу обращения матриц. Чаще всего, особенно в приложениях, Xfa выражают через так называемые технические упругие константы:

Еа, ^ 7» G/37, а ф /З ф ^ ф а ;

а ,/?,7 =

1,2, 3,

(2.6.15)

следующим образом:

 

 

 

= 1 /Е а, Х'3+а = -иру/Ер,

2Х'6+а =

1/(2Сд7,

(2.6.16)

OL ф /3 ф 7 ф а; а,

/3, 7 = 1, 2, 3,

 

где Еа — модули упругости; z£ 7 — коэффициенты Пуассона; Сд7 — модули сдвига (всего технических констант — девять штук).

Тогда матричное представление

компонент IPJ^ тензора 4П в базисе сц

имеет следующий вид:

 

 

 

 

 

 

П12

П43

0

0

0

 

П22

П23

0

0

0

(4п ) = (1Щ )

сим.

Пзз

0

0

0

 

2П2з

0

0

 

 

 

 

 

 

2П13

0

2П12/

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

109

/ п 1111

д 1122

д 1133

0

0

0

\

 

 

д2222

д2233

0

0

0

 

 

 

 

дЗЗЗЗ

 

 

 

 

0

0

0

 

 

 

СИМ.

 

2П2323

0

0

 

 

 

 

 

 

2П1313

0

 

 

V

 

 

 

2П1212/

 

( \ / Е { - * w Е \ -V \z/E \

0

0

 

0

 

\ / е 2

—1/2з/ е 2

0

0

 

0

 

СИМ.

1/Яз

0

0

 

0

(2.6.17)

 

 

1/( 2G23)

0

 

0

 

 

 

 

 

 

 

 

1/(2G13

 

0

 

 

 

 

 

 

1/( 2а д /

 

Отметим, что иногда кроме v\2 ,

^23 также вводят взаимные коэффи­

циенты Пуассона U2 \, ^зь ^32 по формулам взаимности:

 

 

 

 

vp-y/Ep = v1p/E 1,

/3 ^7 -

 

 

(2.6.18)

Для трансверсально-изотропной линейно-упругой среды (группа сим­ метрии Gs = Т3) инвариантное представление упругого потенциала (2.6.1)

имеет следующий вид:

 

 

 

Ъ = р Фо+ ^(Ai/[3) \ е ) + A2/ f

2 (в)) + A3i; 3)( e ) i f (е) +

 

 

+ 2А41 ^ \ е ) +

(2.6.19)

А2 = А2 + Ai + 2 A3 + 2А5,

A3 = Ai + A3,

А4 = А4 +

А5,

где Ai, ..., А5 — независимые упругие константы, число которых для трансверсально-изотропной среды равно пяти (см. т. 1, п. 4.6.9). Соответству­ ющий тензор модулей упругости 4С можно представить в следующем виде (см. т. 1, п. 4.6.9):

4С —AiE (8) Е + А2С3 С) С3 + Аз(Е (8) С3 + С3 (8) Е)+

-\- A4(O I 0 Oj + О2 0 О2) ~\~2А5Д . (2.6.20)

Матричное представление компонент Сг^к1 тензора (2.6.20) имеет вид

 

/ с 1111 £41122

£41133

0

0

0

 

с 1111

£41133

0

0

0

(4с) = (а д

 

£43333

0

0

0

СИМ.

 

 

2С 1313

0

0

 

 

 

 

 

 

2С 1313

0

V

по

Глава 2. Упругие среды с малыми деформациями

 

 

 

(Х\ + 2А5

Х\

Аз

0

0

° \

 

 

\\

+ 2А5

Аз

0

0

0

 

 

 

 

А2

0

0

0

(2.6.20а)

 

 

СИМ.

 

2А4

0

0

 

 

 

 

 

 

 

 

 

2А4

0

 

2А5/

т. е. в этом случае

^2222 _ £ч1111

^1133 _ ^2233

^1313 _ ^2323

2 С Ш 2 = С 1111 — С П22

Тензор упругих податливостей 4П для трансверсально-изотропной среды имеет аналогичную структуру:

4П —Л^Е Е Т А2С3 (S) с2 Т А3(Е с3 Т с3 (S) Е)Т

+ A4 (O I ® Oi + 0 2 0 0 2) + 2 А5 А = n i j k l Ci 0 C j 0 с к 0 Cl, (2.6.21)

где пять констант А), ..., Ag обычно выражают через пять трансверсально­ изотропных технических констант:

 

Е', Е, v, U, G'

(2.6.22)

следующим образом:

 

А) = - v /E ,

Х\ + Аз + 2Ag + 2Ag = 1 /Е ', А) + Ag =

 

А) + 2X^ =

1 /Е , А4 + А5 = 1 /(2 G'), 2Х'Ъ= 1 /(2 G).

(2.6.23)

Технические константы (2.6.22) можно связать с техническими констан­

тами (2.6.15):

 

 

Е' = Е3,

Е = Е\ = Еъ

 

V = Щ Ъ

v ' = V 13 = V2 3 ,

 

G = G\2,

G' = Gn = G23,

(2.6.24)

их обычно называют следующим образом: Е' — модуль упругости в продоль­ ном направлении (направление оси трансверсальной изотропии с вектором с3); Е — модуль упругости в поперечном направлении (в плоскости транс­ версальной изотропии); v и G — коэффициент Пуассона и модуль сдвига в плоскости трансверсальной изотропии; z/ и G' — коэффициент Пуассона и модуль сдвига в продольном направлении.

Для трансверсально-изотропной среды G выражается через и и Е:

G =

Е

(2.6.25)

2(1+ 1/)*

Матричное представление компонент

тензора (2.6.21) в базисе

имеет следующий вид: