Курс лекций по теоретической механике Динамика (II часть) Бондаренко А.Н. Москва Электронный учебный курс написан на основе лекций, читавшихся автором для студентов, обучавшихся по специальностям СЖД, ПГС и СДМ в НИИЖТе и МИИТе ( гг.). Учебный материал соответствует календарным планам в объеме трех семестров. Для полной реализации анимационных эффектов при презентации необходимо использовать средство просмотра Power Point не ниже, чем встроенный в Microsoft Office операционной системы Windows-ХР Professional. Запуск презентации – F5, навигация – Enter, навигационные клавиши, щелчок мыши, кнопки. Завершение – Esc. Замечания и предложения можно послать по Московский государственный университет путей сообщения (МИИТ) Кафедра теоретической механики Научно-технический центр транспортных технологий

Содержание Лекция 9. Работа, мощность силы. Кинетическая энергия. Теоремы об изменении кинетической энергии для материальной точки и системы. Пример решения задач на использование теоремы об изменении кинетической энергии материальной точки. Лекция 9. Лекция 10. Пример решения задач на использование теоремы об изменении кинетической энергии системы. Потенциальное силовое поле. Силовая функция. Потенциальная энергия системы. Закон сохранения механической энергии. Лекция 10 Лекция 11. Динамика поступательного и вращательного движения твердого тела. Физический маятник.Лекция 11. Динамика плоского движения твердого тела. Принцип Даламбера для материальной точки и механической системы. Приведение сил инерции точек при поступательном и вращательном движениях. твердого тела. Лекция 12. Пример приведения сил инерции при вращательном движении тела. Приведение сил инерции точек при плоском движении твердого тела. Определение динамических реакций подшипников при вращении твердого тела. Балансировка. Лекция 12. Лекция 13. Аналитическая механика. Обобщенные координаты. Уравнения связей. Возможные перемещения. Идеальные связи. Принцип возможных перемещений. Примеры использования принципа возможных перемещений при определении реакций связей. Лекция 13 Лекция 14. Общее уравнение динамики. Пример решения задачи на применение общего уравнения динамики. Обобщенные силы. Лекция 14 Лекция 15. Уравнение Лагранжа II рода. Кинетический потенциал. Пример решения задачи на применение уравнения Лагранжа II рода. Вариационный принцип Гамильтона-Остроградского. Понятие об устойчивости состояния равновесия системы в потенциальном поле. Лекция 15 Лекция 16. Малые колебания систем с несколькими степенями свободы. Общая форма дифференциальных уравнений колебаний. Прямая форма. Обратная форма. Главные координаты. Свободные колебания с учетом сопротивления среды. Лекция 16 Лекция 17. Элементарная теория удара. Общие теоремы теории удара. Удар тела о неподвижную преграду. Случай косого удара. Гидравлический удар в трубах. Прямой центральный удар двух тел. Лекция 17 Рекомендуемая литература 1. Яблонский А.А. Курс теоретической механики. Ч.2. М.: Высшая школа г. 368 с. 2. Мещерский И.В. Сборник задач по теоретической механике. М.: Наука г. 416 с. 3. Сборник заданий для курсовых работ /Под ред. А.А. Яблонского. М.:Высшая школа г. 366 с. 4. Бондаренко А.Н. Теоретическая механика в примерах и задачах. Динамика (электронное пособие ), 2004 г.

Лекция 9 Работа, мощность силы. Кинетическая и потенциальная энергия – механическое движение в результате взаимодействия механических систем может переноситься с одной механической системы на другую: 1. без превращений в другую форму движения, т.е. в качестве того же механического движения, 2. с превращением в другую форму движения материи (потенциальную энергию, теплоту, электрическую энергию и т.д.) Каждый из этих случаев имеет свои измерители (меры) механического движения и механического взаимодействия, отстаиваемые в свое время Декартом и Лейбницем (см. таблицу): Мера механического движенияМера механического взаимодействия ДекартКоличество движенияИмпульс силы ЛейбницКинетическая энергияРабота силы Ф. Энгельс показал существование и равноправность обоих (векторных и скалярных) мер движения, каждой из которых соответствуют свои меры механического взаимодействия. Импульс силы является мерой действия силы при изменении механического движения. Работа является количественной мерой превращения механического движения в какую-либо другую форму движения материи. Работа силы, приложенной к материальной точке – Пусть точка приложения переменной по величине и направлению силы перемещается по некоторой произвольной траектории. На малом (элементарном) перемещении силу можно считать постоянной и элементарная работа силы равна проекции силы на направление перемещения (касательную к траектории движения), умноженной на элементарное перемещение : M T Знак элементарной работы определяется величиной угла и знаком cos : Поскольку часто более удобно работать с острыми углами, то в этом случае используют острый угол и знак присваивают по следующему простому правилу: если сила и перемещение совпадают по направлению, то присваивается знак +, если противоположны по направлению, то знак. Элементарная работа может быть записана в виде скалярного произведения: и в проекциях: Работа на конечном перемещении M M 1 получается суммированием или интегрированием: Частные случаи: 1. Сила постоянная по величине (F = const) и направлению ( =const): 2. Сила постоянная по величине (F = const) и параллельна перемещению ( =0): 3. Сила перпендикулярна перемещению: 1

Лекция 9 ( продолжение – 9.2 ) Можно доказать следующие теоремы и утверждения: Работа равнодействующей на некотором перемещении равна алгебраической сумме работ составляющих сил на том же перемещении: Работа постоянной сил по величине и направлению на составном перемещении равна алгебраической сумме работ этой силы на каждом из составляющих перемещений: Работа внутренних сил неизменяемой системы равна нулю: Работа силы тяжести не зависит от вида траектории и равна произведению силы тяжести на разность высот: Работа силы, приложенной к твердому телу, вращающемуся вокруг неподвижной оси. Запишем выражение для элементарной работы силы, приложенной к точке, и выразим элементарное перемещение через угол поворота тела: ω R d dsds T h h z -работа силы, приложенной к вращающемуся твердому телу, выражается через момент силы относительно оси. В частном случае постоянного значения момента силы относительно оси работа равна произведению момента силы на угол поворота: Работа силы, приложенной к вращающемуся твердому телу, для конечного угла поворота: Мощность – величина, характеризуемая количеством работы, произведенной в единицу времени: Мощность силы, приложенной к точке: Мощность силы, приложенной к вращающемуся твердому телу: Работа линейной силы упругости (реакции пружины) при перемещении из состояния равновесия: 2

Лекция 9 ( продолжение – 9.3 ) Кинетическая энергия – характеризует способность механического движения превращаться в эквивалентное количество другого движения: Кинетическая энергия материальной точки: Кинетическая энергия системы материальных точек: Кинетическая энергия твердого тела при поступательном движении: Кинетическая энергия твердого тела при вращательном движении: Кинетическая энергия твердого тела при плоском движении: Теорема об изменении кинетической энергии материальной точки – Изменение кинетической энергии точки равно работе сил, действующих на точку на том же перемещении: Запишем основной закон динамики точки: Выразим ускорение через скорость и умножим левую и правую части соотношения скалярно на дифференциал радиуса-вектора : Проинтегрируем полученное соотношение: После подстановки пределов получаем: Теорема об изменении кинетической энергии системы – Изменение кинетической энергии системы равно работе сил, действующих на систему на соответствующих перемещениях точек системы: Запишем теорему об изменении кинетической энергии для произвольной точки системы, при этом выделим работу внешних и внутренних сил, приложенных к данной точке: Просуммируем левые и правые части соотношений: В левой части получили разность кинетических энергий системы: Для неизменяемой системы: 3

Лекция 9 ( продолжение – 9.4 ) Пример решения задачи на применение теоремы об изменении кинетической энергии для материальной точки – Снаряд массы m выбрасывается пружинным устройством из канала под углом к горизонту. Длина нерастянутой пружины жесткостью c равна длине канала l 0. Перед выстрелом пружина сжимается на величину d. Определить скорость снаряда при вылете из канала, а также максимальную высоту полета. Дано:, c, d, m, l 0 Найти: v 1, H d 1. Выбираем объект - снаряд 2. Отбрасываем связи – ствол, пружину 3. Заменяем связи реакциями – N, R 4. Добавляем активные силы – G 5. Записываем теорему об изменении кинетической энергии для точки: Начальная скорость снаряда равна нулю: Работа сил, приложенных к объекту, равна: Работа нормальной реакции равна нулю (направление реакции перпендикулярно перемещению): Работа силы тяжести: Работа упругой реакции пружины (направление реакции совпадает с перемещением): Подставляем определенные величины в теорему: Отсюда величина скорости вылета снаряда: Определяем максимальную высоту полета (повторяем шаги 1-5): H Вертикальная скорость снаряда в наивысшей точке траектории равна нулю : Горизонтальная скорость снаряда постоянная (из закона сохранения проекции на ось x количества движения точки) и равна: Работа силы тяжести: Подставляем определенные величины в теорему: После некоторых сокращений и преобразований: Отсюда максимальная высота полета: Заметим, что предыдущее выражение можно более быстро получить, записывая теорему об изменении кинетической энергии только для вертикальной скорости движения точки, поскольку горизонтальные силы отсутствуют и горизонтальная скорость не изменяется.. 4

