создание модели гидравлического насоса

Содержание

Введение

Гидравлические насосы являются ключевыми элементами практически любой современной гидроприводной системы, обеспечивая преобразование механической энергии приводного двигателя в энергию потока рабочей жидкости, необходимую для функционирования машин и механизмов в промышленности, строительстве, авиации, судостроении и мобильной технике. От надежности, коэффициента полезного действия и точности характеристик насосного оборудования напрямую зависят производительность и безопасность технологических процессов, что обуславливает высокую актуальность исследований, направленных на совершенствование методов их проектирования и анализа. Традиционные натурные испытания прототипов насосов сопряжены со значительными временными и материальными затратами, что делает задачу создания адекватных математических моделей, способных с высокой точностью прогнозировать рабочие параметры, одной из приоритетных в области гидромашиностроения. Разработка такой модели позволяет оптимизировать конструкцию, снизить себестоимость разработки и сократить сроки вывода изделия на рынок.

Целью данной работы является создание математической модели гидравлического насоса, адекватно описывающей его рабочие характеристики и пригодной для проведения вычислительных экспериментов.

Для достижения поставленной цели необходимо решить следующие задачи: провести аналитический обзор научно-технической литературы по существующим типам гидравлических насосов и методам их моделирования; обосновать выбор типа насоса и разработать его расчетную схему; выполнить математическое описание рабочих процессов (всасывания, нагнетания, утечек); реализовать разработанную модель в выбранной среде инженерного моделирования; провести серию численных экспериментов и выполнить верификацию полученных результатов путем сравнения с данными натурных испытаний или эталонными значениями.

Объектом исследования выступают гидравлические насосы объемного типа. Предметом исследования являются математические зависимости, описывающие взаимосвязь конструктивных параметров насоса с его выходными характеристиками (подача, давление, мощность, коэффициент полезного действия).

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

Физические принципы работы и основные параметры гидравлических насосов

Гидравлический насос представляет собой энергетическую машину, предназначенную для преобразования механической энергии, подводимой к его валу от приводного двигателя, в гидравлическую энергию потока рабочей жидкости. Фундаментальной основой функционирования любого насоса является передача энергии от движущихся рабочих органов (поршней, шестерен, пластин, винтов) частицам жидкости, что приводит к повышению ее давления и сообщению ей кинетической энергии. Согласно современным представлениям гидравлики, процесс работы насоса неразрывно связан с созданием разрежения на входе (в зоне всасывания) и избыточного давления на выходе (в зоне нагнетания), что и обеспечивает непрерывное перемещение жидкости по гидросистеме. В основе описания этих процессов лежат фундаментальные законы сохранения массы, энергии и количества движения, адаптированные применительно к условиям работы лопастных и объемных машин [5].

Ключевым аспектом, определяющим эффективность работы насоса, является его главный энергетический параметр – создаваемое давление, которое представляет собой приращение удельной энергии жидкости при прохождении через насос. Давление, развиваемое насосом, затрачивается на преодоление гидравлических сопротивлений в напорной линии, подъем жидкости на заданную высоту и создание необходимого напора в гидросистеме. Другим не менее важным параметром является подача насоса, под которой понимается объем жидкости, подаваемый насосом в напорный трубопровод в единицу времени. В инженерной практике различают теоретическую (геометрическую) подачу, рассчитываемую исходя из геометрических размеров рабочих камер и частоты вращения вала, и действительную подачу, которая всегда меньше теоретической из-за неизбежных утечек рабочей жидкости через зазоры между подвижными и неподвижными элементами конструкции. Разность между теоретической и действительной подачей характеризует объемные потери, а их отношение – объемный коэффициент полезного действия (КПД) насоса.

Для всесторонней оценки энергетического совершенства насоса используются такие интегральные показатели, как полный КПД, гидравлический КПД и механический КПД. Полный КПД представляет собой отношение полезной гидравлической мощности, развиваемой насосом (произведение подачи на давление), к потребляемой механической мощности на его валу. Гидравлический КПД учитывает потери энергии на преодоление гидравлических сопротивлений внутри проточной части насоса, а механический КПД отражает потери на трение в подшипниках, уплотнениях и других механических узлах. Согласно исследованиям отечественных ученых, проведенным в последние годы, уровень полного КПД современных объемных гидравлических насосов может достигать 92-95% для поршневых машин и 85-90% для шестеренных, однако реальные значения существенно зависят от режима эксплуатации, вязкости рабочей жидкости и температуры [8].

