![](/user_photo/_userpic.png)
книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf422 |
Глава 4. Вязкоупругие среды |
||
ленного времени |
и |
- — |
— rf — гЛ, следовательно, инте- |
дт\ |
1 |
д т \ д т 2 |
|
гралы по ( и t в формуле (4.5.16) |
можно вычислить независимо. |
Подставляя (4.5.17) и (4.5.21) в (4.5.13), получаем асимптотические раз
ложения уравнения теплопроводности:
р с е / д в |
д в {х) |
(4.5.23) |
t\ Vd t |
- а т ( ё - в е ) + w * {0) + х О ( 1). |
|
д £ |
|
|
4 .5 .7 . |
Осредненное уравнение теплопроводности |
|
Осредняя уравнение (4.5.23) по периоду колебаний согласно (4.5.9) и
учитывая периодичность функции в ^ (т. е. ( д в ^ / д £ ) = 0 ^^(1) — ^^(О ) = 0 ),
приходим к окончательному виду уравнения теплопроводности для темпера туры в:
|
р с £ ^ = - а Т ( в - в е ) + Гй*- t = 0: в = |
в 0 , |
(4.5.24) |
где |
= гс* — осредненная функция рассеивания. |
При выводе |
(4.5.24) |
сделана обратная замена безразмерного времени на размерное: t = t t \ .
Если собрать в уравнении (4.5.23) члены при более высоких степенях х
(при х , х 2 и т. д.), то получим уравнение для вычисления функций в^х\ в ^
и т. д. Однако вклад этих членов в значение истинной температуры в мал, согласно формуле (4.5.17), поэтому для рассматриваемых задач достаточно ограничиться только нулевым приближением для определения температуры в .
Осредняя (4.5.22), находим выражение для гс* — функции рассеивания,
средней за цикл колебаний:
W* = (го*(0)) = ae ^2(q(0)(f 2a{k\)) - ( / |
Q(4 ))2 (<?(0) - g(2F))). (4.5.25) |
||
а = \ |
|
|
|
Если брус упругий, то q ( t ) = 0, гс* = |
0 |
и задача |
теплопроводности при |
6Q = ве имеет элементарное решение: 6 (t) |
= |
во = ве = |
const, т. е. температура |
упругого бруса при циклическом деформировании не изменяется. Для вяз коупругого бруса гс* ^ 0 , что описывает появление теплового источника в
уравнении (4.5.24), |
поэтому дв/d i ^ 0, |
т. е. теплоизолированный брус всегда |
будет разогреваться |
(при во = ве и ат = |
0). Этот разогрев обусловлен только |
рассеиванием (д и с с и п а ц и е й ) энергии, поэтому его называют д и с с и п а т и в н ы м
р а з о г р е в о м , или с а м о р а з о г р е в о м .
4.5.8. Температура саморазогрева при симметричном цикле
Процесс деформирования называют с и м м е т р и ч н ы м ц и к л о м , если среднее
(п) |
(п) |
значение функции f |
\ (к \ ) за цикл колебаний равно нулю: ( / i(fci)) = 0 . |
![](/html/65386/197/html_PhmaX2xzBv.Kjra/htmlconvd-3PPgnP423x1.jpg)
![](/html/65386/197/html_PhmaX2xzBv.Kjra/htmlconvd-3PPgnP424x1.jpg)
|
|
§ 4.5. Диссипативный саморазогрев бруса |
425 |
||||
Если <20(0 ) = |
1, то возможен явный вид решения задачи теплопроводности |
||||||
(4.5.27) (упр. |
1 |
к §4 .5): |
|
|
|
|
|
|
|
в = во + g(0)/ 2pCe |
( l |
- exp ( |
- ^ |
) ) . |
(4.5.33) |
|
|
°LT |
V |
V |
рс£ 4/ |
|
|
В случае, |
когда <20(0 ) неограничена |
при 0 —►+ о о |
(такова экспоненциаль |
ная функция (4 .5.31)), реализуется температурный режим 5, при котором кривая саморазогрева имеет три участка: а — начальный, b — установивший
ся, где 0 |
рз const, и с — неустановившийся, на котором функция 0(7) выпукла |
||
вниз (см. рис. 4.5.3) и является неограниченной при t |
—►+ 0 0 . |
||
Если <20(0 ) |
зависит от температуры, но является |
ограниченной при 0 —► |
|
—►+ о о |
и при |
0 = в у имеет точку перегиба (такова функция (4 .5.30)), то |
реализуется температурный режим 6 саморазогрева (см. рис. 4 .5.3), который
характеризуется наличием четырех участков: начального а, установившегося b , неустановившегося с, за которым снова следует установившийся участок d .
В этом случае температура 0(£) остается ограниченной при t —►+ 0 0 .
Режим 7 имеет место при тех же условиях, что и режим 6 . Но в этом
случае на неустановившемся участке с температура саморазогрева 0 (£) дости гает такого критического значения 0&, что происходит тепловое разрушение материала (по аналогии с режимом 3), и участок d не успевает реализоваться.
4.5.11. Экспериментальные и расчетные данные
осаморазогреве вязкоупругих тел
Экспериментальные и расчетные кривые саморазогрева полиуретанового бруса при циклическом деформировании по гармоническому закону (4.5.3)
показаны на рис. 4.5.4.
Расчетное изменение температуры А 0 = 0 — 0О найдено путем конечно разностного решения задачи (4.5.24), (4.5.25) по неявной разностной схеме:
|
|
|
д ^ + 1 |
_ А0г + w*At/рс£ |
|
|
(4.5.34) |
|||
|
|
|
|
|
|
1 -\- OLT/ рс£ |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
где |
в г = в ( и ) |
— |
температура |
в |
г -м |
узле в момент |
времени |
tp, A t |
— |
шаг |
по времени. |
|
|
|
|
|
|
|
|
|
|
|
Значения |
ядра |
релаксации |
q ( 2 t f) |
и q ( 0), входящие в выражение |
(4.5.25) |
||||
для |
функции |
рассеивания, определяли по (4.5.5), а |
В ^ ) и |
в этой |
фор |
муле — по кривым релаксации методом, изложенным в п. 4.4.5. Значения констант В ( ч ) и для полиуретана приведены в табл. 4.4.1.
Функция а е { 9 ) аппроксимировалась выражением (4.5.30), а константы а\
и <22 в этой формуле приняты равными 21 и 208 К соответственно. Значения
![](/html/65386/197/html_PhmaX2xzBv.Kjra/htmlconvd-3PPgnP426x1.jpg)
§ 4.6. Модели вязкоупругих сред с малыми деформациями |
427 |
§ 4 . 6 . М о д е л и в я з к о у п р у г и х с р е д с м а л ы м и д е ф о р м а ц и я м и
4.6.1. Общее представление моделей вязкоупругих сред с малыми деформациями
Теорию вязкоупругих сред, изложенную в §§ 4 .1 -4 .4 , можно использовать и для моделей твердых сред с малыми деформациями, общ ее определение которых, применимое для любых моделей (идеальных и неидеальных), было
дано в п. 2.1.1. В § 2.1 изложен метод получения соотношений для моделей
тел с малыми деформациями из соответствующих соотношений для моделей с произвольными конечными деформациями. Этот метод можно применить и для вязкоупругих сред, заменив градиент деформации F по формуле (2.1.2) с учетом следствий А — Ж из этой формулы (см. п. 2.1.3).
Тогда, поскольку в модели малых деформаций все энергетические тензоры
( п )
деформации С совпадают с тензором малых деформаций е и все энергетиче-
(п)
ские тензоры напряжений Т совпадают между собой и обозначаются сг, то из
формул (4.1.61) и (4.1.64) получаем общий вид определяющих соотношений
для |
с т а б и л ь н ы х в я з к о у п р у г и х |
с р е д с |
м а л ы м и |
д е ф о р м а ц и я м и : |
|
|||||||||||
|
|
|
|
|
|
|
|
|
Т00, |
|
|
ОО |
|
|
|
|
а = £(е,в) = |
|
|
|
|
д(Р0 |
|
Е |
д р т |
ре d/r\ . . . d /Trip , |
|||||||
р{дф/де) = p Y i P l S ) 1^ (t) + |
|
■ |
||||||||||||||
|
|
|
|
|
|
|
|
' д р 3 |
|
|
д Л р |
|
|
|||
|
|
|
|
|
|
|
7—1 |
|
|
771=1 |
|
|
|
|
||
|
|
p f o |
E |
д р т |
7 |
7 |
|
* |
|
|
|
d p |
|
|||
V |
= |
" д в |
ШГ) |
' |
m- |
|
= |
- ? E |
|
&Т\ . . . CLTm , |
||||||
|
|
|
|
|
|
|||||||||||
|
|
|
777=1 |
|
|
|
|
|
|
|
|
777=1 о |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
Ф = (po(pa)(e(t)),9(t)) + |
|
|
Lpm d r 1 . . . d r m , |
|
(4.6.1) |
||||||||
|
|
|
|
|
|
|
|
|
m =i 0 |
|
0 |
|
|
|
|
|
Pm = Pm{t-T\, . . . , t |
- T m |
, 9 ( t ) , 9 { T \), |
. . . , |
9 ( т т |
), |
J ^ ( e ( 7 ) , |
e (r i), |
. . . , e(r m))), |
||||||||
где |
I j J = |
d i p |
1' / d e { t) \ |
J p j |
= |
d |
J p 1'/ d e ( t) - , |
J p ^ |
|
— |
совместные |
инварианты m |
||||
тензоров |
S (T \), |
. . . , s ( r m ) |
относительно рассматриваемой |
группы |
симметрии |
G s (7 = 1, . . . , 2 ^ 6 m ).
Для всех моделей А п определяющие соотношения (4.6.1) очевидно одина ковы, поэтому далее не будем различать эти модели.
Аналогичным образом из соотношений, представленных в пп. 4.2.1 и 4.2.2,
можно получить к в а д р а т и ч н у ю м о д е л ь т е р м о в я з к о у п р у г о й с р е д ы р а з н о с т
н о г о т и п а с м а л ы м и д е ф о р м а ц и я м и :
430 |
Глава 4. Вязкоупругие среды |
4.6.3. Постановки задач в теории линейной вязкоупругости
Сформулировать постановку задачи теории линейной вязкоупругости мож
но двумя способами:
1)применить общую постановку динамической 0?7-задачи (4.4.20) в ма териальном описании, проведя ее линеаризацию для случая малых
деформаций;
2 ) использовать постановку динамической задачи (2 .2 .2 ) для малых дефор
маций, которая была сформулирована для общего случая твердых сред (как идеальных, так и неидеальных).
В результате применения обоих этих способов получаем следующую по
становку с в я з а н н о й |
д и н а м и ч е с к о й |
з а д а ч и |
л и н е й н о й |
т е р м о в я з к о у п р у г о с т и : |
||||||
' °p(d2 u / d t 2) |
= V |
• (4R |
• • V |
<g> u) + |
p i ' |
в V , |
|
|
|
|
h e ^ = V • (Л • V 0 ) - е Г ( а . . |
+ °p qm + w * в V t |
|
||||||||
< n • (4R - - V |
0 u )E<T = |
t'n e , |
u |Eu = u e, |
|
|
(4.6.9) |
||||
- n • (Л • V 0 )E , |
= Qe, |
0\Е в = в ё , |
|
|
|
|
|
|||
t = 0: u = UQ, d u / d t = V Q , |
в = 6 Q . |
|
|
|
|
|||||
Здесь 4R и гс* определяют |
по формулам (4.6.8); |
f' |
и t nef |
— обобщенные |
||||||
массовые силы и поверхностные силы (см. (2.6.62)): |
|
|
|
|||||||
f ' = |
f - |
( l/p ) V • ф |
д ) , |
|
t'n e = |
t n e + |
n |
• f i d , |
|
|
|
|
д = в - в 0, |
fid = 4R • • ad. |
|
|
(4.6.10) |
||||
Н а б л а - о п е р а т о р |
V , как всегда в теории малых деформаций, вычисляется |
|||||||||
|
|
|
|
|
|
о |
|
|
|
о |
по отношению к отсчетной конфигурации /С, т. е. совпадает с V ; знак «°» при |
||||||||||
этом не ставится. |
|
|
|
|
|
|
|
|
|
|
Подобным образом формулируется постановка с в я з а н н о й |
к в а з и с т а т и ч е - |
|||||||||
с к о й з а д а ч и т е р м о в я з к о у п р у г о с т и : |
|
|
|
|
|
|
||||
' V |
• (4R • • V <g> u) + |
p f ’ |
= 0 |
в V , |
|
|
|
|||
h e ^ = V • (Л • V 0 ) + h m + W * В V , |
|
|||||||||
< п |
• (4R • -V |
0 u ) s CT= |
t'ne, |
u |Eu = |
u e, |
(4.6.11) |
||||
—n • (Л • V 9 ) s q = |
qe, |
0 \ ^q |
= e e , |
|
|
|
||||
t = 0: |
U = U Q, |
в = 6Q . |
|
|
|
|
||||
При изотермических |
процессах |
имеем |
постановку |
д и н а м и ч е с к о й з а д а ч и |
л и н е й н о й в я з к о у п р у г о с т и :