
книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 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