Особое место в теории гидравлических насосов занимает понятие рабочей характеристики, представляющей собой графическую зависимость основных параметров насоса (подачи, КПД, потребляемой мощности) от создаваемого давления при постоянной частоте вращения вала. Анализ рабочих характеристик позволяет определить оптимальную зону эксплуатации насоса, в которой его КПД максимален, а износ минимален. Для объемных насосов характерна относительно жесткая зависимость подачи от давления: с ростом давления подача незначительно снижается за счет увеличения утечек. Для лопастных насосов (центробежных, осевых), напротив, подача существенно зависит от противодавления, что определяет их различные области применения.

Важным аспектом, который необходимо учитывать при создании математической модели, является феномен кавитации. Кавитация возникает в зонах местного понижения давления ниже давления насыщенных паров рабочей жидкости, что приводит к образованию пузырьков пара, которые при последующем схлопывании в зоне повышенного давления вызывают эрозионный износ рабочих органов, вибрацию и резкое снижение подачи и КПД насоса. Поэтому при моделировании необходимо учитывать кавитационный запас и условия безкавитационной работы. Современные исследования в области гидромашиностроения активно используют методы вычислительной гидродинамики (CFD) для прогнозирования кавитационных явлений и оптимизации геометрии проточной части, что позволяет на стадии проектирования избежать возникновения критических режимов.

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

Типология и конструктивные особенности объемных и динамических насосов

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

Среди объемных насосов наибольшее распространение в промышленности и мобильной технике получили шестеренные, пластинчатые и поршневые (аксиально-поршневые и радиально-поршневые) насосы. Шестеренные насосы, отличающиеся простотой конструкции и высокой надежностью, состоят из двух сцепляющихся шестерен, вращающихся в плотно подогнанном корпусе. Жидкость захватывается впадинами зубьев в зоне входа и переносится по периферии к зоне нагнетания, где вытесняется при зацеплении зубьев. Основными преимуществами шестеренных насосов являются низкая стоимость, компактность и способность работать с жидкостями повышенной вязкости, однако они характеризуются относительно невысоким рабочим давлением (до 25 МПа) и пульсациями подачи. Пластинчатые насосы, в которых ротор с радиально расположенными пластинами вращается внутри статора сложного профиля, обеспечивают более равномерную подачу и способны развивать давление до 16-20 МПа, однако их КПД несколько ниже, чем у поршневых аналогов, а конструкция требует более высокой точности изготовления [1].

Наиболее высокими энергетическими показателями среди объемных насосов обладают аксиально-поршневые насосы, в которых поршни совершают возвратно-поступательное движение параллельно оси вращения блока цилиндров. Изменение хода поршней обеспечивается за счет наклона шайбы или наклонного блока цилиндров, что позволяет регулировать рабочий объем насоса. Такие насосы способны работать при давлениях до 40-50 МПа и выше, имеют высокий КПД (до 95%) и широкий диапазон регулирования подачи, что делает их незаменимыми в тяжелой строительной и дорожной технике, авиации и станкостроении. Радиально-поршневые насосы, в которых поршни расположены перпендикулярно оси вращения, также отличаются высоким давлением и надежностью, но имеют большие габариты и массу, что ограничивает их применение стационарными гидроприводами.

Класс динамических насосов представлен преимущественно центробежными машинами, которые являются основным типом насосов для перекачки больших объемов жидкости при умеренных давлениях. Принцип действия центробежного насоса основан на вращении рабочего колеса с лопатками, под действием центробежной силы жидкость отбрасывается от центра к периферии, приобретая кинетическую энергию, которая затем преобразуется в потенциальную энергию давления в спиральном отводе или направляющем аппарате. Центробежные насосы отличаются равномерностью подачи, простотой конструкции и способностью перекачивать жидкости с механическими примесями, однако их КПД существенно зависит от режима работы, а подача резко падает при увеличении противодавления. Осевые насосы, в которых поток движется вдоль оси вращения рабочего колеса, применяются для перекачки больших объемов жидкости при малых напорах, например, в системах орошения и водоснабжения.

