Математическое моделирование рабочего процесса в вальцово-матричном пресс-грануляторе с торцевым ограничением клиновидного пространства

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

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

Одним из целесообразных путей совершенствования процесса прессового гранулирования является разработка новых и оптимизация существующих конструкций пресс-грануляторов. Однако обзор научных документов дает основание полагать, что за последнее двадцатилетие конструкции отечественных [1], [3] и зарубежных [1], [4], [5], [6] прессующих механизмов с кольцевой матрицей хотя и достигли более высокого технического уровня, но по структурным признакам изменились незначительно.

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

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

Процесс прессования без бокового выдавливания возможен только в условиях, когда контактные поверхности рабочих органов пресс-гранулятора образуют замкнутый контур поперечного сечения слоя продукта. Это может быть обеспечено путем создания дополнительных торцевых контактных поверхностей, позволяющих реализовать в зонах отставания и выдавливания в фильеры матрицы схему плоского деформированного состояния прессуемого материала. Эффективным конструктивным решением является предлагаемое нами устройство (патент РФ № 2 412 819) в виде двух колец, ограничивающих торцевые поверхности слоя прессуемого материала. Боковые поверхности сменных ограничительных колец, обращенные к торцам прессующих роликов и контактируемые с продуктом, а также рабочая поверхность матрицы образуют кольцевую полость с входящими внутрь нее прессующими роликами. Важной задачей исследования многопараметрического нелинейного процесса прессового гранулирования в вальцово-матричном пресс-грануляторе является его адекватное математическое описание. Математическая модель процесса прессования должна отражать энергетическое состояние системы «пресс-гранулятор — обрабатываемый материал», которое выступает в качестве ее внутренней характеристики. В пресс-грануляторах с торцевым ограничением клиновидного пространства между рабочими органами механическая энергия передается прессуемому материалу через контактные поверхности матрицы, прессующих роликов и ограничительных колец. Мерой передаваемой энергии целесообразно считать напряженное состояние элементарного объема продукта на этих поверхностях, которое является основой для определения параметров эффекта, характеризующих процесс гранулирования. Таким образом, необходимо описать напряженное состояние прессуемого материала в клиновидном пространстве между рабочими органами пресс-гранулятора. Прессуемый материал в клиновидном пространстве будем считать сплошной изотропной средой, обладающей свойствами упругопластического тела. Возможность распространения данной модели прессуемого материала на процесс гранулирования растительного сырья и преимущества ее применения обоснованы различными авторами [8], [9], [10]. При этом, рассматривая пластическую деформацию продукта, можно исходить как из теории малых упругопластических деформаций (деформационной теории пластичности), так и из теории пластического течения.

Так как обрабатываемый продукт является уплотняемым (сжимаемым) материалом, то дифференциальное уравнение неразрывности движения прессуемого материала в клиновидном пространстве имеет вид:

(1).

Где — плотность продукта в некоторой точке клиновидного пространства, заданной радиусом-вектором (напряженное состояние в данной точке известно), кг/м3;

— время, с;

— поле скоростей движения частиц прессуемого материала.

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

Уже при относительной плотности порядка 0,8 и сравнительно небольших значениях среднего давления (что справедливо для процесса прессового гранулирования) поведение уплотняемого материала мало отличается от поведения несжимаемых тел, что допускает использование методов теории пластичности несжимаемых тел [12].

Основываясь на вышесказанном, для описания напряженного состояния прессуемого материала примем допущение:

(2).

где — средняя плотность прессуемого материала в активной части клиновидного пространства, кг/м3.

В таком случае будет справедливо условие постоянства элементарного объема пластически деформируемого материала. С учетом (2) полная производная. Тензор градиента скорости обращается в нуль, и уравнение неразрывности (1) запишется следующим образом:

. (3).

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

. (4).

Учитывая, что и, из уравнения (4) получаем в общем виде дифференциальное уравнение движения продукта в клиновидном пространстве:

прессование гранулятор торцевой клиновидный.

(5).

где — тензор напряжений в прессуемом материале (составленный для схемы плоского деформированного состояния).

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

(6).

описывающее параметры напряженного состояния в прессуемом материале во всех точках замкнутого клиновидного пространства.

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

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

Рисунок 1 — Схема взаимодействия прессуемого материала с рабочими органами пресс-гранулятора в условиях торцевого ограничения клиновидного пространства: 1 — кольцевая матрица; 2 — прессующий ролик; 3 ограничительное кольцо

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

(7).

Где — межосевое расстояние, м;

— эксцентриситет, м;

— относительный радиус прессующего ролика (- радиус цилиндрической контактной поверхности обечайки прессующего ролика, м).

При этом угол, отсчитываемый от линии (соединяющей центры матрицы и эксцентриковой оси) в направлении, противоположном вращению матрицы, находится в пределах:

Выделим в зоне отставания элементарный объем прессуемого материала, ограниченный цилиндрическими поверхностями матрицы 1 и прессующего ролика 2, контактными поверхностями ограничительных колец 3, а также двумя радиальными сечениями с элементарным углом между ними. Границу, на которой прекращается контакт торцевых поверхностей ограничительных колец со слоем материала, зададим полярным уравнением:

