Механика сплошной среды
..pdfследует, что ненапряженное анизотропное тело (а*/=0) может за
счет равномерного нагревания |
получать не |
только объемные, на |
|||||
и сдвиговые деформации. |
|
|
динамической |
термоупругости |
|||
Замкнутая система уравнений |
|||||||
анизотропного тела получается для и и Т из соотношений |
|||||||
доИ |
|
v |
р — -L и |
dp0s |
div q , |
||
1 рХ*= |
dt |
||||||
dxj |
г |
1 |
dt2 |
|
То |
|
|
причем последнее |
уравнение |
(теплопроводности) |
имеет вид |
||||
|
|
дТ |
а |
д2Т |
де; |
(15.34) |
|
|
Р°с |
ut |
|
uX.£uX.j |
-Т0оР$i f |
dt |
Граничные условия совпадают с условиями изотропного случая. Применение этих уравнений возможно как к собственно анизо тропным телам (кристаллам), так и к конструктивно анизотроп ным. В последнем случае сетка арматуры должна быть достаточ но густой, и все рассматриваемые величины (температура Ту по ток тепла q, деформации, напряжения) являются средними в не котором смысле. Понятия средних могут быть уточнены на осно ве опытов с образцами, в которых создаются «однородные» усло вия, или из теоретических соображений, которые специфичны для конкретных моделей тела. Конструктивно анизотропные тела на
зываются композитами.
§16. НЕЛИНЕЙНАЯ ТЕОРИЯ УПРУГОСТИ
Внелинейной теории упругости сохраняются все основные пред положения линейной теории упругости, за исключением предполо жения о малости деформаций; последние могут быть произволь ными.
Компоненты тензора деформации 6 (А) в лагранжевых коор динатах (хг=х*) выражаются через компоненты метрического тен
зора § (Д ) соотношениями
ъ = 1 + 2l, g t,= 6i/ + 2e!7) |
(16.1) |
dxL |
dxi |
Следовательно, как уже отмечалось, сам метрический тензор так же является тензором деформации, т. е. деформация тела вполне определяется его компонентами. Метрический тензор, как и тен
зор <§ определяется законом движения частицы х= х(х, t).
Инварианты тензора ^ относительно ортогональных преобра зований начальной системы координат х‘=х/ определяются как коэффициенты кубического уравнения |gij—Лб//|=0, их обозна чим здесь через а, Ь, с:
—h3 + ah2— bh + с= 0 ,
a = g tfi‘i\ b= 1/2 (а2—gijgij), c= g = (p0/p)2. |
(16.2) |
Уравнение \gu—h8ij\=0 получается при определении .главных де формаций и в обозначениях § 4 получаем h=gx,
а= 1йхУ b = Igt, с = Iёз, |
(16.3) |
Тело называется упругим, если все входящие в табл, б и в ос новное термодинамическое тождество (10.30) функции являются параметрами состояния, причем рассеяние w* равно нулю, так что функционал энтропии совпадает с энтропией s (10.20), r\=\s. Лю бая пара параметров таблицы (я, г, V) из реакции r(t) в момент t представляет вместе с Vf=V функции параметров процесса я (0 в этот же момент. Принимая, например, в качестве процесса (10.24), так что свободная энергия ф — функция деформации и температуры, получим уравнения состояния упругого тела в лагранжевых координатах (Л)
± S ‘I = |
Ро |
Ti |
g - . Ч»=Ч»(е,/. Т). |
(16.4) |
Р |
|
дТ |
|
Соотношения (16.4) можно записать и короче:
^ - 5 = ^ - , |
(16.5) |
|
Ро |
3 9 |
* > Т ’ |
если такие производные действительно суметь эффективно вычис лить. Это можно сделать просто для тела изотропного и по меха ническим свойствам, и по температурным, т. е. когда все термо механические свойства тела определяются скалярными константа ми и одним только единичным тензором /= (б/у).
Следовательно, для анизотропного тела соотношения (16.4) имеют вид
Emnpq, OCm |
Т), |
или |
(16.6) |
Ф = Ф ( | , ё , а , Г ),
где £, а — четырех- и двухиндексные тензоры константы упруго сти и теплового линейного расширения. Аргументами функции (16.6) будут различные свертки указанных в скобках тензоров; фактическое построение таких сверток для кристаллов основыва ется на свойствах симметрии кристаллической решетки. В конст руктивно анизотропных телах они зависят от свойств симметрии неоднородной структуры таких материалов, называемых компози тами. Наименьшее число алгебраически независимых между со бой сверток тензорных аргументов функции (16.6) свободной энер
гии или других скалярных функций состояния V1, каждая из ко торых представляет собой инвариант преобразования системы ко» ординат, допустимого термоупругой симметрией тела, образуют базис скалярных инвариантов 3 \, Sf^ •••, &п, так что (16.6) пре» образуется к виду
Ф =Ф (#1. |
&п; Т), tfh= # h{i, I , а, Т). |
(16.7) |
Среди этих инвариантов, конечно, инварианты любых ортогональ ных преобразований (16.3) и другие типа атпгтп, Emnpq(CLtnp&nq"Ь + anpQmq), но их не больше шести (/2<6); остальные зависимы*
так как <§ — симметричный 6-компонентный тензор. Из (16.7) находим
|
д |
+ JS _ dT= -?$- |
d l + |
or, |
|
||
|
|
|
|
|
д% |
дТ |
|
и из основного тождества |
p(d-^+r\dT)=SdS получаем уравнения |
||||||
состояния |
(сумма по k= \, |
2, ..., |
п): |
|
|
|
|
|
Р |
|
|
99 |
ч = — gL |
(16.8) |
|
|
|
|
|
дТ |
|
||
Тензоры |
(А=1, |
2, ..., |
/г<6) |
называются тензорным |
базисом: |
для напряжений рассматриваемого анизотропного тела.
Для изотропного упругого тела, которое будем теперь рассмат»
ривать, |
инварианты могут быть^ взяты |
в |
виде |
(16.2), |
свободная |
||
энергия |
(16.6) выражена через 9=1 + 26 |
и Т: |
|
|
|||
|
Л Т )=ф (а, |
6, с, |
Т). |
|
(16.9> |
||
Базис, получающийся дифференцированием |
инвариантов а, ft* |
||||||
с (16.2) |
по gij или по е//: |
|
|
|
|
|
|
|
°gij |
|
|
|
|
|
|
|
B($ = -^ -= a & ij- g lh |
В(2)= |
(Я!;>). |
(16.10) |
|||
|
Bfi)= - jZ - = g g 4 , |
g $ ~ l= (Blp). |
(16.11) |
||||
|
dgij |
|
|
|
|
|
|
Тензор напряжения в этом базисе представйм в виде |
|
||||||
|
S=2p% B{k), |
P = Y = ~ . |
|
(16.12) |
|||
|
(Ф« Ф*. Фс)= (Фг. ФгФз)= (Фй). |
k=\, 2, 3, |
|||||
|
|
или через ковариантные компоненты 5
5г/=2р% В ^ |
(16.13) |
где фа, фь, фс — некоторые скалярные функции, зависящие от ин вариантов а, Ь, с и температуры Г, определяемые физическими •свойствами твердого тела. Эти три функции считаем заданными на основании опытов:
г М а , 6, с, Г) ( * = 1 , 2 , 3 ) . |
(16 .14) |
Следовательно, (16.13) полностью определяет связь между тен зорами напряжений и деформаций. В действительности для обла сти больших деформаций упругих тел, таких как высокоэластич ные полимеры, эти функции изучены еще слабо, и их определение
лредставляет всегда сложную задачу для экспериментатора.
В случае, когда существуют потенциал и свободная энергия ф (16.9), экспериментально определяемые функции Ф* (16.14) долж ны удовлетворять условиям
Ъ = оу/г (*=1, 2, 3), |
(16.15) |
т. е. между ними должны существовать три соотношения, доказы вающие существование потенциала:
дфд _ |
дурь |
дура _ |
дурс |
дурь _ |
дурс |
,jg jgx |
db |
да |
дс |
да |
дс |
дЬ |
|
Они могут быть проверены в опытах.
Действительное приращение работы внутренних сил в изменяю щемся объеме V тела постоянной массы (dM=pdV) за время dt равно
Оно и при условии (16.15) не будет полным дифференциалом, как и приращение работы в объеме единичной массы
— SdS= ^ d£ij = ,tykcl<ifh>
Рдгц
поскольку полное приращение урза dt будет
dV= % dffh + - Q d.T=— SdH—sdT. dl p
Но в случае изотермического процесса (r=const) это будут пол ные дифференциалы
— Sd§= d\1) |
6'W=dWT, WT= [pydV + const. (16.17) |
p |
t |
Аналогично в случае адиабатического процесса-, при этом 8'Q=0, и потому s=const, вследствие чего
— Sd$ = d(q + sT)=du, |
8'W=8WS, Ws= J pudV + const. (16.18) |
P |
V |
Работа внутренних напряжений в единице массы и во всем объеме тела будет за время dt полным дифференциалом во всех тех слу чаях, когда в процессе деформации обеспечивается какое-нибудь соотношение между функциями состояния, приводимое к виду
f(s, Т)= О
т. е. дающее определенную связь между энтропией s и темпера турой Т, так как при этом
■j-Sdg = dy + s{T)dT=d(^> + |
J s d r ) , |
|
|
|
|
|
(16.19) |
8'W=dWh Wf= |
(* + jsdTJ dV + const. |
||
Если существует потенциал напряжений |
(16.5), |
то существует |
|
и потенциал деформаций. Обозначая (§ 9) * |
|
|
|
— S = o, |
— 5'-'=а'7(х, |
t), |
(16.20) |
РР
имеем из (10.30) и таблицы (я, г, V) при ri=s
8ф+ s8T= о1'18еи-= 8 (o‘iег/)—ег/8о‘>,
для термодинамического потенциала ф (я ) при п—(Т, о) получим
уравнение
6ф-НбТ+18ст=0, |
(16.21) |
||
откуда |
|
|
|
£ = |
----s= ----------- |
(16.22) |
|
|
до |
дТ |
|
Коэффициент теплоемкости cv при постоянных деформациях (6е/у= 0 ) определяется из условия 6'Q=pcv8T=pdu(T, §). Сле довательно, теплоемкость
* Отличие от компонент а*/(х, t) в Э мы при необходимости отмечаем аргу ментами (х, /).
__ |
д и (Т у У) |
У) |
(16.23) |
|
* — |
дТ |
дТ2 |
||
|
Аналогично можно получить ср=ди(Т, а)/дТ. Отсюда следует, что
|
д2ф < 0 , |
дгф (S, |
Т) |
0 |
|
(16.24) |
||
|
~дТ2 |
|
дТ2 |
|
^ |
|
|
|
Уравнения движения изотропного упругого тела в перемеще |
||||||||
ниях получаются |
подстановкой |
закона |
упругости |
(16.13) при за |
||||
данной функции |
ф(<?, Т) |
в уравнения |
(8.11). Система |
(8.11) |
при |
|||
заданной ф(#, Т) и заданном поле температуры Т(х, |
t) замкну |
|||||||
та относительно |
перемещения и или вектора х=х+и. |
|
|
|||||
Е сли поле тем пературы |
неизвестно, то |
испол ьзуется |
за ко н |
те п |
||||
лопровод ности, опред ел яю щ и й п о то к тепла |
q через |
гр а д и е н ты |
те м |
|||||
пературы , и уравнение б аланса эн тр о п и и |
|
|
|
|
|
|||
q‘ = —A£l дТ |
р Г — = — |
- £ г ( а " - % |
г ) - |
<16-25> |
||||
|
dxi |
dt |
А |
причем в изотропном случае A‘i=Xg‘i.
Граничные условия для основных задач в перемещениях или
напряжениях имеют вид (§ 12) |
|
||
|
х — Хг или S'vi= ^ )2 ) при Ф = 0, |
|
|
а для поля температуры |
|
|
|
|
Т ~ T z |
или q ‘v t = qs при Ф = 0 . |
|
Если функция ф вида |
(16.9) представима полиномом по 6, то |
||
напряжение |
не только |
представимо трехчленной |
формулой |
(16.12), но с |
помощью интерполяционного полинома |
Лагранжа |
(9.14) упругий потенциал представим явно через |
корни характе |
||||
ристического |
уравнения, определяющего |
главные |
значения |
g, |
|
метрического |
тензора .^ = /+ 2 § |
(4.55), |
которые выражаются |
че |
|
рез главные кратности деформаций (§4) |
|
|
|
||
|
h = V g i |
(1 = 1, 2, |
3). |
(16.26) |
Эти три инварианты взаимно однозначно связаны с инварианта ми (а, Ь, с) и потому сами для изотропных тел принимаются в ка честве аргументов потенциала я|г.
Ч>='|>(*1 К К Т). |
(16.27) |
Они удобны также в экспериментальных исследованиях. Истин ное главное напряжение о/ист, соответствующее деформации Хч
согласно термодинамическому тождеству p(d^ + sdT)= SdS, в ко тором
SdS = GTTdKh |
(16.28) |
непосредственно находящееся через ф
а?ст= р |
дф |
Ро |
дф |
(16.29) |
|
дкс |
|
|
|
и измерения Огист и кратностей Ki позволяют строить ф по опытным данным, а также проверять потенциальность напряжений:
= ~S~ < |
• |
(16-29') |
Задача построения ар (К, Т) на основании |
(9.14) |
упрощается, так |
как нужно построить функцию SF(ka, Т) только одного аргумента, после чего для тензора напряжений получится выражение
(8 — Ьр) |
(g — tff) |
S = F (i, Т )= |
(16.30) |
a,p,v |
(*.*-*$) |
|
Однако экспериментальные исследования по конечным дефор мациям твердых тел очень сложны: за длительное время установ лено немного потенциалов, главным образом для эластомеров, к которым в первую очередь относятся синтетические и натуральные каучуки (резины). Укажем несколько потенциалов, приведенных в [53]:
(1940): |
+&)(/g. - 3 ) |
+ ( l - b ) ( / g-l- 3 ) ] > |
(16.31') |
где / gi= g j + £2 + £з совпадает с а (16.3) |
в главных осях, / |
-1 то же |
|
с заменой |
gt на g~'(gi= lf); |
|
|
|
3 |
3 |
|
(1977): |
Ф ,= |* [(1 + Ь )2 (Я*-1) + ( 1 - р ) 2 ( ^ ‘-1 )1 ; |
(16.31") |
|
|
k—\ |
*=1 |
|
|
3 |
|
|
(1976): |
£ ( t f - 1) + -£ - |
(XS— 1)«] +Ф (Л), |
(16.31"') |
|
fc=i |
|
|
причем Л = ‘^ ^ = Я 1Х.2Яя. |
пренебрегают, |
полагая |
|
Сжимаемостью эластомеров часто |
А—1=А,|Я2Аз—1=0. При вычислении напряжений на основе потен циалов Ф (16.31) рассматривают функцию Ф'=Ф —рАДгА з, при чем в главных осях
(16.32)
дк<х |
\ дХа |
множитель Лагранжа р остается неопределенной функцией (х, /). Изотермическую сжимаемость учитывают, например, соотноше нием
1), - з p= a1 + o.2 + a3= S ‘igil. (16.33')
П
В общем случае соотношение Ми — Грюнайзена отражает свойст ва многих тел, включая металлы:
Р(А, Т ) - Р н(А )= ^Г (А )[и (А , Т ) - и и(А)], |
зз> |
Р0/Р= А , Г(А)=Г0— а1(1— А).
Здесь и — внутренняя энергия, РН(А), ип(А )— функции, опреде ляемые на фронтах ударных волн; при больших давлениях р фор мула (16.33) сохраняет достаточную точность.
Отметим, что при больших давлениях р (16.33') отношение максимального касательного напряжения тт а х = |о а—сгТ|/2 к дав лению становится малым:
> ах КР) ) |
->-0. |
(16.34) |
\Р /р-+°°
При больших давлениях термомеханические свойства твердых тел приближаются к свойствам идеальных жидкостей, как, впрочем, и свойства жидкостей вязких.
Заметим, что главные значения симметричного тензора £, за данного в декартовых координатах компонентами z*/, в главных осях имеют компоненты zk (6=1, 2, 3), выражающиеся через ин
варианты |
|
|
3z=Zifiih Z = V Z ifri, |
Izu I = det ztj |
|
с помощью угла вида напряженного состояния |
|
|
cos<p=(det ?j/)/z3 3*1^6 |
(16.35) |
|
формулами |
|
|
2ft—2=1/273 z cos |
|
|
rtii=0, m2= 2, |
m3=4. |
(16.36) |
Поэтому интерполяционный полином Лагранжа представляет со бой эффективное квадратичное представление произвольной функ ции тензора 8Г{2), удобное для экспериментального определения вида функций ЗГ ..) по простейшим опытам.
Глава V
СРЕДЫ СО СЛОЖНЫМИ СВОЙСТВАМИ
Механические свойства жидкостей и твердых тел, не обладаю щих совершенной упругостью и вязкостью, настолько переплета ются, что для тех и других нередко используются одни и те же со отношения между напряжениями и деформациями, и в этих слу чаях основные дифференциальные уравнения МСС для них совпа дают. Важный пример таких сред представляют полимерные ма териалы (смолы, каучук, ...). Технология их производства охва тывает область жидкого и твердого состояния, причем упругие и вязкие свойства являются существенными. Поведение металлов в технологических процессах и конструкциях в зависимости от диа пазона температур определяется вязкими, вязкопластическими, уп ругопластическими или упругими свойствами.
Установление связей между напряжениями и деформациями и замыкание системы уравнений производятся методами, изложен ными в гл. И. В ряде случаев поле температуры Т предполагает ся известным, и потому уравнения МСС становятся замкнутыми только на основании определяющих соотношений.
§ 17. ПРОСТЕЙШИЕ МАТЕМАТИЧЕСКИЕ МОДЕЛИ СРЕД СО СЛОЖНЫМИ СВОЙСТВАМИ
В этом параграфе рассмотрим формально близкие между собой модели сред со сложными свойствами; конечные деформации бу дем рассматривать только в ортогональных эйлеровых координа тах хц малые же деформации — в начально ортогональных лагранжевых координатах (х/). Компоненты тензора напряжений попрежнему в Э и Л будем обозначать о//, S*7, 5 Г/, скорости дефор
мации — Vij, Vij=Eijy деформации — Eih e2-j, девиаторы отмечать волной сверху
ви = о ц — оЬи, а = 1/3 <тг/6(7 = ~ |
tfigif, |
|
•■"'—И |
dvi |
dvj |
|
dxi |
|
Vi/=Vii— - divv6(/. |
|
|
dxj
Существенную роль будут играть также вторые инварианты девиа-
торов о,/, vij, причем квадратные корни из этих инвариантов бу дем называть модулями девиаторов и обозначать
a = V дери =
= у | - Y (ou —СГ22)2 + (a22—<T33)2 + (о3л— >ou)2 + 6 (o\2 + e223+ a2,),
v = V V iP tf = Y ~ Y Vu=~Y$ |
^22)*+ |
+ 6 (®i2 + .. .)• |
В излагаемых ниже теориях разработаны эффективные методы решения задач и вариационные принципы (§ 12).
1. Нелинейно-вязкие стабильные жидкости в простейшем слу чае отличаются от рассмотренной ранее (§ 14) классической жид кости тем, что коэффициенты вязкости зависят от тензора ско рости деформации и температуры. Для изотропной нелинейной вязкой несжимаемой жидкости, как и для классической, девиаторы напряжений и скорости деформаций пропорциональны:
оц=2\ша=2роц. (17.1)
Возводя левые и правые части первого равенства в квадрат, получим
v |
3 vu • |
( 1 7 ’ Г> |
и, следовательно,
°U=o8tf + -irL vif, так как vhh= 0. |
(17.2) |
зVu
В классическом случае ньютоновской несжимаемой жидкости ко эффициент (I при r=const постоянен. В рассматриваемом здесь случае коэффициент вязкости р есть некоторая функция инвари
антов v и сг однако такая, что а-^0 при 5->0:
а = Ф (vyсг Г), v\i (v)-+ 0 при о-*-0. |
(17.3) |
Эта функция находится из опытов на сдвиг и обычно не зависит от а, а только от 0 и Т Среднее напряжение о, как и в других несжимаемых средах, находится в процессе решения краевой за дачи.
Соотношение (17.1), справедливое как для классической, так и нелинейно-вязкой жидкости, можно трактовать как условие сов
падения направлений тензоров Oij и иг;- и потому оно называет ся векторным свойством среды.