Научная статья на тему 'ОБОСНОВАНИЕ ПАРАМЕТРОВ УСТАНОВКИ ЗАЖИМНОГО ТРАНСПОРТЕРА И ДЕКИ ОБМОЛАЧИВАЮЩЕГО УСТРОЙСТВА ЛИНИИ ПЕРВИЧНОЙ ПЕРЕРАБОТКИ'

ОБОСНОВАНИЕ ПАРАМЕТРОВ УСТАНОВКИ ЗАЖИМНОГО ТРАНСПОРТЕРА И ДЕКИ ОБМОЛАЧИВАЮЩЕГО УСТРОЙСТВА ЛИНИИ ПЕРВИЧНОЙ ПЕРЕРАБОТКИ Текст научной статьи по специальности «Механика и машиностроение»

CC BY
29
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЛЕН / ОБМОЛАЧИВАЮЩЕЕ УСТРОЙСТВО / ЭЛАСТИЧНЫЙ РАБОЧИЙ ОРГАН / ЛЕНТА СТЕБЛЕЙ ЛЬНА / ЗОНА ОБМОЛОТА / ЗАЖИМНОЙ ТРАНСПОРТЕР / СЕПАРИРУЮЩАЯ РЕШЕТКА

Аннотация научной статьи по механике и машиностроению, автор научной работы — Шаршунов В. А., Курзенков С. В., Левчук В. А., Цайц М. В.

В разработанном в УО БГСХА обмолачивающем устройстве с эластичным рабочим органом линии первичной переработки льна технологический процесс обмолота определяют конструктивное исполнение и взаимное расположение основных рабочих органов таких, как зажимной транспортер, сепарирующая решетка и эластичный бич. Обмолот лент стеблей льна происходит в пространстве между бичом и декой, при этом сама лента стеблей льна удерживается зажимным транспортером. Смещение зажимного транспортера относительно траектории движения бича приводит к изменению как зоны обмолота, так и «мертвой зоны». Кроме того, воздействие бича на обрабатываемую ленту льна может привести к ее изгибу на величину больше предельного угла, что приведет к повреждению самих стеблей, а как следствие к потере урожая. Получены аналитические зависимости позволяющие определить параметры взаимного расположения рычажного механизма и зажимного транспортера разрабатываемого обмолачивающего устройства с учетом их конструктивных параметров и технологических параметров процесса обмолота. Определена геометрическая форма сепарирующей решетки и параметры ее установки. Установлена связь параметров установки рабочих органов обмолачивающего устройства с технологическими параметрами процесса обмолота. Приведен пример расчета в математическом пакете MathCad, позволяющий определить параметры взаимного расположения зажимного транспортера, сепарирующей решетки и рычажного механизма предлагаемого обмолачивающего устройства при заданных параметрах рычажного механизма. Аналитические зависимости, рассмотренные в данной статье, войдут в основу методики проектирования разрабатываемого обмолачивающего устройства с эластичным рабочим органом линии первичной переработки льна и моделирования процесса обмолота.In the threshing device with an elastic working body of the flax primary processing line developed at the BE BSHA, the technological process of threshing is determined by the design and relative position of the main working bodies of such a cassette conveyor, a separating grid and an elastic scourge. The threshing of the flax stalk strips takes place in the space between the whip and the deck, while the flax stalk belt itself is held by the clamping conveyor. The displacement of the clamping conveyor relative to the trajectory of the movement of the whip leads to a change in both the threshing zone and the «dead zone». In addition, the impact of the whip on the processed flax tape can lead to its bending by an amount greater than the limiting angle, which will damage the stems themselves, and as a result, yield loss. Analytical dependences are obtained that allow determining the parameters of the relative position of the lever mechanism and the clamping conveyor of the developed threshing device, taking into account the design and technological parameters of the lever mechanism. The geometric shape of the separating grating and the parameters of its installation are determined. The relationship between the installation parameters of the working bodies of the threshing device and the technological parameters of the threshing process has been established. An example of a calculation in the mathematical package MathCad is given, which makes it possible to determine the parameters of the relative position of the clamping conveyor, the separating grid and the lever mechanism of the proposed threshing device for given parameters of the lever mechanism. The analytical dependencies considered in this article will form the basis of the design methodology for the developed threshing device with an elastic working body of the flax primary processing line and modeling the threshing process.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Шаршунов В. А., Курзенков С. В., Левчук В. А., Цайц М. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «ОБОСНОВАНИЕ ПАРАМЕТРОВ УСТАНОВКИ ЗАЖИМНОГО ТРАНСПОРТЕРА И ДЕКИ ОБМОЛАЧИВАЮЩЕГО УСТРОЙСТВА ЛИНИИ ПЕРВИЧНОЙ ПЕРЕРАБОТКИ»

