close

Вход

Забыли?

вход по аккаунту

?

5847.Прикладные задачи динамики твердого тела.

код для вставкиСкачать
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Министерство образования и науки Российской Федерации
Федеральное государственное бюджетное образовательное учреждение
высшего профессионального образования
«Оренбургский государственный университет»
Л.И. Кудина, Ю.Л. Власов
ПРИКЛАДНЫЕ ЗАДАЧИ ДИНАМИКИ
ТВЕРДОГО ТЕЛА
Рекомендовано Ученым советом федерального государственного
бюджетного
образовательного
учреждения
высшего
профессионального образования «Оренбургский государственный
университет» в качестве учебного пособия для студентов,
обучающихся
по
программам
высшего
профессионального
образования по направлениям подготовки 190109.65 Наземные
транспортно-технологические средства, 190600.62 Эксплуатация
транспортно-технологических машин и комплексов
Оренбург
2014
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
УДК 531.39(075.8)
ББК 22.213я73
K88
Рецензент – профессор, доктор технических наук М. И. Филатов
Кудина, Л.И.
Прикладные задачи динамики твердого тела : учебное пособие / Л.И. Кудина,
Ю.Л. Власов; Оренбургский гос. ун-т. – Оренбург: ОГУ, 2014. – 118 с.
ISBN
К88
Учебное пособие содержит 5 разделов, соответствующих основным
разделам элективной дисциплины «Прикладные задачи динамики твердого
тела»: устойчивость положения равновесия механических систем, малые
колебания систем с одной и двумя степенями свободы вблизи положения
устойчивого
равновесия,
элементарная
теория
удара,
динамика
плоскопараллельного движения твердого тела.
Учебное пособие предназначено для самостоятельной работы студентов по
направлениям подготовки 190109.65 «Наземные транспортно-технологические
средства», 190600.62 «Эксплуатация транспортно-технологических машин и
комплексов» всех форм обучения.
УДК 531.39(075.8)
ББК 22.213я73
ISBN
© Кудина Л.И.
Власов Ю.Л., 2014
© ОГУ, 2014
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Содержание
Введение ............................................................................................................................... 5
1 Устойчивость положения равновесия ............................................................................ 6
1.1 Некоторые сведения из курса теоретической механики ........................................... 6
1.2 Понятие об устойчивости положения равновесия ................................................... 10
1.3 Потенциальная энергия консервативной системы с конечным числом
степеней свободы .............................................................................................................. 13
1.4 Теорема Лагранжа – Дирихле. Критерий Сильвестра ............................................. 15
1.5 Устойчивость равновесия неконсервативных механических систем .................... 21
2 Малые колебания системы с одной степенью свободы ............................................. 23
2.1 Общие сведения о колебаниях ................................................................................... 23
2.2 Дифференциальное уравнение малых колебаний системы с одной степенью
свободы в общем случае ................................................................................................... 28
2.2.1 Кинетическая энергия системы с одной степенью свободы
в обобщенных координатах.............................................................................................. 29
2.2.2 Обобщенная сила для системы с одной степенью свободы ................................ 30
2.3 Свободные колебания механической системы с одной степенью свободы ......... 35
2.4 Влияние сил сопротивления на свободные колебания системы с одной степенью
свободы ............................................................................................................................... 39
2.5 Вынужденные колебания системы при гармоническом возбуждении.................. 47
2.5.1 Способы возмущения колебаний. Определение обобщенной силы QqB ( t ) ....... 48
2.5.2 Вынужденные колебания при отсутствии сопротивления .................................. 52
2.5.3 Вынужденные колебания при наличии линейного сопротивления .................... 58
2.6 Вынужденные колебания в случае периодической возмущающей силы ............. 62
2.7 Вынужденные колебания при произвольной возбуждающей силе ....................... 65
2.8 Основы теории регистрирующих приборов ............................................................. 68
2.9 Основы виброзащиты.................................................................................................. 72
3 Малые колебания системы с двумя степенями свободы ........................................... 76
3
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
3.1 Дифференциальное уравнение малых колебаний системы с двумя степенями
свободы ............................................................................................................................... 76
3.2 Свободные колебания системы с двумя степенями свободы ................................. 80
3.2.1 Дифференциальные уравнения свободных колебаний ........................................ 80
3.2.2 Общее решение дифференциальных уравнений свободных колебаний системы
с двумя степенями свободы.............................................................................................. 81
3.2.3 Главные координаты ................................................................................................ 85
3.3 Влияние линейного сопротивлении на собственные колебания ............................ 89
3.4 Вынужденные колебания системы без учета сопротивления ................................ 94
3.5 Влияние линейного сопротивления на вынужденные колебания .......................... 97
4 Теория удара ................................................................................................................... 99
4.1 Основное уравнение элементарной теории удара ................................................... 99
4.2 Общие теоремы теории удара .................................................................................. 101
4.2.1 Теорема об изменении количества движения системы при ударе .................... 101
4.2.2 Теорема об изменении главного момента количества движения
системы при ударе ........................................................................................................... 102
4.3 Коэффициент восстановления при ударе о неподвижную поверхность ............. 103
4.4 Прямой центральный удар двух тел ........................................................................ 103
4.5 Теорема Карно ........................................................................................................... 105
4.6 Удар по твердому телу, вращающемуся вокруг неподвижной оси ..................... 106
5 Динамика плоского движения..................................................................................... 108
5.1 Дифференциальные уравнения плоского движения тела ..................................... 108
5.2 Задача о качении диска по наклонной плоскости .................................................. 111
5.2.1 Качение диска по наклонной поверхности при отсутствии скольжения ......... 111
5.2.2 Качение диска по наклонной поверхности при наличии скольжения .............. 114
5.2.3 Качение диска по наклонной поверхности с учетом трения качения .............. 117
Список использованных источников ............................................................................ 118
4
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Введение
Содержание настоящего учебного пособия соответствует курсу лекций по
дисциплине «Прикладные задачи динамики твердого тела» для студентов,
обучающихся по направлениям подготовки 190109.65 «Наземные транспортнотехнологические
средства»
и
190600.62
«Эксплуатация
транспортно-
технологических машин и комплексов».
Элективная дисциплина «Прикладные задачи динамики твердого тела»
вариативной части математического и естественнонаучного цикла основных
образовательных программ (ООП) указанных направлений подготовки является
логическим продолжением курса теоретической механики.
Изложение материала пособия предполагает знакомство читателя с основами
аналитической механики и, прежде всего, с использованием дифференциальных
уравнений
Лагранжа
второго
рода
для
получения
уравнений
движения
механических систем с несколькими степенями свободы.
В
первом
разделе
приведены
краткие
теоретические
сведения
из
предшествующего курса теоретической механики, необходимые для дальнейшего
изложения материала. Для удобства определения расположены в алфавитном
порядке. Основное содержание первого раздела посвящено изложению вопросов,
связанных с исследованием характера положения равновесия механических систем.
Во втором и третьем разделе рассматриваются малые колебания систем
твердых тел с одной и двумя степенями свободы вблизи положения устойчивого
равновесия.
Четвертый раздел содержит изложение основ элементарной теории удара.
Пятый раздел учебного пособия посвящен решению основных задач динамики
в случае плоскопараллельного движения твердого тела.
Завершает учебное пособие список литературы, рекомендуемой для более
глубокого самостоятельного изучения вопросов, не отраженных в настоящем
учебном пособии.
5
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
1 Устойчивость положения равновесия
1.1 Некоторые сведения из курса теоретической механики
Приведем краткие сведения о некоторых понятиях и определениях,
встречающихся в тексте данного учебного пособия и относящихся к различным
разделам курса теоретической механики.
Возможное перемещение точки механической системы – воображаемое
бесконечно малое перемещение, допускаемое наложенными на эту точку связями.
Геометрические связи – связи, налагающие ограничения на положения
(координаты) точек механической системы.
Голономная система – механическая система материальных точек, все связи
которых являются голономными.
Голономные связи – геометрические связи, уравнения (неравенства) которых
содержат только координаты точек системы и не содержат производные от этих
координат по времени (проекции скоростей точек).
Диссипативные силы – силы, вызывающие убывание (диссипацию) полной
механической энергии системы за счет перехода этой энергии в немеханические
формы, например, в теплоту. К диссипативным силам относятся силы вязкого и
сухого трения, силы аэродинамические сопротивления и т.п.
Идеальные связи – связи, сумма элементарных работ сил реакций которых на
любом возможном перемещении равна нулю:
n
∑ Rk δ rk = 0 ,
(1.1)
k =1
где Rk - вектор реакции k-ой связи;
δ rk - вектор возможного перемещения точки приложения k-ой реакции связи.
Консервативная
система
–
механическая
система
с
идеальными,
голономными, стационарными связями, на точки которой действуют только
потенциальные (консервативные) силы.
6
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Обобщенные координаты – независимые параметры любой размерности,
однозначно определяющие положение всех точек механической системы. Для
систем с голономными связями число обобщенных координат равно числу степеней
свободы системы. В дальнейшем для обобщенных координат принято обозначение
q1 , q2 , ..., qS .
Приведем несколько примеров, поясняющих выбор обобщенных координат
для механических систем с различным числом степеней свободы.
х
φ
а)
б)
Рисунок 1
Например, при прямолинейном движении твердого тела его положение в
любой момент времени можно однозначно определить заданием одной обобщенной
координаты q = x , в качестве которой выбрано отклонение тела от некоторого
начального положения (рисунок 1, а).
φ
R
2
1
х
Рисунок 2
7
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Положение твердого тела, вращающегося относительно неподвижной оси
(рисунок 1, б), можно определить заданием одной обобщенной координаты q = ϕ , в
качестве которой принят угол отклонения тела от вертикального положения
равновесия.
Положение механической системы с одной степенью свободы, изображенной
на рисунке 2, определяется одной обобщенной координатой, в качестве которой
можно выбрать либо отклонение груза 1 от положения равновесия q = x , либо угол
поворота барабана 2 q = ϕ , так как данные параметры связаны между собой
соотношением x= ϕ ⋅ R .
φ
2
1
х
Рисунок 3
В случае если тела 1 и 2 (рисунок 3) соединены пружиной, механическая
система будет иметь две степени свободы, а ее положение будет определяться двумя
обобщенными координатами q1 = x и q2 = ϕ , так как если запретить вращение тела
2, то тело 1 все равно сможет перемещаться за счет растяжения (сжатия) пружины.
На рисунке 4 изображена система с тремя степенями свободы, положение
которой определяется тремя обобщенными координатами q1 = x1 , q2 = x2 и q3 = ϕ .
Обобщенные скорости – производные по времени от обобщенных
координат, в дальнейшем обозначаются q1 , q2 , ..., q S .
Потенциальные силы – силы, работа которых зависит только от начального
и конечного положения точки приложения силы и не зависит от траектории, по
8
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
которой происходит перемещение этой точки, например, силы тяжести, силы
упругости.
Рисунок 4
Потенциальная энергия механической системы в данном положении –
сумма работ потенциальных сил, действующих на точки системы, при
перемещении системы из рассматриваемого (данного) положения в начальное
(нулевое).
Стационарные связи – связи, уравнения (неравенства) которых не содержат
время t в явном виде.
Удерживающие (двухсторонние, неосвобождающие) связи – связи, действие
которых сохраняется во все время движения системы. Удерживающие связи
описываются уравнениями в отличие от неудерживающих, описываемых
неравенствами.
Уравнение (неравенство) связи – математическая запись ограничения,
налагаемого связью на точку (точки) механической системы.
Число степеней свободы – число независимых возможных перемещений
точек механической системы.
9
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
1.2 Понятие об устойчивости положения равновесия
Положением равновесия называется такое состояние механической системы, в
котором она может находиться сколь угодно долго, если в начальный момент
времени система была приведена в это положение с нулевыми скоростями.
Положение равновесия механической системы может быть устойчивым,
неустойчивым или безразличным.
Рассмотрим виды положений равновесия на примере одного твердого тела
(рисунок 5).
N
N
N
Р
P
P
a)
б)
в)
Рисунок 5
Положение равновесия механической системы называется устойчивым, если
при сколь угодно малом отклонении от этого положения система стремится
вернуться в исходное состояние и в последующем движении весьма мало
отклоняется от него (рисунок 5, а).
Положение равновесия механической системы называется неустойчивым,
если при сколь угодно малом отклонении от этого положения система все более
удаляется от него и не проявляет стремления вернуться в исходное состояние
(рисунок 5, б).
10
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Положение равновесия механической системы называется безразличным,
если при отклонении от него система и в новом положении может оставаться в
равновесии (рисунок 5, в).
В общем случае кроме начальных отклонений точкам системы следует
сообщить и некоторые начальные скорости. Тогда безразличное равновесие следует
отнести к неустойчивому, так как точки системы будут удаляться от него по
инерции.
Наибольший практический интерес представляет устойчивое равновесие, так
как в таком положении система может находиться длительное время.
Строгое определение устойчивости положения равновесия было дано в
работах выдающегося русского ученого А.М. Ляпунова.
Рассмотрим механическую голономную систему материальных точек с S
степенями свободы. Для того чтобы определить положение такой механической
системы в любой момент времени, требуется выбрать S обобщенных координат
q1 , q2 , ..., qS . Условимся отсчитывать обобщенные координаты от положения
равновесия системы, тогда в положении равновесия все обобщенные координаты
будут одновременно равняться нулю: q1= q2= ...= qS= 0.
Выведем систему из положения равновесия, сообщив ее точкам небольшие по
модулю начальные возмущения. Начальные возмущения системы в общем случае
состоят из начальных значений обобщенных координат q10 , q20 , ..., qS0 и начальных
значений обобщенных скоростей q10 , q20 , ..., qS0 .
Обозначим значения обобщенных координат и обобщенных скоростей при
дальнейшем движении системы через q1( t ), q2 ( t ), ..., qS ( t ) и q1( t ), q2 ( t ), ..., qS ( t ) .
Равновесие системы называется устойчивым, если для всякого сколь
угодно малого числа ε > 0 существует такое малое число η ( ε ) > 0 , что при
начальных возмущениях, удовлетворяющих условиям:
qi0 < η ( ε ); qi0 < η ( ε ) ( i =
1, 2, ..., S ) ;
(1.2)
11
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
в дальнейшем движении механической системы выполняются неравенства:
qi ( t ) < ε ; qi ( t ) < ε .
(1.3)
Если при устойчивом положении равновесия все обобщенные координаты и
обобщенные скорости с течением времени стремятся к нулю:
lim q=
( i 1,2,...,S ),
i ( t ) 0; lim q=
i ( t ) 0=
t →∞
t →∞
(1.4)
то такое положение равновесия называется асимптотически устойчивым.
Приведенному определению устойчивого положения равновесия можно дать
наглядную геометрическую трактовку (рисунок 6). Обозначим равновесные
положения точек механической системы через А1, А2, …, Аn. Опишем вокруг этих
точек сферы малого радиуса ε. Тогда если положение равновесия системы
устойчиво, то как бы ни было мало ε, всегда можно найти такое значение начальных
возмущений η(ε), что во все время дальнейшего движения каждая точка
механической системы останется внутри этой описанной сферы. Разумеется, чем
меньше ε, тем меньше η(ε).
А2
ε
ε
Аk
ε
Аn
А1
ε
Рисунок 6
12
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
1.3 Потенциальная энергия консервативной системы с конечным числом
степеней свободы
Потенциальная энергия П консервативной системы с S степенями свободы
является функцией обобщенных координат этой системы:
П = П( q1 ,q2 ,...,qS ).
(1.5)
Разложим потенциальную энергию в окрестности положения равновесия в ряд
Маклорена по степеням обобщенных координат q1 ,q2 ,...,qS :
 ∂П 
 ∂П 
+
+
+
П( q1 ,q=
,...,q
)
П(
0,0,...,0
)
q
...

 qS +
2
S

 1
q
∂
S
 ∂q1 0

0
 ∂2 П  2
 ∂2 П 
1  ∂ 2 П  2
q
...
q
2
+ 
+
+
+

 1

 q1q2 + ... +
S
2
2  ∂q12 0
q
q
∂
∂
q
∂
 1 2 0
 S 0

(1.6)

 ∂2 П 
+ 2
 qS −1qS  + ...

 ∂qS −1∂qS 0
Здесь и далее обозначение ( ...)0 показывает, что соответствующее выражение
вычисляется в положении равновесия при q1= q2= ...= qS= 0.
Из курса теоретической механики известно, что для равновесия механической
системы с идеальными, стационарными и голономными связями необходимо и
достаточно, чтобы все обобщенные силы системы равнялись нулю:
=
Qi 0=
( i 1,2,...,S ).
(1.7)
13
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
В случае если все действующие на систему силы потенциальны, то
обобщенные силы равны взятым с противоположным знаком частным производным
от потенциальной энергии системы по соответствующим обобщенным координатам:
∂2 П
Qi =
−
(i =
1,2,...,S ).
∂qi
(1.8)
Следовательно, для консервативной системы в положении равновесия должно
выполняться условие:
 ∂2 П 
( i 1,2,...,S ).
 =
 0=
∂
q
 i 0
(1.9)
Будем считать, что в положении равновесия потенциальная энергия системы
равна нулю:
=
П0 П(
=
0,0,...,0 ) 0.
(1.10)
Введем обозначения:
 ∂2 П
cij = 
 ∂qi ∂q j

Коэффициенты
(1.11)
называются

 .
0
обобщенными
(1.11)
коэффициентами
жесткости (квазиупругими коэффициентами) и вычисляются в положении
равновесия при q1= q2= ...= qS= 0. Из выражения (1.11) следует, что
=
cij const;
=
cij c=
( i, j 1,2,...,S ).
ji
(1.12)
14
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Очевидно, что единицы измерения обобщенных коэффициентов жесткости cij
определяются единицами измерения выбранных обобщенных координат q1 ,q2 ,...,qS .
C учетом (1.10), (1.11) выражение для потенциальной энергии (1.6) примет
вид:
П
=
1
 c11q12 + ... + cSS qS2 + 2c12 q1q2 + ... + 2cS −1,S qS −1qS  + ...

