ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ КИНЕМАТИЧЕСКИХ И ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ШНЕКОВОГО РАБОЧЕГО ОРГАНА ЛЕСОПОЖАРНОЙ ГРУНТОМЕТАТЕЛЬНОЙ МАШИНЫ
Аннотация и ключевые слова
Аннотация (русский):
В работе представлены результаты исследований рабочих процессов шнековых рабочих органов в различных отраслях производства. Установлено, что исследований технологических параметров шнековых рабочих органов при удалении напочвенного покрова с горючими растительными остатками проведено недостаточно. Предложена новая конструкция лесопожарной машины, включающая шнековые рабочие органы для удаления напочвенного покрова с горючими растительными остатками, вырезные диски для образования почвенного вала перед ротором-метателем, регулируемый кожух-направитель с возможностью изменять направленный поток грунта. Рабочий процесс гидропривода шнековых рабочих органов лесопожарной грунтометательной машины описан системой дифференциальных уравнений поступательного, вращательного движений шнекового рабочего органа и уравнением расходов рабочей жидкости. Рассмотрена задача Коши на интервале t∈(t_0;t_end ). В силу нелинейности исходная система сводится к системе рекуррентных соотношений, заменяя производные искомых функций их конечными аналогами. Эта система неразрешенная относительно старшего члена. Благодаря свойству фредгольмовости операторного коэффициента перед старшим членом можно применить метод расщепления системы на системы в подпространствах уменьшающихся размерностей и решить систему в этих подпространствах. Аналитическое решение математической модели удаления напочвенного покрова шнековыми рабочими органами лесопожарной грунтометательной машины при прокладке противопожарных полос и тушении низовых лесных пожаров позволяет по конечным формулам определить кинематические и динамические параметры шнековых рабочих органов. Проведено имитационное моделирование рабочего процесса шнека на лабораторном виртуальном стенде. Получены зависимости динамических характеристик шнековых рабочих органов при преодолении препятствий. Установлено, что максимальное среднее значение крутящего момента наблюдается при преодолении препятствия высотой 100 мм и скорости вращения рабочего органа 400 об/мин и составляет 1468,49 Н∙м, при котором не нарушается его работоспособность

Ключевые слова:
моделирование, напочвенный покров, шнековый рабочий орган, грунтометательная машина
Текст
Текст произведения (PDF): Читать Скачать

Введение

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

В работе Лобанова Г. Л. [1] математическая модель взаимодействия стеблей растений с ножами шнековых рабочих органов включает уравнения проекций скоростей крайних точек ножа на оси координат в параметрическом виде. Определены численные значения  параметров режущего аппарата в уравнения в программе Microsoft Excel и выявлена траектория  движения ножа при  срезании стеблей растений. Однако, предложенная математическая модель - это система уравнений первого порядка. Она не содержит производных второго порядка, т.е. не учитывает второй закон Ньютона.

В расчетах Хайбуллина Р.Р. [3] на основе анализа математической модели выкопки траншей под фундамент рабочим органом землеройной машины получены зависимости сопротивления грунта, скорости бурения, фрезерования, размыва грунта, а также мощности и энергоемкости. Следует отметить, что данная математическая модель является линейной по искомым функциям x,φ . Линейная модель сильно огрубляет описание реального процесса: для угла φ , близкого к 0, в силу первого замечательного предела допустима  замена sinkφ  на φk , k=0,1,2,…,  но при больших значениях φ  такая замена влечет накопление погрешности в вычислениях.

 В работе Согина А.В. [5] на основании математической модели, включающей уравнения движения шнекового рабочего органа при разработке донных отложений в водоемах, обоснованы его оптимальные параметры. Однако, математическая модель шнекового рабочего органа хотя и является нелинейной по искомой функции скорости V , но VdVdx=dV2 , поэтому сводится к линейной, и функция V  находится в явном виде. Кроме того шнек движется  параллельно своей оси, а в составе лесопожарных грунтометательных машин  ось шнеков расположена  перпендикулярно направлению движения трактора..

