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

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

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

Учитывая уравнение (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, Фікіт (tt’)

■ядра последействия

или ползучести.

 

Некоторые конкретные примеры решения задач по рассеянию энергии колебаний в монокристаллах будут рассмотрены в гл. 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)

" 1

V

 

В качестве численного примера рассмотрим сталь. Для нее п = 2, [X= 18,6. Учитывая, что р.а0 <С 1, рассмотрим лишь a ja Ö.Приведен­ ная на рис. 48, а частотная зависимость относительной величины амплитудной добавки в очень узкой области резонанса первой гар­ моники достигает 20%. Для упругих частот эта добавка значительно меньше и принимает положительные и отрицательные значения.

Рис. 48, б иллюстрирует частотную зависимость фазовой.добавки

Фх, которая меняет знак

в области резонанса первой

гармоники.

- 2.. Гипотеза Д. Ю. Панова [199]. В наших обозначениях

® = |xsignx (l —

(197)

■В этом случае коэффициенты разложения равны:

 

с = сх^ 0,

ß = -f-“ n '’ D1= — ~ .

(198)

89

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