2
(1.13)
Будем рассматривать только малые смещения точек системы из положения
равновесия. В этом случае обобщенные координаты qi ( i = 1,2,...,S ) , отсчитываемые
от положения равновесия, можно считать величинами первого порядка малости.
Отбросим в выражении (1.13) слагаемые выше второго порядка малости. В
результате выражение для потенциальной энергии консервативной системы примет
вид:
1 S S
П( q1 ,q2 ,...,qS ) = ∑ ∑ cij qi q j .
2=i 1 =j 1
(1.14)
Функция вида (1.14) называется квадратичной формой. Следовательно,
потенциальная энергия механической системы вблизи положения равновесия
является однородной квадратичной формой обобщенных координат с точностью до
величин второго порядка малости включительно.
Такое представление потенциальной энергии значительно упрощает решение
многих технических задач, а с другой стороны, обеспечивает точность, достаточную
для практических целей.
1.4 Теорема Лагранжа – Дирихле. Критерий Сильвестра
Критерий, дающий простое достаточное условие устойчивости равновесия
механической системы материальных точек, был впервые указан Ж.Л.Лагранжем.
15
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Приведенное Лагранжем недостаточно строгое доказательство теоремы уточнил
П.Г.Л. Дирихле, в связи с чем она и получила название теоремы Лагранжа –
Дирихле:
Если в положении равновесия консервативной системы потенциальная
энергия имеет строгий минимум, то это положение равновесия устойчиво.
Поясним понятие строгого минимума. Функция нескольких переменных
V ( q1 ,q2 ,...,qn ) имеет минимум при значениях аргументов q*1 ,q*2 ,...,q*n в том случае,
если при
qi − q*i < ε ( i =
1,2,...,n ) , где ε > 0 – малая величина, имеет место
неравенство
V ( q1 ,q2 ,...,qn ) ≥ V ( q*1 ,q*2 ,...,q*n ) .
(1.15)
Если при тех же условиях выполняется строгое неравенство
V ( q1 ,q2 ,...,qn ) > V ( q*1 ,q*2 ,...,q*n ) ,
(1.16)
то говорят о строгом, или изолированном, минимуме.
Например, изображенная на рисунке 7 функция П( q ) имеет строгий минимум
при значениях аргумента равного q = q3 и q = q5 . При значении аргумента q = q6
функция имеет нестрогий минимум.
Перейдем к доказательству теоремы Лагранжа – Дирихле. Предположим, что
потенциальная энергия механической системы материальных точек M 1 ,M 2 ,...,M n в
положении равновесия A1 ,A2 ,...,An (рисунок 6) имеет минимальное значение П min .
Выясним, какие значения может принимать потенциальная энергия системы, когда
одна из точек системы mk будет находиться на поверхности сферы радиуса ε ,
окружающей соответствующую точку Ak . Очевидно, что возможно n таких
значений и каждое из них больше, чем П min . Обозначим меньшее из всех
возможных значений потенциальной энергии
( Пmin + η ) .
Очевидно, что величина
16
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
η > 0 и тем меньше, чем меньше выбранный радиус сферы ε . Таким образом, если
хотя
бы
одна
точка
механической
системы
находится
на
поверхности
соответствующей сферы радиуса ε , то потенциальная энергия системы П должна
удовлетворять условию
П ≥ П min + η .
(1.17)
Дадим точкам механической системы незначительные начальные отклонения
от положений равновесия A1 ,A2 ,...,An и сообщим им малые начальные скорости.
Покажем, что при начальных отклонениях и начальных скоростях, не
превышающих некоторой малой величины η , точки системы во все время
последующего движения не выйдут за пределы соответствующих сфер радиуса ε ,
как бы мала ни была эта величина. Для этого достаточно доказать, что во все время
последующего движения потенциальная энергия системы будет оставаться меньше,
чем ( П min + η ) .
Обозначим значения кинетической и потенциальной энергий системы в
некоторый момент времени Т и П, соответственно, а их значения в начальный
момент времени – Т0 и П0. Для рассматриваемой консервативной системы будет
справедлив закон сохранения энергии:
T + П = T0 + П0 ,
откуда
T =T0 + П0 − П .
Так как кинетическая энергия системы T > 0 , то потенциальная энергия
системы П во все время движения должна удовлетворять неравенству:
П < T0 + П0 .
(1.18)
17
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Примем начальные отклонения и начальные скорости настолько малыми,
чтобы выполнялись неравенства:
П0 < П min +
η
2
;
T0 <
η
2
.
(1.19)
С учетом (1.19) неравенство (1.18) примет вид:
П < П min + η .
(1.20)
Однако, полученное выражение (1.20) противоречит соотношению (1.17).
Таким образом, во все время последующего движения точки системы будут
оставаться внутри соответствующих сфер, окружающих положения равновесия
A1 ,A2 ,...,An . Следовательно, положение равновесия системы устойчиво.
Согласно теореме Лагранжа – Дирихле для доказательства устойчивости
положения равновесия консервативной системы достаточно убедиться, что
потенциальная энергия имеет в рассматриваемом положении минимум.
Пусть потенциальная энергия П(q) некоторой консервативной механической
системы с одной степенью свободы зависит от обобщенной координаты q так, как
показано на рисунке 7. Положение равновесия будет устойчивым при значениях
обобщенной координаты q = q3 и q = q5 , так как только при этих значениях
аргумента потенциальная энергия системы имеет строгий минимум.
Для консервативных систем с одной степенью свободы определение
минимума выполняется достаточно просто. В этом случае производная второго
порядка от потенциальной энергии по обобщенной координате, вычисленная в
положении равновесия, должна быть положительна (предполагается, что она
существует).
 ∂2 П 
 2  > 0.
 ∂q 0
(1.21)
18
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Для консервативной системы, имеющей S степеней свободы, устойчивость
рассматриваемого положения равновесия также определяется из условия минимума
потенциальной энергии. Но в этом случае решение вопроса о существовании
минимума потенциальной энергии существенно усложняется.
Рисунок 7
Ранее условлено, что в положении равновесия потенциальная энергия системы
П0 = 0 . Тогда наличие минимума потенциальной энергии должно одновременно
означать, что вблизи положения равновесия потенциальная энергия системы должна
быть положительной. Следовательно, на основании теоремы Лагранжа-Дирихле
потенциальная энергия системы (1.14) должна представлять собой определенно
положительную квадратичную форму обобщенных координат.
Квадратичная форма называется определенно положительной, если она
положительна при всех значениях аргументов, одновременно не равных нулю.
Вопрос о знаке любой квадратичной формы устанавливается теоремой
Сильвестра:
19
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Для того чтобы квадратичная форма была определенно положительной
необходимо и достаточно, чтобы все главные диагональные миноры матрицы
квадратичной формы были положительны.
Матрица квадратичной формы (1.14) имеет вид:
с11 с12
с21 с22
с31 с32
... ...
сS1 сS 2
с13
с23
с33
...
сS 3
...
...
...
...
...
с1S
с2S
с3S
...
сSS
(1.22)
Заметим, что эта матрица симметричная. Главные диагональные миноры
матрицы (1.22) имеют вид:
с11 с12 с13
с11 с12
; ∆3 с21 с22 с23
=
∆ 1 с=
=
11 ; ∆2
с21 с22
с31 с32 с33
с11 с12
с21 с22
;...;
=
∆S с31 с32
... ...
сS1 сS 2
с13
с23
с33
...
сS 3
...
...
...
...
...
с1S
с2S
с3S . (1.23)
...
сSS
Таким образом, критерий определенной положительности квадратичной
формы имеет вид:
∆1 > 0; ∆2 > 0;...,∆S > 0.
(1.24)
Заметим особо, что дискриминант определенно положительной формы не
равен
нулю.
Кроме
того,
нумерация
обобщенных
координат
выбирается
произвольно, т.е. любую из них можно сделать первой, следовательно, все элементы
матрицы (1.22), стоящие на главной диагонали, должны быть строго больше нуля:
cii > 0 ( i =
1,2,...,S ).
(1.25)
20
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Теорема Лагранжа-Дирихле дает критерий, позволяющий утверждать, что
равновесное положение консервативной системы устойчиво, если ее потенциальная
энергия имеет минимум. Однако эта теорема не указывает, каким будет равновесие
системы, если потенциальная энергия в положении равновесия не имеет минимума.
Ответ на этот важный для весьма обширного класса задач вопрос содержится
в теоремах А.М. Ляпунова [7]:
Теорема 1. Равновесие системы неустойчиво, если отсутствие минимума
потенциальной энергии узнается по членам второго порядка малости в
разложении потенциальной энергии без необходимости рассмотрения членов
высшего порядка.
Теорема 2. Равновесие системы неустойчиво, если потенциальная энергия
имеет максимум и наличие этого максимума может быть установлено из
рассмотрения членов наименее высокого порядка малости в разложении
потенциальной энергии.
Последнюю теорему применяют в том случае, если невозможно определить
наличие или отсутствие минимума по членам второго порядка малости, например в
случае, когда эти члены в разложении потенциальной энергии просто отсутствуют
(равны нулю).
1.5 Устойчивость равновесия неконсервативных механических систем
В реальной механической системе всегда существуют силы сопротивления
движению, возникающие благодаря трению или вязкости среды. Такие силы
принято называть диссипативными.
При наличии в системе диссипативных сил для оценки устойчивости
положения равновесия можно дополнительно воспользоваться тремя теоремами
Кельвина [3, 4].
Теорема
1.
Если
положение
равновесия
консервативной
системы
устойчиво при одних только потенциальных силах, то оно будет оставаться
21
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
устойчивым и при добавлении диссипативных сил.
Теорема
2.
Устойчивое
положение
равновесия
становится
асимптотически устойчивым при добавлении диссипативных сил с полной
диссипацией.
Понятие полной диссипации имеет существенное значение для систем с двумя
и более степенями свободы и определено в подразделе 3.3.
Теорема 3. Изолированное и неустойчивое при одних потенциальных силах
положение равновесия
не может быть стабилизировано диссипативными
силами.
Первые две теоремы Кельвина показывают, что наличие диссипативных сил
не может нарушить устойчивость положения равновесия. Третья теорема,
утверждает, что диссипативные силы не могут привести к переходу системы из
неустойчивого положения равновесия в устойчивое. Кроме того, из приведенных
теорем следует, что для выяснения характера положения равновесия реальной
механической
системы
с
диссипативными
силами
ее
можно
заменить
консервативной моделью.
В заключение заметим, что реальные механические системы всегда имеют
бесконечно большое число степеней свободы, т.к. такие системы включают в себя
упругие (деформируемые) тела, отдельные точки которых могут перемещаться
относительно друг друга. Следовательно, для определения положения реальных
механических систем требуется задание бесконечно большого числа обобщенных
координат, определяющих положение каждой отдельной материальной точки
системы. Реальную механическую систему можно рассматривать как систему с
конечным числом степеней свободы только в том случае, если все инерционные
(имеющие массу) элементы системы считать абсолютно твердыми, а все упругие
элементы (например, пружины) безынерционными, т.е. считать, что их массы
пренебрежимо малы по сравнению с массами твердых тел.
22
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2 Малые колебания системы с одной степенью свободы
2.1 Общие сведения о колебаниях
Механическими колебаниями называется движение, в процессе которого
происходит многократное чередование возрастания и убывания некоторых величин,
например обобщенных координат точек механической системы.
Механические колебания можно классифицировать по целому ряду
признаков: по характеру закона изменения обобщенных координат, по
характеру сил, действующих на точки системы и т.п.[3, 4].
По кинематическому признаку, т.е. по характеру закона изменения во
времени обобщенной координаты, все встречающиеся на практике механические
колебания можно разделить на периодические, почти периодические и непериодические
(рисунок 8).
Механические колебания
Периодические
Почти периодические
Непериодические
Рисунок 8 – Классификация колебаний по кинематическому признаку
Периодическим называется такое движение, при котором обобщенная
координата q( t ) , взятая в любой момент времени t, через определенный отрезок
времени Т имеет то же самое значение:
q( t + T ) =
q( t ).
(2.1)
Непериодическим называются движение, не удовлетворяющее сформулированному
условию периодичности (2.1).
23
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Промежуточное положение между периодическим и непериодическим
занимает почти периодическое движение [3, 4].
Обобщенная координата q( t ) при почти периодическом движении в любой момент
времени t удовлетворяет условию:
q( t + T ) − q( t ) ≤ ε .
(2.2)
где T и ε – некоторые постоянные величины.
Если значение ε мало по сравнению со средним значением обобщенной
координаты q( t ) за время T , то можно считать, что почти периодическое движение на
некотором промежутке времени практически близко к периодическому.
Среди периодических огромную роль играют прежде всего гармонические
колебания, при которых обобщенные координаты во времени изменяются по закону
синуса или косинуса.
Непериодические колебания гораздо разнообразнее периодических. Наибольший
практический интерес представляют колебания, происходящие по закону «затухающей
(нарастающей) синусоиды» (рисунки 9, 10), и лимитационные движения (рисунок 11).
q(t)
t
Рисунок 9 – Затухающие гармонические колебания
24
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Колебания
по
закону
«затухающей
синусоиды»
или
затухающие
гармонические колебания (рисунок 9) описываются уравнением вида:
=
q( t ) Ae − nt sin( kt + α ),
(2.3)
где A, n, k , α – постоянные величины.
Очевидно, что закон движения при нарастающих гармонических колебаниях
(рисунок 10) отличается от уравнения (2.3) только знаком величины n .
q(t)
t
Рисунок 10 – Нарастающие гармонические колебания
Обобщенная координата при лимитационном движении изменяет знак
ограниченное число раз и с течением времени стремится к некоторому постоянному
значению.
Математически
лимитационные
движения
могут
быть
описаны
уравнением вида:
=
q( t )
( Ae
− kt
)
+ Bekt e − nt ,
(2.4)
где A, B, k , n – действительные числа.
Примеры лимитационных движений показаны на рисунке 11.
25
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
q(t)
t
Рисунок 11 – Лимитационные движения
Перечисленные признаки механических колебаний являются чисто внешними
и не учитывают всех особенностей происходящих при этом процессов, поэтому
принято классифицировать колебания по целому ряду признаков [3, 4].
По числу учитываемых в расчетной схеме степеней свободы различают:
а) колебания в системах с одной степенью свободы;
б) колебания в системах с конечным числом степеней свободы;
в) колебания в системах с распределенными параметрами (с бесконечным
числом степеней свободы).
По совокупности действующих на механическую систему сил:
а) свободные колебания;
б) вынужденные колебания;
в) параметрические колебания;
г) автоколебания.
Свободные колебания происходят в изолированных механических системах
после некоторого внешнего возмущения под действием только восстанавливающих
сил, в качестве которых чаще всего выступают силы упругости. Характеристики
колебательного процесса в этом случае определяются внутренними параметрами
самой системы, зависящими от ее физического строения.
26
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Вынужденные
вынуждающими
колебания
вызываются
(возмущающими)
и
силами
или
поддерживаются
внешним
внешними
кинематическим
возбуждением. Характер колебаний при этом определяется не только физическими
свойствами самой механической системы, но и существенно зависит от характера
внешнего воздействия.
Параметрические колебания – механические колебания, вызываемые и
поддерживаемые некоторым параметрическим возбуждением, т.е. изменением во
времени
одного
или
нескольких
параметров
системы,
например,
массы,
коэффициента жесткости, момента инерции и т.п.
Автоколебания – незатухающие периодические колебания механической
системы, поддерживаемые за счет внешнего постоянного (непериодического)
воздействия. Отличительной особенностью автоколебательных систем является
способность регулировать поступление энергии от источника движением самой
системы.
По
типу
дифференциальных
уравнений
движения,
описывающих
механические колебания, различают:
а) линейные колебания, при которых движение механической системы
описывается линейными дифференциальными уравнениями;
б)
нелинейные
колебания,
которые
описываются
нелинейными
дифференциальными уравнениями.
Изучение и анализ нелинейных колебаний представляет собой математически
сложную задачу в связи с отсутствием общих методов интегрирования нелинейных
дифференциальных уравнений. Вместе с тем для решения многих практических
задач оказывается возможной линеаризация уравнения движения, т.е. замена точных
нелинейных
уравнений
приближенными
линейными.
Естественно,
что
линеаризованные уравнения не могут в точности отобразить движение системы и
дают искаженную картину явлений, присущих только нелинейным колебаниям.
Вместе с тем, эти искажения тем менее существенны, чем меньше отброшенные
члены уравнений по сравнению с оставшимися.
27
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2.2 Дифференциальное уравнение малых колебаний системы с одной
степенью свободы в общем случае
Рассмотрим механическую систему n материальных точек c одной степенью
свободы, на которую наложены голономные, стационарные, удерживающие связи.
Пусть система находится в положении устойчивого равновесия. Выведем систему из
положения равновесия, сообщив ее точкам малые начальные скорости и отклонения.
Положение системы в любой момент времени t характеризуется одной обобщенной
координатой q( t ) , которую условимся отсчитывать от положения устойчивого
равновесия.
Для составления дифференциального уравнения дальнейшего движения
системы воспользуемся уравнением Лагранжа ІІ рода:
d  ∂T  ∂T
−
=
Qq ,
d  ∂q  ∂q
(2.5)
где T – кинетическая энергия системы;
q – обобщенная координата;
q – обобщенная скорость;
Qq – обобщенная сила, соответствующая выбранной обобщенной координате q( t ) .
В общем случае уравнение (2.5) будет нелинейным, но как показано в разделе 1
при принятых допущениях во все время дальнейшего движения обобщенная
координата q( t ) и обобщенная скорость q( t ) будут оставаться малыми величинами,
и в уравнении (2.5) можно отбросить нелинейные слагаемые. Кроме того
дифференцирование понижает порядок малости на единицу, поэтому в выражении
для кинетической энергии можно ограничиться членами не выше второго порядка
малости.
28
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2.2.1 Кинетическая энергия системы с одной степенью свободы в обобщенных
координатах
Из курса теоретической механики известно, что кинетическая энергия
системы определяется как сумма кинетических энергий всех точек системы:
1 n
T = ∑ mk vk2 ,
2 k =1
(2.6)
где vk – скорость точки массой mk .
Все связи системы являются стационарными, следовательно, радиус-векторы
точек системы зависят только от обобщенной координаты r = r ( q ) , тогда
vk =
drk ∂rk dq ∂rk