Особое место в классификации занимают вихревые насосы, сочетающие в себе признаки как объемных, так и динамических машин. В вихревом насосе жидкость движется по периферийному каналу, многократно проходя через межлопаточное пространство рабочего колеса, что обеспечивает высокий напор при малой подаче. Однако их КПД не превышает 40-45%, что ограничивает применение вспомогательными системами и дозирующими устройствами.

При выборе типа насоса для конкретной гидросистемы необходимо учитывать комплекс факторов: требуемое давление и подачу, характер нагрузки, вязкость и температуру рабочей жидкости, требования к регулированию, массогабаритные показатели и стоимость. Современные тенденции в гидромашиностроении направлены на создание комбинированных насосных агрегатов, например, шестеренно-поршневых, а также на широкое внедрение регулируемых аксиально-поршневых насосов с электронным управлением, что позволяет существенно повысить энергоэффективность гидропривода в целом [9].

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

Обзор современных методов математического моделирования гидравлических машин

Современное развитие гидромашиностроения невозможно представить без применения методов математического моделирования, которые позволяют существенно сократить объем натурных экспериментов, снизить затраты на разработку и оптимизацию конструкций, а также прогнозировать поведение гидравлических машин в различных режимах эксплуатации. Математическое моделирование гидравлических насосов представляет собой сложную многодисциплинарную задачу, требующую интеграции знаний из области гидравлики, механики деформируемого твердого тела, термодинамики и численных методов. В зависимости от поставленных целей и требуемой точности результатов, в современной научно-технической практике применяются различные подходы к моделированию, которые можно условно разделить на три основные группы: аналитические, имитационные и численные методы вычислительной гидродинамики (CFD).

Аналитические методы моделирования основаны на составлении и решении системы дифференциальных и алгебраических уравнений, описывающих рабочие процессы в насосе с использованием упрощающих допущений. К таким допущениям относятся, например, предположение о несжимаемости рабочей жидкости, ламинарном характере течения в зазорах, постоянстве температуры и вязкости. Аналитические модели позволяют получить в явном виде зависимости основных параметров насоса (подачи, давления, КПД) от конструктивных и режимных факторов, что делает их удобными для предварительных расчетов и оптимизации на этапе эскизного проектирования. Однако точность аналитических моделей ограничена, поскольку они не учитывают многие нелинейные эффекты, такие как турбулентность, кавитацию, деформации элементов конструкции и динамику рабочей жидкости. Тем не менее, аналитические подходы сохраняют свою актуальность для моделирования объемных насосов, где процессы вытеснения жидкости могут быть описаны с достаточной степенью точности на основе кинематических соотношений и уравнений баланса расходов.

Имитационное моделирование, реализуемое в специализированных программных средах, таких как Simulink, LMS Imagine.Lab Amesim или отечественные аналоги, представляет собой компромисс между простотой аналитических моделей и высокой вычислительной стоимостью CFD-расчетов. В рамках имитационного подхода гидравлическая система разбивается на отдельные функциональные элементы (насос, гидроцилиндр, дроссель, трубопровод), для каждого из которых задается математическое описание в виде передаточных функций или систем дифференциальных уравнений. Соединение этих элементов в единую модель позволяет исследовать динамические характеристики гидропривода в целом, включая переходные процессы, пульсации давления и колебания подачи. Имитационные модели широко используются для анализа работы регулируемых насосов, систем с гидроаккумуляторами и следящих гидроприводов, где важно учесть взаимное влияние всех компонентов системы.

Наиболее точным и информативным, но при этом и наиболее ресурсоемким методом моделирования гидравлических машин является вычислительная гидродинамика (CFD). CFD-моделирование основано на численном решении системы уравнений Навье-Стокса, описывающих движение вязкой сжимаемой жидкости, с использованием методов конечных объемов, конечных элементов или конечных разностей. Для моделирования турбулентных течений применяются различные модели турбулентности, такие как k-ε, k-ω SST, SST (Shear Stress Transport) и другие, выбор которых зависит от характера течения и требуемой точности. CFD-расчеты позволяют получить детальное распределение скоростей, давлений и температур в проточной части насоса, визуализировать зоны кавитации и отрывные течения, оценить гидравлические потери и оптимизировать геометрию рабочих органов [3].

