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

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

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

§ 3.2. Замкнутые системы в материальном описании

281

Система (3.2.1), как и (3.1.1), замыкается после присоединения к ней

определяющих соотношений, состоящих из двух уравнений:

 

q = -А • V0, г] = -дф/дО,

(3.2.2)

и универсальных определяющих соотношений (3.1.3), (3.1.4), которые для тензора Пиолы — Кирхгофа имеют вид

Р = % ( F , в ), = ф {1^{C G (F)), в) = ф{¥, в ),

(3.2.3)

где тензорная функция

 

S ^ (F , в) = (р/р)F - 1• 5F G (F, в).

(3.2.4)

Система (3.2.1)—(3.2.3) содержит 16 скалярных неизвестных (плотность р

полагают выраженной через F):

 

0, u, v, F || Х \ £,

(3.2.5)

являющихся функциями лагранжевых координат Х г и времени t,

и состоит

из 16 скалярных уравнений (после подстановки (3.2.2) и (3.2.3) в (3.2.1)). Эту систему называют 6 U VF -системой уравнений термоупругости в ма­ териальном описании.

Для твердых сред чаще отдают предпочтение именно этой системе, по­

скольку ее область определения для искомых функций (3.2.5) известна —

о

это V х [0,£тах]. Исключение составляют задачи с фазовыми превращениями,

в которых V изменяется со временем и определяется в процессе решения. Однако и в этом случае система (3.2.1)—(3.2.3) чаще является более предпо­ чтительной, чем соответствующая система (3.1.1)—(3.1.6) в пространственном

описании. ( ,

( п )

Замечание 3.2.1. Хотя в определении (3.2.4) тензорной функции T°G присут­ ствует F -1, она, тем не менее, может быть представлена в виде зависимости от аргумента F. Для этого следует использовать представление F -1 в соб­

ственном базисе (см. т. 2, (1.3.35)),

а также

представление (3.1.4)

функции

(п)

 

 

 

 

 

 

 

 

Gi тогда получим

 

(п)

 

г

(п)

 

 

 

 

 

(3-2.6)

 

 

 

P = ^ ( F , 0 )

= y >

4E G° --

4<E G =

> )

 

3

О

Т'=1

/ 7

°

 

h ~ l Е

G

 

Еа/ЗРа ® Рр ® Ыорр >Ра + (1 ~ Ьс)Р/3®Ра),

G

р

 

а,(3=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р

G = А, В, С, D.

 

 

 

 

Е ар = ~

 

РА 0

Аа)

Вид тензоров 1 ^ совпадает с (3.1.9в). Тензоры 4Е£) называют лагранжевыми тензорами энергетической эквивалентности, они связывают тензор

282

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

 

 

 

( п )

 

 

 

 

 

напряжений Пиолы — Кирхгофа с Т ^:

 

 

 

 

 

 

P

, (п)

 

(п)

 

 

(3.2.7)

 

= 4E ^ - - T g .

 

 

Например, выбирая тензорные функции

(п)

(п)

для моделей А п в

= IFG ( C G ,0)

 

 

 

 

( п )

 

 

 

форме (3.1.8), получаем тензорную функцию PF°G в следующем виде:

 

Р

(п)

(п)

(п)

(п)

(п)

(3.2.8)

= 4Е° • • (4М --

С + £ C 2 + 6L ----

(X) С)),

(п)

(п)

 

 

 

 

 

 

где 4Е° = (р/р)F - 1• 4Е . □

 

 

 

 

 

 

Отметим, что в систему уравнений (3.2.1) входит тензор напряжений Пиолы — Кирхгофа Р; однако поскольку для заключительного (после реше­ ния системы (3.2.1) с граничными и начальными условиями) анализа поля напряжений в упругом теле обычно нужно знать тензор напряжений Коши Т, в материальном описании, кроме определяющих соотношений (3.2.6), исполь­ зуют также и соотношения (3.1.4).

