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

книги из ГПНТБ / Постников В.С. Внутреннее трение в металлах

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.97 Mб
Скачать

Следуя способу Фурье [221], найдем решение

уравнения (346)

в виде

 

СО

(347)

W= Ъ Xn{z)Yn{t),

Л= 1

 

где Хп (г) — функция только координаты z, а

У„ (t) — функция

только времени t. При этом каждый член ряда должен удовлетво­ рять граничным условиям задачи (по два на каждом конце стержня/,

Рис. 70. К задаче об нзгнбных колебаниях стержня с распределенными парамет­ рами (схема):

а — деформации стержня; 6 деформации элемента стержня; в — распределение напряжений в элементе

а вся сумма — также начальным условиям. Подставляя выражение (347) в уравнение (346) и требуя, чтобы равенство удовлетворялось для любого номера я, получим

£7

Уп

(348)

ps

YnJr~]TYn

 

 

Здесь точка обозначает дифференцирование по времени, а штрих — по координате.

Так как левая часть уравнения (348) не зависит от времени t, а правая — от координаты г, то для тождественного выполнения равенства необходимо, чтобы обе части равнялись одной и той же постоянной.

Обозначая ту. постоянную через

а>1, приходим к

уравнениям

Х2-~Ъ\Хп =

0,

(349)

Yn+T\-4-Yn+<£Yn = Q-

(350)

130

В уравнении (349) обозначено

л4 _ Ps

(351)

Оп— £Y ®/г-

Интеграл уравнения (349) записывается в форме

*n(z) =

•ln(2) + -B,Än (z) + CnZ3n (z)

(352)

Выражения

 

21

(ch bnZn-f- cos bnz):

 

 

Z\,i —

 

 

Z,n =

21

(sh bnzn+ sin bnz)

(353)

 

Z$n ~

21

