Примеры выполнения курсовой работы

Запишем уравнения Лагранжа второго рода

.

Кинетическая энергия рассматриваемой механической системы относительно неподвижной системы отсчета: Т = Т ₁+ Т ₂. Представим ее как функцию времени t, обобщенных координат q ₁ = j, q ₂ =S и обобщенных скоростей , а именно:

T = T ( q ₁ = j, q ₂ =S, , t).

Кинетическая энергия основания, совершающего вращательное движение относительно оси OZ,

Кинетическая энергия цилиндра, совершающего плоско-параллельное движение относительно неподвижной системы отсчета,

Формулы преобразования координат и поворотная матрица относительно оси OZ в соответствии с формулами (3.18) и (3.19) имеют следующий вид:

– для центра масс основания

– для центра масс цилиндра

Матрица скоростей:

После приведения подобных членов относительно обобщенных скоростей, получаем T= T ₁ +T ₂,

Уравнения Лагранжа второго рода для данной системы имеют следующий вид:

– для q ₁

– для q ₂

Вычислим производные от кинетической энергии системы:

– для q ₁ =φ ,

окончательно

где

Для q ₂ = S ,

.

В левой части уравнений Лагранжа, как правило, оставляют слагаемые со вторыми производными от обобщенных координат. Все остальные переносятся в правую часть. Таким образом, обозначив слагаемое в последнем выражении как - окончательно получаем уравнение Лагранжа второго рода

.

Обобщенные силы для механических систем с числом степеней свободы i = s = 2 , , отвечающие выбранным обобщенным координатам, целесообразно вычислять последовательно, учитывая, что обобщенные координаты, а значит, и их вариации не зависят друг от друга. Системе всегда можно сообщить такое виртуальное перемещение, при котором изменяется только одна обобщенная координата, а другие при этом не варьируются. В этом случае из (2.15) получаем:

. (6.6)

Для определения обобщенной силы дадим системе такое приращение, что 1) , и найдем виртуальную работу от всех заданных активных сил:

,

, ,

Сравнивая множитель (в квадратных скобках) в выражении виртуальной работы перед вариацией dj с формулой (6.6), получаем выражение для первой обобщенной силы

Определяя , будем полагать, что d S ¹ 0,а для угла jповорота цилиндра 2, будем считать dj =0 (j = const), т.е.

2)d S ¹ 0, dj =0 (j = const); .

Сравнивая множитель (в квадратных скобках) в выражении полученной виртуальной работы перед вариацией d S с формулой (6.6), получаем выражение для второй обобщенной силы

.

Составим дифференциальные уравнения в виде матриц

.

Для

Для

окончательно

1)

2) ,

где , если , , если ,

т.е. ;

;

.

Пример 2. Несвободная система с тремя степенями свободы. Дано. Рекомендуемые значения физических величин (рис. 6.5): m ₁ = 3 700 кг – масса тела; = 10 300 кг·м² – момент инерции относительно оси О ₁ z ₁; l ₁ = 3,8 м – расстояние от точки С до точки К; координаты центра масс тела О ₁в системе координат СX ₁ Y ₁ Z ₁: х ₁₀ = 3,8 м, y ₁₀ = 0,8 м, y ₂₀ = 1,0 м – расстояние от линии действия силы до оси СХ ₁; С ₁ = 2,0·10⁶ Н/м – жесткость пружины КМ; С ₂ = 2,0·10⁷ Н/м – жесткость пружины CD; С ₃ = 2,0·10⁷ Н/м – жесткость пружины СЕ; С _j = 2,0·10⁶ Нм/рад – жесткость спиральной пружины; j₀ = 6°, 30° – начальный угол возвышения основания.

Рис. 6.5

За обобщённые координаты приняты следующие параметры:

q ₁ = x_C, q ₂ = y_C, q ₃ = j. q ₁ = x_C – перемещение шарнира С по горизонтальной оси, отсчитанное от положения статического равновесия пружины CD; q ₂ = y_C – перемещение шарнира С вдоль вертикальной оси, отсчитанное от положения статического равновесия пружины CЕ; q ₃ = j – угол поворота тела, отсчитанный от горизонтали СХ ₁.

