Воронеж, Воронежская область, Россия
Воронеж, Воронежская область, Россия
Воронеж, Воронежская область, Россия
Россия
с 01.01.2018 по 01.01.2019
Воронеж, Россия
Аэросев лесных семян с использованием беспилотных летательных аппаратов (БПЛА) представляет собой эффективный метод быстрого и недорогого лесовосстановления, что особенно актуально в условиях современных экологических вызовов. Обзор исследований и технической оснащенности в этой области показал, что сейчас отсутствуют отечественные высевающие аппараты для БПЛА способные осуществлять точный посев дражированных лесных семян с их внедрением в поверхностный слой почвы. Цель данного исследования заключается в разработке и испытании высевающего аппарата точного высева для аэросева дражированных лесных семян с БПЛА, а также изучение его рабочих режимов для определения параметров эффективного заглубления семян. Материалы и методы исследования включают 3D-моделирование геометрии высевающего аппарата с использованием систем автоматизированного проектирования (САПР). Моделирование выполняется с использованием метода дискретных элементов (DEM), представляющего почву в виде сферических частиц разной степени связности, а индентор в виде твердой сферы с заданными физическими свойствами (масса, положение, скорость, сила). Взаимодействие частиц описывается моделью контакта с пружиной и демпфером. Ударное взаимодействие осуществляется при скоростях индентора 25, 50 и 75 м/с с несвязной, срендесвязной и связной почвами. Лабораторный эксперимент проводится на несвязных почвах при скорости движения индентора 75 м/с. Результаты показали, что, согласно 3D-модели, масса высевающего аппарата составила 1140 г, а масса изготовленного экспериментального образца 1074 г. По результатам DEM-моделирования установлено, что проникновение индентора наблюдается только на несвязных почвах при скоростях 75 м/с, однако для исключения выброса индентора из образовавшейся воронки достаточно скорости 25 м/с. В случае среднесвязных и связных почв наблюдается только незначительный боковой охват индентора почвенными частицами или нахождение на поверхности в пределах образованного углубления. Экспериментальная проверка выявила ряд расхождений с данными моделирования, так глубина проникновения при моделировании составила 17,5 мм, а в лабораторном эксперименте 24,6 мм. При этом средний диаметр образуемых воронок в верхней части напротив был больше для данных моделирования и составил 53,5 мм, а для экспериментальных данных 19,3 мм. Для уменьшения расхождений в дальнейшем необходимо модифицировать модель с учетом влияния влажности, гранулометрического состава и объемной плотности естественных почв.
высевающий аппарат, аэросев, БПЛА, имитационное моделирование, системы автоматизированного проектирования (САПР), испытания, метод дискретных элементов (DEM)
Текст (PDF): Читать Скачать
Текст (PDF): Читать Скачать
Текст (PDF): Читать Скачать
Текст (PDF): Читать Скачать
Введение
Введение Интенсификация процессов глобального обезлесения и климатических изменений актуализирует необходимость разработки высокотехнологичных методов лесовосстановления [1]. Технология аэросева с использованием беспилотных летательных аппаратов (БПЛА) представляет особый интерес благодаря операциональным преимуществам и экономической эффективности - прямой посев в 10-30 раз дешевле традиционных методов [1]. Однако эффективность технологии ограничена комплексом факторов, снижающих всхожесть семян до уровня менее 10% [1].
Проведенный анализ оборудования ведущих компаний (DroneSeed, AirSeed, Flash Forest, Dronecoria) выявил четыре основных типа высевающих аппаратов. Наибольший потенциал демонстрируют аппараты точного высева с ускорением дражированных семян потоком воздуха, используемые компаниями Dendra Systems, Flash Forest, AirSeed Technologies и Ghaf Tree Project.
Согласно лесоводческим рекомендациям для Российской Федерации, аэросев может применяться на участках, пострадавших от пожаров, с супесчаными и хорошо дренированными суглинистыми почвами, где огонь полностью уничтожил лесную подстилку [2], а также на транспортно-удаленных лесных участках и склонах, куда затруднен доступ тракторной техники [3].
Для моделирования процессов ударного взаимодействия объектов сферической формы с дискретной средой широко используются экспериментальные и DEM-методы. Фундаментальные основы процессов ударного проникновения в сыпучие гранулированные слои были заложены в работе Hou M. и др. [6] (рисунок 1а), где методами высокоскоростной съемки детально исследована динамика формирования полости и выброса частиц.
Значительный вклад в понимание влияния механических свойств индентора внесло исследование Ye X. и Van Der Meer D. [7], продемонстрировавшее особенности кратерообразования и поведения при ударе для деформируемых гидрогелевых сфер в зависимости от их модуля Юнга.
Особую актуальность для моделирования реальных почвенных условий представляют работы Zhang X. и коллектива авторов [8-10], в которых последовательно исследуется влияние капиллярных эффектов и всасывания в влажных гранулированных материалах на динамику удара, силу сопротивления и глубину проникновения сферического индентора.
Важным аспектом при моделировании поведения семян в почве является учет вращательного движения. Исследование Carvalho D.D. и др. [11] специально посвящено анализу проникновения вращающейся сферы в гранулированную среду, где установлено существенное влияние вращения на траекторию и характеристики внедрения.
Проведенный анализ современных исследований, выполненных методом дискретных элементов (DEM), демонстрирует широкие возможности данного подхода для изучения динамики ударного взаимодействия в гранулированных средах.
В исследовании Ye X. et al. [12] (рисунок 1б) методом DEM изучен наклонный удар снаряда о гранулированную среду. Авторы установили количественные зависимости глубины проникновения от угла атаки и детально проанализировали динамику процесса.
Работа Shen W. et al. [13] посвящена анализу влияния межчастичного трения и демпфирования на динамику удара сферического снаряда о грунтовое основание. Авторы продемонстрировали, что изменение коэффициента трения между частицами на 30% приводит к значительному изменению глубины проникновения (до 25%).
В исследовании Liang S. et al. [14] подробно изучены механизмы передачи и диссипации энергии при ударе индентора о гранулированный материал. Установлено, что до 60% начальной кинетической энергии рассеивается за счет работы сил трения между частицами.
Значительный интерес представляет работа Zhang X. et al. [15], в которой исследовано влияние формы частиц на динамику ударного взаимодействия. Показано, что использование несферических частиц в моделировании позволяет более точно воспроизводить экспериментальные данные.
Важным направлением развития DEM является применение гибридных методов, что продемонстрировано в работе Lin J. et al. [16]. Авторы использовали CFD-DEM подход для моделирования проникновения снаряда в насыщенные жидкостью гранулированные слои.
Исследование Wada K. et al. [17] демонстрирует возможности DEM для моделирования высокоскоростных ударов (100-900 м/с) и кратерообразования. В работе Ghazi Alshanti W. [18] изучено проникновение вращающегося снаряда в гранулированные системы.
Исследования Sunday C. и соавт. [19, 20] представляют систематический анализ влияния гравитационных условий на процессы ударного взаимодействия в гранулированных средах. В работе [19] проведено сравнительное исследование производительности различных реализаций метода дискретных элементов (DEM) для моделирования ударных процессов в условиях различной гравитации. В исследовании [20] разработана масштабируемая модель столкновений для медленных взаимодействий, учитывающая гравитационные эффекты.
Исследование Ye X. и Zhang C. [21] посвящено анализу влияния геометрии индентора на силовое воздействие и реологические свойства гранулированных сред при ударном взаимодействии. Авторы выявили, что в квазистатическом режиме сила сопротивления не зависит от скорости и может быть описана через расширенную модель локального трения, учитывающую площадь контакта.
Работа Wang F. и соавт. [22] представляет комплексное исследование модели силы сопротивления при вертикальном проникновении в гранулированную среду, основанное на комбинированном использовании DEM-моделирования и экспериментальных данных. Особенностью работы является системный учет таких ключевых параметров, как форма носовой части индентора, скорость воздействия, диаметр частиц гранулированной среды и коэффициент трения между ними. Проведенный анализ выявил, что результаты существующих исследований нельзя использовать напрямую для описания процессов ударного взаимодействия дражированных семян с почвенной средой из-за значительных различий в свойствах сред, инденторов и режимов взаимодействия.
Рисунок 1. Методы исследования процессов ударного
взаимодействия объектов с дискретной средой
Figure 1. Methods for studying the processes of impact interaction of objects with a discrete medium
Источник: А – Hou M. и др. [6], Б – результаты моделирования авторов
Source: A – Hou M. et al. [6], B – authors' own modeling results
Таким образом, существует потребность в специализированном исследовании ударного взаимодействия дражированных семян с почвенными средами.
Целью данного исследования является разработка и испытание высевающего аппарата точного высева для аэросева дражированных лесных семян с БПЛА, а также изучение его рабочих режимов для определения параметров эффективного заглубления семян.
Материалы и методы
Объектом исследования является пневматический высевающий аппарат точного высева для аэросева лесных семян с БПЛА.
В данной работе спроектирована параметрическая 3D-модель высевающего аппарата точного высева дражированных семян с ускорением потоком воздуха (рисунок 2), который состоит из направляющего канала 1, корпуса 2, крышки корпуса 3, двигателя 4, замкового узла 5, семенного бункера 6, заслонки семенного бункера 7, высевающего барабана 8, подшипника 9, датчика положения 10, воздушного распределителя 11, нагрузочной пружины 12 и контактной колодки 13. 3D-модель создавалась адаптированной для быстрого прототипирования с использованием 3D-печати, что позволяет в кратчайшие сроки создавать функциональные прототипы из инженерных полимерных материалов с минимальными массовыми характеристиками [5].
Виртуальный стенд, созданный для исследования процессов взаимодействия дражированного семени с почвенной средой, представляет собой фрагмент почвенного пласта размером 100х100х50 мм (рисунок 3). Плотность частиц составляет 2500 кг/м3 что соответствует плотности твердой фазы почвы, масса частицы равна 0,011 г, диаметр 2 мм. Полученный фрагмент почвенного пласта включает 64547 частиц общей массой 0,676 кг с насыпной плотность 1352 кг/м3 что соответствует плотности лесных почв.
Сферический индентор (далее по тексту дражированное семя) диаметром 6 мм генерируется по центру почвенного пласта на высоте 200 мм от его поверхности. Плотность дражированного семени принята равной 1900 кг/м3 что соответствует усредненной плотности дражировочных смесей, разрабатываемых в данный момент авторами статьи и семян сосны.
Рисунок 2. Высевающий аппарат точного высева с ускорением дражированных семян потоком воздуха
Figure 2. Precision seeding device with acceleration of pelleted seeds by an air flow
Источник: собственная 3D-модель авторов | Source: authors' own 3D model
Рисунок 3. Виртуальный стенд для исследования внедрения дражированного семени в почвенный пласт
Figure 3. Virtual test bench for studying the introduction of pelleted seed into the soil layer
Источник: собственная композиция авторов| Source: authors' own composition
В таблице 1 представлены основные параметры DEM-модели позволяющие воспроизводить свойства несвязных, среднесвязных и связных почв. Для определения этих параметров были приняты величины удельного сцепления (cn) свойственные для рассматриваемых типов почв. Далее при постоянном коэффициенте жесткости (kzh) подбирался порог растяжения ts необходимый для получения заданного удельного сцепления и
коэффициент демпфирования (kdamp) для обеспечения стабильности системы.
Движение дискретных элементов под действием упругих сил (FY) и сил сухого и вязкого трения (FС и FВ) рассчитывается по законам классической динамики (рисунок 4). Моделирование производится в трехмерном декартовом пространстве (x, y, z). Состояние каждого элемента Ei задается шестью переменными: координатами его центра (xi, yi, zi) и компонентами скорости (vxi, vyi, vzi). Более подробное описание математической модели приводится в другой работе авторов [24].
Таблица 1.
Параметры DEM-модели
|
Тип почвы |
Порог растяжения, ts |
Коэффициент демпфирования kdamp, Н·с/м |
Коэффициент жесткости kzh, Н/м |
Получаемое удельное сцепление cn, кПа |
|
Несвязная |
1,05 |
30 |
6,0 ⋅ 104 |
0 |
|
Среднесвязная |
1,12 |
40 |
6,0 ⋅ 104 |
15000 |
|
Связная |
1,18 |
50 |
6,0 ⋅ 104 |
30000 |
Рисунок 4. Схемы силового контакта дискретных элементов:
А – контакт двух дискретных элементов; Б – контакт дискретного элемента и твердой поверхности
Figure 4. Diagrams of force contact of discrete elements
A – contact of two discrete elements; B – contact of a discrete element and a solid surface
Источник: собственная композиция авторов| Source: authors' own composition
Мониторинг процесса моделирования производился с использованием цветовых поверхностных эпюр скоростей движения частиц. Синий цвет соответствует нулевой скорости, красный скорости в 1м/с и более, оттенки зеленого – промежуточные значения (рисунок 5, а). Контроль параметров производился в центральном срезе в продольно-вертикальной плоскости, проходящем через ось движения частицы (рисунок 5, б).
фиксировалась по нижней точке семени. Также измеряется наибольшая глубина воронки. Для лабораторных исследований процесса взаимодействия дражированных семян с почвенной средой бы создан испытательный стенд (рисунок 6). Он состоит из основания 1 с кронштейнами крепления камеры и экрана, подвижного кронштейна с механизмом регулировки угла отклонения от вертикали 2, высевающего аппарата 3, крепления для фиксации аппарата 4, хронографа 5, камеры для скоростной съемки 6, емкости для почвы 7.
Рисунок 5. Мониторинг процесса моделирования:
А – цветовые поверхностные эпюры; Б – центральный срез
Figure 5. Monitoring the modeling process:
A – color surface diagrams; B – central section
Источник: собственная композиция авторов| Source: authors' own composition
Для измерения скорости использовался хронограф WoSporT WST-X01. Он позволяет исследовать объекты диаметром 0,01-20 мм, с частотой до 120 шт/с, диапазон измеряемых скоростей 0-800 м/с при погрешности до 0,3 %. Для съемки процесса применялась USB камера со скоростью съемки 260 fps. Разгон дражированного семени диаметром 5,95мм и весом 0,22 г осуществлялся страйкбольным воздушным компрессором. Скорость полета регулировалась подбором жесткости пружины компрессора.
Рисунок 6. Лабораторный стенд
Figure 6. Laboratory stand
Источник: собственная композиция авторов| Source: authors' own composition
Для проведения опыта использовалась несвязная песчаная почва, взятая непосредственно с участка подлежащего лесовосстановлению (рисунок 7). Гранулометрический состав типичный для песчаных почв сосновых лесов Воронежской области: песок 0,063…2 мм от 70% до 95%; пыль 0,002…0,063 мм от 5 до 25%; глина < 0,002 мм менее 5%. Влажность почвы 10 % по массе.
Поверхность выравнивалась вровень с боковыми стенками короба. В опытах выполнялась серия отстрелов, в которой для исключения взаимного влияния выдерживалось расстояние между соседними воронками и стенками не менее 50 мм.
Для фиксации параметров образованных кратеров и глубин проникновения семени применялся штангенинструмент. Измерение диаметра кратера в верхней части выполнялось с использованием штанги штангенглубиномера ШГ, укладываемой на бортики ящика. Для измерения глубин проникновения дражированного семени применялся глубиномер штангенциркуля ШЦ-1, устанавливаемый опорной поверхностью на штангу штангенглубиномера. Сначала измерялся диаметр кратера в двух взаимоперпендикулярных направлениях и вычислялось среднее значение, а затем выполнялся частичное раскапывание дражированного семени и контроль глубины проникновения.
Регулировка скорости высева при проверке стабильности работы высевающего аппарата осуществлялась изменением скорости вращения шагового двигателя высевающего барабана.
В качестве показателя однородности данных был принят коэффициент вариации (CV). Число повторов опытов при моделировании и эксперименте определялось пилотной серией экспериментов, проводившейся до достижения CV менее 10%.
Рисунок 7. Исследование проникновения сферического индентора в почву
Figure 7. Study of penetration of a spherical indenter into soil
Источник: собственное фото авторов| Source: authors' own photo
1. Mohan M. et al. UAV‐Supported Forest Regeneration: Current Trends, Challenges and Implications. Remote Sens. 2021, 13: 1-30. – DOI: https://doi.org/10.3390/rs13132596.
2. Приказ Министерства природных ресурсов и экологии Российской Федерации от 04.12.2020 № 1014 «Об утверждении Правил лесовосстановления, состава проекта лесовосстановления, порядка разработки проекта лесовосстановления и внесения в него изменений»: зарегистрировано в Мин-ве юстиции РФ от 18.12.2020 г., № 61556.
3. Каляшов В. А., До Т.А., Хитров Е. Г., Григорьева О. И., Гурьев А.Ю., Новгородов Д.В. Современные системы машин и технологии заготовки древесины и лесовосстановления в условиях горных лесосек. Resources and Technology. 2022; 19: 1-47. DOI:https://doi.org/10.15393/j2.art.2022.6163.
4. Lysych M.N., Chernyshev V.V., Nagaytsev V.M. Design and Simulation of Seed Metering Device for Aerosowing of Forest Pelleted Seeds. J. Phys. Conf. Ser. 2021; 2032: 1-7. – DOI: https://doi.org/10.1088/1742-6596/2032/1/012062.
5. Lysych M., Bukhtoyarov L., Druchinin D. Design and Research Sowing Devices for Aerial Sowing of Forest Seeds with UAVs. Inventions. 2021; 83: 1-26. – DOI: https://doi.org/10.3390/inventions6040083.
6. Hou M., Peng Z., Liu R., Liu K., Li J., Li H. Projectile Impact and Penetration in Loose Granular Bed. Science and Technology of Advanced Materials. 2005; 6(7): 855-859. – DOI:https://doi.org/10.1088/1468-6996/6/7A08.
7. Ye X., Van Der Meer D. Hydrogel Sphere Impact Cratering, Spreading and Bouncing on Granular Media. Granular Matter. 2021; 23(4): 80. – DOI:https://doi.org/10.1007/s10035-021-01139-w.
8. Zhang X., Li J., Wang Y., Chen H., Liu R., Zhao T. Dynamic Characteristics of Sphere Impact into Wet Granular Materials Considering Suction. Powder Technology. 2023; 415: 118147. – DOI:https://doi.org/10.1016/j.powtec.2022.118147.
9. Zhang X., Wang C., Li H., Yang J., Zhou W., Liu Y. The Force and Dynamic Response of Low-Velocity Projectile Impact into 3D Dense Wet Granular Media. Physical Review E. 2024; 109(2): 024901. – DOI:https://doi.org/10.1103/PhysRevE.109.024901.
10. Zhang X., Li J., Wang Y., Chen X., Liu M. Resistance Force Scaling and the Solution for Penetration Depth of Impacting Wet Granular Materials. Journal of Fluid Mechanics. 2025; 987: 45. – DOI:https://doi.org/10.1007/s40571-024-00867-z.
11. Carvalho D.D., Bertho Y., Franklin E. Penetration of a Spinning Sphere impacting a Granular Medium. Physical Review E. 2024; 110(1): 014902. – DOI:https://doi.org/10.1103/PhysRevE.110.014902.
12. Ye X., Zhao T., Wang H., Liu K., Zhang R., Chen L. Projectile Oblique impact on Granular Media: Penetration Depth and Dynamic Process. Powder Technology. 2021; 378: 144-152. – DOI:https://doi.org/10.1016/j.powtec.2020.09.065.
13. Shen W., Zhang J., Liu K., Wang X., Li M., Yang H. Influence of Inter-Particle Friction and Damping on the Dynamics of Spherical Projectile Impacting Onto a Soil Bed. Journal of Terramechanics. 2022; 101: 35-45. – DOI:https://doi.org/10.1016/j.jterra.2022.03.002.
14. Liang S., Xu Y., Wang J. On Energy Transfer and Dissipation of Intruder Impacting Granular Materials based on Discrete Element Simulations. Granular Matter. 2022; 24(3): 88. – DOI:https://doi.org/10.1007/s10035-022-01254-2.
15. Zhang X., Wang Y., Li J., Chen H., Liu R., Zhao T. The Effect of Particle Shape on The Dynamics of Spherical Projectile Impacting into Granular Media. Physical Review E. 2024; 110(3): 034901. – DOI:https://doi.org/10.1103/PhysRevE.110.034901.
16. Lin J., Zhao T., Jiang M. Investigating Projectile Penetration into Immersed Granular Beds Via CFD-DEM Coupling. Powder Technology. 2023; 415: 118151. – DOI:https://doi.org/10.1016/j.powtec.2022.118151.
17. Wada K., Senshu H., Matsui T. Numerical Simulation of Impact Cratering on Granular Material. Icarus. 2006; 180(2): 528-545. – DOI:https://doi.org/10.1016/j.icarus.2005.10.029.
18. Ghazi Alshanti W. Discrete Element Modeling of a Projectile Impacting and Penetrating into Granular Systems. Ballistics. 2019; 1: 345-356. – DOI:https://doi.org/10.5772/intechopen.75550.
19. Sunday C., Umbanhowar P.B., Goldman D.I. The Influence of Gravity on Granular Impacts: I. A DEM Code Performance Comparison. Powder Technology. 2021; 377: 910-919. – DOI:https://doi.org/10.1016/j.powtec.2020.09.043.
20. Sunday C., Umbanhowar P.B., Goldman D.I. The Influence of Gravity on Granular Impacts: II. A Gravity-Scaled Collision Model for Slow Interactions. Physical Review E. 2022; 105(4): 044901. – DOI:https://doi.org/10.1103/PhysRevE.105.044901.
21. Ye X., Zhang C. Impact Granular Media for Intruders with Different Geometries: Force and Rheology. Journal of Fluid Mechanics. 2023; 958: A34. – DOI:https://doi.org/10.1017/jfm.2023.98.
22. Wang F., Zhang H., Li X., Chen J., Liu Y., Yang K. A Drag Force Model of Vertical Penetration into a Granular Medium Based on DEM Simulations and Experiments. Acta Mechanica Sinica. 2024; 40(1): 123156. – DOI:https://doi.org/10.1007/s10409-023-23156-x.
23. Alshanti W.G., Wu Y. Numerical Simulation of Vertical Penetration into Granular Beds. Australian Journal of Basic and Applied Sciences. 2011; 5(1): 182-188.
24. Бухтояров Л. Д. Лысыч М. Н., Малюков С. В. Моделирование взаимодействия почвы c комбинированным орудием – ножевой каток и пара дисковых плугов. Лесотехнический журнал. 2025; 3 (59): 186-200. DOI:https://doi.org/10.34220/issn.2222-7962/2025.3/12.