В работе Закирова М.Ф. [2] проводились исследования влияния шага шнека на мощность привода питателя малогабаритного шнекороторного снегоочистителя. Шнековый рабочий орган выполнен в виде барабана диаметром 0,3 м, к которому приварена ленточная лопасть шириной 0,7 м и шагом 0,3 м с углом подъема лезвия шнека около 17,8°. Установлено, что мощность привода шнекового питателя расходуется на вырезание снега из массива и на перемещение снега питателем.

       В работах Бартенева И.М. и Nunez-Regueira. L. [4, 8] выполнено математическое моделирование рабочих процессов шнековых рабочих органов лесопожарной машины, с использованием метода конечных элементов, при котором лесной грунт представляется множеством (порядка 6000) шарообразных элементов.    Определена оптимальная высота ленты шнекового барабана в пределах 8-10 см.  При метании грунта доля напочвенного покрова составила 11-15 %, мощность на привод шнека   11-16 кВт. Однако кинематические и динамические характеристики шнековых рабочих органов авторами не изучались.

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

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

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

Материалы и методы

Нами предлагается конструктивно-техноло-гическая схема лесопожарной грунтометательной машины с многофункциональными модулями:  шнековых рабочих органов c зубьями или черенковыми ножами, которые возможно менять модульно с винтовой металлической нарезкой на щеточную, в зависимости от участка и вида почвы в зависимости от участка и вида почвы, ротора метателя  с возможностью замены на роторы-метатели с различной формой и расположением лопаток, осуществлять привод фрез – метателей и шнековых рабочих органов от вала отбора мощности трактора, так и от гидромоторов, регулируемого кожуха-направителя с возможностью изменять направленный поток грунта [10].

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

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

1) Mш+Mгрd2xdt2=Px-Gш+Gгр+Pzf+1+f2cosφ+arcctgfPр+Pтр+1+f2cosφ+arcctgf ;

2) Iпрd2φdt2=qM2πp+Gш+Gгр+PzPf-R+r2+PfcosφPр+Pтр-R+r2+Pfcosφ ;

3) qнnн=qMπdt+ayp+Кр∙dpdt,

где Мш, Мгр – масса шнека и  напочвенного покрова на шнеке, кг;

Gш, Gгр – силы тяжести шнека и  напочвенного покрова на шнеке, Н; 

f – коэффициент трения скольжения;

R – радиус шнека, м;

L – длина шнека, м;

h – глубина удаляемого покрова (0,10…0,15), м;

t – время, с;

Iпр – приведенный к валу гидромотора момент инерции, кг·м²;

b – ширина рабочей грани черенкового ножа, м;

b – угол заточки черенкового ножа, град;  Zн-количествоножей,находящихсявконтактеспочвой;       

α – угол наклона винтовой лопасти от перпендикуляра к оси шнека, град;

φ  – угол поворота вала шнека, рад;

qн,, qм  – рабочие объемы насоса и гидромотора, м3/об;

р – давление рабочей жидкости, Па;

aY – коэффициент утечек, м3/(с.Па);

nH – частота вращения насоса, с-1;

kP – коэффициент податливости упругих элементов гидропривода, м5/(Н.с);

К – удельное сопротивление резанию черенкового ножа;

τ  – сопротивление сдвигу напочвенного покрова, Па;

φтр  – угол трения почвы о винтовую поверхность шнека, град;

s  – шаг винтовой ленты, м.

 

 

Рисунок 1. Расчетная схема шнекового рабочего органа (собственные разработки). Vтр – скорость трактора, Gш, Gгр – силы тяжести шнека и грунта на шнеке, Рz – сила пружины предохранителя, Fтр – сила сопротивления перекатыванию
шнека,
N – реакция поверхности почвы, Мк – крутящий момент гидромотора, Fc – сила сопротивления резанию,

Px – тяговая сила, Pи – сила инерции