=
⋅
=
q.
∂q dt ∂q
dt
(2.7)
С учетом тождества ( vk ) = vk2 и соотношения (2.7) представим выражение
2
для кинетической энергии системы (2.6) в виде:
2
 ∂rk   2 1
1 n
 ∑ mk 
 q
=
=
T
A( q ) ⋅ q 2 ,

2  k =1  ∂q  
2


(2.8)
где A( q ) является некоторой функцией обобщенной координаты q( t ) .
Воспользуемся разложением функции A( q ) в ряд Маклорена по степени
обобщенной координаты q( t ) в окрестности положения равновесия:
 ∂ 2 A  q2
 ∂A 
A( q )= A( 0 ) +   ⋅ q +  2  ⋅ + ...;
 ∂q 0
 ∂q 0 2
(2.9)
29
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Здесь как и ранее обозначение
( ...)0
показывает, что соответствующее
выражение вычисляется в положении равновесия при q = 0.
Т.к. в выражении для кинетической энергии (2.8) следует ограничиться
слагаемыми не выше второго порядка малости, то с требуемой степенью точности
можно считать, что
A( q )= A( 0 )= a= const.
Величина
характеризует
называется
a
инертность
обобщенным
механической
(2.10)
коэффициентом
системы.
Единицы
инерции
и
измерения
обобщенного коэффициента инерции a определяются единицами измерения
обобщенной координаты
q( t ) . Поскольку кинетическая энергия – величина
положительная, то обобщенный
коэффициент инерции
может быть только
положительным ( a > 0 ).
Окончательно
для
системы
с
одной
степенью
свободы
выражение
кинетической энергии в обобщенных координатах примет вид:
1 2
 .
T = aq
2
(2.11)
2.2.2 Обобщенная сила для системы с одной степенью свободы
В общем случае сила, действующая на k-ую точку системы, может быть
функцией положения точки rk , ее скорости vk и времени t:
Fk = Fk ( rk ,vk ,t ).
С учетом малости колебаний представим Fk в виде
Fk = FkП ( rk ) + FkR ( vk ) + FkВ ( t ),
(2.12)
30
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где FkП ( rk ) – потенциальные силы, зависящие только от положения точки;
FkR ( vk ) – силы сопротивления (диссипативные силы), линейно зависящие от
скорости и вызывающие уменьшение полной механической энергии системы;
FkВ ( t ) – вынуждающие (возмущающие) силы, являющиеся функцией времени t .
По аналогии с (2.12) обобщенную силу Qq , соответствующую выбранной
обобщенной координате q( t ) , представим в виде:
Qq = QqП + QqR + QqВ .
(2.13)
где QqП – составляющая обобщенной силы, обусловленная действием потенциальных
сил FkП ( rk ) ;
QqR – диссипативная составляющая обобщенной силы;
QqВ – вынуждающая составляющая обобщенной силы соответственно.
Из курса теоретической механики известно [1, 5], что составляющая обобщенной
силы от действия потенциальных сил равна взятой с противоположным знаком частной
производной от потенциальной энергии по обобщенной координате:
QqП = −
∂П
,
∂q
(2.14)
где П( q ) – потенциальная энергия системы.
Для системы с одной степенью свободы потенциальная энергия в
соответствии с выражением (1.14) примет вид:
П( q ) =
1 2
cq ,
2
(2.15)
31
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где c = c11 – обобщенный коэффициент жесткости, удовлетворяющий условию (1.25).
Получим выражение для составляющей обобщенной силы от действия
диссипативных сил QqR . Предположим, что силы сопротивления, действующие на
точки механической системы, линейно зависят от скорости этих точек:
FkR ( vk ) = − µk vk ,
(2.16)
где µk – коэффициент постоянный коэффициент пропорциональности.
Из (2.7) имеем:
d vk ∂rk
=
.
∂q
dq
Тогда
n
n
dv
∂rk
∂rk
=
− ∑ µ k vk
=
− ∑ µ k vk k =
∑ Fk ∂q =
dq
∂q
k 1=
k 1 =
k 1
=
QqR
n
R
d  n µk vk vk 
d  n µk vk2 
= − ∑
.
 = − ∑


dq
2
dq
2
=
 k 1=

k 1

.
(2.17)
Обозначим
1 n
Ф = ∑ µk vk2 ,
2 k =1
(2.18)
где Ф – диссипативная функция Рэлея или функция рассеивания [3, 7], по форме
аналогичная выражению для кинетической энергии (2.6).
Составляющая обобщенной силы от действия сил сопротивления можно
равна взятой с противоположным знаком частной производной от диссипативной
функции по обобщенной скорости:
32
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
QqR = −
∂Ф
.
∂q
(2.19)
Повторяя рассуждения, приведенные в пункте 2.2.1 для кинетической
энергии, представим диссипативную функцию с требуемой степенью точности по
аналогии с выражением (2.11) в виде:
2
2
 ∂rk   2 1
1 n
1 n
1 2
 drk 
2


 ∑ µk 

µk  =
=
=
⋅
=
Ф
q
B(
q
)
q
bq ,
∑


2
dt
2
q
2
2
∂


 k 1 
 
k 1=
=
(2.20)
где b – постоянный неотрицательный коэффициент, называемый обобщенным
коэффициентом диссипации [4].
Единицы измерения обобщенного коэффициента диссипации, также как и
обобщенных коэффициентов инерции и жесткости, определяются единицами
измерения выбранной обобщенной координаты.
Диссипативная функция Рэлея имеет простой механический смысл: скорость
уменьшения (рассеивания) полной механической энергии системы равна удвоенному
значению диссипативной функции Рэлея [7]:
d( T+Π ) dE
=
= −2Φ ,
dt
dt
(2.21)
где E — полная механическая энергия системы.
Составляющую обобщенной силы от действия вынуждающих сил FkВ ( t ),
зависящих от времени, можно получить известным из курса теоретической
механики стандартным способом, как отношение суммарной возможной работы
вынуждающих сил FkВ ( t ) к вариации обобщенной координаты δ q :
33
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
n
QqВ
=
∑ δ Ak
k =1
δq
.
(2.22)
С учетом полученных соотношений (2.14) и (2.19) выражение для обобщенной
силы (2.13) примет вид:
∂П ∂Ф
Qq =
−
−
+ QqВ .
∂q ∂q
(2.23)
Подставив (2.23) в уравнение движения (2.5), получим
∂П ∂Ф
d  ∂T  ∂T
−
=
−
−
+ QqВ .


d  ∂q  ∂q
∂q ∂q
(2.24)
Как следует из полученного выражения (2.11), кинетическая энергия системы
с точностью до величин второго порядка малости является функцией только
обобщенной скорости q( t ) . Тогда частная производная от кинетической энергии
системы по обобщенной координате тождественно равна нулю
∂T
= 0,
∂q
и уравнение движения (2.24) примет вид:
d  ∂T  ∂Ф ∂П
+
+
=
QqВ .


d  ∂q  ∂q ∂q
(2.25)
С учетом выражений (2.11), (2.15) и (2.20) имеем
34
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
∂T
∂Ф
∂П
d  ∂T 
 =


= aq,
= bq,
=
aq,
cq,


∂q
d  ∂q 
∂q
∂q
(2.26)
Подставляя (2.26) в (2.25), получим
 + bq
 + cq =
aq
QqВ ,
или
1
a
q + 2nq
 + k 2q =
QqВ ,
(2.27)
где введены обозначения
=
2n
b 2 c
=
,k
.
a
a
Выражение (2.27) представляет собой дифференциальное уравнение малых
колебаний системы с одной степенью свободы в общем случае.
2.3 Свободные колебания механической системы с одной степенью
свободы
В случае отсутствия диссипативных и вынуждающих сил уравнение (2.27)
примет вид:
q + k 2 q =
0.
(2.28)
Решение однородного линейного дифференциального уравнения второго
порядка (2.28) с постоянными коэффициентами будем искать в виде
35
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
q( t ) = eλt .
(2.29)
Подставляя (2.29) в (2.28), найдем, что параметр λ должен удовлетворять
характеристическому уравнению
λ2 + k2 =
0.
(2.30)
Характеристическое уравнение (2.30) имеет два мнимых корня:
λ1,2 = ± ki.
В этом случае общее решение дифференциального уравнения (2.29) имеет вид:
=
q C1 cos ( kt ) + C2 sin ( kt ) .
(2.31)
Откуда
q =
− kC1 sin ( kt ) + kC2 cos ( kt ) .
(2.32)
Произвольные постоянные C1 и C2 определяются из начальных условий
движения:
=
q( 0 ) q0 =
, q( 0 ) q0 ,
при t0 0=
где q0 , q0 – начальные значения обобщенной координаты и обобщенной скорости
соответственно.
Подставляя начальные условия в (2.31) и (2.32) получим:
=
C1 q=
0 , C2
q0
.
k
36
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
С учетом найденных значений постоянных интегрирования общее решение
(2.31) примет вид:
=
q q0 cos ( kt ) +
q0
sin ( kt ) .
k
(2.33)
Представим общее решение (2.31) в амплитудной форме:
=
q A sin( kt + α ),
(2.34)
где A и α – произвольные постоянные, определяемые из начальных условий
движения.
Для определения A и α воспользуемся общеизвестной формулой синуса
суммы в выражении (2.34):
q= A  sin (α ) ⋅ cos ( kt ) + cos (α ) ⋅ sin ( kt ) ,
или
q = Asin (α ) ⋅ cos ( kt ) + Acos (α ) ⋅ sin ( kt ) .
(2.35)
Сопоставляя выражения (2.35) и (2.33), получим
 q0 = A sin (α ) ;

 q0
 = Acos (α ) .
k
(2.36)
Откуда
q k
C 
 q 
+  0  ; α = arctg  1  = arctg  0  .
k 
 C2 
 q0 
2
A=
C12
+ C22
=
q02
(2.37)
37
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Из (2.34) следует, что свободные колебания в линейной консервативной
системе с одной степенью свободы являются гармоническими.
На рисунке 12 показан график зависимости обобщенной координаты q от
времени t при свободных колебаниях системы с одной степенью свободы.
Перечислим основные характеристики такого движения:
– круговая, или циклическая частота k , определяющая число полных
колебаний за 2π секунд;
– фаза колебаний ( kt + α ) ;
– начальная фаза колебаний α ;
– амплитуда
колебаний A , определяющая максимальное отклонение
обобщенной координаты
q( t ) от положения равновесия, соответствующего
значению q = 0 ;
– период колебаний T , представляющий собой время совершения одного
полного колебания.
q(t)
А
q0
t
O
Т
Рисунок 12
38
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Период колебаний T и круговая частота колебаний k связаны между собой
соотношением:
=
T
2π
a
.
= 2π
k
c
(2.38)
В инженерной практике часто используют величину, обратную периоду
колебаний и называемую частотой колебаний (число колебаний в единицу времени):
N=
1
k
.
=
T 2π
(2.39)
Как видно из (2.38) круговая частота k и период колебаний T не зависят от
начальных условий (свойство изохронности колебаний). Свойство изохронности
является
следствием
допущения
о
малости
колебаний
и
линейности
дифференциального уравнения (2.28).
2.4 Влияние сил сопротивления на свободные колебания системы с одной
степенью свободы
Свободные
происходят
при
колебания
наличии
механической
сил
системы
сопротивления,
в
реальных
вызывающих
условиях
рассеивание
(диссипацию) механической энергии системы.
В случае наличия сил сопротивления и отсутствия вынуждающих сил
дифференциальное уравнение движения системы (2.27) примет вид:
q + 2nq
 + k 2q =
0,
где n =
(2.40)
b
– коэффициент затухания.
2a
Будем искать решение дифференциального уравнения (2.40) в виде (2.29).
39
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Подставляя (2.29) в (2.40), найдем, что параметр λ должен удовлетворять
характеристическому уравнению
λ 2 + 2nλ + k 2 =
0.
(2.41)
λ1,2 =−n ± n 2 − k 2 .
(2.42)
Откуда
В зависимости от соотношения между величинами n и k возможны три
различных случая:
– случай малого сопротивления n < k ;
– случай критического сопротивления n = k ;
– случай большого сопротивления n > k .
Рассмотрим каждый случай в отдельности.
1. Случай малого сопротивления n < k :
В случае малого сопротивления корни характеристического уравнения (2.42)
являются мнимыми
λ1,2 =−n ± ik ∗ ,
где введено обозначение=
k∗
k 2 − n2 .
Общее решение дифференциального уравнения (2.40) имеет вид:
=
q e − nt ( C1 cos k ∗t + C2 sin k ∗t ) ,
(2.43)
а выражение для обобщенной скорости представляется в виде:
q =
−ne − nt ( C1 cos k ∗t + C2 sin k ∗t ) + e − nt ( −C1k1 sin k ∗t + C2 k1 cos k ∗t ) .
(2.44)
40
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
t0 0,=
q( 0 ) q0=
, q( 0 ) q0 в (2.43) и (2.44),
Подставив начальные условия
получим значения постоянных интегрирования:
q0 + nq0
=
C1 q=
C2
0;
k∗
.
(2.45)
С учетом найденных значений постоянных (2.45) общее решение (2.43)
примет вид:
=
q e − nt ( q0 cos k ∗t +
q0 + nq0
k
∗
sin k ∗t ) .
(2.46)
Общее решение (2.43) можно представить в амплитудной форме:
=
q Ae − nt sin( k ∗t + α ) ,
(2.47)
где амплитуда A и начальная фаза α могут быть определены через постоянные
интегрирования (2.45):
q0 k ∗
C1
 q0 + nq0 
.
+
= arctg
 ; α = arctg
q0 + nq0
C2
 k∗ 
2
A=
C12
+ C22
=
q02
График движения, соответствующего общему решению (2.47), приведен на
рисунке 13. Он представляет собой синусоидальную кривую, ограниченную двумя
кривыми
f1 ( t ) = Ae − nt
и
f 2 ( t ) = − Ae − nt . Движение такого вида называется
затухающими колебаниями.
Не являясь периодическим движением, затухающие колебания вместе с тем
сохраняют некоторые свойства периодичности. Общее решение (2.47) представляет
41
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
собой произведение двух функций – экспоненты и синусоиды с периодом T ∗ =
2π
,
k∗
что приводит к чередованию через равный промежуток времени T ∗ нулей и
максимумов q( t ) (рисунок 13). Поэтому можно считать затухающие колебания
условно-периодическими.
∗
Величина T=
2π
=
k∗
k∗
колебаний, величина =
2π
k −n
2
2
называется условным периодом затухающих
k 2 − n 2 – условной частотой затухающих колебаний, а
функция f1 ( t ) = Ae − nt – условной амплитудой затухающих колебаний.
Очевидно, что условный период затухающих колебаний больше периода
свободных колебаний, определяемого формулой (2.38), т.е. T ∗ > T .
При малых значениях коэффициента затухания n  k условная частота
затухающих колебаний и условный период затухающих колебаний приближенно
равны соответствующим характеристикам свободных колебаний при отсутствии
сопротивления.
q(t)
f1(t)= Аe-nt
q0
t
O
f2(t)= -Аe-nt
Т1
Рисунок 13
42
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Решение (2.47) показывает, что в рассматриваемой системе с линейно-вязким
сопротивлением колебания затухают только при t → ∞ , что не соответствует
действительному движению реальных механических систем, которое всегда
заканчивается за конечный промежуток времени. Это противоречие является
результатом того, что в принятой расчетной схеме не были учтены другие виды
действующих сопротивлений, например силы сухого трения, которые и приводят к
прекращению колебаний через конечный промежуток времени.
1
Введем в рассмотрение величину τ 0 = , называемую постоянной времени
n
затухающих колебаний [3, 4], и найдем отношение двух условных амплитуд
затухающих колебаний, отличающихся друг от друга на постоянную времени τ 0 =
1
n
.
A( t1 )
Ae − nt1
nτ 0
=
= e=
e
− n( t1 +τ 0 )
A( t1 + τ 0 ) Ae
Из полученного отношения видно, что за каждый промежуток времени,
равный τ 0 , условная амплитуда затухающих колебаний уменьшается в e ≈ 2,72 раз.
Через промежуток времени 3τ 0 условная амплитуда уменьшится в e3 ≈ 20 раз, а
затухающие колебания можно условно считать прекратившимися.
Отношение двух последовательных амплитудных значений обобщенной
координаты, взятых через условный период
T ∗ , называется декрементом
колебаний:
A( ti )
Ae − nti
nT ∗
=
=
η =
e
.
∗
A( ti + T ∗ ) Ae − n( ti +T )
Логарифмическим
декрементом
колебаний
называют
(2.48)
натуральный
логарифм от декремента колебаний η :
43
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
δ ln=
(η ) nT ∗ .
(2.49)
Логарифмический декремент колебаний используется для характеристики
медленно затухающих колебаний ( n  k ), когда изменение амплитудных значений
обобщенной координаты за условный период мало:
∆ A=i Ai − Ai +1 << Ai .
(2.50)
Тогда
 Ai 



Ai
∆ Ai
ln 
− ln  1 −
=
=
Ai
 Ai +1 
 Ai − ∆ Ai 

ln 
δ=
 ∆ Ai
.
≈
A
i

