МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ
DOI: 10.34031/article_5cb1e65b635e41.79569989 1 *Хуртасенко В.А., 1Шрубченко И.В.
белгородский государственный технологический университет им. В.Г. Шухова Россия, 308012, Белгород, ул. Костюкова, д. 46 *E-mail: [email protected]
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДЛЯ ОПТИМИЗАЦИИ ПАРАМЕТРОВ ОБРАБОТКИ ПОВЕРХНОСТЕЙ КАЧЕНИЯ ТЕХНОЛОГИЧЕСКИХ АГРЕГАТОВ
МОБИЛЬНЫМ ОБОРУДОВАНИЕМ
Аннотация. Восстановительная обработка крупногабаритных деталей вращающихся технологических агрегатов, к которым относятся бандажи и ролики обжиговых печей, в настоящее время выполняется с использованием встраиваемых станочных модулей. При этом периодическая обработка ведется в процессе работы агрегата без его остановки. Такая обработка требует особых технических и технологических подходов. Это связано не только с большими габаритами самих деталей, но и с непостоянством положения оси вращающейся детали, установленной на два опорных ролика и имеющей погрешность формы в поперечном сечении. Используемые технологии характеризуются недостаточной изученностью процесса формирования поверхности при механической обработке и произвольно задаваемыми режимами обработки. Актуальной задачей является разработка методик для поиска рациональных режимов, при которых будет обеспечиваться требуемая точность и производительность обработки. Особое место при этом занимают вопросы моделирования обработки для поиска оптимальных технологических режимов с учетом технических и физико-механических особенностей процесса. Разработка адекватных математических моделей при этом требует глубокого анализа функциональных особенностей оборудования и технологических особенностей процесса. Установлено, что для оптимизации режимов обработки необходимо ввести ряд ограничений, связанных с инструментом, приводом подач, жесткостью технологической системы, качеством обработки. При этом в качестве целевой функции используется производительность обработки, а в качестве варьируемых параметров глубина и подача.
Ключевые слова: моделирование обработки, крупногабаритные детали, мобильное оборудова-
ние, технологические параметры, оптимизация.
Введение. В современной индустрии применяются вращающиеся технологические агрегаты, непрерывного производственного цикла. К таким агрегатам относятся вращающиеся печи. Они используются в цементной, металлургической, силикатной и химической, сахарной и пищевой промышленностях, а также для утилизации отходов [1, 2]. В конструкциях вращающихся печей особое внимание уделяется прямолинейности оси вращения, настройка которой во многом зависит от точности формы деталей опор - бандажей и опорных роликов [3]. Так как в процессе эксплуатации вследствие высоких динамических и термических нагрузок происходит потеря точности формы как корпуса печи [4], так и поверхностей качения деталей опор, то возникает необходимость их периодической обработки с целью восстановления требуемой точности. В настоящее время для периодической обработки поверхностей качения в технологических агрегатах по производству строительных материалов, таких как сушильные барабаны,
вращающиеся печи и т.п., широко применяют технологии с использованием мобильного оборудования. Это позволяет восстановить их форму без остановки основного производственного процесса, что существенно влияет на технологичность процесса. При реализации известных схем обработки с использованием существующих образцов оборудования [5, 6, 7, 8, 9] глубина резания обеспечивается настройкой вылета режущего инструмента. Однако данный технологический параметр на различных участках обрабатываемой поверхности может существенно отличаться даже в процессе выполнения одного рабочего хода. В силу таких особенностей используемые технологические процессы, основанные на применении как лезвийной, так и абразивной обработки, не всегда позволяют обеспечивать требуемую точность поверхностей, особенно на поверхностях качения бандажей. Основной причиной здесь является недостаточная изученность процесса формирования поверхностей при обработке и произвольно зада-
ваемые режимы обработки. С учетом этого актуальной задачей является разработка методик поиска рациональных режимов, при которых будет обеспечиваться требуемая точность и производительность обработки. Наиболее применяемым видом мобильного оборудования для лезвийной обработки поверхностей качения бандажей и опорных роликов, является универсальный встраиваемый станок УВС-01М [9],
пН.нВпЖ
Прибод ручного I предельного I перемещения
разработанный на кафедре технологии машиностроения Белгородского государственного технологического университета им. В.Г. Шухова, элементы конструкции которого защищены патентами [5, 6].
Станок имеет опорные стойки, в которые установлена направляющая, несущая продольный и поперечный суппорта и оснащается приводом продольных подач (рис. 1).
А
Разец
Резцедержатель
Сменное техно/югичесхиеУу-
Опорный ро/ш
Рис. 1. Универсальный встраиваемый станок УВС-01М
Поперечный суппорт !
Рукоятка поперечного I су&зор&а
Для установки станка на различные конструкции восстанавливаемых опор, применяют сменные технологические наладки, соответствующие конструктивным особенностям этих опор (рис. 2), и как бы встраивается в саму конструкцию опоры.
Наличие таких сменных технологических наладок позволяет обеспечить возможность применения станка модели УВС-01М на различных типоразмерах опор восстанавливаемых агрегатов, обеспечивая тем самым его универсальность.
Рис. 2. Конструкции сменных технологических наладок для установки станка на корпуса подшипников
Основная часть. Использование такого технологического подхода с использованием мобильных станочных модулей и специальной технологической оснастки позволяет выполнить восстановительную обработку бандажей и роликом практически всех типоразмеров существующих опор. Однако, следует отметить, что различные технологические агрегаты имея широкий спектр типоразмеров, могут осуществлять еще и вращение с различной частотой. Так как
основным движением для условий обработки поверхностей качения будет являться вращение самого агрегата, то очевидно, и линейная скорость на различных типоразмерах агрегатов окажется различной. Например, для ряда типоразмеров вращающихся печей, линейная скорость вращения, а, следовательно, и скорость резания, будет иметь следующие значения: V = 11,62; 12,25; 13,53; 15,23; 17,17; 19,15; 21,35; 26,53 м/мин. Процесс обработки поверхностей, с
учетом их размеров и небольших скоростей, требует значительных затрат времени. Поэтому в качестве целевой функции следует выбрать производительность обработки. А в качестве варьируемых параметров: глубину резания и подачу. Для условий оптимизации режимов обработки поверхностей качения следует ввести ряд технических ограничений: по стойкости режущего инструмента; по допустимой силе резания, определяемой тяговой силой привода подач; по минимальной и максимальной величине продольной подачи по характеристике привода подач; по минимальной и максимальной величине снимаемого за рабочий ход припуска, определяемого условиями формообразования для бесцентровой схемы обработки бандажа; по минимальной и максимальной величине снимаемого за один рабочий ход припуска, определяющего точность обрабатываемой поверхности в продольном сечении; по допустимой прочности режущего инструмента; по жесткости режущего инструмента; по прочности механизма подач; по допустимой шероховатости обрабатываемой поверхности. Далее приведены обоснования и расчетные зависимости для рассматриваемых ограничений.
Ограничение по требуемой стойкости инструмента. Связь варьируемых параметров с периодом стойкости инструмента, можно получить из следующей зависимости [10]:
Су
V = -
Tm t*
■■ kv ,
(1)
где Су = 340 , Ку = 0,83 , ху = 0,15 , Уу = 0,45 т = 0,2 - табличные значения для условий наружного продольного точения инструментом, оснащаемым твердосплавной пластиной; Т -требуемый период стойкости инструмента, мин.
Табличное значение периода стойкости, при скорости резания - 150 м/мин, для токарных резцов составляет 60 мин. При таких условиях, путь резания составит - 9000 м. Для условий обработки бандажей и роликов мобильным оборудованием, при скорости резания - 20 м/мин, чтобы резец прошел путь резания - 9000 м, потребуется 450 мин. С учетом преобразований, получаем техническое ограничение в виде:
^ . ^ < <СУ^ (2)
Тт ■ V
Логарифмируя это выражение получим:
C ■ к
*v ■ ln( t) + yv ■ ln( 5) < ln(—m—-)
Tm ■ V
(3)
И, приняв ln( t) = x1, ln(s) = x2, ln(
C ■ к
^ v ^v Tm ■ V
) = b^
получаем техническое ограничение, приведенное к линейному виду:
XV ■ Х1 + УУ ■ х2 < Ъ1 (4)
Ограничение по допустимой силе резания, определяемой тяговой силой привода подач. Для установления связи допустимой тяговой силы с варьируемыми параметрами, воспользуемся следующей зависимостью [10]:
тх и (5)
Px = 10 ■ Cx ■ tx* ■ sy* ■Vn* ■ кр,
где Сх = 339, кр = 0,74, хх = 1, Ух = 0,5, пх = -0,4 - табличные значения коэффициентов, для условий обработки поверхностей качения.
После преобразований и логарифмирования, получаем:
** ■ ln(t) + Ух ■ ln(s) < ln(
Px
10 ■ C* -Vn* ■ кр
) (6)
Производя соответствующие замены, получаем:
хх • х1 + Ух • х2 < Ь2 (7)
Ограничение по минимальной и максимальной величине продольной подачи зависящее от характеристик привода подач. Это ограничение будет определяться характеристиками самого привода продольных подач, т.е.
х2 < ¿3, х2 > ¿4 (8)
Ограничение по минимальной и максимальной величине припуска, снимаемого за рабочий ход, определяемого условиями формообразования для бесцентровой схемы обработки бандажа. Для условий бесцентровой обработки поверхностей качения бандажей, глубина резания будет оказывать существенное влияние на процесс формообразования [11]. При неправильном задании величины снимаемого припуска, возможно даже дальнейшее искажения профиля, а не его исправление. Поэтому на каждый рабочий ход должна задаваться требуемая величина снимаемого припуска, которую мы можем получить по результатам предварительного моделирования процесса формообразования. Такое ограничение по минимальной и максимальной величине припуска, мы можем задать так же в виде двух неравенств:
х1 < ¿5 , х1 < ¿6 (9)
Ограничение по максимально допустимой величине припуска, снимаемого за один рабочий ход, определяющего точность обрабатываемой поверхности в продольном сечении будет также зависеть от величины отжатий в технологической системе. Для условий технологической системы с применением мобильного оборудования, максимальное влияние на величину отжатий, будет оказывать жесткость элементов самого применяемого станка. Таким образом, для определения максимальной и минимальной глу-
v
v
s
бины резания и подачи, которые обеспечат точность формы в продольном сечении, в пределах допусков, мы можем применить предварительное моделирование напряженно деформированного состояния элементов технологической системы с применением метода конечных элементов в CAE приложениях [12].
Ограничение по минимально и максимально допустимой величине подачи, определяемой характеристиками привода продольных подач станка УВС-01. Данное ограничение также можно представить в виде двух неравенств:
* 2 < ¿7 , * 2 < b8 (10)
Ограничение по прочности режущего инструмента. Это ограничение мы можем получить из зависимости по определению предела прочности материала державки резца, на изгиб:
M ■ к
1У± изг. "'з.и.
W
(11)
где М изг - величина изгибающего момента при вылете резца - 1в.р. =150мм, Мизг =р '1в.р.; кз п = 2,4 - коэффициент запаса прочности; W
3
- момент сопротивления державки резца, мм (для прямоугольного сечения державки резца,
шириной Вд = 25 мм и высотой Нд = 40 мм,
Вд • Н2
W =
6
-). После соответствующих подста-
новок и преобразований, получаем:
240 ■ Вд ■ Нд
t*z ■ syz <-
60 ■ C ■ Vnz ■ к ■ l ■ к
W V !\,z Ig.p. !^з.п.
(12)
Логарифмируя и производя соответствующие замены, получаем следующее техническое ограничение:
х • х+У. х < ¿9 (13)
Ограничение по жесткости инструмента. Это ограничение устанавливает взаимосвязь между глубиной резания, подачей и жесткостью режущего инструмента. Максимальная нагрузка, допускаемая жесткостью резца, может быть определена по известной зависимости [13]:
3 • f • Е. 1М
Рж
ж.доп.
l
(13)
е. р.
где f = 0,1 мм - допустимая стрела прогиба для чернового точения, мм; е = (2...2,5) 105МПа ; I -
4
момент инерции сечения державки резца, мм (для прямоугольного сечения державки шири-
Вд = 25 мм и высотой Нд = 40 мм
ной
I =
д
вд ■ Нд
12
Из условия соотношения окружной составляющей Рz и максимальной нагрузки, допускаемой жесткостью резца, имеем:
Pz < Рж .доп., (14)
• ¡¡У, < —- ■/--о "О , (15)
3 ■ f ■ E ■Вд ■ Нд 120 ■ Cz ■ Vnz ■ к, ■ ер.
х2 • у • х2 < ¿10 (16)
Ограничение по допустимой шероховатости поверхности. Это ограничение устанавливает взаимосвязь между варьируемыми параметрами и шероховатостью поверхности [10]. Ожидаемую шероховатость поверхности мы можем вычислить по следующей зависимости:
Ra =
„ 0,31 0,58 0,4 0,4
0,85 ■ t ■ s ■ q> ■ ф1
т,0,06 0,65 TTr> 0,05
V ■ r ■ HB
(17)
где р = 45°, р1 = 45 ° , Г = 2 мм - параметры геометрии режущей части инструмента; НВ = 145 - твердость обрабатываемой поверхности.
После соответствующих преобразований выражения, с учетом обеспечения требуемого параметра шероховатости, получаем:
.0,31 0,58 Ra•V0,06 • Г°,65 • НВ0,05 0,85 -р0'4 .р^4 После соответствующей замены, получим:
0,31-х[+0,58. хг < Ъ11 (19)
Представление полученных ограничений в виде системы неравенств, в совокупности с оценочной функцией, дает математическую модель процесса лезвийной обработки бандажей и роликов вращающихся печей.
' ху • х1 + Уу • х2 < Ъ1 хх • х1 + Ух • х2 < ¿2 х2 < ¿3 х2 > ¿4 х2 > ¿5 х2 > ¿6
х2 > ¿7 (20)
х2 > ¿8 Xz • х1 + Уz • х2 < ¿9 xz • х1 + Уz • х2 < ¿10
к4 • х1 + к3 • х2 < ¿11
_ = (х1 + х2 )тах Задача определения оптимальных параметров технологической системы сводится к отысканию среди всех возможных неотрицательных значений х1 и системы таких значений
).
х1опт и х2опт , при которых линейная функция принимает максимальное значение (./0max)- После определения координат х1опт и х2опт можно вычислить и оптимальные значения варьируемых параметров:
tопт = еХ1опт ,(21)
с — еХ2опт
лопт~е (22)
Выводы: Полученная модель, позволяет определять оптимальные значения технологических параметров для условий бесцентровой лезвийной обработки крупногабаритных деталей, базирующихся на два опорных ролика непосредственно на работающих агрегатах, на основе применения встраиваемых станочных модулей. Используемые в модели зависимости учитывают основные значимые факторы, которые будут оказывать наиболее существенное влияние на качество поверхности для рассматриваемых условий обработки. Применение такой модели возможно, как на этапе проектирования элементов обрабатывающего оборудования и технологической оснастки, так и на этапе реализации восстановительной обработки для настройки технологических параметров оборудования и управления технологическими режимами при многопроходной обработке. Принципиальный подход к построению модели позволяет адаптировать её также для условий других видов обработки резанием, например, для шлифования, применяемого на чистовых операциях восстановительной обработки.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Vijayan S.N., Sendhilkumar S. Industrial Applications of Rotary Kiln in Various Sectors - A Review // International Journal of Engineering Innovation & Research. 2014. Vol. 3. Pp. 342-345.
2. Boateng A.A. Rotary Kilns. Elsevier Inc. Publ., 2015, 390 p.
3. Mogilny S., Sholomitskii A. Precision Analysis of Geometric Parameters for Rotating Machines during Cold Alignmen // Procedia Engineering. Vol. 206. 2017. Pp. 1709-1715.
4. Ramanenka D., Stjernberg J., Jonsen P.. FEM investigation of global mechanisms affecting
brick lining stability in a rotary kiln in cold state // Engineering Failure Analysis. 2016. Vol. 59. Pp. 554-569.
5. Шрубченко И.В, Мурыгина Л.В., Рыбалко В.Ю. Станок для обработки бандажей. Пат. № 125499, Российская Федерация, опубл. 10.03.2013.
6. Санин С.Н., Бондаренко В.Н., Погонин А.А. Следящий суппорт. Пат. № 58420, Российская Федерация, МПК B23Q1/76 (2006.01). опубл. 27.11.2006.
7. Шрубченко И.В., Мурыгина Л.В., Рыбалко В.Ю., Черняев А.С. Станок для обработки бандажей и роликов. Пат. № 118235, Российская Федерация, МПК B23B 5/00 (2006.01), опубл. 20.07.2012.
8. Санин С.Н., Оникиенко Д.А. Разработка концепции мобильного стенда для механической обработки бандажей вращающихся печей с базированием по торцовой поверхности и отверстию // Вестник БГТУ им. В.Г. Шухова. 2016. №2. С. 104-109.
9. Универсальный встраиваемый станок УВС-01 / Федеральный каталог высокотехнологичного оборудования и объектов научного потенциала России. [Электронный ресурс]. Режим доступа: URL: https://каталог-нп.рф/project/281.
10. Справочник технолога-машиностроителя. В 2-х т. Т. 1. Под ред. А.М. Дальского, А.Г. Косиловой, Р.К. Мещерякова, А.Г. Суслова. 5-е изд., перераб. и доп. М.: Машиностроение^, 2001. 912 с.
11.Тимофеев С.П., Хуртасенко В.А., Шруб-ченко И.В. Расчет глубины резания при обработке поверхностей крупногабаритных деталей с нестационарной осью вращения // Вестник БГТУ им. В.Г. Шухова. 2018. № 2. С. 68-74.
12. Гончаров М.С., Шрубченко И.В., Мурыгина Л.В., Щетинин Н.А. Моделирование технологической системы мобильного оборудования // Наукоемкие технологии инновации (XXI Научные чтения): Междунар. науч.-практ. конф., (Белгород, 9-10 октября 2014 г.) Изд-во БГТУ, 2014, Ч.4. С. 223-228
13. Справочник конструктора-инструментальщика: Под общ. ред. В.И. Баран-чикова. М.: Машиностроение, 1994. 560 с.
Информация об авторах
Шрубченко Иван Васильевич, доктор технических наук, профессор кафедры технологии машиностроения. E-mail: [email protected]. Белгородский государственный технологический университет им. В.Г. Шухова. Россия, 308012, Белгород, ул. Костюкова, д. 46.
Хуртасенко Владислав Андреевич, аспирант кафедры технологии машиностроения. E-mail: [email protected]. Белгородский государственный технологический университет им. В.Г. Шухова. Россия, 308012, Белгород, ул. Костюкова, д. 46.
Поступила в феврале 2019 г. © Хуртасенко В.А., Шрубченко И.В., 2019
}'*Khurtasenko V.A., 1Shrubchenko I.V.
1Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.
*E-mail: [email protected]
MATHEMATICAL MODEL FOR OPTIMIZATION OF PROCESSING PARAMETERS OF ROLLING SURFACE OF TECHNOLOGICAL UNITS BY MOBILE EQUIPMENT
Abstract. Restorative processing of large-sized parts of rotating technological units, which include bandages and kiln rollers, is currently performed using embedded machine modules. At the same time, periodic processing is carried out during the operation of the unit without stopping. This processing requires special technical and technological approaches. This is due to the large dimensions of the parts and the inconstancy of the position of the axis of the rotating part, installed on two support rollers and having an error in shape in cross section. The technologies used are characterized by insufficient study of surface formation process during machining and arbitrarily set processing modes. An urgent task is to develop methods for finding rational modes in which the required accuracy and processing performance will be provided. A special place is occupied by the issues of processing modeling to find the optimal technological modes, taking into account the technical and physico-mechanical features of the process. The development of adequate mathematical models requires a deep analysis of the functional features of the equipment and technological features of the process. It is established that for optimization ofprocessing modes it is necessary to enter a number of restrictions related to the tool, the driving feed, the rigidity of technological system, the quality of processing. In this case, the processing capacity is used as the target function; the depth and feed are used as varied parameters.
Keywords: processing modeling, large-sized parts, mobile equipment, technological parameters, optimization.
REFERENCES
1. Vijayan S.N., Sendhilkumar S. Industrial Applications of Rotary Kiln in Various Sectors - A Review. International Journal of Engineering Innovation & Research. 2014. Vol. 3. Pp. 342-345.
2. Boateng A.A. Rotary Kilns. Elsevier Inc. Publ., 2015, 390 p.
3. Mogilny S., Sholomitskii A. Precision Analysis of Geometric Parameters for Rotating Machines during Cold Alignment. Procedia Engineering. Vol. 206. 2017. Pp. 1709-1715.
4. Ramanenka D., Stjernberg J., Jonsen P. FEM investigation of global mechanisms affecting brick lining stability in a rotary kiln in cold state. Engineering Failure Analysis. 2016. Vol. 59. Pp. 554-569.
5. Shrubchenko I.V., Murygina L.V., Rybalko V.U. The machine for processing of bandages. Patent RF, № 125499, 2013. (rus)
6. Sanin S.N., Bondarenko V.N., Pogonin A.A. The watching support. Patent RF, № 58420, 2006. (rus)
7. Shrubchenko I.V., Murygina L.V., Rybalko V.U., Chernyaev A.S. The machine for processing of bandages and rollers. Patent RF, № 118235, 2012. (rus)
8. Sanin S.N., Onikienko D.A. Development of the concept of a mobile stand for the mechanical processing of rotary kiln bands with basing on the end surface and the hole [Razrabotka koncepcii mo-
bil'nogo stenda dlya mekhanicheskoj obrabotki bandazhej vrashchayushchihsya pechej s baziro-vaniem po torcovoj poverhnosti i otverstiyu]. Bulletin of BSTU named after V.G. Shukhov. 2016. No. 2. Pp. 104-109. (rus)
9. Universal Embedded Machine UEM-01 / Federal catalog of high-tech equipment and objects of scientific potential of Russia [Universal'nyj vstraivaemyj stanok UVS-01 / Federal'nyj katalog vysokotekhnologichnogo oborudovaniya i ob"ektov nauchnogo potenciala Rossii]. Access mode: URL: https://KaTa^or-Hn.p$/project/281 .(rus)
10.Reference technologist-mechanical engineer [Spravochnik tekhnologa-mashinostroitelya]. T. 1. Ed. A.M. Dalsky, A.G. Kosilova, R.K. Mesh-cheryakov, A.G. Suslov. Moscow: Mashi-nostroenie-1 Publ., 2001, 912 p. (rus).
11.Timofeev S.P., Khurtasenko V.A., Shrubchenko I.V. Calculation of the depth of cut when machining surfaces of large parts with non-stationary axis of rotation [Raschet glubiny rezaniya pri obrabotke poverhnostej krupnogabaritnyh de-talej s nestacionarnoj os 'yu vrashcheniya]. Bulletin of BSTU named after V.G. Shukhov. 2018. No. 2. Pp. 68-74. (rus)
12.Goncharov M.S., Shrubchenko I.V., Murygina L.V., Schetinin N.A. Simulation of the technological system of mobile equipment [Mod-elirovanie tekhnologicheskoj sistemy mobil'nogo oborudovaniya]. Naukoemkie tekhnologii innovacii
(XXI Nauchnye chteniya): Mezhdunar. nauch.-prakt. konf., (Belgorod, 9-10 oktyabrya 2014 g.) Pp. 223-228. (rus)
13.Reference Tool Designer. Ed. V.I. Baranchikov [Spravochnik konstruktora-instrumental'shchika]. Moscow: Mashinostroenie Publ., 1994, 560 p. (rus)
Information about the authors
Shrubchenko, Ivan V. DSc, Professor. E-mail: [email protected]. Belgorod State Technological University named after V.G. Shukhov.Russia, 308012, Belgorod, st. Kostyukova, 46.
Khurtasenko, Vladislav A. Postgraduate student. E-mail: [email protected]. Belgorod State Technological University named after V.G. Shukhov. Russia, 308012, Belgorod, st. Kostyukova, 46.
Received in February 2019 Для цитирования:
Хуртасенко В.А., Шрубченко И.В. Математическая модель для оптимизации параметров обработки поверхностей качения технологических агрегатов мобильным оборудованием // Вестник БГТУ им. В.Г. Шухова. 2019. № 4. С. 144-150. DOI: 10.34031/article_5cb1e65b635e41.79569989
For citation:
Khurtasenko V.A., Shrubchenko I.V. Mathematical model for optimization of processing parameters of rolling surface of technological units by mobile equipment. Bulletin of BSTU named after V.G. Shukhov. 2019. No. 4. Pp. 144-150. DOI: 10.34031/article_5cb1e65b635e41.79569989