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

Курс механики сплошных сред

..pdf
Скачиваний:
20
Добавлен:
15.11.2022
Размер:
20.58 Mб
Скачать

в) Функции запаздывания. Так же просто определяются для та­ кой среды функции запаздывания или функции вязкоупругой пол­ зучести.

Положим:

М* = (3/(*)-\ ./* = (2р*)“\ С*= (£•)-*.

(111)

Ассоциированные' функции М (t), J (t) и C(t) являются функциями запаздывания при сжатии, при сдвиге и при чистом растя­ жении. Например, соотношение (101) примет такой вид:

eu = J*{D)sif.

(112)

Предполагая, что J* (D) можно разложить на простые элементы

J* (D) = J 0 -)- £

ip

J _

у

,

Q»D

 

1+6,0

 

р - 1

 

f?i

P l+ QpD>

имеем

 

 

 

 

 

 

п

 

 

 

 

'

 

/ > хР 1 -

 

II

+

яъ*

Р=1

i

 

8

о

IW

 

 

 

 

 

и в развернутом виде

"-^ -(Х )5„(г-ь)с1\.

(ИЗ)

(114)

Эта формула обобщает

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

указанные ранее в VIII.7.2

для

случая

трехпараметрической модели. Для того чтобы функция

J (/)

была

возрастающей,

предположим,

что величины J 0 jp, Qp—

неотрицательны. Параметры Qp—времена запаздывания (или вре­

мена ползучести)

при сдвиге;

J 0 и Jx,

равные соответственно

(2р0)_1 и (2|i„)_1,

называются

мгновенной

и длительной податли­

востью при сдвиге.

 

 

 

Очевидно, аналогичные рассуждения могут быть проведены и для других функций запаздывания M(t) и С (t). В последующем придется пользоваться функцией запаздывания N (t), ассоцииро­ ванной N* (D) оператору, соответствующему упругой константе v £ -‘, входящей в формулировки законов поведения. С учетом фор­ мулы (VI,33) имеем простую зависимость

#*=■ (-^-),== у (С* - М ‘) - J*- С*.

Х.3.3. Чистое растяжение цилиндрического стержня. Рассмотрим снова эту очень простую задачу, изученную в X. 1.4. Среду на этот

раз считаем вязкоупругой, а величину F — зависящей (возможно) от времени. Очевидно, что формулы

083 — F,

<Jji = (Т2« — — стаа

^31 — 0|

(115)

е83 “ C*F,

= е22 =

X F,

е12 ^

е22 е21 =ж=О,

X a = xaC*F, X ^

— xfl*?,

X, = — xtN*F

,

 

определяют,

с одной стороны, поле

напря­

 

жений,

удовлетворяющее уравнениям

равно­

 

весия и граничным условиям, а с другой —

 

поле деформаций и поле

перемещений,

по­

 

лученные из данного поля напряжений

через

 

законы поведения. Таким

образом, формулы

 

(115) дают решение задачи единственное с

 

точностью до движения

тела

как

абсолютно

 

твердого.

 

 

 

 

 

 

 

 

 

 

Применим эти результаты для различ­

 

ных видов нагрузки F (<); наиболее важными

 

будут

следующие зависимости:

 

 

 

 

 

 

 

а=‘Е*г,

е =

С*ог,

 

 

 

 

 

в которых

введены

обозначения:

<х3, = а,

 

е88 = е;

а представляет

собой

равномерно

Рис. 15.

распределенную

поверхностную

нагрузку?,

приложенную к торцу 21( е —относительное

а-деформация при постоян­

осевое

удлинение.

 

 

 

 

 

 

 

ной скорости; б-изменение

 

 

 

Испытание при

тельной в начале координат;

а)

 

 

 

 

 

 

напряжения; наклон Е 0 каса­

 

 

 

 

 

 

 

 

 

 

 

наклон Е оо асимптоты

деформации. Предположим,

что в опыте от­

 

носительное удлинение е (/) равно

нулю

при

/ < 0 и пропорционально времени

t

при t

> 0,

иными

словами,

е = at.

 

 

 

 

 

 

 

 

 

нахо­

Используя зависимость (109), интегрированием по частям

дим, что для положительных значений

t

 

 

 

 

 

 

 

o{t) = a \\E ( l) dk.

Если измерить o(t), можно получить интеграл от модуля £ (/); результаты приведены на рис. 15.

Если положить а =‘0 , то будем иметь

в((0) = г0 ^° Е(Х) dX.

График функции Е (t) представляет собой кривую, обращенную выпуклостью вниз, поэтому напряжение о (/0), соответствующее данной деформации е0, достиг­ нутой при постоянной скорости деформации за время t0, будет тем меньшим, чем большим будет интервал времени t0 (или чем меньше будет скорость деформации). Точно так же можно показать, что соответствующая заданному напряжению о0 деформация за время t0 при условии постоянства скорости увеличения нагрузки

^а = о0 -у-^ будет тем большей, чем больше отрезок времени t0.

б) Постоянная нагрузка, приложенная в течение конечного от-

резка времени. Итак, по предположению, нагрузка o(t) равна нулю везде, кроме интервала 0 < i < /0, в течение которого а(/) = а0. Тогда е = а0С(/) при 0 < t < t01 а при t < t 0

е=а°S L . w d x = a °

как это видно из формулы (114) и определения функции запаздыва­ ния C(t). Результаты вычислений представлены На рис. 16.

в) Колебания. Комплексный модуль. По определению, предпо­ лагаем, что е(/) = 0 при t < 0, а при t > О

е (/) = е0 sin со/ «=* Im {е0 ехр (/со/)}, (116)

где введено обозначение Im { } для мнимой части.

Для простоты воспользуемся трехпараметрической моделью, Е* при этом определим по формуле (ПО).

Из соотношения (109) имеем

0 (0 = е0 Im { £ 0 ехр {to*}—

Рис. 16.

а-постоянное напряжение в течение ко­ нечного отрезка времени; б-деформация

exp {to (* ) — £ } dA,} .

Интеграл в правой части вычисляется без труда:

-

Ы Ч } - е > Р <«>)■

Таким образом, имеем следующий результат:

0 (О = ео I m { ^ ^ i e x p ( t o / } + ^ = ^ e x p { - l } } . (117)

Вторая из фигурирующих здесь экспонент со временем убывает и описывает переходной процесс; первая экспонента —синусоидальная функция времени. Через некоторое достаточно долгое время с мо­ мента начала возбуждения значение ст(*) становится практически

равным периодической функции сг(*), если пренебречь переходным процессом. Учитывая выражение (110), можно написать, что

о (*) = е0 Im * (to) ехр (to*)},

(118)

Е* (to) называется комплексным модулем при простом растяжении.

Положим

Е* (to) = E i (со) + iE2(со) = А (со) ехр {*<5 (со)},

отделив действительную и мнимую части, модуль и аргумент вели­ чины £*(to). Теперь имеем в явном виде зависимость

о (*) = Е 0А (со) sin (со* + 6 (со)}.

(119)

Итак, А (со) — коэффициент усиления амплитуды,

б (со) — разность

фаз, что видно из сравнения найденного решения с формулой (116).

 

Для

трехпараметрического

 

модуля

имеем

 

 

Ew-{- Е0(й2т;а

 

 

 

1+<Ла

 

 

_ (£0 —

 

( 120)

 

1 + ( 0 2Та

 

 

 

Таким

образом,

и £ а

 

положительны и,

следова­

Рис. 17. Изменение тангенса сдвига фаз 6

тельно,

б

представляет со­

в зависимости от безразмерной частоты ота

бой угол,

содержащийся ме­

жду 0 и п/2 (при со—>0 и со—►оо угол б бесконечно мал, т. е. для коротких и длинных волн). График tg б представлен на рис. 17.

Заключение. Для реализации периодических колебаний (116) в вязкоупругой среде (тЛ^=0) при пренебрежении переходными процессами требуется приложить периодическое возбуждение (119). При этом деформации отстают по фазе от усилий (сдвиг фаз зависит от частоты возбуждения) и наблюдается усиление амплитуды (также

зависящее от частоты).

так как для

Такое явление в упругих средах не наблюдается,

них Ei) = Eoo= Ei

та = б = 0. Принято говорить, что

вязкоупругие

волны являются

диспергирующими.

 

Рассуждения проведены на примере трехпараметрической модели, однако формулы (118) и (119) и качественные выводы остаются

справедливыми и для более

общих

моделей, для

которых

модуль

релаксации

Е (t) описывается выражениями (107)

и (108).

 

Опыты с

вибрацией легко осуществить и получить тем

самым практический

метод определения вязкоупругих

свойств среды, комплексного модуля

Е* (ио) и

затем в итоге

модуля релаксации

Е (/).

 

 

 

Напомним, однако, что в связи с тем что использована квазистатическая

гипотеза, приведенные результаты не будут

верными для очень высоких частот.

Заметим,

наконец, что если представляется возможным

проводить

измерения

в переходном

режиме, то можно определить

время релаксации та.

 

Х.3.4. Другие задачи теории вязкоупругости. Можно вновь рас­ смотреть задачи, связанные с кручением и изгибом цилиндрического стержня или о раздувании сферического резервуара и решить их приведенным выше методом. Заменяя в решении упругие константы соответствующими операторами, получаем решение (в операторных обозначениях) вязкоупругой задачи. Справедливость полученного решения подтверждается тем, что уравнения равновесия, законы поведения и граничные условия выполняются. Для получения реше­ ния в явном виде достаточно применить правила операционного исчисления, установленные в Х.3.1. Этот метод решения называют

принципом соответствия.

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

В задаче об изгибе балки операторная запись формулы, связы­ вающей изгибающий момент (относительно главной оси инерции) с кривизной дуги, в соответствии с решением (26) такова:

 

M W ^ I E U D ) - ^ ,

= \ c * ( D ) M ( t ),

где

/ — момент инерции поперечного

сечения относительно главной

оси,

коллинеарной приложенному моменту.

Взадаче кручения крутящий момент М (t) и угол закручивания

всоответствии с формулой (51) связаны соотношениями:

М (/) = ф* (D) а (0,

а (0 = 4

J* (D) М (/),

где d константа, определяемая

по формуле (51).

(114) имеем

Если, например, дан момент М (t),

то

согласно

da (t) = 2J (0) М (t) + 2 J "

(k) M ( t - k )

dk,

где J (/) функция запаздывания

при

простом сдвиге.

В сферическом резервуаре, у которого внешняя поверхность сво*

бодна от нагрузок, а внутренняя

испытывает давление

p(t), поле

напряжений снова

дается выражениями (37), а операторная запись

поля перемещений

такова:

 

 

&(0 = - й£ ту( р

+

(121)

Очевидно, что во всех этих задачах можно рассмотреть частные формы нагрузок, приведенные в предыдущем параграфе. В случае гармонического возбуждения, например, если пренебречь переходным процессом, получаем решение также в гармонической форме, в ко­ тором коэффициент усиления и сдвиг по фазе зависят от частоты и определяются соответствующим комплексным модулем £*(/со), р*(ко), /С* (/со), или комбинацией этих модулей. Для случая, опи­ сываемого формулой (121), если при / > 0

Р (0 = Ро sin со/, гармоническая реакция дается формулой

^ (0 = R3°—r* Irn { (зК*(ito)"f“ 4(1* (ICO) р2 ) ехр

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

чистом растяжении заданное

значение напряжений сг0 достигается

в момент tQ за время Т, то

деформация в момент t0 будет тем

 

 

 

 

 

 

 

больше,

чем

больше

промежуток

време­

 

 

 

 

 

 

 

ни Т, т. е. чем медленнее

растет

нагрузка.

 

 

 

 

 

 

 

 

Эти

явления—характерные

признаки

 

 

 

 

 

 

 

того,

что

среда

является

вязкоупругой.

 

 

 

 

 

 

 

В частности,

металлы

при высоких

тем­

 

 

 

 

 

 

 

пературах,

бетон,

пластмассы

обладают

 

 

 

 

 

 

 

именно

такими

особенностями,

т. е. явля­

 

 

 

 

 

 

 

ются

вязкоупругими

средами.

 

 

 

ве­

 

 

 

 

 

 

 

 

Проведенные

выше

исследования

 

 

 

 

 

 

 

 

лись

в

 

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

что

поведение

 

 

 

 

 

 

 

среды

было

линейным.

С помощью тео­

 

 

 

 

 

 

 

рии

вязкоупругости

хорошо описывается

 

 

 

 

 

 

 

поведение

некоторых

пластических

мате­

 

 

 

 

 

 

 

риалов, состоящих

 

из макромолекул

(на­

 

 

 

 

 

 

 

пример,

поликарбютен,

полиизобутилен,

 

, ,

,

, ,

,

,

, ,

политен,

полиуретан

и

т. п.).

Эти

мате-

 

риалы

 

используются,

в частности,

для

7-6 -5 -1

-3 -2

- 1

0 1

2 3 1 log i

твердого

топлива

ракетных

двигателей,

и

 

 

 

 

 

 

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

Рис.

18.

Модули

релаксации

чительной

степени

обязана

своим

разви­

полиизобутилена

 

при рас­

тием

необходимости

знать

механическое

ширении (К) и при сдвиге (р):

обе

шкалы—логарифмические;

поведение шашек твердого топлива при пе

 

<-с;

К и

д -П а

ременных нагрузках

 

во

время горения

и

 

 

 

 

 

 

 

при внешних

воздействиях,

возникающих

во время работы двигателей на марше.

 

описанные в Х.3.3,

могут

Простые

испытания,

как,

например,

помочь в поисках наиболее простой формы выражения закона пове­ дения среды, при которой число неизвестных параметров минимально (например, модель с тремя параметрами или еще более частный случай—тело Максвелла или тело Кельвина—Фойхта). Более Общие формулы, как, например (108), с большим числом параметров также могут быть полезными.

На рис. 18 приведены для наглядности некоторые эксперимен­ тальные значения модулей релаксации при расширении и при сдвиге для полиизобутилена *.

Х.4. ЗАКЛЮЧЕНИЕ

В заключении предлагаем читателю изучить рио. 19, в котором сведены для случая разных сред (упругой, упруго-идеально-пласти-

ческой, вязкоупругой) графики кривых

изменения напряжения а

при чистом растяжении под действием

деформации, приложенной

по определенному закону [относительное удлинение к тождественно

равно нулю, кроме интервала (0, (0),

когда е = е0].

напряжение а

При этом предполагается, что

в

момент

t = 0

принимает одно и то же значение

сг0.

Кроме

того,

предполагается,

* См. работу Хуана, Ли и Роджерса tStanford University Technical Report», 1963, n° 140.

что для упруго-идеально-пла­

 

 

 

стической среды напряжение &

 

 

а0 достаточно велико и по­

 

 

 

зволяет достичь

предела

те- е0

 

 

кучести; в

последующие мо­

 

 

 

менты

(t >

tо)

в

материале

 

 

 

имеется

постоянное по вели­

 

 

 

чине сжатие. В

случае вяз­

 

 

 

коупругой среды

также

на­

Рис. 19. Сравнение эволюции напряжений

блюдается сжатие

в моменты

для различных материалов для одной и той

t > /0,

но оно постепенно па­

же предыстории деформаций:

дает до нуля. В упругой сре­

-------------- упругая ср еда;

упруго-идеаль­

но-пластическая

среда;

— вязкоупругая

де реакция а всегда пропор­

 

среда

 

циональна возбуждению е, и в

 

 

 

моменты t >

t0материал мгновенно, а не в конечное время, как в слу­

чае вязкоупругой

среды,

возвращается в

начальное

естественное

состояние. Пластические среды в естественное состояние не возвра­ щаются.

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

 

ПРИЛОЖЕНИЯ

 

 

 

 

 

 

ПРИЛОЖЕНИЕ I

 

ВЕКТОРЫ И ТЕНЗОРЫ

 

 

 

 

П1.1. АФФИННОЕ ЕВКЛИДОВО ПРОСТРАНСТВО

 

Пространство, которым обычно оперирует классическая меха­

ника,— это

трехмерное аффинное евклидово пространство <£, эле­

менты которого —точки. Паре

упорядоченных

точек

(Р, Q) из $

соответствует элемент X векторного евклидова трехмерного прост­

ранства Е‘

(Ру Q)*X=PQ. Это отображение обладает известными

свойствами

PQ = — QP, PQ = PM + MQ. Если

точку

О из при­

нять за начало координат, то существует биективное

отображение,

ставящее в

соответствие любой

точке М из £

вектор

V из £, для

которого (О, М) V.

пространство Е является векторным

П 1.1.1. Векторы. Линейное

евклидовым

пространством, если каждой паре векторов

Ху Y из Е

можно поставить в соответствие действительное число,

называемое

скалярным произведением вектора X на К и обозначаемое X

Y Отобра­

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

 

(aX + bY) Z = a(X-Z) + b(Y

Z),

 

 

 

 

X

Y = Y Ху

 

 

 

 

 

 

Х Х ^ О .

 

 

 

Неотрицательный скаляр

(X2)1/2f где Х 2= Х Ху

называют моду­

лем вектора

X и

обозначают [ А" | . Модуль

равен

нулю

тогда и

только тогда,

когда

X — нулевой вектор из пространства

Е. Если

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

нулю.

 

 

 

 

 

 

 

 

означает,

 

Утверждение о том, что Е —трехмерное пространство,

что существует базис (или бесконечное число базисов):

базис обра­

зуется

из трех векторов

е19 e 2f е8, для

которых

всякий

вектор X

из

Е

будет

однозначно

представим

в

виде линейной

комбинации

e lt

е 2у

е3. Известно, что в пространстве Е может существовать бес­

конечное число ортонормированных

базисов.

Базис е 19 e 2f е3

или

e;(i= 1, 2, 3)

будет ортонормированным тогда и только тогда, когда

 

 

 

в гв / = б,/ (/, У= 1,2,3),

 

 

 

 

(1)

где

 

—символы Кронекера. Символ Кронекера равен

1,

если

оба

индекса совпадают (6n = 6j2 = 63s = 1), и

равен

0,

если

индексы не

равны

Цф})

(6,g = 6t, =

= 6 3J = 0).

Иными

словами,

величины

б/у-

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

книге,

если

не

оговорено противное,

будут использоваться только

ортонорми-

рованные базисы.

 

 

 

 

Выбрав в пространстве Е ортонормированный базис, можно

представить любой вектор

X в виде

 

 

 

X = Х ге^ + Х 2е2+ Х 8е8 = Х (е(.

 

(2)

Величины X i—составляющие вектора X в базисе et.

Здесь

будем

использовать соглашение о суммировании по повторяющемуся

(немому) индексу. В соответствии

с этим соглашением по индексу,

который повторяется два раза в каком-либо

одночлене,

должно

быть произведено

суммирование,

иными

словами, одночлен должен

быть заменен суммой всех членов,

которая получается, когда индексу

придают значения

1, 2, 3.

Индекс i в

(2) называется немым, так

как буква, которой он обозначен, не

имеет,

очевидно,

никакого

значения: Х {е( и Xkek обозначают один

и тот же вектор X. С уче­

том (1) и в соответствии с соглашением о суммировании

 

X

У= Х ^ + Х 2У2+ Х3У3 =

 

(3)

Линейная форма 2 (Л),

значения которой —скалярные

величины,

есть линейное отображение пространства Е на множество дейст­

вительных чисел. Таким образом,

 

* (X) = х гш (ег) + Х 2% (е2) + Х32 (е8).

(4)

Линейные формы являются элементами трехмерного

векторного

пространства Е, называемого сопряженным пространству Е\ формы

А\, Х2, Х8 составляют базис этого

пространства,

называемый сопря­

женным

базису

(е() пространства

Е. Величины

2 (е()—координаты

формы 2

в этом

базисе. Если задать вектор А из Е его составля­

ющими Л, = 2(е/) в базисе

et пространства

£, то будем иметь

 

 

%(Х) = А - Х= А;Х(.

 

(5)

Отображение

2 —►Л,

которое

ставит в

соответствие форме 2

вектор А из пространства £, реализует изоморфизм пространства Е на пространство Е и позволяет, таким образом, отождествить сопря­

женное пространство Ё с самим пространством Е. Ниже будем это делать систематически. Практические упрощения, которые из этого следуют, в конечном счете оправдывают появляющиеся при этом теоретические сложности. Так, вектор Л из £ можно рассматри­ вать либо непосредственно как один из элементов £, либо как операцию, определяющую линейную форму:

2(Х) = А-Х.

В ортонормированном базисе е( имеем одновременно

 

А = А{еь %{ X) ^AiXi.

(6)

Замечание. Тот факт, что здесь используются только ортонормированные базисы, позволяет оперировать только с нижними индексами и применять правило суммирования по немым индексам в нижнем положении. Пусть, например, базис

в{ не ортонормированный, тогда составляющие вектора X будут

обозначаться X1

н уравнение (2) перепишется в такой форме:

 

X = Xieh

(7)

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

внизу.

Равенство (1) перепишется следующим образом:

er ej = 8i/>

(8)

где gij— элементы симметричной матрицы

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

и называемой фундаментальной матрицей базиса в[.

 

Скалярное произведение (3) запишется в виде

 

X - Y = g uXiYJ,

(9)

правая часть в силу соглашения о суммировании по повторяющимся (немым) индексам представляет собой сумму:

gliX m + g i2XiY*+g21X*Yi+...

Квадратичная

форма gijXlYl положительно определена, а определитель мат­

рицы g/y, det (g/y)

отличен от

нуля.

Линейная форма % (X) на

пространстве Е в базисе Х'Е имеет выражение

*(*) = *'■*(*/)•