Замечание 3.2.2. Достаточно сложный и, вообще говоря, не имеющий анали­ тического выражения вид тензорной функции (3.2.6) значительно упрощается для двух исключительных моделей Ау и By, поскольку

Р

= (р/р) F -1 • Т

= (р/р) F - 1

• (F • Т

• F T) = (р/р) Т • F T.

(3.2.9)

 

 

 

 

V

в

Подставляя выражение (3.1.4) для тензорной функции Т =

(3.2.9),

получаем

 

 

 

 

 

P

= :F ^ F ,0 )

= y ^

7l W - F T,

(3.2.10)

 

 

 

7—1

 

 

Г = к д ф /д 1ф), ф = ф(1ф(С), в ) ,

1\з )= д 1 ф /д С, С = (1/2)(FT• F —Е),

— тензорную функцию, имеющую явное аналитическое представление.

V V

Проведя в (3.2.10) замену С —►G, находим тензорную функцию F°B( ¥ ,6 ) для модели By.

Вычислять собственные базисы и собственные значения для моделей Ау и By при построении функций (3.2.10) не требуется. Указанные преимущества выделяют модели Ау и By среди других при использовании материального описания (также, как и модели А\ и В\ в пространственном описании (см. за­ мечание 3.1.3 в § 3.1)). Эти модели наиболее широко используют на практике при численном решении задач теории упругости с конечными деформациями.

Следующими по уровню сложности в материальном описании являются модели А\ и В\. Так, определяющие соотношения (3.2.10) в модели А\ заме­ няют выражениями

§ 3.2. Замкнутые системы в материальном описании

283

P

= ^ ( F , 0 ) = J 2 ^ F ~ l - F - 1T - I ^

• F _1,

(3.2.11)

 

7—1

 

 

^= °р (д ф /д р з)),

ф = ф (рв\А ) ,в ) , 1 ^ = д р 8р д А ,

Л = ^(Е —F -1 • F -1 т).

В этой модели, хотя можно и не вычислять собственные значения и соб­ ственные векторы, однако дополнительно по сравнению с моделями Ау и By, требуется обращать тензор F -1.

Наиболее сложными для численной реализации, как и в пространственном описании, являются модели Ац, В\\ и А\у и В\у.

Уравнение баланса энтропии в системе (3.2.1) можно преобразовать по­ добно тому, как это было сделано в п. 3.1.1, используя при этом основное термодинамическое тождество в материальном описании (т. 2, (3.3.19)):

о dib

о d6

_

°

т ~

n

7

+ '”' s

- p

- - v e v

=° -

(3'2' 12>

Дифференцируя это тождество по в, а затем используя второе соотношение в (3.2.2), получаем выражение для скорости изменения плотности энтропии в материальном описании:

о др _ о д2ф 80 дР

V ® v T.

(3.2.13)

Pdt = ~ PW ~dt ~ ~дв

 

 

Подставляя (3.2.2), (3.2.3) и (3.2.13) в уравнение баланса энтропии в си­ стеме (3.2.1), получаем искомое уравнение теплопроводности упругой среды

в материальном описании:

 

 

дв

»

н

(3.2.14)

рс£-

= V - (A - V 0) + (p/p)(4E G

Т $G) ■■V (8) V T + pqm,

где теплоемкость се определяется по (3.1.11а). Здесь учтено, что из (3.2.6) следует выражение для производной по в:

я р

Я (п)

,(п)

E

= Ср/ р/ ъ Ъ • • т А ,

(3.2.14a)

— =

— J=°r(F,0)

= 4 B°r

дв

дв Gy ’ '

G

7=1

 

 

 

( n )

 

 

 

где тензор

 

 

 

 

Т QG введен по формуле (т. 2, (3.8.147)):

 

 

 

т

= 2 ^ ж

10)

(3.2.146)

 

 

т «

1■4G-

 

 

 

 

7=1

 

 

3.2.2. 0U V - и 0U -системы уравнений термоупругости в материальном описании