Figure 1. Design scheme of the screw working body (own developments). Vtr is the tractor speed, Gsh, Ggr is the gravity force of the auger and the soil on the auger, Pz is the spring force of the safety device, Ftr is the resistance to rolling of the auger, N is the reaction of the soil surface, Mk is the torque of the hydraulic motor, Fc is the cutting resistance force, Px is traction force, Pi – force of inertia

 

2.1. Постановка задачи

Рассматривается задача Коши на интервале tt0;tend :

Mш+Mгрd2xdt2=Px-

Gш+Gгр+Pzf+

Fс1+f2cosφ+arcctgf,

 

Iпрd2φdt2=qМ2πp-FсR+r2+        (1)

Gш+Gгр+Pz-FсcosφRf,

qнnн=qМπdt+aуp+Kpdpdt,

xt0=x0,x't0=x1,

φt0=φ0,φ't0=φ1,

pt0=p0,                                (2)

 

Для решения задачи аппроксимируем производные конечными разностями, для чего отрезок t0;tend  разобьём n  узловыми точками ti  с шагом h : ti+1=ti+h ; так, величина шага равна h=tend-t0n . Выбираем шаг так, чтобы h≠Kpaу-1.  В обозначениях xi=xti , φi=φti , pi=pti  имеем

d2xdt2xi+2-2xi+1+xih2,

dtφi+1-φih2,d2φdt2φi+2-2φi+1+φih2,

dpdtpi+1-pih2.

Тогда систему можно переписать, выделив xi+2  в первом уравнении и φi+2  во втором и умножив на h2 :

Mш+Mгрxi+2=

2Mш+Mгрxi+1-Mш+Mгрxi+

+h2Px-h2Gш+Gгр+Pzf+

h2Fс1+f2cosφi+arcctgf,

Iпрφi+2=2Iпрφi+1-Iпрφi

+h2qM2πpi-h2FсR+r2+              (3)

h2RfGш+Gгр+Pz-Fсcosφi,

0=qMπφi+1-qMπφi+Kppi+1

+haу-Kppi-hqнnн.

 

Рассмотрим отдельно второе и третье соотношения системы. Их можно записать в векторном виде

Ayi+2=Byi+1+Fi,

 (4)

с операторами A,B:R2R2

A=Iпр000,B=2Iпр0qMπKp,

                               (5)

и вектор-последовательностями

yi=φipi,Fi=Fi1Fi2

 (6)                             

в обозначениях:

Fi1=-Iпрφi+h2qM2πpi-h2FсR+r2

+h2RfGш+Gгр+Pz-Fсcosφi,

Fi2=-qMπφi+haу-Kppi-hqнnн.

Необходимые сведения

Далее, рассматривается векторное линейное рекуррентное соотношение (далее, ЛРС) второго порядка

Ayi+2=Byi+1+Fi,   i=0,1,2,…,

(2.1)

с заданными начальными значениями y0 , y1 , где A,B  – операторы, задаваемые квадратными числовыми матрицами, причем detA=0 .

При разрешении ЛРС (2.1) используется факт, что линейный оператор, задаваемый вырожденной квадратной матрицей, фредгольмов [8].

Пусть оператор A  обладает одномерным ядром. Для него вводится проектор Q  в подпространство Coker A , полуобратный оператор A-=A-1(I-Q) , где A=ACoim Adom A ;   фиксируются элементы eKer A , e≠0 , φCoker A .

Вводится в Coker A  скалярное произведение <,>  так, что <φ,φ>=1 .

Используется следующий результат [9].

Лемма 1  Линейное уравнение

Av=w,   vE1dom A,   wE2,  

 

равносильно системе

v=A-w+ce   для любых cC,

<Qw,φ>=0.

Пусть выполнено следующее условие.

Теорема 1  Пусть выполнены условия

1. Выражение <QFi,φ>  определено при всех  i=0,1,… .

2. d=des<QBe,φ>≠0 .

Тогда ЛРС (3.2.1) равносильно системе

yi+2=Kyi+1+Φi,

(2.2)

<QByi+1,φ>+<QBFi,φ>=0,

(2.3)