(2.51)
Как видно из (2.51), логарифмический декремент колебаний при медленно
затухающих колебаниях характеризует относительное изменение максимальных
отклонений обобщенной координаты за условный период.
Выясним энергетический смысл логарифмического декремента. Найдем
изменение полной механической энергии системы за условный период колебаний. В
положениях максимальных отклонений обобщенная скорость q( ti ) = 0 , а полная
механическая энергия системы равна ее потенциальной энергии, вычисляемой по
формуле (2.15):
cAi2
Е=
П=
;
i
i
2
c(Ai − ∆ Ai )2
.
Еi =
Пi =
+1
+1
2
С учетом (2.50) для медленно затухающих колебаний имеем:
c(Ai2 − 2 Ai ⋅ ∆ Ai )
.
=
Еi +1 Π=
i +1
2
44
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Тогда относительное изменение полной механической энергии системы за
условный период колебаний составит
Еi − Еi +1 2 Ai ⋅ ∆ Ai
∆ Ai
=
=
2
= 2=
δ ψ,
Еi
Ai
Ai2
где ψ – коэффициент поглощения энергии за один период (цикл) колебаний.
Следовательно, удвоенное значение логарифмического декремента колебаний
равно коэффициенту поглощения энергии за один условный период колебаний.
2. Случай критического сопротивления n = k :
В случае критического сопротивления характеристическое уравнение (2.41)
имеет два кратных действительных корня:
λ1,2 = −n .
Общее решение дифференциального уравнения (2.40) имеет вид:
=
q e − nt ( C1 + C2t ) ,
(2.52)
а выражение для обобщенной скорости представляется в виде:
q = e − nt ( − nC1 + C2t (1 − n ) ) .
(2.53)
=
t0 0,=
q( 0 ) q0=
, q( 0 ) q0 в (2.52) и (2.53),
Подставив начальные условия
получим значения постоянных интегрирования
C=
q0 ;
1
q0 + nq0 ,
C=
2
а закон движения в виде
45
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
q e − nt [ q0 + ( q0 + nq0 )t ] .
(2.54)
Общий характер движения в этом случае определяется тем, что при
бесконечном возрастании времени t обобщенная координата q( t ) стремится к нулю.
На рисунке 14 представлено решение (2.54) при различных начальных
условиях.
Очевидно, что движение не имеет колебательного характера. Такое движение
называют в общем случае апериодическим, а в рассматриваемом случае
критического сопротивления – предельно апериодическим.
q(t)
q0
t
0
Рисунок 14
3. Cлучай большого сопротивления n > k :
В случае большого сопротивления характеристическое уравнение (2.41) имеет
два отрицательных действительных корня:
λ1,2 =−n ± n 2 − k 2 .
46
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Общее решение дифференциального уравнения (2.40) в этом случае имеет
вид:
q e − nt  C1et
=

n2 −k 2
+ C2 e − t
n2 −k 2
.


Взяв производную по времени t , определим выражение для обобщенной
скорости:
q e − nt C1et
=

n2 −k 2
(
)
− n + n 2 − k 2 − C2 e − t
n2 −k 2
(n +
)
n 2 − k 2  .

=
t0 0,=
q( 0 ) q0=
, q( 0 ) q0
Подставив начальные условия
в полученные
уравнения, найдем значения постоянных интегрирования C1 и C2 :
q0 + nq0 
q0 + nq0 
1
1
C1 =
 q0 +
 ; C2 =
 q0 −
.
2
k2 
2
k2 
Таким образом, в случае большого сопротивления движение также имеет
апериодический характер, аналогичный представленному на рисунке 14, причем с
увеличением коэффициента n графики растягиваются вдоль оси абсцисс, т.к. с
повышением вязкого сопротивления при прочих равных условиях скорость
движения убывает.
2.5 Вынужденные колебания системы при гармоническом возбуждении
В
случае
действия
возмущающих
сил
дифференциальное
уравнение
вынужденных колебаний линейной системы с одной степенью свободы имеет вид
(2.27).
47
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Предположим, что имеет место гармоническое возбуждение колебаний, т.е.
обобщенная сила QqВ изменяется по гармоническому закону:
=
QqB ( t ) Q0 sin( pt + β ) ,
(2.55)
где Q0 – амплитуда обобщенной силы;
p, β – частота и начальная фаза обобщенной силы соответственно.
2.5.1 Способы возмущения колебаний. Определение обобщенной силы QqB ( t )
Выделяют три характерных способа возбуждения колебаний [3, 4]: силовой,
кинематический и инерционный.
Рассмотрим особенности определения обобщенной силы в каждом из трех
случаев.
1. Силовое возмущение. Пусть математический маятник, связанный с
пружиной (рисунок 15), находится под воздействием внешней силы F( t ) ,
изменяющейся
по
F( t ) F0 sin( pt + β ).
гармоническому закону=
Положение
системы будем характеризовать обобщенной координатой q , в качестве которой
выберем угол отклонения маятника ϕ от вертикального равновесного положения.
Дадим системе возможное перемещение, при котором обобщенная координата
q = ϕ получит положительное приращение δϕ > 0 , и найдем возможную работу
возмущающей силы F( t ) . Обобщенную силу QqB ( t ) найдем, разделив возможную
работы силы F( t ) на приращение обобщенной координаты δϕ :
=
QqB ( t )
F( t )l cos( ϕ ) ⋅ δϕ
= F0 l sin( pt + β )cos ϕ .
δϕ
48
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
O
l
φ
δφ
А
c
Рисунок 15
Полагая с учетом малости угла
cos ϕ = 1 , получаем выражение для
обобщенной силы в виде:
QqB ( t ) F0 l sin( pt + β )
=
2. Кинематическое возмущение. Пусть вынужденные колебания возникают в
результате
задаваемого
извне
перемещения
точки
крепления
пружины
=
S( t ) S0 sin( pt + β ) (рисунок 16).
Изменение
потенциальной
энергии
пружины
при
перемещении
ее
закрепленного конца с учетом малости угла ϕ будет равно:
1 
1
2
2
2
− c  2lϕ S( t ) − ( S( t ))  .
П′ =
c ( lϕ − S( t )) − ( lϕ )  =


2 
2 
Тогда
∂П ′
−
=
QqB ( t ) =
clS( t ) =
clS0 sin( pt + β ) .
∂ϕ
49
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
O
φ
l
S(t)
А
c
Рисунок 16
3. Инерционное возмущение. Рассмотрим два возможных случая.
Вынужденные относительные колебания (рисунок 17). Пусть маятник
находится на подвижном основании, которое перемещается поступательно и
S( t ) S0 sin( pt + β ) . При составлении дифференциального
=
прямолинейно по закону
уравнения относительных вынужденных колебаний необходимо добавить к реально
действующим
силам
переносную
силу
инерции
Feu = −mae ,
направленную
противоположно переносному ускорению ae = S( t ) . Обобщенную силу QqB ( t )
найдем по аналогии со случаем силового возмущения как отношение возможной
работы переносной силы инерции к приращению обобщенной координаты:
Feи l cos( ϕ ) ⋅ δϕ
B
Qq ( t ) =
δϕ
 t )l cos( ϕ ) ⋅ δϕ
−mS(
=
=
−mp 2 S0 l sin( pt + β )cos ϕ ,
δϕ
или с учетом малости угла ϕ
QqB ( t ) =
−mp 2 S0 l sin( pt + β )
50
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
O
φ
l
А
c
s(t)
Рисунок 17
Вынужденные колебания, вызываемые вращающимся эксцентриком (рисунок
18). Пусть с движущимся поступательно телом массой M скреплен эксцентрик,
имеющий массу m , эксцентриситет l и вращающийся с постоянной угловой
скоростью p . Обозначив через ϕ угол отклонения эксцентрика от вертикали,
выразим QqB через проекцию на горизонталь центробежной силы mp 2l :
QqB ( t )
=
mp 2l sin( pt + β )δ x
= mp 2l sin( pt + β ) .
δx
Из полученных выражений следует, что при инерционном возмущении
колебаний в отличие от силового и кинематического возбуждения амплитуда
обобщенной силы QqB ( t ) пропорциональна квадрату частоты силы p 2 . Полагая
частоту возмущающей силы постоянной p = const , выражение для составляющей
обобщенной силы, обусловленной действием возмущающих сил, независимо от
способа возмущения, можно привести к виду:
51
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
QqB ( t ) Q0 sin( pt + β ),
где Q0 – некоторый постоянный коэффициент.
m
l
с
φ(t)
x
Рисунок 18
2.5.2 Вынужденные колебания при отсутствии сопротивления
При отсутствии сил сопротивления дифференциальное уравнение движения
системы при гармоническом возбуждении в соответствии с (2.27) примет вид:
1
a
q + k 2 q =
QqВ ,
(2.56)
или
q + k 2 q =
h0 sin( pt + β ),
где h0 =
(2.57)
Q0
.
a
Общее решение линейного неоднородного уравнения (2.57) ищем в виде:
52
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
q = q + q ∗ ,
(2.58)
где q – общее решение соответствующего однородного дифференциального
уравнения вида (2.28);
q ∗ – частное решение исходного неоднородного дифференциального уравнения
(2.57), определяемое видом его правой части.
Общее решение однородного уравнения в соответствии с (2.31) и (2.34) имеет
вид
q C1 cos kt + C2 sin kt ,
=
(2.59)
q A ( sin kt + α ) .
=
(2.60)
или
Частное решение q ∗ зависит от соотношения частоты возмущающей силы p
и частоты собственных колебаний k . Рассмотрим два возможных случая.
1. Отсутствие резонанса p ≠ k .
Частное решение исходного неоднородного дифференциального уравнения
(2.57) будем искать в виде:
=
q∗ H sin( pt + β ) ,
(2.61)
где H – постоянный коэффициент.
Продифференцировав (2.61) дважды по времени, получим
q∗ = − Hp 2 sin( pt + β ) .
(2.62)
Подставляя (2.61) и (2.62) в уравнение (2.57), получим
53
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
β ) h0 sin( pt + β ),
− Hp 2 sin( pt + β ) + Hk 2 sin( pt +=
(
)
H k 2 − p 2 sin( pt +=
β ) h0 sin( pt + β ).
Откуда следует
(
)
H k 2 − p2 =
h0 ,
или
H=
(k
h0
2
− p2
)
(2.63)
.
B соответствии с (2.59), (2.62) и (2.63) общее решение (2.58) будет иметь вид:
q = C1 cos kt + C2 sin kt +
h0
k − p2
2
sin( pt + β ) ,
(2.64)
а закон изменения обобщенной скорости по времени:
q =
−C1k sin kt + C2 k cos kt +
h0 p
k 2 − p2
cos( pt + β ) ,
(2.65)
или с учетом (2.60)
q A sin( kt + α ) +
=
h0
k 2 − p2
sin( pt + β ) .
(2.66)
Произвольные постоянные C1 , C2 получим, подставив начальные условия:
=
t0 0,=
q( 0 ) q0=
, q( 0 ) q0 в уравнения (2.64) и (2.65):
54
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
h
q0 =
C1 + 2 0 2 sin β ;
k −p
h p
q0 =
C2 k + 2 0 2 cos β .
k −p
Отсюда
h
C1 =
q0 − 2 0 2 sin β ;
k −p
q0
h0 p
C2 =
cos β .
−
k k k 2 − p2
(
)
Постоянные A и α , входящие в (2.66), при необходимости можно определить
через C1 , C2 по формулам:
A =+
C12 C22 ;
α=
arctg
C1
.
C2
Как следует из (2.65) и (2.66), в рассматриваемом случае система совершает
сложное движение, состоящее из двух гармонических колебаний с частотами k и p
соответственно. Первое слагаемое с частотой k можно по аналогии со случаем
отсутствия возмущающей силы условно назвать свободными колебаниями, а второе
с частотой p – вынужденными колебаниями. Условность названия «свободные
колебания» связана с тем, что определяющие их произвольные постоянные зависят
q( 0 ) q=
не только от начальных условий=
0 , q( 0 ) q0 , но и от параметров
возмущающей силы h0 , p1 , β .
Из-за наличия в реальных системах сил сопротивления, свободные колебания
с частотой k с течением времени затухают, и в дальнейшем система совершает
вынужденные колебания с частотой p , уравнение которых имеет вид:
q
=
h0
k − p2
2
sin( pt + β ) .
55
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Отметим особо, что параметры этого движения не зависят от начальных
условий. Если p < k , то установившиеся вынужденные колебания будут совпадать
по фазе с возмущающей силой, если же p > k , то вынужденные колебания будут
находиться с ней в противофазе (сдвинуты по фазе на π по отношению к
возмущающей силе).
Введем в рассмотрение амплитуду вынужденных колебаний Aвын :
Aвын =
h0
k 2 − p2
.
Тогда уравнение установившихся вынужденных колебаний можно будет
представить в виде:
=
q Aвын sin( pt + β − γ ) ,
где γ – сдвиг по фазе вынужденных колебаний от колебаний возмущающей силы,
причем
0
γ =
π
npu p < k;
npu p > k.
Представим Aвын в виде
Aвын
=
где Aстат
=
h0
1
1
,
А
А
=
=
стат
стат
2
2
k 2 − p2
p
1− z
1− 2
k
h0 Q0
– статическое смещение системы от положения равновесия под
=
c
k2
действием постоянной силы, совпадающей по величине с ее амплитудой,
56
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
z=
p
k
– относительная частота возмущающей силы, называемая также
коэффициентом расстройки.
Величину
η=
1
1 − z2
называют
коэффициентом
динамичности.
Коэффициент динамичности показывает во сколько раз амплитуда вынужденных
колебаний Aвын при гармоническом воздействии больше статического смещения
системы Aстат от действия постоянной силы Q0 .
2. Резонанс p = k . В случае совпадения частоты возмущающей силы с
частотой
свободных
колебаний
(собственной
частотой)
возникает
явление
резонанса, отличительной особенностью которого является нарастание амплитуды
вынужденных колебаний с течением времени.
Будем искать частное решение уравнения (2.57) при p = k в виде
=
q∗ Ht cos( pt + β ) .
(2.67)
Продифференцировав (2.67) дважды по времени, получим
q2 =
− Hp sin( pt + β ) − Htp 2 cos( pt + β ) .
(2.68)
После подстановки (2.67) и (2.68) в (2.57), получим
−2Hp sin( pt +=
β ) h0 sin( pt + β ) ,
H= −
q∗ = −
h0
,
2p
h0 t
ht
π
cos( pt + β ) = 0 sin( pt + β − ) .
2p
2p
2
(2.69)
57
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Как следует из (2.69) вынужденные колебания при резонансе происходят с
нарастающей пропорционально времени амплитудой и запаздывают по фазе от
вынуждающей силы на
π
2
.
В реальной колебательной системе всегда присутствуют силы сопротивления,
а при достижении больших размахов колебаний нарушается допущение о малости
колебаний, поэтому амплитуда колебаний при резонансе хотя и достигает больших
значений, но все же не является неограниченно возрастающей.
Резонанс может стать причиной возникновения опасных напряжений в
материале конструкций и, следовательно, причиной их разрушения или сокращения
срока эксплуатации.
2.5.3 Вынужденные колебания при наличии линейного сопротивления
При наличии линейного сопротивления дифференциальное уравнение (2.27)
при гармоническом возбуждении примет вид:
2
q + 2nq
 + k=
q h0 sin( pt + β ) .
(2.70)
Ограничимся рассмотрением случая отсутствия резонанса, т.е. будем считать,
что n ≠ p . Общее решение уравнения (2.70) будем искать в виде суммы общего
решения q соответствующего однородного уравнения (2.40) и частного решения
q ∗ исходного неоднородного уравнения:
q = q + q ∗ .
(2.71)
Как было показано в подразделе 2.4 общее решение однородного уравнения в
зависимости от соотношения между частотами n и k может иметь вид:
58
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
q e − nt ( C1 cos k ∗t + C2 sin k ∗t ) при n < k ;
=
q e − nt ( C1 + C2t ) при n = k ;
q e − nt  C1et
=

n2 −k 2
+ C2 e − t
n2 −k 2
(2.72)
(2.73)
 при n > k ;


(2.74)
где C1 , C2 – постоянные интегрирования.
Частное решение исходного неоднородного дифференциального уравнения
(2.70) будем искать в виде (2.61):
=
q∗ H sin( pt + β − ε ) .
(2.75)
Продифференцировав (2.75) дважды по времени, получим
q∗ = Hp cos( pt + β − ε ),
q∗ = − Hp 2 sin( pt + β − ε ) .
(2.76)
Подставляя (2.75) и (2.76) в уравнение (2.70), получим
− Hp 2 sin( pt + β − ε ) + 2nHp cos( pt + β − ε ) + Hk 2 sin( pt + β=
− ε ) h0 sin( pt + β ).
Преобразуем правую часть полученного уравнения:
h0 sin( pt=
+ β ) h0 sin( pt + β − ε=
+ε )
h0 sin( pt + β − ε )cos ε + h0 cos( pt + β − ε ) sin ε .
После подстановки в предыдущее уравнение и группировки слагаемых
получим уравнение вида:
59
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
( H (k
2
)
)
− p 2 − h0 cos ε sin( pt + β − ε ) + ( 2nHp − h0 sin ε ) cos( pt + β − ε ) =
0.
Приравнивая коэффициенты к нулю, придем к соотношениям:
(
)
H=
k 2 − p 2 h0 cos=
ε , 2nHp h0 sin ε .
Откуда
Н=
(k
h0
2
−p
tgε =
2
) + 4n
;
2 2
(2.77)
p
2np
;
k 2 − p2
(2.78)
Подставляя найденные значения Н и ε , получаем частное решение в виде:
q∗
=
(k
h
2
−p
2
)

 2np  
sin  pt + β − artg  2
2 

.
2 2
k
p
−


+ 4n p

(2.79)
Общее решение уравнения (3.70) будет иметь вид:
1) при n < k
=
q e − nt ( C1 cos k ∗t + C2 sin k ∗t ) +
(k
h
2
−p
2
) + 4n
2 2
sin ( pt + β − ε ) ;
(2.80)
p
2) при n = k
=
q e − nt ( C1 + C2 t ) +
(k
h
2
−p
2
) + 4n
2 2
sin ( pt + β − ε ) ;
(2.81)
p
3) при n > k
60
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
q e − nt  C1et
=