Возвращаясь к общей системе (3.2.1)—(3.2.3), следует отметить, что гра­ диент деформации F можно исключить из этой системы с помощью соотно­ шения (т. 2, (1.2.10))

284 Глава 3. Упругие среды с конечными деформациями

 

F = Е + V <g>u T,

(3.2.15)

тогда получаем OUV-систему уравнений термоупругости в материальном

описании:

, ч

 

 

о B v

° Щ)

°

о

Р ^ = V - ^ ( ( E + V < g > u T) , 0 ) + p f,

= V • (Л • V0) + | ^ ( Е

+ V ® и т, в) . • V ® v T+ pqm,

 

du/dt = v,

(3.2.16)

состоящую из семи уравнений относительно семи скалярных неизвестных:

в, u, v || Х \ t.

(3.2.17)

Плотность р явным образом не входит в эту систему и всегда может быть вычислена из уравнения неразрывности (3.1.17):

р — °р det (Е + V ® и т).

(3.2.18)

Поскольку в материальном описании скорость v связана с вектором пере­ мещений и явным кинематическим соотношением, то можно исключить ско­ рость v из числа неизвестных, в результате получим OU-систему уравнений термоупругости в материальном описании:

 

о В^ЛХ

°

(п)

 

°

 

о

 

 

PW

= V

• ^ ( ( Е + V СХ) ит),

в) + pH,

 

о ВО

° ° °

в

( ( Е

°

 

°

я и

о

pce°g = V . (Л . V0) + | ^

+ v

® и т), 0) • . V ® ^

+ pqm, (3.2.19)

относительно четырех скалярных неизвестных:

 

 

 

 

 

 

0, и

||

Х \ t.

 

 

(3.2.20)

В частности, для исключительной модели Ау 07/-система (3.2.19) с учетом

формул (3.2.10)

и (т. 2, (1.2.11)) записывается в виде

 

 

 

оВ2и

°

о

 

 

 

 

р э ?

= v ■Р + pf•

 

 

<9(9

 

г

 

 

 

 

Bt = V

• (А •

VO) + ( ] Г ¥>7*1%) • • (F T• V

сх) V T) + pqm,

(3.2.21)

 

 

7—1

 

 

 

 

^ = ° » щ й у

дф

р = й «

' ' '

Ф = а Г ‘Н0 ),е),

 

7

 

7

7=1

 

 

 

I ^ d Z ^ / d C ,

С =

^ ( V ® u + V

® u T+ V

0 и • V

® u T), F =

Е + V ® и т.

§ 3.2. Замкнутые системы в материальном описании

285

3.2.3. T O U V F -система уравнений термоупругости в материальном описании

( п )

Если использовать соотношения (3.2.7) между тензорами Р и Т<з, то, добавив к системе (3.2.1) определяющие соотношения (т. 2, (3.8.166)) «в ско­ ростях», можно записать следующую T 6 UVF-систему уравнений термоупру­ гости в материальном описании:

 

о длг

°

 

,( п)

(п)

о

 

 

 

p ^

= V - ( 4Е& ..

Те ) + pi,

 

 

д в

 

 

 

,(п)

 

( п )

 

V 0 v T+ pqm,

рс£^ = V - (A . V 0 ) + (p/p)(4E b

 

Т вс)

 

du/dt = v,

5F /d t = V 0 v 1

 

(3.2.22)

( n )

»

 

 

 

 

(n)

(n)

 

( п ) вв

 

 

 

 

 

 

“ 5 Г -

(F . 0 )"

V ® v T + ZGh- T G + T G - ZGh+ T 0G —

 

 

 

 

 

 

 

 

d t

относительно 22 скалярных неизвестных:

 

 

 

 

 

 

( n )

 

u,

V, F

 

Х г,

t.

 

(3.2.23)

 

T Q , в,

 

 

Здесь тензор четвертого ранга

 

 

 

 

 

 

 

 

»

_

/4хэ (1243)

F

— 1 T )(1243)e

(3.2.24)

 

G h

=