Пример решения задачи на применение теоремы об изменении кинетической энергии для системы – Массивный бумажный рулон радиуса R, приведенный в движение толчком, катится без проскальзывания по инерции вверх по наклонной шероховатой плоскости под углом к горизонту с некоторой начальной скоростью. Коэффициент трения качения f k. Определить начальную скорость рулона, необходимую для того, чтобы он мог перевалить через вершину высотой H от начального положения. Дано:, f k, H, R Найти: v 0 1. Выбираем объект - рулон 2. Отбрасываем связи – опорную плоскость 3. Заменяем связи реакциями – N, F тр, M к 4. Добавляем активные силы – G 5. Записываем теорему об изменении кинетической энергии для твердого тела: Кинетическая энергия на вершине равна нулю: Работа сил, приложенных к объекту, равна: Работа нормальной реакции равна нулю: Работа силы тяжести: Работа момента сопротивления качению: Подставляем определенные величины в теорему: Заметим, что выражение для начальной скорости не зависит от массы рулона. Масса рулона, как мера инертности, будет влиять на величину усилия, которое должно быть приложено к телу, чтобы сообщить ему указанную начальную скорость. Кинетическая энергия в начальный момент времени равна: Момент инерции массы сплошного цилиндра равен: Угловая скорость равна: Тогда кинетическая энергия в начальный момент времени: Работа силы трения скольжения равна нулю (приложена в МЦС): Момент сопротивления качению: Разность углов поворота рулона: После некоторых сокращений и преобразований получаем: Лекция 10 Потенциальное силовое поле Силовое поле – пространство, в каждой точке которого на материальную точку действуют силы, зависящие от координат точки. Стационарное силовое поле – действующие силы которого не зависят от времени, F = F(x, y,z) (поле силы тяжести, поле силы упругости). Нестационарное силовое поле - действующие силы которого зависят от времени, F = F(x, y,z, t) (электромагнитное поле). 5

Потенциальное силовое поле – в котором существует функция, в каждой точке пространства удовлетворяющая соотношениям: где U = U(x, y, z) – силовая функция. Лекция 10 ( продолжение – 10.2 ) Силовая функция определяется с точностью до постоянной: Основные свойства силовой функции: 1. Элементарная работа силы потенциального поля равна полному дифференциалу силовой функции: 2. Полная работа силы потенциального поля не зависит от траектории перемещения точки и равна разности значений силовой функции в конечном и начальном положениях: Следствие: Работа силы потенциального поля при перемещении точки по замкнутой траектории равна нулю: x y z Потенциальная энергия системы – функция, характеризующая запас энергии (потенциальной энергии) в данной точке потенциального силового поля. Потенциальная энергия равна работе сил потенциального поля, действующих на материальную точку, при ее перемещении из данного положения в начальное (нулевое). Для системы материальных точек потенциальная энергия равна сумме работ сил потенциального поля на всех перемещениях точек системы в начальное положение. Величина потенциальной энергии в начальном положении принимается равной нулю: П(x 0,y 0,z 0 ) = 0. В произвольной точке потенциальная энергия является функцией координат: П(x,y,z). Тогда по определению: – связь потенциальной энергии с силовой функцией. С учетом: П(x 0,y 0,z 0 ) = 0 соотношение связи можно записать как разность: Поскольку потенциальная энергия также определена с точностью до постоянной, то работа силы потенциального поля на перемещении из точки M 0 в точку M равна: Таким образом, изменение потенциальной энергии равно и обратно по знаку изменению силовой функции. Тогда: и 6

Лекция 10 ( продолжение 10.3 ) Примеры потенциальных силовых полей Поле силы тяжести. Сила тяжести, работа которой не зависит от траектории, является примером силы, имеющей потенциал – геометрическое место точек пространства, в которых потенциальная энергия постоянна. Проекции силы тяжести на координатные оси равны: Последнее выражение есть дифференциальное уравнение, которое легко решается разделением переменных и интегрированием левой и правой частей: x y z x y z Эквипотенциальные поверхности (П = const) представляют собой горизонтальные плоскости. Сила тяжести направлена перпендикулярно к этим плоскостям в сторону уменьшения значений потенциальной энергии. Работа силы тяжести на перемещении из точки M 1 в точку M 2 : Поле центральной силы притяжения. Силы тяжести могут считаться параллельными и постоянными по величине только в небольшой области пространства в поле тяготения Земли и эквипотенциальные поверхности могут считаться плоскими только в пределах этой области. В случае рассмотрения силы притяжения к центру величина силы прямо пропорциональна массе и обратно пропорциональна квадрату расстояния между материальной точкой и центром тяготения O: x y z x y z Проекции силы притяжения на координатные оси равны: Элементарная работа силы притяжения: Дифференциал потенциальной энергии: Полученное выражение есть дифференциальное уравнение, которое легко решается интегрированием левой и правой частей: Эквипотенциальные поверхности (П = const) поля центрального тяготения представляют собой сферические поверхности с центром в точке O. Сила притяжения направлена по нормали к этим поверхностям в сторону уменьшения значений потенциальной энергии. O Закон сохранения механической энергии – При движении механической системы в стационарном потенциальном поле полная механическая энергия системы остается постоянной. По теореме об изменении кинетической энергии системы: Отсюда: Сумму кинетической и потенциальной энергий называют полной механической энергией системы. 7

Динамика поступательного и вращательного движений твердого тела – рассмотренные теоремы динамики системы дают дифференциальные уравнения, описывающие эти два типа движения твердого тела. Дифференциальные уравнения поступательного движения твердого тела – из теоремы о движении центра масс системы: Дифференциальные уравнения вращательного движения твердого тела вокруг неподвижной оси – из теоремы об изменении момента количества движения системы: Физический маятник – твердое тело, имеющее неподвижную горизонтальную ось вращения, не проходящую через его центр тяжести, и находящийся под действием только силы тяжести. При отклонении физического маятника от положения равновесия возникает возвращающий момент от силы тяжести, наличие которого является условием колебательного движения (качания) относительно положения равновесия. x y z С O 1. Выбираем объект (маятник): 2. Отбрасываем связи (цилиндрические шарниры): 3. Заменяем связи реакциями (суммарные реакции двух шарниров): 4. Запишем дифференциальное уравнение вращения оси x : Представим уравнение в виде: - дифференциальное уравнение качаний физического маятника. Рассмотрим математический маятник длиной l : x y z O1O1 O Подставим момент инерции - дифференциальное уравнение качаний и представим уравнение в виде: математического маятника. Поскольку полученные уравнения отличаются лишь коэффициентами, то всегда можно поставить в соответствие физическому маятнику математический маятник, период качаний которого равен периоду данного физического маятника. Для этого достаточно приравнять коэффициенты: Отсюда можно определить приведенную длину физического маятника: Последнее неравенство легко доказывается: Точка O 1 физического маятника, находящаяся на расстоянии l по прямой OC называется центром качаний маятника. O1O1 В случае малых колебаний sin φ φ : Период колебаний: Используя формулу для периода колебаний можно определять опытным путем моменты инерции тел сложной формы (положение центра тяжести можно найти методом подвешивания). 8 Лекция 11

