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

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

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

§ 4.4. Постановки задач в теории вязкоупругости

401

Зам ечание 4 .4 .3 . Как и постановки задач термоупругости в пространствен­ ном описании (см. п. 3.3.3), представленные выше постановки задач термо­ вязкоупругости являются сильно связанными, поскольку они не могут быть разделены на постановки задач теплопроводности и задачи вязкоупругости, даже если пренебречь «энтропийным фактором» связанности (т. е. членом

( п )

(ос/р) • • Т в уравнении теплопроводности (4.4.9)).

К перечисленным в замечании 3.3.4 шести факторам связанности, в за­ дачах термовязкоупругости добавляется еще один: функция рассеивания в уравнении баланса энтропии (4.4.1в) или в уравнении теплопроводности

(4.4.9), что является следствием неидеальности вязкоупругих сред. Вклад функции рассеивания в уравнение теплопроводности (4.4.9) во многих задачах оказывается весьма существенным, и пренебрегать им при неизотер­ мических процессах нельзя.

Эффект повышения температуры в вязкоупругих средах без подвода тепла к телу извне, а только за счет внутреннего тепловыделения при деформирова­ нии (обусловленного функцией рассеивания w*), называют диссипативным разогревом или саморазогревом тел (см. § 4.6). □

Обратим внимание также на шестой фактор связанности по классифи­ кации из замечания 3.3.4 к п. 3.3.3: для вязкоупругих сред зависимость определяющих соотношений (4.4.2в) от температуры можно подразделить на три составляющие:

1)зависимость от тепловой деформации е (4.1.66), когда применяется модель Дюгамеля — Неймана;

2 ) зависимость упругих свойств от температуры 0{t)\

3) зависимость «вязких» свойств, т. е. интегральной части соотношений (4.4.2в), от предыстории температуры вг(т).

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

пературы, чем упругие. Поскольку функция рассеивания

зависит именно

от вязких свойств, то она также явным образом зависит

от температуры;

в модели Ап с экспоненциальными ядрами (4.4.11) эта зависимость прояв­ ляется в виде функции ao(0(t)). Температурная зависимость w*(0) приводит к интенсификации саморазогрева вязкоупругих сред и может вызывать при определенных условиях эффект типа теплового взрыва (см. п. 4.6.9).

4 .4 .2 . Постановки динамических задач в материальном описании

Если использовать постановку динамической 0[ЛЛР-задачи термоупруго­ сти в материальном описании (см. пп. 3.2.1 и 3.3.4) и заменить участвующие в ней определяющие соотношения (3.2.3) на соотношения вязкоупругости в

402

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

форме (4.1.38), (4.3.16), то для моделей А п и Вп получим постановку ди­ намической eUVF-задачи теории термовязкоупругости в материальном описании, которая состоит из системы уравнений

 

р = °р det F -1,

p(dv/dt) = V • Р + р{,

 

 

p6(dr]/dt) = V • (Л • V0) + °pqm + w*,

(4.4.14)

о

1dF T/dt = V (8 ) v,

9u/9t = v

о

в области V х (0,£тах), а также определяющих соотношений в V х (0,£тах):

 

 

 

P

 

(п)

(п)

 

 

(4.4.15а)

 

 

 

= 4 E ^ - - T g ,

 

 

(п)

(n)

 

 

t

(п)

 

(п)

 

(4.4.156)

Т G = р(дф/дСа {€)) = J=G( C G(t), 9(t),

С %(т), в\т)),

 

 

rj

= дф/дв,

= —р5ф,

 

(4.4.15в)

 

 

 

г=0

 

 

 

 

 

 