(4P

Gh

 

 

 

Остальные обозначения такие же, как и в п. 3.1.3. При записи уравнения теплопроводности в системе (3.2.22) использовано соотношение (3.2.14).

Компонентное представление систем уравнений теории упругости в мате-

о

риальном описании чаще всего используют в базисе г* отсчетной конфигура­ ции К, (см. упр. 2 и 3 к § 3.2).

3.2.4. Система уравнений термоупругости для квазистатических процессов в материальном описании

Модель квазистатических процессов, введенная в п. 3.1.5, может быть рассмотрена и для материального описания. В этом случае система уравнений

(3.2.1)—(3.2.3) с учетом допущения (3.1.34) принимает вид

 

V • Р + pf = 0,

(3.2.25)

рСе^ = V ■Cx-ve) + °pqm,

(3.2.26)

( п )

О

 

Р = ^ G (E + V ® и, в).

(3.2.27)

286

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

Она состоит из четырех скалярных уравнений относительно четырех скаляр­ ных неизвестных функций:

0, и || Х \ t,

(3.2.28)

и называется квазистатической системой уравнений термоупругости в материальном описании.

Динамическое уравнение совместности и кинематическое соотношение в системе (3.2.1) при этом также удовлетворяются лишь приближенно и не рассматриваются. Точное выполнение этих уравнений обеспечивается для статических процессов, начиная с некоторого £Q-

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

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

Р = ( ° р / р ) ( АЫ • С) • F T, Т = F • (4М • • С) • F T,

или через градиент вектора перемещении:

Р = (р/р)(4М • • (Е + i v О и) • V О ит) • (Е + V ® и).

Упражнение 2. Используя результаты упр. 1 к § 3.1, показать, что для изотропных сред соотношения (3.2.10) имеют вид

 

Р = (Д е + ф 2с

+ V>3C2) • F T,

Т = JF • (У, + ф2С + ф3С2) • F T,

°i

° о

о

°

о

о

°

о о

о д ф

Т т

Л ~ n

Ф1 —

^ 1 + ^ 2 ^ 1

+ ^ 3 ^ 2 » —

^ 2 =

^ 2

+ Д ^ 3 ’

Фз=

Га

= P~Q]~ '

Д * = Д * ( С ) ,

OL= 1 , 2 , 3 .

Показать, что для линейной модели Ау изотропных сред эти соотношения принима­ ют вид

Р = l \ I \ (C)FT+ 212С • F t, Т = J { l \ I \ (C)F • F T+ 2/2F • С • F T),

или через градиент вектора перемещений:

P = / i ( V- u+^V( g) u - - V(g)uT)(E + V(g) u)+

+(8)u + V (8)u T+ V(8) u - V ( 8) u T) • (E + V ® u).

Упражнение 3. Используя результат упр. 1, показать, что компонентное представ­ ление определяющих соотношений квазилинейной модели Ау в базисе гi имеет вид

Р \ = (°р/р) M jkms e m s F l k g u ,

для модели Ау изотропных сред:

Р \ = (ф

+ ф 2р т °дк ° е т з + ф 3д * т °дк в д ', '’е т л е 81, ) Р 1к д и

§ 3.2. Замкнутые системы в материальном описании

287

Упражнение 4. Показать, что компонентное представление 9UVF-системы (3.2.1)- (3.2.3), (3.2.6) в материальном описании в базисе гц имеет вид

°p(dvl/dt) = WjPji + °p°f\

, Кде/dt) =

+ fc V ifr g v + °pqm,

d b p d t = V ju\ dul/dt = vl;

(n), , о

pi j = 1 е = ф - в(дф/дв),

 

Л = ффо(С°а Ф \ )). 9) = ф ф \ , о)-

^poij _

Щоijkl°j(s)

Х G —

1 jG k l

 

7=1

LjG k l

= д 1М'jG/д с ъ * .

G

= к д ф / д М ,

T ( s )

