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

книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости

.pdf
Скачиваний:
38
Добавлен:
22.10.2023
Размер:
10.86 Mб
Скачать

§ 2.7. Лример состояния гармонических колебаний

67

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

ораспространении вязкоупругих волн.

§2.7. Пример стационарного состояния

гармонических колебаний

Составив систему основных уравнений теории изо­ термической вязкоупругости, можно перейти к решению конкретных краевых задач. В качестве первого примера рассмотрим стационарные гармонические крутильные колебания прямого кругового цилиндра. Соответствую­ щее уравнение движения упругого цилиндра в цилиндри­ ческих координатах г , Ѳ и г имеет вид

£

4 ,

 

 

 

 

 

 

/2 59ч

д г2

г

дг

д г2

г 2

р dt2

где «0= и ѳ

(г, 2,

/),

а р — массовая плотность. Положим

 

ие =

rO{z,t),

Ф (z,t) =

f(z )e iat,

(2.60)

где и — частота

гармонических колебаний. Комбинируя

(2.59) и (2.60), легко найти

 

 

 

 

 

f(z)

=

A sin {Q.zjh) +

В cos (£2z//i),

(2.61)

где

 

 

 

 

 

 

 

 

 

 

 

 

Q2 =

p©«А*/,*.

 

(2.62)

Из получающегося поля перемещений находим, что единственной ненулевой компонентой напряжения явля­ ется компонента сгѲг. Интегрируя авг по всему попереч­

ному сечению, определяем крутящий момент М в виде

М (г, t) = (я/2) (Ь* — а4) (рю2Л/й) ІА cos (Qz/h) —

*— Bsin(Qzjh)]ei(ai, (2.63)

5*

68

Гл. 2. Изотермические краевые задачи

где а

и b — внутренний и внешний

радиусы соответст­

венно,

h — длина цилиндра, а Л и

В —-произвольные

постоянные, которые должны быть определены из усло­

вий на концах.

 

 

 

 

 

 

 

 

Принимая,

что конец z = h закреплен и не допускает

перемещений,

тогда

как к концу 2 =

0 приложен крутя-

(\

tat

,

находим, что отношение крутящего

щии момент Me

 

момента на конце 2 =

0 к углу закручивания дается фор­

мулой

 

 

 

 

 

 

 

 

 

м

(0 ,0

 

 

М

_

я ( Ь * - о * )

2. ctgQ

(2.64)

Ф (0 ,0 - / ( 0 ) ~

2

Q

 

Следуя методу, изложенному в § 2.4, превратим уп­ ругое решение (2.64) в гармоническое вязкоупругое ре­ шение стационарного состояния, заменив в выражении (2.62) для П р на р*(іа>). Такая процедура дает

л

 

М __

п (64 — а4)

р

ctg Q*

 

(2.65)

 

По) ~

 

2

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

Q* (іа)

= [рсо%2/р* (ш )і,г.

 

 

Е сли известно р*(і'(о),

уравнение

(2.65)

дает

возмож­

ность определить угол закручивания /(0)

при 2 = 0 для

заданного значения

приложенного крутящего

момен-

л

согласно

(2.65), величина

л

та М. Заметим, что,

М рас­

сматривается как действительная,

а f(0)

как комплекс­

ное

число, которое определяется абсолютной величиной

Ф (0,

t) и фазовым

углом, на

который /(0) отстает от

Аі(0, t). В другом случае, когда рассматривается экспе­ риментальная установка анализируемого здесь типа, соотношение (2.65) позволяет найти комплексный мо­ дуль сдвига путем наблюдения за приложенным крутя­ щим моментом и угловым перемещением. Эта последняя процедура была проделана Готтенбергом и Кристенсе­ ном [2.9], и ее результаты будут подробно приведены в § 7.3.

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

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

69

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

Приведем теперь пример использования метода ин­

тегральных

преобразований для получения решений.

В задаче,

рассмотренной в § 2.7, исследовались ста­

ционарные крутильные колебания цилиндра и было по­ казано, как можно использовать этот результат для оп­ ределения комплексного модуля р* (іи ). При этом неяв­ но подразумевалось, что любые начавшиеся переходные процессы «отмирают», уступая место условиям стацио­ нарного процесса. Та же задача будет рассмотрена вновь, но теперь будут учтены начальные переходные процессы, связанные с вязкоупругостью материала. Пусть прямой круговой цилиндр с внутренним радиусом а, наружным радиусом b и длиной h считается закреп­ ленным на конце z = h, тогда как его конец 2 = 0 испыты­ вает закручивание с заданным углом, определяемым за­ висимостями

ф (0,і) =

0,

і < 0 ,

(2.66)

Ф (0, t) =

&sino)i,

0 < i < o o ,

 

где k — заданная амплитуда колебаний. Задача состоит в том, чтобы найти крутящий момент M(z, t) |г = о, требу­ емый для того, чтобы вызвать колебания (2.66). В этой задаче частота и считается достаточно малой, чтобы можно было пренебречь инерционными членами.

Соответствующее соотношение статической теории упругости между крутящим моментом М(0, t ) , действую­

щим на цилиндр, и углом

закручивания Ф (0, t) конца

2 = 0 имеет вид

 

 

ЛҢ0, t)

_

п (Ь* — ар р

Ф ( 0 , 0

~

2h

Следуя методу интегральных преобразований из § 2.5, это соотношение теории упругости можно рассматривать как преобразование Лапласа соответствующего вязко-

упругого соотношения, если заменить р на sp(s),

М (0, t) на М (s), а Ф (0, t) на Ф (0, s ), где р (s) — преоб­ разование Лапласа вязкоупругой функции релаксации

70 Гл. 2. Изотермические краевые задачи

\ x { t ) .

Это приводит к формуле

 

 

 

М (s) =

я Ф4 — «д) ^

Ф) kü>

(2.67)

 

К ’

2 h (s2 +

0)2)

 

где

использовано преобразование Лапласа

от Ф (0, і ) ,

определяемое условиями (2.66).

Чтобы обратить зави­

симость

(2.67), необходимо задать p(s). Предположим,

что рф)

задано в виде

 

 

 

N

 

V-(t) = G0 +

(2.68)

 

 

/=1

где Go, Gj и tj (/ =1, ..., N) — константы, которые нужно выбрать так, чтобы они представляли любую представ­ ляющую интерес частную функцию релаксации. Умно­ женное на s преобразование Лапласа формулы (2.68) после приведения к общему знаменателю дается выра­ жением

Sli ( s ) = n Л(5)------

,

(2.69)

П ^ -И ф 1) /=1