Кинетическая энергия рассматриваемого тела, совершающего плоскопараллельное движение относительно неподвижной системы отсчета и связанной с ней системой координат XOYZ:

.

Формулы преобразования координат от СX ₁ Y ₁ Z ₁ к XOYZ и поворотная матрица относительно оси СZ в соответствии с формулами (3.18) и (3.19) для центра масс тела 1 имеют следующий вид:

Матрица скоростей:

.

Квадрат абсолютной скорости центра масс тела в обобщённых координатах:

Кинетическую энергию тела представим как функцию времени t, обобщенных координат q ₁ = x_C, q ₂ = y_C, q ₃ = j и обобщенных скоростей :

T = T (q ₁= x_C, q ₂ = y_C, q ₃ = j, .

Уравнения Лагранжа второго рода для данной системы имеют вид:

– для ;

– для ;

– для .

Частные производные по обобщённым координатам:

;

Частные производные по обобщённым скоростям:

Полные производные по времени:

Запишем окончательно левые части уравнений Лагранжа второго рода, оставив только слагаемые со вторыми производными от обобщенных координат. Все остальные слагаемые перенесем в правую часть уравнения, обозначив их как :

–

где

–

где

–

, где .

Для нахождения обобщенных сил рассмотрим приложенные к системе силы. Определим первоначальные сжатия пружин КМ и СЕ в положении статического равновесия.

Уравнение равенства моментов при t =0 относительно оси CZ ₁

,

где - момент силы тяжести в момент времени t= 0; – момент силы упругости пружины КМ. Отсюда

- первоначальное сжатие пружины КМ.

Уравнение равенства сил по оси OY

,

где – сила упругости пружины КМ; – сила упругости пружины СЕ, отсюда

- первоначальное сжатие пружины СЕ

Силы, действующие на тело в текущий момент времени t:

– сила тяжести;

– переменная нагрузка, прикладываемая к телу; где Р ₁ = 2,37·10⁶ Н – максимальная нагрузка;
а ₁ = 6,68·10¹⁰ Н/с² – коэффициент; t ₁ = 0,005 с – временной параметр;

– сила упругости пружины КМ,

где – координата точки К в момент времени t; – координата точки К в начальный момент времени; – константа;

– сила упругости пружины CD; –

сила упругости пружины СЕ.

Моменты, действующие на тело, относительно оси CZ ₁:

– момент силы тяжести;

– момент переменной нагрузки;

– момент силы упругости пружины КМ;

– момент спиральной пружины.

Виртуальная работа согласно (2.15)

. (6.7)

Так как обобщённые координаты – независимые друг от друга параметры, то и их вариации тоже независимы. Поэтому, используя принцип замораживания, находим виртуальные работы, соответствующие виртуальным перемещениям поочередно:

1) ,

;

2) , (6.8)

;

3) ,

.

Сравнивая множители в (6.7) перед вариациями обобщённых координат и в формуле (6.8), находим обобщённые силы, соответствующие обобщённым координатам:

;

;

Окончательный вид уравнений Лагранжа второго рода или дифференциальные уравнения движения рассматриваемой системы в обобщенных координатах:

–

,

где

–

где

–

где ,

или в матричной форме ,

где – инерционная матрица, , если , , если ,

,

=

– матрица-столбец обобщенных сил;

– матрица-стол-бец слагаемых, перенесенных из левых частей уравнений Лагранжа, не содержащих обобщенных ускорений.

Решение системы уравнений Лагранжа второго рода, описывающих движение указанной механической системы (рис. 6.5), на интервале времени от t = 0 до t = 0,01 с выполнено в системе Mathcad 11. Распечатка результатов расчёта в системе Mathcad 11 для начального угла j₀ = 30° приведена на рис. 6.6, а для начального угла j₀ = 6° на рис. 6.7.

Рис. 6.6

Рис. 6.6. Окончание

Пример 3. Несвободная система с тремя степенями свободы. Дано. Полый цилиндр 2 массой m ₂ (рис. 6.8) скользит по основанию 1 массой m ₁, опирающемуся на цилиндрический шарнир О и поддерживаемому вертикальной пружиной КМ. Жесткость пружины C ₁, длина недеформированной пружины l₂, расстояние от шарнира О до точки К опоры пружины l ₁. Коэффициент трения цилиндра об основание f = 0,12. K шарниру C прикреплена горизонтальная пружина CD жесткостью C ₂. В начальном положении пружина CD не деформирована, а пружина КМ поддерживает систему в положении статического равновесия.