( \

(n )

okl \

+ (i -

 

74SJ /

 

h e y p f c 'g O ik O f l) ,

< ДЗ = hGG > { c% 1)

(n)

1

3

 

 

о о

c°rG¥ = — —

Y . K

 

l(hGQ \ Q i a + ( l - h G) Q \ Q i a) -

hp ( n - ПГ

Е°ры = E Ea0QiaQj0(hGQkl3Qla + (1 - hG)Qk0Qla).

OL,(3=\

Здесь pa = Q 7 rb pa = Q\pi.

Упражнение 5. Показать, что компонентное представление 6U-системы (3.2.21) для модели Ау в базисе гi имеет вид

'°p(d2uz/dt2) = VjPb+pf,,

°p(de/dt) = VjfaViO) + Р\\7Лдик/ д ф к + °pqm,

р \

=

t 0^{ d I (^ / d s jk)Flkgu.

 

 

 

 

7=1

 

 

 

< 7

= р ( д ф / д 1 ^ ,

ф = ф(1^\е^к),

в),

е = ф - в(дф/дв),

£ j k

=

( 1 / 2 ) ( V J IX/J

Х7k Uj -р V j4 X m

k Ui у

),

Flk = Si + V kum°g™\

 

 

где £jk = Cjk — компоненты тензора деформации С = е^тъ® rJ (см. т. 2, (1.2.1)).

Упражнение 6. Учитывая, что компоненты тензора Т в базисе гц совпадают с компонентами тензора Коши Т в базисе iy (см. т. 2, упр. 1 к § 3.2), показать, что для модели Ау компоненты Tli можно легко вычислить с помощью определяющих соотношений (3.1.3а):

д1(Л

7=1

§ 3.3. Постановки задач

289

/

[и]= О,

п• [Т] = О,

(3.3.2)

' т = о ,

кп • [q] = О.

Граница £ области V фактически также является поверхностью разрыва, поэтому на £ необходимо формулировать граничные условия, используя для этого указанные соотношения на скачках (3.3.1). Выделим еще несколько основных частных случаев граничных условий.

3. При наличии перехода материальных точек через часть поверхности £ i граничные условия на Хц для твердой среды имеют следующий вид:

 

М \ у \

+ п • [Т] + С2Е = 0,

 

 

М [ е

+ у ]

+ п • [Т • v q] + СЗЕ = 0,

(3.3.

3)

где [v] = ve —v и т. д.,

причем функции без индексов

соответствуют рас­

сматриваемой области V

и являются неизвестными, а функции с индексом

е

соответствуют внешности области V (внешняя среда) и считаются известны­

ми. Оставшиеся соотношения на скачках в системе (3.3.1) не используются при формулировке граничных условий — они накладывают ограничения на параметры внешней среды.

Если рассматривают модель квазистатических процессов в упругой среде,

то, согласно

определению

(3.1.34),

в граничных условиях (3.3.3)

полагают

v = 0 , в этом случае они принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.3.3а)

4. Если внешняя среда на части поверхности £ 2

является идеальной ж ид ­

костью, то для нее тензор

напряжений

Коши (Те = — р е )

является

шаровым

(см. т. 2, (3.8.168а)) и соотношения (3.3.3)

можно

записать в виде

(см. т. 3,

(1.2.416),

(1.2.41B), (1.2.41 г))

 

 

 

 

 

 

 

 

 

 

 

' M [ v n ] — Т п = р е — Спе ,

 

 

 

 

 

 

 

< M [ V T I ] - T

T I =

О,

I

=

1,2,

 

 

 

(3.3.4)

 

 

^ М [ в +

(t>^/2)] —T n v n

=

p e v n e —

С3 2

 

 

Здесь е,

Т п =

п Т n, v n

= п

- v ,

Т Т1

=

т /

Т п — неизвестные

функции,

относящиеся

к твердому

телу; ео,

р е ,

v T i e , v ne, qne

заданные

функции,

относящиеся

к жидкости;

C sЕf

определяется формулой

(т. 3, (1.2.40а))

/=1