в обозначениях

K=A-B-d-1<QBA-B,φ>e,

Φi=A-Fi-d-1<QBA-Fi,φ>+<QFi+1,φ>e.

Доказательство. В силу леммы 3.2.1 соотношение (3.2.1) равносильно системе

yi+2=A-Byi+1+A-Fi+cie,

        (2.4)

и равенству (3.2.3), где ci  – искомая последовательность. В равенстве (3.2.3) заменим i  на i+1 :

<QByi+2,φ>+<QFi+1,φ>=0

и вместо yi+2 подставим выражение

(3.2.4):

<QBA-Byi+1,φ>+<QBA-Fi,φ>+ci<QBe,φ>+

+<QFi+1,φ>=0.                     (2.5)

Из (3.2.5) выразив ci  и подставив в (3.2.4), получим соотношение (2.2).

Теорема доказана.

Отметим, что равенство (3.2.3) при i=0  – это условие разрешимости начальной задачи для соотношения (3.2.1):

<QBy1,φ>+<QBF0,φ>=0.

(2.6)

Далее рассматривается скалярное ЛРС второго порядка

ui+2+aui+1+bui=fii=0,1,2,…,

(2.7)

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

λ2++b=0

(2.8)

D=a2-4b  –дискриминант, λ1,2=-aD2  – корни.

Лемма 2 Пусть D>0 . Тогда общее решение ui  соотношения (3.2.7) равно

ui=c1(λ1)i+c2(λ2)i+1Dl=0i-2(λ2)i-1-l-(λ1)i-1-lfl,   i=2,3,….,

где c1 , c2  - произвольные постоянные.

Лемма 3 Пусть D=0 , λ  – корень уравнения (3.2.8). Тогда общее решение ui  соотношения (3.2.7) равно

ui=c1λi+c2iλi-1+l=0i-2(i-l-1)λi-2-lfl,   i=2,3,….

Обе леммы доказываются непосредственной подстановкой решения в (3.2.7).

Об операторе

A=Iпр000 .

Результаты и обсуждение

Решение начальной задачи для выражений (4), (5), (6)

Для вычисления xi  получено соотношение (3), откуда

xi+2=2xi+1-xi+h2 ·

Mш+Mгр-1Px

-h2Mш+Mгр-1 ·

Gш+Gгр+Pzf

+h2Mш+Mгр-1Fс1+f2 ·

cosφi+arcctgf.                 

 

 

(7)                    

Имеем:

d=(QBe,φ)=Kp≠0.

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

φi+2=2φi+1-φi+h2Iпр-1qM2πpi

-h2Iпр-1FсR+r2+

+h2Iпр-1Rf

Gш+Gгр+Pz-Fсcosφi                                                             

pi+2=1-hKp-1aуpi+1-

Kp-1qMπφi+1+hKp-1qнnн.    

   

 

  (8)

 

 

  (9)                       

Из разложения функций x,φ  в ряд Тейлора в окрестности точки t=t0  при h→0

xt0+h=xt0+x't0h+Oh2,

φt0+h=xt0+x't0h+Oh2,

и начальных значений (2) вытекают формулы для нахождения значений x1 , φ1 :

x1=x0+x1h,

                            (10)

φ1=φ0+φ1h.

                          (11)

Далее, из соотношения (10) вытекает формула для вычисления p1 :

p1=1-hKp-1aуp0-

Kp-1qMπφ1-φ0+         

hKp-1qнnн.                                                   

   (12)                    

Теорема 2.

 1. Решение xi  с начальными значениями x0=x0  и x1 , определяемым формулой (10), равно

xi=x0+x1ih+

+h2Mш+Mгр

l=0i-2i-l-1Px

-Gш+Gгр+Pzf+

+1+f2Fсcosφl+arcctgf,

i=2,3,….                                   

                       (13)

2. Значения угла поворота шнека φi  определяются из итерационного процесса (8), где φ0=φ0 , а φ1  из (11).

