Вычисление производных от скоростей отнесённых к угловой скорости кривошипа
Вычисление производных Ω'k (k=1,2,3)
Формирование производных Ω'k (k=1,2,3) в виде векторов
Блок вычисления приведённого момента внешних сил и углового ускорения звена
Вычисление производной Jпр'
Формирование векторов сил тяжести звеньев механизма
Вычисление коэффициентов в выражении приведённого момента внешних сил
Отобразим изменение вычисляемых величин на графике за один оборот кривошипа
Вычисление углового ускорения кривошипа ОА
Конечный момент времени
Процедура интегрирования дифференциальных уравнений
Конечный момент времени
Количество узловых точек
Формирование системы дифференциальных уравнений
Вывод результатов вычислений
Применение процедуры Radau
Формирование векторов искомых величин
Вычисление средней угловой скорости
График изменения угловой скорости ω=ω(t) и величины ωср
Вычисление среднего углового ускорения εср
График углового ускорения ε и εср
Вычисление MД и его среднего значения
График Мд и Мдср
Процедура вычисления реакций внешних и внутренних связей
Вычисление угловых ускорений звеньев механизма
Формирование вектора правых частей системы уравнений
Нахождение угловых ускорений звеньев
Вычисление ускорения узловых точек и центров масс звеньев
Формирование матрицы коэффициентов V при неизвестных реакциях и вектора правых частей H для каждого момента времени задаваемого вектором решений t
Графики реакций внешних и внутренних связей
5. Анализ результатов вычислений
Анализ результатов вычислений позволяет сделать следующие выводы
1. Время неустановившегося движения механизма невелико и составляет
около 1.3 с.
2. В установившемся режиме движение кривошипа близко к равномерному вращению, средняя угловая скорость которого порядка
Максимальные и минимальные значения угловой скорости в установившемся режиме приблизительно равны
и , а его период - 0.162 с. Таким образом, коэффициент неравномерности движения механизма приблизительно равен3. В установившемся режиме среднее угловое ускорение маховика приблизительно равно
. Амплитуда изменения углового ускорения значительна и составляет около , а коэффициент динамичности в этом случае