Динамика плоского движения твердого тела – Плоское движение может быть представлено как совокупность поступательного движения тела со скоростью центра масс и вращательного движения вокруг центра масс. Это представление было использовано ранее при вычислении кинетической энергии: Первые два дифференциальные уравнения, описывающие поступательное движение тела, найдем, используя теорему о движении центра масс: Здесь найдем дифференциальные уравнения движения твердого тела. В кинематике уравнения движения плоской фигуры были получены при использовании в качестве полюса любой произвольной точки. В динамике в качестве полюса выбирают центр масс: Третье дифференциальное уравнение, описывающее вращательное движение тела вокруг центра масс, найдем, используя теорему об изменении кинетического момента системы: Подобным образом с использованием различных теорем динамики могут быть получены соответствующие уравнения при рассмотрении динамики других более сложных видов движения твердого тела (сферическое, свободное, ударное взаимодействие). Принцип Даламбера (Германа, Эйлера) – общий метод, при помощи которого уравнениям динамики по форме придается вид уравнений статики. Благодаря простоте этот метод получил широкое применение во многих прикладных дисциплинах. Принцип Даламбера для материальной точки. Основное уравнение динамики точки: Перенесем произведение массы на ускорение в правую часть: Получившееся дополнительное слагаемое имеет размерность силы и принимается за силу инерции, направленную в сторону противоположную ускорению: С введением силы инерции уравнение динамики точки принимает вид уравнения равновесия: Таким образом, геометрическая сумма приложенных к точке сил и силы инерции этой точки равна нулю. Сила инерции условно добавляется к действующим на точку силам, образуя взаимно уравновешенную систему сил. Пример 1: Кабина лифта весом G поднимается тросом с ускорением a. Определить натяжение троса. 1. Выбираем объект (кабина лифта). 2. Отбрасываем связь (трос) и заменяем реакцией R. 3. Добавляем к действующим силам силу инерции: 4. Составляем уравнение равновесия: Определяем реакцию троса: Определяем натяжение троса: Сравните ход решения и результат с примером 1 в лекции 2 (стр.3) 9 Лекция 11 ( продолжение – 11.2 )

Лекция 11 ( продолжение – 11.3 ) Пример 2: Груз весом G подвешен на тросе длиной l и движется по круговой траектории в горизонтальной плоскости с некоторой скоростью. Угол отклонения троса от вертикали равен. Определить натяжение троса и скорость груза. 1. Выбираем объект (груз). 2. Отбрасываем связь (трос) и заменяем реакцией R. 3. Добавляем к действующим силам силу инерции: 4. Составляем уравнение равновесия: x y Из первого уравнения определяем реакцию троса: Определяем натяжение троса: Подставляем значение реакции троса и силы инерции во второе уравнение и определяем скорость груза: Принцип Даламбера для несвободной механической системы. Принцип Даламбера для k-той точки: Сложим все n уравнений: Здесь P k – равнодействующая задаваемых сил, приложенных к точке, R k – равнодействующая реакций связей, приложенных к точке, Ф k = - ma k – сила инерции точки. Здесь P * – главный вектор задаваемых сил, приложенных к точке, R * – главный вектор реакций связей, приложенных к точке, Ф * – главный вектор сил инерции точек системы. Таким образом, геометрическая сумма главных векторов задаваемых сил, реакций связи и сил инерции материальных точек равна нулю. Умножим уравнение, выражающее принцип Даламбера на радиус-вектор, проведенный из центра O к точке: Сложим все n уравнений: Здесь M O P – главный момент задаваемых сил относительно центра O, M O R – главный момент реакций связей относительно центра O, M O Ф – главный момент сил инерции точек системы относительно центра O. Таким образом, геометрическая сумма главных моментов задаваемых сил, реакций связи и сил инерции материальных точек относительно любого центра равна нулю. Приведение сил инерции точек твердого тела к простейшему виду - В динамике за центр приведения принимается обычно центр масс системы. В результате приведения сил инерции в общем случае получается главный вектор сил инерции и главный момент сил инерции относительно центра масс: 10 Сравните ход решения и результат с примером 3 в лекции 2 (стр.4).

Лекция 11 ( продолжение – 11.4 ) Главный вектор сил инерции твердого тела – не зависит от выбора центра приведения и для всех типов движения равен: Приведение сил инерции точек твердого тела при поступательном движении – В случае поступательного движения ускорения всех точек одинаковы и главный момент сил инерции относительно центра масс равен нулю: т.к. радиус-вектор центра масс обращается в нуль, если центр приведения совпадает с центром масс. Таким образом, силы инерции приводятся к равнодействующей силе, приложенной в центре масс, равной по модулю произведению массы тела на модуль ускорения его центра масс и направленной противоположно этому ускорению. Приведение сил инерции точек твердого тела при вращательном движении вокруг неподвижной оси – Рассмотрим тело, имеющее плоскость материальной симметрии, перпендикулярной оси вращения. В этом случае ось вращения является главной осью инерции тела в точке O. Ускорения всех точек, лежащих на одной прямой, параллельной оси вращения геометрически равны. Поэтому силы инерции симметричных точек относительно плоскости материальной симметрии также равны и их равнодействующая будет лежать в этой плоскости (в точке M к ). 11 z ω ε O MkMk C Главный вектор сил инерции равен: O C MkMk Ускорение центра масс: При выборе центра приведения в точке O главный вектор сил инерции должен быть приложен в этой точке параллельно вектору ускорения центра масс в противоположную сторону: В произвольной точке M к ускорение равно: Так как линии действия центробежных сил инерции проходят через центр вращения, то главный момент сил инерции вычисляется как сумма моментов только вращательных сил инерции: Главный момент сил инерции равен произведению углового ускорения на момент инерции тела относительно оси вращения и направлен в сторону противоположную угловому ускорению: Как и в статике, силу и пару можно заменить одной силой, равнодействующей, приложенной в новом центре приведения. Можно показать, что равнодействующая сил инерции будет приложена в центре качаний. В частном случае, если центр масс лежит на оси вращения, главный вектор сил инерции обращается в нуль и силы инерции приводятся к паре: Таким образом, силы инерции приводятся к главному вектору и главному моменту сил инерции.

Приведение сил инерции точек твердого тела при плоском движении – Рассмотрим тело, совершающее плоское движение, и имеющее плоскость материальной симметрии параллельную плоскости движения. Это движение может быть разложено на поступательное движение с центром масс тела C и вращательное вокруг подвижной оси z C, проходящей через центр масс тела перпендикулярно плоскости движения xCxC C yCyC zCzC ε В соответствии с этим силы инерции поступательного движения приводятся к главному вектору сил инерции, приложенному в центре масс, и главному моменту сил инерции (паре сил, лежащей в плоскости движения): Пример: Однородный стержень OA массы M длиной l, шарнирно подвешенный в точке O к вращающейся оси со угловой скоростью ω, находится в относительном равновесии под углом к оси вращения. Определить силы инерции и угол. ω O C x y z A Силы инерции в каждой точке стержня пропорциональны осестремительному ускорению, величина которого в свою очередь пропорциональна расстоянию точки до оси вращения (треугольная эпюра распределения). Элементарная сила инерции, приложенная к элементарной массе длины ds, расположенной на расстоянии s от точки О, равна: s ds dФdФ Главный вектор сил инерции находится интегрированием dФ по длине стержня: Этот же результат можно гораздо проще получить используя ускорение центра масс: Главный момент сил инерции нельзя найти по формуле M Ф 0 = I x x, т.к. стержень OA находится в относительном равновесии и x = x = 0. Однако силы инерции от вращения стержня относительно оси z создают момент сил инерции: Таким образом, силы инерции приводятся к главному вектору, приложенному в центре приведения O, и главному моменту относительно этого центра. Полученную систему силы и пары можно заменить одной силой, равнодействующей сил инерции, приложенной в точке, отстоящей от центра приведения по перпендикуляру к направлению силы на расстоянии: d Угол можно определить из уравнения относительного равновесия: Таким образом, равнодействующая сил инерции приложена в центре тяжести эпюры распределения сил инерции. Такому расположению соответствует центр качаний O 1. Лекция 12

Определение динамических реакций подшипников при вращении твердого тела вокруг неподвижной оси – Тело произвольной формы вращается вокруг неподвижной оси z с угловой скоростью ω и угловым ускорением ε. ω ε A C B x y z O 1. Освобождаем объект движения от связей и заменяем их реакциями: 2. Показываем внешние задаваемые (активные) силы: 3. Добавляем к системе сил силы инерции. В каждой точке тела прикладываем вращательную и центробежную силы инерции: MkMk K 4. Полученная система сил удовлетворяет уравнениям равновесия: 5. Силы инерции в общем случае приводятся к главному вектору и главному моменту сил инерции. Определим величину и направление компонент главного вектора сил инерции: O C x y Выбираем центр приведения в начале координат и прикладываем компоненты главного вектора сил инерции в этом центре: 6. Определим компоненты главного момента сил инерции относительно центра приведения как осевые моменты непосредственным суммированием по каждой из компонент сил инерции, приложенных в точке: K x y MkMk 6. Запишем уравнения равновесия полученной системы сил с учетом динамических добавок от возникающих сил инерции: Из первых пяти полученных уравнений можно определить реакции подшипников. Последнее уравнение является дифференциальным уравнением вращения. Динамические добавки зависят от положения центра масс тела и значений центробежных моментов инерции этого тела. zkzk xkxk ykyk 1313 Лекция 12 ( продолжение – 12.2 )

