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

Механика сплошной среды

..pdf
Скачиваний:
33
Добавлен:
15.11.2022
Размер:
13.71 Mб
Скачать

следует, что ненапряженное анизотропное тело (а*/=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=— SdHsdT. 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 и иг;- и потому оно называет­ ся векторным свойством среды.