УДК 631.36

ОБОСНОВАНИЕ ПАРАМЕТРОВ УСТАНОВКИ ЗАЖИМНОГО ТРАНСПОРТЕРА И ДЕКИ ОБМОЛАЧИВАЮЩЕГО УСТРОЙСТВА ЛИНИИ ПЕРВИЧНОЙ ПЕРЕРАБОТКИ

В. А. ШАРШУНОВ

УО «Белорусский государственный университет пищевых и химических технологий», г. Могилев, Республика Беларусь, 212027, e-mail: [email protected]

С. В. КУРЗЕНКОВ, В. А. ЛЕВЧУК, М. В. ЦАЙЦ

УО «Белорусская государственная сельскохозяйственная академия», г. Горки, Республика Беларусь, 213407, e-mail: [email protected]

(Поступила в редакцию 15.03.2022)

В разработанном в УО БГСХА обмолачивающем устройстве с эластичным рабочим органом линии первичной переработки льна технологический процесс обмолота определяют конструктивное исполнение и взаимное расположение основных рабочих органов таких, как зажимной транспортер, сепарирующая решетка и эластичный бич. Обмолот лент стеблей льна происходит в пространстве между бичом и декой, при этом сама лента стеблей льна удерживается зажимным транспортером. Смещение зажимного транспортера относительно траектории движения бича приводит к изменению как зоны обмолота, так и «мертвой зоны». Кроме того, воздействие бича на обрабатываемую ленту льна может привести к ее изгибу на величину больше предельного угла, что приведет к повреждению самих стеблей, а как следствие к потере урожая.

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

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

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

In the threshing device with an elastic working body of the flax primary processing line developed at the BE BSHA, the technological process of threshing is determined by the design and relative position of the main working bodies of such a cassette conveyor, a separating grid and an elastic scourge. The threshing of the flax stalk strips takes place in the space between the whip and the deck, while the flax stalk belt itself is held by the clamping conveyor. The displacement of the clamping conveyor relative to the trajectory of the movement of the whip leads to a change in both the threshing zone and the «dead zone». In addition, the impact of the whip on the processed flax tape can lead to its bending by an amount greater than the limiting angle, which will damage the stems themselves, and as a result, yield loss.

Analytical dependences are obtained that allow determining the parameters of the relative position of the lever mechanism and the clamping conveyor of the developed threshing device, taking into account the design and technological parameters of the lever mechanism. The geometric shape of the separating grating and the parameters of its installation are determined. The relationship between the installation parameters of the working bodies of the threshing device and the technological parameters of the threshing process has been established. An example of a calculation in the mathematical package MathCad is given, which makes it possible to determine the parameters of the relative position of the clamping conveyor, the separating grid and the lever mechanism of the proposed threshing device for given parameters of the lever mechanism.

The analytical dependencies considered in this article will form the basis of the design methodology for the developed threshing device with an elastic working body of the flax primary processing line and modeling the threshing process.

Key words: flax, threshing device, elastic working body, line of flax stalks, threshing zone, clamping conveyor, separating grid.

Введение

Любой технологический процесс определяется его конструкторским исполнением и взаимным расположением рабочих органов. Это утверждение справедливо можно отнести и к обмолачивающему устройству лент льна линии первичной переработки, предложенному сотрудниками УО БГСХА. Основными узлами данного обмолачивающего устройства являются: зажимной транспортер, сепарирующая решетка и эластичный бич. Их состав и принцип работы подробно рассмотрены в источниках [1-3].

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

ширине. Ее неоднородность зависит от длины стеблей льна. Также установлено [5], что участок расположения семенных коробочек льна в ленте по ее ширине изменяется в широком диапазоне.

Основная часть

Лента льна, подаваемая на обмолот, представляет собой слой стеблей, длина которых находится в пределах 0,7...1,2 м. Стебли льна в ленте смещены друг относительно друга, это смещение характеризуется растянутостью ленты льна [6]. Кроме того, в ленте имеется относительный перекос стеблей.