Рис.6.8

Введем неподвижную систему координат OXYZ, ось OX которой горизонтальна, и связанную с основанием систему координат OX ₁ Y ₁ Z ₁, ее ось OX ₁ параллельна направляющей основания, по которой скользит цилиндр. Положение центра масс O ₁ основания задается координатами x_{1 O} и y_{1 O} точки O ₁ в системе координат OX ₁ Y ₁ Z ₁, причем y_{1 O} = 0,6 м. Начальное положение цилиндра на основании определяется начальными значениями координат x _{2 O} и y _{2 O} точки O ₂ в системе координат OX ₁ Y ₁ Z ₁. Точка O _н совпадает с начальными положением центра масс цилиндра O ₂.

В начальный момент времени к внутренней поверхности дна цилиндра прикладывается переменная нагрузка P (t), определяемая по зависимости (6.1). При этом цилиндр 2 начинает скользить по основанию 1, вызывая вращение последнего вокруг оси шарнира OZ. Движение цилиндра тормозится реакцией R (t) тормозного устройства 3, приложенной к внешней поверхности дна цилиндра. Величина реакции определяется по формуле (6.2). Моменты инерции основания и цилиндра относительно осей O ₁ z ₁ и O ₂ z ₂ равны J ₁ и J ₂ соответственно.

Составить дифференциальные уравнения движения системы и рассчитать конкретное движение на ЭВМ.

Указание. В качестве обобщенных координат выбрать:
1) перемещение точки C по горизонтальной оси , отсчитанное от положения статического равновесия пружины CD;
2) угол поворота j основания , отсчитанный от горизонтали, и 3) перемещение S центра масс цилиндра 2 по направляющей основания 1 , отсчитанное от его начального положения O _н. Начальное значение j = j₀. Длину недеформированной пружины l ₂ определить из условия статического равновесия системы в начальный момент времени. Рекомендуемые значения величин приведены в табл. 6.3.

Исходные данные: m ₁ – масса основания; m ₂ – масса цилиндра. Число степеней свободы: i = s = 3. За обобщенные координаты приняты следующие параметры: ; ; . Обобщенные скорости: ; ; . В рассматриваемую механическую систему входят: основание 1 и цилиндр 2, совершающие плоскопараллельное движение относительно неподвижной системы координат OXYZ,

; ;

; , i= 1,2,

где

Координаты центра массы основания в матричной форме:

(1), или

, где , т.е.

Скорости центра массы основания в матричной форме:

, или

;

;

Кинетическая энергия основания в функции обобщенных координат и обобщенных скоростей:

Координаты центра массы цилиндра в матричной форме:

.

Скорости центра массы цилиндра 2 в матричной форме:

, или

;

.

их квадраты:

Кинетическая энергия цилиндра в функции обобщенных координат и обобщенных скоростей:

Кинетическая энергия системы

Уравнения Лагранжа второго рода для данной системы:

Для :

;

окончательно левая часть уравнения по

;

Для :

окончательно левая часть уравнения по

;

Для :

окончательно

Перепишем уравнения Лагранжа второго рода, оставив в левой части только слагаемые со вторыми производными от обобщенных координат. Все остальные слагаемые перенесем в правую часть уравнения, обозначив их как

Уравнения Лагранжа второго рода в матричной форме:

,

где – инерционная матрица, , если , , если ,

или ,

инерционные коэффициенты:

;

;

;

;

.

Виртуальная работа сил, действующих на рассматриваемую систему:

;

.

Так как обобщённые координаты - независимые друг от друга параметры, то и их вариации тоже независимы. Поэтому, используя принцип замораживания, находим виртуальные работы, соответствующие виртуальным перемещениям поочередно:

1) ,

;

2) ,

3) ,

.

Сравнивая множители в выражениях виртуальных работ перед вариациями соответствующих обобщённых координат и в формуле (6.7), находим обобщённые силы, соответствующие обобщённым координатам:

;

.