(8).

Где — угол поворота кольцевой матрицы, рад;

— радиальная высота слоя прессуемого материала в клиновидном пространстве, м;

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

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

Если, то:

(9).

Где — угловая координата радиального сечения, где начинается зона отставания и возникает контакт слоя материала с прессующим роликом, рад;

— угловая координата радиального сечения, где заканчивается зона опережения и прекращается контакт слоя материала с прессующим роликом, рад.

Если, то:

(10).

где — угол поворота матрицы, при котором образуется замкнутое клиновидное пространство, рад.

. (11).

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

. (12).

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

В случае, когда торцевая контактная поверхность ограничительного кольца перпендикулярна контактной поверхности матрицы,.

. (13).

Когда торцевая контактная поверхность ограничительного кольца расположена под углом, рад, к контактной поверхности матрицы, находится по формуле:

. (14).

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

(15).

Где — радиус производящего круга циклоиды, м;

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

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

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

(16).

Где — нормальное радиальное напряжение в продукте, Па;

, — касательные тангенциальные напряжения соответственно на контактных поверхностях матрицы, прессующих роликов и ограничительных колец, Па;

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

— рабочая ширина кольцевой матрицы, м;

— коэффициент, определяемый конфигурацией ограничительных колец: контактная поверхность вертикальная —, наклонная —, циклоидальная — .

Для зоны опережения дифференциальное уравнение будет отличаться от уравнения (16) знаками контактных касательных напряжений.

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

Поскольку где фактор трения равен единице, то для плоского деформированного состояния цилиндрическое условие пластичности, которое является точным при, запишется в выбранной системе координат следующим образом [13], [15]:

. (17).

Тогда уравнение (16) можно преобразовать к виду:

(18).

где — для зоны отставания; - в нейтральном сечении зоны выдавливания в фильеры (где касательные тангенциальные напряжения равны нулю); - для зоны опережения.

Ранними исследователями [9], [10] установлено, что предел текучести при сдвиге при гранулировании растительного материала представляет собой функцию от напряжения всестороннего сжатия (гидростатического напряжения), Па:

(19).

Где — предел текучести сжатия продукта при атмосферном давлении, Па;

— коэффициент всестороннего давления, Па-1.

Параметры и зависят от состава продукта, температуры, влажности и определяются экспериментально.

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

. (20).

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

Рисунок 2 — Поперечное сечение слоя прессуемого материала в замкнутом клиновидном пространстве пресс-гранулятора (,): 1 — кольцевая матрица; 2 — прессующий ролик; 3 — ограничительное кольцо Из выражений (17) и (20) имеем:

. (21).

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

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

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

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

(22).

Где — плотность продукта в зоне выдавливания в фильеры матрицы, кг/м3;

— число прессующих роликов;

— расход материала через фильеры матрицы, м3/с.

В условиях незамкнутого клиновидного пространства между рабочими органами пресс-гранулятора расход может быть представлен в виде:

(23).

Где — расход материала через радиальное сечение (высотой) клиновидного пространства, соответствующее началу зоны отставания, м3/с;

— расход материала через радиальное сечение (высотой) клиновидного пространства, определяющее границу между зоной выдавливания в фильеры матрицы и зоной опережения, м3/с;

— расход материала через сечение зоны выдавливания в фильеры матрицы плоскостью, совпадающей со свободной торцевой поверхностью слоя, м3/с;

— расход материала через сечение зоны отставания плоскостью, совпадающей со свободной торцевой поверхностью слоя, м3/с.

Из выражения (23) следует, что общая производительность пресс-гранулятора, определяемая расходом, может быть увеличена, во-первых, за счет предотвращения бокового выдавливания материала и, во-вторых, за счет увеличения разности расходов и .

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

. (24).

Для замкнутого клиновидного пространства на основании уравнения (3) можно считать, что:

(25).

где — расход материала через радиальное сечение (высотой) клиновидного пространства, соответствующее началу зоны выдавливания в фильеры матрицы, м3/с.

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

Объемные расходы и при пластическом течении материала в клиновидном пространстве можно определить интегрированием элементарных потоков. Тогда без учета величины зазоров и конфигурации ограничительных колец имеем:

(26).

где — скорость элементарного слоя материала в граничных сечениях зоны выдавливания в фильеры матрицы, м/с.

Считаем, что окружные скорости и, м/с, точек контактных поверхностей матрицы и прессующего ролика равны между собой:

. (27).

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

(28).

Где — угловая скорость кольцевой матрицы относительно водила прессующих роликов, рад/с;

— координата, определяющая положение элементарного слоя прессуемого материала, м.

Преобразовав выражение (26) с учетом (28) и подставив его в формулу (22), получим:

(29).

где — средняя высота зоны выдавливания в фильеры матрицы, м.

. (30).

Переменные и (см. рисунок 1) можно выразить через угловые координаты верхней, рад, и нижней, рад, границ зоны выдавливания:

(31).

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