Существование скалярного

произведения позволяет

поставить

в

соответствие

всякому вектору А линейную форму:

 

 

 

 

 

 

 

 

 

A-X = gijA‘X J = gljA!Xi = AiXi,

 

 

 

 

(10)

если положить

 

 

Ai = giJAi.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(11)

Таким образом, можно снова

привести в соответствие £ и £ в

силу

изомор­

физма g — >-Л,

но при заданных

g составляющие

Ai

получаются

из

решения

уравнений

(11)

приi4/ = g (e/),

что всегда

возможно,

так

как det (g/y) Ф 0.

Если

обозначить

через

g'V матрицу,

обратную

g/y, решение запишется

в таком

виде:

 

 

 

 

A* = gVAj.

 

 

 

 

 

(12)

Таким образом, вектор А пространства Е в случае неортонормированного

базиса е/ имеет два разложения

на составляющие— составляющие А

называемые

контрвариантными, и составляющие Л/, называемые ковариантными, которые следует различать. Формулы (6) запишутся тогда в виде

A = Aieh g (X) = А(Хг\

(13)

Можно к тому же убедиться, что Л,- действительно являются составляющими

вектора А в базисе, определяемом векторами

 

ei^g'J'ej

(14)

и называемом взаимным базису е(.

 

А — А ^ ^ А ^ ,

(15)

Можно теперь оценить упрощения, вытекающие из рассмотрения только ортонормированных базисов: g/y, gU — единичные матрицы, базисы е,- и ei тож­ дественны, как и составляющие Л* и Л/ одного и того же вектора Л. Положение

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

П1.1.2. Использование соглашения о немом индексе. Следует сделать несколько замечаний об употреблении немого индекса,