Лекция 12 ( продолжение – 12.3 ) Балансировка вращающегося тела – Добавление сил инерции к действующим силам есть результат введения подвижной системы координат, связанной с вращающимся телом. Полученные уравнения представляют собой относительные уравнения равновесия, записанные для этой подвижной системы координат: 1414 ω ε A C B x y z O MkMk K zkzk xkxk ykyk Из 4 и 5 уравнений системы определим эти составляющие: Реакции подшипников разложим на две составляющие, условно статическую и дополнительную (динамическую), например: Модуль полной дополнительной динамической реакции равен: Центробежные моменты инерции, вычисляемые относительно подвижных осей, не изменяются при вращении тела. Модуль полной дополнительной реакции зависит от угловой скорости и ускорения. В случае установившегося вращения ( = 0 ), величина полной дополнительной реакции пропорциональна квадрату угловой скорости. В современных машинах угловые скорости по величине могут быть значительными, так что дополнительные динамические реакции могут во много раз превышать статические реакции. Кроме того, направления этой динамической реакции и условно статической составляющей изменяются по отношению к неподвижной системе координат. Это вызывает знакопеременное нагружение опорных узлов, приводящее к их усталостному разрушению. Если центр масс находится на оси вращения, но главная центральная ось инерции тела не совпадает с осью вращения, то условно статические реакции при консервативных заданных силах не будут изменяться по величине и по направлению (см. уравнения 1,2), но дополнительные динамические реакции образуют пару сил: Плоскость действия пары и направления дополнительных реакций вращаются вместе с рассматриваемым телом со всеми вытекающими отсюда последствиями. При конструировании машин и механизмов необходимо исключить возникновение дополнительных динамических реакций. Для этого вращающиеся элементы должны иметь нулевые центробежные моменты инерции, т.е. главная центральная ось должна совпадать с осью вращения. При изготовлении таких элементов возможны незначительные отклонения в размерах, которые приводят к неуравновешенности возникающих при вращении сил инерции и возникновению дополнительных динамических реакций на опорные устройства. Для предотвращения этого выполняется балансировка - введение дополнительных масс, силы инерции которых уменьшают неуравновешенность сил инерции тела, или, напротив, высверливание части материала, что изменяет распределение сил инерции. т.к. по определению условно статических реакций.

15 Лекция 13 Аналитическая механика – устанавливает общие, единые методы изучения движения и равновесия любых самых сложных материальных систем средствами математического анализа. Для этого вводятся новые понятия и обобщаются старые. Связи – рассматриваются теперь как некоторые условия, налагаемые на систему, которые должны удовлетворяться в процессе движения системы. Они содержат соотношения (уравнения или неравенства) между координатами, компонентами скоростей и ускорений и, возможно, времени. Классификация связей: По интегрируемости: Голономные (геометрические) – выражаются конечными уравнениями относительно координат или интегрируемыми дифференциальными уравнениями относительно координат: Неголономные (кинематические) - выражаются неинтегрируемыми дифференциальными уравнениями относительно координат, т.е. уравнениями, содержащими не только координаты точек системы, но и их производные по времени: Неинтегрируемость состоит в том, что их нельзя привести к виду уравнений голономной связи. По зависимости от времени: Склерономные (стационарные) – не зависящие от времени: Например, уравнение траектории, полученное для некоторой точки шатуна кривошипно-шатунного механизма: рассматривается как уравнение cклерономной голономной связи: Реономные (нестационарные) – зависящие от времени. Например, кинематическое возбуждение колебаний. По освобождаемости: Неосвобождающие (удерживающие или двухсторонние) – описываются уравнением, исключающим возможность покидания точкой траектории или поверхности, описываемой уравнением. Этому соответствует, например, жесткая связь в виде шарнирного стержня. Освобождающие (неудерживающие или односторонние) – выражаются неравенством, регламентирующим связь лишь в одном направлении, например, гибкая нить или гладкая поверхность. Обобщенные координаты – независимые параметры, однозначно определяющее положение механической системы при ее движении. Обобщенность состоит в том, что они могут иметь различную природу (линейные или угловые перемещения относительно некоторого начального положения или какие-либо другие величины). Общее обозначение – q i (i = 1,…,n). Число степеней свободы – число независимых обобщенных координат, через которые можно выразить декартовые координаты всех точек системы. Например: A x y yAyA xAxA O Здесь положение любой точки стержня (например, А) однозначно определяется значением всего одной величины – угла, который является обобщенной координатой (q = ). Число степеней свободы равно n = 1. Уравнение связи для рассматриваемой точки A: Если на систему N точек в пространстве наложено m голономных связей, то декартовые координаты всегда могут быть выражены конечными соотношениями: Число обобщенных координат равно n = 3N – m.

Лекция 13 ( продолжение – 13.2 ) Возможные перемещения – бесконечно малые перемещения, допускаемые наложенными на систему связями. С точностью до бесконечно малых приращения радиуса-вектора лежат в касательной плоскости к поверхности связи и представляют собой возможные перемещения. В случае нестационарной голономной связи f(x,y,z,t) = 0 возможные перемещения рассматриваются для положения и формы поверхности связи, соответствующих данному моменту времени. Возможные перемещения не зависят от приложенных к системе сил. Действительные перемещения – бесконечно малые (элементарные) перемещения, действительно (фактически) происходящие за время dt, допускаемые наложенными на систему связями. Действительные перемещения зависят от сил, приложенных к системе, от вида связей (стационарных, нестационарных, голономных, неголономных) и начальных условий. Таким образом, возможные перемещения являются более общим понятием, чем действительные перемещения. T бsбs ds Поскольку вектор положения точки системы можно выразить через обобщенные координаты, то возможные перемещения выражаются через приращения обобщенных координат как полный дифференциал: или Вычисление возможных перемещений: Геометрический способ - в силу малости возможных перемещений при повороте твердого тела любая его точка может рассматриваться движущейся не по дуге, а по перпендикуляру к радиусу вращения в сторону угла поворота: бyAбyA бxAбxA A x y yAyA xAxA O x A O бxAбxA бyA=бsAбyA=бsA Для малых углов cos 1, sin, тогда: Например, для наклонного стержня: A y O x бsAбsA Аналитический способ – вычисляется вариация от координат: В отличие от геометрического способа знаки возможного приращения координат получаются автоматически. При использовании геометрического способа в дальнейших вычислениях, например, работы, необходимо учитывать направление полученного приращения (перемещения) Возможная работа силы – элементарная работа силы на том или ином возможном перемещении: В координатном виде: В естественном виде:

Лекция 13 ( продолжение – 13.3 ) 1717 Примеры использования принципа возможных перемещений для определения реакций связей: Пример 1. Определить реакцию балки в правой опоре: A B a l Балка неподвижна и не имеет ни возможных, ни действительных перемещений. Отбросим связь, реакция которой отыскивается, и заменим ее реакцией: Без правой опоры балка может поворачиваться под действием активных сил, реакцию R B причисляем к активным силам. Зададим малое возможное перемещение: Идеальные связи – связи, при которых сумма элементарных работ сил реакций связи на любом возможном перемещении равна нулю: Примеры идеальных связей: абсолютно гладкая поверхность (при скольжении), абсолютно твердая поверхность (при качении без скольжения). Любую неидеальную связь можно рассматривать как идеальную, если соответствующие реакции связи (совершающие работу на возможных перемещения) причислить к задаваемым (активным) силам. Принцип возможных перемещений – Для равновесия материальной системы, подчиненной голономным, стационарным, двухсторонним и идеальным связям, необходимо и достаточно, чтобы сумма элементарных работ всех активных сил на любом возможном перемещении из предполагаемого положения равновесия равнялось нулю: Доказательство необходимости: Система находится в равновесии и для каждой точки удовлетворяется уравнение равновесия: Умножим скалярно на вектор возможного перемещения точки и сложим: = 0 Доказательство достаточности: Дано: Предположим, что равновесия нет. Тогда каждая из точек под действием активных сил придет в движение, переместится за время dt на малое действительное перемещение dr. Рассматривая эти перемещения, как возможные, вычислим работу и просуммируем: = 0 Получили противоречие с исходным равенством. Значит предположение об отсутствии равновесия неверно. б бsPбsP бsBбsB Запишем сумму работ: Вычислим возможные перемещения: Пример 2. Определить опорный момент многопролетной составной балке в левой опоре: б бsPбsP бsBбsB Отбросим в жесткой заделке связь, препятствующую повороту балки, и заменим ее парой сил M A : MAMA A B a l l C D E b b бsDбsD Вычислим возможные перемещения: Запишем сумму работ: Заметим, что 1. для нахождения опорного момента M A из уравнений статики потребовалось бы решить как минимум три уравнения равновесия; 2. эпюра возможных перемещений пропорциональна линии влияния усилия; 3. если задать возможное перемещение для искомой реакции равным 1, например, б =1, то эпюра перемещений будет полностью тождественна линии влияния поскольку