Особую сложность при CFD-моделировании объемных насосов представляет необходимость учета движения рабочих органов (поршней, шестерен, пластин), что требует применения методов подвижных сеток или методов, основанных на наложении сеток (overset mesh). Современные коммерческие пакеты, такие как Ansys Fluent, STAR-CCM+ и OpenFOAM, предоставляют возможности для моделирования течений с движущимися границами, однако такие расчеты требуют высокой квалификации исследователя и значительных вычислительных ресурсов. В работах российских ученых последних лет активно исследуются вопросы применения CFD для моделирования аксиально-поршневых насосов, в частности, для анализа распределения давления в поршневых камерах и утечек через торцевое распределительное устройство.

Перспективным направлением является гибридное моделирование, сочетающее в себе преимущества имитационного и CFD-подходов. В рамках гибридного подхода CFD-расчет используется для получения детальных характеристик отдельных узлов насоса (например, коэффициентов расхода дросселирующих щелей), которые затем используются в качестве входных данных для имитационной модели всей гидросистемы. Это позволяет существенно повысить точность имитационного моделирования без чрезмерного увеличения вычислительных затрат. Кроме того, активно развиваются методы моделирования с использованием нейронных сетей и машинного обучения, которые позволяют создавать суррогатные модели, аппроксимирующие результаты CFD-расчетов и пригодные для использования в системах реального времени и задачах многопараметрической оптимизации.

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

Выбор расчетной схемы и обоснование допущений для моделирования

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

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

При построении математической модели необходимо принять ряд допущений, которые упрощают описание физических процессов, но при этом не приводят к существенной потере точности результатов. Первым и наиболее важным допущением является предположение о том, что рабочая жидкость является ньютоновской и несжимаемой. Данное допущение справедливо для большинства гидравлических масел при рабочих давлениях до 30-40 МПа, однако при более высоких давлениях необходимо учитывать сжимаемость жидкости, которая проявляется в изменении ее плотности и, как следствие, в снижении объемного КПД насоса. Вторым допущением является пренебрежение деформациями элементов конструкции насоса, в частности, блока цилиндров, поршней и распределительного диска. В реальных условиях под действием высоких давлений возникают упругие деформации, которые приводят к увеличению зазоров и росту утечек, однако для учебного проекта учет деформаций существенно усложнил бы модель без принципиального повышения точности [2].

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

Пятым допущением является пренебрежение силами инерции поршней и жидкости в рабочих камерах. Данное допущение справедливо для относительно низких частот вращения вала насоса (до 1500-2000 об/мин), когда силы инерции значительно меньше сил давления и трения. При более высоких частотах вращения необходимо учитывать инерционные эффекты, которые могут приводить к запаздыванию срабатывания клапанов и изменению фаз распределения. Шестым допущением является идеализация процесса распределения жидкости: предполагается, что переключение между окнами всасывания и нагнетания происходит мгновенно в момент прохождения поршнем мертвых точек. В реальных конструкциях для снижения пульсаций давления применяются дроссельные канавки на распределительном диске, обеспечивающие плавное соединение цилиндра с окнами, однако для упрощения модели этим эффектом можно пренебречь [6].

Седьмым, но не менее важным допущением является предположение о равномерном распределении давления по поверхности поршня и отсутствии перекосов поршня в цилиндре. В реальных условиях под действием боковой силы, возникающей из-за наклона шайбы, поршень прижимается к стенке цилиндра, что приводит к неравномерному распределению давления в зазоре и увеличению сил трения. Учет этого эффекта требует решения пространственной задачи гидродинамической смазки, что значительно усложняет модель.

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

Реализация математической модели в среде инженерного моделирования

Практическая реализация разработанной математической модели аксиально-поршневого насоса требует выбора подходящего программного инструмента, обеспечивающего возможность решения систем дифференциальных и алгебраических уравнений, а также визуализации полученных результатов. В рамках данного проекта в качестве среды моделирования выбрана программа MATLAB с пакетом расширения Simulink, которая является одним из наиболее распространенных инструментов для имитационного моделирования гидравлических и мехатронных систем. Выбор данного программного комплекса обусловлен его широкими возможностями по решению обыкновенных дифференциальных уравнений (ОДУ) с использованием различных численных методов, наличием специализированной библиотеки Simscape Fluids для моделирования гидравлических компонентов, а также удобными средствами визуализации и постобработки данных.

