Научная статья на тему 'Методология создания единого вибровысевающего аппарата для семян и удобрений'

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

CC BY
101
43
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВЫСЕВ СЕМЯН И УДОБРЕНИЙ / ВИБРАЦИОННОЕ ПОЛЕ / ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКОГО ВЕСА МАТЕРИАЛОВ / РАВНОМЕРНОСТЬ ВЫСЕВА / ВИБРАЦИОННЫЙ АППАРАТ ДИСКРЕТНОГО ДЕЙСТВИЯ / SEED AND FERTILIZER / VIBRATORY FIELD / IDENTIFICATION OF DYNAMIC WEIGHT MATERIALS / UNIFORM SOWING / VIBRATION UNIT DISCRETE ACTION

Аннотация научной статьи по технологиям материалов, автор научной работы — Беспамятнова Наталья Михайловна, Беспамятнов Юрий Алексеевич, Реутин Виталий Валентинович

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

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

Похожие темы научных работ по технологиям материалов , автор научной работы — Беспамятнова Наталья Михайловна, Беспамятнов Юрий Алексеевич, Реутин Виталий Валентинович

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

The article discusses the application of the theory of vibration in the design and creation of devices seeding. The ways of solving the problem of creating a single universal device for seed and fertilizer using vibration processes influencing the change flow rates when sowing seed from the hopper. The ability to transfer the balance of volume of crops and fertilizers at the expiry of their bunker in order to bring the different flow rates for identical simplifies management processes noncontact sowing materials of different physical and mechanical properties

Текст научной работы на тему «Методология создания единого вибровысевающего аппарата для семян и удобрений»

Научный журнал КубГАУ, №113(09), 2015 года

1

УДК 631.17:631.331:631.816.1

05.00.00 Технические науки

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

Беспамятнова Наталья Михайловна д-р техн. наук, профессор,

SPIN-код: 6075-5920

Беспамятнов Юрий Алексеевич SPIN-код: 3197-7431

Реутин Виталий Валентинович канд. техн. наук,

SPIN-код: 6701-6938

Федеральное государственное бюджетное научное учреждение «Северо-Кавказский научноисследовательский институт механизации и электрификации сельского хозяйства»

(ФГБНУ СКНИИМЭСХ),

347740, г. Зерноград Ростовская область, Россия, ул. им. Ленина, д. 14., e-mail: vniptim@gmail. com

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

Ключевые слова: ВЫСЕВ СЕМЯН И УДОБРЕНИЙ, ВИБРАЦИОННОЕ ПОЛЕ, ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКОГО ВЕСА МАТЕРИАЛОВ, РАВНОМЕРНОСТЬ ВЫСЕВА, ВИБРАЦИОННЫЙ АППАРАТ ДИСКРЕТНОГО ДЕЙСТВИЯ

UDC 631.17:631.331:631.816.1 Technical sciences

METHODOLOGY OF THE COMMON VI-BROSOWING MACHINE FOR SEEDS AND FERTILIZER

Bespamyatnova Natalya Mikhailovna Dr.Sci.Tech., professor,

RSCI SPIN-code: 6075-5920

Bespamyatnov Yuri Alekseevich RSCI SPIN-code: 3197-7431

Reutin Vitali Valentinovich Cand.Tech.Sci.,

RSCI SPIN-code: 6701-6938 Federal State Budgetary Scientific Institution Northern-Caucasian Scientific Research Institute of Mechanization and Electrification of Agriculture (FSBSINCSRIMEA)

347740, ZernogradRostov region, Russia,

Lenina st. 14., e-mail: vniptim@gmail. com

The article discusses the application of the theory of vibration in the design and creation of devices seeding. The ways of solving the problem of creating a single universal device for seed and fertilizer using vibration processes influencing the change flow rates when sowing seed from the hopper. The ability to transfer the balance of volume of crops and fertilizers at the expiry of their bunker in order to bring the different flow rates for identical simplifies management processes noncontact sowing materials of different physical and mechanical properties

Keywords: SEED AND FERTILIZER, VIBRATORY FIELD, IDENTIFICATION OF DYNAMIC WEIGHT MATERIALS, UNIFORM SOWING, VIBRATION UNIT DISCRETE ACTION

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