1818 Лекция 14 Общее уравнение динамики – Принцип возможных перемещений, дающий общий метод решения задач статики, можно применить к решению задач динамики, а именно: 1. Применить принцип Даламбера, сводящий задачу динамики с задаче статики: 2. Применить принцип возможных перемещений, решающий эту статическую задачу: Просуммируем по всем точкам: = 0 – для идеальных связей Получим общее уравнение динамики: В любой момент времени сумма работ всех задаваемых сил и сил инерции несвободной механической системы с двухсторонними идеальными связями на любом возможном перемещении равна нулю. Более короткие записи общего уравнения динамики: или Или еще короче:где бA – возможная работа всех задаваемых сил и сил инерции на любом возможном перемещении. Пример. Центробежный регулятор вращается вокруг вертикальной оси с постоянной скоростью. При = 0 пружина не деформирована. Жесткость пружины c. Длина каждого из стержней l. Плечо подвески a. Вес каждого из шаров G, вес муфты G 1. Определить угловую скорость установившегося вращения для данного угла. a 1. Покажем заданные силы: 2. Добавим силы инерции: 3. Упругая связь (пружина), не являющаяся идеальной (совершает работу на возможных перемещениях), должна быть отброшена и заменена реакцией, которая включается в число заданных сил: Модуль реакции пружины пропорционален изменению длины (укорочению) пружины: 4. Определим проекции возможных перемещений (вариации координат) точек приложения сил: A B x y 5. Составим общее уравнение динамики: Подставим значения сил инерции и реакции пружины: Отсюда после некоторых сокращений и упрощений:

Лекция 14 ( продолжение – 14.2 ) Обобщенные силы – следующий шаг к обобщению, а именно, механического действия заданных сил на систему, после введения обобщенных координат (обобщения задания движения системы). Пусть механическая система имеет s степеней свободы, ее положение определяется s обобщенными координатами q 1, q 2,…, q s совокупность этих перемещений представляет одно из возможных перемещений системы. Элементарная работа всех заданных сил системы на этих перемещениях равна: Сообщим некоторой обобщенной координате q j бесконечно малое приращения, оставляя остальные обобщенные координаты неизменными, т.е. бq 1 =бq 2 = … = бq j-1 = 0, бq j 0, бq j+1 =…= бq s = 0. В результате все N точек системы получат какие-то бесконечно малые перемещения: Поставим в соответствие ко всем заданным силам системы некоторую одну (воображаемую) силу, которая совершает такую же работу на данном возможном (обобщенном) перемещении бq j, что и все силы системы: Отсюда величина этой силы определяется как: -обобщенная сила Q j, соответствующая обобщенной координате q j – скалярная величина, равная отношению элементарной работы заданных сил на всех перемещениях системы, вызванных элементарным приращением бq j 0 координаты q j, к величине этого приращения. 1. Размерность этой силы определяется размерностью обобщенной координаты. Например, если q j есть линейная обобщенная координата, то размерность обобщенной силы Q j соответствует силе (Н). Если q j есть угловая обобщенная координата, то размерность обобщенной силы Q j соответствует паре сил или моменту (Нм). 2. Число обобщенных сил равно числу обобщенных координат. Размерность каждой из обобщенных сил определяется размерностью соответствующей обобщенной координаты. Другие формулы для вычисления обобщенной силы: В векторной форме: Радиус-вектор k-той точки есть функция всех обобщенных координат: Вариация радиуса-вектора по обобщенным координатам при бq 1 =бq 2 = … = бq j-1 = 0, бq j 0, бq j+1 =…= бq s = 0 : Отсюда: В координатной форме: В случае потенциальных сил:

Лекция 14 ( продолжение – 14.3 ) Пример вычисления обобщенных сил – Для механической системы трех грузов с двумя неподвижными и одним подвижным блоками определить обобщенные силы Q j Число обобщенных сил равно числу обобщенных координат. Число обобщенных координат равно количеству степеней свободы, которое можно определить последовательным наложением связей: Ограничим горизонтальное перемещение груза 1, грузы 2 и 3 могут вертикально перемещаться. Ограничим дополнительно вертикальное перемещение, например, груза 3. Итак, n = 2. Выбираем обобщенные координаты q 1 = s 1 и q 2 = s 2 : s1s1 s2s2 2. Для определения Q 1 задаем произвольное малое перемещение бq 1 = бs 1 (бq 2 = бs 2 =0). бs1бs1 бs1бs1 бs 1 /2 Вычисляем возможную работу заданных сил: 3. Для определения Q 2 задаем произвольное малое перемещение бq 2 = бs 2 (бq 1 = бs 1 =0). Вычисляем возможную работу заданных сил: бs2бs2 бs2бs2 бs 2 /2 Груз 2 перемещаться не может (связи считаем двухсторонними). Q1Q1 Q2Q2 Уравнения равновесия в обобщенных силах – Согласно принципа возможных перемещений при равновесии системы: Зададим возможные перемещения точек системы, вызванные бесконечно малыми приращениями всех обобщенных координат: Вычислим возможную работу заданных сил: Перегруппируем суммы произведений: QjQj Приращения обобщенных координат произвольны и независимы друг от друга. Поэтому в полученном уравнении все коэффициенты при них (обобщенные силы) должны быть равны нулю: - условия равновесия сил в обобщенных силах. В рассмотренном выше примере, для равновесия системы необходимо, чтобы Q 1 и Q 2 равнялись нулю. Видно, что Q 1 0 и равновесия нет. Равновесие этой системы возможно лишь при наличии силы трения определенной величины между грузом 1 и опорной плоскостью. Тогда эта сила войдет в выражение для Q 1 : Теперь уравнения равновесия для данной системы определяют соотношения между силами и имеют вид:

2121 Лекция 15 Уравнение Лагранжа II рода – Уравнения представляют собой дифференциальные уравнения движения системы относительно обобщенных координат системы. Воспользуемся общим уравнением динамики: где бA – возможная работа всех задаваемых сил и сил инерции на любом возможном перемещении. 1. Зададим возможные перемещения точек системы, вызванные бесконечно малыми приращениями всех обобщенных координат: Вычислим возможную работу заданных сил и сил инерции: Перегруппируем суммы произведений: или QjQj QjФQjФ Приращения обобщенных координат произвольны и независимы друг от друга. Поэтому в полученном уравнении все коэффициенты при них (обобщенные силы) должны быть равны нулю: - уравнения движения системы, эквивалентные общему уравнению динамики. 2. В обобщенные силы инерции Q j Ф входят массы и ускорения точек системы. Попытаемся выразить эти силы через скорости точек и в конечном итоге через кинетическую энергию: Добавим к этому выражению два одинаковых слагаемых разного знака и следующего вида: Вычислим частную производную кинетической энергии системы по обобщенной координате: Вычислим производную по времени от частной производной кинетической энергии системы по обобщенной скорости: Производная по обобщенной скорости имеет аналогичное выражение: Таким образом: Подставим в уравнение движения: Отсюда: - уравнения Лагранжа II рода. Для консервативных (потенциальных) сил:

