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

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

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

§ 4.5. Диссипативный саморазогрев бруса

421

4.5.6. Асимптотические разложения по малому параметру

Поскольку функция k \ ( t ) периодическая, то любые алгебраические функ-

(п)

ции от k 1, в частности, функции f а ( к i), определенные по (4.4.36а), будут

(п)(п)

также периодическими: / а = f а ( к i(£ )). Тогда температура 0, являющая­

ся решением уравнения (4.5.13), будет, вообще говоря, квазипериодической функцией.

Решение уравнения (4.5.13) ищем в виде асимптотического разложения

(ряда) по малому параметру н \

 

 

 

 

 

 

9

=

6 (F) +

О +

х 20 ( 1).

(4.5.17)

Первый член этого ряда в зависит только от медленного времени 7.

Согласно правилу (4.5.7), производную д в / d t от ряда (4.5.17)

вычисляем

следующим образом:

 

1( дб

дв ^ \

 

Л/1,

- 1оч

 

 

дв

 

 

 

Э1 = й \ Я + ~ д Г ) + х 0 ^ '

<4 -5 1 8 )

Функцию а о ( 6 )

и приведенное время t 1 и г '

после подстановки в них ряда

(4.5.17) также можно представить в виде асимптотического разложения:

а # (в ) =

CLQ( 6 )

+ я а ^

+

я 2 0 { 1),

i! = i!

+

+ я 2 0 ( 1).

(4.5.19)

Здесь, как и ранее,

0 (1 )

означает члены, сравнимые по порядку величины с

первыми членами ряда.

 

 

 

 

 

 

Функции а <9 = CLQ( 6 ) и

V определяются теми же формулами, что и исходные

функции а о и

если в них сделать замену в

—►в:

 

 

 

 

 

 

t

 

 

 

 

 

а е ( в )

= а е ( в ) ,

t f

а о (в ) d t .

(4.5.20)

 

 

 

 

 

о

 

 

 

Подставляя в формулу (4.5.16) разложения (4.5.17) и (4.5.19), а также учитывая правило (4.5.8) интегрирования квазипериодических функций, на­

ходим асимптотическое разложение для функции рассеивания ш*:

 

где

w * = w * ^

+

 

+ х 20 (1 ),

(4.5.21)

аЩ<г(0)'/;;№(?)) -C

 

 

f a(kt))(g(0) -,(?))+

 

° * ( ° )

f j k

, m

 

 

а =\

 

 

 

 

 

+

( / а ( к l)) 2 (^(0 ) - 2 ?(i') + q ( 2 t'

(4.5.22)

При

выводе этого выражения учтено,

что

подынтегральные функции пред-

 

 

 

 

(п)

 

ставляют собой произведения функций быстрого времени f а(кi(£)) и мед-

422

Глава 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 .

 

 

§ 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 К соответственно. Значения

§ 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,

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

н о г о т и п а с м а л ы м и д е ф о р м а ц и я м и :

428

Глава 4. Вязкоупругие среды

& =

w

I е

7=1

Ф = Р’о Р Г ^ в Ф ф Л ф ) )

 

 

t t

 

d r ) , W * = p

— ip2 dr\ d r,

 

9 t W

d T ~ P

 

 

о 0

