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

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

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

$ 3.6. Термореологически простые материалы

147

случае отсутствует принцип соответствия или зависи­ мость между упругими и вязкоупругими решениями краевых задач. Однако есть два частных случая, кото­ рые заслуживают внимания. Это случай, когда задан­ ное температурное поле обладает пространственной од­ нородностью, но зависит от времени, и случай, когда температурное поле зависит от координат, но постоянно во времени. В последнем случае соответствующая крае­ вая задача подобна краевой задаче при изотермических условиях, но механические свойства при этом зависят от координат. Такой тип задач относится, следователь­ но, к задачам для неоднородной среды, все еще доволь­ но сложным. В другом случае, когда температурные по­ ля зависят от времени, но не зависят от координат,, второй член в уравнении равновесия (3.95) тождествен­ но равен нулю и можно применять методы интеграль­ ных преобразований по приведенной переменной време­ ни подобно тому, как это делалось в изотермическом случае.

Пример

Приведем пример, который иллюстрирует использо­ вание зависимостей (3.93) и (3.94). Этот пример, так же как и некоторые из приведенных ранее формул, при­ надлежит Муки и Стернбергу [3.16]. Рассматривается задача о бесконечной изотропной вязкоупругой плите тол­ щиной 2а, подвергнутой нагреванию, которое приводит к известным распределениям температуры по толщине. Используются декартовы координаты, в которых свобод­ ные от напряжений грани плиты имеют уравнения z — = ± а . Распределение температуры T(z, t) считается из­ вестным либо из экспериментальных наблюдений, либо из теоретического анализа. Задача состоит в отыскании получающегося в результате напряженного состояния.

Поле перемещений принимается в виде

 

и х = и у = 0, U z — U z { z , t) .

(3.96)

Все получающиеся в результате деформации равны ну­

лю, кроме еzz— dUzldz, откуда следует,

что

Shk = 6z z i ехх == вуу =■

Y36ZZ,

вгг — 2/зеи, (3.97)

10*

148

Гл. 3. Термовязкоупругость

а остальные компоненты деформации равны нулю. Не­ нулевые компоненты напряжения соответственно равны

= а,уу»

syy

Ѵз(^zz

&хх)>(3.98)

}kk = 2

S z z = 2lз(<7г

®xx)-

 

Девиаторное соотношение между напряжениями и деформациями (3.93) для компонент с индексами zz в силу (3.97) и (3.98) принимает вид

(Z, I) - охх (z, S) = \ G, ( l - l ' ) (dlkk (z, S')/öS') d t . (3.99)

о

Подобным же образом уравнение (3.94) с учетом (3.97)

и(3.98) принимает вид

£« (*. і) + 2oxx(z,l) =

= j G3 (S - S') j - r h kk (z, t ) - a (z, t ) W .

(3.100)

о

Единственным нетривиальным уравнением равновесия оказывается уравнение

dozz(z, t)/dz = 0.

(3.101)

Граничные условия таковы: .

Ozz{z, t) = 0 при г =

(3.102)

Согласно (3.101) и (3.102), требуется, чтобы выполня­ лось условие

<Jzz(z, t) = Gzz{z, I) = 0.

(3.103)

Применим теперь к (3.99) и (3.100) преобразование Лапласа по |. Оно дает

 

л

 

 

 

л

 

(3.104)

 

G x x { z ,

S ) =

s G

i S k h ( Z , s)

 

и

_

 

 

_

 

 

 

 

л

=

_

л

л

s)],

(3.105)

 

2oxx{z, s)