n2 −k 2
+ C2 e − t
n2 −k 2

+

(k
h
2
−p
2
) + 4n
2 2
sin ( pt + β − ε ) ; (2.82)
p
где ε – величина, определяемая выражением (2.78);
C1 , C2 – постоянные интегрирования, определяемые из начальных условий.
Структура полученных общих решений (2.80) – (2.72) показывает, что система
совершает сложное движение, но т.к. первое слагаемое с увеличением времени t
стремится
к
нулю,
то
определяющим
становится
последнее
слагаемое,
соответствующее q∗ . В этом случае система будет совершать установившиеся
вынужденные колебания.
Полученные результаты позволяют сформулировать основные свойства
установившихся вынужденных колебаний при гармоническом возбуждении.
Во-первых, вынужденные колебания, даже при наличии сил сопротивления,
являются незатухающими и длятся до тех пор, пока действует возмущающая сила.
Во-вторых, амплитуда и частота вынужденных колебаний не зависят от начальных
условий. Наличие сил сопротивления движению приводит к уменьшению
амплитуды
вынужденных
колебаний.
В-третьих,
вынужденные
колебания
совпадают по частоте и отстают по фазе от возмущающей силы на величину
0 ≤ε ≤π .
В заключение отметим, что свободное движение, затухающее по истечении
некоторого промежутка времени, возникает каждый раз, когда изменяются
параметры возмущающей силы (амплитуда, частота и т.д.). При этом система
осуществляет плавный переходный процесс от одного установившегося режима
вынужденных колебаний к другому. Для того, чтобы получить уравнение движение
системы во время переходного процесса, требуется определить постоянные
интегрирования в общем решении (2.80) – (2.82) при начальных условиях, равных
значениям отклонения и скорости в момент изменения параметров возмущающей
силы.
61
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2.6 Вынужденные колебания в случае периодической возмущающей силы
Рассмотрим случай действия возмущающей силы, являющейся некоторой
периодической функцией времени с периодом τ (рисунок 19). Подобные силы часто
встречаются при решении самых различных технических задач.
Q(t)
t
τ
Q(t)
t
τ
Рисунок 19
Воспользуемся способом определения периодического решения с помощью ряда
Фурье, который целесообразно применять в тех случаях, когда ряд быстро сходится и
можно ограничиться несколькими первыми членами.
Пусть функция QqВ ( t ) удовлетворяет условиям Дирихле, т. е. она
ограничена, имеет разрывы первого рода и конечное число экстремумов на
конечном интервале. Тогда QqВ ( t ) можно разложить в ряд Фурье:
∞
(
)
QqB ( t ) =
Q0 + ∑ a j cos ( jpt ) + b j sin ( jpt ) ,
j =1
(2.83)
62
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где p = 2π – основная частота возмущающей силы.
τ
Постоянные коэффициенты Q0 , a j , =
b j ( j 1,2,3,...,∞ ) , входящие в (2.83),
равны:
τ
τ
τ
1
2
2
B
B
=
Q0 = ∫ QqB ( t )dt; an =
∫ Qq ( t )cos( jpt )dt; bn ∫ Qq ( t ) sin( jpt )dt.
τ
τ
0
τ
0
0
Тригонометрический ряд (2.83) можно представить в ином виде:
∞
B
Qq ( t ) =
Q0 + Q j
j =1
∑
(
)
sin jpt + β j ,
arctg
где введены обозначения Q j =a 2j + b 2j ; β j =
aj
bj
(2.84)
.
Отдельные члены ряда (2.83) или (2.84) называются гармониками j-го
порядка.
Дифференциальное
уравнение
вынужденных
колебаний
(2.70)
в
рассматриваемом случае примет вид:
∞
q + 2nq
 + k q =h0 + ∑ h j sin ( jpt + β j ) ,
2
(2.85)
n =1
где
=
h0
Qj
Q0
=
; hj
.
a
a
Общее решение уравнения (2.85) можно представить в виде суммы общего
решения однородного уравнения (2.40), характеризующего свободное движение
системы,
и
частного
решения
исходного
неоднородного
уравнения,
характеризующего вынужденные колебания. При установившемся движении
63
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
практическое значение имеют только незатухающие вынужденные колебания,
определяемые уравнением
=
q
∞
∑ A j sin ( jpt + β j − ε j ) ,
(2.86)
j =1
где A j и
(β j − ε j )
– соответственно амплитуда и сдвиг фазы вынужденных
колебаний, соответствующие частоте=
p j jp
=
( j 1,2,3,...,∞ ) и определяемые
соотношениями
=
Aj
hj
2njp
=
; tgε j
.
2
2 2
2
2 2
2 2 2
−
k
j
p
k − j p + 4n j p
(
)
(2.87)
Если частота одной из гармоник возмущающей силы совпадает с частотой
свободных колебаний системы, т.е.
=
jp k=
( j 1,2,3,...,∞ ) ,
то говорят о наступлении резонанса j-го порядка, а соответствующие частоты
называются критическими. Как следует из (2.87), амплитуда j-той гармоники и
сдвиг фаз получаются равными
Aj
=
hj
π
.
; εj
=
2nk
2
Стоит заметить, что вследствие сдвига фаз амплитуда A j имеет максимум не
при резонансе p=
jp
= k , а при другом значении частоты p j .
j
Исследовав выражение для амплитуды гармоники (2.87) на наличие
64
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
максимума, получаем, что при n < k
2
2
максимума
Aj
не существует, и,
следовательно, резонанса в этом случае быть не может.
2.7 Вынужденные колебания при произвольной возбуждающей силе
Рассмотрим движение механической системы с одной степенью свободы в
случае действия возбуждающей силы, меняющейся во времени по произвольному
закону.
При отсутствии сил сопротивления дифференциальное уравнение движения
системы (2.27) примет вид
1
a
q + k 2 q =
QqВ ( t ),
(2.88)
Воспользуемся методом вариации произвольных постоянных и представим
общее решение дифференциального уравнения (2.88) в виде
=
q C1( t )cos kt + C2 ( t ) sin kt ,
(2.89)
где C1 ( t ) и C2 ( t ) – некоторые функции времени.
Продифференцировав (2.89) по времени, получим
q =
dC1 ( t )
dC ( t )
cos kt + 2
sin kt − C1 ( t )k sin kt + C2 ( t )k cos kt .
dt
dt
Т.к. мы ввели две неизвестные функции, то можно связать их некоторым
условием, например,
65
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
dC1 ( t )
dC ( t )
cos kt + 2
sin kt =
0.
dt
dt
(2.90)
При этом условии выражение для обобщенной скорости будет иметь такой же
вид, как и в случае постоянных коэффициентов C1 и C2 , т.е.
q =
−C1 ( t )k sin kt + C2 ( t )k cos kt .
(2.91)
Продифференцируем (2.91) по времени:

q =
k − 1
dC ( t )
dC ( t )

sin kt + 2
cos kt  − k 2 ( C1 ( t )cos kt + C2 ( t ) sin kt ) .
dt
dt


(2.92)
Подставив (2.89) и (2.92) в уравнение (2.88), имеем
−
dC1 ( t )
dC ( t )
1 B
sin kt + 2
cos kt =
Qq ( t ) .
dt
dt
ak
Решая совместно уравнения (2.90) и (2.93) относительно
(2.93)
dC1 ( t )
dC2 ( t )
и
,
dt
dt
получим
dC1 ( t )
1
dC2 ( t ) 1 B
=
− QqB ( t ) sin kt;
=
Qq ( t )cos kt .
dt
ak
dt
ak
Откуда
t
1
C1( t ) =
D1 − ∫ QqB ( τ ) sin kτ dτ ;
ak 0
t
1
C2 ( t ) =
D2 +
QqB ( τ )cos kτ dτ ,
∫
ak 0
(2.94)
66
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где τ – переменная интегрирования.
Подставляя
найденные
значения
(2.94)
в
решение
(2.89),
после
преобразований получим
τ
1
q = D1 cos kt + D2 sin kt +
QqB ( t ) sin k ( t − τ )dτ ,
∫
ak 0
(2.95)
где D1 и D2 – постоянные интегрирования, определяемые из начальных условий
=
t0 0,=
q( 0 ) q0=
, q( 0 ) q0 и не зависящие от параметров возмущающей
движения
силы.
=
C1( 0 ) D=
Из (2.94) следует, что
1 ; C2 ( 0 ) D2 , тогда подставляя начальные
условия в выражения (2.89) и (2.94) получим
=
D1 q=
0 ; D2
q0
.
k
С учетом найденных значений постоянных D1 , D2 общее решение (2.95)
примет вид
τ
q0
1
q = q0 cos kt + sin kt +
QqB ( t ) sin k ( t − τ )dτ .
∫
k
ak 0
(2.96)
При наличии вязкого сопротивления дифференциальное уравнение движения
(2.27) имеет вид:
1
а
q + 2nq
 + k 2q =
QqB ( t ) .
(2.97)
Ограничимся наиболее важным для инженерной практики случаем малого
67
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
сопротивления ( n < k ) . Представим решение этого уравнения в виде
(
)
=
q e − nt C1 ( t )cos k ∗t + C2 ( t ) sin k ∗t ,
(2.98)
где C1 ( t ) и C2 ( t ) – некоторые функции времени;
=
k∗
k 2 − n2 .
Повторяя
приведенные
выше
рассуждения,
после
преобразований
и
подстановки начальных условий найдем
q =e
− nt
τ
q0 + nq0
1

∗
∗ 
+
+
q
cos
k
t
sin
k
t
Q B ( t ) sin k ∗ ( t − τ )dτ .
 0

∗
∗∫ q
k

 ak 0
(2.99)
Нетрудно установить, что при n = 0 выражение (2.99) приводится к виду
(2.96).
Полученные результаты можно применять для исследования вынужденных
колебаний систем с одной степенью свободы при действии каких угодно
возмущающих сил, особенно в случае возмущающих сил, резко меняющих величину
и направление, а также заданных графически или численно.
2.8 Основы теории регистрирующих приборов
Рассмотрим принципы работы приборов, предназначенных для регистрации
различных переменных во времени физических величин. Ограничимся случаем,
когда
измеряемая
величина
является
периодической
или,
в
частности,
гармонической.
Практически все механические приборы в качестве одного из основных
упругих элементов имеют пружину того или иного типа, деформация которой,
вызываемая воздействием измеряемой величины и определяет показания прибора.
68
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Рассмотрим два типа приборов.
Первый тип приборов предназначен для непосредственного измерения
некоторой физической величины. Простейшая схема такого прибора представлена
на рисунке 20, где F – сила, определяемая измеряемой величиной х; x1 = f(x) –
показания прибора, пропорциональные деформации пружины под действием силы
F. Выясним, какие ограничения должны накладываться на измеряемые процессы
для того, чтобы показания прибора соответствовали измеряемым величинам, т.е.
выполнялось условие x1(t) = Ax(t), где А – некоторая постоянная величина.
x
F
x1
Рисунок 20
Рассмотрим прибор как некоторую колебательную систему с одной степенью
свободы, подверженную силовому воздействию измеряемой величины. Вследствие
того, что стрелка прибора и некоторые другие его подвижные части обладают
массой, необходимо учитывать инерцию подвижных частей и возникающие при
движении силы вязкого сопротивления.
Величину Д, равную отношению частоты собственных колебаний k к
удвоенному коэффициенту затухания n , называют добротностью колебательной
системы:
69
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Д=
Величина,
обратная
k
.
2n
добротности,
носит
название
безразмерного
коэффициента затухания:
d=
2n
.
k
Добротность Д представляет собой коэффициент динамичности при резонансе и
показывает, во сколько раз амплитуда вынужденных колебаний Aвын при резонансе
превышает статическое смещение Астат . Как указывалось ранее в пункте 2.5.2 в
реальных колебательных системах при наличии вязкого сопротивления амплитуда при
резонансе имеет конечное значение. Если частота изменения возмущающей силы p
значительно меньше частоты свободных колебаний k , т. е. p  k , то амплитуда
вынужденных колебаний Aвын близка к статическому смещению Астат , а коэффициент
динамичности η близок к единице. Если же частота изменения возмущающей силы p
значительно больше частоты свободных колебаний k , т. е. p  k , то колебательная
система ведет себя как фильтр, практически не воспринимая возмущения с частотами,
существенно превышающими собственную частоту.
Для того чтобы исключить искажения показаний прибора, необходимо
предусматривать как можно большую его добротность и, как следствие, как можно
меньший безразмерный коэффициент затухания.
Однако если использовать приборы с большой добротностью, то возникающие
свободные колебания будут затухать очень медленно и в значительной степени
искажать показания прибора. Поэтому на практике приходится выбирать достаточно
малые значения добротности, близкие к единице.
Второй тип приборов, используется для замера вибраций подвижных
объектов (железнодорожных вагонов, автомобилей, самолетов и т. д.). Эти приборы
70
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
иногда называют сейсмическими, поскольку они построены по такому же принципу,
что и сейсмографы – приборы для регистрации колебаний земной коры.
Такие приборы должны измерять абсолютные колебания изучаемого объекта
относительно некоторой реально несуществующей «неподвижной» системы
отсчета, не имея возможности опереться на нее. Следовательно, в приборе должно
быть тело, покоящееся относительно «неподвижной» системы отсчета.
Таким телом отсчета может служить тело 1, упруго соединенное с
исследуемым объектом 2 (рисунок 21).
S(t)
h
с
1
m
2
S1(t)
Рисунок 21
В этом случае коэффициент динамичности η ≈ 0 и собственная частота
прибора k существенно меньше частоты низшей регистрируемой гармоники, т. е.
тело 1 остается практически неподвижным, и запись движения объекта относительно тела с точностью до знака будет соответствовать колебаниям объекта.
Таким образом, в сейсмических приборах соотношение между собственной
частотой и частотами регистрируемых гармоник обратно тому, что имеет место в
приборах первого типа.
71
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2.9 Основы виброзащиты
Проблемы, связанные с защитой от негативного воздействия вибраций
(виброзащитой), рассмотрим на простейших примерах системы с одной степенью
свободы, представляющей собой тело массой m , совершающее вынужденные
прямолинейные колебания.
Необходимость виброзащиты возникает в двух случаях:
1. Необходимость уменьшения воздействия вибраций, возникающих в
процессе работы, на фундамент какой-либо машины.
2. Необходимость защиты какого-либо устройства (экипаж, прибор и т. д.) от
вредного
воздействия
вибраций,
возникающих
при
транспортировке
или
являющихся результатом работы находящихся рядом машин.
В первом случае возмущающая сила приложена непосредственно к телу
(рисунок 22, а) и имеет место силовое возбуждение колебаний, во втором (рисунок
22, б) – кинематическое возбуждение колебаний вследствие вибрации основания.
Воздействие на фундамент, при силовом возбуждении по возможности
должно быть малым. Станки, двигатели или другие аналогичные машины,
создающие в процессе работы вибрации, устанавливают на фундамент, как правило,
через различные амортизирующие подушки, изготовленные из специальных сортов
резины,
обладающей
помимо
упругости
большим
внутренним
неупругим
сопротивлением. Эти подушки можно аппроксимировать пружиной жесткостью c и
демпфером с коэффициентом вязкого сопротивления µ (рисунок 22, а). Тогда
динамическое воздействие на фундамент будет равно
R(=
t ) cxc + µ xc ,
(2.100)
где xc – перемещение тела относительно положения статического равновесия.
При кинематическом возбуждении колебаний необходимо, чтобы абсолютные
перемещения xk тела были как можно более малыми (рисунок 22, б).
72
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
При силовом возбуждении, сила тяжести компенсируется статической
деформацией пружины, и уравнение движения тела имеет вид:
mxc + µ xc + cxc =
F (t ) ,
(2.101)
а при кинематическом возбуждении
 t )) + c( x − s( t )) =
mx + h( x − s(
0,
k
k
k
или
mxk + µ xk + cxk = cs(t ) + µ s(t ) ,
(2.102)
где µ – коэффициент вязкого сопротивления материала пружины или
специально установленного демпфера (например, амортизатора);
c – жесткость пружины подвески.
m
m
с
с
а)
б)
Рисунок 22
73
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Обозначим отношение амплитуд динамического R(t ) и силового F (t )
воздействий на фундамент через β c , а отношение амплитуд xk (t ) и s(t ) при
кинематическом воздействии β k . Как показано в [3, 4], условие виброзащиты не
зависит от природы возмущения и имеет вид:
k 4 + 4n 2 p 2
=
β β=
β=
c
k
( k − p ) + 4n p
2
2 2
2 2
< 1.
(2.103)
Разделив числитель и знаменатель (2.103) на k 2 , преобразуем условие виброзащиты к виду:
=
β
1 + d 2z2
(1− z ) + d z
2 2
2 2
< 1,
где z = p k – коэффициент расстройки;
d = 2n k – безразмерный коэффициент затухания.
Отсюда, имеем:
( 1 − z 2 )2 + d 2 z 2 > 1 + d 2 z 2 ,
или
z2( z2 − 2 ) > 0 .
Это неравенство выполняется при z > 2 или p > k 2 .
Следовательно, независимо от способа возмущения и значения вязкого
сопротивления для осуществления виброзащиты необходимо, чтобы собственная
74
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
частота колебаний системы была значительно ниже (по крайней мере, в
2 раза)
частоты возбуждения.
При p > k 2 демпфирование будет играть отрицательную роль, поскольку,
чем меньше демпфирование, тем больше эффект виброзащиты. Вместе с тем
уменьшение
демпфирования
не
всегда
приводит
к
увеличению
эффекта
виброзащиты.
Необходимо учитывать, что при силовом возбуждении машина при пуске
проходит режим раскрутки, а при остановке – торможения. Частота возбуждения
при этом изменяется от нуля до p и наоборот. При этом система проходит через
резонанс,
что
вынуждает
конструктора
в
ущерб
виброзащите
повышать
демпфирование, чтобы уменьшить амплитуду резонансных колебаний. При
кинематическом возбуждении возможно скачкообразное перемещение основания
(наезд на препятствие, попадание колеса в яму и т. д.), что при отсутствии
амортизаторов может привести к недопустимым перемещениям и, следовательно,
перегрузкам защищаемого объекта, а также к длительному процессу затухания
возникающих свободных колебаний.
75
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
3 Малые колебания системы с двумя степенями свободы
3.1 Дифференциальное уравнение малых колебаний системы с двумя
степенями свободы
Рассмотрим
механическую
систему
c
голономными,
стационарными,
удерживающими связями, состоящую из n материальных точек и имеющую две
степени свободы. Предположим, что система имеет устойчивое положение
равновесия, от которого будем отсчитывать обобщенные координаты q1 и q2 .
Выведем систему из этого положения, сообщив ее точкам малые начальные
возмущения (малые отклонения от положения равновесия и скорости). Вследствие
устойчивости
положения
равновесия
при
дальнейшем
движении
системы
отклонения и скорости всех ее точек будут также малы по модулю, т.е. система
будет совершать движение вблизи этого положения равновесия.
Дифференциальные уравнения движения механической системы с двумя
степенями свободы в обобщенных координатах (уравнения Лагранжа II рода)
имеют вид:
 d  ∂T  ∂T