2

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

В СКНИИМЭСХ проведена работа по созданию вибрационного аппарата дискретного действия для высева всех видов сельскохозяйственных культур и удобрений. В основе создания вибрационного аппарата принято положение о том, что сыпучее тело при истечении из бункера находится под воздействием дополнительных вибрационных внешних сил, изменяющих как граничные, так и объемные силы этого тела. В частности, объемный вес высеваемого материала изменяется по величине и направленности [1, 2, 3].

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

При воздействии на истекающий поток материала из бункера дополнительными ускорениями jd тело приводится в некоторое новое поле тяготения, как бы отличное от поля тяготения земли. Объемный вес у сыпучего тела в этом новом поле тяготения будет меньше, чем обычный, принято называть его динамическим объемным весом Уд [3].

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

Динамический объемный вес определяется по зависимости:

Гд = g

( .Л2

1-

J_ g )

(1)

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

3

где у - объемный вес, т/м3;

j - внешнее ускорение, придаваемое сыпучему телу, м/с2; g - ускорение свободного падения, м/с2.

Практически изменять динамический вес можно, изменяя силу и направленность внешнего ускорения вибратором аппарата вибродискретного действия [1, 2].

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

Ранее было установлено, что расширение технологических возможностей высевающей системы достигается в результате динамического равновесия между давлением семян, накопленных на дозаторе, и упругими силами затвора. Частота колебаний затвора зависит от частоты датчика скорости движения сеялки, длительность постоянного импульса определяет амплитуду отключения затвора [1].

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

4wi

k

2 2

Щ _4_ Щ 4 ,2

>

(2)

k

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

4

„1

30

20

10

-10

-20

-30

2®? + -Г 4h

у S 1 0 1 2 14

2а£ ^

В)

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

5

мерности высева семян заключается в обеспечении устойчивого процесса массового расхода семян (удобрений).

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

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

Среднее давление семян или удобрений на пластинку-вибратор можно определить по зависимости:

P

Scp = - = Щ

1

hmf 2 R j

где P - сила удержания материала на пластине затвора кг/м ; S - площадь основания выгрузного отверстия, м ; [4]

H - глубина погружения, м;

у - объемный вес, г/м ;

f - коэффициент внутреннего трения;

(3)

R = W L

гидравлический радиус отверстия, равный в исследуемом

случае единице;

m - масса выделенного элемента (норма высева в единицу времени, массовый расход), кг/с:

m =

SDhg

g

(4)

Ah - высота элементарного объема высеваемого материала, м.

http://ej.kubagro.ru/2015/09/pdf/04.pdf

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

Научный журнал КубГАУ, №113(09), 2015 года

6

В исследуемом случае глубина погружения H ограничивается требованиями АТТ не проводить высев семян из бункера, если в нем не остается менее 15% от общего объема бункера.

Среднее значение напряжения на пластину вибратора можно опре-

делить [3]:

°ср = Rc, (5)

где х - коэффициент, равный:

Х = у + 2/ % -/2 . (6)

На рисунке 2 представлено среднее значение теоретического давления на пластину-вибратор для семян пшеницы, льна и удобрений (селитры).

Оф

2.40 2,20 2,00 1,80 1,60

1.40 1,20

у = 29,246х2 - 42,783х 4-17,303 ^

/

/

/

у = 7,414хг - 7,331х+3,289

у = 6 1 Д31хг- 6,117 1 х 4- 2,929

♦ Пшеница

■ Лён

А Селитра

^“Полиномиальная

(Пшеница)

■Полиномиальная

(Лён)

‘Полиномиальная

(Селитра)

0,2

0,4

0,6

0,8

Рисунок 2 - Среднее теоретическое давление на пластину-вибратор для семян пшеницы, льна и удобрений (селитра)

Из графика (рис. 2) видно, что среднее теоретическое давление на пластину выпускного отверстия для семян пшеницы с высотой слоя растет

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

7

незначительно, а для удобрений давление возрастает почти в 2 раза. Это обстоятельство объясняется тем, что коэффициент трения f селитры составляет 0,85 ... 1,0, а для пшеницы - 0,47 ... 0,73, объемный вес селитры составляет 0,8 ... 1,0 т/м ; а пшеницы - 0,65 ... 0,81 т/м .

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

