THEORETICAL STUDY OF THE KINEMATIC AND DYNAMIC CHARACTERISTICS OF THE AUGER WORKING BODY OF A FOREST FIRE GROUND-THROWING MACHINE
Abstract and keywords
Abstract (English):
The paper presents the results of studies of the working processes of screw working bodies in various industries. Studies of the technological parameters of screw working bodies during the removal of the ground cover with combustible plant residues have been insufficiently carried out. A new design of a forest fire machine is proposed, including screw working bodies for removing the ground cover with combustible plant residues, cut-out discs, for the formation of a soil shaft in front of the rotor-thrower, with an adjustable guide casing with the ability to change the directional flow of soil. Working bodies of the forest fire soil-throwing machine are described by a system of differential equations of the translational, rotational movements of the auger working body and the equation of the flow rates of the working fluid. The Cauchy problem is considered on the interval t∈ (t_0; t_end). Original system is reduced to a system of recurrence relations due to its nonlinearity, replacing the derivatives of the desired functions with their finite analogs. This system is unresolved with respect to the senior member. One can apply the method of splitting the system into systems in subspaces of decreasing dimensions and solve the system in these subspaces due to the Fredholm property of the operator coefficient in front of the leading term. The analytical solution of the mathematical model for removing the ground cover by auger working bodies of a forest fire soil-throwing machine when laying fire strips and extinguishing grassland forest fires makes it possible to determine the kinematic and dynamic parameters of the auger working bodies using the final formulas. Simulation modeling of the screw working process was carried out on a virtual laboratory bench. The dependences of the dynamic characteristics of the auger working bodies when overcoming obstacles have been obtained. It was found that the maximum average value of the torque is observed when overcoming an obstacle with a height of 100 mm and the rotation speed of the working body of 400 rpm is 1468.49 N ∙ m, at which its performance is not impaired.

Keywords:
modeling, ground cover, screw working body, soil throwing machine
Text
Publication text (PDF): Read Download

Введение

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

В работе Лобанова Г. Л. [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 Н∙м.

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

 

Заключение

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

References

1. Lobanov G. L., Pogorov T. A. Matematicheskaya model' traektorii dvizheniya nozhey shnekovogo rezhuschego apparata v zone rezaniya stebley rasteniy. Nauchnyy zhurnal Rossiyskogo NII problem melioracii. 2017; 1(25). URL: https://cyberleninka.ru/article/n/matematicheskaya-model-traektorii-dvizheniya-nozhey-shnekovo go-rezhuschego-apparata-v-zone-rezaniya-stebley-rasteniy (data obrascheniya: 14.05.2021).

2. Zakirov M. F. Issledovanie vliyaniya shaga shneka na moschnost' privoda pitatelya malogabaritnogo shnekorotornogo snegoochistitelya. Intellektual'nye sistemy v proizvodstve. 2015;2: 56-57. URL: https://www.elibrary.ru/item.asp?id=40999966/ (data obrascheniya: 27.04.2021).

3. Haybullin R.R. Teoreticheskie osnovy razrusheniya grunta pri frezerovanii i burenii rabochimi organami stroitel'nyh mashin: avtoref. dis. ... d-ra tehn. nauk : 05.05.04. Moskva, 2008. URL: http://www.scholar.ru/speciality.php?spec_id=276 (data obrascheniya: 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. Sogin A. V. Sredstva gidromehanizicii dlya ochistki vodnyh ob'ektov ot donnyh otlozheniy: dis. … d-ra tehn. nauk. Moskva, 2011. 360 s. URL: https://www.dissercat.com/content/sredstva-gidromekhanizatsii-dlya-ochistki-vodnykh-obektov-ot-donnykh-otlozhenii (data obrascheniya: 16.05.2021).

6. Bartenev I. M., Pozdnyakov A. K. Analiz rabochih organov tehnicheskih sredstv dlya tusheniya lesnyh pozharov. Aktual'nye napravleniya nauchnyh issledovaniy XXI veka: teoriya i praktika. Voronezh, 2020. № 1 (48). S. 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 (data obrascheniya: 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. Bartenev I. M. (i dr.) Kombinirovannyy lesopozharnyy gruntomet i rekomendacii po ego primeneniyu. Politematicheskiy setevoy elektronnyy nauchnyy zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta. 2012;84. URL: https://cyberleninka.ru/article/n/kombinirovannyy-lesopozharnyy-gruntomet-i-rekomendatsii-po-ego-primeneniyu (data obrascheniya: 16.05.2021).

10. Patent № 2684940 Rossiyskaya Federaciya, MPK E02 F 3/18. Pozharnyy gruntomet-polosoprokladyvatel' : № 2016104672; zayavl. 09.07.2018; opubl. 16.04.2019 / I. M. Bartenev, P. I. Popikov, S. V. Malyukov, S. V. Zimarin, N. A. Sherstyukov; zayavitel' i patentoobladatel' FGBOU VO «VGLTU im. G.F. Morozova». URL: https://patents.google.com/patent/RU2541987C1/ru (data obrascheniya: 01.02.2021).

11. Drapalyuk, M. V. (i dr.) Matematicheskaya model' processa podachi i vybrosa grunta rabochimi organami kombinirovannoy mashiny dlya tusheniya lesnyh pozharov. Politematicheskiy setevoy elektronnyy nauchnyy zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta. 2012;84. URL: https://cyberleninka.ru/article/n/matematicheskaya-model-protsessa-podachi-i-vybrosa-grunta-rabochimi-organami-kombinirovannoy-mashiny-dlya-tusheniya-lesnyh-pozharov (data obrascheniya: 12.03.2021).

12. Gnusov M. A., Popikov P. I., Malyukov S. V., Sherstyukov N. A., Pozdnyakov A. K. Povyshenie effektivnosti preduprezhdeniya i tusheniya lesnyh pozharov s pomosch'yu lesopozharnoy mashiny. IOP: Materialovedenie i inzheneriya. 2020;919(3): 032025. URL: https://cyberleninka.ru/article/n/vliyanie-rezhimov-raboty-lesopozharnoy-gruntometatelnoy-mashiny-s-gidroprivodom-na-pokazateli-effektivnosti (data obrascheniya: 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/ (data obrascheniya: 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. Orlovskiy S. N. Opredelenie energeticheskih i dinamicheskih parametrov traktorov, rezhimov rezaniya aktivnyh rabochih organov mashinno-traktornyh agregatov: monogr. Krasnoyarsk: KrasGAU, 2011. 376 s. URL: http://lesnoizhurnal.ru/issuesarchive/ ELEMENT_ID=294450 (data obrascheniya: 21.05.2021).


Login or Create
* Forgot password?