=
Q1 ;
 
−
 dt  ∂q1  ∂q1

 d  ∂T  − ∂T =
Q2 ,
 dt  ∂q  ∂q
2
  2
(3.1)
где T – кинетическая энергия системы;
Q i – обобщенная сила, соответствующая выбранной обобщенной координате
qi ( i = 1,2 ) .
С учетом (2.13), уравнения (3.1) примут вид:
76
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
 d  ∂T  ∂T
−
=
Q1П + Q1R + Q1В ;
 

 dt  ∂q1  ∂q1

 d  ∂T  − ∂T =
Q2П + Q2R + Q2В ,


 dt ∂q
  2  ∂q2
(3.2)
где QqП – потенциальная составляющая обобщенной силы;
QqR – диссипативная составляющая обобщенной силы;
QqВ – вынуждающая составляющая обобщенной силы соответственно.
Кинетическая энергия механической системы материальных точек в общем
случае определяется выражением (2.6):
T=
1 n
mk vk2 ,
∑
2 k =1
где vk – скорость точки массой mk .
Все связи системы являются стационарными, следовательно, радиус-векторы
точек системы зависят только от обобщенных координат q1 ( t ), q2 ( t ) , которые в
свою очередь являются функциями времени, т.е. ri = ri ( q1 ( t ),q2 ( t )) . Тогда,
vk =
drk ∂rk dq1 ∂rk dq2 ∂rk
∂r
=
⋅
+
⋅
=
q1 + k q2 .
∂q1
∂q2
dt ∂q1 dt ∂q2 dt
(3.3)
Подставляя (3.3) в выражение кинетической энергии (2.6) с учетом тождества
( vk )2 =
vk2 , получим
2
 ∂rk

∂rk
1 n
1 n
2


T= ∑ mk vk
m
q
q
=
+
∑ k  ∂q 1 ∂q 2  ,
2
2
2
 1

=
k 1=
k 1
77
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2
  ∂r  2
 ∂rk   ∂rk 
 ∂rk  2 
1 n
2
k
T = ∑ mk 
 q1 + 2 

 q1q2 + 
 q2  ,
2 k =1  ∂q1 
q
q
q
∂
∂
∂

 2
 1  2 

2
2
n
n
 ∂rk  2
 ∂rk   ∂rk 
 ∂rk  2 
1 n
 ∑ mk 
T=
 q1 + 2 ∑ mk 

 q1q2 + ∑ mk 
 q2  .
2  k 1=
∂q1 
∂q=
∂q2 
∂q2 

1





=
k
1
k
1

Введем обозначение
n
 ∂r   ∂r 
Aij = ∑ mk  k   k  .
 ∂qi   ∂q j 
k =1
(3.4)
тогда выражение для кинетической энергии системы примет вид:
=
T
1
q12 + 2 A12q1q2 + A22q22  ,
A
11

2
(3.5)
где Aij являются, в общем случае, функциями обобщенных координат q1 ( t ), q2 ( t ) .
Разложим Aij в ряд Маклорена по степеням обобщенных координат q1 , q2 в
окрестности положения равновесия, приняв в положении равновесия q=
q=
1
2 0.
 ∂Aij 
 ∂Aij 
+
+
q
 1 
 q2 + ...
0
∂
q
∂
q
 1 0
 2 0
( )
Aij ( q1 ,q2 ) =
Aij
Здесь, как и ранее, обозначение
( ...)0
(3.6)
показывает, что соответствующее
выражение вычисляется в положении равновесия при q=
1 q=
2 0.
78
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Т.к. в выражении для кинетической энергии (3.5) следует ограничиться
слагаемыми не выше второго порядка малости, то с требуемой степенью точности
можно считать, что
Aij ( q1 ,q2=
)
( Aij )=0
a=
const .
ij
Величины aij называются обобщенными коэффициентами инерции, причем
согласно (3.4) aij = a ji .
Окончательно
для
системы
с
двумя
степенями
свободы
выражение
кинетической энергии в обобщенных координатах примет вид:
1
T =  a11q12 + 2a12q1q2 + a22q22  ,
2
(3.7)
т.е. кинетическая энергия механической системы является однородной
квадратичной формой обобщенных скоростей.
Так как кинетическая энергия всегда положительна и равняется нулю только
при нулевых значениях обобщенных скоростей, то ее коэффициенты должны
удовлетворять критерию Сильвестра (1.24):
a11 > 0, a22 > 0,
a11
a12
a12
> 0.
a22
Потенциальная энергия системы с двумя степенями свободы в соответствии с
выражением (1.14) имеет вид:
1
П( q1 ,q2 ) =  c11q12 + 2c12 q1 q2 + c22 q22  ,
2
(3.8)
79
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где сij ( i, j = 1,2 ) – обобщенные коэффициенты жесткости, характеризующие
упругие свойства механической системы, причем сij = с ji .
Выражение (3.8) показывает, что потенциальная энергия механической
системы является однородной квадратичной формой обобщенных координат.
Обобщенные коэффициенты жесткости сij в соответствии с (1.24) должны
удовлетворять условию:
c11 > 0, c22 > 0,
c11
c12
c12
> 0.
c22
3.2 Свободные колебания системы с двумя степенями свободы
3.2.1 Дифференциальные уравнения свободных колебаний
В случае отсутствия диссипативных и возмущающих сил, уравнения (3.2)
примут вид
 d  ∂T  ∂T
−
=
Q1П ;
 

 dt  ∂q1  ∂q1

 d  ∂T  − ∂T =
Q2П .


 dt ∂q
  2  ∂q2
(3.9)
Составляющая обобщенной силы от действия потенциальных сил равна взятой с
противоположным знаком частной производной от потенциальной энергии по
соответствующей обобщенной координате:
∂П
QiП =
−
(i =
1,2 ),
∂qi
(3.10)
80
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где П( q1 ,q2 ) – потенциальная энергия системы, определяемая выражением (3.8).
Тогда,
Q1П =
− ( c11q1 + c12 q2 ) , Q2П =
− ( c12 q1 + c22 q2 ) .
(3.11)
С учетом выражения (3.7) для кинетической энергии системы, получим
∂T ∂T
= = 0;
∂q1 ∂q2
∂T
=
a11q1 + a12 q2 ;
∂q1
d  ∂T 
a11q1 + a12 q2 ;

=
dt  ∂q1 
∂T
=
a12 q1 + a22 q2 ;
∂q2
d  ∂T 
a12 q1 + a22 q2 .

=
dt  ∂q2 
(3.12)
Подставляя выражения (3.11) и (3.12) в уравнения (3.9), получим уравнения
движения системы в случае свободных колебаний без учета сил сопротивления:
a11q1 + a12 q2 + c11q1 + c12 q2 = 0;

a12 q1 + a22 q2 + c12 q1 + c22 q2 = 0.
(3.13)
Коэффициенты a12 и c12 называют коэффициентами инерционной и упругой
связи соответственно. Если в колебательной системе коэффициент a12 = 0 , ее
называют системой с упругой связью, а если c12 = 0 – системой с инерционной связью.
3.2.2 Общее решение дифференциальных уравнений свободных колебаний
системы с двумя степенями свободы
Решение системы уравнений (3.13) будем искать в виде q1 = C1eλ t , q2 = C2 eλ t .
В случае малых колебаний в окрестности устойчивого равновесия все значения λ
81
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
мнимые и, следовательно, частные решения можно представить в следующем виде:
q1 = A1 sin( kt + α ); q2 = A2 sin( kt + α ) ,
(3.14)
где A1 , A2 , k , α – неизвестные постоянные, подлежащие определению.
Продифференцируем (3.14) дважды по времени:
q1 =
−k 2 A1 sin( kt + α ); q2 =
−k 2 A2 sin( kt + α )
(3.15)
и подставим полученные выражения в систему уравнений (3.13). После
преобразований получим систему двух алгебраических однородных линейных
уравнений относительно A1 и A2 :
(
(
)
)
(
(
)
)
 c11 − a11k 2 A1 + c12 − a12 k 2 A2 =
0;


2
2
0.
 c12 − a12 k A1 + c22 − a22 k A2 =
