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

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

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

§ 2.13. Итоги и выводы

107

Фурье. Таким приближенным методам обращения посвя­ щена обширная литература (см., например, [2.6] и [2.30]). Пример такого приближенного метода дан в при­ ложении В. Следует также отметить, что в некоторых случаях краевые задачи вязкоупругости можно свести к интегрированию интегральных уравнений стандартных типов. Хотя для получения аналитических решений уже описанным путем может использоваться метод интег­ ральных преобразований, другой путь получения реше­ ния состоит в непосредственном численном интегрирова­ нии интегральных уравнений.

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

Наконец, для тех задач,

в

которых часть

границы

В и, где заданы перемещения,

и

часть границы

В а, где

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

108

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

Задачи

2.1. Сформулируйте процедуру решения первой и второй краевых задач для несжимаемого вязкоупругого материала с применением метода разделения переменных. При этом предположите, что имеют место квазистатические условия.

2.2. Длинный вязкоупругий цилиндр с внутренним радиусом а

и

внешним радиусом b подвергается

действию внутреннего давления

p

(t), тогда как его внешняя граница

свободна от нагрузки. Исполь­

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

2.3. Сфера с внутренним радиусом а и внешним радиусом b подвергается действию внутреннего давления p(t) = p 0h (t), тогда как ее внешняя граница свободна от нагрузки. Получите преобразование Лапласа решения для компонент перемещения и напряжения для вязкоупругого материала общего вида в квазистатических условиях. Выполните обращение при дальнейшем допущении о том, что функ­ ция релаксации при объемном расширении является упругой, тогда как функция релаксации при сдвиге состоит из одного экспоненциального члена.

2.4. Опишите квазистатическое поведение вязкоупругого полу­ пространства, на которое действует внезапно приложенная сосредо­ точенная нормальная нагрузка p(t) = p 0h(t) в некоторой точке, в то время как остальная поверхность остается свободной от нагрузки (это соответствует задаче Буссинеска из теории упругости). Считай­ те, что функция релаксации при объемном расширении является уп­ ругой, тогда как функция релаксации при сдвиге та же, что и для обобщенной модели Максвелла. Получите результаты для частного случая, когда задана функция релаксации при сдвиге, соответствую­ щая модели Максвелла.

2.5. Упругая балка с изгибной жесткостью Е еІе имеет вязкоупру­ гое покрытие, которое сообщает ей дополнительную жесткость. Пре­ образование Лапласа от вязкоупругого вклада в изгибную жесткость

балки равно E v(s ) I v. Модифицируйте анализ задачи о вдавливании штампа из § 2.12 таким образом, чтобы он охватывал балки такого типа. Найдите перемещение центра балки, когда нагрузка на штамп

задана в виде p ( t ) — pah(t) и когда функция Е ѵ отвечает одному убывающему экспоненциальному члену в функции релаксации.

2.6. Проведите анализ задачи о вдавливании штампа, подобный проведенному в § 2.12, но вместо вязкоупругой балки рассмотрите круглую вязкоупругую пластинку, для которой классическое уравне­ ние, связывающее нагрузку и прогибы упругой пластинки, модифици­ ровано так, чтобы описывать поведение вязкоупругой пластинки. Считайте края пластинки свободно опертыми.

2.7. Предлагается измерить комплексный модуль путем исполь­ зования эксперимента со свободными колебаниями. Рассматривается

Задачи

109

вязкоупругая консольная балка с грузом на свободном конце, испы­ тывающая свободные колебания. Если пренебречь массой балки по сравнению с сосредоточенной массой на конце, то можно считать, что система имеет одну степень свободы. Определите действительную и мнимую части комплексного модуля Е * (т ) в зависимости от на­ блюдаемой частоты колебаний ш и логарифмического декремента за­ тухания. Пусть 7 — момент инерции балки, L — ее длина, т — масса на конце, и пусть используется классическая теория балок; можно принять Е"/Е'<^ 1. Получите модификации, которые следует провести в этой процедуре для случая упругой балки с вязкоупругим покры­ тием. В этом случае для вязкоупругого материала не используйте предположения Е " /Е '^ :і.

