книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости
.pdf§ 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}} |
=\