3. Решение для определения давления рабочей жидкости в гидромоторе pi  с начальными значениями p0=p0   и p1 , определяемым формулой (12), равно

-1-hKp-1aу-1

-Kp-1qMπφi+hKp-1qнnн+

+1-hKp-1aу-1

l=0i-11-hKp-1aуi-1-l

-Kp-1qMπφl+1+hKp-1qнnн .

   (14)

Тем самым решение xt , φt , pt  задачи (1), (2) определяется приближенно в виде значений xi , φi , pi  в узловых точках ti

При участии доц. Лысыча М.Н. проведено имитационное моделирование рабочего процесса шнековых рабочих органов. На рис. 2 представлены зависимости изменения средних значений крутящего момента от скорости вращения шнекового рабочего органа.

Рисунок 2. Зависимости изменения средних значений крутящего момента от скорости вращения шнекового рабочего органа [11]

Figure 2. Dependences of changes in the average values on the rotational speed of the auger working body [11]

 

Анализ полученных зависимостей показал, что крутящий момент гидромотора незначительно возрастает с увеличением скорости вращения шнекового рабочего органа практически по линейному закону. При возрастании частоты вращения от 200 об/мин до 400 об/мин крутящий момент гидромотора увеличивается от 1400 Нм до 1490 Нм.

Характер изменения средних значений крутящего момента от высоты препятствия, преодолеваемого шнеком, представлен на рис. 3.

Рисунок 3. Зависимости изменения средних значений крутящего момента от высоты препятствия [9]

Figure 3. Dependences of the change in the average values on the height of the obstacle [9]

 

С увеличением высоты препятствия крутящий момент заметно возрастает. При этом интенсивность роста увеличивается с увеличением высоты препятствия. Так среднее минимальное значение крутящего момента для высоты препятствия 0 мм и скорости вращения рабочего органа 200 об/мин составляет 880,20 Н∙м. Максимальное среднее значение крутящего момента наблюдается для высоты препятствия 100 мм и скорости вращения рабочего органа 400 об/мин, это 1468,49 Н∙м.

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

 

Заключение

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

Список литературы

1. Лобанов Г. Л., Погоров Т. А. Математическая модель траектории движения ножей шнекового режущего аппарата в зоне резания стеблей растений. Научный журнал Российского НИИ проблем мелиорации. 2017; 1(25). URL: https://cyberleninka.ru/article/n/matematicheskaya-model-traektorii-dvizheniya-nozhey-shnekovo go-rezhuschego-apparata-v-zone-rezaniya-stebley-rasteniy (дата обращения: 14.05.2021).

2. Закиров М. Ф. Исследование влияния шага шнека на мощность привода питателя малогабаритного шнекороторного снегоочистителя. Интеллектуальные системы в производстве. 2015;2: 56-57. URL: https://www.elibrary.ru/item.asp?id=40999966/ (дата обращения: 27.04.2021).

3. Хайбуллин Р.Р. Теоретические основы разрушения грунта при фрезеровании и бурении рабочими органами строительных машин: автореф. дис. ... д-ра техн. наук : 05.05.04. Mосква, 2008. URL: http://www.scholar.ru/speciality.php?spec_id=276 (дата обращения: 15.05.2021).

4. Bartenev I. M., Popikov P. I., Malyukov S. V. Research and development of the method of soil formation and delivery in the form of a concentrated flow to the edge of moving ground forest fire. IOP Conference Series: Earth and Environmental. 2019;226(1): 012052. DOI:https://doi.org/10.1088/1755-1315/226/1/012052.

5. Согин А. В. Средства гидромеханизиции для очистки водных объектов от донных отложений: дис. … д-ра техн. наук. Москва, 2011. 360 с. URL: https://www.dissercat.com/content/sredstva-gidromekhanizatsii-dlya-ochistki-vodnykh-obektov-ot-donnykh-otlozhenii (дата обращения: 16.05.2021).

