книги из ГПНТБ / Постников В.С. Внутреннее трение в металлах
.pdfСледуя способу Фурье [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Т/ |
ГуСОбу |
|
|
|
|
|
|
|
|
т,= - |
іум ту |
8У |
(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