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

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

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

Задачи

197

[ѴхЬ _

— Р (VsQ— 0

х

üPo { (TsO— l) 2 + 4"^/0Vso}

 

 

xexD(_

2-y^o (Я0 +

2ц0)

y)

 

(4.166)

l vy h _ _

_

?PVfo_ _

i

i

+

0(4.167)

_ _

ex_p

J _

 

 

( ? 5 0

ü

Hfl {(VsO —

l)2 +

4V/o Vso}

1

 

 

2VsoHo

 

J

и

 

 

 

 

 

 

 

 

 

 

 

Ы і

PVfQ (VsO -

0

 

 

 

 

 

 

 

 

V

Ho {(Vso —

1)2 +

Yso}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xexp

 

(Vfo +

Q

( h

+

2Ht)

 

(4.168)

 

 

 

2Vfo (^o ~b 2,u0)

У

 

 

 

 

 

 

 

Подобные выражения

можно

найти

и

для

разрывов

в компонентах напряжения. Соотношения

(4.165) —

(4.168)

аналогичны результатам решения соответствую­

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

Подобная задача о нагрузке, движущейся по вязкоупругому слою конечной толщины, с учетом отражения волн рассмотрена в работе [4.7].

Задачи

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

198

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

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

4.3. В задаче о течении сдвига в § 4.2 проведите полный вывод решения для напряжения для вязкоупругого материала обще­ го вида. Дайте частный случай этих результатов для функции релак­ сации, представленной моделью Максвелла, и сравните результаты с (4.61). В заключение получите полное решение соответствующей упругой задачи, для которой модуль сдвига равен начальному значе­ нию функции релаксации при сдвиге.

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

4.5.Распространите метод § 4.4 на учет отражения и преломле­ ния падающих гармонических волн на поверхности раздела между двумя вязкоупругими средами (см. работу [4.17]).

4.6.Пренебрегая инерционными членами, получите квазистати­ ческое решение задачи в § 4.5, вводя движущуюся нагрузку на вяз­ коупругом полупространстве (при плоской деформации).

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

4.1.Achenbach J. D., Vibrations of a Viscoelastic Body, АІАА J., 5, 1213 (1967). Русский перевод: Ракетная техника и космонавти­ ка, 5, № 6, 221 (1967).

4.2.Achenbach J. D„ Reddy D. P., Note on Wave Propagation in Linearly Viscoelastic Media, Z. Angew. Math. Phys., 18, 141

(1967).

4.3.Berry D. S., A Note on Stress Pulses in Viscoelastic Rods, Phil. Mag. [8], 3, 100 (1958).

4.4.Berry J. G., Naghdi P. M., On the Vibration of Elastic Bodies Having Time-Dependent Boundary Conditions, Quart. Appl. Math., 14, 43 (1956).

4.5.Christensen R. M., Application of the Method of Time-Dependent

Boundary Conditions in Linear Viscoelasticity, I. Appl. Mech., 34, 503 (1967).

4.6.Chu В. T., Stress Waves in IsotropicLinear Viscoelastic Mate­ rials, J. Alec., 1, 439 (1962).

4.7.Chu В. T., Response of Various Material Media to High-Velocity Loadings; I. Linear Elastic and Viscoelastic Materials, J. Mech.- Phys. Solids, 13, 165 (1965).

4.8.Cooper H. F., Jr., Reflection and Transmission of Oblique Plane Waves at a Plane Interface between Viscoelastic Media, I. Acoust. Soc. Amer., 42, 1064 (1967).

4.9.Cooper H. F., Jr., Reiss E. L., Reflection of Plane Viscoelastic Wa-

Список литературы

 

 

 

199

ves from Plane Boundaries, J. Acoust.

Soc.

Amer.,

39,

1133

(1966).

 

 

the

Li­

4.10. Fisher G. M. C., Gurtin M. E., Wave Propagation in

near Theory of Viscoelasticity, Quart.

Appl.

Math.,

23,

257

(1965).

Rods, J. Appl. Mech.,

4.11. Herrmann G., Forced Motion of Elastic

21, 221 (1954).

 

 

 

 

4.12. Hunter S. C., Tentative Equations for the Propagation of Stress,

Strain, and Temperature Fields in Viscoelastic Solids,

J.

M ech,

Phys. Solids, 9, 39 (1961).

 

 

4.13. Hunter S. C., The Transient Temperature Distribution

in

a Se­

mi-Infinite Viscoelastic Rod Subject to Longitudinal Oscillations,

Int. ]. Eng. Sei., 5, 119 (1967).

4.14. Knauss W. G., Uniaxial Wave Propagation in a Viscoelastic Mate­ rial Using Measured Material Properties, /. Appl. Mech., 35,

449 (1968).

4.15.Lee E. H., Kanter L, Wave Propagation in Finite Rods of Vis­ coelastic Material, J. Appl. Phys., 24, 1115 (1953).

4.16.Lee E. H., Morrison J. A., A Comparison of the Propagation of

Longitudinal Waves in Rods of Viscoelastic Materials, J. Polytn. S ei, 19, 93 (1956).

4.17.Lockett F. J., The Reflection and Refraction of Waves at an In­ terface between Viscoelastic Materials, J. Mech. Phys. Solids, 10, 53 (1962).

4.18.Mindlin R. D., Goodman L. E., Beam Vibrations with Time-Depen­ dent Boundary Conditions, J. Appl. Mech., 17, 377 (1950).

4.19.Petrof R. C., Gratch S., Wave Propagation in a Viscoelastic Ma­ terial with Temperature-Dependent Properties and Thermomecha­ nical Coupling, I. Appl. Mech., 31, 423 (1964).

4.20. Sackman J.

L., Kaya I., On the Propagation of Transient Pulses

in Linearly

Viscoelastic Media, J. Mech. Phys. Solids, 16, 349

(1968).

 

4.21.Thomas T. Y., Plastic Flow and Fracture in Solids, New York, 1961. Русский перевод: Томас T., Пластическое течение и раз­ рушение твердых тел, «Мир», М., 1964.

4.22.Tsai Y. М., Kolsky Н., Surface Wave Propagation for Linear Viscoelastic Solids, /. Mech. Phys. Solids, 16, 99 (1968).

4.23.Valanis К- C., Propagation and Attenuation of Waves in Linear Viscoelastic Solids, J. Math. Phys., 44, 227 (1965).

4.24.Valanis К- C., Chang S., Stress Wave Propagation in a Finite Viscoelastic Thin Rod with a Constitutive Law of the Hereditary Type, in «Developments in Theoretical and Applied Mechanics» (Huang T. C., Johnson M. W., Jr., eds.), New York, 1965.

4.25.Valanis К- C., Sun С. T., Axisymmetric Wave Propagation in a Solid Viscoelastic Sphere, Int. J. Eng. Sei., 5, 939 (1967).

4.26.Wolosewik R. M., Gratch S., Transient Response in a Viscoelastic Material with Temperature-Dependent Properties and Thermome­ chanical Coupling, J. Appl. Mech., 32, 620 (1965).

Глава 5

ОБЩИЕ ТЕОРЕМЫ И ФОРМУЛИРОВКИ

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

§ 5.1. Единственность решения краевой задачи

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

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

Соответствующие линейные уравнения из гл. 3 для динамической неизотермической анизотропной связанной теории вязкоупругости сводятся к уравнениям (3.45)— (3.52), причем уравнение движения (3.45) учитывает инерционные члены. Предполагается, что существует преобразование Лапласа этих соотношений с нулевыми начальными условиями при £ < 0 . Преобразованные со­ отношения даются зависимостями (3.53)—-(3.57), кото­ рые для удобства повторяются ниже. Уравнения движе­ ния, соотношения между деформациями и перемещени­ ями и определяющие соотношения для напряжений соот­ ветственно таковы:

Oij,j -\-Fi = s 2pui,

( 5 . 1 )

§ 5.1. Единственность решения краевой задачи

201

е(7 — 1!"і fa ij I ч-j.i)

(5.2)

и

 

Ои = sGiiKi Чі — Sfpа Ѳ,

(5.3)

где s — переменная преобразования, а Ѳ(х{, t) — инфини­ тезимальное изменение температуры по отношению к ис­ ходной температуре Т0. Уравнение теплопроводности имеет вид

(ku/To) Ѳ.</ = s2mQ -! s2(р(/8/;,

(5.4)

где kn и функции релаксации т (і) и cpij(t) определяют­ ся механическими свойствами материала. Граничные ус­ ловия записываются как

a ii пі ~

S,

щ =

А,

CDI 1!

CD >|

kijQ jn j =

0

на Ва ,

на Ви>

на В!

на в%,

(5.5)

(5.6)

(5.7)

(5.8)

где В а , В и, В\ и В%— соответствующие подобласти гра­

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

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

К раевая задача изотропной связанной термовязкоуп­ ругости, определяемая уравениями и граничными усло­ виями (3.45) — (3.52), имеет единственное решение, если существует преобразование Л апласа всех переменных поля с преобразованными уравнениями поля и гр а ­ ничными условиями (5.1) — (5.8) и если начальные зна­

чения функций релаксации являются положительно оп­

ределенными, а коэффициент теплопроводности поло­

жительно полуопределенным, т. е. если

6цы(0)уцум > 0,

т ( 0 ) > 0,

к ц у і у ^ б ,

где уц произвольный

симметричный

тензор второго

порядка,

 

 

202 Гл. 5. Общие теоремы и формулировки

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

переменные поля имеют вид

 

 

«і1»,

е<;>,

а<)>,

ѲО>

 

и

 

 

u f ),

е<2),

а<2>,

Ѳ<2>

и считаются непрерывными функциями времени и коор­

динат.

Обозначим разности этих решений через ud., e f ,

и 0d, так что

 

 

 

 

 

 

“і М ) М ’ М М Р М ) ’

 

 

4

M

) - - М М М ) ?

М ) >

(5.9)

 

4

К . О =

<4' М

М ) ?