(3.16)
Система алгебраических однородных линейных уравнений (3.16) имеет
отличное от нуля решение, если определитель системы равен нулю:
c11 − a11k 2
c12 − a12 k 2
c12 − a12 k
c22 − a22 k
2
2
= 0.
(3.17)
Раскрывая определитель (3.17), получаем уравнение частот:
(c
)(
) (
2
c22 − a22 k 2 − c12 − a12 k 2
11 − a11k
)
2
=
0.
(3.18)
Уравнение частот (3.18) в общем случае имеет два корня k12 и k22 . Если
82
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
квадратичные формы для кинетической и потенциальной энергии удовлетворяют
условиям определенной положительности, то оба корня k12 и k22 являются
действительными и положительными. В этом случае обобщенные координаты q1 и
q2 выражаются синусоидальной зависимостью от времени. В противном случае
движение системы не является колебательным.
Извлекая из k12 и k22 квадратные корни и отбрасывая отрицательные значения
частот, как не имеющие физического смысла, получаем две собственные частоты
колебаний системы k1 и k2 . Меньшую из частот обозначают k1 и называют
низшей, большую k2 называют высшей. Частоты k1 и k2 не зависят от начальных
условий и полностью определяются инерционными и упругими свойствами самой
системы.
В случае равенства частот k=
1 k=
2 k каждое из уравнений системы является
тождеством,
справедливым
при
любых
значениях
A1
и
A2 .
Система
дифференциальных уравнений (3.13) распадается на два независимых уравнения –
одно для q1 , другое для q2 . Их решения имеют вид:
q1 = A1 sin ( kt + α1 ); q2 = A2 sin ( kt + α 2 )
(3.19)
Амплитуды и фазы колебаний A1 , A2 ,α1 , α 2 определяются независимо для
каждой координаты из начальных условий:
)
=
t0 0,=
q1 q1( 0 )=
, q2 q2( 0 )=
, q1 q1( 0 =
, q2 q2( 0 ) .
(3.20)
Система совершает гармонические колебания, при которых обобщенные
координаты q1 и q2 изменяются по синусоидальному закону независимо друг от
друга с одинаковыми частотами.
83
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
В случае разных частот k1 и k2 каждой частоте соответствуют определенные
значения A1 , A2 ,α . Обозначим амплитуды и фазу для частоты k1 через A1( 1 ) , A2( 2 ) , α1
и A1( 2 ) ,A2( 2 ) ,α 2 – для частоты k2 .
В соответствии с этим получим два независимых частных решения,
соответствующих частотам k1 и k2 :
 q1( 1 )
=
 (1)
=
 q2
 (2)
=
 q1
 (2)
=
 q2
A1( 1 ) sin ( k1t + α1 );
A2( 1 ) sin ( k1t + α1 );
(3.21)
A1( 2 ) sin ( k2t + α 2 );
A2( 2 ) sin ( k2t + α 2 );
где q1( 1 ) и q2( 1 ) составляют главные колебания для частоты k1 , а q1( 2 ) и q2( 2 ) – для
частоты k2 .
Каждое из главных колебаний является гармоническим для обеих обобщенных
координат.
Система
однородных
линейных
уравнений
(3.16)
дает
возможность
определить только отношение амплитуд A2 и A1 . Для первого и второго главных
колебаний соответственно получаем:

A2(1)
c11 − a11k12
c12 − a12 k12
=
−
=
−
,
 β1 =
(1)
2
2
−
−
A
c
a
k
c
a
k

12
12 1
22
22 1
1

A2(2)
c11 − a11k22
c12 − a12 k22
β =
=
−
=
−
.
2
2
 2 A(2)
−
−
c
a
k
c
a
k
12
12 2
22
22 2

1
Отношения
амплитуд
в
главных
колебаниях
β1
и
(3.22)
β2 ,
называют
коэффициентами формы:
84
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
A2(1) q2(1)
A2(2) q2(2)
=
β1 =
β2 =
; =
A1(1) q1(1)
A1(2) q1(2)
(3.23)
Общее решение системы уравнений (3.13) получается путем суммирования
частных решений:
 q1 = q1( 1 ) + q1( 2 ) = A1( 1 ) sin(k1t + α1 ) + A1( 2 ) sin(k2t + α 2 ),

(1)
(2)
(1)
(2)
 q2 = q2 + q2 = β1 A1 sin(k1t + α1 ) + β 2 A1 sin(k2t + α 2 ).
(3.24)
Четыре постоянных A1( 1 ) , A1( 1 ) , α1 , α 2 определяются из начальных условий:
)
=
t0 0,=
q1 q1( 0 )=
, q2 q2( 0 )=
, q1 q1( 0 =
, q2 q2( 0 ) .
(3.25)
Полученные результаты показывают, что собственные колебания системы с
двумя степенями свободы не являются простым гармоническим движением и
представляют собой сумму двух главных гармонических колебаний с частотами k1 и
k2 , которые содержатся в каждой обобщенной координате q1 и q2 .
3.2.3 Главные координаты
2
Выберем в качестве новых обобщенных координат системы q1( 1 ) и q1( ) , тогда
главное колебание с частотой k1 будет характеризовать только обобщенная
координата q1( 1 ) , а главное колебание с частотой k2 , – обобщенная координата q1( 2 ) .
Обобщенные координаты, которые совершают только по одному главному
колебанию, называются главными координатами системы. Произвольные
обобщенные координаты можно представить линейными комбинациями главных:
85
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
β1q1( 1 ) + β 2 q1( 2 ) .
q1 =+
q1( 1 ) q1( 2 ) ; q2 =
(3.26)
Очевидно, что в случае равенства частот k1 = k2 , любые обобщенные
координаты будут главными.
Как и в случае равных частот, система уравнений (3.13) для главных
координат должна распадаться на два отдельных, независимых уравнения. Это
означает, что выражения потенциальной (3.7) и кинетической (3.8) энергий не
a12 0,=
c12 0 .
должны содержать членов с произведениями переменных, т.е.=
Найдем из этого условия главные координаты. Обозначим q1 и q2 – произвольные
обобщенные координаты, q1∗ и q2∗ – главные координаты.
Из (3.26) получим
q1 =
q1∗ + q1∗ ; q2 =
β1q1∗ + β 2 q1∗ ,
(3.27)
где β1 , β 2 – коэффициенты формы, определяемые соотношениями (3.22).
Выражения
для
кинетической
и
потенциальной
энергий
в
главных
обобщенных координатах имеют вид:
( )
1
П = с*11 q*1
2
2
( )
+ с*22 q*2
2
( )
1 * *
 ; T = a11 q1
2

2
( )
+ a*22 q*2
2
.

(3.28)
Подставляя (3.27) в выражения потенциальной (3.7) и кинетической (3.8)
энергий, получим:
86
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»

=
П




T
=




(
1
с11 + 2с 12 β1 + с22 β12
2 
(
)( q )
∗ 2
1
+ с11 + 2с 12 β 2 + с22 β 22
(
(
)
+ 2 с11 + с 12 ( β1 + β 2 ) + с22 β1 β 2 q1∗q2∗ +
)( q )  ;
∗ 2
2
)( )
1
a11 + 2a 12 β1 + a22 β12 q1∗

2
2
(3.29)
(
)
+ 2 a11 + a 12 ( β1 + β 2 ) + a22 β1 β 2 q1∗q2∗ +
)( )
(
2
+ a11 + 2a 12 β 2 + a22 β 22 q2∗  .

Сопоставляя выражения (3.28) и (3.29), получим следующие уравнения:
с11 + с ( β1 + β 2 ) + с22 β1 β 2 =
0,
12

0,
a11 + a 12 ( β1 + β 2 ) + a22 β1 β 2 =
(3.30)
и соотношения между старыми и новыми коэффициентами инерции и жесткости:
2
=
a*11 a=
a*22
11 + 2a12 β1 + a22 β1 ;
*
2
=
c11
c=
c*22
11 + 2c12 β1 + c22 β1 ;
a11 + 2a12 β 2 + a22 β 22 ;
c11 + 2c12 β 2 + c22 β 22 .
(3.31)
Решая систему (3.30) относительно коэффициентов формы β1 и β 2 , найдем:
β1 β 2
=
a11c12 − a12 c11
a22 c11 − a11c22
=
; β1 + β 2
.
a12 c22 − a22 c12
a12 c22 − a22 c12
(3.32)
Заметим, что соотношения (3.32) являются корнями квадратного уравнения:
(3.33)
0.
( a11c22 − a22c12 ) β 2 − ( a22 c11 − a11c22 ) β + ( a11c12 − a12c11 ) =
87
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Таким образом, коэффициенты формы β1 и β 2 можно определять двумя
способами: как по формулам (3.22), так и как корни квадратного уравнения (3.33).
Система уравнений малых колебаний (3.13) для главных координат q1∗ и q2∗ с
*
*
учетом a=
c=
12
12 0 распадется на два независимых уравнения:
a11q1∗ + c11q1∗ = 0;

∗
∗
a22 q2 + c22 q2 = 0.
(3.34)
Общее решение системы дифференциальных уравнений (3.34) будет иметь
вид:
 q1∗ A1( 1 ) sin ( k1∗t + α1 ),
=
 ∗
(2)
∗
=
 q2 A2 sin ( k2t + α 2 ),
(3.35)
где k1∗ и k2∗ – частоты главных колебаний, вычисляемые по формулам:
=
k1*
c*11
=
; k*2
*
a11
c*22
.
a*22
(3.36)
Отметим особо, что т.к. частоты главных колебаний являются физическими
характеристиками системы, то они не могут зависеть от выбора обобщенных
координат, следовательно,
*
=
k1* k=
k2 .
1 ; k2
Каждая из главных координат системы изменяется по гармоническому закону
определенной частоты, амплитуды и фазы, как это наблюдается в случае системы с
одной степенью свободы. Полученные результаты будут справедливы и для систем
88
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
с любым конечным числом степеней свободы. Использование главных координат,
не
упрощает
вычислений,
однако
главные
координаты
имеют
большое
теоретическое значение, например, при исследовании вынужденных колебаний
системы.
3.3 Влияние линейного сопротивлении на собственные колебания
Пусть на точки системы действуют линейные силы сопротивления Rk ,
пропорциональные скорости точек vk , т. е.
FkR ( vk ) = − µk vk ,
(3.37)
где µk – постоянные коэффициенты, характеризующие сопротивление.
Повторяя рассуждения, приведенные в параграфе 2.2, получим выражение для
обобщенной силы сопротивления Q Rj , соответствующей обобщенной координате q j
:
n
n
n
dv
∂rk
R ∂rk
=
−
=
−
µ
v
QqRj =
F
∑ k ∂q
∑ k k ∂q
∑ µk vk dq k =
j
j
j
=
k 1=
k 1 =
k 1
.
d  n µ k vk v k 
d  n µk vk2 
dФ
= − =
,
∑
 −
∑
 −=



dq
2
dq
2
d
q
j  k 1=
j k 1
j
=


(3.38)
где Ф – диссипативная функция (2.18).
Следовательно,
обобщенная
сила
сопротивления
равна
взятой
с
противоположным знаком частной производной от диссипативной функции по
соответствующей обобщенной координате.
Учитывая,
что
выражение
для
диссипативной
функции
аналогично
выражению для кинетической энергии, из (3.7) получим:
89
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
1
2
Φ = ( b11q12 + 2b12q1q2 + b22q22 ).
Постоянные
величины
называются
b11 , b12 , b22
(3.39)
обобщенными
коэффициентами сопротивления.
Так же как и кинетическая энергия, диссипативная функция является
однородной квадратичной формой обобщенных скоростей.
В
общем
случае
диссипативная
функция
является
неотрицательной
квадратичной формой [3, 4]. Если диссипативная функция является положительно
определенной квадратичной формой, и ее коэффициенты удовлетворяют условию (1.24):
b11 > 0, b 22 > 0,
b11 b12
> 0,
b12 b 22
(3.40)
то диссипация называется полной.
Пусть q1 и q 2 – главные координаты системы. Тогда выражения для
кинетической и потенциальной энергий примут вид:

=
T

П
=

1
( a11q12 + a22q22 ),
2
1
( с11q12 + с22 q22 ).
2
(3.41)
Диссипативная функция содержит слагаемое с произведением переменным и в
случае главных координат:
Φ=
(
)
1
b11q12 + 2b12q1q2 + b22q22 .
2
(3.42)
90
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Подставив (3.41) и (3.42) в уравнения Лагранжа (3.1), получим с учетом (3.38)
уравнения собственных колебаний системы с двумя степенями свободы при
наличии линейного сопротивления:
( a11q1 + b11q1 + c11q1 ) + b12q2 =
0,

b12q1 + ( a22q2 + b22q2 + c22 q2 ) = 0.
(3.43)
Общее решение системы (3.43) будем искать в виде
λt
=
q1 C=
C2 e λ t
1e , q2
(3.44)
Подставляя (3.44) в систему уравнений (3.43), после преобразований получим:
(
)
C1 a11λ 2 + b11λ + c11 + C2b12λ = 0,


2
0.
C1b12λ + C2 a22λ + b22λ + c22 =
(
)
(3.45)
Для того чтобы эта система однородных линейных уравнений относительно
неизвестных C1 , C2 имела ненулевое решение, определитель системы должен равняться нулю:
a11λ 2 + b11λ + c11
b12λ
b12λ
a22λ + b22λ + c22
2
= 0,
(3.46)
или
2
a11a22λ 4 + ( a11b22 + a22b11 )λ 3 + ( a11c22 + a22 c11 + b11b22 − b12
)λ 2 +
+( b11c22 + b22 c11 )λ + c11c22 =
0.
(3.47)
91
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Каждому из четырех корней
уравнения (3.47) λ1 ,λ2 ,λ3 ,λ4 соответствуют
определенные значения постоянных C1 и C2 .
Т.к. определитель системы (3.46) равен нулю, то уравнения (3.45) линейно
зависимы, поэтому их решение позволяет определить только отношения C1 и C2
каждого корня λi ( i = 1,2,3,4 ) :
a11λi2 + b11λi + c11
b12λi
C2i
1, 2, 3, 4 ) . (3.48)
=
−
=
−
(i =
2
i
b
λ
C1
a22λi + b22λi + c22
12 i
Возможны три варианта корней характеристического уравнения (3.47) и,
соответственно, общего решения системы дифференциальных уравнений (3.43).
1 Корни комплексные.
В случае малого сопротивления корни характеристического уравнения имеют
вид:
λ1,2 =
−n1 ± k1i; λ3,4 =
− n2 ± k 2 i ,
(3.49)
где n1 , n2 и k1 ,k2 – некоторые положительные числа.
Подставляя значения (3.49) в (3.44), после несложных преобразований
получим:
=
q1 e − n1t ( B1( 1 ) cos k1t + D1( 1 ) sin k1t ) + e − n2t ( B1( 2 ) cos k2t + D1( 2 ) sin k2=
t)
A1( 1 )e − n1t sin( k1t + α1( 1 ) ) + A1( 2 )e − n2t sin( k2t + α1( 2 ) );
=
q2 e
− n1t
( B2( 1 ) cos k1t + D2( 1 ) sin k1t ) + e − n2t ( B2( 2 ) cos k2t + D2( 2 )
A2( 1 )e − n1t sin( k1t + α 2( 1 ) ) + A2( 2 )e − n2t sin( k2t + α 2( 2 ) ),
sin k2=
t)
(3.50)
где B,D, A,α – постоянные величины.
Как следует из (3.50), каждая главная координата является суммой двух
затухающих колебаний. Четыре из восьми постоянных A1 , A2 ,α1 ,α 2 определяются
из начальных условий; остальные – через постоянные C1 ,C2 , определяемые из (3.48).
92
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2 Корни действительные.
В случае большого сопротивления корни характеристического уравнения
(3.47) можно представить в виде:
λ1,2 =
−n1 ± k1 ; λ3,4 =
− n2 ± k 2 ,
причем n1 > k1 ; n2 > k2 . Главные координаты (3.44) в этом случае примут вид:
− n1t
q=
( C1( 1 )ek1t + C1( 2 )e − k1t ) + e − n2t ( C1( 3 )ek2t + C1( 4 )e − k2t );
1 e
− n1t
q=
( C2( 1 )ek1t + C2( 2 )e − k1t ) + e − n2t ( C2( 3 )ek2t + C2( 4 )e − k2t ).
2 e
(3.51)
Как следует из (3.51) система будет совершать непериодическое, затухающее
движение.
3 Два корня действительные, два – комплексные.
λ1,2 =
−n1 ± k1 ; λ3,4 =
−n2 ± k2i.
Тогда главные координаты(3.44) примут вид:
− n1t
q=
( C1( 1 )ek1t + C1( 2 )e − k1t ) + A1( 2 )e − n2t sin( k2t + α1( 2 ) );
1 e
− n1t
q=
( C2( 1 )ek1t + C2( 2 )e − k1t ) + A2( 2 )e − n2t sin( k2t + α 2( 2 ) ).
2 e
В этом случае система будет совершать сложное движение, представляющие
собой наложение затухающих колебаний на апериодическое, затухающее движение.
Т.к., произвольные обобщенные координаты линейно выражаются через главные в
соответствии с (3.27), то в этом случае каждая из обобщенных координат является
линейной комбинацией трех рассмотренных случаев.
93
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
3.4 Вынужденные колебания системы без учета сопротивления
Рассмотрим вынужденные колебания системы с двумя степенями свободы без
учета сил сопротивления под действием гармонических возмущающих сил.
Получим уравнение движения системы в главных обобщенных координатах. В этом
случае кинетическая и потенциальная энергии определяются выражениями (3.41).
Пусть обобщенные возмущающие силы изменяются по синусоидальному закону
одинаковой частоты и фазы:
Q1B = H 1 sin( pt + δ ); Q2B = H 2 sin( pt + δ ) ,
где H 1 и H 2 – амплитуды обобщенных возмущающих сил;
p – частота возмущающей силы;
δ – начальная фаза возмущающей силы.
Из уравнений Лагранжа (3.2) получим следующую систему уравнений
вынужденных колебаний системы:
a11q1 + c11=
q1 H 1 sin( pt + δ ),

a22q2 + c22 q2 = H 2 sin( pt + δ ),
или
q1 + k12 q=
h1 sin( pt + δ ),
1

2
q2 h2 sin( pt + δ ).
q2 + k2 =
(3.52)
где введены следующие обозначения:
=
k12
c11
c22
H1
H2
=
; k22 =
; h1 =
; h2
.
a11
a22
a11
a22
(3.53)
94
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Общее решение системы дифференциальных уравнений (3.52) складывается
из общего решения системы соответствующих однородных уравнений (3.34) и
частного решения системы исходных неоднородных уравнений. Общее решение
системы (3.34) получено в пункте 3.2.3. Остановимся на определении частного
решения этой системы уравнений, характеризующего вынужденные колебания
системы.
Будем искать частное решение в виде:
q1B = A1B sin( pt + δ ); q2B = A2B sin( pt + δ ) ,
(3.54)
где A1B и A2B – постоянные коэффициенты, подлежащие определению.
Продифференцируем (3.54) дважды по времени:
q1B =
− p 2 A1B sin( pt + δ ); q2B =
− p 2 A2B sin( pt + δ ) .
(3.55)
Подставляя (3.54) и (3.55) в уравнения (3.52), получим после преобразований
систему алгебраических уравнений для определения A1B и A2B :
( k12 − p 2 ) A1B =
h1 ;
 2
2
B
h2 ;
( k2 − p ) A2 =
откуда
h1
 B
;
A
=
1

2
2
−
k
p
1


h2
 A2B =
.

k22 − p 2

(
)
(
)
(3.56)
С учетом найденных значений (3.56) вынужденные колебания определяются
выражениями:
95
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
q1B
=
(
h1
k12
−p
2
)
sin( pt + δ );=
q2B
(
h2
k22
−p
2
)
sin( pt + δ ) .
Если частота возмущающей силы совпадает с одной из частот собственных
колебаний, например p = k1 (резонанс на данной частоте), частное решение для
вынужденных колебаний ищем в виде
q1B = A1B t cos( pt + δ ); q2B = A2B sin( pt + δ ) .
(3.57)
где A1B и A2B – постоянные коэффициенты, подлежащие определению.
Подставляя (3.57) в систему уравнений (3.52), после преобразований получим:
h
h
A1B =
− 1 ; A2B =
− 2 2 2.
2k1
k2 − p
(3.58)
С учетом (3.58) вынужденные колебания в главных координатах определяются
выражениями:
π
h1t
h1t
 B
=
−
+
=
−
+
−
δ
δ
q
cos(
pt
)
sin(
pt
);
 1
2
2k1
2k1


h2
 q2B
=
sin( pt + δ ).

k2 − p 2
Произвольные обобщенные координаты связаны с главными соотношениями
(3.27), следовательно, при резонансе обе произвольные обобщенные координаты
будут неограниченно возрастать во времени.
В случае совпадения обеих собственных частот с частотой возмущающей
силы k=
p (резонанс по обеим координатам), обе главные координаты будут
1 k=
2
неограниченно увеличиваться во времени. Таким образом, для наступления
96
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
резонанса в системе с двумя степенями свободы достаточно совпадения частоты
возмущающей силы с одной из двух частот собственных колебаний.
3.5 Влияние линейного сопротивления на вынужденные колебания
Рассмотрим вынужденные колебания системы с линейным сопротивлением
при гармоническом возмущении. Уравнения движения системы в произвольных
координатах примут вид:
(
(
) (
) (
)
)
 



H 1 sin( pt + δ );
2
 a11q1 + b11q1 + c11q1 + a12 q2 + b12 q2 + c12 q=
(3.59)

H 2 sin( pt + δ ).
 a12q1 + b12q1 + c12 q1 + a22q2 + b22q2 + c22 q=
2
Общее решение системы уравнений (3.59) для каждой координаты является
суммой собственных q1C , q2C и вынужденных колебаний q1B , q2B :
q1 =
q1C + q1B ; q2 =
q2C + q2B .
В
зависимости
от
вида
корней
характеристического
(3.60)
уравнения
соответствующей системы однородных уравнений, получаемой из (3.59) путем
отбрасывания правых частей, общее решение (3.60) будет либо линейной
комбинацией затухающих колебаний и затухающего непериодического движения,
либо каждым из этих движений в отдельности.
Вынужденные колебания q1B , q2B являются частным решением системы
уравнений (3.59), которое следует искать в виде
 q1B B1 sin( pt + δ ) + D1 cos( pt=
=
+ δ ) A1B sin( pt + δ − ε 1 );
 B
=
+ δ ) A2B sin( pt + δ − ε 2 ).
 q2 B2 sin( pt + δ ) + D2 cos( pt=
(3.61)
97
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Постоянные A1B , A2B ,ε 1 ,ε 2 связаны с постоянными B1 ,B2 ,D1 ,D2 соотношениями:
 D1 
;
B
1 

A1B =
B12 + D12 ;
arctg  −
ε1 =
A2B =
B22 + D22 ;
 D 
arctg  − 2  .
ε2 =
 B2 
(3.62)
где A1B , A2B , ε 1 , ε 2 – соответственно амплитуды и сдвиги фаз.
Подставив
преобразований
(3.61)
получим
в
систему
систему
уравнений
четырех
(3.59),
после
алгебраических
несложных
уравнений
для
определения неизвестных коэффициентов B1 ,B2 ,D1 ,D2 :
 B1 ( c11 − a11 p 2 ) − D1b11 p + B2 ( c12 − a12 p 2 ) − D2b12 p =
H1 ;

2
2
0;
 B1b11 p + D1 ( c11 − a11 p ) + B2b12 p + D2 ( c12 − a12 p ) =

2
2
H2 ;
 B1 ( c12 − a12 p ) − D1b12 p + B2 ( c22 − a22 p ) − D2b22 p =

2
2
0.
 B1b12 p + D1 ( c12 − a12 p ) + B2b22 p + D2 ( c22 − a22 p ) =
(3.63)
Решая систему уравнений (3.61), определим постоянные B1 ,B2 ,D1 ,D2 , а затем
по формулам (3.62) постоянные A1B , A2B , ε 1 , ε 2 .
Как и в случае системы с одной степенью свободы, вынужденные колебания
являются незатухающими гармоническими колебаниями, амплитуда и фаза которых
не зависят от начальных условий, а частота совпадает с частотой возмущающей
силы. При резонансе амплитуды вынужденных колебаний остаются постоянными в
отличие от случая отсутствия сопротивления.
98
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
4 Теория удара
4.1 Основное уравнение элементарной теории удара
Явление, при котором за ничтожно малый промежуток времени скорости
точек тела изменяются на конечную величину, называется ударом.
Силы, действующие на тела, можно подразделить на силы, изменяющие
скорости точек в течение некоторого конечного промежутка времени, – конечные
силы (например, сила тяжести) и силы, изменяющие скорости точек тела в течение
весьма малого промежутка времени (порядка десятой и менее доли секунды), –
ударные силы.
Ударной называется такая сила, которая действуя в течение ничтожно малого
промежутка времени, достигает при этом настолько больших значений, так что ее
импульс есть величина конечная.
Импульс ударной силы определяется выражением:
τ
S = ∫ Fdt ,
(4.1)
0
где S – импульс ударной силы, H ⋅ c ;
F – ударная сила, Н;
τ – время удара, с.
Рассмотрим действие ударной силы на материальную точку за время удара.
Импульсы неударных сил за время удара будут малыми величинами и ими по
сравнению с ударными можно пренебречь. Тогда, на основании теоремы об
изменении количества движения материальной точки имеем:
mu − mv =
S,
(4.2)
99
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где m – масса точки, кг;
u – вектор скорость точки после удара, м/с;
v – вектор скорость точки до удара, м/с;
S – суммарный импульс ударных сил, H ⋅ c .
Из выражения (4.2) следует, что изменение количества движения
материальной точки за время удара равно импульсу действующих на точку
ударных сил.
Уравнение (4.2) называется основным уравнением теории удара.
Найдем из уравнения (4.2) скорость материальной точки в конце удара:
u= v +
S
.
m
(4.3)
На рисунке 23 показаны скорости точки до и после удара и импульс ударной
силы в соответствии с формулой (4.3). Здесь через 1 обозначена траектория
движения точки до удара, 2 – траектория движения точки после удара.
1
2
Рисунок 23
Таким образом, в элементарной теории удара вводятся следующие допущения:
1. Перемещениями точек тела за время удара можно пренебречь, считая тело
за время удара практически неподвижным.
2. Действием неударных сил за время удара можно пренебречь.
3. Изменения скоростей точек тела за время удара определяются основным
уравнением теории удара (4.2).
100
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
4.2 Общие теоремы теории удара
4.2.1 Теорема об изменении количества движения системы при ударе
Изменение количества движения системы за время удара равно сумме всех
внешних импульсов ударных сил, действующих на систему:
Q1 − Q0 =
∑ Ske ,
(4.4)
где Q1 – количество движения системы после удара;
Q0 – количество движения системы до удара;
∑ Ske
– сумма внешних импульсов ударных сил.
Если геометрическая сумма всех внешних импульсов ударных сил равна нулю
( ∑ Ske = 0 ), то количество движения системы за время удара не изменится, т.е.
Q1 = Q0 .
(4.5)
В проекциях на любую неподвижную координатную ось, например ось х,
уравнение (4.4) примет вид:
Q1x − Q0 x =
∑ Skxe .
(4.6)
Если сумма проекций всех внешних импульсов ударных сил на какую-либо
e
неподвижную ось, например ось х, равна нулю ( ∑ Skx
= 0 ), то проекция количества
движения системы на эту ось за время удара не изменится, т.е.
Q1x = Q0 x .
(4.7)
101
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
4.2.2 Теорема об изменении главного момента количества движения системы
при ударе
Изменение
главного
момента
количества
движения
системы
относительно какого-либо неподвижного центра за время удара равно сумме
моментов внешних импульсов ударных сил относительно того же центра:
KO ( u ) − KO ( v ) =
∑ M O ( Ske ) ,
(4.8)
где KO ( u ) – главный момент количества движения системы относительно центра О
после удара;
KO ( v ) – главный момент количества движения системы относительно центра О
до удара;
∑ M O ( Ske )
– сумма моментов внешних импульсов ударных сил относительно
центра О.
Если сумма моментов внешних импульсов ударных сил относительно какоголибо неподвижного центра О равна нулю ( ∑ M O ( Ske ) = 0 ), то главный момент
количества движения системы относительно этого же центра за время удара не
изменится, т.е.
KO ( u ) = KO ( v ) .
(4.9)
В проекциях на любую неподвижную ось х равенство (4.8) примет вид
Kx( u ) − Kx( v ) =
∑ M x ( Ske ) .
(4.10)
Если сумма моментов внешних импульсов ударных сил относительно какой102
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
либо неподвижной оси х равна нулю ( ∑ M x ( Ske ) ), то главный момент количества
движения системы относительно этой же оси х за время удара не изменится, т.е.
Kx( u ) = Kx( v ) .
(4.11)
4.3 Коэффициент восстановления при ударе о неподвижную поверхность
При прямом ударе тела о неподвижную поверхность величина k, равная
отношению модуля скорости тела в конце удара к модулю скорости в начале удара,
называется коэффициентом восстановления:
k=
u
.
v
(4.12)
Коэффициент восстановления изменяется в пределах 0 ≤ k ≤ 1 и определяется
опытным путем. В случае, если k = 0, удар называется абсолютно неупругим, если k
= 1 – абсолютно упругим. В практических расчетах коэффициент восстановления
обычно принимают зависящим только от материала соударяющихся тел, однако, как
показывают опытные данные, коэффициент восстановления зависит также от
соотношения масс соударяющихся тел и их формы.
4.4 Прямой центральный удар двух тел
Общая
нормаль
к
поверхности
соударяющихся
тел
в
точке
их
соприкосновения называется линией удара.
Удар называется центральным, если центры масс соударяющихся тел лежат
на линии удара (рисунок 24).
Центральный удар называется прямым, если скорости центров масс
103
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
соударяющихся тел в начале удара направлены по линии удара.
Обозначим массы соударяющихся тел через m1 и m2 , скорости их центров
масс до удара – v1 и v2 (рисунок 24,а), а в конце удара – u1 и u2 (рисунок 24,б).
Линия удара
х
х
а)
б)
Рисунок 24
Очевидно, что в рассматриваемом случае при ударе должны выполняться два
условия: v1 > v2 и u1 < u2 .
Ударные силы, действующие между телами, будут внутренними, т.е. сумма
внешних импульсов ударных сил равна нулю:
∑ Ske = 0 .
В результате уравнение
(4.6) примет вид:
Q1x = Q0 x ,
или
m1v1 + m2 v2 = m1u1 + m2u2 .
(4.13)
Коэффициент восстановления k в случае прямого центрального удара двух
тел определяется выражением:
k=
u2 − u1
.
v1 − v2
(4.14)
104
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Из (4.14) следует, что коэффициент восстановления при прямом центральном
ударе двух тел равен отношению модулей относительной скорости тел после
удара и до удара.
Для определения скоростей тел после удара u1 и u2 необходимо решить
систему уравнений (4.13) и (4.14).
В случае абсолютно неупругого удара ( k = 0 ):
u=
u=
1
2
m1v1 + m2 v2
.
m1 + m2
(4.15)
В случае абсолютно упругого удара ( k = 1 ):
2m2 ( v1 − v2 )