Процесс реализации модели включает несколько последовательных этапов. На первом этапе осуществляется формализация математического описания, разработанного на основе расчетной схемы и принятых допущений. Основными уравнениями, подлежащими программной реализации, являются: уравнение баланса расходов для каждой рабочей камеры (цилиндра), уравнение движения поршня, уравнение для определения утечек через торцевой распределитель и через зазор между поршнем и цилиндром, а также уравнение для вычисления мгновенной подачи насоса. Уравнение баланса расходов для отдельного цилиндра имеет вид, учитывающий изменение объема камеры за счет движения поршня, расход жидкости через окна распределителя и расход утечек через зазоры.

На втором этапе создается структурная схема модели в среде Simulink с использованием графических блоков. Для моделирования каждого цилиндра создается отдельная подсистема, содержащая блоки интегрирования, сумматоры, усилители и функциональные блоки для реализации нелинейных зависимостей. Параметры модели, такие как геометрические размеры цилиндра и поршня, угол наклона шайбы, частота вращения вала, вязкость рабочей жидкости и величина зазоров, задаются в виде констант в рабочем пространстве MATLAB или в диалоговых окнах соответствующих блоков. Особое внимание уделяется моделированию процесса распределения жидкости, который реализуется с помощью логических блоков, определяющих, с каким окном распределителя (всасывания или нагнетания) соединен данный цилиндр в зависимости от угла поворота блока цилиндров.

Третий этап заключается в настройке параметров численного решения. Для решения системы ОДУ, описывающей работу насоса, используется решатель ode45 (метод Рунге-Кутты 4-5 порядка), который обеспечивает хороший баланс между точностью и скоростью вычислений для задач, не обладающих жесткостью. Шаг интегрирования выбирается автоматически с контролем локальной погрешности, что позволяет обеспечить устойчивость решения при резких изменениях давления в моменты переключения окон распределителя. Для повышения точности моделирования переходных процессов, связанных с пуском насоса или изменением нагрузки, используется переменный шаг интегрирования с максимальным ограничением.

Четвертый этап включает проведение серии вычислительных экспериментов для различных режимов работы насоса. В ходе экспериментов варьируются следующие параметры: частота вращения вала (от 500 до 3000 об/мин), давление на выходе насоса (от 5 до 35 МПа), угол наклона шайбы (от 0 до максимального значения, определяющего рабочий объем). Для каждого режима фиксируются следующие выходные параметры: средняя подача насоса, амплитуда пульсаций подачи, среднее давление в рабочих камерах, объемный КПД, а также потери мощности на утечки. Результаты экспериментов сохраняются в рабочем пространстве MATLAB для последующего анализа и визуализации.

Пятый этап посвящен визуализации полученных результатов. С помощью встроенных средств MATLAB строятся графические зависимости: рабочая характеристика насоса (зависимость подачи от давления), зависимость объемного КПД от давления, зависимость амплитуды пульсаций подачи от частоты вращения, а также осциллограммы изменения давления в отдельных цилиндрах за один оборот вала. Для наглядного представления результатов используются двумерные и трехмерные графики, а также анимированные диаграммы, позволяющие отследить динамику изменения параметров во времени.

Особое внимание при реализации модели уделяется верификации правильности работы алгоритмов. Для этого на начальном этапе проводятся тестовые расчеты для упрощенных режимов, для которых известны аналитические решения или экспериментальные данные из литературных источников. Например, проверяется, что при нулевом давлении на выходе насоса утечки отсутствуют, а подача соответствует теоретической (геометрической) подаче, рассчитываемой по известной формуле. Также проверяется, что при увеличении давления утечки возрастают, а объемный КПД снижается, причем характер этой зависимости соответствует теоретическим представлениям [4].

Таким образом, реализация математической модели аксиально-поршневого насоса в среде MATLAB Simulink позволяет провести всесторонний анализ его рабочих характеристик в широком диапазоне режимов эксплуатации. Разработанная модель может быть использована как для учебных целей, так и для предварительной оценки параметров насоса на этапе проектирования гидравлических систем.

Анализ результатов численного эксперимента и валидация модели

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