слоя т - — = Rg, а вертикальные давления в слоях не будут изменяться по

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

_ Rg т/

ношению mn = f, или / f, следовательно, при большом коэффициенте

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

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

Скорость истечения материалов может быть получена по зависимости [3]:

V = 12 gcR -

2 gP

Sg

(7)

где к - коэффициент истечения, для пшеницы 0,67, для суперфосфата - 0,44.

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

8

Нами аналитически и практически установлено, что все виды сельскохозяйственных культур в вибрационном поле приобретают собственные скорости истечения из бункеров, зависимые не от механических размеров, а скорее - физических свойств: динамического объемного веса уд, коэффициента трения f, начального сопротивления сдвигу т0 [1, 2].

Исходя из полученных результатов исследований, было установлено, что физические свойства семян и удобрений следует разделять по принципу их скоростей истечения из бункеров. Полученные по представленным зависимостям данные для бункера сеялки показали, что скорости истечения зерновых культур находятся в пределах 3,15 ... 3,66 м/с; ячменя

4,37 ... 5,16; гречихи и льна - 5,16; проса - 5,3; риса - 4,7 м/с. Для высева удобрений эти значения следующие: для суперфосфата 3,5 м/с; селитры -4,19 м/с.

Семена льна и частицы удобрений обладают большим коэффициентом сдвига (6 ... 15 кг/м ), но у семян льна коэффициент внутреннего трения очень мал (0,34 ... 0,47), а у селитры - значительно выше (0,85 ... 1,0). При этом объемный вес льна составляет 0,65 ... 0,75 т/м , а для селитры и суперфосфата до 1,4 т/м . Скорость истечения для льна составляет 5,16 м/с, для селитры и суперфосфата - 3,5 ... 4,19 м/с.

По результатам факторного эксперимента в соответствии с матрицами ортогонального плана второго порядка были получены адекватные математические модели в виде уравнений регрессии параметров оптимизации для скорости истечения V0, объемного веса у и коэффициента трения f.

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

9

V0=7,2992-3,7892*f-2,7232* у + 15,6082*f*f-25,1197*f* у +11,9688* у * у

а)

V0=6,9666-6,9551*f+1,2107*Y-

26,5415*f*f+114,0648*f*y-111,9591*y*y

в)

V0=9,104+0,0988*f-15,0914* у + 1,9434*f*f-8,7237*f* у +17,5157* у * у

V0=36,9782+59,5545*f-133,5142*y+

49,0479*f*f-163,0152*f*y+154,9588*y*y

г)

Рисунок 3 - Графики поверхностей отклика и изолинии, характеризующие изменение скоростей истечения V0 (м/с) семян и удобрений в зависимости от объемного веса у (т/м ) и коэффициента внутреннего трения f а) пшеница; б) ячмень; в) подсолнечник; г) лен; д) селитра

V0=182,536-906,3217*f+466,9792*y+ 1127,2446*f*f-1141,6733*f*Y+283,8956*Y*Y

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

10

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

Для высева семян пшеницы в выбранных параметрах центр (ядро) изолиний отличается при скоростях 3,2 ... 3,5 м/с, а градиент скорости неравномерен и сильно сужается по краям, что свидетельствует о неравномерности высева семян за границами центра (рис. 3 а).

Для высева семян ячменя центр поверхности отклика уходит заметно ближе к влиянию коэффициента трения, скорость истечения составляет до 5 м/с, что, в общем-то, соответствует отмеченному в практике «плохому» высеву семян ячменя при использовании катушечных аппаратов. Градиент изолинии значительно увеличивается к осям обоих факторов (рис. 3б).

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

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

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

3,5 ... 4,2 м/с (рис. 3д). При изменении скорости истечения градиент скоростей на изолиниях значительно снижается, что вероятно и приводит к уплотнению частиц в бункере и быстрому возникновению сводов в слое.

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

11

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

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

Поскольку изменить физические свойства материалов невозможно, то технически вполне допустимо реальное изменение динамического объемного веса при обеспечении соотношения j/g в заданном направлении. Согласно нашим исследованиям для сельскохозяйственных культур соотношение j/g не превышает 1,3 ... 2,0 [2]. В конкретном случае это соотношение достигалось при частоте вибратора 6 Гц.

