Вестник Дагестанского государственного технического университета. Технические науки. Том 45, №1, 2018 Heraldof Daghestan State Technical University.TechnicalSciences. Vol.45, No.1, 2018 _http://vestnik.dgtu.ru/ISSN (Print) 2073-6185 ISSN (On-line) 2542-095Х_
Для цитирования: Мурадова С.Ш., Федосеева Л.В. Влияние статических сводов на эффективность технологических процессов в бункере. Вестник Дагестанского государственного технического университета. Технические науки. 2018; 45 (1):30-48. DOI:10.21822/2073-6185-2018-45-1-30-48
For citation: Muradova S. Sh., Fedoseeva L.Vl. The infuence of static arches on the efficiency of technologicalprocesses in a hopper. Herald of Daghestan State Technical University. Technical Sciences. 2018; 45 (1):30-48 (In Russ.) DOI: 10.21822/2073-6185-2018-45-1-30-48
ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ МЕХАНИКА
УДК 62-405.6
DOI: 10.21822/2073-6185-2018-45-1-30-48
ВЛИЯНИЕ СТАТИЧЕСКИХ СВОДОВ НА ЭФФЕКТИВНОСТЬ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ В БУНКЕРЕ
1 2 Мурадова С.Ш. , Федосеева Л.В.
Южный университет (Институт управления бизнеса и права), 344068, г. Ростов-на-Дону, пр. Михаила Нагибина, 33а/47, Россия, Донской государственный технический университет, 344000,г. Ростов-на-Дону, пл. Гагарина, 1, Россия, 1e-mail: [email protected],2e-mail:[email protected]
Резюме. Цель. Как показывает практика, процесс истечения сыпучих материалов из бункеров в различных технологических цепочках часто нарушается, в результате чего снижается эффективность переработки содержащихся сыпучих материалов. Одной из распространенных причин остановки истечения является образование динамических сводов в бункере. С этих позиций целью исследования является оценка влияния статических сводов на эффективность технологических процессов в бункере. Метод. Основу исследования составляет система дифференциальных уравнений, описывающая стационарное и равновесное состояние сыпучей среды, которое интегрируется в граничных условиях, соответствующих круглым си-лосам и бункерам. При этом выводятся аналитические выражения для компонент давления в сыпучем материале по различным направлениям. Результат. Полученные решения проверяются на логическую непротиворечивость, иллюстрируются с помощью графиков, построенных в среде символьной математики, для модельных и реальных параметров силосов и бункеров. Производится сравнение результатов математического моделирования с экспериментальными данными для типовых силосов. Разброс экспериментальных значений связан с процессами уплотнения сыпучего материала со временем его выдержки в силосе, с направлениями микровибрации сейсмического или техногенного характера. Сыпучий материал в силосе все время находится в промежуточном положении между стационарными и равновесными состояниями. Параметры сыпучего материала (насыпная плотность, углы внешнего и внутреннего трения) являются некоторыми обобщенными, мало соответствующими конкретному сыпучему материалу и, к тому же, меняющимися с течением времени. Вывод. Подтверждено, что ещё в статике в бункере образуются своды, которые при истечении материала переходят в динамические, что препятствует процессу истечения, понижая экономическую эффективность бункеров.
Ключевые слова: давление, сыпучий материал, дно и стенки бункеров, статические и динамические своды, истечение сыпучего материала
PHYSICAL-MATEMATICAL SCIENCE MECHANICS
THE INFLUENCE OF STATIC ARCHES ON THE EFFICIENCY OF TECHNOLOGICAL PROCESSES IN A HOPPER
1 2 Safura Sh. Muradova , Lyudmila Vl. Fedoseeva
Southern University (Institute of Management, Business and Law),
33a/4 Nagibina Ave., Rostov-on-Don 344068, Russia,
Don State Technical University,
1 Gagarina Sguare, Rostov-on-Don 344000, Russia,
1e-mail: [email protected],2e-mail:[email protected]
Abstract Objectives. In practice, the process of loose material flowing from hoppers in different technological processes is often disrupted, resulting in reduced processing efficiency of contained loose materials. One of the most common reasons for the interruption of flows is the formation of dynamic arches in the hopper. From this perspective, the aim of the study is to assess the effect of static arches on the efficiency of technological processes in the hopper. Methods. The study was based on a system of differential equations describing the stationary and equilibrium state of the loose medium, which is integrated in the boundary conditions corresponding to round silos and hoppers. In this case, the analytical expressions for the pressure components of loose material in different directions are derived. Results. The obtained solutions are checked for logical consistency and illustrated with the help of graphs constructed in the area of symbolic mathematics for modeled and real parameters of silos and hoppers. The results of mathematical modeling are compared with experimental data for typical silos. The scatter of the experimental values is related to the processes of compacting the loose material, the time of its storage in the silo, as well as the microvibrational directions of a seismic or technogenic nature. The loose material in the silo resides in an intermediate position between stationary and equilibrium states at all times. The parameters of the loose material (bulk density, external and internal friction angles) represent some generalised values that do not correspond to a specific loose material, which, in addition, can vary over time. Conclusion. It is confirmed that arches formed in the hopper during the static state acquire a dynamic character during material flow, hampering the flow process and consequently reducing the efficiency of the hoppers.
Keywords: pressure, loose material, hopper bottom and walls, static and dynamic arches, loose material flow
Введение. Конические бункеры, силосы, находят широкое применение в различных сферах человеческой деятельности. Чтобы адекватно рассчитать давление, испытываемое стенками емкости со стороны сыпучего материала, необходимо правильно выбрать модель этого сыпучего материала. При этом очевидна важность аналитического, а не численного, решения, позволяющего прогнозировать изменения давления с изменением, как параметров самих емкостей, так и параметрами сыпучего материала, находящегося в них.
К сожалению, до сих пор нет законченной, общепризнанной модели сыпучих сред. В связи с этим рассматривают в основном два направления.
Первое - сыпучее тело моделируют сплошной средой. Основоположник такой модели французский физик и инженер Coulomb [1] еще в 1783 году сформулировал основные положения так называемой теории предельного равновесия. При этом предусматривается создание строгой теории предельного равновесия. Здесь необходимо отметить работы М.В. Малышева [2] и Г.К. Клейна [3]. Обобщением и развитием работ по статике сплошной среды занимался В.В. Соколовский [4]. Дальнейшее развитие модель сплошной среды получила в работах [5-8] и др. Второе направление - это, так называемая, дискретная модель сыпучего материала. Дискретную модель зернового материала впервые предложил английский инженер С. Дженкин [9] в 1931 г.
Отечественный ученый Л. В. Гячев развил идеи Coulomb и C. Дженкина, предложив дискретную модель сыпучего тела с сухим трением [10-11]. Развитием модели, предложенной Л.В. Гячевым, занимались и за рубежом [12-14].
Постановка задачи. Здесь необходимо отметить ещё одну проблему, как оказывается тесно связанную с данной - это проблема истечения сыпучего материала из бункеров. В большинстве случаев закономерности такого истечения не укладываются в общую схему, разработанную в [10-11]. Очень часто процесс истечения вообще останавливается по непонятным, казалось бы, причинам. Нашим учёным Богомягких В. А. и его школой [15-17] было показано, что всё дело в образовании динамических сводов, которые тормозят и даже останавливают процесс сам истечения. Это направление получило дальнейшее развитие в трудах [18-20].
Методы исследования. В данной работе под сыпучим материалом будем понимать вещество, состоящее из частиц различной величины и формы, взаимодействующих между собой и стенами ограждающей емкости посредством электромагнитных сил, сил сухого и вязкого трения, а также сил сцепления. Однако в отличие от модели Л.В. Гячева, здесь не вводится никаких предварительных предположений о форме частиц и вариантах их укладки в материале.
Если электромагнитные силы, силы вязкого трения и силы сцепления отсутствуют, а проявляются лишь силы сухого трения, то такое состояния вещества в нашей модели будем называть идеальным сыпучим материалом.
Дифференциальные уравнения, описывающие такую модель сыпучей среды, получены нами ранее и опубликованы в [21].
1. Давление сыпучего материала в круглом силосе.
Рассмотрим сыпучий материал в глубоком круглом силосе (рис. 1). Первое уравнение дифференциальной системы (система (11) нашей работы [21]) справедливо и для данного случая:
Рис. 1. Сыпучий материал в круглом силосе и цилиндрическая система координат Fig. 1. The bulk material in a round silo and a cylindrical coordinate system
Проинтегрируем это уравнение по координате р в пределах от 0 до Я:
0 = п• У • д • К2 • бг- п • К2 • бг • (1 + У ?) • дРг/дг + С (1)
Отметим, что уравнение (1) справедливо на цилиндрической поверхности р = Я. Очевидно, что константа интегрирования С по физическому смыслу представляет собой силу трения, действующую на боковую поверхность элементарного слоя со стороны стенки силоса:
Здесь использована связь между вертикальной и горизонтальной компонентами давления через боковой коэффициент. Подставляя (2) в (1), получим
0=n y • g • R2 • dz - n • R2 • dz • (1+ ji,2) • d Pz/dz-j ek-Pz-2 n - Rdz (3)
Преобразуем уравнение (3) к стандартному виду:
дР2 2 /лекР2 уд _
+
dz Я (1 +íif)
Интеграл уравнения (4) имеет вид:
14
(4)
Р, =
YS
R
+ С ■ ехр
(2 ■ fig - к - ,z\
R V • дгу
- 2 - fie -к г ^ Я - (1 + ¡11
Из условия, что на поверхности силоса давление равно нулю, найдем константу интегрирования С. Выражение для давления при этом примет вид:
(5)
Напомним, что выражение (5) справедливо на поверхности р = R.
Соответственно, выражение для горизонтальной компоненты давления будет иметь вид:
Рр =
Y-д -R
R<(6) соответственно решение (5) для равновесного со-
2 - ¡Л е
В случае равновесия ц^ 0, 1 стояния примет вид:
Ря у-д -Я
3 ■ \ V Л // (7)
Очевидно, также, что при этом давление по горизонтали равно давлению по вертикали. Графики зависимости давлений, рассчитанных по формулам (5) и (7) приведены на рис. 2. Из рис. 2 следует, что вертикальное давление в стационарном состоянии (сплошная линия) больше давления в равновесном состоянии (линия из точек). В то же время горизонтальное давление в стационарном состоянии меньше давления в равновесном состоянии (пунктирная линия). Следовательно, в равновесном состоянии стены силоса берут на себя большую нагрузку, в результате чего давление на дно уменьшается.
40000
30000
20000
10000
Р (Па) Àf
jf / Y f J t y s f/ Ji 7.'
10
Z
12 14 16 18 20
Рис. 2. Вертикальное (сплошная линия) и горизонтальное (пунктирная линия) давление в стационарном состоянии и давление (линия из точек) в равновесном состоянии для круглого силоса с сыпучим материалом (R = 3,5 м, у = 160, ф = 200, х = 300, у = 800 кг / м3) Fig. 2. Vertical (solid line) and horizontal (dashed line) pressure in the steady state and the pressure (line of dots) to the equilibrium condition for a circular silo with the bulk material (R =3,5m, у = 160, ф = 200, х = 300, у = 800 kg / m3) Таким образом, в стационарном состоянии (при движении сыпучего материала) на стены силоса оказывается меньшая нагрузка, чем в равновесном состоянии, когда сыпучий материал неподвижен.
При нахождении сыпучего материала в силосе формально он не перемещается и, следовательно, находится в равновесном состоянии. При этом давление в нем необходимо рассчиты-
вать по формуле (7). Но в действительности, в сыпучем материале все время происходят процессы его уплотнения, т.е. он находится в состоянии микродвижения. Процессы уплотнения, слеживания происходят скачкообразно, под действием микроколебаний сейсмического или техногенного характера. Соответственно, сыпучий материал часть времени «проводит» в стационарном состоянии, а часть в равновесном. При этом и давление также меняется скачкообразно. Однако, в решение (5) нет зависимости от горизонтальной координаты. При определении зависимости от координаты Р будем искать ее аналогично тому, как мы это делали для щелевого бункера:
(8)
В таком виде решение (8) удовлетворяет второму уравнению системы (11) работы [21]. Действительно, найдем производные от (8) по координатам:
2 -fie-к
И - ■ t-r vr Р II — —
Эр
2 ■ ■ к
■ ехр
^iVJL (
dz
{1 + fin-R
yg-R
■ [z - <R - p) ■ tg x\
/-leinet a u-ivi.
)2 ■ fig ■ к (l+tf)-*"
¿SX
• exp
(2 ■ ■ к
"(1 + ßf)R
[z-(R-p)tgX]
(1 + rf)-*
2 ■ ■ к
Подставив найденные значения производных во второе уравнение системы (11, [21]),
■
найдем, что оно удовлетворяется при К-'^ВХ — Р-г, откуда следует, что £дх
Кроме того, решение (8) удовлетворяет дифференциальному уравнению (4). Однако решение (8) не удовлетворяет граничному условию на поверхности г = 0.
O
Рис. 3. Деление силоса на области Fig. 3. The division of the silage on the field
Поэтому, для нахождения более точного решения, разобьем весь объем силоса, занятый сыпучим материалом, на области, как показано на рис. 3.
В области I влияние стенок силоса еще не будет сказываться, т.к. эта область расположена выше плоскости естественного откоса, с углом естественного откоса Поэтому в этой области решение будет определяться формулой:
1 I ft
(9)
Соответственно решение в области II можно было бы взять в виде (8). Оно удобно тем, что на боковой поверхности силоса (при р = R) переходит в (5), которое является решением дифференциального уравнения (4).
На поверхности естественного откоса, определяемой уравнениями:
уравнение (8) обращается в ноль.
Однако на этой границе - (10), необходимо, чтобы искомое решение переходило бы в (9). Поэтому решение в области //представим в виде:
(11)
Как видим, уравнение (11) при г = 0 равно нулю, при р = R (т.е. на боковой стенке силоса) переходит в уравнение (5). А при р = Rb (на поверхности откоса, т.е. на границе областей) оно переходит в уравнение (9).
Таким образом, согласование решений в областях I и II выполнено.
Посмотрим, как будет вести себя решение (11) при предельных переходах.
Рассмотрим случай, когда сыпучий материал находится в силосе с абсолютно гладкими стенами (/¿е 0 ).
и™ ^т^'У^-т-рУ^дхи^Ц- -г
1 +Mf
1+iq
l+lif
То есть, сыпучий материал в этом случае ведет себя как безграничная насыпь.
В данном случае влиянием стенок можно пренебречь, что не противоречит физическому смыслу задачи.
В области I решение имеет такой же вид. Если теперь и внутреннее трение стремится к нулю, то решение переходит в решение для идеальной жидкости, что вполне естественно.
Таким образом, решение (11) удовлетворяет граничным условиям.
Оно также удовлетворяет и предельным переходам. Из уравнения (11) вытекает частное решение для случая р = R.
Следовательно, оно описывает условие равновесия элементарного слоя, но не элементарного объема.
На рис. 4 показана зависимость вертикального и горизонтального давления в силосе, R = 3,5 м, в котором находится сыпучий материал (у = 800 кг/м3, у = 160, ф = 200, х = 300) непосредственно вдоль его стенки и по центру силоса.
Как видим, давление в центре несколько выше, чем непосредственно у стенки. Кроме того, и горизонтальное давление меньше вертикального.
Рис. 4. Зависимость вертикального и горизонтального давления от координаты z при р = R = 3,5 м (по стенке силоса) и р = 0 (по осевой линии силоса, у = 800 кг/м3, у = 160, ф = 200, х = 300) Fig. 4. The dependence of the vertical and horizontal pressure on the coordinate z at р = R = 3.5 m (along the wall of the silo) and р = 0 (centerline of the silo, у = 800 kg/м3, у = 160, ф = 200, х = 300)
На рис. 5 представлены изолинии вертикального давления от координат при различных значениях координаты z в областях I и II.
Из данных рис. 5 также следует, что давление в центре силоса несколько выше, чем у стенок. Здесь также можно видеть тенденцию к образованию статических сводов.
Рис. 5. Изолинии вертикального давления на разных горизонтах
3 0 0 0
R = 3,5 м, Y = 800 кг/м , у = 16 , ф = 20 , X = 30 Fig. 5. Contours of vertical pressures at different horizons
3 0 0 0
R = 3,5m, у = 800 kg/m , у = 16 , ф = 20 , х = 30
Кроме того, из рис. 5 следует, что в круглом силосе изолинии вертикального давления не являются горизонтальными. Следовательно, пласты сыпучего материала, формирующиеся под действием вертикального давления, также не будут горизонтальными. В данном случае они уже больше напоминают куполообразные поверхности сводов в сыпучем материале. Образовавшиеся своды уже будут оказывать некоторое дополнительное сопротивление движению сыпучего материала в круглой вертикальной трубе.
На рис. 6. представлен трехмерный график зависимости вертикального давления для того же сыпучего материала в областях I и II . Из графика (рис.6) следует, что при г = 0 вертикальное давление также равно нулю. Кроме того, наглядно видна «сшиваемость» решений на границе областей, а также видно, что давление в центре силоса (р = 0) несколько превышает давление у его стен (р = Я = 3,5 м).
Рис. 6. Трехмерный график вертикального давления в силосе
3 0 0 0
R = 3,5 м, у = 800 кг/м , у = 16 , ф = 20 , х = 30 Fig. 6. Three-dimensional graph of the vertical pressure in the silo
3 0 0 0
R = 3,5m, y = 800 kg/m , у = 16 , ф = 20 , х = 30
На рисунках 7 - 8 представлены графики сравнения теоретических кривых (непрерывные линии) и экспериментальных значений. При этом экспериментальные данные взяты из [22].
На рис. 7 представлена зависимость горизонтального давления на стену круглого железобетонного силоса высотой H = 21 м и радиусом R = 3,25 м.
30000 25000 20000 15000 10000 5000
0|.....................
0 5 10 15 20
г (м)
Рис. 7. Теоретическое горизонтальное стационарное (сплошная линия) и равновесное (линия из точек) давление в круглом железобетонном силосе (Н = 21 м R = 3,25 м х = 300, ф = 200, у = 200). Вертикальные линии - разброс экспериментальных значений Fig. 7. Theoretical horizontal stationary (solid line) and equilibrium (line of dots) pressure in a circular, reinforced concrete silos (Н = 21 m, R = 3,25 m, х = 300, ф = 200, у = 200). Vertical lines the spread of experimental values На рис. 8 представлена зависимость силы давления на дно стального круглого силоса, высотой Н = 6 м и радиусом R = 0,5 м.
Рис. 8. Теоретическая сила давления на дно стального круглого силоса (Н = 6 м R = 0,5 м) от высоты засыпки. Сплошная линия - стационарное состояние, линия из точек - равновесное (х = 30°, ф = 20, у = 20°). Вертикальные линии - разброс экспериментальных значений Fig. 8. The theoretical pressure force to the bottom steel round silo (H = 6 m R = 0.5 m) to the height of backfill. The solid line is the steady state, a dotted outline - the equilibrium (х = 300, ф = 200, у = 200). Vertical lines - the spread of experimental values
На рис. 9 представлена зависимость силы трения, действующей на боковую стенку стального круглого силоса, высотой Н = 6 м, радиусом R = 0,5 м. Как видно из рисунков, разброс экспериментальных данных достаточно велик. Кроме того, как отмечается в цитируемом источнике [22], само горизонтальное давление разное по разным направлениям на стороны света.
Рис. 9. Зависимость силы трения на боковую стенку круглого, стального силоса (Н = 6 м, R = 0,5 м, х = 300, ф = 200, ф = 160) от высоты засыпки.
Стационарное состояние - сплошная линия, равновесное - линия из точек. Вертикальные линии - разброс экспериментальных значений Fig. 9. The dependence of the friction force on the side wall round steel silo (H = 6 m, R = 0.5 m, х = 300, ф = 200, ф = 160) of the height of the backfill.
Stationary state - a solid line, the equilibrium line of dots.
Vertical lines the spread of experimental values
Очевидно, что этот разброс экспериментальных значений связан с процессами уплотнения сыпучего материала, со временем его выдержки в силосе, а также с тем, с какого направления идут микровибрации сейсмического или техногенного характера. Сыпучий материал в силосе все время находится в промежуточном положении между стационарными и равновесными состояниями. Кроме того, сами параметры сыпучего материала (насыпная плотность, углы внешнего и внутреннего трения) являются некоторыми обобщенными, мало соответствующими конкретному сыпучему материалу и, к тому же, меняющимися с течением времени. Поэтому более точного согласия с экспериментальными данными ожидать трудно. Но из сравнения с экспериментальными данными, все же можно сделать вывод, что состояние сыпучего материала ближе к стационарному состоянию, чем к равновесному.
2. Давление сыпучего материала в коническом бункере.
Исследуем давление в сыпучем материале, находящемся в коническом бункере (рис. 10).
Рис. 10. Параметры конического осесимметричного бункера Fig. 10. Parameters of conical axisymmetric hopper
Уравнение (1) для силоса будет справедливо и в данном случае. Проинтегрируем это уравнение по координате р в пределах от 0 до Rd:
0 = m y • g • R2d • dz -n • R2d • dz • (1 +jj,2) • ^ + C
(12)
Отметим еще раз, что уравнение (12) справедливо на поверхности р = Rd. Очевидно, что константа интегрирования - С (по аналогии со щелевым бункером) по физическому смыслу представляет собой силу трения, действующую на боковую поверхность элементарного слоя со стороны стенки силоса. То есть, выражение для константы будет иметь аналогичный
вид:
В выражениях (12) и (13) Rd - точка на боковой поверхности бункера (рис. 10), связанная с координатой г соотношением:
Здесь, как и раньше, для сокращения записи введен параметр b = tga Подставим теперь (14) и (13) в (12):
0=н-y -g -R2d • dz - п -R2d • dz • (1 + щ2) -2n • Rd • \e • Rz (1 - (1 - k) • cosa) -dz
Преобразуем полученное выражение:
(14)
дР2
dz
+
2 ■ /ле ■ (1 - (1 - k) ■ cosa) p y-g
(15)
Если стенки вертикальны - а = 0 , то уравнение (15) переходит в уравнение (4).
Решение уравнения (15) имеет вид:
п г , г, и У -9 ■ ( Я + Ь
Р2 = С ■ (-/? +Ь-г) ь + ' в ^ *
где А, для сокращения записи, имеет вид:
Ь = tga А =
(1+/I?) ( 2 А+Ь) ¡ie -(1 - (1 - fcD- cosoiï
1+ДГ
Из условия, что на поверхности бункера давление равно нулю, найдем константу интегрирования С:
у д-Я за
С =
(Я) Ъ
(2 А-ы-а+рр
С учетом значения константы С, выражение для вертикального давления в бункере будет иметь вид:
■л г-д-я /■.!?-(-" I''-ю
Преобразовывая это выражение, получим:
B R Л- ^Ук ^
(2 A-by[l+tf)
Pz <2Э
ЧгГ)
2-А-Ь V \ к ) ) (16)
Выражение (16) описывает вертикальную компоненту давления в коническом бункере. Оно подобно аналогичному выражению, описывающему давление в щелевом бункере. Поэтому можно сказать, что оно описывает давление в стационарном состоянии.
В равновесном состоянии, при ¡и^ 0, х ^ 0, к^ 1,А ^ ¡ие , давление в коническом бункере будет выражаться как:
Рв сю =
у-д-R
2 Це-Ъ \ Д / \ V К ) / (17)
Из выражения (17) видно, что при Ъ^ 2ие, а ^аг^(2аг^ф) возникает неопределенность, типа 0/0. Раскрыв ее по правилу Лопиталя, получим:
р _ Y9R п
"okp---— (1
bz
R
) In (1
bz
R
)
(18)
На рисунке 11 приведены кривые равновесного давления в коническом бункере (сплошные линии), рассчитанные по формуле (17), и кривая критического равновесного давления (линия из точек Ъ = 2ие), рассчитанная по формуле (18).
Мы видим, что никаких разрывов не наблюдается.
14000
12000
10000
16000
18000
2000
8000
6000
4000
0
0 2 4
G
8
z
Рис. 11. Кривые равновесного давления в коническом бункере (сплошные линии) при разных углах наклона стенки и кривая критического равновесного давления (b = 2^e) в том же бункере (кривая из точек).
R = 3,5 м, у = 160, ф = 200, х = 300, у = 800 кг / м3 Fig. 11. Curves of equilibrium pressure in a conical hopper (solid lines) at different angles of inclination of the wall and the curve of the critical equilibrium pressure (b = 2^e) in the same bunker (curve of dots).
R = 3.5 m, у = 160, ф = 200, х = 300, у = 800 kg / m3
На рис. 12 приведены кривые веритикального (сплошная линия) и горизонтального (пунктирная линия) давления в стационарном состоянии и давление в равновесном состоянии (линия из точек) для конического бункера.
Это указывает на то, что в стационарном состоянии давление на стены меньше, чем в равновесном состоянии. Отметим что здесь, также как и в силосе, вертикальное давление в стационарном состоянии больше давления в равновесном состоянии. В то же время горизонтальное давление в стационарном состоянии меньше давления в равновесном состоянии.
Рис. 12. Вертикальное (сплошная линия) и горизонтальное (пунктирная линия) давление в коническом бункере в стационарном случае и давление в равновесном состоянии (линия из точек) R = 3,5 м, у = 160, ф = 200, х = 30°,у = 800 кг / м3 Fig. 12. Vertical (solid line) and horizontal (dotted line) pressures in a conical hopper in the stationary case, and the
pressure in the equilibrium state (line of dots) R = 3.5 m, у = 160, ф = 200, х = 300, у = 800 kg / m3
G
z
На рис. 13 представлены графики зависимости вертикального давления вдоль стенки от координаты z для щелевого (сплошная линия) и конического (линия из точек) бункеров. При этом параметры бункеров одинаковы: L = R = 3,5 м, одинаков и угол наклона стенки а = 300, одинаков и сыпучий материал в бункерах.
Заметим, что графики подобны, хотя есть различия и количественные (давление в щелевом бункере больше давления в коническом) и качественные (максимум давления в щелевом бункере сдвинут ближе к выходному отверстию). Однако здесь различий меньше, чем в случае прямоугольной траншеи и круглого силоса.
Рис. 13. Вертикальное давление в щелевом (сплошная линия) и в коническом (линия из точек) бункерах L = R = 3,5 м, а = 30°у = 800 кг / м3, у = 160, ф = 200, х = 300) Fig. 13. Vertical pressure in the slot (solid line) and conical (line of dots) in bins of L = R = 3.5 m, а = 30°, у = 80° kg / m3, у = 16°, ф = 20°, х = 30°)
Из уравнения (16) следует, что также как и в случае щелевого бункера, при а ^ 0 решение (16) переходит в решение (5) для круглого силоса.
Если стенка бункера абсолютно гладкая, то выражение (16) преобразуется к виду: То есть, найдем значение выражения:
Рассчитаем составляющие уравнения (19).
Давление на стены конического бункера определяется выражениями (16):
а - В - К [л Ьг\(л Ь ■ гУ^'1]
[{ Н ){ Н ) \ (20)
Рп сю =
2 А-Ь
где а = 1 - (1 - Ю ■ cosa .
lim pz сгэ =
А—»n
Y ■ g - R b - (1 +
_ (R h z\
n'y R У
Y ■ 9 -R z —= Y ' 9 b-( 1+aî?) R 1+MÏ
При ц^ 0 выражение (16) переходит в формулу гидростатического давления. Посмотрим, выполняется ли условие равновесия сыпучего материала в коническом бункере. Для этого, как обычно, найдем разность А между силой тяжести FT сыпучего материала в бункере и суммой сил вертикальной компоненты силы трения о боковую стенку FTP.в и вертикальной компонентой силы реакции стенок бункера FB .
Сила тяжести FT сыпучего материала в коническом бункере будет определяться как:
где h - мощность пласта (рис. 14).
Рис. 14. Пласт сыпучего материала в коническом бункере Fig. 14. The layer of granular material in a conical hopper
Вертикальная компонента cIFb силы реакции стенки бункера будет равна:
dz
dFB = Р„ ■___ ■ 27Г ■ г ■ sina
п ■
cosa
Значение текущего радиуса определяется выражением: Силы реакции стенки конического бункера, действующая на весь пласт, будет равна:
Произведя интегрирование, получим:
27Г -а-В ■ R3
Fr =
i-(
1 -(1 -b -ху +-,
3 -b
2-Л. 1 - Ь ■ X) ÍJ
3 ■ (2 А-Ъ) I v ' 2 ■ (А + й)
где х - безразмерная переменная: х = h / R.
Соответственно, сила трения Ftp, также, как и для щелевого бункера, будет равна:
FTP = FB-f
В
в
Следовательно, полная сила сопротивления, действующая на пласт, будет равна:
Ктг)
В этом случае, условие (19) запишется в виде:
д= FT -FB -FjP£ = \у -д-п■ R3fi-(1 -Ъ ■ х)
: /1 \\ т-а
1 -(гТз^
B-R' L п , 3-й 1 -(1 -Ьх) +-,
¡■Ь-».^-)Wuf)
г ° - - \Ь ^ 1 и? )) 3-{2 А-Ъ) I 1 ' 2-{А + Ъ) Максимальное значение безразмерной переменной хт = 1 / Ъ. При этом условие (19)
3-й
примет вид:
д= FT - FB - FTP B
i n, i 2тг ■ a ■ В ■ R3
З-к-R' b-JT^TJ
^lYi- 3± Yfi +Щ
-b)\ 2-(A + b)) \ b)
вид:
В равновесном состоянии B^ ygA^ /ла ^ 1 .В этом случае уравнение (19) примет
д= FT
р р -1v « 77 я*1 (л 2'№е + Ь)\г-уе-Ь)\_
Обсуждение результатов. Таким образом, получили обычный результат: в равновесном состоянии сила тяжести в точности равна сумме сил трения о стенки бункера плюс сила реакции со стороны этих стенок. Во всех остальных случаях А > 0.
Однако необходимо учитывать, что это чисто формальный вывод. Действительно, если начнется движение вниз пласта сыпучего материала, то возникнет смещение частиц сыпучего
материала друг относительно друга (ведь бункер суживается к низу). Следовательно, возникнут силы внутреннего, сухого трения, против которых будет совершаться работа за счет силы тяжести, которая не учитывалась в данном выводе.
Для нахождения общего решения, разобьем весь объем конического бункера, также, как и для щелевого, на такие же две области / и II (рис. 15).
Рис. 15. Области бункера Fig. 15. Regions of the hopper
Это необходимо сделать для удовлетворения граничному условию на поверхности бункера и дальнейшего введения зависимости от горизонтальной координаты р.
При этом решение в области I будет иметь такой же вид, как и для силоса в той же области I - (9):
/у сю = В ■ z
Решение в области II представим в виде:
Рис. 16. Вертикальное давление в коническом бункере вдоль осевой и образующей боковую поверхность
линий при различных углах наклона стенки бункера а (у = 800 кг / м3, у = 160, ф = 200, х = 300, R = 3,5 м) Fig. 16. Vertical pressures in a conical hopper along the axis and forming the side surface lines at various angles of inclination of the hopper walls а (y = 800 kg / m3, у = 160, ф = 200, х = 300, R = 3.5 m)
В таком виде решение (21) удовлетворяют граничному условию на поверхности бункера. Решение (21) на границе областей, при р = Яъ переходит в решение (9) и при р = Я^ переходит в решение (16). На рис. 16 приведены графики зависимости вертикального давления
от координаты z для конического бункера с углом наклона 300 и 150. При этом сплошная кривая отображает ход давления по центру бункера, а кривая из точек - по его боковой поверхности. Выделенная темным цветом прямая соответствует решению в области I. Как видно из графиков (рис. 16), угол наклона влияет на величину давления и его характер зависимости от координаты. Кроме того, из рис. 16 следует, что вертикальное давление в центре конического бункера превышает давление на его поверхности.
Исследуем теперь характер расположения изолиний давления в коническом бункере. При этом, во-первых, можно рассмотреть давление действующее на площадку, ори-
ентированную перпендикулярно направлению на точку симметрии конического бункера (точка Ц на рис. 17).
H
ц '^.„Центр симметрии , Center of symmetry
Z
Рис.17. Ориентация давления и угол ß Fig. 1T. The orientation of the pressure Pц and the angle ß
Это направление выбрано из следующих соображений. При истечении сыпучего материала из бункера, осредненные траектории частиц могут представлять собой непересекающиеся прямые направленные в точку симметрии конического бункера. Вследствие этого, можно рассмотреть именно это направление давления. В этом случае именно оно будет создавать форму поверхности движущегося пласта. При этом можно сказать, что изолиния этого давления будет в значительной степени совпадать с формой поверхности движущегося пласта. Если изолиния будет иметь форму свода, то и поверхность пласта будет иметь форму свода, и возможна остановка истечения сыпучего материала.
Соответственно давление, действующее на площадку, ориентированную вдоль направления на точку симметрии бункера, обозначим как.
Формулы для расчета давлений P^{p,z} и запишем в виде:
Pц (p, z) = Pz II (p, z) (1 - (1 - ^ ) sin//; Pn, (p, z) = Pz II (p, z) (1 - (1 - ^ ) cos// (22)
Здесь, согласно рис. 17, тригонометрические функции угла ß определяются выражениями:
Рис. 18. Трехмерный график зависимости давления Рц(р^).(у = 800 кг / м3, у = 160,
ф = 200, х = 300, а = 300, R = 3,5 м) Fig. 18. Three-dimensional graph of pressure RC(p,z). (y = 800 kg / m3, у = 160, ф = 200, х = 300, а = 300, R = 3.5 m)
На рис. 18 представлен трехмерный график давления Рц (p,z) (давления, направленного к центру симметрии) для сыпучего материала в коническом бункере в областях I и II . Как видно, давление в центре превышает давление вдоль стен. Также достаточно видна хорошая «сши-ваемость» решений. Но осредненные траектории движения частиц могут иметь форму прямых, параллельных боковой поверхности бункера - пересекающиеся траектории. В этом случае удобно рассматривать давления параллельное - и перпендикулярное боковой стенке
бункера - Они будут рассчитываться по формулам:
На рис. 19 представлены изолинии давленийPn(ß-Z) и рассчитанные по
формулам (24).
(24)
Рис. 19. Изолинии давления Рц(р^) (а) и Р„д(р) (б) в коническом бункере R = 3,5 м; у = 160;
ф = 20°; х = 300 ; у = 800 кг / м3 ; а = 30°. Fig. 19. Contours of pressure Rc(p,z) (a) and Rnc(p) (b) in a conical hopper R = 3.5 m ; ф = 160 ;
ф = 200 ; x = 300 ; у = 800 kg / m3 ; a = 300
Рис. 20.а.
Давление на дно бункера,r = 2см,
а
Рис. 20.б.
350. Пшеница ф = 200,
кг/м3, Равновесное давление - линия из точек,
V = 160, X = 300, y = 770
008 0.1
Давление на боковую стенку бункера,
r = 2см, а = 350. Пшеница ф = 200 , v = 16° , X
3
стационарное давление - сплошная линия. Экспериментальные значения - выделенные, вертикальные линии.
Fig. 20. a. The pressure on the bottom of the hopper, r = 2cm, a = 350. Wheat ф = 200, ф = 160, x = 300, y = 770 kg/m3, the Equilibrium pressure of a line of points, a stationary pressure - solid line. Experimentalvalues - highlightedverticalline.
30 , у = 770 кг/м , Равновесное давление -линия из точек, стационарное давление сплошная линия. Экспериментальные значения - выделенные, вертикальные линии. Fig. 20.b. The pressure on the side wall of the hopper, r = 2cm, a = 350. Wheat ф = 200 ,ф = 160 , x = 300 , y = 770 kg/m3 , the Equilibrium pressure of a line of dots, stationary pressure solid line. Experimental values - highlighted vertical line.
Мы видим, что в коническом бункере изолинии давления также представляют собой кривые, представляющие собой части окружности. Следовательно, и здесь происходит образование сводчатых структур, тормозящих истечение сыпучего материала из бункера. Результаты сравнения теоретических расчетов с экспериментальными данными, взятыми из [22] для модели бункера малых размеров, представлены на рис. 20.
Как видно из рис.20, в пределах точности экспериментальных данных, а также усредненных параметрах сыпучего материала, согласие можно считать вполне удовлетворительным, гораздо лучшим, чем на рис. 7 - 9. Вероятно, это можно объяснить тем, что в модели бункера малых размеров процессы уплотнения сыпучего материала завершаются за гораздо меньший промежуток времени, чем в реальном бункере. Кроме того, из данных рис. 20 следует, что в этом случае, состояние сыпучего материала ближе к равновесному состоянию, чем к стационарному, что вполне естественно.
Вывод. Из всего изложенного вытекает, в сыпучем материале, находящемся в некоторой ёмкости, образуются статические своды, которые переходят в динамические своды при его истечении. Эти своды тормозят истечение, уменьшая его скорость, тем самым понижая экономическую эффективность различных технологических процессов, имеющих в своём составе бункерные устройства.
Библиографический список:
1. Coulomb. Application des regles de maximisetminimis a quelquesproblemes de statiqurelatifs a J'architecture. Mem des Savants Etrangers de J'dcad d / Coulomb // Sc. de Paris, 1773.
2. Малышев М.В. Об определении угла внутреннего трения и сцепления предельно напряженной сыпучей среды / М.В. Малышев // Изв. АН СССР. Сер. ОТН. - 1954.- № 7. -С. 18-21.
3. Клейн Г.К. Строительная механика сыпучих тел / Г.К. Клейн. - М.: Госсторйиздат, 1956. - 252 с.
4. Соколовский В.В. Статика сыпучей среды / В.В. Соколовский. - М.: Физматгиз, 1960. - 186 с.
5. Соловых С.Ф. О связи напряженного состояния сыпучего тела с формой движения его в силосе / С.Ф. Соловых
// Изв. вузов. Строительство и архитектура. - 1962. - № 5. -С. 22-25.
6. Ехансон. А.И. Поля напряжений и скоростей при гравитационном течении масс / А.И. Ехансон // Прикладная механика. - 1964. - № 3. - С. 149.
7. Гениев Г.А. Об одном варианте течения сыпучей среды / Г.А. Гениев // Строительная механика. - 1965. - № 6. - С. 23-26.
8. Дерсевич Г. И. Механика зернистой среды / Г.И. Дер-севич // Проблемы механики. - М., 1961. - Вып. III. - 368 с.
9. Jenkin C.F. Pressure Exerted by Granular Material an Application of the Principles о!" Dilatancy / C.F. Jenkin // Proceedings pf Royal Society of London. Ser. A. - 1931. - Vol. 131. -P. 53-89.
10. Гячев Л.В. О механической модели сыпучего тела / Л.В. Гячев // Механика сыпучих материалов: тез. докл. Всесоюз. конф. - Одесса, 1975. - С. 3-4.
11. Гячев Л.В. Движение сыпучих материалов в трубах и бункерах / Л.В. Гячев. - М.: Машиностроение, 1968. - 184 с.
12. Keller H. Das Schuttgutmodell von L.V.GjachevzurBeschreibung der Gesetzmabigkeiten der Bun-kerungkohasionsloserSchuttgutter und seine experimentelleU-berprufung / Н. Keller . - Weimer, 1982. - 76 s.
13. Keller H. KorrekturformelnzurTheorie des Schuttgutaus-flussesausBunkern von L.V.Gjachev / Н. Keller // Wissen-schaftliceZeitschrift der Hochschule furArchitectur und Bau-wessen. - Weimar, 1987. - RaieB.h. 5/6. - S. 292-295.
14. Keller H. BeitragezumSchuttgutausflubausBehaltern: Dissertationen / Н. Keller. - Weimar, 1989. - 238 s.
15. Богомягких В.А., Ялтанцев В.Г. Зыков В.А. Процесс образования сводов в силосах и бункерах при истечении сыпучих материалов // Механика деформируемых систем в сельхозмашиностроении. - Ростов - н/Д, 1974. С. 115 -119.
16. Богомягких В.А., Пепчук А.П. Интенсификация разгрузки бункерных устройств в условиях сводообразования зернистых материалов. - Зерноград, 1997, - 161 с.
17. Богомягких В.А., Кунаков В.С. Расчёт параметров и режимов работы сводоразрушающих устройств. // Труды
References:
1. Coulomb C.-A. Essai sur une application des règles de maximis et de minimis à quelques problèmes de Statique relatifs à l'Architecture. Mémoires de l'Académie royale des sciences de Paris (Savants étrangers). Paris; 1773.
2. Malyshev M.V. Ob opredelenii ugla vnutrennego treniya i stsepleniya predel'no napryazhennoi sypuchei sredy. Izv. AN SSSR. Ser. OTN. 1954;7:18-21. [Malyshev M.V. On the determination of the angle of internal friction and the adhesion of an extremely stressed loose medium. Izv. AN SSSR. Ser. OTN. 1954;7:18-21. (In Russ.)]
3. Klein G.K. Stroitel'naya mekhanika sypuchikh tel. M.: Gosstoriizdat; 1956. 252 s. [Klein G.K. Construction mechanics of loose bodies. M.: Gosstoriizdat; 1956. 252 p. (In Russ.)]
4. Sokolovskii V.V. Statika sypuchei sredy. M.: Fizmatgiz; 1960. 186 s. [Sokolovskii V.V. Statistics of the loose medium. M.: Fizmatgiz; 1960. 186 p. (In Russ.)]
5. Solovykh S.F. O svyazi napryazhennogo sostoyaniya sypuchego tela s formoi dvizheniya ego v silose. Izv. vuzov. Stroitel'stvo i arkhitektura. 1962;5:22-25. [Solovykh S.F. On the relationship between the stressed state of a loose body and the form of its movement in a silage. News of higher educational institutions. Construction. 1962;5:22-25. (In Russ.)]
6. Ekhanson A.I. Polya napryazhenii i skorostei pri gravi-tatsionnom techenii mass. Prikladnaya mekhanika. 1964;3:149. [Ekhanson A.I. Fields of stresses and velocities in the gravitational flow of masses. International Applied Mechanics. 1964;3:149. (In Russ.)]
7. Geniev G.A. Ob odnom variante techeniya sypuchei sredy. Stroitel'naya mekhanika. 1965;6:23-26. [Geniev G.A. About a version of the loose medium flow. Structural Mechanics and Analysis of Constructions. 1965;6:23-26. (In Russ)]
8. Dersevich G. I. Mekhanika zernistoi sredy. Problemy mek-haniki. Vyp. III. M.: 1961. 368 s. [Dersevich G. I. Mechanics of a granular medium. Issues of Mechanics. Vol. III. M.: 1961. 368 p. (In Russ.)]
9. Jenkin C.F. Pressure Exerted by Granular Material an Application of the Principles оf Dilatancy. Proceedings of Royal Society of London. Ser. A. 1931;131:53-89.
10. Gyachev L.V. O mekhanicheskoi modeli sypuchego tela. Tezisy dokladov Vsesoyuznoy konferentsii "Mekhanika
УМеждународной конференции по динамике технологических систем. Ростов - н/Д, 1997, с. 85 - 87.
18. Федосеев В.Б., Зацаринная И.А. К определению законов моментов образования и разрушения динамических сводов при установившемся режиме истечения зернистых материалов. // Вестник аграрной науки Дона: теоретический и научно-практический журнал. - Зерноград: ФГОУ ВПО АЧГАА, 2010. - Вып. 3. - С. 64 - 68.
19. Федосеев В.Б., Зацаринная И.А., Титученко А.А. Марковский процесс неустановившегося режима истечения сыпучего материала. Труды Кубанского государственного аграрного университета. - Краснодар: ФГОУ ВПО КубГАУ, 2011. - Вып. 2 (29). С. 188 - 192.
20. Федосеев В.Б., Зацаринная И.А. Стохастический характер образования динамических сводов при установившемся режиме истечения сыпучих материалов из бункеров. Вестник Мичуринского государственного аграрного университета. - Мичуринск: ФГОУ ВПО МичГАУ, 2011. - №1, ч. 1. С. 196 - 199.
21. Мурадова С.Ш., Федосеева Л.В. Воздействие сыпучих материалов на стенки щелевых бункеров. - Вестник Дагестанского государственного технического университета. Технические науки. - том 44, №3, 2017. - С. 24 - 38.
22. Ким В.С. Давление зерна и совершенствование конструкций силосов зерновых элеваторов / В.С. Ким. - М.: Хлебоиздат, 1959. - 55 С.
sypuchikh materialov". Odessa; 1975. S. 3-4. [Gyachev L.V. On the mechanical model of a loose body. Materials of AllUnion conference "Mechanics of loose materials". Odessa; 1975. P. 3-4. (In Russ.)]
11. Gyachev L.V. Dvizhenie sypuchikh materialov v trubakh i bunkerakh. M.: Mashinostroenie; 1968. 184 s. [Gyachev L.V. Motion of loose materials in pipes and hoppers. M.: Mashinostroenie; 1968. 184 p. (In Russ.).
12. Keller H. Das Schuttgutmodell von L.V.Gjachev zur Beschreibung der Gesetzmabigkeiten der Bunkerung kohasionslo-ser Schuttgutter und seine experimentelle Uberprufung. Weimer; 1982. 76 p.
13. Keller H. Korrekturformeln zur Theorie des Schuttgutausflusses aus Bunkern von L.V.Gjachev. Wissenschaftlice Zeitschrift der Hochschule fur Architectur und Bauwessen. Raie B.h. 5/6. Weimar; 1987. P. 292-295.
14. Keller H. Beitrage zum Schuttgutausflub aus Behaltern: Dissertationen. Weimar; 1989. 238 p.
15. Bogomyagkikh V.A., Yaltantsev V.G. Zykov V.A. Protsess obrazovaniya svodov v silosakh i bunkerakh pri istechenii sy-puchikh materialov. Mekhanika deformiruemykh sistem v sel'khozmashinostroenii. Rostov-na-Donu; 1974. S. 115 - 119. [Bogomyagkikh V.A., Yaltantsev V.G. Zykov V.A. The process of formation of arches in silos and hoppers at the expiration of bulk materials. Mechanics of deformable systems in agricultural machinery. Rostov-on-Don; 1974. P. 115 - 119. (In Russ.)]
16. Bogomyagkikh V.A., Pepchuk A.P. Intensifikatsiya razgruzki bunkernykh ustroistv v usloviyakh svodoobrazovani-ya zernistykh materialov. Zernograd;1997. 161 s. [Bogomyagkikh V.A., Pepchuk A.P. Intensification of unloading of hopper devices under the conditions of arch formation in granular materials. Zernograd; 1997. 161 p. (In Russ.)]
17. Bogomyagkikh V.A., Kunakov V.S. Raschet parametrov i rezhimov raboty svodorazrushayushchikh ustroistv. Trudy V Mezhdunarodnoi konferentsii po dinamike tekhnologicheskikh sistem. Rostov-na-Donu; 1997. C. 85 - 87. [Bogomyagkikh V.A., Kunakov V.S. Calculation of parameters and modes of operation of the arch-breaking devices. Proceedings of the V International Conference on the dynamics of technological
systems. Rostov-on-Don; 1997. C. 85 - 87. [(in Russ.)] (In Russ.)]
18. Fedoseev V.B., Zatsarinnaya I.A. K opredeleniyu zakonov momentov obrazovaniya i razrusheniya dinamicheskikh svod-ov pri ustanovivshemsya rezhime istecheniya zernistykh mate-rialov. Vestnik agrarnoi nauki Dona. 2010;3:64 - 68. [Fedoseev V.B., Zatsarinnaya I.A. On the determination of the laws of the formation moments and destruction of dynamic arches under the steady regime of the outflow of granular materials. Don agrarian science bulletin. 2010;3:64 - 68. (In Russ.)]
19. Fedoseev V.B., Zatsarinnaya I.A., Tituchenko A.A. Mar-kovskii protsess neustanovivshegosya rezhima istecheniya sypuchego materiala. Trudy Kubanskogo gosudarstvennogo agrarnogo universiteta. 2011;2(29):188 - 192. [Fedoseev V.B., Zatsarinnaya I.A., Tituchenko A.A. Markov process of the unsteady outflow mode of loose material. Work collection of the Kuban State Agrarian University. 2011;2(29):188 - 192. (In Russ.)]
20. Fedoseev V.B., Zatsarinnaya I.A. Stokhasticheskii kharak-ter obrazovaniya dinamicheskikh svodov pri ustanovivshemsya
rezhime istecheniya sypuchikh materialov iz bunkerov. Vestnik Michurinskogo gosudarstvennogo agrarnogo universiteta. 2011;1:196-199. [Fedoseev V.B., Zatsarinnaya I.A. Stochastic character of the formation of dynamic arches under steady-state conditions for the flow of loose materials from hoppers. Vestnik Michurinskogo gosudarstvennogo agrarnogo universiteta. 2011;1:196-199. (In Russ.)]
21. Muradova S.Sh., Fedoseeva L.V. Vozdeistvie sypuchikh materialov na stenki shchelevykh bunkerov. Vestnik Dage-stanskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhnicheskie nauki. 2017;44(3):24 - 38. [Muradova S.Sh., Fedoseeva L.V. Impact of loose materials on track hopper walls. Herald of Daghestan State Technical University. Technical Sciences. 2017;44(3):24 - 38. (In Russ.)]
22. Kim V.S. Davlenie zerna i sovershenstvovanie konstruktsii silosov zernovykh elevatorov. M.: Khleboizdat; 1959. 55 s. [Kim V.S. The grain pressure and the improvement of silo grain silos. M.: Khleboizdat; 1959. 55 p. (In Russ.)]
Сведения об авторах.
Мурадова Сафура Шиховна - кандидат педагогических наук, доцент. Федосеева Людмила Владимировна - старший преподаватель. Information about the authors.
Safura S. Muradova - Cand. Sci. (Pedagogical), Assoc.Prof. Lyudmila V. Fedoseeva - Senior Lecturer.
Конфликт интересов. Conflict of interest.
Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.
Поступила в редакцию 11.01.2018. Received 11.01.2018.
Принята в печать 20.02.2018. Accepted for publication 20.02.2018.