4> ( с с ( Ц

0(i),

(c i ( r ) ,

04(т»,

G = A,B,

(4.4.15r)

 

r = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( n )

( n )

 

которые дополняют выражениями для тензоров 4 Е и и С^:

 

 

4 ( п ) 0

 

^

°

о

 

о

о

 

 

Е

=

2 ^

 

 

 

 

 

 

 

 

 

а,/5=1

 

 

 

 

 

 

 

 

 

ск=1

^

Ра -

^еЕ).

(4.4.16)

 

 

 

 

 

 

 

 

о

^СК5 Рек’ Рек II Е,

граничными условиями (3.3.10)—(3.3.17), которые при отсутствии фазовых превращений имеют следующий вид:

О

=

и

О О

 

О

 

 

е ц

...,

^

^

 

п • Р

t ne,

n - q

= qne

на Е ь

Е4,

Е7;

 

 

 

u = и е,

в = ве

 

о

о

 

 

(4.4.17)

 

 

 

на Е5,

Еб;

 

о

 

о

0,

о

о

=

0,

о

т ск =

°

 

и • п = 0,

v • п =

q

n

п - Р

0 на Eg,

 

и начальными условиями

 

 

 

 

 

 

 

 

 

 

t = 0:

v = VQ ,

 

U =

U Q ,

F = E,

6 = 6Q.

(4.4.18)

После подстановки определяющих соотношений (4.4.15) и (4.4.16) в (4.4.14) получим систему для 16 неизвестных функций — компонент следую­ щих векторов и тензоров:

в, u, v, F

Х \ £.

§ 4.4. Постановки задач в теории вязкоупругости

403

Плотность р, в силу уравнения неразрывности, обычно исключают из числа неизвестных функций.

Как и для

задач термоупругости в материальном

описании, задачу

(4.4.14)—(4.4.19)

о

что значительно

формулируют для известной области V,

облегчает ее решение.

Для конкретных моделей вязкоупругих сред формулы (4.4.156)—(4.4.15г) заменяют на какие-либо из соотношений, представленных в §§ 4.1-4.3.

Если

рассматривать модели

А п вязкоупругих

сред

с разностными яд­

рами и

моделью Дюгамеля —

Неймана (4.1.66),

то

плотность энтропии

г] выражается по формуле (4.1.69). Полагая д'ус^/дв зависящей только от температуры, уравнение баланса энтропии в системе (4.4.14) в этом случае можно переписать в форме уравнения теплопроводности вязкоупругой среды в материальном описании:

о дв ° Д ° лч

оп д (

(п)

о *

 

/АЛЛ г\\

Т \ О

,

рсе— = V • (Л • V0) -

рб»— рх "

—J + РЧт + w

(4.4.19)

причем вторым «энтропийным» членом в правой части уравнения, как прави­ ло, можно пренебречь по сравнению с w*.

В отличие от постановки задачи термоупругости в материальном опи­ сании, сформулированной в и. 3.3.4, задача термовязкоупругости (4.4.14)- (4.4.18) является сильно связанной даже при отсутствии фазовых превраще­ ний — ввиду наличия функции рассеивания w*. В общем случае неизотер­ мических процессов вклад функции в уравнение теплопроводности может быть весьма существенным и пренебрегать им нельзя (см. в и. 4.4.1).

Используя постановки динамических 0UV-, T0UVF- и 0?7-задач термо­ упругости в материальном описании (см. и. 3.3.4), можно сформулировать соответствующие динамические задачи термовязкоупругости. Так, постановка

динамической OU-задачи термовязкоупругости в материальном описании

состоит из системы уравнений (3.2.19):

од2м

°

_

о

о дО

° Д °

о

о *

(4.4.20)

р-^ 2 =

V

Р

+ pf,

рс£ —

= V • (Л • V0) + pqm + w

о

в области V х (0,£тах), в которые подставляют определяющие соотношения

( п )

( п )

кинематическое соот­

(4.4.15), выражения для тензоров 4 Е и и

(4.4.16),

ношение

и

(4.4.20а)

F = Е + V

а также граничные условия (4.4.17) и начальные данные (4.4.18). Эта задача рассматривается относительно четырех скалярных функций — компонент вектора перемещений и и температуры 0.

404

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

4.4.3. Постановки квазистатических задач теории вязкоупругости в пространственном описании

Формально постановки квазистатических задач теории вязкоупругости мо­ гут быть образованы из соответствующих постановок квазистатических задач теории упругости при конечных деформациях (см. п. 3.3.5) путем замены входящих в них определяющих соотношений упругости на какие-либо соот­ ношения вязкоупругости из моделей, представленных в §§ 4.1-4.3. Так, по­ становка связанной квазистатической задачи термовязкоупругости в про­ странственном описании для линейных механически-детерминированных моделей Ап термореологически простых сред имеет следующий вид:

 

V • Т + p f = 0,

 

 

 

 

 

 

 

 

 

(4.4.21)

 

pc£(de/dt) = V - ( \ - V e )

+ pqm + w*

в V\

 

 

( п )

*

 

• •

(п)

 

 

F U S ,

 

 

 

 

 

 

 

Т

J / 4R ( £ ' - T ')

с1Св(т)

B

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

< W* = -J (a ,/2 )JJd C ,(T 1) • • ^7 4R(2t'

-

 

 

 

• •

( п )

(4.4.22)

т[ -

т'2)

( Ю д(т

 

 

0 0

 

 

 

Ul

 

 

 

 

2),

 

(п)

в

_

_

 

 

(*.т)

 

 

^

^

 

( п )

 

(t’,т’) =

 

 

 

С 0 = С —£,

e = ja(0)de,

J

ад(6(т))(1т’,

 

 

 

 

во

 

 

 

 

о

 

 

 

 

 

 

 

,(п)

(п)

 

 

 

 

 

 

 

 

 

 

 

 

Т = 4Е • •

т ,

 

 

 

 

 

 

 

 

 

 

 

 

( n )

|

3

—III °

 

 

 

 

 

 

 

 

 

С =

-----

 

 

 

 

 

 

 

 

 

 

 

 

п - Ш —

 

Р“

 

 

 

 

 

 

 

(4.4.23)

 

,(п)

3

1а = \

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 т?. —

Е а/ЗРа <£>Р/3 ® Рр (8) р 0

 

 

 

 

 

 

Е =

Е

 

 

 

 

 

 

 

а ,(3 = \

 

 

 

 

 

 

 

 

 

 

 

 

 

kAa , Pa , Pa || F,

F = E —V ® u T B EU E ;

 

n ‘T | ECT= t ne,

u |Eu= u e,

n - q | E(j= g e,

 

в\Ев=ве,

 

< u

n = 0,

n • T • r Q = 0 на Eg,

 

 

 

 

 

 

(4.4.24)

Vt = 0 : в = 90.

Здесь (4.4.21) — система уравнений равновесия и теплопроводности; (4.4.22) — определяющие соотношения; (4.4.23) — система кинематических соотношений и соотношений энергетической эквивалентности; (4.4.24) — гра­ ничные и начальное условия, где Х а = Hi U Х2 U Х3 U Х4 U Х7, X w = Х5 U Не,

Tig = Hi U X2 U Хз U Х4 U Xg, Хб> = Х5 U Х7.

После подстановки соотношений (4.4.22) и (4.4.23) в (4.4.21) получим систему четырех скалярных уравнений относительно четырех скалярных

§ 4.4. Постановки задач в теории вязкоупругости

405

неизвестных — компонент вектора перемещений и температуры:

U, в || X, t.

(4.4.25)

Если рассматривать модель Ап с экспоненциальными ядрами, то опреде­ ляющие соотношения (4.4.22) заменяют на соотношения (4.4.11).

В частном случае могут рассматриваться изотермические процессы, при которых поле температуры в области V не изменяется: 0(х, t) = 6Q = const, тогда уравнение теплопроводности исключают из системы (4.4.21) и получают следующую постановку квазистатической задачи вязкоупругости в про­ странственном описании для линейных моделей Ап\

V - T + p{ = 0

в V,

(п)

Г .

(п)

T = J

4R(t — т) • • dC(r);

 

 

 

 

о

 

.(п)

(п)

 

 

 

Т = 4Е

• •

Т ,

 

 

 

с = ( 1/ ( п - 1П ) ) Е л г ш ра ® ра ,

^

3

а=1

 

(4.4.26)

.(п)

Еа/ЗРа

Р/3

0

1 7

Е =

X

Р/3

Р

а ,(3 = \

 

p - 1 = E —V ® u T

Aaa PaОРек

II F -1,

u . T|

i u -3 II

=

p

U = ue, n •1T • r a = 0 на

относительно трех компонент вектора перемещений и(х, £).

Если рассматривают линейные модели Ап и Вп несжимаемых изотропных сред (4.3.6) и (4.3.19), то квазистатическую постановку задачи вязкоупруго­ сти в пространственном описании записывают следующим образом:

 

V - T + pf = 0, det F - 1

= 1,

Т = —рЕ + (m +

 

(n)

(n)

(n)

r i{ t - T )

d /i(C G(r)))E + 4E ..2 Г2 Ч

 

 

G = A,B\

 

 

C G = (l/(n - Ш )) X ( А Г Ш -

~hG)Pa ® Pa .

 

 

a=l

 

 

< 4 (П )

Д

о

о

 

E =

X Ea/3pa

® pf3® p^ 0

pa,

 

O L , ( 3 = 1

 

 

 

,Аа ,р а ,р а II F -1,

F _ 1 = E —V ® u T;

( n )

r) d C G(r),

(4.4.27)

406

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

Гп-т|

= t ne, uL = ue,

H s 8 -n = 0’ n ' T ls8 - T“ = 0 -

Напомним, что согласно (3.1.4a): TIQ = 1, если G = А, и ho = 0, если G = B.

(n)(n)

Здесь также обозначен тензор: Е = 4Е • • Е.

Поскольку область V в пространственном описании является неизвестной, то для ее определения к системам (4.4.21)-(4.4.24), а также к (4.4.26) и (4.4.27) следует добавить соотношение (3.3.22) для определения неизвестной геометрии области V.

4.4.4. Постановки квазистатических задач для моделей вязкоупругих сред в материальном описании

Аналогичным образом получаем постановки квазистатических задач в ма­ териальном описании. Так, постановка связанной квазистатической задачи термовязкоупругости в материальном описании для линейных механиче­ ски детерминированных моделей Ап термореологически простых сред имеет вид

V • Р + pi = 0,

О дв

° /° ° лч

О

о *

(4.4.28)

рс£ —

= V • (А • V0) +

pqm + w

в V;

( п ) t

( п )

о

 

о

 

 

Р = 4 E ° J 4R ( £ / т ') • • <1Сд(т) в У US,

 

 

0

 

 

 

 

 

 

W* = -{ад/2 ) J J dCg{r\) ■

^ 74R (2 1' -

т[ - т£) • • сЮд(т2),

(4.4.29)

ооаг

О в _ _

(tr)

С - е . , в =

 

 

(п)

1

 

г—Illi

 

 

 

 

С =

п -Ш , Е

 

 

 

 

Ла" "Ф а «Р а-

 

 

 

(п)

3

Е,0.(3 1

Р/3 ® Р/з ® Ро

(4.4.30)

 

4 Е° := Е

Ап

 

 

 

 

 

а,(3=1

 

 

 

 

Аа , ра , ра || F, F = Е + V ® и т в V U Е;

 

п • Р

= tr;

uL

=

ue, —n - A - V 0 |o = qe, ОI о

= 0 ,

U о • n =

0,

n - P |o

-та = 0,

 

(4.4.31)

E8

 

 

E8

 

 

 

t = 0 :

0 =

OQ.

 

 

 

 

 

После подстановки определяющих соотношений (4.4.29) и кинематических соотношений (4.4.30) в (4.4.28) получим систему четырех скалярных уравне­

§ 4.4. Постановки задач в теории вязкоупругости

407

нии равновесия и теплопроводности с граничными и начальным условиями (4.4.31) относительно четырех скалярных неизвестных — компонент вектора перемещений и температуры:

и, в || Х \ t.

(4.4.32)

Если рассматривать изотермические процессы,

при которых 6 ( X \ t ) =

= const, то, в силу (4.4.28)-(4.4.31), получаем постановку квазистатической задачи вязкоупругости в материальном описании для линейных моделей Ап\

о

 

0,

р =

( п )

,

(п)

 

V - P + pf =

4е °

lR ( t - T )

• • dC(r);

 

( n )

 

1

X А2"ШРв 0

 

 

 

С =

 

Ро

 

 

n - l l l

 

 

 

а =\

 

 

 

 

4 Е 0 =

£

{Еар/ЮЬа 0 Р/3 0 P/J 0 Ро

(4.4.33)

 

 

 

а ,(3=1

 

 

 

 

 

Ра-Ра II F,

F = Е +

 

 

 

 

О

^ 1

О

 

 

 

 

П • Р о

= tп е , U о

 

 

 

<

1

'Ха

 

Sи

 

 

 

 

°

 

П • Р | о

 

 

 

U О • п = 0,

 

 

 

 

%

 

 

%

 

 

 

относительно трех компонент вектора перемещений u

Если рассматривают модели Вп несжимаемых изотропных сред (4.3.19), то определяющие соотношения в (4.4.33) заменяют на следующие:

t

 

 

(n)

t

(n)

( n )

( n ) .

 

 

 

>4

тр 0

r2(t -

r) rfG(r),

P = - PF~[ + {m+ r\{t T) d /i(G (r)))E U+ 2

4

 

о

 

 

 

 

0

det F =

1,

 

 

 

(4.4.34)

а задачу (4.4.33), (4.4.34) рассматривают относительно четырех неизвестных функций: и || X \ t .

4.4.5. Задача об одноосном растяжении бруса

В качестве примера рассмотрим классическую задачу об одноосном рас­ тяжении бруса, которая была подробно изучена в и. 3.4.2.

Брус полагаем вязкоупругим, изотропным и несжимаемым, а его определя­ ющие соотношения соответствуют простейшим линейным моделям Ап или Вп с экспоненциальными ядрами (см. (4.3.13), (4.3.25) и упр. 2 к §4.3). Общая постановка квазистатической задачи в пространственном описании имеет вид (4.4.27).

408

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

Закон движения бруса при растяжении не зависит от особенностей его механических свойств, он одинаков как для упругих, так и вязкоупругих сред, и определяется формулой (3.4.1). Следовательно, все кинематические

(п)(п)

характеристики: тензоры F, С и G — такие же, как и для упругих сред, и определяются в соответствии с результатами т. 2, упр. 1 к § 1.2, упр. 2 к § 1.3

и упр.

13 к § 3.2:

 

4

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

F = J 2 ka(t)e2a,

= 6 a,

Aa = ka,

( 4 .4 . 3 5 )

 

 

 

P a = P a

 

с Ч - Ш Ь г ш М.

(n)

3

(n)

 

 

 

d e V C =

 

 

( 4 .4 . 3 6 )

 

(п)

 

 

 

а=\

 

 

a=\

 

 

 

 

-

 

1

3

 

 

 

 

 

 

 

Г

п -

1

V Ьп~ш с2

с - 1= ( п - т ) ^ к Г ~ пё 1

 

^

 

 

III Щ

е«’

 

 

 

 

 

а=\

 

 

 

a=\

 

 

где обозначены функции от к\\

 

 

 

 

 

 

(Df

( и . ) -

 

2

 

fun - III

p m -n)/2 ,

( f _ ( ? _

 

(4.4.36a)

!

-

3 ( n - i n ) (fel

 

)’

1 2

/ з

2

1

ka{t) — коэффициенты кратности удлинения бруса по соответствующим коор­ динатным направлениям, причем из условия несжимаемости = 1/у/к\ .

Подставляя выражения (4.4.36) в дифференциальное уравнение определя­

ющих соотношений для

(см. упр. 2 к § 4.3), находим, что тензоры

имеют диагональный вид:

 

 

 

w w =

(4-4-37)

 

 

а =\

Функции Wa ^ одинаковы для моделей Ап и Вп, они удовлетворяют

следующим дифференциальным уравнениям:

 

 

d

 

 

 

1

(П)

/ ч

(4.4.37а)

 

 

dt

 

 

 

 

 

которые имеют решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

( п )

 

 

 

 

и 47) =

(

1 ~

т \

/

а (к [{т ))с1т

(4.4.376)

ехР (

- ^

)

----- ^

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

Подставляя выражения

(4.4.36)

в определяющие соотношения из упр. 2

к § 4.3, получаем, что энергетические тензоры напряжений

( п )

Т также имеют

диагональный вид:

 

 

Н

 

з

Н

_9

 

 

 

 

(4.4.38)

 

 

 

Т

 

Т ааеа,

а =\

§ 4.4. Постановки задач в теории вязкоупругости

409

(п)

 

 

 

m

( ^ - III + 2fciIII- n)/2 - 3 ) +

 

 

Т аа = - p k la ~ n + m + ^

 

 

 

 

 

 

 

i

 

(7 n - l l l

 

 

 

N

 

 

 

 

 

 

 

^ 2

 

1) -

^

 

(4.4.39)

 

 

 

 

 

н-------

 

 

 

 

— для моделей An и

 

 

n - I I I r a

 

 

7=1

 

 

 

 

 

 

 

 

 

 

 

 

Т аа = -pfc™"n + Ai(n - III)(1 + /3 + (1 -

/3 )(^ -ш + 2к(Н ~ п)/2-

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

~ К ~ Ш) ~

Wa )Bb)

(4.4.40)

— для моделей Bn.

 

 

 

 

 

 

 

^

7=i

 

 

 

 

 

 

 

 

 

 

 

 

 

Как было показано в п. 3.4.3, тензоры Т

и

Т

в данной

задаче

связаны

следующими соотношениями:

 

 

 

 

 

 

 

 

 

 

 

 

( п )

3

 

П(Т ё2

Т —

3

 

 

 

(4.4.41)

 

 

т

= £ С

о;=1

 

 

 

 

 

 

а = \

 

 

 

тТТ(п)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_ -уп—III

гр

 

 

 

 

 

(4.4.42)

 

 

 

 

 

'*

 

 

-1 ск*

 

 

 

 

Там же было

показано, что

при одноосном

растяжении

бруса

а\ ф 0,

а (Т2 = сгз = 0.

Тогда, подставляя

(4.4.39) и

(4.4.40)

в (4.4.42), получаем

следующую систему двух уравнений:

 

 

 

 

 

 

 

 

СГ\ = —р +

(га + —-Ц— (/сТ—111 + 2к\fclI I I _ n ) / 2

- з ) +

 

 

 

 

 

п —III

1

 

 

1

N

 

 

 

 

 

 

212

 

 

 

 

 

(7)

;

ип

 

 

 

+

{кпх~ш -

1) -

X

 

 

 

 

п - III

w \1]BO))kЧ

 

(4.4.43а)

0 = - р Щ т + ^ Ч Ш (* ;* Ш

 

 

7=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Щ~ш + 2/С|1 П _п )/2 —3)+

 

 

 

 

+

J

h

{k(m-n),2 _

1} _ g

w Ws ( 7 ))fc(ra-n)/2

 

 

п

— I I I

 

 

 

 

 

 

 

 

 

 

 

для моделей Ап;

 

 

 

 

 

 

 

 

 

 

 

 

ст1 р

(p(n -

Ш ) ( 1

+ /3 + 2 (1

(3)к™~ш )—

 

 

 

 

 

 

 

 

 

 

 

 

х

W ^ B ^ k ,711 -111

 

0 = - р + {р(п -

Ш)(1 + /3 + (1 -

7=1

 

 

 

п)/2)У

 

(4.4.436)

/З Щ ~ Ш+ к?11

 

 

 

 

 

 

 

 

 

N

w ^ ]B O ))kfll~n)/2

 

 

 

 

 

 

 

 

X

 

— для моделей Вп.

 

 

 

 

7=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исключая из этих систем р,

получаем соотношения между а\ и к\\

(п)

 

 

(п)

 

(п)

 

(п)

 

N

W\1]В™

(4.4.44а)

а\ = fhQ(k\) + 1\М{к\) + hN {k\) -

L {к\) ^

.

 

 

 

 

 

 

 

 

7=1

 

 

 

— для моделей

А п\

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

и

 

н

 

 

н

Щ

w

м

(4-4.446)

<71 =/x(l+/?)Z(fci) + /x(l -/?)tf(fci)-

L ( k i) J 2

Wi

B l)

 

 

 

 

 

 

 

 

7=1

 

 

 

— для моделей Вп, где обозначены функции от к\\

 

 

 

 

_

ттт

(П)

 

 

 

(n—III)/2

III—п\

 

п-Ш

(III—п)/2

Я

=

(n -

 

 

Q = к \

 

к1

III) (/с1

 

к1

)

 

( п )

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М =

п - Ш {кпх~ш + 2к?п~п)/2 -

3)(к»~ш -

к ^ - п]/2),

 

N

=

 

_ Х)кп-ш _

(fc(ni-n)/2 _ 1)fc(III-n)/2))

(4 .4 .4 5 )

L

= Щ~ш +

( l/ 2 ) / C | I II _ n ) / 2 ,

 

Z

= (п -

Ш ) ( А ; ? _ 1П -

/с[Ш _ п ) / 2 ).

Здесь учтено, что, согласно (4.4.36а)

и (4.4.41),

W^

= —(1/2)

 

При

= 0 из (4.4.446) следуют упругие соотношения (3.4.19а).

Кривые релаксации рассчитывают в несколько этапов.

 

 

А. Сначала вычисляют константы

 

и т ^\

 

 

 

Рассмотрим ступенчатое деформирование бруса, при котором коэффици­

ент кратности удлинения k\(t)

задается в виде

 

 

к\(£) = k®h{t).

(4.4.46)

Подставляя (4.4.46) в (4.4.376), получаем

 

w ^ \ t ) =

/ 1(£;?)( 1 - exp ( - t / r (7))).

(4.4.47)

С учетом (4.4.47) из (4.4.44а) и (4.4.446) находим выражение для Aa(t) — напряжения релаксации:

Aa(t) = CTJ (0) - <71 (t) = L (fc°) / ! (/c?)(r(0) - r(t))

(4.4.48)

— для всех моделей An и Bn, где r(t) — функция релаксации, которая в соответствии с (4.2.83) имеет экспоненциальный вид:

 

N

 

 

N

 

r(t) = г00 +

В ^

exp (—t/т ^), г(0) = г°° +

В ^ = 2^-

(4.4.49)

 

7=1

 

 

7=1

 

Уравнение (4.4.48), записанное в форме

 

 

 

 

 

N

 

 

 

(П)

}

= Е д(7)(1 ~ ехР Н

/ ^ ) ) ,

(4.4.50)

 

£ (*?)/i(*?)

7=1

 

 

можно использовать для вычисления констант В ^

и

 

Если известна экспериментальная кривая релаксации Асгэ(£), полученная при некотором фиксированном значении fcj, то с помощью соотношения (4.4.50) эту кривую можно аппроксимировать, выбирая параметры В^) и