6. Бартенев И. М., Поздняков А. К. Анализ рабочих органов технических средств для тушения лесных пожаров. Актуальные направления научных исследований XXI века: теория и практика. Воронеж, 2020. № 1 (48). С. 119-122. DOI:https://doi.org/10.34220/2308-8877-2020-8-1-119-122.

7. Whitehouse N. J. Forest fires and insects: palaeoentomological research from a subfossil burnt forest. Palaeogeography, Palaeoclimatology, Palaeoecology. 2000;164(1-4): 231-246. URL: https://www.sciencedirect.com/science/article/pii/S0031018200001887 (дата обращения: 15.05.2021).

8. Nunez-Regueira L. [et al.] Calculation of forest biomass indices as a tool to fight forest fires. Thermochimica Acta. 2001;378(1-2): 9-25. DOI:https://doi.org/10.15372/SJFS20150609.

9. Бартенев И. М. (и др.) Комбинированный лесопожарный грунтомет и рекомендации по его применению. Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2012;84. URL: https://cyberleninka.ru/article/n/kombinirovannyy-lesopozharnyy-gruntomet-i-rekomendatsii-po-ego-primeneniyu (дата обращения: 16.05.2021).

10. Патент № 2684940 Российская Федерация, МПК Е02 F 3/18. Пожарный грунтомет-полосопрокладыватель : № 2016104672; заявл. 09.07.2018; опубл. 16.04.2019 / И. М. Бартенев, П. И. Попиков, С. В. Малюков, С. В. Зимарин, Н. А. Шерстюков; заявитель и патентообладатель ФГБОУ ВО «ВГЛТУ им. Г.Ф. Морозова». URL: https://patents.google.com/patent/RU2541987C1/ru (дата обращения: 01.02.2021).

11. Драпалюк, М. В. (и др.) Математическая модель процесса подачи и выброса грунта рабочими органами комбинированной машины для тушения лесных пожаров. Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2012;84. URL: https://cyberleninka.ru/article/n/matematicheskaya-model-protsessa-podachi-i-vybrosa-grunta-rabochimi-organami-kombinirovannoy-mashiny-dlya-tusheniya-lesnyh-pozharov (дата обращения: 12.03.2021).

12. Гнусов М. А., Попиков П. И., Малюков С. В., Шерстюков Н. А., Поздняков А. К. Повышение эффективности предупреждения и тушения лесных пожаров с помощью лесопожарной машины. IOP: Материаловедение и инженерия. 2020;919(3): 032025. URL: https://cyberleninka.ru/article/n/vliyanie-rezhimov-raboty-lesopozharnoy-gruntometatelnoy-mashiny-s-gidroprivodom-na-pokazateli-effektivnosti (дата обращения: 22.04.2021).

13. Kruggel-Emden H., Wirtz S., Scherer V. (2008) A study on tangential force laws applicable to the discrete element method (DEM) for materials with viscoelastic or plastic behavior. Chemical Engineering Science, 63(6), 1523-1541. URL: https://www.edemsimulation.com/papers/a-study-on-tangential-force-laws-applicable-to-the-discrete-element-method-dem-for-materials-with-viscoelastic-or-plastic-behavior/ (дата обращения: 27.05.2021).

14. Tsunazawa Y., Shigeto Y., Tokoro C., Sakai, M. (2015) Numerical simulation of industrial die filling using the discrete element method. Chemical engineering science, 138, 791-809. DOI:https://doi.org/10.14419/ijet.v7i2.23.11876.

15. Guo Y., Curtis J. S. (2015) Discrete element method simulations for complex granular flows. Annual Review of Fluid Mechanics, 47, 21-46. URL: https://www.annualreviews.org/doi/abs/10.1146/annurev-fluid-010814-014644.

16. Орловский С. Н. Определение энергетических и динамических параметров тракторов, режимов резания активных рабочих органов машинно-тракторных агрегатов: моногр. Красноярск: КрасГАУ, 2011. 376 с. URL: http://lesnoizhurnal.ru/issuesarchive/ ELEMENT_ID=294450 (дата обращения: 21.05.2021).


Войти или Создать
* Забыли пароль?