книги из ГПНТБ / Кристенсен, Р. Введение в теорию вязкоупругости
.pdf§ 2.9. Цилиндр под действием внутреннего давления |
77 |
m
°zz(r,f) = — р о ' £ |
[p ln r (cij'fB (tZj) -f E (fl/)], |
“ |
Hm [Д (s)/(s — а/)] |
/= 1 |
|
|
(2.87) |
где комплексные корни а/ уравнения F ( s ) = 0 можно оп ределить на ЭВМ с помощью стандартной программы, когда в представлении для р(^) заданы характеристики какого-либо конкретного материала.
Квазистатическое поведение сжимаемого цилиндра
Для сжимаемого вязкоупругого цилиндра при осевой симметрии единственное нетривиальное уравнение квазистатического равновесия в цилиндрических координа тах г, Ѳи z имеет вид
+ |
_ L |
^ |
_ J L = o, |
(2.88) |
д г 2 |
г |
д г |
г |
’ |
где u = u(r, s ) — преобразование Лапласа от радиально
го перемещения. Решение уравнения (2.88) |
будет |
и = Cr D/r, |
(2.89) |
где C(s) и D(s) — функции переменной преобразования, подлежащие определению.
Используя зависимости между деформациями и пере мещениями и преобразование соотношений между на пряжениями и деформациями в виде интеграла релакса ции, получаем
(2.90)
(2.91)
(2.92)
где
е = (ди/дг) -f- (и/г). |
(2.93) |
78 |
Гл. 2. Изотермические краевые задачи |
Граничные условия для этой задачи те же, что и для предыдущей, и задаются зависимостями (2.78) с уче
том (2.79). Постоянные С и D определяются из преобра зования Лапласа от (2.78) и (2.79). Оказывается, что
|
|
o rr = |
C i — (b2 г2) С2, |
|
|
|
оѳѳ - |
C ,+ (b */r*)C 2, |
|
|
|
azz = |
2svCu |
|
где |
|
|
|
|
п |
— |
|
— Р(К —sp) |
|
^ 1 |
(К-— sp) + |
(b2/a2)[K (1 — 2sv) + |
sp] |
|
|
|
|||
г |
— |
p [/C (1 — 2 s v )+ s p ] |
|
|
и 2 |
(К — sp) + |
(b2/a2)[K (1 — 2sv) + |
sp] |
|
|
|
|||
причем |
|
|
|
|
|
|
~ |
2 * ( 1 — V *). ‘ |
|
(2.94)
(2.95)
(2.96)
(2.97)
(2.98)
(2.99)
Если ц(7) и k(t) принять по-прежнему в виде (2.68),
то р. и V будут отношениями полиномов относительно па раметра преобразования s. Примем, например, p(t) в ви де ступенчатой функции времени (2.85) и подставим ее
преобразование и преобразования p,(s) и v(s) в (2.94) — (2.99). Это приводит к преобразованиям напряжений, которые содержат отношения полиномов по s и могут быть обращены с помощью метода, изложенного в при ложении Б. Выполнив эту процедуру, считая вязкоупру гий коэффициент Пуассона действительной константой
и принимая sp(s) в форме (2.80), получим для агг вы ражение вида
Orr = Po[C(s) — (b2/r2)D (s )]/F (s ), |
(2.100) |
где
С (s) = — КВ (s) + А (s), D (s) = (1 — 2v) KB (s) + А (s),
F (s) = s {KB (s) — А (s) + (bW )[( 1 — 2v) KB (s) + А (s)]}.
Функцию F (s) можно разложить на множители соглас-
§ 2.10. Действие давления на границу сферической полости 79
но (2.86). Обратное преобразование Лапласа зависимо сти (2.100) дает
|
fC (о/) - (ЬѴг») D (aj) 1 еV |
|||
Огг (П 0 |
p°S |
lim [Д (s)/(s-— fly)] |
(2.101а) |
|
|
/=1 |
|
|
|
Другие две компоненты |
напряжений |
в |
соответствии |
|
с (2.95) и (2.96) |
оказываются равными |
|
fl;/ |
|
|
|
|
|
|
°ѳѳ 'Г>0 |
~ Po У ' |
[С (g /)+ (fca/r*) D (fl/)) e ; |
||
lim [F (s)/(s — a.j)] |
|
|||
|
/=1 |
s->ay |
|
(2.1016) |
|
|
а |
< |
|
|
|
|
||
|
|
2vC (a,j e |
1 |
|
<Дг (/■. 0 = Po 5 ] lim [F (s)/(s — aj)] |
|
l—‘ s->a-
Примеры вычисления напряжений по формулам (2.101) дали Кристенсен и Шрейнер [2.5]. В этой работе функ ция релаксации р(^) представлялась рядом (2.68) убы вающих экспонент и находились соответствующие корни F (s ). В работу [2.5] включен также один частный спо соб учета влияния границы, которая изменяется вслед ствие сгорания или коррозии материала.
Задаче о нагружении вязкоупругих цилиндров внут ренним давлением посвящены многочисленные работы. Один подход, типичное описание которого содержится в работе [2.16], заключается в сведении задачи к реше нию интегрального уравнения, или уравнений, что мож но легко провести численно с использованием непосред ственно измеренных механических свойств. Другое исследование по этому вопросу проведено в работе [2.34], где также учитывается возможное влияние корро зии границы цилиндра.
§ 2.10. Действие давления на границу сферической полости
В качестве последнего примера использования мето да интегральных преобразований рассмотрим действие давления на границу сферической полости в бесконечной
80 |
Гл. 2. Изотермические краевые задачи |
среде. В отличие от предыдущего примера эта задача используется для иллюстрации использования преобра зования Фурье. При этом учитываются инерционные чле ны; таким образом, рассматривается задача о динами ческом поведении, включающем эффект распространения волн. Этот эффект не будет, однако, рассматриваться здесь детально; основное внимание будет уделено мето ду решения задачи. Специальное исследование вопроса о распространении волн будет проведено в гл. 4. Реше ние, приведенное здесь, взято из работы Локкета [2.23].
В сферических координатах г, Ѳ, ср при условии сфе рической симметрии единственное нетривиальное урав нение движения упругой среды имеет вид
д2и |
. 2 ди |
2и _ |
р |
дЧі |
|
|
dr2 |
г дг |
г2 _ А + 2ц ä 2 ’ |
' |
' |
где и = и (г, /) — радиальная |
компонента перемещения, |
а X и ц — постоянные Ламе. |
Следуя методу, изложенно |
му в § 2.5, превратим это уравнение движения в преобра зование Фурье уравнения движения вязкоупругой среды
заменой и (г, t) |
на и (г, со), К на К* (іа) и ц на ц*(гсо). Это |
|||||
дает уравнение |
|
|
|
|
|
|
а* 5 |
, |
2 |
d“ |
2ц = —Й2со2«, |
|
(2.103) |
дг2 |
^ |
г |
дг |
г2 |
|
|
где |
|
|
|
|
|
|
Q2 = |
р/[Я*(г<о) + 2ц* (гео)]. |
|
(2.104) |
|||
Решение уравнения |
(2.103) |
имеет вид |
|
|
||
|
д |
' |
С (со) еШаг + D (со) e - iQ(ür |
' |
(2.105) |
|
и = ~я~ |
|
|
Г |
|
||
|
дг |
|
|
|
|
где функции С (со) и D (a) подлежат определению. По скольку. внешней границы не существует, не может быть и отражения, вызывающего приходящие волны; в силу этого могут существовать только исходящие волны и С((о) должно обращаться в нуль. Это приводит (2.105) к виду
« = — (Цт2)(1 ! tQcor) е~ іаслг- |
(2.106) |
Ненулевыми являются только компоненты напряжения огг и сгѳѳ= а фф. Используя преобразование соотношения
§ 2.10. Действие давления на границу сферической полости 81
между напряжением и деформацией и зависимость меж ду деформацией и перемещением, их можно получить из
(2.106) в виде |
|
|
|
Ъ„ - (Dir3) [4fx* (1 + iQ w ) — pcoV2J е~ійш |
(2.107) |
||
и |
|
|
|
■ Стѳѳ = — (Dir3) [2ц* (1 + a im ) |
+ |
Q2coV21*\ е~ аыг. |
(2.108) |
На границе полости г = а |
краевое условие имеет вид |
||
огг(г, а ) = р(ы ), |
г = а, |
(2.109) |
|
где р — преобразование Фурье |
заданного давления на |
поверхности полости. Используя (2.107) и (2.109), имеем
дЗ П |
|
D = 4ц* (1 + iQ m ) — рш2 а2 • |
(2.110) |
Получающиеся после подстановки (2.110) в (2.106) — (2.108) формулы удобны для обращения. Формулу обращения преобразования Фурье (1.73) можно перепи сать в иной форме
оо
f (t) = — |
Re Г / (м) еш dar, |
(2.111) |
ЗХ |
,} |
|
|
о |
|
где через Re обозначена действительная часть. Напри
мер, применяя (2.111) |
к обращению перемещения грани |
|||
цы полости и (a, t), |
находим его в виде |
|
||
и (a, t) = |
— |
R i' |
(1 і®ша) реш da |
(2. 112) |
|
я |
^ е .) |
4ц* ( 1 -f-ißcoa)— pco2a2 |
|
Интегрирование в (2.112) можно выполнить, если из вестно давление p(t) и величины ц* (ш) =ц '(ш ) -Нц'^®) и К* (too) = А/(<о) -[-ІА," (со) заданы как функции частоты ш для конкретного вязкоупругого материала. В общем случае провести это интегрирование аналитически невоз можно, однако практически можно проделать его чис ленно. В качестве примера зададим такой закон измене ния давления, для которого численное определение (2.112) упрощается.
6—851
82 Гл. 2. Изотермические краевые задачи
Примем давление на границу сферической полости
равным |
|
|
ап |
|
|
p ( t ) = - PЛ± J |
(О |
(2.113) |
и, |
|
|
Используя определение интегрального синуса |
|
|
U |
|
|
Si (и) = I* |
dx, |
(2.114) |
о |
|
|
можно переписать (2.113) в виде |
|
|
p {t) = (Ро/п) [Si((o20 |
— Si (cöi^)]. |
(2.115) |
Преобразование Фурье от этого выражения, определяе мое с помощью (1.72), дается формулами
|
О |
при |
0 < ; со < |
од, |
|
р(ю) |
Po/t© |
при |
ОД ДС О) < |
ОД, |
(2.116) |
|
О |
при |
о д ^ ® ^ 00- |
|
Получив p ( t ) из (2.115), можно найти интегральные си нусы по любым стандартным таблицам, например таб лицам Янке и Эмде [2.20]. Схематический график p (t) в зависимости от t дан на рис. 2.3, где од = 5, ©2= 1000.
Р и с . 2.3. Давление на границе сферической полости.
§ 2.10. Действие давления на границу сферической полости 83
Из него мы видим, что давление p (t), вычисленное по (2.115), является приближением к разрывному измене нию давления. Преимущество использования (2.115) вместо ступенчатой функции изменения давления состо ит в том, что характеристики, получаемые из (2.116), позволяют нам использовать интеграл обращения на ко нечном интервале частот вместо обычного бесконечного.
Теперь, используя выражение для р из (2.116), прида дим (2.112) вид
|
|
ш2 |
|
|
и (а, t) = |
Re |
(1 -|- t'Qcoa) ieim da» |
(2.117) |
|
4м|х* (1 4- ійсоа) — раз2 а2 |
||||
|
и,
Интегрирование в (2.117) производится в конечных пре делах и может быть осуществлено непосредственно, ког да p,*(t'co) и Х*(ш) заданы как функции частоты.
Вэтой динамической задаче, исследованной с помо щью преобразования Фурье, процесс обращения слиш ком труден, чтобы его можно было выполнить аналити чески. Однако в квазистатических задачах метод преоб разования Фурье может использоваться для получения аналитических решений. Процедура в таких квазистати ческих задачах совершенно аналогична уже описанной для преобразования Лапласа, где, как было показано, можно практически получать аналитические решения для механических свойств общего вида. Кроме того, ис пользование метода преобразования Фурье как в квази статических, так и в динамических задачах часто ведет
крешению, выраженному в конечной форме, требующей простого численного интегрирования.
Вобщем случае, когда мы ищем аналитические ре шения в переходных динамических вязкоупругих зада чах, используя методы интегральных преобразований, процесс обращения связан с точками ветвления. В дан ном примере такие точки ветвления вносит П в (1.117). В силу этой трудности аналитические решения удается получить редко. Более полно эта ситуация будет рас смотрена в гл. 4.
6*
84 |
Гл. 2. Изотермические краевые задачи |
§2.11. Свободные колебания
Взадаче, рассмотренной в предыдущем параграфе, задавались граничные условия, содержащие изменения во времени, и определялось динамическое поведение вяз коупругого тела. Другим типом задач, в которых учиты ваются инерционные члены, полностью определяющие
характер поведения материала, является |
класс |
задач |
о свободных колебаниях. Определением |
этого |
класса |
служит то, что работа всех сил на границе там, где они существуют, тождественно равна нулю; таким образом, граничные условия имеют вид
Ui==n |
На |
а«' |
(2.118) |
а ц П/ = 0 |
на |
Во. |
’ |
Кроме того, тождественно равными нулю должны счи таться внешние массовые силы, исключая, возможно, их временное приложение в самом начале движения. Та ким образом, эти задачи принадлежат к типу задач с на чальными значениями. Анализ поведения вязкоупругих тел, на которых наложены такие условия, и будет пред метом нашего изучения.
Преобразование Лапласа уравнений движения дает ся зависимостью
Щ,П + (1 — 2sv) 1 Ufji = - J — 2 |
д 2 щ |
(2.119) |
|
dt2 |
|||
Sfl (s) |
|
Мы будем искать решение уравнения (2.119), удовлетво ряющее граничным условиям (2.118) при некоторых на чальных условиях. Прежде всего рассмотрим частный случай, в котором коэффициент Пуассона будет считать ся постоянным, поскольку этот случай приводит к боль шим упрощениям.
Частный случай постоянного коэффициента Пуассона (ѵ = const)
Сначала напомним анализ свободных колебаний для упругих тел. Уравнения движения имеют вид
§ 2.11. Свободные колебания |
85 |
и по-прежнему применяются граничные условия (2.118). Решение ищется в форме
ul {xh t) = U f{x i) e M . |
(2.121) |
Функции координат Uf(x.) определяются таким обра зом, чтобы удовлетворялись уравнения движения (2.120), а удовлетворение граничным условиям (2.118) приводит к задаче о собственных значениях. Соответст вующие собственные значения и собственные векторы определяют собственные частоты и отвечающие им фор мы колебаний, которые обозначаются через
С0„ И U f ’n(Х;).
Спектр частот со» может быть дискретным или непре рывным в зависимости от того, конечны размеры упру гого тела или нет. Здесь и далее в этом параграфе мы будем рассматривать только случай дискретных собст венных значений; однако результаты легко можно обоб щить и на непрерывный случай. Полное решение при этом представляется суперпозицией собственных векто ров, и, согласно (2.121), это дает
= |
• |
(2-122) |
П
Произвольные постоянные, отвечающие каждому собст венному вектору, определяются таким образом, чтобы удовлетворялись начальные условия. С формальной сто роны процедура определения этих комплексных постоян ных использует ортогональность собственных векторов таким же образом, как это делается для определения постоянных при представлении функции комплексным рядом Фурье.
Рассмотрим теперь соответствующую вязкоупругую задачу. Предположим, что вязкоупругое решение име ет вид
’"(*,) П 0> |
(2.123) |
где fn(t) — неизвестные функции времени, а Uf-n (х{) ~
собственные векторы для соответствующей упругой за дачи, для которой коэффициент Пуассона таков же, как
86 |
Гл. 2. Изотермические краевые задачи |
|
|
и для вязкоупругой задачи. Подстановка (2.123) |
в урав |
||
нения движения (2.119) дает |
|
|
|
№ + (1 |
■ 2»Г'ир\ Г = |
|
|
|
= (р; sp) [s2 г - f n(0) - |
s f (0)j £/f |
(2.124) |
где |
|
|
|
/«(0) |
=/»(*)|<=о и /"(0) = |
гі/"(ОА«І*=о, |
(2.125) |
и величины /™(0) и /” (0) определяются из начальных условий. В то же время из соответствующего уравнения (2.120) для упругой задачи следует, что
U f f , + (1 - 2ѵ)—1U f:ü = { - со2 p/pj u f ' n. (2.126)
Подставим (2.126) в (2.124) и найдем соответствующие значения fn (s). Имеем
-п /"(0) + s/” (0)
(2.127)
sa + (cü^sp/p) '
Поскольку (о2п из соответствующей упругой задачи про
порционально р, решение, определяемое формулой (2.127) , в действительности не зависит от принятого уп ругого модуля сдвига р. Принимая вязкоупругую функ цию релаксации при сдвиге p(f) в виде ряда из убываю щих экспонент, как в (2.68), видим, что (2.127) имеет только простые полюсы, соответствующие которым корни легко найти, исполызуя стандартные программы для
ЭВМ, и обращение /(s) можно провести непосредствен но. В точности та же ситуация возникала и в задачах о квазистатическом поведении, где мы сочли полезным оставить в представлении (2.68) для р(^) столько чле нов, сколько требуется для достижения точности. Про стой пример, включающий обращение зависимости (2.127) , дан в конце этой главы в качестве задачи.
Полное решение задачи о свободных колебаниях име ет вид
М * ц* ) = 5 Х ,л(*і) П 0 |
(2.128) |
п |
|
при начальных условиях