Из выражений (29) и (31) следует, что производительность пресс-гранулятора зависит от границ зоны выдавливания в фильеры матрицы. Определим границы зоны выдавливания в условиях торцевого ограничения клиновидного пространства.

Перепишем уравнение (18) относительно переменной :

(32).

В уравнении (32) сумму, заключенную в квадратные скобки, обозначим как функцию конструктивных параметров рабочих органов пресс-гранулятора. Тогда с учетом (19) и (21), разделив переменные, получим:

. (33).

Решение уравнения (33) для зоны отставания () и для упругопластической части зоны опережения () (нормальным радиальным напряжением в сечении с высотой слоя пренебрегаем) может быть записано в виде функции с промежуточным переменным параметром — определенным интегралом с переменным верхним пределом. При этом должны быть учтены условия (9), (10). Решение для зоны отставания имеет вид:

. (34).

Решение для зоны опережения:

. (35).

Для нахождения и (численными методами) в выражения (34), (35) необходимо подставить граничные значения пределов интегрирования:

(36).

. (37).

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

Анализ уравнений (36), (37) и других зависимостей, полученных на основе уравнения (32), показывает, что за счет введения дополнительных контактных поверхностей ограничительных колец (функция) при том же сопротивлении фильер обеспечивается смещение зоны отставания в направлении, противоположном вращению рабочих органов, с одновременным расширением границ зоны выдавливания. Это означает, что увеличивается разность и в формуле (29) и, соответственно, разность расходов и в правой части выражения (24). Следовательно, увеличивается производительность пресс-гранулятора.

На основании изложенного можно сделать вывод, что боковое выдавливание материала является важным недостатком в рабочем процессе всех существующих пресс-грануляторов с кольцевой матрицей и его необходимо учитывать с целью повышения производительности и снижения энергоемкости вновь проектируемых прессующих механизмов. Теоретически боковое выдавливание характеризуется различной интенсивностью в зависимости от напряженного состояния, вида и реологических свойств обрабатываемого материала, конструктивных параметров пресса. В то же время при гранулировании растительного сырья (различных рецептов комбикормов, подсолнечного жмыха и шрота) в ходе экспериментальных исследований значительное боковое выдавливание наблюдалось нами как на лабораторном пресс-грануляторе ПГМ-05, так и на промышленной установке Б6-ДГВ. Таким образом, устройство в виде элементов рабочих органов — сменных ограничительных колец, предотвращающих поперечную деформацию продукта и улучшающих условия его поступления в зону отставания, — является эффективным конструктивным решением, учитывающим особенности прессования в замкнутом клиновидном пространстве.

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

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

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

• 1. Технологическое оборудование предприятий отрасли (зерноперерабатывающие предприятия) / Л. А. Глебов [и др.]. — М.: ДеЛи принт, 2006. — 816 с.
• 2. Handbook of powder technology Vol. 11 Granulation / edited by A.D. Salman, M.J. Hounslow, J.P.K. Seville. — Oxford: Elsevier, 2007. — 1402 p.
• 3. Жислин, Я. М. Дробильное и прессующее оборудование комбикормового завода / Я. М. Жислин, Б. И. Пикус. — М.: Агропромиздат, 1987. — 118 с.
• 4. Шестерина, С. А. Гранулирование комбикормов за рубежом / С. А. Шестерина. — М.: ЦНИИТЭИ Министерства хлебопродуктов СССР, 1989. — 28 с.
• 5. Пат. 2 412 819 Российская Федерация, МПК8 B 30 B 11/20, B 28 B 3/18. Пресс-гранулятор / Ковриков И. Т., Кириленко А. С.; заявитель и патентообладатель Государственное образовательное учреждение высшего профессионального образования «Оренбургский государственный университет». — № 2 009 145 789/02; заявл. 09.12.2009; опубл. 27.02.2011, Бюл. № 6. — 10 с.: ил.
• 6. Винников, Г. А. Исследование процесса гранулирования комбикорма в прессах с вращающейся кольцевой матрицей: дис. … канд. техн. наук / Г. А. Винников. — М., 1970. — 164 с.
• 7. Технологическое оборудование предприятий по хранению и переработке зерна / под ред. А. Я. Соколова. — М.: Колос, 1984. — 445 с.
• 8. Полищук, В. Ю. Исследование процесса гранулирования комбикормов и путей повышения долговечности рабочих органов прессов-грануляторов: дис. … канд. техн. наук: 05.02.14 / В. Ю. Полищук. — М., 1979. — 192 с.
• 9. Друянов, Б. А. Прикладная теория пластичности пористых тел / Б. А. Друянов. — М.: Машиностроение, 1989. — 168 с.
• 10. Унксов, Е. П. Инженерная теория пластичности: методы расчета усилий деформирования / Е. П. Унксов. — М.: Машгиз, 1959. — 328 с.
• 11. Теория пластических деформаций металлов / под ред. Е. П. Унксова, А. Г. Овчинникова. — М.: Машиностроение, 1983. — 598 с.
• 12. Сторожев, М. В. Теория обработки металлов давлением / М. В. Сторожев, Е. А. Попов; под ред. Е. П. Унксова. — М.: Высш. шк., 1963. — 390 с.