Лекция 15 ( продолжение – 15.2 ) 2 Кинетический потенциал – функция, определяемая выражением: L = T – П - функция Лагранжа, где и Кинетический потенциал L будет также функцией обобщенных координат, обобщенных скоростей и времени: Определим кинетическую энергию через кинетический потенциал как T = L + П и вычислим необходимые частные производные, участвующие в уравнении Лагранжа II рода: = 0, т.к не зависит от или a 2. Покажем заданные силы: 3. Упругую связь (пружину) заменяем реакцией и включаем ее в число заданных сил: 4. Определим проекции возможных перемещений (вариации координат) точек приложения сил: A B x y 5. Определим обобщенную силу Q: 6. Вычислим кинетическую энергию: Составляем уравнение Лагранжа II рода: Пример 1. Центробежный регулятор вращается вокруг вертикальной оси с постоянной скоростью. При = 0 пружина не деформирована. Жесткость пружины c. Длина каждого из стержней l. Плечо подвески a. Вес каждого из шаров G, вес муфты G 1. Определить угловую скорость установившегося вращения для данного угла. 1.Система имеет 2 степени свободы (поворот вокруг оси и изменение угла наклона стержней подвески). При установившемся вращении рассматриваем только изменение угла наклона и выбираем его в качестве обобщенной координаты q =. После подстановки R находим :

Лекция 15 ( продолжение – 15.3 ) Обобщенные силы Q 1, Q 2 были вычислены в примере: 3. Вычислим кинетическую энергию: Пример 2. Для механической системы трех грузов с двумя неподвижными и одним подвижным блоками определить ускорения грузов. 1.Система имеет 2 степени свободы (см. пример вычисления обобщенных сил Лекция 14, стр.19).Лекция 14, стр.19 Уравнения Лагранжа имеют следующий вид при выборе обобщенных координат q 1 = s 1 и q 2 = s 2 : s2s2 s1s1 4. Вычислим частные производные кинетической энергии: 5. Вычислим производные по времени: 6. Подставим полученные выражения и обобщенные силы в уравнения Лагранжа: Или: Неизвестные (независимые) ускорения: При равенстве масс M 1 =M 3 =M, как например, в задаче М [2]:

Лекция 15 ( продолжение – 15.4, дополнительный материал ) 24 Вариационный принцип Гамильтона-Остроградского – устанавливает, какому соотношению удовлетворяет действительное движение механической системы в некотором интервале времени в отличие от всех иных возможных движений (перемещений) – кривых сравнения. Кривая сравнения соответствует движению, допускаемому существующими связями, бесконечно близкому к действительному. Общее уравнение динамики имеет вид: или Первое слагаемое – работа задаваемых сил на возможном перемещении системы (бA). Попробуем представить второе слагаемое в виде совокупности членов, содержащих скорости и в конечном счете кинетическую энергию: Таким образом, общее уравнение динамики принимает вид: Потребуем, чтобы на границах интервала времени [t 1,t 2 ] действительная траектория совпадала с кривой сравнения: t=t1t=t1 t=t2t=t2 Данное требование эквивалентно отсутствию вариации перемещений в начале и в конце интервала [t 1,t 2 ], например, при рассмотрении свободных колебаний можно задаться формой решения: Такая вариация координаты обращается в нуль на концах интервала [0,T/2], где T – период колебаний (2π/k). Полученное общее уравнение динамики в дифференциальной форме справедливо в любой момент времени рассматриваемого интервала. Умножим его на dt и проинтегрируем по всему интервалу: = 0 Отсюда получаем вариационный принцип Гамильтона-Остроградского: -Только для действительного движения системы с голономными, двухсторонними и идеальными связями для данных условий интеграл по рассматриваемому интервалу времени суммы вариаций работы заданных сил и кинетической энергии равен нулю. В случае потенциальных (консервативных) сил: где, L = T – П - функция Лагранжа. Введем подобно импульсу действия силы интеграл вида - действие по Гамильтону. Тогда получаем вариационный принцип Гамильтона-Остроградского в другой форме: -Только для действительного движения консервативной системы с голономными, двухсторонними и идеальными связями для данных условий вариация интеграла S по рассматриваемому интервалу времени равна нулю или действие по Гамильтону имеет стационарное значение. (uv)=uv+uv или -uv=-(uv)+uv бAбA бT – d(…)/dt

Лекция 15 ( продолжение – 15.5, дополнительный материал ) 25 Вывод уравнения Лагранжа из вариационного принципа Гамильтона-Остроградского – установленный принцип Гамильтона- Остроградского позволяет получить уравнение Лагранжа II рода: Вариационный принцип Гамильтона-Остроградского: Подставим вариации работы и кинетической энергии в принцип Гамильтона-Остроградского: Возможная работа, выраженная через обобщенные силы: Вариация кинетической энергии в каждый момент времени: или Проинтегрируем члены с вариацией обобщенной скорости по частям: U dV U V V dU = 0 На границах интервала бq j = 0 Таким образом, принцип Гамильтона-Остроградского принимает вид: или Равенство нулю подынтегральной функции должно выполняться при любых значения вариации обобщенной координаты. Следовательно, для каждой из обобщенных координат коэффициенты при вариациях должны быть равны нулю: Понятие об устойчивости равновесия механической системы – Состояние покоя (равновесия) механической системы в консервативном (потенциальном) поле может быть устойчивым, неустойчивым и безразличным: Устойчивое – если система, выведенная из положения равновесия, возвращается в это положение и совершает колебания около него. Неустойчивое – если система, выведенная из положения равновесия при сколь угодно малом отклонении от него, не возвращается в положение равновесия и не совершает колебания около него. Безразличное - если система, выведенная из положения равновесия, занимает новое положение равновесия.

Лекция 15 ( продолжение – 15.6, дополнительный материал ) 26 Уравнения равновесия в обобщенных координатах (в обобщенных силах): По определению обобщенной силы для консервативной системы: Отсюда следует, что положению равновесия (покоя) консервативной системы соответствуют экстремальные значения потенциальной энергии системы. При этом по обращению в нуль частной производной потенциальной энергии нельзя судить об устойчивости состояния покоя (равновесия) в этих положения системы. Условие устойчивости состояния покоя устанавливается критерием Лагранжа-Дирихле: Те положения покоя консервативной системы, в которых потенциальная энергия достигает минимума, являются ее устойчивыми состояниями покоя: - условие минимума потенциальной энергии. Если вторая производная потенциальной энергии меньше нуля, то это соответствует случаю неустойчивого положения равновесия. Если вторая производная потенциальной энергии равна нулю, то она не может служить критерием минимума потенциальной энергии и для решения вопроса об устойчивости положения равновесия необходимо последовательно исследовать знаки производных более высокого порядка: Если первая по порядку, ненулевая производная потенциальной энергии, имеет четный порядок и больше нуля, то это соответствует минимуму потенциальной энергии (положение равновесия устойчивое). Если первая по порядку, ненулевая производная потенциальной энергии, имеет нечетный порядок, то потенциальная энергия не имеет экстремума (нет ни минимума, ни максимума), соответствует безразличному состоянию равновесия. Пример. Метроном представляет собой маятник с двумя грузами: A – неподвижный, весом G A, B – перемещаемый, весом G B. Определить условия устойчивого и неустойчивого положения равновесия. A B LALA LBLB x y Выберем в качестве обобщенной координаты угол отклонения стержня метронома от вертикали, φ : φ Потенциальная энергия системы грузов: Условие равновесия системы грузов: Исследуем устойчивость равновесия системы грузов при выполнении условий (a) и (b): При выполнении условия (a) (G A L A =G B L B ) вторая производная, как и все последующие, обращается в нуль. Это соответствует безразличному состоянию равновесия. Если G A L A >G B L B и φ =0 вторая производная оказывается больше нуля и это соответствует устойчивому состоянию равновесия. Если G A L A >G B L B и φ =180 о вторая производная оказывается меньше нуля и это соответствует неустойчивому состоянию равновесия.

27 Лекция 16 Малые колебания упругих систем с конечным числом степеней свободы – Рассмотрим упругую, например, двух опорную балку, распределенную массу которой далее будем считать сосредоточенной в нескольких точках, расположенных по длине балки, например, с некоторым равным шагом: Считаем, что точечные массы соединяются жесткими шарнирными стержнями, а упругое (изгибное) взаимодействие отдельных частей исходной балки представим в виде пружин определенной жесткости, обеспечивающих возвращение системы в положение равновесия при ее отклонении. m1m1 m2m2 m3m3 При этом перемещение одной из точек вызывает перемещение других соответственно некоторой гладкой кривой прогиба балки и значит определенную величину реакции каждой из пружин. Реакции упругих связей линейно связаны с перемещениями точек и могут быть вычислены как сумма реакций от единичных перемещений: Заметим, что такие связи не являются идеальными (совершают возможную работу) и должны быть заменены соответствующими реакциями: q 1 =1 q 2 =1 q 3 =1 В матричном виде: Здесь матрица C – симметричная матрица жесткости, элементы которой c ij – реакция по направлению i-той упругой связи от единичного смещения по направлению j-той связи. Аналогично можно представить обратную зависимость (перемещения точек от действия единичной силы): R 1 =1 q 11 q 21 q 31 R 2 =1 R 3 =1 q 22 q 23 q 33 q 12 q 13 В матричном виде: Здесь матрица D – симметричная матрица податливости, элементы которой d ij – перемещение по направлению i-той упругой связи от единичной силы, действующей по направлению j-той связи. Уравнения Лагранжа II рода для рассматриваемой системы: Обобщенная сила: Сумма полных работ всех сил: Приращение работы на вариациях перемещений: Возможная работа: Обобщенная сила: q 32

