- •Е.А.Коган
- •и естественно-научным дисциплинам
- •Москва 2010
- •Настоящее пособие посвящено двум важным и тесно связанным между собой разделам математики: обыкновенным дифференциальным уравнениям и вариационному исчислению.
- •Подробное изложение теории дифференциальных уравнений и вариационного исчисления содержится во многих учебниках и пособиях, некоторые из которых, рекомендуемые студентам втузов для самостоятельного изучения, приведены в списке литературы.
- •1.1. Основные понятия
- •Обыкновенное дифференциальное уравнение n-го порядка в самом общем виде записывается так:
- •Для дифференциальных уравнений первого порядка различают общее, частное и особое решения, а также общий, частный и особый интегралы.
- •Поэтому общее решение дифференциального уравнения можно определить как множество всех частных решений уравнения.
- •Особым решением дифференциального уравнения называется решение, которое не может быть получено из общего решения ни
- •при одном частном значении произвольной постоянной.
- •Часто при интегрировании уравнения первого порядка не удается найти его общее решение в явном виде, а получается конечное (не дифференциальное) соотношение вида
- •Знание изоклин позволяет во многих случаях даже для не интегрируемых явно дифференциальных уравнений получить графическое решение задачи Коши и выявить характер интегральных кривых.
- •Пример. Построить методом изоклин интегральную кривую уравнения
- •Таким образом, изоклинами данного уравнения являются концентрические окружности с центром в начале координат, причëм угловые коэффициенты касательных к искомым интегральным кривым равны радиусам этих окружностей С.
- •Очевидно, это уравнение с разделëнными переменными. Интегрируя его, получим
- •Следовательно, общий интеграл уравнения будет
- •Тогда
- •Тогда уравнение примет вид
- •Дифференциальное уравнение вида
- •называется приводящимся к однородному. В частности, к этому классу относится уравнение вида
- •В результате уравнение приводится к однородному
- •Возвращаясь к старой переменной, получим
- •Общее решение уравнения примет вид
- •После разделения переменных получим
- •Тогда уравнение (1.18) примет вид
- •Следовательно,
- •Подставляя (1.19) и (1.20) в (1.16), получим общее решение уравнения (1.12) в виде (1.15).
- •Уравнение вида
- •Следовательно, его общий интеграл, а значит, и общий интеграл уравнения (1.24) имеет вид
- •В результате общее решение уравнения будет иметь вид
- •После интегрирования получим
- •Следовательно, общее решение уравнения (2.18) будет
- •Далее ищем решение уравнения (2.17) в форме, аналогичной по структуре выражению (2.19), но произвольную постоянную в (2.19) заменяем неизвестной функцией
- •Подставляя (2.20) в (2.17), получим
- •В качестве практического примера применения метода понижения порядка рассмотрим задачу об осесимметричном изгибе упругих круговых пластин.
- •Определителем Вронского или вронскианом называется функциональный определитель вида:
- •Поэтому теорема о структуре общего решения линейного однородного дифференциального уравнения n-го порядка может быть переформулирована так:
- •а общее решение запишется в виде
- •Если известно общее решение однородного линейного дифференциального уравнения, соответствующего заданному неоднородному, то его частное решение может быть найдено методом вариации произвольных постоянных (методом Лагранжа).
- •Рассмотрим его реализацию для линейных дифференциальных уравнений второго порядка:
- •Интегрируя дифференциальные уравнения первого порядка (2.66), находим
- •Пример. Решить уравнение
- •Определитель этой системы
- •Интегрируя уравнения (2.71), находим
- •Подставляя (2.72) в (2.70), получим общее решение уравнения в виде
- •Рассмотрим линейное однородное дифференциальное уравнение второго порядка
- •Пример. Найти собственные значения и собственные функции однородной краевой задачи для уравнения (2.90) при граничных условиях
- •Общее решение уравнения согласно (2.92) имеет вид
- •Вычисляем
- •Раскрывая определитель, получим уравнение относительно параметра
- •корни которого являются собственными значениями задачи:
- •Уравнением Эйлера называется линейное уравнение вида
- •Пример. Решить уравнение
- •Уравнение (2.103) есть уравнение Эйлера второго порядка. Применим замену независимой переменной и производных по формулам (2.101), (2.102). Тогда уравнение примет вид
- •Его общее решение
- •Эти решения находятся методом подбора (см. выше) и имеют вид
- •Поэтому общее решение уравнения (2.104) будет
- •Пример. Решить задачу Коши для уравнения
- •Общее решение уравнения имеет вид
- •Убедимся в этом на примере системы двух уравнений относительно двух неизвестных функций
- •Будем искать частное решение однородной системы в виде
- •Пример. Решить методом Эйлера систему
- •Характеристическое уравнение системы
- •Общее решение однородной системы (3.18) запишется в виде
- •Этот метод применим к решению систем неоднородных линейных уравнений n-го порядка. Ограничимся для простоты нормальной системой двух линейных уравнений с постоянными коэффициентами
- •Пусть общее решение однородной системы уравнений известно:
- •Пример. Решить систему
- •Принимаем частное решение системы (3.31) в виде
- •Интегрируя эти уравнения, получим
- •Общее решение системы запишется в виде:
- •Пусть требуется решить краевую задачу для уравнения
- •Будем искать решение задачи в виде
- •Задача сводится к решению дифференциального уравнения
- •Выбранные функции удовлетворяют всем перечисленным выше требованиям. Они линейно независимы, непрерывно дифференцируемы и удовлетворяют граничным условиям (4.12).
- •Решение методом Бубнова
- •Как видим, ошибка, полученная при решении методом Бубнова в двучленном приближении, составляет сотые доли процента.
- •Решение методом наименьших квадратов
- •Решение методом коллокаций
- •В результате находим
|
|
|
|
113 |
|
|
|
|
|
|
|
y (x) |
|
|
y (x) |
|
y |
(x) |
|
||
|
11 |
|
|
12 |
|
|
1n |
|
|
|
|
y21(x) |
|
|
y22 (x) |
|
y2n (x) |
|
|||
Y1 |
|
|
Y 2 |
|
|
|
|
|
|
|
(x) = . |
, |
(x) = . |
|
|
, ,Y n (x) = . |
|
|
|||
|
. |
|
|
. |
|
|
. |
|
|
|
|
. |
|
|
. |
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yn1 (x) |
|
yn 2 |
(x) |
yn n |
(x) |
||||
|
|
|
|
|
|
|
|
|
|
|
образуют фундаментальную систему решений этой системы. Следовательно, общее решение однородной системы диф-
ференциальных уравнений запишется в виде
Y (x)=C1 Y1 (x)+C 2Y 2(x)+…+C n Y n (x)
или
y1 (x) = C1 y11 (x) + C2 y12 (x) + + Cn y1n (x),
y2(x) = C1 y21(x)+C2 y22(x)+ +Cn y2n (x),
…………………………………………….
yn (x) = C1 yn1 (x) + C2 yn 2 (x) + + Cn yn n (x),
где C1, C2,…,Cn |
- произвольные постоянные. |
|
|
||||||||||||
Пример. Решить методом Эйлера систему |
|
||||||||||||||
|
|
y′ |
= −4 y |
+ y |
|
, |
|
|
(3.18) |
||||||
|
|
|
1 |
|
1 |
|
|
|
2 |
|
|
|
|
||
Принимаем |
|
y2′ = −2y1 − y2. |
|
|
|
||||||||||
y =α |
|
e k x, |
y |
|
=α |
|
e k x |
|
|
||||||
|
1 |
2 |
2 |
. |
(3.19) |
||||||||||
Подставляя (3.19) |
1 |
|
|
|
|
|
|
|
|
|
|
|
|||
в (3.18), |
получим систему |
|
|
||||||||||||
|
(−4 − k)α1 |
+α2 = 0, |
|
|
(3.20) |
||||||||||
|
|
|
+ (−1 − k)α2 = 0. |
|
|||||||||||
|
− 2α1 |
|
|
||||||||||||
Характеристическое уравнение системы |
|
|
|
|
|
|
|||||||||
|
|
− 4 − k |
1 |
|
|
= 0. |
|
|
|||||||
|
|
|
|
|
|
||||||||||
|
|
− 2 |
−1 − k |
|
|
|
|
|
|
|
|||||
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(−4 −k)(−1−k)+2 = 0 |
|
или |
|
|
|
k 2 |
+ 5k + 6 = 0. |
|
Корни уравнения k1 = −2, k2 = −3 - действительны и различны.
|
|
|
|
|
|
|
|
|
|
114 |
|
|
|
|
|
|
|
|
|
|
|
||
Подставляя k1 |
= −2 в (3.20), получим |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
− 2α1 +α2 |
= 0, |
|
|
|
(3.21) |
|||||||||||
|
|
|
|
|
|
|
|
2α1 +α2 |
= 0. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
− |
|
|
|
|
||||||||||||
Система (3.21) совместная, |
|
|
но |
неопределённая. |
Полагая |
||||||||||||||||||
α1 =α11 =1, находим α2 |
=α21 |
= 2. |
|
Поэтому |
|
|
|
|
|||||||||||||||
|
|
|
y |
|
|
= e−2x , |
|
|
|
y |
21 |
= 2e−2x . |
|
|
|
(3.22) |
|||||||
|
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Подставляем теперь в систему (3.20) |
|
k2 = −3. Получим |
|
||||||||||||||||||||
|
|
|
|
|
|
|
−α |
1 |
+α |
2 |
|
= 0, |
|
|
|
(3.23) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
− |
2α1 + 2α2 = 0. |
|
|
|
|
|||||||||||
Решение этой системы принимаем в виде α1 =α12 |
=1, α2 |
=α22 =1. |
|||||||||||||||||||||
При этом |
|
|
y |
|
|
= e−3 x , |
|
|
|
|
y |
|
=e−3 x . |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
22 |
|
|
|
(3.24) |
||||||||||
|
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Общее решение однородной системы (3.18) запишется в виде |
|||||||||||||||||||||||
y |
|
(x) = C y |
|
(x)+C |
2 |
y |
|
(x) |
= С e−2x +C |
e−3 x , |
|
||||||||||||
1 |
1 |
11 |
|
|
|
|
|
12 |
|
|
|
|
|
1 |
2 |
|
|
|
|||||
y |
2 |
(x) = C y |
21 |
(x)+C |
2 |
y |
22 |
|
(x) = 2С e−2x +C |
e−3 x . |
(3.25) |
||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
2 |
|
|
3.3. Метод вариации произвольных постоянных
Этот метод применим к решению систем неоднородных линейных уравнений n-го порядка. Ограничимся для простоты нормальной системой двух линейных уравнений с постоянными коэффициентами
dy1 |
= a y |
+ a y |
2 |
+ f |
(x), |
|
|
11 1 |
12 |
1 |
|
||
dx |
|
|
|
|
(3.26) |
|
dy2 |
= a21y1 + a22 y2 + f2 (x). |
|||||
|
|
|||||
dx |
|
|
|
|
|
Пусть общее решение однородной системы уравнений известно:
y10 |
= C1 y11 |
+C2 y12 |
, |
(3.27) |
|
y20 = C1 y21 +C2 y22, |
|||||
|
где С1 , C2 - произвольные постоянные, а y11 , y21, y12 , y22 - част-
ные решения однородной системы, соответствующие различным корням характеристического уравнения.
115
В соответствии с методом вариации частное решение неоднородной системы отыскивается в форме, аналогичной по структуре общему решению однородной системы, но произвольные постоянные в ( 3.27) заменяются неизвестными функциями, то есть принимается
y1 = C1 (x) y11 + C2 (x) y12, |
|
y2 = C1 (x) y21 + C2 (x) y22. |
(3.28) |
Подстановка (3.28) в ( 3.26) приводит к следующей системе двух дифференциальных уравнений первого порядка относитель-
но производных от неизвестных функций |
C1 (x), C2 (x): |
|
|||||||||
|
C′(x) |
y |
+C′ |
(x) y |
= |
f |
(x), |
(3.29) |
|||
|
|
1 |
11 |
2 |
12 |
|
|
1 |
|
|
|
|
C1′(x) y21 + C2′(x) y22 = f2 (x). |
|
|||||||||
Разрешая систему (3.29) по формулам Крамера, получим |
|||||||||||
два дифференциальных уравнения первого порядка |
|
||||||||||
С1′(x) = |
f1 y22 |
− f2 y12 |
, С2′(x) = |
|
f2 y11 − f1 y21 |
. |
(3.30) |
||||
y11 y22 |
− y12 y21 |
|
|
||||||||
|
|
|
|
y11 y22 − y12 y21 |
|
Интегрируя эти уравнения, находим варьируемые функции С1 (x), C2 (x) и подставляем их в (3.28). Общее решение системы
(3.26) запишется в виде
y1 = y10 + y1 ,
y2 = y20 + y2.
Пример. Решить систему |
|
|
|
|
|
y′ = −4 y |
+ y |
|
+ x, |
(3.31) |
|
1 |
1 |
|
2 |
|
|
y′ = −2y1 − y2 + 3x. |
|
||||
Общее решение однородной системы, согласно |
(3.25), |
||||
имеет вид |
|
|
|
|
|
y10 = C1 y11 +C2 y12 = С1 e−2x +C2 e−3x ,
y |
20 |
= C y |
21 |
+C |
2 |
y |
22 |
= 2С e−2x +C |
2 |
e−3x . |
(3.32) |
||||||
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|||||
Принимаем частное решение системы (3.31) в виде |
|
||||||||||||||||
|
|
y |
= C (x)e−2x + C |
2 |
(x)e−3 x , |
|
|
|
|||||||||
|
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
2 |
= 2C (x)e−2x + C |
2 |
(x)e−3 x . |
|
(3.33) |
|||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
116
Подставляя (3.33) в ( 3.31), получим после элементарных преобразований
C1′(x)e−2x + C2′(x)e−3 x = x,2C1′(x)e−2x + C2′(x)e−3 x = 3x.
Решение этой системы приводит к дифференциальным уравнениям первого порядка с разделяющимися переменными:
|
|
|
x |
e−3x |
|
|
|
|
|
|
e−2x |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
C1′(x) = |
|
|
3x e−3x |
|
|
= 2xe2x , |
C2′(x) = |
|
|
2e−2x |
3x |
|
= −xe3x . |
||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
e−2x |
e−3x |
|
|
|
e−2x |
e−3x |
|
|
||||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||
|
|
2e−2x e−3x |
|
|
|
|
|
2e−2x e−3x |
|
|
|
Интегрируя эти уравнения, получим
|
|
|
|
С1 (x) |
= 2∫xe |
2x |
dx |
= e |
2x |
|
|
|
− |
1 |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
x |
2 |
, |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
С2 (x) |
= −∫xe |
3 x |
dx |
= |
e3 x |
1 |
− x |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
3 |
|
3 |
. |
|
|
|
|
|
|
||||||||||||||||||||
Поэтому |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 = C1 (x)e |
−2x |
|
+ C2 (x)e |
−3 x |
= x − |
1 |
+ |
1 |
1 |
− x |
|
= |
2x |
|||||||||||||||||||||||
|
|
|
|
|
|
2 |
3 |
|
3 |
|
|
3 |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
y2 = 2C1 (x)e |
−2x |
+ C2 (x)e |
−3 x |
|
|
|
|
|
|
1 |
+ |
|
1 |
|
1 |
|
|
|
|
|
= |
|||||||||||||||
|
|
|
|
|
|
|
= 2 x − |
2 |
|
|
3 |
|
3 |
− x |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Общее решение системы запишется в виде: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
y |
= y |
|
+ y |
|
= C e−2x +C |
|
e−3 x |
+ 2x − |
7 |
, |
|
|||||||||||||||||||||||||
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
1 |
|
|
|
10 |
|
1 |
|
|
1 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
3 |
|
|
18 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
y |
2 |
= y |
20 |
|
+ y |
2 |
|
= 2C e−2x +C |
2 |
e−3 x + 5x − 8. |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
9 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−187 , 53x − 89.
3.4. Задача о свободных вертикальных колебаниях автомобиля
Эта задача подробно рассмотрена в [18]. Динамическая схема автомобиля показана на рис. I.15.
117
Рис. I.15
Будем предполагать, что вертикальные колебания автомобиля происходят параллельно его плоскости симметрии, деформациями кузова по сравнению с осадками опор можно пренебречь, поэтому на динамической схеме автомобиля горизонтальный стержень, моделирующий кузов автомобиля, можно считать абсолютно жестким.
Введëм следующие обозначения: m- масса кузова,
C1,C2 - приведëнные жесткости передней и задней подвесок автомобиля,
I p = mr 2 - момент инерции кузова относительно поперечной оси,
проходящей через центр тяжести (r - радиус инерции),
a,b - расстояния от передней и задней опор до центра тяжести, q1 - вертикальное перемещение центра тяжести (q1 > 0, если на-
правлено вниз),
q2 - угол поворота стержня относительно центра тяжести (q2 > 0,
если поворот происходит против хода часовой стрелки). Рассматриваемая система обладает двумя степенями свобо-
ды.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
118 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Так как осадки упругих опор равны |
|
q1 + aq2 |
|
(для передней |
||||||||||||||||||||||||||||||||||||||||||||||||||||
опоры) и |
q1 −bq2 (для задней опоры), то кинетическая и потен- |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
циальная энергия системы равны |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
T |
= |
1 |
|
|
dq |
|
|
2 |
+ |
1 |
I |
|
dq |
|
|
2 |
|
|
|
|
П = |
C (q + aq |
|
)2 |
+ |
C |
(q −bq |
|
)2 |
. |
||||||||||||||||||||||||||||
2 |
m |
1 |
|
|
|
2 |
|
|
|
|
|
2 |
|
, |
|
|
|
|
1 1 |
2 |
|
|
2 |
|
|
2 |
|
|
1 |
2 |
|
|
2 |
|
||||||||||||||||||||||||
|
|
|
|
|
dt |
|
|
|
|
|
|
|
p |
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂T |
|
|
|
|
|
|
|
dq |
|
|
|
|
|
d |
|
dq |
|
|
|
|
|
d |
2q |
|
|
|
|
|
|
|
||||||||||||
|
|
Так как |
|
|
|
|
|
|
|
|
|
|
|
|
|
= m |
|
|
|
1 |
, |
|
|
|
|
m |
|
1 |
= m |
|
|
|
1 |
, |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
dq1 |
|
|
|
|
|
|
|
|
|
dt2 |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
dt |
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
∂T |
|
|
|
|
|
|
|
dq |
|
|
|
|
|
|
|
d |
|
|
dq |
|
|
|
|
|
|
d |
2q |
|
|
|
|
|
∂T |
|
|
|
∂T |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
= |
I |
|
|
|
|
2 |
, |
|
|
|
|
|
|
|
I |
|
|
|
2 |
|
= I |
|
|
|
2 |
, |
|
|
|
|
|
= |
|
|
|
= |
0, |
|
|||||||||||
|
|
|
dq2 |
p |
|
|
|
|
|
|
|
dt |
p |
|
|
p dt2 |
|
|
∂q |
|
∂q |
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
∂ |
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
∂П |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂П |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
= C |
|
(q + aq |
|
) |
+C |
|
(q −bq |
|
), |
|
|
= C a(q + aq |
|
)−C |
b(q −bq |
|
), |
|||||||||||||||||||||||||||||||||||||||||
∂q |
|
2 |
2 |
2 |
|
|
∂q |
|
2 |
2 |
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
2 |
|
1 |
|
1 |
|
|
|
|
|
|
|
2 |
|
|
1 |
|
|
|
|
||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q1 |
|
то уравнения Лагранжа относительно обобщëнных координат |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
и q2 [22] |
|
|
|
|
|
d |
|
|
|
∂T |
|
|
|
|
|
∂T |
|
|
|
∂П |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(i =1,2) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
= − ∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
dq |
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
q |
i |
|
|
|
|
q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
примут вид:
m d 2q1 |
+(C |
+C |
)q +(C a −C |
b)q |
2 |
= |
||||||||
|
|
dt |
2 |
1 |
2 |
|
1 |
1 |
2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2q2 |
|
|
|
|
|
|
|
|
|
|
||
I |
p |
+(C a −C |
b)q |
+(C a2 |
+C |
|
b2 |
|||||||
|
dt |
2 |
|
1 |
2 |
1 |
1 |
|
2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
0,
(3.34)
)q2 = 0.
В результате приходим к системе двух обыкновенных линейных однородных дифференциальных уравнений второго порядка с постоянными коэффициентами.
Будем искать частные решения этой системы в виде
q1 = A1 sin(pt +ϕ), |
q2 = A2 sin(pt +ϕ). |
(3.35) |
Тогда подставляя (3.35) в (3.34), получим систему линейных однородных алгебраических уравнений относительно амплитуд A1 и A2 :
119
(C |
+C |
2 |
−mp2 )A +(C a −C |
b)A = 0, |
|
|
|||||||
|
1 |
|
|
|
1 |
1 |
2 |
|
|
2 |
|
(3.36) |
|
|
|
|
|
|
b)A +(C a2 |
+C |
b2 − I |
|
p2 )A |
= 0. |
|||
(C a −C |
p |
|
|||||||||||
|
1 |
|
|
2 |
1 |
1 |
2 |
|
|
2 |
|
|
|
Система (3.36) имеет тривиальное решение |
A1 = A2 = 0, ко- |
торое соответствует отсутствию колебаний и реализуется только при однородных начальных условиях
q1 = q2 = dqdt1 = dqdt2 = 0.
Это решение, очевидно, не представляет физического инте-
реса.
Ненулевые (нетривиальные) решения системы (3.36) существуют при условии равенства нулю определителя системы
|
|
|
C +C |
2 |
−mp2 |
|
|
C a −C |
b |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
= 0. |
(3.37) |
|
|
|
|
C a −C |
b |
C a2 |
+C |
b2 − I |
p |
p2 |
||||||||||||
|
|
|
1 |
|
|
|
2 |
|
1 |
|
|
2 |
|
|
|
|
|
|
|||
Раскрывая этот определитель, получим частотное уравне- |
|||||||||||||||||||||
ние: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(a +b)2 |
|
||
p |
4 |
C |
+C |
|
C a2 |
+C |
b2 |
|
2 |
+ |
C C |
|
(3.38) |
||||||||
|
− |
1 |
|
2 |
+ 1 |
2 |
|
p |
|
1 2 |
|
|
= 0. |
||||||||
|
|
|
|
m |
|
|
|
|
I p |
|
|
|
|
|
|
|
mI p |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Каждому из его корней соответствует определëнное соо т- ношение между амплитудами, определяющее собственную форму колебаний. Это соотношение можно найти из любого из уравнений системы (3.36), например, из первого
|
A |
= |
mp2 −(C |
+C |
|
) |
. |
(3.39) |
|||
|
2 |
1 |
|
|
2 |
|
|||||
A |
|
b |
|
||||||||
|
|
C a −C |
|
|
|
|
|||||
1 |
|
1 |
|
2 |
|
|
|
|
|
||
Рассмотрим различные частные случаи. |
|
||||||||||
1) Первый частный случай. |
|
|
|
|
|
|
|||||
Пусть |
|
r = |
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
(3.40) |
|
|
|
|
ab |
|
|
|
|
|
то есть радиус инерции равен среднему геометрическому величин a и b. Это соотношение реализуется, если расстояние a +b между осями колес заметно меньше общей длины автомобиля.
При выполнении условия (3.40) частотное уравнение (3,38) принимает вид
|
|
|
120 |
|
|
|
|
|
|
p4 − (C1a +C2b)(a +b) p2 |
+ C1C2(a +b)2 = 0. |
(3.41) |
|||||
|
mab |
|
|
m2ab |
|
|
||
Дискриминант этого биквадратного уравнения равен |
|
|
||||||
D = |
(C1a +C2b)2(a +b)2 |
− |
4C1C2(a +b)2 |
= |
(a +b)2(C1a −C2b)2 |
. |
||
|
|
|||||||
|
(mab)2 |
m2ab |
|
|
(mab)2 |
|
|
Следовательно, уравнение (3.41) будет иметь действительные и различные корни
p2 |
= |
C1 (a +b) |
, |
p2 |
= |
C2(a +b) |
. |
(3.42) |
|
||||||||
|
|
|||||||
1 |
|
mb |
2 |
|
ma |
|
||
|
|
|
|
|
Подставляя эти корни поочередно в соотношение (3.39), получим:
- для первой формы колебаний (при p2 = p2 ) |
|
|
|
|
||||||
|
|
|
|
C1 (a +b) |
|
1 |
|
|
|
|
|
A21 |
|
m |
−(C +C |
) |
1 . |
|
|||
|
|
mb |
|
|||||||
|
= |
|
|
1 2 |
|
= |
(3.43) |
|||
|
A |
|
C a −C |
b |
|
|||||
|
|
|
|
|
b |
|
||||
11 |
|
1 |
2 |
|
|
|
|
|
Здесь вторые индексы обозначают номер рассматриваемой собственной формы;
- для второй формы колебаний (при p2 |
= p22 ) |
|
|
|||||||
|
A22 |
|
m C1 (a +b) |
−(C +C |
) |
1 . |
|
|||
|
= |
ma |
|
|
1 2 |
|
= − |
(3.44) |
||
|
A |
C a −C |
b |
|
|
|||||
|
|
|
|
|
a |
|
||||
12 |
|
1 |
2 |
|
|
|
|
|
|
Эти формы показаны на рис. 1.16.
Из соотношений (3.43), (3.44) следует, что A11 = bA21 - вертикальное перемещение сосредоточенной массы при том, что пра-
вая пружина недеформирована, а A21- угол поворота стержня; при
колебаниях по второй форме недеформирована левая пружина, из (3.44) следует A12 = −aA22, поэтомуA22 теперь угловое перемеще-
ние стержня.
Таким образом, особенностью рассматриваемых форм колебаний является неподвижность одной подвески автомобиля при колебаниях другой подвески. Практически это означает, что подвески автомобиля являются независимыми, то есть колебания одной подвески не передаются другой. При этом свободные коле-
121
бания каждой из подвесок характеризуются различными частотами.
Рис. I.16
Свойства такой системы наглядно изображаются схемой, показанной на рис. I.17 (здесь масса m распределяется обратно пропорционально расстояниям от центра масс автомобиля, а моменты масс относительно центра тяжести одинаковы).
Рис. I.17 |
|
2) Второй частный случай: |
|
C1a = C2b. |
(3.45) |
Тогда уравнения (3.36) разделяются и принимают вид
A (C +C |
2 |
−mp2) = 0, |
(3.46) |
||||||
1 |
1 |
|
|
|
|
|
|
|
|
A (C a2 |
+C |
b2 |
− I |
p |
p2) = 0. |
(3.47) |
|||
2 |
1 |
|
|
2 |
|
|
|
|
Независимость уравнений (3.46) и (3.47) означает, что при соблюдении условия (3.45) могут происходить чисто вертикальные колебания без поворотов (подпрыгивание) (рис. I.18а) или чисто угловые колебания при неподвижном центре тяжести (га-
лопирование) (рис I.18б).
122
Рис. I.18
В первом случае (подпрыгивание без поворотов) имеем A1 ≠ 0, A2 = 0. Тогда из (3.46) следует C1 +C2 −mp2 = 0.
Поэтому квадрат первой частоты будет равен p2 = p12 = C1 m+C2 .
Во втором случае (угловые колебания при неподвижном центре тяжести) A1 = 0, A2 ≠ 0. Поэтому из уравнения (3.47) сле-
дует выражение для квадрата второй частоты p2 = p22 = C1a2 I+p C2b2 .
Таким образом, частоты подпрыгивания и галопирования различны. Формы таких колебаний показаны на рис. I.18.
3) Случай равенства обеих частот.
Частотное уравнение (3.38) будет иметь равные корни, если его дискриминант равен нулю. Это приводит к следующему условию:
|
|
|
|
|
2 |
|
|
2 |
2 |
|
|
2 |
|
|
|
|
|
|
C1a |
+C2b |
|
|
4C1C2(a +b) |
|
|
||||||
C +C |
2 |
+ |
|
|
|
|
− |
= 0. |
(3.48) |
||||||
|
|
|
|
|
|
||||||||||
|
1 |
|
|
|
r |
2 |
|
|
|
|
r |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Легко убедиться в том, что условие (3.48) выполняется, если одновременно выполняются условия (3.40) и (3.45).
Итак, если параметры динамической схемы автомобиля удовлетворяют условиям (3.40) и ( 3.45), то единственная собственная частота колебаний будет определяться формулой
p2 |
= p2 |
= 1 |
|
+C2 + C1a |
2 |
+C2b |
2 |
|
|
= C1 +C2 , |
|
C1 |
|
|
|
|
|||||||
1 |
2 |
2 |
|
m |
mr |
2 |
|
r2 |
|
m |
|
|
|
|
|
|
=ab |
||||||
|
|
|
|
|
|
|
|
|
C a=C |
b |
|
|
|
|
|
|
|
|
|
1 |
2 |
|
|
123 |
|
|
|
следовательно, |
|
|
|
|
p = |
|
C1 +C2 |
. |
|
|
|
|||
|
|
m |
Заметим, что схема автомобиля, у которого частоты p1 и p2
равны, считается более рациональной с точки зрения комфортности езды, так как в этом случае сильные колебания возникают лишь при некоторой единственной частоте толчков, тогда как в случае p1 ≠ p2 сильные колебания могут происходить при двух
разных значениях частот.