sG2[efefe(2, s) — a(z,

где

использована

зависимость

(3.103) и

s — параметр

§ 3.6. Термореологически простые материалы

149

л

преобразования. Исключив еьь. из (3.104) и (З.Г05), по­ лучим

а хх(-Z, s) = —sG lG2a(z, s )/(2G i + G2). (3.106)

С помощью обратного преобразования Лапласа можно

определить функцию R(£,)

в виде

 

 

./?(!) = S ? - 1[G1(s )G 2(s )/(2G 1(s ) +

G2(s ))] .

(3.107)

Используя (3.107), обращаем (3.106) и находим

 

Ъ х (гЛ ) =

- \ Я ( 1 - Ъ ' ) № ( г Л ,) № < % .

(3.108)

 

 

о'1

 

 

 

 

Сведем теперь

(3.108)

с

помощью

(3.91) и

(3.92)

к зависимости

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

физического времени t.

В результате имеем

 

 

 

 

ѵхх (г, 0 = ауу (г, Ц =

 

 

 

 

= —j

R (z, f) —I' (г, т)] (да (г , т)/дт) dr.

(3.109)

о

 

 

 

 

 

 

Тем самым завершено решение задачи для напряжений. Если задано температурное поле, которое предполагает­ ся совместным с начальными условиями, то £ находится из (3.86), R — из (3.107), а затем находятся напряже­ ния из (3.109). Более общее исследование вместе с под­ робными численными результатами этой задачи дали Муки и Стернберг [3.16]. Ли и Роджерс [3.13] указа­ ли, что эту задачу можно решить без использования преобразования Лапласа, которое применялось здесь; затем ту же задачу рассматривал Стернберг [3.19]. Близкую по смыслу, но более трудную задачу об опре­ делении температурных напряжений в цилиндре реша­ ли Локкет и Морленд [3.14].

Ранее уже было указано, что связанная теория тер­ мовязкоупругости, в которой учитывается зависимость функций механических свойств от переменного поля тем­ пературы, непременно должна быть нелинейной. Эта связь между тепловыми и механическими явлениями вы-

150

Гл. 3. Термовязкоупругость

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

Задачи

3.1.Какие модификации следует провести, чтобы учесть несжи­ маемость материалов при термодинамическом выводе линейной тео­ рии в § 3.1? Считайте условия изотермическими и материалы изотроп­ ными.

3.2.Следуя процедуре § 3.2 и используя диссипативное неравен­ ство, выведите необходимые ограничения для изотропных функции интегрирования в выражении для свободной энергии, где эти огра­

ничения следуют из последовательного применения двух ступенчатых функций изменения деформации и температуры.

3.3. Рассуждения в § 3.1 основаны на представлении свободной энергии, включающем функции релаксации. Можно использовать дру­ гой способ вывода, если применить функции ползучести. Проделайте это, вводя преобразование Лежандра между свободной энергией А, отнесенной к единице массы, и функционалом pW в виде

00

со

р А (ец (t — s)) =

о а EU (t) — pW (оц (t s) ) ,

s = 0

s—О

где предполагаются изотермические условия. Далее примите pW в виде

t t

pW (t) = Va j' j Jtikl (2t — г — T)) (00;/ (х)/дх) (дои O l)/*l) dxdr\,

ои

где Jijh i(t) — изотермические функции ползучести, которые принима­ ются в частном виде, при котором суммируются два возможных не­ зависимых аргумента времени. Используя вывод, подобный тому, который использовался в § 3.1, покажите, что окончательное опреде­ ляющее уравнение для деформации имеет вид

Задачи

151

t

 

ъц (0 = J Jlik l (t x) (dakt (x)/dx)

dt,

о

 

адиссипативное неравенство представляется в форме

г(0 > о,

где

i t

 

 

 

г ( 0 оfоf а*- Jtlkl

^

 

д0» ft) dxdx\.

дх

dr]

 

 

При переходе отсюда к соотношениям теории упругости функционал pW становится просто упругой дополнительной энергией. Можно ли в вязкоупругом случае дать функционалу pW интерпретацию в тер­ минах энергии?

 

3.4. В условиях стационарного состояния гармонических колеба­

ний

получите

выражение для локальной энергии, рассеивающей­

ся

за один

цикл,

путем рассмотрения скорости работы

(ReOij) (Ree, j ) ,

если а ц

и е^ -— обычные комплексные гармониче­

ские функции напряжений и деформаций.

3.5. Произведите анализ температурных напряжений для термо­ реологически простой вязкоупругой сферы подобно тому, как это бы­ ло сделано в § 3.6 для плиты. Считайте температурное поле симмет­ ричным относительно точки и заданным по всей сфере в виде Т (г)

(см. работу [3.16]).

3.6. Покажите, что постулат о термореологически простом по­ ведении эквивалентен утверждению, что времена релаксации т,- в функции релаксации вида

G (t) = G0 + Е G( ë ~ t/Xi l= 1

связаны с температурой зависимостью

Хі = XiF(T),

где г , — времена релаксации при исходной температуре Г0, а F(T ) — функция от температуры, такая, что

F (T 0) = 1.

Чтобы подтвердить такую интерпретацию, начните с того, что выбе­ рите F(T ) в виде

F (T ) = 10-ftx).

Какие выводы отсюда следуют относительно механических моделей при такой интерпретации термореологически простого материала?

152

Гл. 3. Термовязкоупругость

 

СПИСОК ЛИТЕРАТУРЫ

3.1.Biot М. А., Linear Thermodynamics and the Mechanics of Solids, Proc. 3rd U. S. Nat. Congr. Appl. Mech., 1, 1958.

3.2.Breuer S., Onat E. T., On the Determination of Free Energy in Linear Viscoelastic Solids, Z. Angew. Math. Phys., 15, 184 (1964).

3.3.Christensen R. M., Naghdi P. M., Linear Non-Isothermal Visco­

elastic Solids, Acta M echanica, 3, 1 (1967).

3.4.Crochet M. J., Naghdi P. M., A Class of Simple Solids with Fading Memory, Int. J. Eng. Sei., 7, 1173 (1969).

3.5.Crochet M. J., Naghdi P. M., On «Thermorheologically Simple» Solids, Proc. IUTAM Symp. Thermoelasticity, p. 59, New York, 1970.

3.6.Eringen A. C., Irreversible Thermodynamics and Continuum Me­ chanics, Phys. Rev., 117, 1174 (1960).

3.7.Ferry J. D., Mechanical Properties of Substances of High Mo­ lecular Weight; VI. Dispersion in Concentrated Polymer Solu­

tions and Its Dependence on Temperature and Concentration,

J. Amer. Chem. Soc., 72, 3746 (1950).

3.8.Ferry J. D., Viscoelastic Properties of Polymers, 2nd ed., New York, 1970. Русский перевод 1-го издания: см. стр. 50.

3.9.Gurtin М. Е. Herrera I., On Dissipation Inequalities and Linear Viscoelasticity, Quart. Appl. Math., 23, 235 (1965).

3.10.Hunter S. C., Tentative Equations for the Propagation of Stress, Strain, and Temperature Fields in Viscoelastic Solids, J. Mech. Phys. Solids, 9, 39 (1961).

3.11.Kovacs A. J., La Contraction Isotherme du Volume des Polyméres Amorphes, J. Polym , Sei., 30, 131 (1958) .

3.12.Leaderman H., Elastic and Creep Properties of Filamentous Ma­ terials and Other High Polymers, p. 175, Washington, 1943.

3.13.Lee E. H., Rogers T. G., Solution of Viscoelastic Stress Analysis

Problems Using Measured Creep or Relaxation Functions,

J. Appl. Mech., 30, 127 (1963).

3.14.Lockett F. J., Morland L. W., Thermal Stresses in a Viscoelastic,

Thin-Walled Tube with Temperature Dependent Properties, Int.

J. Eng. Sei., 12, 879 (1967).

3.15.Morland L. W., Lee E. H., Stress Analysis for Linear Viscoelastic

Materials with Temperature Variation, Trans. Soc. Rheol., 4,

233 (1960).

3.16.Muki R., Sternberg E., On Transient Thermal Stresses in Visco­ elastic Materials with Temperature Dependent Properties, J. Appl. Mech., 28, 193 (1961).

3.17.Schapery R. A., Application of Thermodynamics to Thermome­ chanical Fracture and Birefringent Phenomena in Viscoelastic Media, /. Appl. Phys., 35, 1451 (1964).

3.18.Schwarzl F., Staverman A. J., Time-Temperature Dependence of Linear Viscoelastic Behavior, J. Appl. Phys., 23, 838 (1952).

3.19.Sternberg E., On the Analysis of Thermal Stresses in Viscoelastic Solids, Proc. 3th Symp. Nav. Struct. Mech., 348, New York, 1964.

3.20 Truesdell C.,

Toupin

R. A., in

«Handbuch der Physik» (Flügge

S., ed.), Band

3, No.

1, Berlin,

1960.

Глава 4

РАСПРОСТРАНЕНИЕ ВОЛН

Проблема распространения волн в вязкоупругих ма­ териалах имеет много различных важных аспектов; не­ которые из них рассматриваются здесь для изотропной среды. В первом параграфе этой главы исследуется изо­ термическое распространение возмущений в полубесконечных стержнях. Некоторые из результатов такого ис­ следования обобщены на случай распространения пло­ ских возмущений в трехмерной среде. В частности, най­ дены две возможные скорости распространения плоско­ стей разрыва напряжений и деформаций в бесконечной среде. Исследование особенностей распространения гар­ монических волн рассмотрено в § 4.3 и 4.4.

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

Другая важная задача распространения волн отно­ сится к исследованию распространения плоских гармо­ нических волн в неограниченной среде. Этот случай ана­ лизируется, и рассматриваются эффекты распростране­ ния связанных термических волн с применением связанной теории термовязкоупругости, развитой в §3.1. Приложение этих результатов для частного случая за­ дачи с изотермическими условиями дало возможность получить в первом параграфе скорости распространения

154

Гл. 4. Распространение волн

и скорости затухания изотермических гармонических волн в трехмерной среде. Для аналогичной задачи тео­ рии упругости импульсы могут распространяться в сре­ де без изменения формы и без затухания. Результаты этого параграфа показывают, что скорости распростра­ нения гармонических волн в вязкоупругой среде зависят от частоты. Следовательно, импульсы, составленные из некоторого спектра гармонических волн, испытывают влияние дисперсии, и форма импульса не сохраняется. Этот эффект, вместе с эффектом затухания вязкоупру­ гих волн, обнаруживает глубокое различие между осо­ бенностями распространения волн в вязкоупругой и упругой средах.

В§ 4.3 рассматриваются только плоские волны; в этом случае направления распространения и затухания сов­ падают. В § 4.4 показывается, что это ограничение не всегда имеет место. Здесь рассматривается задача от­ ражения гармонических волн от свободной поверхности полупространства и показывается, что один из двух ти­ пов отраженных волн действительно имеет различные направления распространения и затухания.

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

§4.1. Изотермическое распространение волн

Одноосный случай

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

§ 4.1. Изотермическое распространение волн

155

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

На сингулярной поверхности, при переходе через ко­ торую напряжение претерпевает разрыв, уравнение дви­ жения теряет силу. Однако условие для скачка, которое справедливо на сингулярной поверхности, можно выве­ стииз глобальной формы линейного уравнения количе­ ства движения (см., например, Томас [4.21]). Такое ус­ ловие для скачка в данном случае имеет вид

[<*(Х, t)] = — рѵ[ди(х, t)/dt];

(4.1)

здесь о и и — одноосные напряжение и перемещение со­ ответственно, р — плотность, V— скорость распростра­ нения сингулярной поверхности, а [h ] означает скачок переменной h(x, t) при переходе через сингулярную по­ верхность, определяемый формулой

\h ] = lim h (х+, /+) — lim h (x~ , t~),

* 4"-УХ,

X ->x

i+^t

t—^t

где x+, t+ и xr, t~ относятся к противоположным сторо­

нам сингулярной поверхности в точке

с

координата­

ми X, t.

 

 

 

 

 

Одноосное

вязкоупругое определяющее

соотношение

имеет вид

 

t

 

 

 

 

 

 

 

 

а (x ,t)=

J Е (t — т) (de (х, x)/dx)dx.

 

(4.2)

Интегрируя (4.2) по частям и подставляя

результат в

(4.1), имеем

 

 

 

 

 

Е ( 0) ди (X, t)

± E ( t - x ) ^ J l dx

 

 

 

дх ч

dt

дх

 

 

 

 

 

= — ру '

ди (х, t)

(4.3)

 

 

 

dt

 

где использована зависимость е — ди/дх между дефор­ мацией и перемещением. Рассмотрим распространение

156 Гл. 4. Распространение волн

разрывов напряжения и деформации, которые приводят к распространению разрыва в первой производной пе­ ремещения. Такие типы волн обычно называют ударны­

ми волнами. Предполагая,

что функция dE (t)/dt

огра­

ничена и непрерывна при

0

получаем, что

ин­

тегральный член в (4.3)

не вносит вклада в скачок,

так

как функция ди (х ,х )/дх

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

т = t , и интеграл не зависит от du(x,t)/dx . Таким

обра­

зом, уравнение (4.3) приводится к виду

 

 

Е ( 0) [ди (х, і)/дх ]

=

рѵ[ди(х, t)/dt],

 

(4.4)

Из работы Томаса [4.21] известны кинематические условия, связывающие скачки в первых производных перемещения:

[ди(х, t)/dt] =

ѵ\ди(х, t)/dx ].

(4.5)

Комбинируя зависимости (4.4) и

(4.5),

получаем

 

Е (0) [ди(х, і)/дх\

= рѵ

2[ди(х,

t)/dx];

 

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

о = (£(0)/р)7*.

(4.6)

Соотношение (4.6) определяет скорость распростране­ ния одноосной вязкоупругой ударной волны, если из­ вестны начальное значение одноосной функции релакса­ ции и плотность.

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

Для решения задач о распространении волн может использоваться как преобразование Фурье, так и преоб­ разование Лапласа. Если используется преобразование