К -О

 

и

 

 

 

 

 

 

 

 

Q d ( X i ,

t) =

0(9 ( X i ,

t) — Ѳ<2)(Хі, / ) .

 

Заметим, что разности решений 4

(хг-, ^),

(лгг-, /),

od(X{,

t) и 0d(Xi,

удовлетворяют однородным гранич­

ным условиям. В силу этого можно записать следующее соотношение:

J

ad. п. ud. da = 0

(5.10)

в

 

 

для преобразований

Лапласа переменных.

Применив

к (5.10) теорему Гаусса — Остроградского и использо­ вав уравнения (5.1) и (5.2), отвечающие разности ре­ шений, находим

j [ps2«d ud +

crf ed.] dv =-- 0.

(5.11)

V

 

 

 

 

Подставляя в (5.11) разностную форму (5.3),

получаем

J [Рs2ud ud +

sG.jkj edf

kld — sep.j ed 0d] dv =

0.

(5.12)

V

 

 

 

 

Последний член в

(5.12) можно исключить. С этой

це­

лью используем уравнение

теплопроводности

(5.4)

для

§ 5.1. Единственность решения краевой задачи

203

разности решений, умножив обе его части на 0d. Тогда получим

[ktiIT0) Щ

, = s2in (№ f +

№.

(5.13)

Первый член в (5.13) необходимо

выразить

в другой

форме. Чтобы сделать это, укажем тождество

f k . fi d№. n. da =

0,

(5.14)

в

 

 

 

которое справедливо,

поскольку

0d= O

на В х и

k iß dt( tij= 0 на В 2, где В у и В 2— дополняющие друг дру­

га части границы В. Применяя к (5.14)

теорему Гаус­

са — Остроградского, находим, что

 

=

(5-15)

Теперь, используя (5.15), разрешим (5.13) относительно' e^0d и подставим результат в (5.12). Это дает

j [pshtfu* +sG .jH l di f 4 , + (ktj/sT0) Щ ,

+ sm (0d)2] dv =

0.

 

 

(5.16)

Рассмотрим действительные значения s,

такие,

что

s > s 0, где so — самая удаленная от

начала

координат

особенность решения справа на действительной оси. На­ помним, что

lim s G ^ (s) = Gijkl (0)

S->c°

иlim sm (s) = tn (0). (5.17)

S-»-co

Теперь потребуем, чтобы начальные значения функций релаксации были положительно определенными:

Оцы{0)уізуы > 0, т ( 0 ) > 0 .

(5.18)

Отсюда следует, что sGijki(s)уі]уы >0 и srn(s) > 0

для

достаточно больших s. Если, кроме того, потребовать, чтобы выполнялось условие k a y iy ^ Q , то при достаточ­ но большом s уравнение (5.16) будет состоять только из одних неотрицательных членов, и, для того чтобы оно удовлетворялось, мы должны иметь

№ — № = 0 при достаточно больших s > s0. (5.19)

204

І'л. 5. Общие теоремы и формулировки

Теорема из книги Дейча [5.4] утверждает, что если пре­ образование Лапласа от непрерывной функции равно пулю на бесконечном множестве точек действительной оси, которые образуют арифметическую прогрессию, то и сама функция тождественно равна нулю. Применяя эту теорему к нашему случаю, обнаруживаем, что раз­ ности функций ud. и 0d должны быть тождественно рав­

ны нулю:

 

ud. = Qd = 0.

(5.20)

Поскольку решение для разностей тождественно рав­ но нулю, решение краевой задачи связанной термовяз­

коупругости должно быть единственным.

Это доказа­

тельство

достаточности следует из

предположений

ЬцУіУі ^О ,

m(0) > 0 и G ijki(0)ynyhi>0.

Примечательно,

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

Использование преобразования Лапласа для уста­ новления теоремы единственности является здесь в не­ котором смысле искусственным. Хотя методы интеграль­ ных преобразований служат весьма ценным инструмен­ том для решения краевых задач, их применение при ус­ тановлении общих теорем менее привлекательно из-за ограничений, с которыми оно связано. В этом парагра­ фе преобразование Лапласа используется только из-за простоты его приложения именно в этом частном дока­ зательстве; оно не является существенным элементом доказательства единственности. Доказательство изотер­ мической теоремы единственности в § 2.2 производилось без использования интегральных преобразований.

В § 2.2 были даны ссылки на несколько различных изотермических теорем единственности. Стернберг и Гёртин [5.14] дали теорему единственности для неизотерми­ ческого случая, применимую в условиях, отличных от приведенных здесь.

§ 5.2. Представление через функции перемещения

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

§ 5.2. Представление через функции перемещения

205

тора Галеркина и функции напряжений Папковича — Нейбера. Обратимся к обобщению этих методов теории упругости на случай вязкоупругости. Такое обобщение связано с использованием введенной в § 1.2 свертки Стильтьеса.

Напомним, что в теории упругости уравнения равно­ весия в перемещениях без учета объемных сил име­ ют вид

цѴ2и + [ц /(1 — 2ѵ) ] V V •и = 0,

(5.21)

где и — вектор перемещения, V — оператор

Лапласа

и Ѵ - и — скалярное произведение. Решение этого урав­ нения с помощью вектора Галеркина записывается как

и = (1/2ц)[2(1 — v )V 2g — V V -g ],

(5.22)

где

 

V 4g = 0.

(5.23)

Такимобразом, решение уравнений равновесия теории упругости сводится к задаче отыскания бигармониче­ ской векторной функции g. Альтернативная процедура связана с использованием функции напряжений Папко­ вича— Нейбера; при этом вектор перемещения име­ ет вид

и = (1/2ц) [Ѵ(<р +

г-ф)

— 4(1 — ѵ)ф],

(5.24)

где г — вектор, соединяющий

рассматриваемую

точку

с началом координат, и

 

 

 

Ѵ 2ф =

Ѵ 2ф = 0.

(5.25)

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

Для задач вязкоупругости уравнения квазистатического равновесия в форме, аналогичной (5.21), имеют вид

V 2u * dGi -f- ѴзѴѴ -u

* d (G x+

2G2) = 0 , (5.26)

где используется обозначение

свертки

Стильтьеса из

§ 1.2. Это представление можно проверить, если исполь­ зовать уравнение ац, j — 0 вместе с (2.3), приняв Gx(t) — = 2 p ,(f), G2( t ) = 3 k ( t ) , а также соотношения между де­ формациями и перемещениями.

206

Гл. 5. Общие теоремы и формулировки

Обобщение упругого вектора Галеркина на вязкоуп­ ругость дается зависимостью

u = 2V 2g *d (2 G i + G2) — V V -g * d(G\ + 2G2),

(5.27)

где

 

V 4g = 0.

(5.28)

Для того чтобы доказать, что такое представление удов­ летворяет уравнениям равновесия, выражение (5.27) можно прямо подставить в (5.26). Чтобы проделать это, получим из (5.27) две следующие зависимости:

Ѵ 2и = 2 [ V 4g

* d(2G i +

G2) ] —

 

 

— V 2[V V -g

* d {G { + 2G2)]

 

(5.29)

и

 

 

 

 

 

 

 

V V •u = V V •[2V 2g * d (2G, +

G2) —

 

— V V -ff *d(G i + 2G2) j.

 

 

(5.30)

Заметим, что первый член в правой части (5.29)

обраща­

ется в нуль в силу зависимости

(5.28). Затем,

подстав­

ляя (5.29) и (5.30) в (5.26), получаем

 

 

 

 

— V 2[V V -g

* d(G i +

2G2) ] *

dGi +

 

+ ѴзѴѴ- [2V 2g*

d (2 G { +

G2)]

. d{G i +

2G2) —

—Ѵ зѴ Ѵ Д Ѵ Ѵ -g * d(G i + 2G2)] *

d(G i +

2G2) = 0 .

 

 

 

 

 

 

 

(5.31)

Воспользуемся известным из векторного

анализа тож­

деством V 2V d iv =V d iv V 2= V d iv Vdiv

и

свойствами

ассоциативности и дистрибутивности свертки Стильтье-

са (см. § 1.2), применив их к (5.31).

Это сводит

(5.31)

к простой форме

 

 

V 2V d iv g * [—d (G i+ 2 G 2) * d G i - f

 

-frf(G j) * d(G\ + 2G2)]

= 0 .

(5.32)

Наконец, воспользуемся свойством коммутативности свертки Стильтьеса, чтобы показать, что зависимость (5.32) удовлетворяется, и доказательство закончено.