Технологическим процессом предусматривается зажатие обмолачиваемой ленты льна зажимным транспортером на расстоянии Ькс = 0,2...0,25 м от комлевой части стеблей (рис. 1), часть ленты непосредственно удерживается транспортером Ьзт = 0,15 м, а это значит, что изменению (в зависимости от урожая и растянутости стеблей в ленте) подвержена часть ленты льна Ьоб, в пределах которой расположена зона семенных коробочек Ьск величина которой изменяется от 0,2 до 0,35 м (по некоторым источником она может составлять 0,2... 0,5 м) [4, 7]. Чем больше зона расположения семенных коробочек в ленте Ьск , тем меньше остается расстояние от зажимного транспортера до нее. А это значит, что зона воздействия бича на ленту должна быть максимально приближена к зажимному транспортеру.

Рис. 1. Общая схема устройства: 1 - кривошип; 2 - шатун; 3 - коромысло, 4 - стойка; 5 - зажимной транспортер: 6 - лента стеблей льна; 7 - сепарирующая решетка (дека)

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

Однако в таком случае возникает опасность излома стеблей льна. Поскольку бич, двигаясь по своей криволинейной траектории, преодолевая сопротивление изгибу стеблей льна, изгибает их в точке

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

Анализируя схему (рис. 1), можно заметить, что при таком расположении бич в точке Б движется вертикально вниз, а угол п стремится к 90°, что существенно превышает предельное значение угла изгиба стебля Ппред = 30...43° (зависит от диаметра стебля) [8]. Решить проблему изгиба стеблей льна можно, подняв траекторию движения бича относительно линии стеблей льна, зажатых в транспортере, или отнеся ее вправо относительно зажимного транспортера. В обоих перечисленных случаях (рис. 2) возникает участок ленты льна, не подверженный воздействию рабочего органа Ьмз «мертвая зона». С точки зрения минимизации величины «мертвой зоны» Ьмз, перспективнее видится сдвиг траектории движения бича вверх (рис. 2). Тем самым при минимальной «мертвой зоне» стеблей можно существенно уменьшить их угол изгиба. При этом максимальное значение угла изгиба стебля льна будет соответствовать точке Р, в которой линия стеблей льна будет располагаться по касательной к траектории движения бича.

Рис. 2. Схема к определению ординаты точки Б и величины «мертвой зоны» - БК На основании изложенного можно сделать вывод о том, что располагать траекторию движения бича необходимо выше линии стеблей льна, зажатых транспортером, таким образом, чтобы угол изгиба п не превышал своего предельного значения (ппред = 30.43°). В тоже время величина «мертвой зоны» не должна превышать 0,15 м (значение величины Ьмз = 0,08.0,15 м) [5].

Установка зажимного транспортера по отношению к механизму обмолачивающего устройства (рис. 2) характеризуется координатами точки Б(хэ;уэ), где хэ иуэ - соответствующие проекции точки Б на декартовую систему координат, центр которой находится в точке О, являющийся центром вращения кривошипа. Заметим, что положение точки Б непосредственно связано с траекторией движения точки С (крайней точки положения шатуна) (рис. 1), которая описывается системой [9].