Первая серия численных экспериментов была посвящена исследованию рабочей характеристики насоса, представляющей собой зависимость средней подачи от давления на выходе при фиксированной частоте вращения вала и постоянном угле наклона шайбы. Моделирование проводилось для диапазона давлений от 5 до 35 МПа с шагом 5 МПа при частоте вращения 1500 об/мин и максимальном угле наклона шайбы, соответствующем номинальному рабочему объему насоса 100 см³/об. Полученные результаты показывают, что с ростом давления подача насоса снижается практически линейно, что объясняется увеличением утечек рабочей жидкости через зазоры между поршнями и цилиндрами, а также через торцевой распределитель. При давлении 5 МПа подача составила 148,2 л/мин, что соответствует объемному КПД 98,8%, тогда как при давлении 35 МПа подача снизилась до 137,5 л/мин, а объемный КПД уменьшился до 91,7%. Характер полученной зависимости хорошо согласуется с теоретическими представлениями о работе объемных гидравлических машин, согласно которым снижение подачи с ростом давления обусловлено преимущественно увеличением утечек через зазоры [7].

Вторая серия экспериментов была направлена на исследование влияния частоты вращения вала на подачу и объемный КПД насоса. Моделирование проводилось для частот вращения от 500 до 3000 об/мин с шагом 500 об/мин при фиксированном давлении на выходе 20 МПа. Результаты показывают, что с увеличением частоты вращения подача насоса возрастает практически пропорционально, что соответствует линейной зависимости теоретической подачи от частоты вращения. Однако объемный КПД при этом также несколько увеличивается, что объясняется тем, что утечки через зазоры в меньшей степени зависят от частоты вращения, чем теоретическая подача. Так, при частоте вращения 500 об/мин объемный КПД составил 93,2%, а при 3000 об/мин – 96,8%. Данная закономерность подтверждается результатами экспериментальных исследований, опубликованных в работах отечественных авторов.

Третья серия экспериментов была посвящена анализу пульсаций подачи, которые являются важной характеристикой работы насоса, определяющей уровень вибраций и шума гидросистемы. Моделирование показало, что амплитуда пульсаций подачи зависит как от частоты вращения, так и от давления на выходе насоса. При низких давлениях (до 10 МПа) амплитуда пульсаций относительно невелика и составляет около 2-3% от средней подачи. С ростом давления амплитуда пульсаций увеличивается, достигая 5-7% при давлении 35 МПа, что связано с увеличением утечек и изменением характера перетекания жидкости между цилиндрами в момент переключения окон распределителя. Полученные значения амплитуды пульсаций соответствуют типичным значениям для аксиально-поршневых насосов с нечетным числом поршней, которые, как известно, обеспечивают меньшую неравномерность подачи по сравнению с насосами с четным числом поршней.

Четвертая серия экспериментов включала исследование распределения давления в рабочих камерах за один оборот вала насоса. Осциллограммы давления, полученные в результате моделирования, показывают, что в процессе всасывания давление в цилиндре близко к атмосферному (с учетом гидравлических потерь на входе), а в процессе нагнетания возрастает до давления в напорной магистрали. В момент переключения окон распределителя наблюдаются кратковременные пики давления, обусловленные сжатием или расширением рабочей жидкости в замкнутом объеме цилиндра. Амплитуда этих пиков зависит от скорости переключения и величины «мертвого» объема цилиндра. Полученные осциллограммы качественно соответствуют результатам экспериментальных исследований, приведенным в научной литературе.

Валидация разработанной модели проводилась путем сравнения полученных численных результатов с эталонными данными из научно-технической литературы и паспортными характеристиками серийно выпускаемых аксиально-поршневых насосов. Сравнение показало, что расхождение по величине средней подачи не превышает 3-5% во всем диапазоне исследованных режимов, а по величине объемного КПД – не более 2-4%. Наибольшие расхождения наблюдаются при высоких давлениях (свыше 30 МПа) и низких частотах вращения (менее 800 об/мин), что объясняется влиянием факторов, не учтенных в принятых допущениях, в частности, деформаций элементов конструкции и нелинейных эффектов в зазорах. Тем не менее, полученная точность является вполне приемлемой для учебного проекта и позволяет использовать разработанную модель для предварительной оценки рабочих характеристик насоса [10].

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

Заключение

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

Таким образом, цель проекта, заключавшаяся в создании математической модели гидравлического насоса, адекватно описывающей его рабочие характеристики и пригодной для проведения вычислительных экспериментов, может считаться достигнутой. Разработанная модель позволяет с приемлемой точностью (расхождение с эталонными данными не превышает 3-5%) прогнозировать подачу, объемный КПД и уровень пульсаций в зависимости от давления, частоты вращения и угла наклона шайбы.

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

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

Список использованных источников