A method for calculation of oil supply to the parts of units of agricultural machinery

Full Text

Введение Двигатели и агрегаты трансмиссий сельхозмашин [1, 2] могут быть оборудованы принудительными системами смазки, в которых масло под действием масляных насосов подается к потребителям, снижая непроизводительные затраты на трение и обеспечивая необходимые тепловые режимы деталей. Для расчета систем смазки при проведении проектировочных и конструкторско-доводочных работ используются математические модели, которые представляют собой сложные системы, состоящие из большого количества нелинейных уравнений. Такие системы требуют разработки специальных методов решения. Цель исследования Цель исследования - разработать метод решения системы линейных и нелинейных уравнений, описывающих движение масла по каналам принудительных систем смазки. Метод должен обеспечивать простоту, однозначность и высокую достоверность расчетов. Материалы и методы В работах [3-5] предложена следующая математическая модель системы смазки: где Здесь для каждого i-го участка подачи масла: - давления на входе и выходе, Па; - давление, обусловленное действием центробежных сил (для радиально вращающегося участка), Па; - коэффициент неравномерности потока масла (коэффициент Дарси); - суммарный коэффициент гидравлических сопротивлений; - среднее арифметическое значение температуры масла, К; - среднее арифметическое значение давления, Па; - массовый расход масла, кг/с; - площадь поперечного сечения, м2; - плотность газомасляной смеси, кг/м3; - частота вращения (для радиально вращающегося участка), рад/с; - плотность масла, кг/м3; - угол наклона к оси вращения (для радиально вращающегося участка), рад; - радиусы входа и выхода масла (для радиально вращающегося участка), м; - коэффициент местных гидравлических сопротивлений; - коэффициент, учитывающий потери на трение масла; - длина участка, м; - гидравлический диаметр участка, м; - коэффициент, учитывающий растворимость газа в масле от температуры, К [5]; - коэффициент, учитывающий растворимость газа в масле от давления, Па [5]; - универсальная газовая постоянная, Дж/(моль·К); - давление равновесной растворимости, Па; - температура равновесной растворимости, К. Данная математическая модель представляет собой сложную систему, состоящую из линейных и нелинейных уравнений, для решения которой необходимо применить специальные методы. Следует отметить, что входящие в уравнения параметры зависят от температуры, но для повышения точности расчета необходимы дополнительные уравнения, позволяющие определять текущие температуры деталей и масла. В предлагаемом методе решения текущие температуры деталей и масла и задаются как средняя температура масла в картере рассчитываемой сборочной единицы, следовательно, все зависящие от нее параметры рассматриваются как постоянные величины. В результате анализа влияния рассчитываемых величин друг на друга выделены три основные зависимости. Во-первых, коэффициент потерь на трение для каждого из участков линий подач представляет собой функцию от объемного расхода масла. Во-вторых, плотность масла представляет собой функцию от давления на входе-выходе участка линии подачи. В-третьих, давление, обусловленное действием центробежных сил, представляет собой функцию от плотности смазочного материала на входе во вращающийся участок. Для решения системы уравнений предлагается применять итерационный метод, включающий три вида итераций: 1) решение системы уравнений с использованием предварительно заданных значений коэффициентов потерь на трение и плотностей смазочного материала для каждого из выделенных участков линий подач смазочного материала; 2) уточнение коэффициентов потерь на трение по полученным значениям подач смазочного материала и повторное решение системы уравнений; 3) уточнение значений плотности смазочного материала по полученным значениям давления на входе-выходе участков линий подач смазочного материала и повторное выполнение пунктов 1 и 2. Итерации выполняются до тех пор, пока не будет достигнута необходимая точность расчета. При выполнении итераций 1 решается система квадратных и линейных уравнений. Для этого предлагается применить метод Гаусса [6], используемый для решения систем линейных уравнений. На первом шаге итерации задаются приближенные величины подач для каждого из выделенных участков линий подач системы смазки, при этом квадратные уравнения принимают вид: где , - предварительно заданная и искомая величины подач. Система уравнений (1) превращается в систему линейных уравнений. После ее решения проводится сравнение заданных и искомых величин подач. Если различия между ними не превышают необходимую точность, то расчет повторяется. Но для последующего шага итерации предварительно заданная величина подачи определяется как среднее арифметическое между заданной величиной для предыдущей итерации и полученной в результате ее выполнения величиной подачи: В результате для решения системы квадратных уравнений удалось применить метод решения системы линейных уравнений. Наиболее удобным методом решения системы линейных уравнений оказался метод Гаусса. Итерацию 2 предложено проводить в следующей последовательности. Коэффициенты потерь на трение определяются исходя из предварительно заданных величин подач. После выполнения итерации 1 коэффициенты потерь на трение определяются по полученным величинам подач и сравниваются с предыдущими. Если их различия превышают заданную точность, то расчет повторяется с использованием полученных коэффициентов. Наибольшую сложность представляет выполнение итерации 3, так как она обусловлена зависимостью плотности смазки не только от давления на входе-выходе участков линий подач, но и от давления, обусловленного действием центробежных сил в радиально вращающихся каналах. В свою очередь, давление, обусловленное действием центробежных сил, зависит от плотности смазки. Проведен анализ методов решения, на основании которого выбран метод, обеспечивающий достоверное решение для рассматриваемой задачи. На рис. 1, а представлены графики 1-3 изменения давления на входе в радиально вращающийся канал в зависимости от плотности смазочного материала для различных частот вращения, а также график изменения плотности газомасляной смеси в зависимости от давления, обусловленного действием центробежных сил, где . Из графиков на рис. 1, а видно, что применение метода простых итераций обеспечивает сходимость только для каналов, вращающихся с относительно малой частотой. По этой причине для выполнения итерации 3 применен метод секущих [6], обеспечивающий сходимость для рассматриваемой системы уравнений. Геометрическая интерпретация данного метода приведена на рис. 1, б. Для реализации метода секущих необходимо задать два исходных значения плотности и . С их использованием определяются величины давлений , и , . Полученные точки дают возможность определить значение плотности , которая служит исходной границей для последующей итерации. Предлагаемый метод позволяет получить искомые величины при минимальном числе итераций. При проведении расчета система смазки разбивается на отдельные участки, не имеющие разветвлений, характеризующиеся близкими значениями диаметров и наличием вращения. Для каждого участка определяются геометрические размеры, выбираются коэффициенты местных гидравлических сопротивлений, задаются частоты вращения для радиально вращающихся участков, определяются углы их наклона к оси вращения. Для расчета, как правило, выбирается около 10-30 участков. Для реализации предложенного метода решения системы уравнений, описывающих работу принудительных систем смазки, требуется применение средств вычислительной техники. Результаты и их обсуждение С целью проверки достоверности предложенного метода выполнен расчет системы смазки планетарной коробки передач, а также проведены измерения суммарного объемного расхода масла в процессе ее функционирования на различных режимах работы. Результаты расчета массовых G и объемных Q расходов масла приведены на рис. 2. Выводы Полученные результаты указывают на высокое совпадение рассчитанных и измеренных значений объемного расхода масла, что подтверждает достоверность предложенного метода решения системы уравнений (1) и возможность его использования для расчета систем смазки при проведении проектировочных и конструкторско-доводочных работ.

About the authors

O. L Malomyzhev

Irkutsk National Research Technical University

Email: olm@bk.ru
PhD in Engineering Irkutsk, Russia

N. E Fedotova

Branch of Irkutsk National Research Technical University in Usolye-Sibirskoye

Email: oip@istu.edu
PhD in Economics Usolye-Sibirskoye, Russia

V. V Skutel'nik

Irkutsk National Research Technical University

PhD in Engineering Irkutsk, Russia

References

Supplementary files