u
=
−
,
v
1
1

m
+
m

1
2

2m ( v − v2 )
u=
v2 + 1 1
.
2

m1 + m2
(4.16)
4.5 Теорема Карно
Потерянной скоростью называется разность скоростей тела (v − u ) до удара и
после удара.
Изменение кинетической энергии системы при упругом ударе равно
 1− k 
умноженной на коэффициент 
 кинетической энергии, которую имела бы
 1+ k 
система, если бы ее тела двигались с потерянными скоростями:
=
T0 − T1
1 − k  m1
m

( v1 − u1 )2 + 2 ( v2 − u2 )2  ,

1+ k  2
2

(4.17)
105
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
где T0 - кинетическая энергия системы в начале удара;
T1 – кинетическая энергия системы в конце удара.
Выражение (4.17) носит название теоремы Карно.
Частичная потеря кинетической энергии при неупругом ударе объясняется
остаточной деформацией и нагреванием соударяющихся тел.
Из (4.17) следует, что при абсолютно упругом ударе ( k = 1 ) потери
кинетической энергии не происходит ( T0 = T1 ).
4.6 Удар по твердому телу, вращающемуся вокруг неподвижной оси
Рассмотрим твердое тело, вращающееся вокруг неподвижной оси z (рисунок
25). Внешний ударный импульс, приложенный к телу в некоторой его точке,
вызывает появление в опорах тела дополнительных ударных реакций, которые
могут вызвать повреждение опор. Сформулируем условия, при которых опорные
устройства вращающегося тела не испытывают действия ударных сил.
Теорема об изменении главного момента количества движения системы (4.10)
в рассматриваемом случае примет вид
J z (ω1 − ω0 ) =
∑ M z ( Ske ) ,
(4.18)
где J z – момент инерции тела относительно оси вращения z , кг ⋅ м 2 ;
ω1 – угловая скорость тела после удара, рад / с ;
ω0 – угловая скорость тела до удара, рад / с .
Для того чтобы при ударе по телу, вращающемуся вокруг неподвижной оси z,
ударные импульсы в подшипниках S A и S B (рисунок 25) обращались в нуль,
должны быть удовлетворены следующие условия:
1 Внешний ударный импульс должен быть перпендикулярен к плоскости,
проходящей через ось вращения z и центр тяжести тела С.
106
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
2 Ось вращения тела z должна быть главной осью инерции в точке
пересечения ее с плоскостью, перпендикулярной к оси вращения и содержащей
J=
0.
точку приложения ударного импульса S , т.е. J=
xz
yz
3 Точка приложения D ударного импульса S должна отстоять от оси вращения
z на расстоянии
OD =
Jz
.
m ⋅ xC
(4.19)
Точка D, в которой приложен внешний ударный импульс S при отсутствии
ударных опорных реакций, называется центром удара.
А
z
xC
C
zC
O
y
D
х
В
Рисунок 25
107
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
5 Динамика плоского движения
5.1 Дифференциальные уравнения плоского движения тела
Рассмотрим
движение
твердого
тела,
масса
которого
распределена
симметрично относительно плоскости Оxy. Предположим, что все внешние силы
F1 ,F2 ,...,Fn , приложенные к телу, действуют в этой же плоскости, а начальные
скорости точек этого тела ей параллельны. При этих условиях тело будет совершать
плоскопараллельное движение, для изучения которого достаточно рассмотреть
движение плоской фигуры, полученной сечением тела плоскостью Оxy (рисунок 26).
Начало системы координат Cx1 y1 помещено в центре масс тела – точке С.
Рисунок 26
Разложим плоское движение тела на поступательное движение вместе с
подвижной системой отсчета Cx1 y1 и вращательное движение вокруг оси,
перпендикулярной плоскости Оxy и проходящей через центр масс С.
Как известно из кинематики уравнения плоского движения тела имеют вид:
108
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
 xC = f1( t ),

 yC = f 2 ( t ),
ϕ = f ( t ).

3
(5.1)
Основными задачами динамики плоского движения твердого тела являются
определение сил по заданным уравнениям движения в форме (5.1) (первая задача
динамики) и определение уравнений движения (5.1) по заданным силам (вторая
задача динамики).
Воспользуемся теоремой о движении центра масс:
MaC = ∑ Fke ,
(5.2)
где M – масса тела;
aC – вектор ускорения центра масс;
∑ Fke
– главный вектор внешних сил, приложенных к твердому телу.
Проецируя векторное равенство (5.2) на оси координат x и y, получим:
 Mx
 C = ∑ Fkxe ,

 C = ∑ Fkye .
 My
Для
составления
дифференциального
(5.3)
уравнения
вращательной
части
движения воспользуемся теоремой об изменении главного момента количества
движения системы:
dKC
= ∑ mC ( Fke ) ,
dt
(5.4)
где KC – главный момент количества движения твердого тела относительно оси,
109
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
перпендикулярной плоскости движения тела и проходящей через центр масс С;
∑ mC ( Fke ) – сумма моментов всех внешних сил относительно оси, проходящей
через центр масс С и перпендикулярной плоскости движения тела.
Как известно, момент количества движения твердого тела, вращающегося
вокруг неподвижной оси, определяется соотношением:
KC = I C ω ,
(5.5.)
где I C – момент инерции тела, относительно оси проходящей через центр масс С и
перпендикулярной плоскости движения тела;
ω – угловая скорость тела.
С учетом (5.5) формула (5.4) примет вид:
d( I Cω )
= ∑ mC ( Fke ) .
dt
Для твердого тела I C = const , тогда с учетом, что
(5.6)
dω
ϕ , получим
= ε= 
dt
дифференциальное уравнение вращательной части плоского движения:
I C
ϕ = ∑ mC ( Fke ) .
(5.7)
Уравнения (5.3) и (5.7) называются дифференциальными уравнениями
плоского движения твердого тела:
 Mx
 C = ∑ Fkxe ,

 C = ∑ Fkye ,
 My

ϕ = ∑ mC ( Fke ).
 I C
(5.8)
110
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
5.2 Задача о качении диска по наклонной плоскости
Рассмотрим
однородный
круглый
диск
с
горизонтальной
осью,
скатывающийся по наклонной плоскости с углом наклона α (рисунок 27).
Пренебрегая трением качения, исследуем движение диска для двух случаев:
при отсутствии скольжения по наклонной плоскости и при наличии скольжения.
5.2.1 Качение диска по наклонной поверхности при отсутствии скольжения
В случае отсутствия скольжения мгновенный центр скоростей диска
находится в точке соприкосновения диска с неподвижной плоскостью.
Тогда скорость центра масс диска
vC =ω ⋅ CP =ω ⋅ R .
(5.9)
Рисунок 27
111
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Располагая неподвижные координатные оси Оху как показано на рисунке 27,
получим:
yC =− R =const; ⇒ yC =0; yC =0.
С учетом (5.9) имеем:
a=
ϕR.
C x=
C 
(5.10)
Дифференциальные уравнения движения (5.8) в рассматриваемом случае
примут вид:
=
ϕ R Mg sin α − Fmp ;
 M

− N + Mg cos α ;
0 =
 I 
 Cϕ = Fmp R.
(5.11)
Решая эту систему уравнений, получим:

ϕ=
Fmp
=
MgR sin α
;
MR 2 + I C
MgR 2 sin α
aC =
;
MR 2 + I C
(5.12)
I C Mg sin α
µ Mg sin α ,
=
MR 2 + I C
(5.13)
где введено обозначение
µ=
IC
.
MR 2 + I C
(5.14)
112
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Для сплошного однородного круглого диска момент инерции относительно
центральной оси равен
MR 2
.
IC =
2
Тогда,
Fmp
=

ϕ=
2g sin α
;
3R
aC =
2g sin α
;
3
Mg sin α
1
.
; µ
=
3
3
Рассматриваемый случай качения диска без скольжения по наклонной
плоскости имеет место при условии, что сила трения при качении Fmp не превышает
*
предельную силу трения при скольжении Fmp
.
Первая определяется формулой (5.13), вторая определяется выражением:
*
F=
mp
fN
=
fMg cos α .
(5.15)
Тогда сформулированное выше условие примет вид:
µ Mg sin α ≤ fMg cos α ,
или
µtgα ≤ f .
(5.16)
Качение диска не будет сопровождаться скольжением, если угол α не больше
некоторой величины α0 , определяемой неравенством (5.16), т.е.
113
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
 f 
arctg   .
α ≤ α0 =
µ
Для сплошного однородного диска полученное выражение примет вид:
α ≤ α0 =
arctg ( 3 f ) .
5.2.2 Качение диска по наклонной поверхности при наличии скольжения
Как показано в пункте 5.2.1, в случае если α > α0 , качение диска по
наклонной плоскости будет сопровождаться скольжением. Дифференциальные
уравнения движения (5.8) примут вид:
*
 Mx
 С Mg sin α − Fmp
=
;

− N + Mg cos α ;
0 =

*
ϕ = Fmp
R.
 I C
С учетом (5.15) имеем:
 С Mg sin α − fMg cos α ;
=
 Mx

ϕ = fMgR cos α ,
 I C
Тогда дифференциальные уравнения движения тела примут вид:
=
xС g ( sin α − f cos α ) ;

fMgR cos α

ϕ=
.

I

C
(5.17)
Для сплошного однородного круглого диска система (5.17) примет вид:
114
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
=
xС g ( sin α − f cos α ) ;


2 fg cos α
ϕ=
.


R
При наличии скольжения диска по наклонной плоскости, скорость точки А
соприкосновения
диска
с
неподвижной
поверхностью
не
равна
нулю.
Следовательно, мгновенный центр скоростей находится в точке P (рисунок 28).
Рисунок 28
Как известно из кинематики скорости точек А и С определяются
соотношениями:
 v A= ω ⋅ AP;

 vC= ω ⋅ CP,
 v A= ω ⋅ AP;

=
ω ⋅ ( R − AP ).
 vC =
115
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Следовательно,
v A + vC =
ωR .
Продифференцировав полученное выражение по времени, выразим ускорение
точки А:
x=
ϕ R −xC .
A 
Откуда с учетом (5.17) имеем:
x A
=
 fMR 2
fMgR 2 cos α
fg cos α g cos α 
− g sin α + =
− tgα +
IC
I
 C

f ,

или

x A g cos α  f
=


 MR 2

+ 1  − tgα  .


 IC

(5.18)
Подставив (5.14) в выражение (5.18), получим ускорение точки А:
 f

x A g cos α  − tgα  ,
=
µ

или
x A g cos α ( tgα0 − tgα ) .
=
Если α0 > α , то выражение в скобках – величина положительная, и скорость
скольжения диска по наклонной плоскости будет возрастать с течением времени.
116
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
5.2.3 Качение диска по наклонной поверхности с учетом трения качения
При наличии трения качения дифференциальные уравнения движения (5.11)
примут вид:
 C Mg sin α − Fmp ;
=
 Mx

ϕ Fmp R − mmp ,
I C
=
N δ Mg cos α – момент трения качения;
где m=
mp δ=
δ – коэффициент трения качения.
После преобразований получим:
a=
C x=
C
Fmp =
MgR ( R sin α − δ cos α )
MR 2 + I C
Mg ( I C sin α + M δ R cos α )
MR 2 + I C
;
.
Условие отсутствия скольжения диска (5.16) примет вид:
µ sin α + δη cos α < f cos α ,
где введено обозначение
η=
MR
,
MR 2 + I C
или
tgα <
f − δη
µ
.
117
Copyright ОАО «ЦКБ «БИБКОМ» & ООО «Aгентство Kнига-Cервис»
Список использованных источников
1 Бутенин, Н.В. Курс теоретической механики: учебник для втузов. В 2 т. Т.2.
Динамика / Н.В. Бутенин, Я.Л. Лунц, Д.Р. Меркин. - 2-е изд., исправл. – М.: Наука,
2000. - 271 с.
2 Добронравов, В.В. Курс теоретической механики: учебник для вузов / В.В.
Добронравов, Н.Н. Никитин, А.Л. Дворников. – 3-е изд. перераб. – М.: Высшая
школа, 1974. – 528 с.
3 Ильин, М.М. Теория колебаний: учебник для вузов / М.М. Ильин, К.С.
Колесников, Ю.С. Саратов. – 2-е изд., стереотип. – М.: МГТУ, 2003. – 272 с.
4 Курс теоретической механики / под ред. К.С. Колесникова. – 2-е изд.,
стереотип. – М.: МГТУ, 2002. – 736 с.
5 Тарг, С.М. Краткий курс теоретической механики: учебник для втузов / С.М.
Тарг. - 18-е изд., стер. – М.: Наука, 2008. – 416 с.
6 Яблонский, А. А. Курс теоретической механики: учебник / А. А. Яблонский,
В. М. Никифорова. - 11-е изд., стер. - CПб.: Лань, 2004. - 768 с.
7 Яблонский, А. А. Курс теории колебаний: учебное пособие для вузов / А. А.
Яблонский, С. С. Норейко. - 5-е изд., стереотип. - СПб. :БХВ-Петербург, 2007. - 336 с.
118
Документ
Категория
Без категории
Просмотров
43
Размер файла
1 095 Кб
Теги
динамика, твердого, 5847, прикладное, тела, задачи
1/--страниц
Пожаловаться на содержимое документа