Динамика плоских шарнирных механизмов (стр. 2 из 5)
где
- ускорение точки А кривошипа ОА, 
- ускорение точки K шатуна АB.Ускорение ползуна B определим дифференцированием уравнения (6)
(11)2.2 Составление дифференциального уравнения движения механизма с помощью теоремы об изменении кинетической энергии
Для составления дифференциального уравнения движения механической системы с одной степенью свободы применим теорему об изменении кинетической энергии в дифференциальной форме
(12)где Т - кинетическая энергия системы;
- сумма мощностей внешних сил; - 
сумма мощностей внутренних сил.Вычислим кинетическую энергию системы как сумму кинетических энергий тел, образующих механическую систему
Т = Т0+Т1+Т2+Т3+Т4,
где
- кинетическая энергия кривошипа ОА, совершающего вращательное движение вокруг оси Oz; 
- кинетическая энергия шатуна АB, совершающего плоское движение; 
- кинетическая энергия шатуна KD, совершающего плоское движение; 
- кинетическая энергия кривошипа О1D, совершающего вращательное движение вокруг оси O1z; 
- кинетическая энергия ползуна B, который движется поступательно.Моменты инерции сплошных однородных стержней, составляющих механизм, относительно осей проходящих через их центры масс равны
Подставляя всё в выражение кинетической энергии системы, окончательно получаем:
(13)где
(14)- приведенный момент инерции механизма, а величины
, (k = 1,4) - скорости точек механизма, отнесенные к угловой скорости ведущего звена 
Для рассматриваемой механической системы, состоящей из абсолютно твердых тел, соединенных идеальными шарнирами сумма мощностей внутренних сил равна нулю
.Сумма мощностей внешних сил будет равна
где
, 
, 
, 
мощности сил тяжести звеньев; 
- мощность момента приводящего механизм в движение; 
- мощность полезной нагрузки.Мощности сил
, 
равны нулю, т.к. реакция опорной плоскости YП и сила тяжести перпендикулярна скорости точки B, а остальные силы приложены к неподвижным точкам.Учитывая выражения для движущего момента МД и полезной нагрузки
, окончательно получим 
(15)где
(16)- приведенный момент внешних сил, а величины
и 
равны 
Подставляя найденные выражение кинетической энергии (13) и мощности внешних сил (15) в теорему об изменении кинетической энергии (12), получим дифференциальное уравнение движения механизма
(17)где
- производная момента инерции механизма по углу поворота ведущего звена.Решив данное дифференциальное уравнение второго порядка с указанными в задаче начальными условиями, найдем закон движения ведущего звена
, его угловую скорость 
и угловое ускорение 
.3. Нахождение реакций внешних и внутренних связей
Для определения реакций внешних и внутренних связей расчленим плоский шарнирный механизм на отдельные звенья и изобразим реакции внешних и внутренних связей каждого звена (рис. 3).
Рис.3 Расчётные схемы звеньев плоского механизма.
Применив к каждому телу, изображенному на расчетной схеме, теорему о движении центра масс (в проекциях на оси координат) и теорему об изменении кинетического момента (для кривошипов относительно осей вращения, для шатунов относительно осей проходящих через центр масс) получим следующую систему уравнений:
Кривошип ОА:
(18)Шатун АB:
(19)Шатун KD:
(20)Кривошип О1D:
(21)Ползун B:
(22)Первое уравнение системы (22) позволяет определить реакцию опорной плоскости YП, а третье из системы (18), после подстановки найденных величин, дифференциальное уравнение движения механизма (17).
Оставшиеся двенадцать соотношений представляют собой систему линейных алгебраических уравнений относительно неизвестных реакций.
4. Результаты расчетов
Решение поставленной задачи сводится к численному интегрированию дифференциального уравнения движения механизма (17) и решению системы двенадцати линейных алгебраических уравнений (18) - (21) относительно неизвестных динамических реакций внешних и внутренних связей.
Задача интегрирования дифференциального уравнения (17) связана с большим количеством предварительных вычислений и может быть условно разбита на пять блоков:
о решение системы уравнений геометрических связей (1) или вычисление геометрических соотношений (2);
о вычисление кинематических соотношений по формулам (3) - (11);
о вычисление приведенного момента инерции механизма и приведенного момента внешних сил.
о вычисление производной от приведенного момента инерции по углу поворота ведущего звена.
о численное интегрирование дифференциального уравнения.
Все это может быть проведено в Mathcad несколькими способами использующими различные встроенные процедуры-функции. Отличие этих способов и методов заключается во времени вычислений, которое требуется для нахождения: решения системы уравнений геометрических связей, приведенного момента инерции механизма и его производной, приведенного момента внешних сил, а также решения дифференциального уравнения движения (17). Ниже рассмотрен алгоритм и приведен документ Mathcad, в котором обеспечивается минимальное время вычислений.
4.1 Алгоритм вычислений
Угловые координаты звеньев
и положение ползуна B 
вычисляются в явном виде по формулам (2).