Лекция 16 ( продолжение – 16.2 ) 2828 Кинетическая энергия. Для механической системы материальных точек кинетическая энергия вычисляется как сумма: m1m1 m2m2 m3m3 Здесь скорость k-той точки: Тогда кинетическая энергия: = m ij В матричном виде: Здесь матрица M – симметричная матрица масс, элементы которой m ij – мера (сила) инерции, соответствующая направлению i-той упругой связи от единичного ускорения по направлению j-той связи: Вычислим необходимые производные кинетической энергии, участвующие в уравнении Лагранжа: Подставим найденные выражения в уравнение Лагранжа: или В матричном виде: - общая форма уравнений малых колебаний При получении уравнений в общей форме матрица масс M и матрица жесткости C являются полностью заполненными, например: Уравнение колебаний в матричном виде может быть представлено и так: Сравните с каноническим уравнением свободных колебаний материальной точки: Специальным выбором обобщенных координат можно добиться диагональности матрицы масс: q 1 =1 q 2 =1 q 3 =1 Тогда кинетическая энергия: Здесь матрица масс M становится диагональной: и система дифференциальных уравнений распадается относительно ускорений: - прямая форма уравнений малых колебаний Физический смысл прямой формы – система расчленяется на отдельные материальные точки, на каждую из которых действуют упругие силы безмассового каркаса.

Лекция 16 ( продолжение – 16.3 ) 29 Обратная форма – получается формально из обратных соотношений обобщенных координат: Здесь силу R j можно рассматривать как силу инерции, приложенной к связи: Тогда: или В матричном виде: Здесь матрица масс M сразу становится диагональной (массы расчленены), а матрица податливости D по-прежнему остается полностью заполненной (упругий безмассовый каркас соединяет точечные массы). Физический смысл – безмассовый каркас нагружен силами инерции со стороны материальных точек. Уравнение колебаний в обратной форме в матричном виде может быть представлено и так: Решение уравнений малых колебаний – получим формально с использованием матричной записи по аналогии с решением дифференциального уравнения свободных колебаний материальной точки: Будем искать решение в виде: - для любого момента времени Это уравнение представляет собой типичную задачу определения собственных значений (и собственных векторов), которое для произвольных значений амплитуд сводится к вычислению определителя следующего (векового) уравнения: Для уравнений малых колебаний в прямой форме матрица масс диагональная и вековое уравнение может быть записано, например, для i = 3 более подробно: Специальные методы и стандартные программы, составленные на их основе, решают данную задачу и в результате получается набор значений собственных частот ( ω i ) и соответствующих им собственных векторов ( φ i ). Последние представляют собой формы собственных колебаний для каждой из собственных частот: Формы собственных колебаний ортогональны между собой, что означает равенство нулю их скалярных произведений: φ 11 φ 12 φ 13 φ 23 φ 33 φ 22 φ 21 φ 31 φ 32 (ω1)(ω1) (ω2)(ω2) (ω3)(ω3) φ1φ1 φ2φ2 φ3φ3 Замечание. При использовании стандартных программ уравнения колебаний приводятся к виду: или Собственные числа и собственные вектора вычисляются для матрицы.

Лекция 16 ( продолжение – 16.4 ) 30 Главные координаты – С целью приведения матрицы масс и матрицы жесткости к диагональному виду можно использовать собственные формы колебаний в качестве новых обобщенных координат. Такие координаты называются главными в силу их ортогональности. В этом случае система дифференциальных уравнений распадается на отдельные совершенно независимые уравнения. Выразим старые координаты q i как линейную комбинацию ортогональных форм с новыми координатами (главными координатами): φ 11 φ 12 φ 13 φ 23 φ 33 φ 22 φ 21 φ 31 φ 32 (ω1)(ω1) (ω2)(ω2) (ω3)(ω3) φ1φ1 φ2φ2 φ3φ3 В матричном виде: Здесь матрица Ф – матрица собственных форм колебаний, элементы которой φ ij – элементы собственных векторов (ординаты эпюр прогибов): Обобщенные скорости выражаются через ортогональные формы и главные координаты аналогично. В матричном виде: Тогда кинетическая энергия системы в матричном виде: Здесь матрица масс – обобщенная матрица масс, которая становится диагональной в силу использованного ортогонального преобразования координат. Точно также потенциальная энергия системы: Здесь матрица жесткости – обобщенная матрица жесткости, которая так же становится диагональной по той же причине. Уравнения движения в главных координатах имеют вид: Подстановка общего решения приводит к равенству, откуда: Свободные колебания с учетом сопротивления среды – Сила вязкого сопротивления движению пропорциональна скорости: Обобщенная сила с учетом силы сопротивления: Уравнения движения в общей форме с учетом силы сопротивления: В матричном виде: Здесь матрица F – диагональная матрица демпфирования

Лекция 16 ( продолжение – 16.5, дополнительный материал ) 31 Примеры вычисления собственных частот для различных схем дискретизации массы – С целью исследования влияния схемы дискретизации массы рассмотрим для одной и той же балки длиной L = 3а, a = 1 м, погонный вес G 1 = 21 кгс/м (206 Н/м), момент инерции поперечного сечения I = 1840 см 4 = м 4, модуль упругости E = кгс/см 2, несколько схем: 1. С одной точечной массой – балка заменяется одной, расположенной посредине, сосредоточенной массой, равной массе балки. 3а3а m =M P=1 Воспользуемся обратной формой: Матрицы имеют первый порядок и представляют собой один всего лишь элемент: Элемент матрицы податливости определим как перемещение от единичной силы с помощью интеграла Мора, используя для его вычисления правило Верещагина: Произведение элементов матриц жесткости и податливости: Обратная величина: Собственная частота колебаний (рад/c): Точное решение (низшая частота колебаний): Ошибка аппроксимации составляет 29.8%. Такая большая ошибка связана с тем, что вся масса балки сосредотачивается в точке максимального прогиба. Вполне логично разбить балку на три части: средняя часть – половина массы балки прикладывается в средней точке, две части – две четверти массы прикладываются в опорных точках и тем самым не участвуют в колебаниях. Таким образом, остается одна точечная масса, равная половине массы балки: Ошибка аппроксимации составляет теперь всего 0.73%. И это очень хороший результат. 2. С двумя точечными массами – балка заменяется двумя, расположенными по третям длины балки, сосредоточенными массами, равными половине массы балки. 3а3а m =0.5M P=1 Матрицы имеют второй порядок: Матрица податливости: Произведение матриц жесткости и податливости: Обратная матрица: Собственные числа: Собственные частоты колебаний (рад/c): Ошибка аппроксимации по нижней частоте 18.4%. Вновь, если принять массу каждого из двух точечных масс, равной одной трети от массы балки при той же геометрии расположения, считая, что оставшаяся треть массы приходится на опоры, то ошибка по нижней частоте снижается до 0.101%. Это не означает, что такая высокая точность получается для высших частот. Для получения достаточной точности по этим частотам приходится увеличивать степень дискретизации (количество точечных масс).