где A (s) — полином степени N относительно s, коэффи­ циенты которого определяются по значениям коэффици­ ентов Gj.

Подставляя (2.69) в (2.67), после обратного преоб­ разования Лапласа получаем

M(t) = я (Ь4 — a4) ka> X

2 h

е-Ш А 1(о)

X

N

2i'ü) П ( — ICD + tj *)

+ 2

+

 

Jat A (i(o)

 

 

 

2І(0 П (ко-{-tj

 

 

/=1

~

^ ( s )

( S + ^ M

 

h ( s + t 7 ')

-

/= i

j

e

— tu

n

\

1

. (2.70)

 

 

 

( ф г 2 +

* 2 )

 

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

71

По поводу метода обратного преобразования Лапласа см. приложение Б.

Первые два члена в (2.70) характеризуют стационар­ ный процесс, а последний — переходную часть процесса.

Готтенберг и Кристенсен [2.10] использовали полное решение (2.70) для предсказания начального переходно­ го поведения образца, который применялся для опреде­ ления комплексного модуля сдвига, как это описано в § 2.7. Функция релаксации при сдвиге для материала образца имеет вид, представленный на рис. 2.1, и опре­ делена по соответствующему комплексному модулю. Процедура установления связи между функцией релак­ сации и комплексным модулем обсуждается ниже

 

Р и с . 2.1.

Функция релаксации при сдвиге.

 

 

П олиуретановая

матрица,

содерж ащ ая кристаллы соли

и алюминиевый

поро­

шок,

при 26,1° С

по Готтенбергу и К ристенсену [2.9].

По

оси

абсц исс: десяти ч ­

ный

логарифм

времени

(в се к у н д а х ); по

оси

ординат:

в

ф унт/дюйм 2

 

 

 

(I ф унт/дюйм 2«0 ,0 7

кг/см 2).

 

 

 

 

в §

7.3. Функция

релаксации на рис. 2.1

в

пределах

точности чертежа

выражается

формулой

(2.68)

при

N = 8 ; значения Gj и tj приведены в табл. 2.1.

 

 

 

Сравнение результатов, предсказываемых формулой

(2.70) с учетом данных табл. 2.1,

с экспериментальными

результатами приведено на рис. 2.2, который взят из работы [2.10]. Соответствие кажется полным.

Таблица 2.1

Данные о функции релаксации

і

°/

ч

0

500

1,5-10—б

1

997

2

538

1,5-10—4

3

494

1 ,5 - ІО -3

4

392

1 .5 - 1 0 - 2

5

306

1 ,5 -1 0 - і

6

154

1,5

7

119

1,5-10

8

20

1 ,5 * ІО2

Р и с . 2.2. Пример переходного поведения при частоте 0,02 Гц и тем­ пературе 26,1° С (по Готтенбергу и Кристенсену [2.10]).

□ — моменты времени в сек у н д ах при

анализе

первого цикла,

О — моменты

времени в сек у н д ах

при анали зе

второго

цикла;

------------ эксперим ентальны е ре­

зул ьтаты . По оси абсц исс:

угол

закручивания

в

р ади ан ах; по

оси ординат:

крутящ ий

момент

в ф унт-дю йм U

ф унт-дю й м «0,01 к гм ).

§ 2.9. Цилиндр под действием внутреннего давления

73

Что касается приближения функции релаксации ря­ дом из убывающих экспонент, как это было проделано здесь, представляются важными следующие замечания. В общем случае требуется один член на каждый поря­ док величины в графических данных. Времена релакса­ ции, относящиеся к этим отдельным членам, должны располагаться в пределах одного порядка так, чтобы самое большое и самое малое время релаксации грубо соответствовали самому большому и самому малому временам, при которых моделировались данные о релак­ сации. Разумеется, приемлемо любое изменение этой процедуры, если аналитическое представление действи­ тельно удовлетворяет экспериментальным данным на требуемом интервале времени, однако изложенный при­ ем всегда можно использовать для окончательной про­ верки. Задача 2.10 в конце этой главы включает модель­ ную процедуру такого типа. Признавая частный харак­ тер этой процедуры и подобных ей, авторы работы [2.11] предложили систематический способ представления функций релаксации и ползучести.

§ 2.9. Цилиндр под действием внутреннего давления

Приведем другой более сложный пример применения преобразования Лапласа для получения вязкоупругих решений. Задача о действии внутреннего давления на вязкоупругий цилиндр, подкрепленный упругой оболоч­ кой, представляет практический интерес. Рассмотрим достаточно длинный цилиндр, для которого можно при­ нять гипотезу о плоском характере деформаций. Внача­ ле анализ проводится для несжимаемого вязкоупругого цилиндра с учетом инерционных членов. После этого со­ ответствующий анализ проводится для сжимаемого вяз­ коупругого цилиндра с пренебрежением инерционными членами.

Динамическое поведение несжимаемого цилиндра

Займемся задачей о действии внутреннего давления на несжимаемый вязкоупругий цилиндр в условиях пло­ ской деформации и осевой симметрии. При этих условия^

74

Гл. 2. Изотермические краевые задачи

зависимость радиальной компоненты перемещения от радиальной координаты находится непосредственно из условия несжимаемости. После того как распределение перемещений становится известным, задача сводится к виду, в котором она имеет только одну независимую переменную — время, т. е. является задачей с одной сте­ пенью свободы. При этих весьма частных условиях, свя­ занных с допущением о несжимаемости, можно учесть инерционные члены, не вызывая чрезмерных ослож­ нений.

Рассуждая вышеизложенным способом, запишем ус­ ловие несжимаемости в виде

ди дг + и'г = О,

(2.71)

где u = u(r, t) — радиальная компонента перемещения в цилиндрических координатах г, Ѳ и г. Решение уравне­ ния (2.71) имеет простой вид

и = С/г,

(2.72)

где C = C ( t ) — функция времени, которую нужно опре­ делить.

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

о „ =

2 \ і^ - +

S,

 

 

дг

 

 

аѲѳ =

2p“ +

S,

(2.73)

= S ,

где S — гидростатическое напряжение, которое является основной неизвестной задачи. Обратим эти зависимости в преобразования Лапласа вязкоупругих соотношений между напряжениями и деформациями:

Ъ„ {г, s) =

2sjx (s) ди

s)

+ s (r, s),

 

0ѲѲ(r, s) =

2щ (s)

+

S (r, s),

(2.74)

öz2 (r, s) =

S (r , s).

 

 

 

§ 2.9. Цилиндр под действием внутреннего давления

75

Единственное нетривиальное уравнение движения име­ ет вид

до гг

I

a tr ~ аѲѲ _

сри_

(2.75)

дг

^

г

Р dt*

 

Подвергая это уравнение преобразованию Лапласа при

условии покоя в начальном состоянии и используя

(2.72)

и (2.74), находим решение

 

5 = ps2C ln г + D,

(2.76)

где D = D (s).

Используя (2.72) и (2.76), можно найти преобразо­

вания напряжений в виде

 

о п (г, s) = — 2sp (C/r2) + ps2C ln г +

D,

оѳѳ (г, s) =

2sp (C/r2) + ps2C \nr + D,

(2.77)

azz =

ps2C ln r Ц-D.

 

Граничные условия задачи принимаются в форме

a f f ( r ,t ) = — p(t)

при

г

=

а,

 

0rr( r , t ) = — q(t)

при

г

— b,

(2.78)

и (г, t) = и0 (t)

при

г

=

Ь,

 

где q(t) и ис(і) — давление и перемещение ограничива­ ющей цилиндр тонкой упругой оболочки, которые свя­ заны зависимостью

ис (0 = q (i) [b2 (1 — v2c)/Ech] ;

(2.79)

здесь Vc и Ec — упругие константы оболочки, h — толщи­ на оболочки, а b — ее радиус. Последние два уравнения (2.78) представляют собой условия непрерывности на границе между вязкоупругим цилиндром и упругой обо­ лочкой. _

Величины C(s) и D(s) находятся из граничных усло­ вий (2.78) с учетом (2.79). Используя то же представле­ ние функции релаксации при сдвиге, которое дается ра­

венством (2.68), можно записать sp, в виде

sp, = A (s)/B (s),

(2.80)

76

 

Гл. 2. Изотермические краевые задачи

 

 

где А (s)

и

B ( s ) — полиномы

относительно

s.

Согласно

(2.80), C(s)

и D(s) выражаются формулами

 

 

 

 

C(s)

= —sp B (s)/F (s),

 

 

(2.81)

 

 

D(s)

= - s t E ( s ) / F ( s ) ,

 

 

(2.82)

где

 

 

 

 

 

 

 

 

F ( s ) = — 2sA ( s ) ( - --------Ц -I- ps2ß (s) ln —-----------^

 

-----,

'

 

Ч а *

b2 ) ' "

w

b

b ^ l - v ^

 

 

 

 

 

 

 

 

(2.83)

 

£ (s) =

ps2ß (s )ln ö ------ £c hR (5)-

.

(2.84)

 

 

&2

 

 

63 ( l - ^ )

 

 

Взяв, например, давление в виде

 

 

 

 

 

 

 

p(t) = p0h(t)

 

 

 

(2.85)

и подставив преобразование Лапласа от этого выраже­ ния в (2.81) и (2.82), получаем полное решение для C(s)

D (s). Знаменатель F(s) в выражениях для

C (s ) и

D(s) можно разбить на множители

 

F (s) = x(s —■ fl;) (s — a2)...(s — am),

(2.86)

где öj — корни, в общем случае комплексные, х — посто­ янная, а целое число пг определяется количеством чле­ нов, включенных в представление типа (2.68) для р(0- После того как проведено обратное преобразование ра­ венств (2.77) методом, изложенным в приложении Б, на­ пряжения выразятся формулами

еа/

 

 

 

X

 

°rr (r,t) = — Ро £ lim [Д (s)/(s — aj)}

 

 

/=1 s->a

 

 

 

 

X ------- A(üj) +

P ln r {cLjf В (fl/) +

E (fl/)

°ѳѳ ('»О = — PoX

 

 

 

 

eV

■ А (а,) +

p ln г (а/)2 В (а/) +

E (fl/)

X

£nlim [F (s)/(s — af}}

=\