(1 +1 )• ХА (()- 1 ■ ХВ (<)

С:

хс (() = -.Ус(() =

I

(1 +1 )• УА (<)-1 • Ув (<)

(1)

I

где I - длина участка шатуна АВ, м; ¡1 - длина участка шатуна АС, м; ^ - время, с. Выделим на этой траектории, представляющей эллипс, характеризующие ее точки 7\, Т2, Тз и Т4, являющиеся вершинами эллипса, а также определим координаты центра эллипса О2.

Угол поворота кривошипа ф рад, в зависимости от времени, определяется выражением [9]:

ф = Ш, (2)

где ю - угловая скорость вращения кривошипа ОА, рад/с.

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

2%

га 172

Тогда вершины траектории движения точки С определяются из (1) как Т:(х«(0); ус(0)), Т2(х<с(Д /4); уе(1к / 4)), Тз(хе(^к / 2); у<(Ъ / 2)), Т4(хс(3^к / 4); ус(3^к / 4)), а центр эллипса будет определятся координатами 02(О;(ус(^к / 4) + Ус(3^ / 4)) / 2) (рис. 2).

С учетом вышесказанного траектория движения точки С с привязкой к системе координат хОу может быть описана уравнением:

х2 (У- Уо2) . ...

— + ,2 2 = 1 , (3) « Ьэ

где аэ = хс(0) - большая полуось эллипса, м;

Ч. Л (з • к

„ УС I 4 J-* ч 4 у

Ь =--—-—-—---— меньшая полуось эллипса, м;

к Л (з • I

Ус [ 4 J + Ус ^ 4 у0г =--—-—-—---— ордината центра эллипса.

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

Уо ^ У02 . (4)

Возможные значения ординаты точки В и величину «мертвой зоны» - ВК определим из схемы (рис. 2).

Примем во внимание условие, что угол изгиба п стеблей в ленте зажимного транспортера при воздействии на них эластичного рабочего органа не должен превышать 30°-43°.

Из рис. 2 видно, что прямая ВР является касательной к траектории движения точки С на участке Т3Т1. Это означает, что угловой коэффициент прямой ВР определяется, как

кт = 1§Сл), (5)

а с другой стороны

Ш(Л) = —. (6)

хС С)

Решив уравнение (6) относительно времени ^ получим момент времени ?р, при котором точка С траектории движения перейдет в точку контакта Р со стеблями. Тогда координаты точки Р можно определить выражениями:

хр = хс ((р), Ур = Ус {Ч). (7)

Составим уравнение прямой ВР:

х - Хо У - Уо

хр ХВ УР У0

Чтобы определить угловой коэффициент этой прямой выразим из (8) переменную у

.Ур -Уо

(8)

У = •{х-хв) + Ув , (9)

где = Ур-Уо = .

Тогда ордината точки В определяется выражением:

Уо = ^(п) •{ хР - хо) + Ур. (10)

Для определения «мертвой зоны» стеблей, т. е. длины отрезка ВК, определим момент времени ^вк при котором точка С переходит в точку К из уравнения:

Ус {о-Уо = 0. (11)

Тогда общая точка К будет равна:

хк = хс {^ок ) , (12) а величина «мертвой зоны» стеблей будет определятся формулой:

А» = ХК - ХП • (13)

Так как эластичный рабочий орган взаимодействует со стеблями на промежутке времени ^эк^к],

то можем определить длину этого участка КТ\ по следующей зависимости:

I-

4з = ХС(Г ))2 +())2 л •

(14)

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

Продемонстрируем пример расчета параметров установки зажимного транспортера и деки в математическом пакете МаШСа^ при фиксированных значениях длины кривошипа (Я = 0,3 м), длины шатуна (I + 1\= 0,808 + 0,4 = 1,208 м), длины коромысла (Я; = 0,95 м), горизонтальное смещение точки крепления коромысла Я; относительно центра вращения кривошипа (а = 0,95 м), вертикальное смещение точки крепления коромысла Я; относительно центра вращения кривошипа (Ь =0,75 м). Предположим, что частота вращения кривошипа ОА п = 3 об/с, а предельное значение угла изгиба стеблей льна

"Ппред 30 .

Вариант оформления исходных данных и расчетной части в математическом пакете MathCad показан на рис. 2.

¡¡¡tail

■ [Обоснование параметров расстановки зажимного транспортера]

jyj зил Правка Вил

□ -ö н

"Л да =

п

100%

fi Xf xfy Xfy ¿ ¿ ™ г I П J H дта JiSV Й- V

Normal

v Arial

V 10 V В JT U Ш Ш Ш

Расчет величины мертвой зоны и длины участка стеблей взаимодействующего с бичом

частота вращения кривошипа ОА п:=3 предельный угол изгиба стеблей льна ц := 30

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

¡aj:= п 2 тт = IS.Sí

ЛЛЛ/

tk:= — = 0.333

uJ

t:= 0,0.002.. tk

t:= ■

tk

-tan

T)-Tt

ISO

-0.577

угол изгиба стеблей

Найдем момент времени tP

-yC(t) , dt i il тт i(t) :=- + tan -—

d .... USO

-xC(t) dt

xP := xC(tP) = -0.274

ÏD := xC| y I = -0.449

tP := root(j(t) ,t) = 0.217 момент времени tP

Найдем координаты точки Р

уР := yC(tP) = -0.631 Найдем координаты точки D

yD := tanl | ■ f хР - хс' ^

координаты точки Р + уР = -0.53 координаты точки D

ISO ) \ \2 Найдем координаты точки К tDK:= root(yC(t) - yD,t) = 0.195 момент времени tDK

Ж. := xC(tDfg = -0.391 уК := yC(tDK) = -0.53 координаты точки К

Определим величину "мертвой зоны" стеблей Ьмз := лК - xD = 0.0575 величина "мертвой зоны" стеблей Определим длину участка стеблей, взаимодействующих с бичем ■tk

1бз :

tDK

-xC(t) I + f-yC(t) I dt =1.049

dt ! I dt !

Рис. 3. Исходные данные и расчетная часть определения величины «мертвой зоны» и длины участка стеблей взаимодействующего с бичом

Полученное в результате расчета значение минимальной величины «мертвой зоны» Lмз = 0,058 м не учитывает особенностей исполнения зажимного транспортера, и в целях обеспечения свободного движения бича, может быть увеличено на величину технологического зазора. Расчетное значение максимальной величины длины участка стеблей, взаимодействующего с бичом, составит 1,049 м, что существенно превысит участок ленты стеблей находящегося в зоне обмолота LCT и обеспечит эффективную работу обмолачивающего устройства.

Заключение

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

По результатам проведенных расчетов установлено, что условию минимального значения величины «мертвой зоны» Lмз = 0,058 м и предельного значения угла изгиба стеблей льна Ппред = 30° удовлетворяют:

- установка зажимного транспортера (точка D) вниз относительно оси вращения кривошипа ОА на величину yD = 0,53 м, с горизонтальным смещением влево на величину xd = 0,449 м;

- установка ближней к зажимному транспортеру части сепарирующей решетки (деки) вниз относительно оси вращения кривошипа ОА на величину ук + А = 0,53 + А м, а горизонтальное смещение влево - хк = 0,391 м без учета технологического зазора предусматривающего свободное движение бича.

ЛИТЕРАТУРА

1. Шаршунов, В. А. Исследование обмолачивающего устройства в линии первичной переработки льна / В. А. Шаршу-нов, В. Е. Кругленя, А. С. Алексеенко, В. А. Левчук, М. В. Цайц // Весщ Нацыянальнай акадэмп навук Беларуси Серыя аграрных навук. - 2015. - № 3. - С. 112-117.

2. Устройство для отделения семенных коробочек от стеблей льна: пат. 7224 Респ. Беларусь, МПК А 01F 11/00 / В. Е. Кругленя, М. В. Лёвкин, В. И. Коцуба, С. Н. Крепочин, В. А. Левчук; заявитель Белорус. гос. с.-х. акад. - №u 20100607; заявл. 06.07.10; опубл. 02.02.11 // Афщыйны бюл. / Нац. цэнтр штэлектуал. уласнасщ. - 2011. - № 2 - С. 154.

3. Устройство для выделения семян из ленты льна: пат. 8183 Респ. Беларусь, МПК А0Ш 45/06 / В. Е. Кругленя, М. В. Лёвкин, В. А. Левчук; заявитель УО «БГСХА». - № u 20110743; заявл. 29.09.11; опубл. 30.04.12 // Афщыйны бюл. / Нац. цэнтр штэлектуал. уласнасщ. - 2012. - №2 - С. 197.

4. Курзенков, С. В. Методика расчета параметров слоя стеблей льна в зоне обмолота / С. В. Курзенков, В. А. Левчук, М. В. Цайц // Вестник Белорусской государственной сельскохозяйственной академии. - 2022. - № 1. - С. 154-159.

5. Курзенков, С. В. Обоснование скорости зажимного транспортера обмолачивающего устройства линии первичной переработки льна / С. В. Курзенков, В. А. Левчук, М. В. Цайц // Агропанорама. - 2022. - № 1(149). - С. 14-19.

6. ГОСТ 33734-2016. Межгосударственный стандарт. Техника сельскохозяйственная. Комбайны и машины для уборки льна. Методы испытаний. - М.: Стандартинформ, 2017. - 73 с.

7. Основы расчета рабочих органов машин и оборудования для производства семян льна / В. А. Шаршунов, В. Е. Круг-леня, А. Н. Кудрявцев [и др.]. - Горки: Белорусская государственная сельскохозяйственная академия, 2016. - 156 с.

8. Ковалев, М. М. Технологии и машины для комбинированной уборки льна долгунца: Дис. ... док. техн. наук. - Тверь. 2010. - 615 с.

9. Шаршунов, В. А. Уравнения траектории движения рабочего органа обмолачивающего устройства колебательного типа линии первичной переработки льна «Van Dommele» / В. А. Шаршунов, А. С. Алексеенко, С. В. Курзенков, В. А. Левчук, М. В. Цайц // Вестник Белорусской государственной сельскохозяйственной академии. - 2017. - № 4. - С. 164-168.

i Надоели баннеры? Вы всегда можете отключить рекламу.