2.8. Обратите соотношение (2.127) для функции р.(г'), представ­ ленной одним убывающим экспоненциальным членом \іе~Ѵх со вре­ менем релаксации т. Интерпретируйте полученное решение для сво­

бодных колебаний в двух случаях о?п ^ 1/(4т2.) Убедитесь в том, что

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

2.9. Используя результаты задачи (2.8), но с заменой р(7) на E (t), найдите решение задачи о свободных колебаниях свободно опертой вязкоупругой балки при равномерном распределена скоро­ стей в начальном состоянии. Используйте классическую теорию ба­ лок с заменой модуля Юнга Е на преобразование функции релакса­

ции sE (s).

2.10. Пусть функция релаксации для некоторого материала мо­ жет быть адекватно описана формулой

G (t) = G/(3 + lg t), 0,01 < t < 10.

Чтобы использовать эту функцию в аналитических исследованиях, желательно представить G(t) в интервале— l ^ l g l ^ l в виде ряда из затухающих экспонент. Примем

А4

G ( 0 / 0 = 0 , 2 0 + ^ G je x p {— < /[1 5 -10(_ 4 + ') ] } ,

;=i

где первый член 0,20 является приближением асимптотического зна-

Л

для длительного времени. Времена релаксации приняты

чения G (t)/G

в интервале

описания G (t). Найдите Gj, і = 1 , 2, 3, 4, приравнивая

 

А

два выражения G (t)jG в точках lg t — —2, — 1, 0, 1, и сравните задан­ ную функцию с приближенной.

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

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

2.1. Alfrey Т., Nonhomogeneous Stresses in Viscoelastic Media. Quart. Appl. Math., 2, 113 (1944).

2.2. Barberan J., Herrera J., Unigueness Theorems and the Speed of Propagation of Signals in Viscoelastic Materials, Arch. Ration. Mech. Anal., 23, 173 (1966).

2.3.Brüll M. A., A Structural Theory Incorporating the Effect of Ti­ me-Dependent Elasticity, Proc. 1st Midwestern Conf. Solid Mech., 141, 1953.

2.4.Calvit IT. H., Numerical Solution of the Problem of Impact of a Rigid Sphere onto a Linear Viscoelastic Half-Space and Compa­ rison with Experiment, Int. J. Solids Structures, 3, 951 (1967).

2.5. Christensen R. M., Schreiner R.

N.,

Response

to Pressurization

of

a Viscoelastic Cylinder with

an

Eroding

Internal

Boundary,

AIAA J., 3, 1451 (1965). Русский

перевод: Ракетная

техника

и

космонавтика, 3, № 8, 99 (1965).

 

 

 

 

2.6.Cost Т. L., Approximate Laplace Transform Inversions in Visco­ elastic Stress Analysis, AIAA J., 2, 2157 (1964). Русский перевод:

Ракетная техника и космонавтика, 2, № 12, 175 (1964).

2.7.Edelstein W. S., Gurtin М. Е., Uniqueness Theorems in the Li­ near Dynamic Theory of Anisotropic Viscoelastic Solids, Arch. Ration. Mech. Anal., 17, 47 (1964).

2.8.Essenburg F., On Surface Constraints in Plate Problems, J. Appl M ech., 29, 340 (1962).

2.9.Gottenberg W. G., Christensen R. M., An Experiment for Deter­ mination of the Mechanical Property in Shear for a Linear Isot­ ropic Viscoelastic Solid, Int. J. Eng. Sei., 2, 45 (1964).

2.10.Gottenberg W. G., Christensen R. M., Prediction of the Transient Response of a Linear Viscoelastic Solid, J. Appl. Mech., 33, 449 (1966) .

2.11. Gradowczyk M. H., Moavenzadeh F., Characterisation of

Linear

Viscoelastic Materials, Trans. Soc. Rheol., 13,

173 (1969).

 

2.12. Graham G. A. C., The Contact Problem in the Linear Theory of

Viscoelasticity, Int. J. Eng. S e i, 3, 27 (1965).

 

 

2.13. Graham G. A. C., The Contact Problem in

the Linear

Theory

of Viscoelasticity When the Time Dependent Contact Area Has

Any Number of Maxima and Minima, Int. J.

Eng. S ei,

5, 495

(1967) .

 

 

2.14.Graham G. A. C., The Correspondence Principle of Linear Vis­ coelasticity Theory for Mixed Boundary Value Problems Invol­ ving Time-Dependent Boundary Regions, Quart. Appl. Math., 26, 167 (1968).

2.15.Gurtin M. E., Sternberg E., On the Linear Theory of Viscoelastici­ ty, Arch. Ration. Mech. Anal., 11, 291 (1962).

2.16.Huang N. C., Lee E. H., Rogers T. G., On the Influence of Visco­ elastic Compressibility in Stress Analysis, in «Proc. 4th Int. Cong. Rheol.» (Lee E. H., ed.), New York, 1965.

2.17.Hunter S. C., The Hertz Problem for a Rigid Spherical Indentor and a Viscoelastic Half-Space, J. Mech. Phtis. Solids, 8, 219 (1960).

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

111

2.18.Hunter S. C., The Rolling Contact of a Rigid Cylinder with a Viscoelastic Half-Space, J. Appl. Mech., 28, 611 (1961).

2.19.Hunter S. C., The Solution of Boundary Value Problems in Linear Viscoelasticity, Proc. 4th Symp. Nav. Structural Mech., 257, Ox­ ford, 1965.

2.20.Jahnke E., Emde F., Tables of Functions, New York, 1945. Рус­ ский перевод: Янке Э., Эмде Ф., Таблицы функций с формула­ ми и кривыми, Гостехиздат, 1948.

2.21.Lee Е. Н., Stress Analysis in Viscoelastic Bodies, Quart. Appl. Math., 13, 183 (1955).

2.22.Lee E. H., Radok J. R. M., The Contact Problem for Viscoelastic Bodies, J. Appl. Mech., 27, 438 (1960).

2.23.Lockett F. J., Interpretation of Mathematical Solutions in Vis­ coelastic Theory Illustrated by a Dynamic Spherical Cavity Prob­ lem, J. Mech. Phys. Solids, 9, 215 (1961).

2.24.Lubliner J., Sackman J. L., On Uniqueness in General Linear Vis­ coelasticity, Quart. Appl. Math., 25, 129 (1967).

2.25.Morland L. W., A Plane Problem of Rolling Contact in Linear Viscoelasticity Theory, J. Appl. Mech., 29, 345 (1962).

2.26.Morland L. W., Exact Solutions for Rolling Contact between Viscoelastic Cylinders, Quart. J. Mech. Appl. Math., 20, 73 (1967).

2.27.Odeh F., Tadjbakhsh L, Uniqueness in the Linear Theory of Vis­

coelasticity, Arch. Ration. Mech. A n al,

18,

244

(1965).

2.28. Onat E. T.,

Breuer S., On Uniqueness in

Linear Viscoelasticity,

in «Progress

in Applied Mechanics»

(Drucker

D. C., ed.) (The

Prager Anniversary Volume), p. 349, New York,

1963.

2.29. Read W. T., Stress Analysis for Compressible Viscoelastic Ma­

terials, /.

Appl. Phys., 21, 671 (1950).

 

2.30. Schapery

R.

A., Approximate Methods of Transform

Inversion

for Viscoelastic Stress Analysis, Proc. 4th U. S. Nat. Cong. Appl.

Mech., 1075

(1962).

 

2.31. Sips R., General Theory of Deformation of Viscoelastic Substan­

ces, J. Polym . Sei., 9, 191 (1951).

 

2.32. Sokolnikoff

I. S., Mathematical Theory of Elasticity,

2nd ed.,

New York, 1956.

2.33.Struik L. С. E., Free Damped Vibrations of Linear Viscoelastic Materials, Rheol. Acta, 6, 119 (1967).

2.34. Ting T. С. T., Remarks

on Linear

Viscoelastic Stress

Analysis

in Cylinder Problems, in «Developments in Theoretical and Ap­

plied Mechanics» (Huang

T. C. and

Johnson M. W.,

Jr., eds.),

New York, 1965.

 

 

 

2.35.Ting T. С. T., Contact Problems in the Linear Theory of Visco­ elasticity, J. Appl. Mech., 35, 248 (1968).

2.36.Ting T. С. T., A Mixed Boundary Value Problem in Viscoelasti­

city with Time-Dependent

Boundary Regions, in «Developments

in Mechanics» (Weis H. J., Young D. F., Riley W. F., Rogge T. R.,

eds.), vol. 5, p. 591, Ames,

1969.

2.37.Tsien H. S., A Generalisation of Alfrey’s Theorem for Viscoelas­ tic Media, Quart. Appl. Math., 8, 104 (1950).

Глава 3

ТЕРМОВЯЗКОУПРУГОСТЬ

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

§ 3.1. Термодинамический вывод определяющих соотношений

В этом разделе мы займемся выводом зависимостей линейной теории термовязкоупругости, основанным на двух фундаментальных термодинамических постулатах: законе сохранения энергии и неравенстве для роста энтропии. Эти два постулата вместе с двумя другими определяющими допущениями приводят к построению полной теории термовязкоупругости. Вывод, который здесь приведен, взят из работы Кристенсена и Нахди [3.3]. Другие способы вывода, также использующие термодинамические соображения, но отличные от приво­ димого здесь, дали Био [3.1], Эринген [3.6], Хантер

[3.10] и Шейпери [3.17].

§ 3.1. Термодинамический вывод определяющих соотношений 113

Установим локальные соотношения, выражающие со­ хранение энергии и неравенство для роста энтропии. Вывод этих локальных соотношений из общих глобаль­ ных дан в работе Трусделла и Тупена [3.20]. При запи­ си названных соотношений мы введем несколько новых полей переменных, включая температуру и энтропию. Хотя в рамках данной работы невозможно дать точные определения этим новым величинам, они нам знакомы в том смысле, что являются обобщениями на случай термо­ динамики необратимых процессов соответствующих из­ вестных понятий теории состояний теплового равновесия.

Локальное условие сохранения энергии для инфини­ тезимальной теории дается уравнением

рг — р [Л + TS + Т 5 ]

+ а

ц е{. QtЛ — 0,

 

(3.1)

где р — плотность, г — функция

притока тепла

на

еди­

ницу массы, А — свободная

энергия Гельмгольца,

отне­

сенная к единице массы,

Т •— абсолютная температура,

S — энтропия, отнесенная

к

единице массы, и

Q, — де­

картовы компоненты вектора теплового потока, отнесен­ ного к единице площади в единицу времени. Функция г в (3.1) дает способ притока или отвода тепла извне. В (3.1) так же, как и в последующем выводе, t обозначает текущее время, если не оговорено иное, а точка сверху обозначает дифференцирование по времени. Форма (3.1) закона сохранения энергии удобна для инфините­ зимальной теории, поскольку включенные в нее напряже­ ние и деформация определяются в том же смысле, кото­ рый они имеют в этой теории. Кроме того, поскольку в силу инфинитезимальной теории плотность с точно­ стью до членов первого порядка остается постоянной, сво­ бодная энергия, отнесенная к единице объема, рА, мо­ жет рассматриваться без связи с плотностью и свобод­ ной энергией, отнесенной к единице массы. Цель вклю­ чения плотности в форму (3.1) и в последующий вывод состоит в том, чтобы добиться общности обозначений с нелинейной теорией, которая будет развита в гл. 6. Ус­ ловие сохранения энергии часто записывается в форме, содержащей внутреннюю энергию вместо свободной энергии. Зависимость между этими двумя величинами

8-851

114

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

устанавливается

преобразованием

Лежандра рU—

= рЛ+ 7р5, где

U — внутренняя

энергия, отнесенная

к единице массы.

 

 

Соответствующее локальное неравенство для роста

энтропии имеет вид

 

p T S - p r ± Q u — Qt [T't/ T ) > 0

(3.2)

ичасто называется неравенством Клаузиуса Дюгема.

Мы покажем, что если постулировать определяющие

соотношения для рА и Qi в соответствующей форме, то в процессе удовлетворения условиям (3.1) и (3.2) мож­ но вывести все прочие определяющие соотношения. Та­ ким образом, отправной точкой является получение вида выражения для свободной энергии рА.

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

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

мы приближения следующим образом.

 

Допустим, что функции Eij{t) и T(t)

непрерывны на

интервале — о о < ^ < о о , и, кроме того,

допустим, что

Eij(t) стремится к нулю и T(t) стремится к Т0 при t-)— оо.

§ 3.1. Термодинамический вывод определяющих соотношений 115

При таких предположениях о непрерывности из тео­ ремы Стоуна — Вейерштрасса следует, что действи­ тельный непрерывный скалярный или тензорный функ­ ционал от еij(r) и Г(т) при —о о < т ^ ^ можно равно­ мерно приблизить полиномом на множестве действи­ тельных непрерывных линейных функционалов от £ц(т) и Т(х). Далее, используя теорему представления Рисса, можно выразить эти линейные функционалы через инте­ гралы Стильтьеса, в которых функции интегрирования имеют ограниченную вариацию. Следуя процедуре из § 1.2, можно записать эти интегралы Стильтьеса в виде (1.5) и в форме, включающей историю инфинитезималь­ ного изменения температуры, начиная с исходной тем­ пературы Tq. Далее мы будем обозначать через Ѳ(?)

разность между текущим значением температуры и

Т0

и предполагать, что ец(%)

и в(х)ІТ0 являются

малыми

величинами порядка 0 ( e ) ,

как это делалось

в §

1.2.

Полиномиальное разложение функции рА относительно этих линейных функционалов дает нам форму

рА = рЛ0 + і Dtj (t — т) дв‘-СФ- d%

{ ß (t — т)

d i +

J

дт

.1

дт

— 00

 

— oo

 

+Т І

----OO — 00

 

-

j

j q > , / ( ^ - T ,

/^ -^TMi d) x^ d r\ —

 

 

----00

— oo

 

---- Г

i

i

m (t — ^ ,t — ^ ) d- ^ - ~ ^ - d x d x \ + 0 ( e z), (3.3)

2

J

J

дт

ÖT]

где Aq— средняя свободная энергия, а функции интег­ рирования, характеризующие механические свойства, полагаются непрерывными по аргументам т ,-^0 и тож­ дественно равными нулю при Т г<0, т. е.

ß(Ti) = 0, Dij(ті) = 0,

Gum(ті, т2)

= О,

 

фи (тъ т2) = 0 , m (ть т2) =

0 при <

0, т2 < 0.

ѵ

8*

116

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

Для рассматриваемой теории членами порядка О (в3) в (3.3) можно пренебречь и функции интегрирования в (3.3) не должны зависеть от деформаций и температуры.

Из (3.1) и (3.2) можно исключить функцию притока тепла г, что дает

—р5Ѳ — рЛ + оцеп Qi(Q,i/T0) ^ 0,

(3.5)

где использована зависимость Г = Г 0+Ѳ и в выражении Т'і/Т сохранены только члены порядка 0 (e ).

Подставим значение рЛ из (3.3) в (3.5) и выполним дифференцирование по t, используя правило Лейбница. Эта операция дает

{ - О ц (0) -

[ Gm (t - х, 0)

 

dx +

 

— COoo —

 

 

 

t

 

 

+

j Фii (0,t — x) d- ^ ~

dx + а и J si} (t) +

 

t

 

 

+

{ß(Ö) + ^ m ( t - x , 0 ) d- ^ - d x +

 

t

 

 

+

j Фі/ (t — T, 0)

 

dt — p5 j è (t) +

+

д е д (т)

dx + j ~ ß (t t )

 

dx

 

 

 

 

 

(3.6)

где

t

 

 

t

 

 

A— І И

S '

 

 

t t

дт 0г]

— OO — oo