(ch bnz ■— cos M ;

 

 

 

Z^n =

21

(sh bnz -—sin bnz)

 

представляют собой функции Крылова, которые связаны между собой дифференциальными соотношениями

Z\п = bn

Zo,n

 

(354)

Z 3 п =

 

Соотношения между постоянными Ап, Вп, Сп, Dn и собственными частотами со„ определяют граничными условиями. Для нашего слу­ чая на закрепленном конце равен нулю прогиб и угол поворота w'\ таким образом, граничные условия для левого конца имеют вид

Хл (0) = 0; Х(0) = 0.

(355)

Свободному концу стержня соответствует отсутствие огибающего момента, пропорционального w" и поперечной силе, пропорцио­ нальной w'". Это приводит к граничным условиям

*„(/) = 0; Х (0 = 0.

(356)

С помощью уравнений (355) и (356) легко найти постоянные Л„, Вп, Сп, Dn и собственные частоты юп, а следовательно, установить вид фундаментальных функций Хп (z). Важно отметить, что фунда­ ментальные функции каждой задачи совершенно не зависят от вязких свойств стержня и образуют ортогональную систему, т. е.

JIХт (г) Хп (z)dz = 0 ,

(357)

о

 

если т Ф п.

 

Возвратимся теперь к уравнению (350), описывающему, оче­

видно, процесс'затухающих колебаний.

 

Его интеграл имеет вид

 

Yn= е Ѵп‘ (Сп sin соnt -f- Dncos <Хпі).

(358)

9*

131

Здесь

 

со=

 

 

 

 

Уп = tio> « / 2 Е,

(359)

Таким

 

j/^l —(I F )2-

(36°)

 

СО*

Ѵ'а>1— уІ ~ (-0,1

 

 

образом, полное решение уравнения задачи

имеет вид

 

w

2 Хп(г)е~~Упі (Спsin co,V-f Dncos со,*0-

(361)

 

 

Л—1

 

 

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

 

 

ш (г,

0) =

h (г);

 

= /а(г),

(362)

где / х (z) и / 2

(г) — известные функции. Тогда при

t = 0, согласно

уравнениям (361)

и (362), имеем

 

 

 

 

 

оо

Dnx n(z);

 

 

со

 

 

 

Ы г) =

X

 

/2 (2 ) =

2, (Cnan — D„yn) x n(z).

 

n = l

 

 

 

 

 

n ~ l

 

 

 

Умножая обе части этих равенств на

Хп (г) и интегрируя в пре­

делах всей длины стержня, получим

 

(z)dz,

 

(C„(ä'n —

t

 

 

1

 

 

 

J

Xl (z)

dz = J

/1 (z)Xn

 

 

 

Dn

 

 

 

 

 

 

о

I

 

0

I

 

 

(363)

 

 

 

 

 

 

 

 

 

 

D„yn)

X2n

 

/ 2 (2 ) Xn

 

Соответственно условию ортогональности фундаментальных функ­

ций (357) все остальные слагаемые,J (z)dz ==Jвходящие в(z)левыеdz.

части этих

равенств, обращаются в нули. Теперь из уравнения (363) нетрудно найти Сп и Dn для любого номера п.

Рассматривая выражение (361), заметим, что каждый его член

описывает затухающие колебания, если — действительное число. Из уравнения (360) видно, что это характерно лишь для нескольких начальных значений п, пока выполняется неравенство

2 Е ^

1

(364)

 

 

 

При достаточно больших значениях п неравенство (364) нару­

шается и величина а>1 становится мнимой. При этом соответствующие члены общего решения (361) уже не будут описывать затухающих колебаний, но будут представлять апериодическое движение стержня. Имея это в виду, можно, как это делали выше, разделить решение (361) на две части: затухающую и апериодическую.

132

Заметим еще, что искусственным выбором начальных условий можно обеспечить монохроматичность процесса колебаний. Пусть, например, f x (2) = Х х (2), f23 (2, = 0), т. е. в начальный момент за­ даны прогибы, точно соответствующие первой фундаментальной функции.

Из уравнения (363) получим

Сх = у]/«:;

С2 = С3 —■• •= 0;

D 1 = 1; D 2 = D 3 = . . .= 0,

 

 

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

колебания 'дол = Х г (2) Ух (2); это

значит, что

заданная в начале

форма f x (2) = Х х (2) устойчива и без всяких искажений сохраняется

в течение всего

колебательного процесса, лишь масштаб ее непре­

рывно меняется

с течением времени.

 

Согласно выражению (359), логарифмический декремент

 

 

6п

Е

(365)

т. е. совпадает с выражением (66), установленным при исследовании модели Фойгта. Следовательно, все, что было сказано относительно формулы (66), остается в силе и для формулы (365).

Исследование поперечных колебаний стержней, материал ко­ торых описывается другими реологическими моделями (Максвелла, Зинера. . .), представляет значительно оолее трудную задачу [222].

Вынужденные колебания

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

F (г, t) = F0 (г) sin cof.

(366)

Для получения уравнения вынужденных колебаний введем на­ грузку (366) в уравнение (345). Тогда вместо (346) придем к урав­ нению

Ш_

.

ці

d*w .

d*w

_

F0(г)

(367)

ps ' dz*

~r

ps

' dz* dt '1_

dt*

~

ps

 

 

 

Далее поступим, как в разделе «Свободные колебания», и будем искать решение разложением в ряд по фундаментальным функциям соответствующей однородной задачи:

w(z, f ) = S Xft(z)K„(f). ,

(368)

n=l

 

В конечном итоге для фундаментальных функций Х п найдем прежнее уравнение (352), а для функции У„ (t) — уравнение

+ со2 J - Yn+ соlYn =

sin <в/,

(369)

где hn — коэффициент в разложении амплитуды F0 (г) в ряд по фундаментальным функциям:

Л ( 2 ) = | K X n{z). -

(370)

133

Для того чтобы определить hn, умножим обе части равенства (370) на Хп (z) и проинтегрируем его по всей длине I стержня.

Благодаря свойству ортогональности фундаментальных функций (357) в правой части сохранится только одно слагаемое, мы получим

F0 (z) Xß(г) dz = hnJ Хп (z)dz.

Отсюда находим

 

 

 

о

 

 

I

 

 

 

 

 

J

I (z)

 

(г) dz

 

(371)

оJ

 

 

 

 

Fo

X n

 

 

 

 

 

\xl(z)dz

 

 

 

 

о

 

 

 

 

Величину hn определяют

конкретным законом изменения ам­

плитуд нагрузки Fо (г).

 

 

 

 

 

 

Пусть, например, нагрузка задана в виде сосредоточенной силы

F 0 sin оЛ, приложенной

в сечении

z — a

стержня. Тогда

после

предельного перехода найдем из

уравнения

(371)

 

 

hП

 

Х„-(а)

 

(372)

 

JIХ2а (г)<Ь

 

 

 

о

 

 

 

 

Стационарная часть решения уравнения (369) имеет вид

 

Уп — \

Sin (tot

cc,j),

 

(373)

где

 

 

 

 

 

 

 

 

 

 

 

 

(374)

 

 

 

 

 

 

(375)

Легко видеть, что при а>п >> ш знаменатель выражения (374) возрастает с увеличением номера п; соответственно происходит убывание амплитуд Ап. Это значит, что высшие гармоники (при условии ш„ >» ш) имеют тем меньшее значение, чем больше номер гармоники. Если значение частоты со возмущающей силы совпадает со значением одной из собственных частот со„, то Ап приобретает значение

Ап

hnß

(376)

 

pSpco® ’

 

Отсюда следует, что резонансу с высшими собственными часто­ тами, как правило, отвечают меньшие резонансные амплитуды.

134

Нетрудно показать, что коэффициент поглощения при со = со„

(377)

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

Тепловая релаксация Зинера

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

Однако положение меняется, если рассматривать поперечные колебания тонких образцов [165]. В этом случае температурная релаксация происходит в результате потока тепла через поверхность образца. Следуя Зинеру [165], рассмотрим эту задачу.

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

г = Ер'о + \АТ.

(378)

В этом выражении Ер — релаксированный

(изотермический)

модуль Юнга; X— коэффициент линейного теплового расширения; АТ — разность температур между растянутыми и сжатыми волок­ нами; е — деформация, определяемая выражением (341).

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

(дЛдифф = —т 1 АТ.

(379)

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

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

135

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

(ДТ)ад = —у&,

(380)

где

 

Объединяя уравнения (379) и (380), получаем

 

(AT) = — x~1AT — yè

(381)

Уравнение, связывающее напряжениеи деформацию, получается исключением величины АТ из уравнений (378) и (381). При этом находим уравнение (67), где величина М 2 = Ен определяется из выражения

Ен = (1+Ку)Ер.

(382)

Следовательно, если температурныефлуктуацииподчиняются уравнению (381), то термоупрутая связь превращает твердое тело, подчиняющееся закону Гука, в тело Зинера, неупругие свойства которого рассматривались подробно выше (см. гл. Ill и IV).

Из уравнения (382) получаем

Ае = Ху.

(383)

Можно вывести и другую формулу для Д£, выразив ее через величины, которые могут быть более непосредственно измерены, чем у. Если рассматривать а и Т как независимые переменные, а в и s — как их функции, то будем иметь

 

de^ET'do + XdT

(384)

 

ds = X da -f- р dT,

(385)

 

( dS\

(386)

 

 

Симметрия коэффициентов в уравнениях (384) и (385) следует

из того факта,

что изменение плотности термодинамического потен?

цила есть полный дифференциал.

 

Поэтому из

уравнения (см. . например, [166])

 

d (Е га TS) = —е da — 5 dT

(387)

вытекает, что

 

 

 

( - І Н г = (-Ц-)т-

<388>

Полагая в уравнении (385) dS=0 и решая его совместно с выраже­ нием (384), получаем для адиабатического модуля

136

и потому

(389)

Теперь вернемся к уравнениям (378) и (379). Уравнение (379)

устанавливает простую релаксационную связь между АТ и АТ. Заменим ее на более общую, которая дается уравнением диффузии. Тогда вместо (379) напишем

(АТ) = D V2 (АТ).

(390)

Объединяя его с уравнением (380), найдем

(АТ) = D V2 (АТ) ук.

(391)

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

D Ѵ2Х + тХ = 0

(392)

при соответствующих граничных условиях.

В рассматриваемом частном случае деформацию в хорошем при­ ближении можно считать функцией только поперечного размера образца (341). Поэтому уравнение (392) можно заменить уравнением

( D - ^ + x-i)X(y) = 0,

(393)

где координата у соответствует поперечному размеру. Граничные условия должны соответствовать отсутствию теплового потока через границы. Если размер образца.а в направлении у (см. рис. 70), то граничное условие имеет вид

-^- = 0 при y = ± - j - .

(394)

Собственные ортонормированные решения и характеристические времена релаксации для этой задачи равны соответственно

X

.=Ѵ=

(2п +

I) 2

а2

(396)

4

sin

' ( 2 n + l ) - ^ - '

(395)

 

—1

 

 

пЮ

 

Разложим теперь е, а и АТ в ряд по ортонормированным функ­ циям X (у):

е = S в/ ( 0 х і\ ° = Е аі ( 0 Xh

 

/

/

 

 

A T = l i Tj (t)Xj. .

(397)

137

Подставим (397) в уравнения (378) и (391), умножим их на Х;- и проинтегрируем по объему образца, используя основное свойство ортонормированных функций [см. (357)]; получим

е/ = Ep'ojXTj,

(398)

Ті = — і у 1Ті - у г і.

(399)

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

Скорость рассеяния энергии образцом выражается формулой

(А№) = I ое dv,

где интеграл берется по всему объему образца. Используя уравне­ ние (397) и условие ортонормированности, приходим к выражению

W) = £ °у6 /.

(400)

і

 

Если теперь предположить, что е, а и Т изменяются периодически со временем, и если взять среднее по времени, то получим, исполь­ зуя уравнение (398),

=

/

/

= —

(401)

 

/

 

так как 8 у8 у = 0

ввиду периодической зависимости

8 у от времени.

Из уравнения

(399)

 

 

 

или

шТ,- = —ту 1Т/

ГуСОбу

 

 

 

 

 

 

т,= -

іум ту

(402)

 

1 -f- /соту

 

 

 

Умножая уравнение (402) на е;- = гсоеу, усредняя

по времени и

взяв вещественную часть произведения Губу, получим

 

"HF

7 ш2ту

а

(403)

 

8 ;I ; = ----- — Б,.

 

“Г 1

1 + (А?

 

 

В том случае, когда Ер/Еѵ лишь немного меньше единицы, вну­ треннее трение может быть определено как [см. (38) ]

0-і =

2 nW

(404)

ѵ

 

1 Величину Т можно рассматривать как термодинамический потенциал, так как система находится в равновесии только в том случае, если Т повсюду имеет оди­ наковое значение.

138

где

 

 

 

 

 

W = EPj гЧѵ

=

Ер%г),

(405)

а потеря энергии за период

 

 

 

 

 

ДГ = — —

А#.

(406)

 

 

 

 

 

Используя уравнения (403)—(406),

а также (383),

получаем

" 1

 

 

 

 

(407)

Q

 

 

 

 

 

Введем величину

 

 

 

 

 

 

 

 

 

 

(408)

Тогда уравнение (407) можно переписать

 

'

1

= Д£

 

СОТ/

(409)

Q

 

 

Т + Щ ’

 

 

 

 

 

 

По определению величины /у удовлетворяют условию

 

 

Е/у =

1.

(410)

Таким образом, определение внутреннего трения, обусловленного релаксацией теплового потока при изгибных колебаниях, сводится

к

проблеме

нахождения

характеристических

времен

релаксации

Ту

и коэффициентов /у.

времена ре­

й'Чо3

 

 

 

Характеристические

 

 

лаксации для

данной

задачи опре­

 

 

 

деляют

уравнением (396).

 

Коэффи­

 

 

 

циенты /у даются выражением

[165]

 

 

 

 

 

fl

=

96л"4 (2/ +

I)"4.

(411)

 

 

 

 

Первые

коэффициенты

равны:

 

 

 

f о =

0,986;

/j

= 0,012;

/ 2

=

0,0016.

 

 

 

 

В этом случае, если предполо­

 

 

 

жить

коэффициент f о

равным

еди­

 

 

 

нице,

а

все

последующие

коэффи­

 

 

 

циенты

равными нулю,

будет допу­

Рис. 71.

Сравнение

теории Зннера

с опытными данными

Беневитца и

щена

очень

небольшая

 

ошибка.

Ретгера

 

 

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

как простую релаксацию со временем релаксации, равным a/n2D, где а — ширина образца, а D = и/£р — коэффициент тепловой

139

Соседние файлы в папке книги из ГПНТБ