32 Лекция 17 Элементарная теория удара. Удар – явление, при котором за ничтожно малый промежуток времени скорости точек изменяются на конечную величину. Ударные силы - силы взаимодействия при соударении тел (удар молота, столкновения экипажей). Время удара – очень малый промежуток времени, в течении которого происходит удар (контакт соударяющихся поверхностей). В силу этого ударные силы могут достигать очень больших значений, при которых возможно изменение скоростей точек на конечную величину. Соотношение между конечным изменением скорости и величиной ударной силы определяется теоремой об изменении количества движения: Здесь - импульс ударной силы. - основное уравнение удара. Импульс ударной силы является конечной величиной не смотря на то, что интегрирование должно выполняться практически на бесконечно малом интервале времени (времени удара). Точный закон изменения ударной силы в течении времени удара, как впрочем и само время удара, как правило, остаются неизвестным и интеграл заменяется произведением некоторого среднего значения силы на время удара: В силу того, что ударные силы много больше по величине других сил (неударных), последними пренебрегают. В силу малости времени удара, возникающие перемещения точек во время удара (v ср При рассмотрении механической системы во время удара из всех теорем динамики используется лишь теорема об изменении количества движения системы и для вращающейся системы ее аналог – теорема об изменении момента количества движения системы В проекции, например, на ось x В проекции, например, на ось z (относительно оси z) Теорема об изменении кинетической энергии использоваться практически не может, поскольку перемещениями во время удара пренебрегается и работа ударных сил не может быть вычислена. Удар шара о неподвижную поверхность – Рассматривается поступательное движение шара массой m со скоростью v перпендикулярно неподвижной массивной поверхности (преграде) – прямой удар. Например, шар падает с высоты h 0 и ударяется о горизонтальную поверхность со скоростью v. h0h0 2.Переход потенциальной энергии в кинетическую при восстановлении первоначальной формы тела за счет упругих сил. Из-за наличия остаточных (пластических) деформаций и нагрева тела кинетическая энергия полностью не восстанавливается и скорость u отделения шара от поверхности будет меньше, чем скорость до удара (u < v). h1h1 Различают две стадии (фазы) удара: 1. Переход кинетической энергии движения в потенциальную энергию деформации. При этом скорость падает до нуля, часть энергии расходуется на нагрев тела. Отношение модуля скорости шара в конце удара к модулю его скорости в начале удара – коэффициент восстановления при ударе : Коэффициент восстановления можно определить опытным путем:

Лекция 17 ( продолжение – 17.2 ) 3 Коэффициент восстановления может изменяться от 0 до 1. При k = 0 – абсолютно неупругий удар (шар не отскакивает от преграды), при k =1 – абсолютно упругий удар (нет потери энергии при деформации, нет нагрева). Реальные материалы всегда имеют такие потери энергии и коэффициент восстановления даже для достаточно упругих материалов лишь приближается в той или иной степени к единице. Кроме того коэффициент восстановления зависит от скорости, при которой происходит удар (k = k(v)). Поэтому сравнение значений коэффициентов восстановления должно выполняться при одной и той же скорости. Например, при скорости v = 3 м/с: стекло – k = 0.94; кость – k = 0.89;сталь – k = 0.56; дерево – k = Можно показать, что коэффициент восстановления определяет так же соотношение между импульсами ударной силы в двух фазах: Основное уравнение удара для первой фазы: для второй фазы: Отсюда, импульс второй фазы и суммарный импульс ударной силы в двух фазах зависят от коэффициента восстановления: Рассмотрим теперь поступательное движение шара массой m со скоростью v, составляющей некоторый угол (угол падения) к нормали неподвижной массивной поверхности (преграде) – косой удар. Запишем основное уравнение удара: Спроецируем на нормаль и касательную к поверхности: n β Коэффициент восстановления: Поскольку коэффициент восстановления k < 1, то угол отражения больше угла падения. Угол отражения равен углу падения только в случае абсолютно упругого удара (k = 1). Модуль скорости после удара: При очень больших углах падения, близких к прямому углу, скорость после удара приближается к скорости до удара (u v). Импульс ударной силы: При очень больших углах падения, близких к прямому углу, импульс ударной силы приближается к нулю (S 0). На этих свойствах, связанных с большими углами падения, основывается эффект запуска блинчиков метанием плоских камней (голышей) под острым углом к водной поверхности.

Лекция 17 ( продолжение – 17.3 ) 34 Гидравлический удар в трубах – при резком закрывании задвижки в магистральных трубопроводах может произойти ударное воздействие на элемент запирающего устройства и саму трубу (резко возрастает давление в системе). Рассмотрим трубу длиной L, по которой движется жидкость под некоторым давлением. L Пусть в некоторый момент поток резко перекрывается опусканием задвижки. Запишем теорему об изменении количества движения для объема жидкости, находящейся под давлением в трубе: Спроецируем уравнение на ось x, совпадающей с осью трубы: Здесь p – разность давлений, A – площадь поперечного сечения потока, t – время, за которое ударная волна пройдет расстояние L. Характер изменения ударной силы (разности давления) неизвестен. Ударный импульс заменим произведением средней разности давлений на площадь поперечного сечения потока и время импульса: Количество движения жидкости в начальный и конечный моменты: Подставим эти выражения в уравнение разности проекций количеств движения: Здесь с – скорость распространения ударной волны. Знак минус означает, что R B > R A и его можно далее опустить. С использованием коэффициента объемного веса получаем формулу Жуковского: Таким образом, повышение давления не зависит от длины трубы и пропорционально скорости течения жидкости. Скорость распространения ударной волны в жидкости: Поскольку при гидравлическом ударе происходит также деформация трубы, влияющая на изменение плотности жидкости Жуковским была получена формула для скорости распространения волны с учетом этого: Здесь E – модуль упругости воды кгс/см2. Здесь E м – модуль материала трубы, d – диаметр трубы, б – толщина стенки трубы. Следовательно, гидравлический удар зависит от скорости и будет сильнее в трубах малого диаметра и трубах, изготовленных из материалов с более высоким модулем упругости (например, модуль упругости стали почти в два раз больше модуля упругости меди).

Лекция 17 ( продолжение – 17.4, дополнительный материал ) 35 Прямой центральный удар двух тел – Рассмотрим соударение двух движущихся тел со скоростями v 1 и v 2 (v 1 > v 2 ) массами M 1 и M 2. 1.В первой фазе удара ударная сила взаимодействия возрастает от нуля до максимального значения (деформация нарастает до момента выравнивания скоростей). Проекция на горизонтальную ось теоремы об изменении количества движения для всей системы дает: M1M1 M2M2 Для определения величины механического взаимодействия (импульса) составим такое же уравнение для одного тела, например, 1: Заметим, что разность скоростей представляет собой относительную скорость (скорость сближения) и поскольку v 2 < v 1, то ударный импульс, приложенный к телу 1, будет направлен в сторону, противоположную движению этих тел. Аналогично можно определить импульс, приложенный к телу 2, но быстрее и проще воспользоваться законом действия и противодействия: 2.Во второй фазе удара ударная сила взаимодействия уменьшается от максимального значения до нуля (упругие деформации восстанавливают полностью и частично форму тел и потенциальная энергия деформации переходит в кинетическую до отделения тел друг от друга). Проекция на горизонтальную ось теоремы об изменении количества движения для одного из тел, например, 1, дает : С использованием коэффициента восстановления можно записать: Поделим это уравнение на уравнение в красной рамке: Подставим выражение для скорости u: Заметим, что разность скоростей (v 2 - v 1 ) опять представляет собой относительную скорость (скорость сближения) и поскольку v 2 < v 1, то скорость тела 1 уменьшается и это уменьшение пропорционально массе тела 2 и относительной скорости. Аналогично можно определить скорость тела 2: Здесь скорость тела 2 увеличивается и это увеличение пропорционально массе тела 1. Замечания: 1. В частном случае равенства масс (M 1 = M 2 ) и абсолютно упругого удара (k =1) скорость тела 1 после удара будет равна скорости тела 2 до удара и наоборот. Т.е. если тело 2, например, как при игре в биллиард, покоилось, то после удара телом 1 тело 2 получит скорость тела 1, а тело 1 остановится. 2. Проверить полученные соотношения можно подставив их в закон сохранения количества движения: 3. Отношения модулей относительных скоростей до и после удара определяют коэффициент восстановления (или наоборот). Для этого подставим вычтем выражение для скорости u 1 (в синей рамке) из аналогичного ему выражение для скорости u 2 : Добавили два одинаковых слагаемых с противоположными знаками для получения еще одной разности скоростей. Определим модуль ударного импульса, приложенного к каждому телу за весь период упругого удара (за две фазы): Используем выражение, полученное для импульса первой фазы:

36 Автор благодарит вас, уважаемые студенты, за то, что вы воспользовались этим материалом для подготовки к экзаменам по рассмотренным разделам теоретической механики. Если представленный материал поможет молодым преподавателям теоретической механики подготовиться к чтению лекций или послужит основой для разработки собственного курса лекций, то автор будет только рад. Успеха всем! Об авторе Список трудов