книги из ГПНТБ / Постников В.С. Внутреннее трение в металлах
.pdfУчитывая уравнение (161) из формулы (160), найдем, что высота пика внутреннего трения равна
sin Я|? |
д _ |
^ OD |
|
Qmax —2 + Д ,, cos ф ’ |
Им -- |
Ѵ м<»мо |
(162) |
В случае стандартного линейного тела из выражения (162) сле дует формула (90).
На рис. 44 сопоставляется экспериментальная кривая Q- 1 (Т) алюминия для зернограничной релаксации с теоретическими кри выми, полученными из уравнения (160) при различных значениях
параметра у. |
Как |
видим, |
соответствующим |
выбором |
параметра, |
||||||||||||
а4- ю* |
|
|
|
который можно назвать параметром |
|||||||||||||
|
|
|
размытости |
спектра, |
можно |
лучше |
|||||||||||
|
|
|
|
|
описать экспериментальную |
кривую |
|||||||||||
|
|
|
|
|
Q~1 (Г), чем это позволяют сделать |
||||||||||||
|
|
|
|
|
экспоненциальные ядра (150). Более |
||||||||||||
|
|
|
|
|
подробный анализ преимуществ ядер |
||||||||||||
|
|
|
|
|
Работнова |
перед |
экспоненциальны |
||||||||||
|
|
|
|
|
ми |
|
ядрами |
при |
описании |
Q“ 1 |
(со) |
||||||
|
|
|
|
|
и |
|
Q- 1 |
(Т) |
дан |
в |
работах |
[87, |
|||||
|
|
|
|
|
180, |
181]. |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
Вернемся еще раз к формуле (160). |
||||||||||
|
|
|
|
|
При М о = 0 эта формула принимает |
||||||||||||
|
|
|
|
|
вид |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
0Г1 = ----— • |
|
|
(163) |
|||||
|
|
|
|
|
|
|
|
|
|
|
соэф + |
(оут)ѵ |
|
, . |
|
||
|
|
|
|
|
|
|
Отсюда |
в случае у = |
1 вытекает |
||||||||
|
|
|
|
|
известное выражение |
(60) |
внутрен |
||||||||||
Рис. |
44. Температурная |
зависимость |
него |
трения |
модели |
Максвелла, |
|||||||||||
внутреннего трения поликрнсталлнче- |
широко используемое для |
описания- |
|||||||||||||||
ского алюминия 99,98%), |
ѵ = 0,6 Гц: |
||||||||||||||||
/ — экспериментальная |
кривая; |
2, |
высокотемпературного фона внутрен |
||||||||||||||
3 .— теоретические, |
найденные по фор |
него трения металлов (см. [182—184, |
|||||||||||||||
муле (160) (2 — у |
— 1; 3 — у = |
0,5) |
|||||||||||||||
|
|
|
|
|
87, |
|
с. |
43]). |
|
|
|
|
|
|
|
|
|
|
Существенное различие внутреннего трения, описываемого мо |
||||||||||||||||
делью Максвелла и выражением (163) при у = |
1, состоит в том, что |
||||||||||||||||
при ю —* 0 в первом случае Q- 1 |
—»оо, тогда как во втором Q" 1 —> |
||||||||||||||||
-> tg ф. |
|
выражение |
(163) |
можно |
записать |
следующим |
|||||||||||
1 |
Если (сот)1’ > 1, |
||||||||||||||||
образом: |
|
|
Q 1 *=« (сот)“Ѵsin ф, |
|
|
|
|
|
|
(164) |
|||||||
|
|
|
|
|
|
|
|
|
|
что с точностью до константы совпадает с формулой Q- 1 = (сот)-"'’, предложенной в работе [185] для описания дислокационного фона внутреннего трения металлов и уточненной (в смысле некоторого обоснования ее вывода) в работе [186]. Вопрос о возможности ис пользования формулы (164) по всей области частот (и температур) будет обсужден позже (см. гл. V).
§0
5. Анизотропные тела
Связь между тензором напряжений о.к и тензором деформации г1т для равновесных состояний монокристалла дается законом Гука [166]:
°ік — |
(165) |
где ХіЫт— тензор модулей упругости (изотермических для инфразвуковых частот и адиабатических для ультразвуковых н гиперзву ковых частот). Для описания неупругнх явлений (без учета теплопро водности) необходимо ввести, как мы говорили выше, тензор релакса ции tlk, который характеризует отклонение состояния систем от равновесного значения. Повторяя рассуждения, изложенные ранее, можно получить
( |
|
aik= \klnfiltn + { FіЫт — О г1т(О dt', |
(166) |
— СО |
|
где FMm — ядра релаксации; по дважды встречающимся значениям производится суммирование. Уравнение (166) впервые получено Вольтерра [178] на основе теории линейных функционалов.
На основе термодинамической теории сходные уравнения полу чены в работе [187]. Там же показано, что при экспоненциальных ядрах релаксация компонент тензора напряжений к своим равновес ным значениям не следует простому экспоненциальному убыванию до равновесного значения, как то характерно для изотропной среды, а определяется всей совокупностью времен релаксации хік. Таким образом, в общем случае релаксация напряжений однородной ани зотропной среды определяется шестью временами релаксации хік, соответствующими шести независимым компонентам тензора напря жений.
Уравнение (168) можно обратить, выразив явно тензор напряже ний через тензор деформации:
eife = |
8іЫггРіт4~ |
— О °іт (0 - ■ |
(167) |
Здесь sikltn— тензор податливости, |
связанный с тензором упругих |
||
постряйных Xiklm |
соотношением |
Sfktm = hklm, Фікіт (t— t’) — |
|
■ядра последействия |
или ползучести. |
|
Некоторые конкретные примеры решения задач по рассеянию энергии колебаний в монокристаллах будут рассмотрены в гл. V. Сзадачами по рассеянию энергии колебаний в поликристаллах можно ознакомиться в работах [188; 9, с. 229, 239 и др. ].
Из всего изложенного следует, что термодинамическая теория позволяет обосновать общую линейную связь между напряжением и деформацией и достаточно точно описать характер температурной и частотной зависимостей внутреннего трения, дефектов модулей и т. д. Однако при помощи уравнений (152), (153) и (166), (167) нельзя опи сать амплитудную зависимость этих величин (в силу их линейности!)
6 р. С. Постников |
81 |
и их зависимость от времени, часто наблюдаемую на опыте. Тот факт, что ядра ползучести и релаксации зависят от разности t — t', ука зывает на неизменность свойств материала во времени; начало от счета времени можно изменять произвольно, результаты, при этом остаются одни и те же.
Для описания свойств стареющих материалов ядрэ должны быть приняты функцией аргументов t и t' по отдельности.
Амплитудную зависимость внутреннего трения можно учесть, например, в рамках нелинейной теории вязко-упругости, которая, к сожалению, находится еще в стадии разработки [189]. Можно ее учесть и путем выяснения характера зависимости времени релакса ции от амплитуды деформации, как это сделано в работе Фастова [92, с. 31 ]. Фастов ввел еще один тензор — тензор остаточных де
формаций еік. Кроме этого, он использовал уравнение для устано вившейся ползучести, которое рассматривается в микроскопической теории ползучести. В результате он получил выражение для времени релаксации, которое для одноосного напряжения имеет вид
т = ■ |
(168) |
2GAsh Wуа
Здесь А и у — константы уравнения ползучести.
Наконец, можно решить вопрос амплитудной зависимости вну треннего трения методами нелинейной теории упругости," опираясь на опытные данные.
Сначала полезно рассмотреть случай линейных колебаний (см. раздел 6 ), а затем нелинейных (раздел 7). Будем рассматривать, для простоты задачу о колебании механической системы с одной сте пенью свободы, к которой могут быть сведены многие реальные системы.
Приемы приведения реальных механических систем к указанно^ расчетной схеме хорошо известны [190].
6. Линейные колебания. Резонансные потери
Свободные колебания
Задачи, которые рассматриваются в теории линейных колебаний, как хорошо известно,-, приводятся к изучению дифференциального уравнения
|
|
пи + |
ßx 4- kx = Р (t), |
, |
• (169) |
где т, ß и |
k — постоянные; |
ßx — вязкая или диссипативная сила; |
|||
kx — восстанавливающая или упругая сила; |
Р (t) — внешняя вы |
||||
нуждающая |
сила. |
|
|
|
|
Рассмотрим случай, когда Р (t) = 0. Введем следующие величины: |
|||||
|
Р |
. 2 |
k |
|
(.170) |
|
т |
С00 = |
~т’ |
|
82
Тогда уравнение (169) примет вііД
X4- sx -j- ©о* = 0. |
(І71) |
Общее решение этого линейного уравнения с постоянными коэф фициентами представляется линейной комбинацией двух показатель ных функций:
X = Ae^i Be^ - ( . |
(172) |
Здесь А к В — произвольные постоянные, а q-yи q2— корни характе-' ристического уравнения
<T + |
s<7 + coo = |
0. |
(173) |
Следовательно, qx и q%определяются |
соотношениями |
|
|
Чх, 2 = |
— - | - ± і ю . |
(174) |
Для нас является важным случай, когда со — величина действи
тельная и не равная нулю (а>о > s2).
В этом случае решение уравнения (170) можно представить в виде
_ _s£_ |
(175) |
х = е 2 (сх cos cot -)- с2sin at) ^ |
|
или |
|
_ st |
(176) |
х = а0е 2 cos (со/ -f- cp). |
Здесь суи с2 или а0и ср — произвольные постоянные, их определяют
по заданным значениям смещения х — х0 и скорости х = х0 |
в неко |
||||
торый момент времени t = t0. |
«временем |
релаксации», так |
|||
Величину |
г = 1/s можно назвать |
||||
как квадрат |
амплитуды колебаний |
_ |
t_ |
за время |
t — х |
а2 = а2е |
т |
||||
уменьшается |
в е раз. |
|
|
|
|
Как видим, возбужденная система, рассеивая энергию, которую она приобрела за счет внешнего воздействия, приближается посте пенно к состоянию, равновесия, причем логарифмический декремент
колебаний |
П77\ |
X---- ж |
не зависит от амплитуды. Вынужденные колебания
Пусть на систему действует внешняя сила Р (t), зависящая от времени. Наиболее важным для дальнейшего является случай, когда сила Р (t) периодическая.
Здесь предполагается, что Р (t) есть простая гармоническая
функция |
Р (t) = b cos со^. |
6 * |
83 |
Решение уравнения (169) представляется суммой решения одно родного уравнения и какого-нибудь частного решения неоднородного, уравнения. Если предположить, что свободные колебания имеют вид (175), то решение уравнения (169), где Р (t) = b cos со£, легко полу чить в виде
2 (сх cos at -|- с2 sin cat) -j- |
bcos(«в*—-Ф) |
(178) |
|
(оусо2) 2 -I-sV |
|||
, у |
|
Таким образом, искомое движение получают сложением свободных колебаний и так называемых вынужденных колебаний, которые воз никают от действия внешней силы. По истечении достаточно боль шого промежутка времени свободные колебания практически пре кратятся (s > 0 ) и будут происходить лишь вынужденные колебания с амплитудой
а = -----= - • |
(179) |
т ~\[( С0о — |
+ s2<ü |
Сдвиг фазы ер вынужденных колебаний относительно фазы внеш ней силы Р (t) определяется соотношениями:
cos ф : |
у« |
|
Sin ф : |
SCO |
(180) |
|
т + |
со2) 4- s2« 2 |
|||||
|
y « |
|
||||
|
|
|
Работа, производимая силой Р (t) в течение периода Т, дается
выражением
т
(181)
С другой стороны, если величина ф мала, можно принять энергию колебаний системы W — ab.
Тогда
sin ф. |
|
|
\Ѵ |
|
|
Имея в виду уравнение (179), запишем |
|
|
so |
(182) |
|
■J/^(o)Q—fl>2) 2 -f-s20)2 |
||
|
Нетрудно показать, что резонансное поглощение (182) как функ ция частоты со имеет максимум
/т- 1 |
1 |
/ 4 со0 — 2s2 |
(183) |
|
У тах _ К |
|
4Ü)Q—s2 ’ |
||
|
|
|||
при резонансной частоте |
"I / |
2 |
S2 |
|
|
(184) |
|||
|
У |
щ |
----J- |
и не зависит от амплитуды внешней силы.
84
На рис. 45 представлена зависимость Q_i от со/со0 при различных s/co0. Как видно, изменение s существенно влияет на форму резонанс ного пика: В случае релаксационного поглощения [см. формулы (78), (1 2 1 )] изменение s приводит к сдвигу максимума, но форма пика не претерпевает значительного изменения. Следовательно, чтобы отли чить релаксационное затухание от резонансного, необходимо ис следовать, как влияет на затухание изменение s, вызванное, напри мер, изменением температуры тела. Изменение температуры сильно влияет на релаксационное затухание [см. формулу (93) ] и мало
влияет на резонансное поглощение в силу медленного падения соо =
= klm с увеличением темпера |
<?" |
||||
туры. |
|
|
|
|
|
Заметим, что, когда коэффи- ^0 |
|||||
циент ß |
и приведенная масса т |
|
|||
стремятся к бесконечности, а их °-8 |
|||||
отношение s |
остается |
конечным, |
|
||
различие между резонансным и-ре- |
|
||||
лаксационным |
поглощением энер |
Oi |
|||
гии |
теряется. |
Тогда |
уравнение |
||
(171) можно переписать в виде |
0.2 |
||||
|
|
. . . |
2 |
||
|
|
х = —sx — СОо*. |
|
||
Разделив обе части |
уравнения |
|
|||
на X, |
получим |
|
|
|
|
|
1 |
|
2 |
X |
* Рис. 45. Зависимость резонансных потерь |
|
——X — — S — CÖ0 |
I • |
от параметра s/ci) 0 |
||
|
X |
j |
|
X |
|
Если при фиксированном s величины т и ß стремятся к бесконеч ности, то CÖQ= klm стремится к нулю, так что — х = —s.
Решение этого уравнения имеет типичный релаксационный вид:
X = x0e~si + хг. |
- |
Обратимся теперь к построению петли гистерезиса для случая,
когда со Ф со0. Сумма сил упругого kx и вязкого ßx сопротивлений для установившиеся вынужденных колебаний [см. (178) ] равна
F = ka cos (at— cp) — ßcoa sin (at— qp). |
(185) |
Исключив время t из уравнения (178), которое для установив шихся вынужденных колебаний будет иметь вид
X = а cos (co^ — cp),
и уравнения (185), получим
' |
F = kx ± ßcoaV ' — £ • |
(186) |
85
Петля динамического гистерезиса имеет форму эллипса (рис. 46). Полуоси, измеренные в направлениях координатных осей, состав ляют а и ßcoa. Соответственно этому площадь петли гистерезиса равна
S = nßcoa2. |
(187) |
Заменив в (181) амплитуду b и фазовый угол ср по формулам (179) и (180), можно легко убедиться в том, что ДЦ7 = S. Следова тельно, поглощение энергии колеба ний можно определить по площади
петли гистерезиса.
7.Нелинейные колебания
Вдействительности вязкие силы и восстанавливающие силы часто
|
могут |
быть |
нелинейными. В общем |
|
|
виде нелинейное уравнение для си |
|||
|
стемы с одной степенью свободы |
|||
|
можно записать так: |
|||
|
|
тх |
ср (х) + |
f (х) = Р (t). (188) |
Рнс. 46. Петля гистерезиса упруго |
Здесь тх — сила |
инерции; ср (х) — |
||
вязкой системы |
вязкая сила; / (х) — восстанавливаю |
|||
|
щая |
сила. |
|
|
Петли гистерезиса для таких систем получаются не эллиптиче скими, а с острыми вершинами и косой симметрией [2, 191 ]. Площадь петель пропорциональна не квадратам, а другим степеням ампли туды деформации и часто почти не изменяется с частотой, начиная от со «=: 0 (статические испытания) до частот в несколько десятков килогерц. Последнее обстоятельство и особенно наличие статических петель гистерезиса как будто указывают на отсутствие прямой связи между внутренним трением и скоростями деформации. Опытные дан ные [192 ] позволяют связать внутреннее трение в этих случаях только с напряжением. Возможность таких случаев вытекает из следующих рассуждений.
Пусть при нагружении каждому напряжению соответствует дефор мация, мгновенно (со скоростью, не превышающей скорости звука) принимающая постоянное значение. После достижения определенной максимальной деформации нагрузка снимается. Предположим, что каждому значению убывающего напряжения соответствует определен ное мгновенно установившееся значение деформации, отличное, однако, от соответствующего значения при нагружении. После снятия нагрузки остается постоянная остаточная деформация, которую можно довести до нуля, только приложив силу в противоположном направлении. Для этого случая соотношение напряжение — дефор мация представляет многозначную функцию (рис. 47). Физически это означает, что при нагружении (и разгрузке) происходит необра тимое изменение состояния тела. Механизмы затухания этого типа называют статическим гистерезисом [83]. Они исключительно раз
86
нообразны и не поддаются такому простому математическому опи санию, как релаксационные явления (см. гл. V). Явления статиче ского гистерезиса возможны при больших (приводящих к микро скопической деформации) и при малых [результат особой атомной или магнитной перестройки) (см. гл. V) ] амплитудах.
Величина затухания характеризуется по-прежнему величиной, выраженной уравнением (2). Очевидно, что эта мера дает интеграль ную оценку относительного рассеяния энергии за цикл и выражает площадь (AW = S) петли гистерезиса, но не описывает ее формы. Из-за малой ширины петли гистерезиса [для металлов площадь петли
составляет десятые доли процента от пло |
|
щади треугольника, определяющего энергию |
|
системы W (см. рис. 2 ) ] опыт дает возмож |
|
ность составить лишь весьма приблизитель |
|
ное представление о ее форме. Если бы мы |
|
знали точно форму петли, то не трудно было |
|
бы решить обратную задачу, т. е. по форме |
|
петли найти закон изменения вязких и вос |
|
станавливающих сил. К сожалению, пока |
|
еще нет надежных экспериментальных и тео |
|
ретических оснований (наши знания о ме |
|
ханизмах гистерезисного внутреннего тре |
|
ния весьма ограниченны, а опытные данные |
рис. 4 7 . схема петли стат»- |
весьма противоречивы), которые ПОМОГЛИ бы |
ческ°го гистерезиса |
нам решить эту очень важную задачу. Это обстоятельство открыло широкий простор для различных умозри
тельных предположений относительно закона вязких и восстанав ливающих сил, поскольку форма петли при решении конкретной задачи о колебании системы выбирается в значительной мере про извольно [192—195].
Гистерезисные потери ч
Уравнение колебаний системы с одной степенью свободы имеет вид (180). Перепишем его в другом виде, более удобном для анализа, при этом имея в виду только моногармоническое воздействие:
X -f- юо (х + аФ) = h cos at. |
, |
(189) |
Здесь X— обобщенная координата; а — амплитуда колебаний; со„ = |
||
-----собственная частота упругих колебаний; |
т — приве |
денные значения коэффициента жесткости и массы соответственно; h — отнесенная к единице массы амплитуда внешней моногармони
ческой силы; Ф (х, X, а) — функция, определяющая петлю гисте резиса.
Для решения уравнения (189) можно использовать два приема. Один из них основан на том предположении, что силы, обусловлива ющие гистерезис, вносят малую поправку в законы*Гука и Ньютона
87
и, следовательно, их можно учитывать как возмущение. Решение задачи находят асимптотическими методами Крылова—Боголюбова разложения по малому параметру. При этом за нулевое приближение принимается точное решение вполне упругой задачи, а все эффекты, связанные с энергетическими потерями, проявляются в высших гар мониках, из которых обычно ограничиваются только'первым. Такое ограничение связано не только с прогрессирующим усложнением выкладок, необходимых для получения высших приближений; глав ное состоит в другом; построение высших приближений имеет смысл при достаточно точно заданном характере возмущающих сил, чего как раз и нет в наших задачах. Возмущающие силы, входящие в урав нения движения, описываются в значительной мере произвольно выбираемыми выражениями. Очевидно, что в этих условиях построе ние высших приближений может дать лишь иллюзию повышения точности решения. Обоснование и эффективное применение этого метода принадлежит Писаренко [192].
Другой прием, предложенный Пановко [193], не требует введе ния малого параметра и основан на условии энергетической эквива лентности равновеликих петель независимо от формы на контуре. Тогда, выбрав эллиптическую петлю, можно легко получить точное решение задачи, которое, однако, не учитывает влияния высших гар моник, неизменно сопровождающих стационарные колебания нели нейной системы.
Так как первое приближение асимптотического метода эквива лентно эллиптизации произвольной петли, то влияние высших гар моник можно учесть,' выбрав в качестве нулевого приближения точ ное решение для эллиптической петли, а малое отклонение ее от заданной формы — как возмущение. Построение нулевого прибли жения в этом случае не вызывает принципиальных затруднений для достаточно широких петель, а для узких совпадает с первым прибли жением известного .метода Крылова—Боголюбова основной гармо
ники. |
Этот |
обобщенный |
прием |
подробно рассмотрен в работах |
||
[196, |
197]. |
этим |
работам, |
будем |
искать решение уравнения |
(189) |
Следуя |
||||||
в виде |
|
|
|
|
|
|
|
|
X |
= а cos Ѳ, |
0 = |
оо/Ң-фҢ-чр (t), |
(190) |
где а и составляющая фазы ср от времени не зависят. Так как для любой симметричной петли Ф (Ѳ) = —Ф (0 + зі), то ряд Фурье для Ф содержит .только нечетные гармоники:
Ф = С cos 0 + D sin Ѳ+ ß £ [C] cos (2/ + 1)0 + D -s(2/ + 1) 0]. (191) i=i
Здесь и в дальнейшем ß —• параметр, учитывающий отклонение петли произвольной формы от эллипса, и, кроме того, принято усло
вие со +> т|>.
В результате Громоздких вычислений для внутреннего трений, за меру которого принята величина, обратная добротности колеба тельной системы, получается следующая формула [196]:
Q-1= Оо1 1 + «I (4 |
дР |
1 |
дС \~ |
(192) |
" ■да0 1 |
+ С (а0) |
да0 / . |
||
где Q^ 1 — внутреннее трение |
нулевого |
приближения; |
а0 и аг — |
|
амплитуды в нулевом (ß = 0 ) |
и в первом (ß = |
1) приближениях. |
||
Из формулы (192) следует два важных вывода: |
1) форма гистере |
зисной петли не влияет на величину внутреннего трения, если коэф фициенты разложения ее в ряд Фурье (С и D) не зависят от.амплитуды колебаний; 2 ) высшие гармоники могут привести к частотной зависимости гистерезисного внутреннего трения, которая опреде ляется (рис. 48) амплитудной добавкой ах!ай.
Для иллюстрации изложенных выше расчетов рассмотрим кон
кретные примеры. |
[198]. В наших обозначениях |
||
1. Гипотеза Н. Н. Давиденкова |
|||
Ф |
■і-а"- signx ( l + |
S i g n x '^ - ) " — 2 " - 1 |
(193) |
Здесь [Xи /г — коэффициенты, значение которых определяется фор
мой петли; sign — знак, равный +1, если х >> Ѳ, и —1, если л; < 0. Для функции Ф, заданной выражением (191), коэффициенты раз
ложения в ряд Фурье имеют вид
С = — 2\шп~1 |
(2л— 1)1! |
г |
_ г (я— 1)(п — 2) . |
(194) |
||
|
(я + 1)1 ’ |
1 |
(п + 2)(п+3) ’ |
|||
D = гца" - 1 |
1 — гг |
DX= D |
л3 — 12л2-)- 17л — 6 |
(195) |
||
|
1 + л |
’ 1 |
3(л— 1)(л + 2)(л + 3) |
|
Для внутреннего трения в этом случае (см. [196]) имеем формулу
- Ң = |
1 + ^ - ( 1 - -рло)-1. |
(196) |
Qö" 1 |
V |
|
В качестве численного примера рассмотрим сталь. Для нее п = 2, [X= 18,6. Учитывая, что р.а0 <С 1, рассмотрим лишь a ja Ö.Приведен ная на рис. 48, а частотная зависимость относительной величины амплитудной добавки в очень узкой области резонанса первой гар моники достигает 20%. Для упругих частот эта добавка значительно меньше и принимает положительные и отрицательные значения.
Рис. 48, б иллюстрирует частотную зависимость фазовой.добавки
Фх, которая меняет знак |
в области резонанса первой |
гармоники. |
- 2.. Гипотеза Д. Ю. Панова [199]. В наших обозначениях |
||
® = |xsignx (l — |
(197) |
|
■В этом случае коэффициенты разложения равны: |
|
|
с = сх^ 0, |
ß = -f-“ n '’ D1= — ~ . |
(198) |
89