CC BY
56
http://vestnik.mrsu.ru ISSN Print 2658-4123
ISSN Online 2658-6525
УДК 631.3 DOI: 10.15507/2658-4123.029.201902.279-294
Методика расчета длины резки растений кормоуборочными комбайнами: роторным и двойного измельчения
М. И. Белов
ФГБОУ ВО «Российский государственный аграрный университет - МСХА имени К. А. Тимирязева» (г. Москва, Россия) B-Mikhael@yandex.ru
Введение. Результаты многочисленных исследований доказывают, что качество корма из измельченных растений зависит от качества измельчения, которое обеспечивается кормоуборочным комбайном и определяется путем взвешивания фракций частиц разной длины. Представляет интерес теоретическое обоснование режимов измельчения трав на сенаж и другие виды кормов. Целью данной статьи является разработка алгоритма и методики оценки качества измельчения растений на основе математической модели измельчения кормоуборочным комбайном с роторным измельчающе-швыря-ющим аппаратом и кормоуборочным комбайном двойного измельчения (с двойным измельчающим аппаратом).
Материалы и методы. Изложены необходимые положения теории роторного измельчающего аппарата и двойного измельчающего аппарата кормоуборочного комбайна, на основе которых разработаны алгоритмы, методика расчета массовой доли частиц в заданных диапазонах длины и компьютерная модель измельчения. Результаты исследования. Представлены алгоритмы и методика расчета массовой доли частиц растений, измельченных роторным аппаратом и двойным измельчающим аппаратом. Установлена связь массовой доли частиц в заданных диапазонах длины с высотой растений, высотой среза, подачей, средней расчетной длиной резки. Обсуждение и заключение. Кормоуборочные комбайны с роторным измельчающе-швыряющим аппаратом не позволяют обеспечить заготовку трав на сенаж с массовой долей частиц длиной 8-19 мм в пределах от 45 % до 65-75 %, комбайны с двойным измельчающим аппаратом дают такую возможность только при регулировании угловой скорости шнека и/или дискового или цилиндрического барабана. Двойной измельчающий аппарат позволяет выполнить требования к качеству измельчения разных кормов, если предусмотрено регулирование угловой скорости шнека и/или дискового или цилиндрического барабана.
Ключевые слова: кормоуборочный комбайн, измельчающий аппарат, закон распределения частиц по длине, качество измельчения, длина резки, длина частиц, роторный комбайн, комбайн двойного измельчения
Для цитирования: Белов М. И. Методика расчета длины резки растений кормоуборочными комбайнами: роторным и двойного измельчения // Инженерные технологии и системы. 2019. Т. 29, № 2. С. 279-294. Ш1: https://doi.org/10.15507/2658-4123.029.201902.279-294
© Белов М. И., 2019
|(с?1 © I Контент доступен по лицензии Creative Commons Attribution 4.0 License.
This work is licensed under a Creative Commons Attribution 4.0 License.
The Method of Calculating Cut Length for Flail and Double Chop Forage Harvesters
M. I. Belov
Russian Timiryazev State Agrarian University (Moscow, Russia)
B-Mikhael@yandex.ru
Introduction. Many studies show that the chopped feed quality depends on the plants chopping quality. The quality of plant chopping, provided by a forage harvester, is defined by particle length distribution. it is important to find the theoretical basis of grass chopping regimes for making haylage and other kinds of feeds. The purpose of this study is to develop an algorithm and techniques to assess the quality of plant chopping based on mathematical models of chopping with the flail and double chop forage harvesters. Materials and Methods. Mathematical models of a flail type unit and double chop unit were presented and used to develop the algorithms and method of calculating the cut length distribution.
Results. The algorithms and computation methods of cut length distribution were presented for a flail type unit and a double chop unit. The effect of the plant height, the height of cut, feed and the average length of cut on the mass fraction of particles within the specified ranges of lengths were investigated.
Discussion and Conclusion. A flail forage harvester is not able to harvest the grass crops for making haylage with 45-75 % forage particles, the length of which is in the range from 8 mm to 19 mm. A double chop forage harvester allows harvesting the grass crops for making haylage only when auger angular velocity and/or flywheel or cylinder drum angular velocity is regulated. A double chop device can be used to chop plants for preparation of various forages if auger angular velocity and/or flywheel or cylinder drum's angular velocity is regulated.
Keywords: forage harvester, chopper unit, particle length distribution, grinding quality, cut length, particle length, flail harvester, double chop harvester
For citation: Belov M.I. The Method of Calculating Cut Length for Flail and Double Chop Forage Harvesters. Inzhenernyye tekhnologii i sistemy = Engineering Technologies and Systems. 2019; 29(2):279-294. DOI: https://doi.org/10.15507/2658-4123.029.201902.279-294
Введение
Измельчение — главная операция, выполняемая кормоуборочным комбайном, и основной критерий эффективности его работы — выполнение агротехнических требований по качеству измельчения растений. Агротехнические требования по качеству измельчения растений регламентируют массовые доли фракций частиц в заданных диапазонах длины.
Результаты многочисленных исследований доказывают, что качество корма, изготовленного с помощью кормоуборочного комбайна, зависит от
качества измельчения растений [1-4], упрощенная методика оценки которого предполагает использование сепаратора частиц, содержащего три решета. Современные рекомендации к качеству измельчения трав на сенаж состоят в обеспечении следующих показателей массовой доли частиц: на верхнем решете с круглыми отверстиями диаметром 19 мм - 10-20 %, на среднем решете с круглыми отверстиями диаметром 8 мм - 45-75 %, на нижнем решете с квадратными отверстиями размером 1,3 мм - 10-20 %, под нижним решетом на нижней панели сепарато-
Vol. 29, no. 2. 2019
ра - менее 5 %1. Представляет интерес определение режимов, обеспечивающих требуемое качество измельчения трав на сенаж. Задача определения качества измельчения сводится при этом к нахождению функции распределения массовой доли частиц по длине и может быть решена построением компьютерной модели измельчения растений.
Целью данной статьи является разработка алгоритма и методики оценки качества измельчения растений на основе математической модели измельчения роторным кормоуборочным комбайном и кормоуборочным комбайном двойного измельчения.
Обзор литературы
Качество измельчения растений, оцениваемое весовым содержанием фракций частиц в разных диапазонах длины, можно рассчитать посредством нахождения функции распределения массовой доли частиц по длине. Анализом разделения растений на части при их резании и измельчении кормоубороч-ными комбайнами занимались многие ученые.
Математическое представление функции распределения по длине массовой доли частиц растений, измельченных кормоуборочными комбайнами, с помощью функций нормального и логнор-мального распределения изучено в ряде исследований [5-7]. Кроме того, отмечалось, что распределение массовой доли частиц по длине не является нормальным; предлагалось представлять плотность распределения в виде суммы функций параболического и экспоненциального типа [8]. Японские ученые предложили математическую модель измельчения растений роторным кор-моуборочным комбайном при низком
и высоком стеблестое кукурузы и сорго [9]. Данная модель не учитывает многократное перерезание наклонных стеблей, сдвинутых передним щитком.
Вопрос о функции распределения массовой доли частиц растений, измельченных кормоуборочными комбайнами двойного измельчения, недостаточно изучен.
Измельчающие и измельчающе-швы-ряющие аппараты современных кормо-уборочных комбайнов традиционно разделяются на две категории: аппараты, обеспечивающие точную (расчетную) длину резки, и не обеспечивающие ее2. Признаком первой категории является тип подачи растений в камеру измельчающего или измельчающе-швыряющего аппарата, при котором растения сжимаются и подаются сжатым слоем между вальцами питающего аппарата. Признаком второй категории служит подача растений в камеру измельчителя не вальцами питающего аппарата, а другим способом (как правило, шнеком).
Классифицируя комбайны на рынке в 1980-е гг., С. Хоу и Т. Бишоп отметили семь групп, главным отличием которых друг от друга был тип измельчающего аппарата [10]. Роторные кормоубороч-ные комбайны и комбайны двойного измельчения занимали и занимают нишу машин, используемых в малых хозяйствах на заготовке зеленых кормов [11; 12].
1. Роторный кормоуборочный комбайн оснащен роторным измельчающе-швыряющим аппаратом со швырковыми ножами (рис. 1, а).
2. Роторный кормоуборочный комбайн со швырялкой оснащен роторным измельчающим аппаратом с режущими
1 From harvest to feed: understanding silage management / C. M. Jones [et al.]. Penn State's College of Agricultural Science. 2004. 34 p. URL: http://pss.uvm.edu/pdpforage/Materials/CuttingMgt/ From_Harvest_to_Feed_Understanding_Silage_management_PennState.pdf; Wiersma D. Theoretical length of cut: Theory and practice // Progressive Dairyman Canada. 2013. URL: https://www. progressivedairycanada.com/topics/feed-nutrition/theoretical-length-of-cut-theory-and-practice
2 Engineering principles of agricultural machines. - 2nd ed. / A. K. Srivastava [et al.]. ASAE. 2006.
604 p. URL: https://fcaib.edu.ng/books/Engineering/%5BAjit_K._Srivastava_Carroll_E._Goering%5D_
Engineer(BookFi.org).pdf
a)
b)
Р и с. 1. Схемы кормоуборочных комбайнов с роторным измельчающе-швыряющим аппаратом (а)
и с роторным измельчающим аппаратом в комплекте со швырялкой (b): 1 - противорежущий брус; 2 - ротор; 3 - швыряющий нож; 4 - выгрузной трубопровод; 5 - ротор; 6 - режущий Г-образный нож; 7 - шнек; 8 - швырялка; 9 - лопасть швырялки; 10 - выгрузной трубопровод F i g. 1. A traditional flail forage harvester (a) and a flail forage harvester with a blower (b): 1 - shearbar; 2 - rotary cutter; 3 - throwing knife; 4 - discharge pipe; 5 - rotary cutter; 6 - L- shape knife; 7 - auger; 8 - blower; 9 - blower blade; 10 - discharge pipe
ножами (рис. 1, Ь). Измельченные растения сначала выгружаются в шнек и затем шнеком подаются в швырялку.
3. Кормоуборочный комбайн двойного измельчения оснащен роторным измельчающим аппаратом и дисковым
измельчающе-швыряющим аппаратом (рис. 2).
4. К кормоуборочным комбайнам, не обеспечивающим точную длину резки, относятся комбайны с цилиндрическим измельчающе-швыряющим
a)
b)
Р и с. 2. Кормоуборочный комбайн двойного измельчения с роторным измельчающим аппаратом и дисковым измельчающе-швыряющим аппаратом при виде слева (а) и сверху (b): 1 - роторный измельчающий аппарат; 2 - режущие Г-образные ножи; 3 - дисковый измельчающе-швыряющий аппарат; 4 - ножи дискового аппарата; 5 - швырковые лопатки дискового аппарата; 6 - шнек; 7 - противорежущие брусья; 8 - кожух дискового измельчающе-швыряющего аппарата и выгрузной трубопровод
F i g. 2. A double chop forage harvester with a rotary cutter and with a flywheel cutterhead on the left side (a) and on the top (b): 1 - rotary cutter; 2 - L-shape blades; 3 - flywheel cutter-head; 4 - blades of a flywheel cutterhead; 5 - throwing paddles of a flywheel cutterhead; 6 - auger; 7 - a shearbars; 8 - flywheel cutterhead casе and discharge pipe
a) b) c)
Р и с. 3. Схема кормоуборочного комбайна с цилиндрическим измельчающе-швыряющим аппаратом (а) и комбайна двойного измельчения с роторным аппаратом и цилиндрическим измельчающе-швыряющим аппаратом при виде сбоку (b) и спереди (с): 1 - подборщик; 2 - шнек; 3 - противорежущий брус; 4 - цилиндрический измельчающе-швыряющий аппарат; 5 - нож роторного измельчающего аппарата; 6 - шнек; 7 - цилиндрический измельчающе-швыряющий аппарат; 8 - выгрузной трубопровод
F i g. 3. A forage harvester with a cylinder cutterhead (a) and a double chop forage harvester with a rotary cutter and with a cylinder cutterhead on the left side (a) and from the front (b): 1 - pick up; 2 - auger; 3 - shearbar; 4 - cylinder type cutterhead with throwing blades; 5 - knife of a rotary cutter; 6 - auger; 7 - cylinder type cutterhead with throwing blades; 8 - discharge pipe
аппаратом, в камеру которого растения подаются шнеком (рис. 3, а) [10].
5. Разновидностью комбайнов двойного измельчения является комбайн с роторным измельчающим аппаратом и цилиндрическим измельчающе-швыряющим аппаратом (рис. 3, b-c) [Там же].
Материалы и методы Методика оценки качества измельчения реализована в виде компьютерной модели измельчения растений роторным кормоуборочным комбайном и кормо-уборочным комбайном двойного измельчения, разработанной в среде Lazarus 1.8.4. Математическая модель и алгоритмы расчета фракционного состава частиц измельченных растений изложены далее. В расчетах использованы данные современных комбайнов (табл. 1). Результаты исследования Математическая модель резания растений ножами роторного измельчающего аппарата
Рассмотрим роторный измельча-юще-швыряющий аппарат. Введем следующие обозначения: sm - наибольшая длина отрезанных частиц растения,
срезанного на корню и перерезаемого в третий раз и далее, мм; Ь - заданная длина отрезка, мм; МЬ - масса частиц растений заданной совокупности, длина каждой из которых не превышает Ь, кг; МТ - масса частиц растений заданной совокупности, кг; тЬ - массовая доля частиц растений заданной совокупности, длина каждой из которых не превышает Ь, %; т - массовая доля частиц растений заданной совокупности, длина каждой из которых находится в пределах от 8 до 19 мм, %. Прочие обозначения указаны на рис. 4.
Будем считать, что ножи роторного комбайна измельчают первоначально прямостоящие вертикальные стебли.
Запишем соотношения между параметрами (Н - И), s, R, k [13]:
Я5 = (2п / k )R, (1)
Asm0 = п /k + 0, (2)
где в - параметр, обозначающий угол, показанный на рис. 4, а.
Примем, что начало О неподвижной системы координат расположено на оси
Таблица 1 Table 1
Модели и параметры роторных кормоуборочных комбайнов и комбайнов двойного измельчения Flail forage harvesters and double chop forage harvesters models and parameters
Компания, модель (и) / Company, model (s) Номер рисунка/ Figure no. Радиус R ротора, мм / Rotor radius R, 111111 Число к ножей на роторе в плоскости вращения / Number к of flail knives in rotation plane Угловая скорость со ротора, об/мин / Angular velocity со of a rotor, rpm Скорость 1' комбайна, км/час / Speed i' of a tractor, km/h Диаметр D дискового измель-чающе-швыря-ющего аппарата по краям ножей, мм / Diameter D of a flywheel cut-terhead, 111111 Угловая скорость соj диска, об/мин / Angular velocity со, of a flywheel cutterhead, rpm Число n ножей дискового аппарата / Number n of flywheel cutterhead knives Ширина W приемной горловины дискового аппарата, мм / Flywheel chopper input throat width W, 111111 с
Staalmeester GH 1303 1 н/д / n.a. 1 н/д / n.a. 8-10 - - - -
Hiniker 57104 1 307* 1 1 503* н/д / n.a. - - - -
CelikelTarim 1300; 15005 1 н/д / n.a. 1 1 350 н/д / n.a. - - - -
Корммаш КИР-1,5M / Komimash KIR-1.5M6 1 286 1 1 500 <8 - - - -
FiMAX RF 132; RF 1527 1 н/д / n.a. н/д / n.a. н/д / n.a. н/д / n.a. - - - -
John Deere 72s 2 н/д / n.a. 2 1 366 н/д / n.a. н/д / n.a. 800 3(6) 246
New Holland 389 2 н/д / n.a. н/д / n.a. 1 373 н/д / n.a. н/д / n.a. 966 3(6) н/д / n.a.
Staalmeester DS 54010 2 н/д / n.a. 2 1 655 н/д / n.a. 810 1 258 6 н/д / n.a.
Elho DC2100 S11 2 н/д / n.a. 2 2 000 н/д / n.a. н/д / n.a. 1 080 8 н/д / n.a.
Fimax 1500.210012 2 н/д / n.a. 2 н/д / n.a. н/д / n.a. н/д / n.a. н/д / n.a. 6 н/д / n.a.
Ростсельмаш Стерх 2000 / Rostselmash Sterh 200013 2 н/д / п. a. 1 1 500 н/д / n.a. н/д / n.a. 1 000 6 н/д / n.a.
■ H I
I g
s =
s
43
Ol й/
Примечание: * - расчетное значение; н/д - нет данных / Note: * - calculated data; n.a. - not available
3 Celmak Forage Harvesters. URL: http://www.staalmeester.co.za/Celmak_forage_harvesters.html
4 5710 Flail Mower/Shredder. URL: https://www.lmiiker.con"i/ag_products%20new/5710failinower.hta"il; Flail chopper shredder Model 5710 (with identification no.'s ending 100 and higher): Operator's manual. URL: https://www.lmiiker.conVag_products%20new/maiuMls/Operator's%20MaiuMls/79202099.pdf
5 Forage harvester. URL: http://www.celikeltariin.con"i/eii/forage-harvester.html
6 Косилка-измельчитель роторная КИР-1,5М. Техническое описание и инструкция по эксплуатации. URL: http://agtz.ru/upload/iblock/420/ Rukovodstvo-po-ekspluatatsii-KIR_l_5.pdf
7Fimaks silaj makineleri. URL: http://fiinaks.conVwp-conteniyuploads/2018/03/Silaj-Makiiieleri-K.pdf
8 Hay and Forage Harvesting Equipment. URL: https://www.deere.com/en/liay-forage/liarvesting
9 Rill-type forage harvesters and flail harvester. URL: http://d3ulquraki94yp.cloudfront.net/iiliag/nar/en-us/assets/pdf/liay-tools/pull-type-forage-harvester-brochure-us-en.pdf
10 Celmak Forage Harvesters. URL: http://www.staalnieester.co.za/Celniak_forage_harvesters.htnil
11 Forage Harvesters. URL: https://www.elho.fi/products/forage-harvesters
12 Fimaks silaj makineleri. URL: http://fimaks.coni/wp-content/uploads/2018/03/Silaj-Makiiieleri-K.pdf
13 Кормоуборочные комбайны. URL: https://rostselniash.coni/products/forage_harvesters
R (lnlk+гЩГ/.
R
ä"T*TT у
\ I \R Al /\
rlnrfff
a)
b)
Р и с. 4. Резание невысоких (а) и высоких растений, высота которых превышает расстояние от оси ротора до их основания (b): 1; 2 - траектории точек лезвий соседних ножей 1 и 2 на одной окружности; AjA2, B1B2, CE, DE, EG - отрезки растений после срезания ножом; Oxy - неподвижная ортогональная система декартовых координат с горизонтальной осью Ox и вертикальной осью Oy, направленной вверх; k- число ножей в плоскости вращения; h - минимальная высота стеблей после срезания, мм; H - высота стеблей до срезания, мм; R - радиус ротора (расстояние от оси ротора до лезвия ножа), мм; d - шаг по оси Ox, мм; s - подача или перемещение ротора вдоль оси Ox за время поворота на угол, равный радиальному углу, образуемому соседними ножами в плоскости вращения, мм; в - угол (параметр), рад; а - угловая скорость ротора, рад/с; v - скорость точки оси ротора вдоль оси Ox, мм/с; X - кинематический параметр (X = Ra/v); xj, yj - координаты x, y точки лезвия ножа 1; x2, y2 - координаты x, y точки лезвия ножа 2; в - угол (параметр); а - угол (параметр); x0, xe - дополнительные значения F i g. 4. Low plants cutting (a) and cutting of high plants, whose height is greater than the distance from
the axis of the rotor to the base of the plants (b): 1, 2 - trajectories of knife tips 1 and 2; AjA2, BjB2, CE, DE, EG - plant cuts after cutting by knife; Oxy - fixed orthogonal Cartesian coordinates with horizontal axis Ox and vertical upward axis Oy; k - number of knives in the plane of rotation; h - minimal height of the stems after cutting; H - height of the stems before cutting; R - radius of a rotor; d - step in the Ox direction; s - distance between neighboring trajectories of knives along an axis Ox; в - angle (parameter); a - angular velocity of the rotor; v - speed of the rotor axis point along an axis Ox; X - kinematic parameter (X = Ra/v); xl, yj - x, y coordinates of the point of the blade 1; x2, y2 - x, y coordinates of the point of the blade 2; в - angle (parameter); aj - angle (parameter); x0, xe - additional values
ротора в тот момент, когда нож 1 находится в низшем положении (рис. 4, а). Уравнения траектории 1 точки лезвия ножа 1 в параметрическом виде запишем следующим образом [Там же]:
= R [(а2 + 2п / k )/Я-y2 = R cosa2,
■sin а
(4)
|xj = r (а / я +sina ) [ y = R cos а,
(3)
где ax - параметр, обозначающий угол, показанный на рисунке 4, а (at > 0).
Уравнение траектории 2 точки лезвия ножа 2 записываются аналогичным образом:
где а2 - параметр, обозначающий угол, показанный на рисунке 4, а (а1 > - 2л/к).
Резание рассматривалось как множество циклов; за цикл принималось резание лезвием ножа 2 растений между траекториями 1 и 2.
Разделим фигуру между траекториями 1 и 2 на две части в пределах двух интервалов по оси Ох:
Х0 £ Х1 £ Хе,
Xg — %2 — Xg + s,
(5)
(6) 285
где х0 - координата нижней точки пересечения траекторий 1 и 2 по оси Ох; х. -координата верхней точки пересечения траектории 1 со стеблем по оси Ох.
Следующее равенство можно использовать для определения х0 (рис. 4, Ь):
x0 = R(в / Я + sin0).
(7)
Величина xe может быть найдена с помощью следующих формул (рис. 4, b):
_ÍR(в /1 + sinв),0 < H-h < R Х _{ R(0,5я-/1 +1),H-h > R, (8)
где cos P=(R + h - H)/R.
В пределах первого интервала (x0 < Xj, x2 < xe) длина u частицы равна длине отрезка AJA2 (рис. 4):
u = У 2 - У1 ,
отрезка EG перерезается следующими ножами в режущей паре «нож - проти-ворежущий брус» на отрезки длиной 5 и последний отрезок, длина которого не превышает 5.
Пусть N - число частей, равное числу резаний высокого растения от нижней точки С с координатами х2, у2 (х. < х2 < х. + 5); и. - длина 1 - й части (1 < 1 < N); LCE - длина отрезка СЕ, где
LcE = У2 + (х. +s - х2)2]1/2.
Найдем число N и длины всех частей высокого растения после резания его ножами с учетом принятых допущений.
Если Н - ^ + h) < LCE, то N =1 и и1 = Н- ^ + И) + У2; иначе:
1) и1 = LCE; к = 2;
2) ик = 5;
3) если
R
- J2 > H, то
(9)
где х1, у1, х2, у2 можно определить из уравнений (3-4) и принять х1 = х2.
В пределах второго интервала х. < х2 < х. + 5 длина и частицы зависит от высоты растений. Короткие растения (Н — h < R) не упираются в передний щиток и не сгибаются им перед срезанием. Длина и частиц коротких растений равна длине отрезка В1В2 (рис. 4, а): 1 2
и = у2 + Н - h - R, (10)
где у2 определяется из уравнений (4).
Примем, что:
1) высокие растения (Н > R + И) упираются в передний щиток и изгибаются перед щитком;
2) срезанное в нижней точке С высокое растение перерезается тем же ножом в режущей паре «нож - проти-ворежущий брус» в точке Е, если длина верхней части растения превышает длину отреза СЕ (рис. 4, Ь);
3) оставшаяся после резания в точке Е часть высокого растения в виде
N = к и щ = Н - ^ + И) + у2 - 1, = 1 и,; иначе к = к + 1. Далее следует переход к пункту 2.
Пусть I - длина верхней части прямого растения; М(1) - масса верхней части растения длиной I. Примем, что плотность р растения на единицу длины в данной точке (сечении) увеличивается прямо пропорционально расстоянию у от вершины растения до точки:
р = а y, М(1) = а I 2/ 2,
где а - коэффициент пропорциональности.
Тогда массу части растения длиной А, расположенной на расстоянии I от вершины растения, можно рассчитать по формуле:
т(1, А) = М(1 + А) — М(1)
или
m(l, A) =a[(l + А)2 - l 2]/2. (11) Согласно определению,
(12)
mL = 100 ML / MT ,
где Мт - сумма вычисляемых по формуле (11) без коэффициента а и при заданных параметрах I, А масс всех частиц; Мь - сумма вычисляемых по формуле (11) без коэффициента а масс тех частиц, длина А которых не превышает заданной величины Ь.
Согласно определению,
т = ШЬ \ь =19тт -т1 11=8шш. (13)
Алгоритм расчета массовой доли частиц в заданных диапазонах длины
Пусть 8иЬА(1, А, Ь, МЬ, Мт) - подпрограмма с входными переменными I, А, Ь и выходными переменными МЬ, Мт, выполняющая три шага.
1. Расчет т(1, А) по формуле (11) без коэффициента а.
2. Увеличение параметра Мт на т(1, А).
3. Увеличение параметра МЬ на т(1, А), если I < Ь.
Следующий алгоритм был реализован для вычисления величины тЬ.
1. Назначение входных переменных R, (Н - И), s, й, Ь, k.
2. Расчет величины X из уравнения (1).
3. Расчет величины 0 из уравнения (2).
4. Расчет величин х0, х, по формулам (7-8).
5. Назначение переменных I, МЬ, Мт: I = 1; МЬ = 0; Мт = 0.
6. Назначение величины х1: х1 = х0+I й.
7. Вычисление величин а1 и у1 из уравнений (3).
8. Назначение величины х2: х2 = х
9. Вычисление величин а2 и у2 из уравнений (4).
10. Вычисление длины и частицы из равенства (9) и выполнение подпрограммы БиЬА(1, А, Ь, М Мт) при I = Н - R - И + у1, А = и.
11. Увеличение переменной I на 1.
12. Выполнение пунктов 6-11, если
х1 < хе.
13. Назначение величины х2: х2 =
= х0 + I й.
14. Если Н - И < R, то:
1) вычисление длины и частицы из равенства (10);
2) выполнение подпрограммы
БиЬА(1, А, Ь, МЬ, Мт) при I = 0, А = и.
15. Если Н- И > R, то:
1) вычисление длины отрезка СЕ: Ьсе = [У22 + (х, +5 - х2)2]1/2;
2) если Н- R - И + у2< ЬСЕ, то:
а) и = Н - ^ + И) + У2;
б) выполнение подпрограммы БиЬА(1, А, Ь, МЬ, Мт) при I = 0, А = и1;
3) если Н- R - И + у2 = ЬСЕ, то:
а) и1 = ЬСЕ; иит = и1; выполнение подпрограммы 8иЪА(1, А, Ь, М, Мт) при I = Н - R - И + у2 - и1, А = и1;
б) если и5ит + 5 + R + И - у2 > Н, то и1 = Н - ^ + И) + у2 - иит; выполнение подпрограммы БиЬА(1, А, Ь, МЬ, Мт) при I = 0, А = и1;
в) если иит + 5 + R + И - у2 < Н, то
и1 = 5 иит = и*ит + 5 выполнение подпрограммы БиЬА(1, А, Ь, МЬ, Мт) при I = Н - R - И + у2 - иит, А = и1 и возвращение к предыдущему шагу 2.
16. Увеличение переменной I на 1.
17. Выполнение пунктов 13-16, если х2 < х, + 5.
18. Вычисление величины тЬ по формуле (12).
Математическая модель резания частиц растений в измельчающем аппарате двойного измельчения
Измельчение растений кормоубо-рочным комбайном двойного измельчения (с двойным измельчающим аппаратом) включает две стадии: измельчение роторным аппаратом на первой стадии и измельчение с выгрузкой растений дисковым измельчающе-швыряющим аппаратом на второй стадии. На выходе из роторного аппарата частицы захватываются шнеком и подаются в приемное окно дискового аппарата. Будем считать, что в шнеке частицы растений представляют собой прямые отрезки, движущиеся вдоль оси шнека с некоторой известной средней скоростью, позволяющей рассчитать среднюю расчетную длину резки.
Введем неподвижную систему декартовых координат с осью О1£, вдоль прямого лезвия противорежущего бруса и осью О1п вдоль оси шнека, а также дополнительные переменные: w - длина частицы в виде отрезка МN в момент резания ее ножом дискового аппарата, мм; тк - масса частицы в виде отрезка МN в момент резания ее ножом дискового аппарата, кг; Жс - ширина приемной горловины дискового измельчающе-швыряющего аппарата, мм; 1с - средняя расчетная длина резки, обеспечиваемая дисковым измельчающе-швыряющим аппаратом, мм.
Пусть в момент перерезания частица MN располагается в плоскости О1£ц, передняя точка М имеет координаты П, отрезок MN образует с осью О1£ угол Ф (рис. 5).
0 <^< Wc/2,
0 <п< ,
Pi ^р^P2,
где величины фх, ф2 зависят от 288
(14)
Pi =
(15)
0,если w < arccos( / w),если w > %,
п,еслиw <ЖС <р2 =\п-arccos[(ЖС-%)/^, (16) еслиw >ЖС
Число п частей, сформированных из частицы МN после всех перерезаний ее ножами дискового аппарата, может быть рассчитано по следующим формулам:
n = <
1,если w sin ф < п, [(wsin<р-п)/ lc] + 2, (17) если w sin ф > r¡.
Р и с. 5. Частица MN в приемном окне дискового аппарата в момент резания ножом F i g. 5. A particle MN at the cutting by a blade of a flywheel cutterhead
Примем, что переменные n, ф ~ случайные величины с равномерным законом распределения на отрезках [0, Wc /2], [0, lc] и [ф1, ф2] соответственно:
Пусть wi и mi - длина и масса части частицы MN с номером i (i = 1, 2,...n). Найдем выражения для длины и массы всех частей частицы, если величины w, mw, ф, n, lc известны, а массы частей прямо пропорциональны их длинам.
Если n = 1, частица MN не перерезается, и ее длина и масса не изменяются:
Wj = w,
mi = mw. (18)
Если n = 2, частица MN перерезается один раз, образуя две части. Длины w1, w2 и массы m1, m2 этих частей можно рассчитать по формулам:
w1 = n /sin ф,
w0 = w - w,
m2 = W2 mw /w
v2' " "j> mj = Wj mw /w,
(19)
Если n > 2, частица MN перерезается с образованием n частей. Длины w w .wn и массы этих частей можно рассчитать по формулам:
w1 = n / sin ф, w = w. = ...= w , = l / sin ф,
2 3 n - 1 c r'
w = w
n
w,
J...
m. = w. m
i = 1,
w
n
/w, n.
(20)
Случайные величины п, ф можно задавать с любой точностью следующими формулами:
iWc / (2Й|),
П = jlc 1
p = p + q (p2-p )
(21.1) (21.2) (21.3)
7. Назначение переменной д: q = 1.
8. Вычисление величины ф по формуле (21.3).
9. Вычисление величины п по формуле (17).
10. Вычисление величин w1 и т1 по формулам (18), если п = 1, или величин w1, w2 и т1, т2 по формулам (19), если п = 2, или величин w1, w2, ... wn и т„ т2,...тп по формулам (20), если п > 2.
11. Увеличение параметра MT на m1, если n = 1, или на m1 и m2, если n = 2, или на m1, m2,... mn, если n > 2.
12. Увеличение параметраML на m1, если w1 < L; увеличение параметра ML на m2, если w2 < L; увеличение параметра ML на mn, если wn < L.
13. Увеличение переменной q на 1.
14. Выполнение шагов 8-13, если Ф < Ф2.
15. Увеличение переменной j на 1.
16. Выполнение шагов 6-15, если П < l,
17. Увеличение переменной i на 1.
18. Выполнение шагов 3-17, если £ < Wc /2.
Алгоритм расчета величины mL для двойного измельчающего аппарата совпадает с алгоритмом расчета mL для роторного измельчающе-швыряющего аппарата, где в пункты 10, 14 и 15 внесены следующие изменения.
10. Вычисление длины u частицы из равенства (9); вычисление массы mw частицы по формуле (11) при a = 1, l = H - R - h + y1, A = u, выполнение подпрограммы SubB(w, mw, L, ML, MT) при w = u.
В пункте 14:
2) Вычисление массы mw частицы по формуле (11) при a = 1, l = 0, A = u, выполнение подпрограммы
SubB(w, mw, L, ML, MT) при w = u. В пункте 15:
2) б) Вычисление массы mw частицы по формуле (11) при a = 1, l = 0, A = u1, выполнение подпрограммы SubB(w, mw, L, Ml, MT) при w = u1.
14 Polyanin A. D., Manzhirov A. V. Handbook of mathematics for engineers and scientists. Taylor & Francis Group, LLC, 2007. 1509 p. URL: https://isidore.co/calibre/get/pdf/6338
где п, пп, пф - заданные натуральные
числа; I, у, д - целые числа: \ = 1, 2,...п^; у = 1 2,...п,; д = 1 2,...пф.
Так как п, ф - случайные величины с равномерным законом распределения, то целые числа I, у, д в равенствах (21) можно также рассматривать как случайные величины в пределах от 1 до п, от 1 до п от 1 до пф соответственно с вероятностью 1/п, 1/п 1/пф соответственно14. Величины ф1, ф2 в формулах (15-16) зависят от /. Число комбинаций, в которых целые числа I, у, д примут все возможные значения, равно произведению (п • пп • пф).
Алгоритм вычисления массовой доли частиц в заданных диапазонах длины
Пусть w, - входные переменные, обозначающие длину и массу частицы до перерезания ножами дискового аппарата. Представим алгоритм подпрограммы SubB(w, т^ L, Мь, МТ), предназначенной для изменения величин переменных М М
mw, L
1. Назначение переменных Wc,
w,
L, ML, Mt.
ф
2. Назначение констант п, п„, п и переменной /: / = 1.
3. Вычисление величины £ по формуле (21.1).
4. Вычисление величин ф1, ф2 по формулам (15-16).
5. Назначение переменной у: у = 1.
6. Вычисление величины п по формуле (21.2).
3) Если Н - R - И + у2 = LCE, то:
а) и1 = ишт = и1; вычисление
массы тн, частицы по формуле (11) при
а = 1, I = Н - R - И + у2 — и1, А = и1, выполнение подпрограммы SubB(w, тц1, L, ML, МТ) при w = и1;
б) если иишп + 5 + R + И - у2 = Н, то и1 = Н - ^ + И) + у2 - и5ит, вычисление массы тм! частицы по формуле (11) при а =1, I = 0, А = и1, выполнение подпрограммы SubB(w, т^, L, М МТ) при
w = uv
в) если usum + 5 + R + h - y2 < H, то
U1 = 5 Usum = Usum + 5 вычисление массы
mw частицы по формуле (11) при a = 1, l =H - R - h + y2 - usum, Л = u1, выполнение подпрограммы SubB(w, m L, M MT) при w = u1 и возвращение к предыдущему шагу 2.
Расчет массовой доли частиц в заданных диапазонах длины
Алгоритм расчета массовой доли частиц растений в заданных диапазонах длины был реализован в среде программирования Lazarus. Исходные данные назначаемых переменных и констант приведены в табл. 2.
Расчеты показали, что распределение по длине массовой доли mL частиц растений, измельченных роторным аппаратом, имеет три участка, на которых функция распределения может быть аппроксимирована линейными зависимостями от аргумента L (рис. 6, а). На первом участке аргумент изменяется в диапазоне от нуля до величины, не превышающей подачи s на нож, на втором - равен подаче s с незначительным варьированием нижней и верхней границы диапазона, на третьем - в диапазоне от величины подачи до максимальной длины частиц. Испытания
роторного кормоуборочного комбайна типа КИР-1,5 на кошении с измельчением люцерны высотой 440 мм15 показали, что, несмотря на значительный диапазон, экспериментальные значения массовой доли частиц длиной до 50 мм и свыше 100 мм практически точно соответствуют теоретически найденным значениям, показанным на графике 1 функции распределения (рис. 6, а). Максимальное согласование между экспериментальными и расчетными данными массовой доли частиц длиной свыше 100 мм наблюдается при кошении с измельчением вики с овсом16 высотой 780 мм (график 3) и эспарцета17 высотой 1 200 мм (график 4).
Двойной измельчающий аппарат повышает качество измельчения в том случае, если частота вращения шнека и/или дискового или цилиндрического аппарата обеспечивают среднюю расчетную длину резки, не превышающую величины подачи роторного аппарата (рис. 6, Ь). При этом массовая доля частиц, длина которых превышает 100 мм, снижается практически до нуля, а наибольшая доля частиц находится в диапазоне увеличенной и уменьшенной на 5 мм расчетной длины резки, изменяясь от 45 % при расчетной длине резки 25 мм до 75 % при расчетной длине резки 5 мм. Таким образом, невозможность регулирования угловой скорости шнека и/или дискового барабана не позволяет использовать двойной измельчающий аппарат на заготовке кормов с различными требованиями к качеству измельчения.
Расчеты показали, что роторные кормоуборочные комбайны не позволяют обеспечить качество измельчения
15 Косилка-измельчитель роторная КИР-1,5М. Протокол испытаний № 11-11-13 (2130042). ФГБУ «Северо-Кавказская МИС». Зерноград. 2013. 6 с.; Косилка-измельчитель роторная КИР-1,85 «Корммаш». Отчет № 11-25-14 (1130122) базовых испытаний. ФГБУ «Северо-Кавказская МИС». Зерноград. 2014. 12 с.
16 Косилка-измельчитель роторная КИР-1,5М. Протокол испытаний. ФГБУ «Подольская государственная зональная МИС». 2015. URL: http://www.podolskmis.ru/2013-12-20-07-09-44/2013-12-27-06-33-09/53-kratkie-protokoly-ispytanij-za-2015-god/250-09-13-15-5130062
17 Косилка-измельчитель роторная КИР-1,5Е. Протокол испытаний № 14-42-2001 (4130232). ГУ «Центрально-Черноземная МИС». 2001. 4 с.
Т а б л и ц а 2 T a b l e 2
Исходные данные Initial data
Параметр / Parameter Значение / Value
R, мм / R, mm 285
R m 10-3, м/с / R m 10-3, m/s 40
d, мм / d, mm 0,5
k 1
W, мм / W, mm 246
nv n 30
£ ч
100 90 80 70 60 50 40 30 20 10
0
-
4 У/
■
А/ 1 // Г
JsV
90 80 70 60 50 40 30 20 10
20 40 60 80100 140 180 L, MM/i,mm
a)
220
20 30 40 50 60 70 80 90 L, мм/ L, mm b)
Р и с. 6. Графики зависимости от L массовой доли mLчастиц растений, длина которых не превышает L, после измельчения роторным кормоуборочным комбайном (а) и комбайном двойного измельчения (b): 1 - H - h = 367 мм, s = 71 мм; • - экспериментальные данные измельчения люцерны; 2 - H - h = 377 мм, s = 80 мм; ■ - экспериментальные данные измельчения
люцерны; 3 - H - h = 660 мм, s = 76 мм; ▲ - экспериментальные данные измельчения вики с овсом; 4 - H - h = 1 001 мм, s = 63 мм; ▼ - экспериментальные данные измельчения эспарцета; 5 - H - h = 660 мм, s = 76 мм, lc = 5 мм; 6 - H - h = 660 мм, s = 76 мм, lc = 10 мм; 7 - H - h = 660 мм, s = 76 мм, l = 25 мм
3 3 c
F i g. 6. The effect of L on mass fraction mL of plant particles, the length of which does not exceed L, after chopping by a rotary forage harveLter (a) and a double chop forage harvester (b): 1 - (H - h) = 367 mm, s = 71 mm; • - experimental data for alfalfa; 2 - (H - h) = 377 mm, s = 80 mm; ■ - experimental data for alfalfa; 3 - (H - h) = 660 mm, s = 76 mm; ▲ - experimental data for alfalfa with oats; 4 - (H - h) = 1001 mm s = 63 mm; ▼ - experimental data for sainfoin; 5 - (H - h) = 660 mm, s = 76 mm, lc = 5 mm; 6 - (H - h) = 660 mm, s = 76 mm, lc = 10 mm;
7 - (H - h) =c 660 mm, s = 76 mm, lc = 25 mm c
трав на сенаж, при котором массовая доля частиц длиной 8-18 мм находится в диапазоне 45-75 % (рис. 7, а). При подаче более 20 мм или при скорости комбайна более 1,6 м/с роторные комбайны обеспечивают массовую долю
таких частиц не более 35 %. Кормо-уборочные комбайны с двойным измельчающим аппаратом позволяют обеспечить заготовку трав на сенаж только при заданной и регулируемой частоте вращения барабана (дисково-
го или цилиндрического) либо шнека, подающего растения в камеру дискового или цилиндрического аппарата.
Обсуждение и заключение
1. Кормоуборочные комбайны с роторным измельчающе-швыряющим аппаратом, осуществляющим одновременно три операции (кошение растений на корню, измельчение растений, выгрузку их в транспортное средство), отличаются простотой конструкции и не позволяют обеспечить заготовку трав на сенаж с массовой долей частиц длиной 8-19 мм в пределах 45-75 %.
2. Разработанная математическая модель измельчения растений роторным измельчающе-швыряющим аппаратом позволяет оценить массовую долю частиц измельченных растений в любом диапазоне их длины в зависимости от высоты растений, высоты среза, подачи. В рамках данной модели в пределах высоты растений 400-1 200 мм и подачи 20-100 мм функция распределения массовой доли частиц по длине представляет собой три участка, на которых линии графика функции близки к прямым. С увеличением высоты растений
возрастает массовая доля частиц, длина которых примерно равна подаче.
3. Кормоуборочные комбайны с двойным измельчающим аппаратом позволяют обеспечить заготовку трав на сенаж с массовой долей частиц длиной 8-19 мм в пределах 45-75 % только при регулировании и соответствующей настройке угловой скорости шнека и/или дискового или цилиндрического барабана.
4. Разработанная математическая модель измельчения растений двойным измельчающим аппаратом позволяет оценить массовую долю частиц измельченных растений в любом диапазоне их длины в зависимости от высоты растений, высоты среза, подачи на нож, средней расчетной длины резки. Закон распределения массовой доли частиц измельченных растений зависит от средней расчетной длины резки. В пределах диапазона длины частиц, границы которого отличаются от расчетной длины резки на 5 мм, содержится 45-75 % частиц. Массовая доля частиц длиной 8-19 мм практически не зависит от подачи и высоты растений при расчетной длине резки менее 10 мм.
2015
20 30 40 50 60 70 80 90 20 30 40 50 60 70 80 90
мм / s, mm s, мм / 5, mm
а) Ъ)
Р и с. 7. Влияние подачи s на массовую долю частиц m длиной 8-19 мм при измельчении растений роторным кормоуборочным комбайном (а) и кормоуборочным комбайном двойного
измельчения (b): 1 - (H - h) = 250 мм; 2 - (H - h) = 600 мм; 3 - (H - h) = 1 000 мм; 4 - (H - h) = 660 мм, lc = 5 мм; 5 - (H - h) = 660 мм, lc = 10 мм; 6 - (H - h) = 660 мм, lc = 25 мм
F i g. 7. The effect of s on mass fraction of plant particles m with lengths from 8 mm to 19 mm for a rotary forage harvester (a) and a double chop forage harvester (b): 1 - (H - h) = 250 mm; 2 - (H - h) = 600 mm; 3 - (H - h) = 1000 mm; 4 - (H - h) = 660 mm, l = 5 mm; 5 - (H - h) = 660 mm, lc = 10 mm; 6 - (H - h) = 660 mm, lc = 25 mm
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Crop processing and chop length of corn silage: effects on intake, digestion, and milk production by dairy cows / M. A. Bal [et al.] // Journal of Dairy Science. 2000. Vol. 83, issue 6. P. 1264-1273. DOI: https://doi.org/10.3168/jds.S0022-0302(00)74993-9
2. Effects of the chop lengths of alfalfa silage and oat silage on feed intake, milk production, feeding behavior, and rumen fermentation of dairy cows / S. K. Bhandary [et al.] // Journal of Dairy Science. 2008. Vol. 91, issue 5. P. 1942-1958. DOI: https://doi.org/10.3168/jds.2007-0358
3. Effect of grass silage chop length when fed alone, or with corn silage, on digestion and metabolism in dairy cows / U. Tayyab [et al.] // Journal of Animal Science. 2018. Vol. 96, issue 3. P. 394. DOI: https:// doi.org/10.1093/jas/sky404.864
4. Addah W., Baah J., McAllister T. A. Effect of silage chop length on feed intake and feeding behaviour of finishing feedlot steers // Acta Agriculturae Scandinavica, Section A - Animal Science. 2016. Vol. 66, no. 2. P. 106-114. DOI: https://doi.org/10.1080/09064702.2016.1267789.
5. Saqib G. S., Finner M. F. Simulated ideal length of cut for forage harvesters // Transactions of the ASAE. 1982. Vol. 25, issue 5. P. 1237-1238. DOI: https://doi.org/10.13031/2013.33704
6. Morgan D. D. V., Osman M. M. A mathematical model of an apparatus for the assessment of the length distribution of chopped forage // Journal of Agricultural Engineering Research. 1984. Vol. 30. P. 157-164. DOI: https://doi.org/10.1016/S0021-8634(84)80017-7
7. O'Dogherty M. J. Chop length distributions from forage harvesters and a simulation model of chopping // Journal of Agricultural Engineering Research. 1984. Vol. 30. P. 165-173. DOI: https://doi. org/10.1016/S0021-8634(84)80018-9
8. Szendro P. Examination of chopping process in self-propelled forage harvesters // Acta Agronomica Hungarica. 1979. Vol. 28, issues 1-2. P. 106-119.
9. Yukueda M, Kawamura N. Studies of a direct-throw flail-type forage harvester (I). Journal of the Japanese society of agricultural machinery. 1976. Vol. 37, issue 4. P. 600-605. DOI: https://doi.org/10.11357/ jsam1937.37.600
10. Howe S. D., Bishop T. Forage harvesters: designs for the 80s // Power Farming. 1982. Vol. 64, issue 4. P. 36-45.
11. Белов М. И., Чепурной А. И. Перспективные измельчающие аппараты кормоуборочных комбайнов // Тракторы и сельхозмашины. 1985. № 5. С. 39-41.
12. Белов М. И., Чепурной А. И. Анализ схем кормоуборочных машин // Тракторы и сельхозмашины. 1986. № 6. С. 26-28.
13. Ерохин М. Н., Белов М. И., Судник Ю. А. Модель и экспериментальное исследование ротационного режущего аппарата // Тракторы и сельхозмашины. 2003. № 12. С. 21-24.
Поступила 08.12.2018; принята к публикации 11.02.2019; опубликована онлайн 28.06.2019
Об авторе:
Белов Михаил Иванович, профессор, кафедра сопротивления материалов и деталей машин, Институт механики и энергетики, ФГБОУ ВО «Российский государственный аграрный университет - МСХА имени К. А. Тимирязева» (127550, Россия, г. Москва, ул. Тимирязевская, 49), доктор технических наук, профессор, ReseacherID: T-5622-2018, ORCID: https://orcid.org/0000-0001-9907-8825, B-Mikhael@yandex.ru
Автор прочитал и одобрил окончательный вариант рукописи.
REFERENCES
1. Bal M.A., Shaver R.D., Jirovec A.G., Shinners K.J., Coors J.G. Crop processing and chop length of corn silage: Effects on intake, digestion, and milk production by dairy cows. Journal of Diary Science. 2000; 83(6):1264-1273. DOI: https://doi.org/10.3168/jds.S0022-0302(00)74993-9
2. Bhandary S.K., Ominsky K.M, Wittenberg K.M., Plaizier J.C. Effects of the chop lengths of alfalfa silage and oat silage on feed intake, milk production, feeding behavior, and rumen fermentation of dairy cows. Journal of Dairy Science. 2008; 91(5):1942-1958. DOI: https://doi.org/10.3168/jds.2007-0358
3. Tayyab U., Wilkinson R., Reynolds C., Sinclair L. Effect of grass silage chop length when fed alone, or with corn silage, on digestion and metabolism in dairy cows. Journal of Animal Science. 2018; 96(3):394. DOI: https://doi.org/10.1093/jas/sky404.864
4. Addah W., Baah J., McAllister T.A. Effect of silage chop length on feed intake and feeding behaviour of finishing feedlot steers. Acta Agriculturae Scandinavica, Section A - Animal Science. 2016; 66(2):106-114. DOI: https://doi.org/10.1080/09064702.2016.1267789
5. Saqib G.S., Finner M.F. Simulated ideal length of cut for forage harvesters. Transactions of the ASAE. 1982; 25(5):1237-1238. DOI: https://doi.org/10.13031/2013.33704
6. Morgan D.D.V., Osman M.M. A mathematical model of an apparatus for the assessment of the length distribution of chopped forage. Journal of Agricultural Engineering Research. 1984; 30:157-164. DOI: https://doi.org/10.1016/S0021-8634(84)80017-7
7. O'Dogherty M.J. Chop length distributions from forage harvesters and a simulation model of chopping. Journal of Agricultural Engineering Research. 1984; 30:165-173. DOI: https://doi.org/10.1016/ S0021-8634(84)80018-9
8. Szendro P. Examination of chopping process in selfpropelled forage harvesters. Acta Agronomy Academy Scientific Hungary. 1979; 28(1-2):106-119.
9. Yukueda M., Kawamura N. Studies of a direct-throw flail-type forage harvester (I). Journal of the Japanese Society of Agricultural Machinery. 1976; 37(4):600-605. DOI: https://doi.org/10.11357/ jsam1937.37.600
10. Howe S.D., Bishop T. Forage harvesters: Designs for the '80s. Power Farming. 1982; 64(4):36-45.
11. Belov M.I., Chepurnoy A.I. [Perspective chopper units of forage harvesters]. Traktory i selkhoz-mashiny = Tractors and Agricultural Machinery. 1985; 5:39-41. (In Russ.)
12. Belov M. I., Chepurnoy A. I. [Analysis of forage harvester schemes]. Traktory i selkhozmashiny = Tractors and Agricultural Machinery. 1986; 6:26-28. (In Russ.)
13. Yerokhin M. N., Belov M. I., Sudnik Yu. A. [Model and experimental study of rotary mower]. Traktory i selkhozmashiny = Tractors and Agricultural Machinery. 2003; 12:21-24. (In Russ.)
Received 08.12.2018; revised 11.02.2019;published online 28.06.2019
About author:
Mikhail I. Belov, Professor, Chair of Strength of Materials and Machine Parts, Institute of Mechanics and Energetics, Russian Timiryazev State Agrarian University (49 Timiryazevskaya St., Moscow 127550, Russia), D.Sc. (Engineering), ReseacherID: T-5622-2018, ORCID: https://orcid.org/0000-0001-9907-8825, B-Mikhael@yandex.ru
The author has read and approved the final version of the paper.