На рисунках представлены графики поверхностей отклика и изолинии в исследуемых параметрах оптимизации при подаче внешней частоты колебаний вибратором в процессе истечения семян из бункера с учетом изменения динамического веса высеваемых материалов в вибрационном поле (рис. 4).

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

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

12

уравнений регрессии показывает, что при использовании вибрации, резко снижается влияние дополнительных факторов на скорость истечения материала, что также свидетельствует об устойчивости процесса высева. В лабораторных и полевых условиях неравномерность высева семян пшеницы в исследуемом поле составила не более 1,8 / [1, 2].

V0=0,2907-27,9016*f+45,6183*Y-

8,2133*f*f+59,2677*f*y-70,2364*Y*Y

V0=15,7757-163,5897*f+170,9764*Y+ 268,1802*f*f-431,816*f*Y+143,1811 *y*Y

в)

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

V0=31,4199+58,4145*f-145,7315*Y+

56,5262*f*f-214,7329*f*Y+227,6813*Y*Y

б)

Рисунок 4 - Графики поверхностей отклика и изолинии, характеризующие изменение скоростей истечения V0 семян и удобрений в виброполе при частоте вибратора 6 Гц:

а) пшеница; б) лен; в) селитра

На графиках 3г и 4б - (высев семян льна) также наглядно видно изменение поверхности отклика: центр (ядро) более четко выражен, а гради-

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

13

енты скоростей значительно меньше, чем без вибрации. В полевых и лабораторных исследованиях неравномерность высева семян льна сократилась от 20 до 3% [1, 2].

Особенно заметно изменение изолиний поверхности отклика для высева в виброполе селитры (рис. 3д и 4в). При введении в виброполе частоты 6 Гц градиенты скоростей более стабильны, поэтому изменение перемещений частиц в потоке также следует ожидать более равномерного. В лабораторных и полевых условиях неравномерность высева селитры не превысила 3%, что значительно ниже допустимых по АТТ 10% [2].

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

Библиографический список

1. Беспамятнова, Н.М. Колебания и вибрации в технологических процессах почвообрабатывающих и посевных машин и агрегатов / Н.М. Беспамятнова. - Зерноград: ВНИПТИМЭСХ, 2008. - 224 с.

2. Беспамятнова, Н.М. Вибровысев сельскохозяйственных культур и удобрений / Н.М. Беспамятнова, Ю.А. Беспамятнов, В.В. Реутин, Ю.А. Семенихина. - Зерноград: СКНИИМЭСХ, 2015. - 234 с.

3. Зенков, Р.Л. Механика насыпных грузов / Р.Л. Зенков. - Москва: Машиностроение, 1952. - 216 с.

4. Семенихина, Ю.А. Оптимизация высоты высевных окон вибрационного дозатора семян / Ю.А. Семенихина // Механизация и электрификация сельского хозяйства. - 2011. - №8. - с. 32.

References

1. Bespamjatnova, N.M. Kolebanija i vibracii v tehnologicheskih processah pochvoobrabatyvajushhih i posevnyh mashin i agregatov / N.M. Bespamjatnova. - Zerno-grad: VNIPTIMJeSH, 2008. - 224 s.

2. Bespamjatnova, N.M. Vibrovysev sel'skohozjajstvennyh kul'tur i udobrenij / N.M. Bespamjatnova, Ju.A. Bespamjatnov, V.V. Reutin, Ju.A. Semenihina. - Zernograd: SKNIIM-JeSH, 2015. - 234 s.

3. Zenkov, R.L. Mehanika nasypnyh gruzov / R.L. Zenkov. - Moskva: Mashino-stroenie, 1952. - 216 s.

http://ej.kubagro.ru/2015/09/pdf/04.pdf

Научный журнал КубГАУ, №113(09), 2015 года

14

4. Semenihina, Ju.A. Optimizacija vysoty vysevnyh okon vibracionnogo do-zatora semjan / Ju.A. Semenihina // Mehanizacija i jelektrifikacija sel'skogo hozjajst-va. - 2011. -

№8. - s. 32.

http://ej.kubagro.ru/2015/09/pdf/04.pdf

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