rj = -(д'ро/дв) + (1 /р)а

a,

 

t

 

( p i ( t - T , j y ) ( e e ( t ) , e e ( T ) ) d T +

t t

 

 

 

+

(p2 ( t ~ T i , t -

r 2, A s\ e e ( T x ) , s e {T2) ) d T \ d T 2 , (4.6.2)

 

о 0

 

 

где

= d(po/dl!f'1', <^i7 =

dipi/dJ^',

9'ipo/90 производная ipo по второму

аргументу.

4.6.2. Модель линейно-вязкоупругих сред

Для линейных моделей вязкоупругих сред с малыми деформациями (их называют также моделями л и н е й н о - в я з к о у п р у г и х с р е д ) справедливы все соотношения, представленные в пп. 4 .2 .3 -4 .2 .1 2 , если формально провести

(п)

(п)

 

 

 

 

 

 

 

 

 

 

 

замену: Т ^ с г , С

—►£,

J

= 1.

В

частности,

представление

определяющих

соотношений

(4.2.59),

(4.2.60), (4.2.64)

и

(4.2.65)

в форме

Больцмана для

м е х а н и ч е с к и

д е т е р м и н и р о в а н н о й

м о д е л и

л и н е й н о - в я з к о у п р у г и х с р е д имеет

следующий вид:

 

 

 

 

 

t

 

 

 

 

 

 

 

<7

=

R

• • £0 =

4R (t

- т)

• • d e e { r ) ,

(4.6.3a)

 

 

 

t t

 

 

 

 

 

 

 

 

 

,

О ,

1

 

d e e (T\)

4R (2 1 -

 

 

• • d e $ (r2),

(4.6.36)

рФ = рФо + 2

 

T \ -

T 2)

 

 

 

о о

 

 

 

 

 

 

 

 

 

 

t

t

 

 

 

 

 

 

 

 

 

 

W

= " 2

d e e ( n )

■■ ^

4 R (2 t

-

т \ ~

т 2) -

■ d e e {T2),

(4.6.3B)

 

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е в =

e

- a

(6

- 6 0 ),

 

(4 .6 .Зг)

где 4R (t) — как

и ранее, тензор функций релаксации, обладающий свойства­

ми (4.2.70), (4.2

.71).

С учетом определения (4.2.62) тензора ядер релаксации соотношение (4.6.3а) можно представить в форме Вольтерры:

t

<т= АС - - е е - (t —т) • ■£в{т) dr.

(4.6.4)

§ 4.6. Модели вязкоупругих сред с малыми деформациями

429

Здесь, как всегда в теории малых деформаций, тензор модулей упругости 4М

(4.2.61) обозначен 4С (см. также § 2.6):

 

*R(0) = 4 С .

(4.6.5)

В соответствии с результатами п. 4.2.15, соотношения (4.6.3а) или (4.6.4)

можно обратить, представив их в виде обратных зависимостей:

 

 

t

 

 

 

ев = 4П • •

<т =

*П(£ — г)

• •

ё с г ( т )

(4.6.6)

— в форме Больцмана,

 

 

 

 

 

ев = 4П • •

о- +

*Г(£ — т)

• •

<т(т)

(4.6.7)

— в форме Вольтерры, где 4Г(£) — тензор ядер ползучести (4.2.101).

Модель линейно-вязкоупругих сред наиболее широко применяется для описания вязкоупругих свойств реальных твердых сред: полимеров, компози­

ционных материалов, биологических сред и др.

Если рассмотреть линейные модели А п термореологически простых вяз­ коупругих сред (см. упр. 7 к § 4.2) при малых деформациях, то придем к сле­

дующей модели т е р м о р е о л о г и ч е с к и

п р о с т о й

л и н е й н о - в я з к о у п р у г о й

с р е д ы :

 

 

t

 

 

 

 

 

t

 

 

 

сг = 4R - - e e =

l R ( t '

- T ') ■■(1ев ( т ) ,

ев =

4П Д — т ') ■■ d a ( r ) ,

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

t

t

 

 

 

 

 

 

 

ф =

ф 0 + ф ь ф \ =

\

йев{т\)

4R (2 1 '~

 

■ ■ d e e ( T 2),

(4.6.8)

 

 

 

2 P

 

 

 

 

 

 

 

 

 

 

 

о о

 

 

 

 

 

 

 

 

t t

 

 

 

 

 

 

 

 

 

 

a e

■47 4R ( 2 t , ~ т [ -

T '2) ••d £ d( T 2),

(t', r') =

а0(в(т))с1т,

w = —-

d £ e ( Ti)

 

0 0

 

 

 

 

 

 

 

 

 

 

где tr

приведенное время; т[ его значения при разных цц ао(в)

функция

температурно-временного сдвига (см. п. 4.1.12).

 

 

 

 

Используя формулу

(4.2.28а),

выражение

для

функции

рассеивания

(4.6.3в) можно представить также в следующей эквивалентной форме, кото­ рая справедлива и для механически недетерминированных моделей:

>д/ф\

d 0

(4.6.8а)

W = & • • £ — р —------ QL • • а -

.

d t

d t

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 .

 

 

 

 

При изотермических

процессах

имеем

постановку

д и н а м и ч е с к о й з а д а ч и

л и н е й н о й в я з к о у п р у г о с т и :