CC BY
20
УДК 631.333.53 05.00.00 Технические науки
ОБОСНОВАНИЕ КОНСТРУКТИВНЫХ ПАРАМЕТРОВ АКТИВНОГО КОЛЬЦА ЦЕНТРОБЕЖНОГО ДИСКА С ПОДАЧЕЙ МАТЕРИАЛА ВДОЛЬ ЛОПАТОК
Припоров Евгений Владимирович к.т.н., доцент.
БРШ-код автора: 9965-6360
Павлов Сергей Николаевич студент
ФГБОУ ВО Кубанский государственный аграрный университет имени И. Т. Трубилина г. Краснодар, Россия
Разработана конструкция однодискового центробежного аппарата с подачей сыпучего материала вдоль лопаток. Рассеивающий диск состоит из лопаток на верхней полки, которых установлен плоский диск с воздухозаборными окнами (дефлектор). Назначение дефлектора - устранить отражение частиц в момент контакта с лопаткой и создать попутный воздушный поток для увеличения дальности полета частиц. Анализ конструктивной схемы рассеивающего диска показал, что на объем поступающего воздуха внутрь рассеивающего диска влияет ширина окна и угол наклона корпуса дефлектора. Определена зависимость ширины окна от угла наклона корпуса дефлектора и высота. На основе анализа конструкции установлены пределы изменения ширины окна и угла наклона корпуса дефлектора. Угол наклона корпуса дефлектора изменяется в интервале 9°...320 Увеличение угла наклона корпуса дефлектора приводит к пропорциональному уменьшению ширины окна. Установлены интервалы варьирования факторов. Выбран симметричный композиционный план второго порядка и составлена матрица планирования эксперимента. В качестве функции отклика принят коэффициент относительной кучности. Реализация эксперимента позволила определить значение коэффициентов регрессии. Адекватность уравнения регрессии полученному уравнению проводилась по критерию Фишера. Установлено, что полученное уравнение адекватно описывает технологический процесс. Для анализа полученного уравнения привели к каноническому виду путем переноса начала координат в новую точку и поворотом на некоторый угол -14,950. Анализ канонического уравнения позволил определить оптимальные значения исследуемых факторов. Ширина окна дефлектора составляет 22.25 мм, угол наклона корпуса дефлектора 8.120
UDC 631.333.53 Technical sciences
JUSTIFICATION OF DESIGN PARAMETERS OF THE ACTIVE RING OF A CENTRIFUGAL DISC WITH THE FILING OF MATERIAL ALONG BLADES
Priporov Evgeny Vladimirovich Cand.Tech.Sci., associate professor SPIN-code: 9965-6360
Pavlov Sergey Nikolaevich student
FGBOU VPO Kuban State Agrarian University named after I. T. Trubilin, Krasnodar, Russia
There has been developed the design of the single disc centrifugal apparatus with the feed of the granular material along the blades. The scattering disc consists of blades, on the top shelf there is a flat disc with air intake windows (deflector). The purpose of the deflector is to eliminate the reflection of particles at the moment of contact with the blade and create a passing air flow to increase the range of flight of the particles. The analysis of the constructive scheme of the scattering disk showed that the amount of air coming inside the scattering disk is affected by the width of the window and the tilt angle of the deflector body. There was determined the dependence of the window width from the angle of the deflector body and height. Based on the analysis of the design there have been determined the limits of setting the window width and the tilt angle of the deflector body. The angle of the deflector body changes in the rate of 90 to 320 . Increasing the tilt angle of the deflector body leads to a proportional reduction of the window width. There have been set the intervals of variation of factors. There have been chosen the symmetrical composite plan of the second order and have been made a matrix of the experimental design. As a function of the response factor there have been taken a relative accuracy coefficient. The implementation of the experiment allowed us to determine the value of the regression coefficients. The adequacy of the regression equation to the equation obtained was performed by the Fisher test. It was determined that the equation adequately describes the technological process. For the analysis of this equation it was turned to the canonical form by moving the origin to a new point and turning at a certain angle of 14,950. The analysis of the canonical equations allowed the authors to determine the optimal values of the studied factors. The window width of the deflector is 22 - 25 mm, the angle of the body of the deflector is 8 - 120
Ключевые слова: РАССЕИВАЮЩИЙ ДИСК, ДЕ- Keywords: SPREADING DISCS, DEFLECTOR, ФЛЕКТОР, ШИРИНА ОКНА, УГОЛ НАКЛОНА WINDOW WIDTH, TILT ENCLOSURE, COMPO-КОРПУСА, КОМПОЗИЦИОННЫЙ ПЛАН, КА- SITE DESIGN, CANONICAL EQUATION, ELLIP-НОНИЧЕСКОЕ УРАВНЕНИЕ, ЭЛЛИПСОИД, ОП- SOID, OPTIMAL VALUE ТИМАЛЬНЫЕ ЗНАЧЕНИЯ
DOI: 10.21515/1990-4665-123-113
На кафедре эксплуатации МТП Кубанского ГАУ проводились работы по улучшению технологических параметров центробежного аппарата с подачей материала вдоль лопаток. Предложено устройство, включающее горизонтальный центробежный диск с радиальными лопатками каждая из которых имеет торцовый отгиб [1,2,3]. На верхней кромки лопаток смонтировано кольцо с воздухозаборными окнами [4,5] . Активное кольцо имеет две основные задачи - устранить отражение частиц от плоскости и второе - создать попутный воздушный поток в направлении движения частиц.
Наличие активного кольца на верхней технологической полке лопатки создает замкнутое пространство внутри диска. Во время вращения диска лопатки выполняют функции центробежного вентилятора, создавая разряжение по оси диска и избыточное давление на выходе. Образовавшееся разряжение в междисковом пространстве нарушает непрерывность подачи материала лопаткой и повышению неравномерности распределения материала рассеивающим диском.
Образующиеся воздушные потоки, во время вращения диска, оказывают влияние на качество рассева. Для нейтрализации этих воздушных потоков требуется обеспечить подачу воздуха внутрь рассеивающего диска. На активном диске над каждой последующей лопаткой в тыльной ее части по всей длине располагаются воздухозаборные карманы навстречу потоку. Эти карманы (дефлекторы) обеспечивают подачу воздуха в рабочую зону каждой предыдущей лопатки. Такая схема подачи воздуха позволяет снизить величину разряжения по оси кольца. Величина подачи воздуха через дефлектор зависит от геометрических размеров воздухозаборного окна
(ширина и длина окна) и угла наклона корпуса к плоскости активного кольца.
На рисунке 1 представлена схема плоского рассеивающего диска с установленным дефлектором для подачи воздуха.
Рисунок1 - Схема диска с дефлектором активного кольца: 1 - активное кольцо; 2 - воздухозаборный карман;
3 - рассеивающий диск (лопатки условно не показаны)
Пусть объемная подача воздуха лопаткой определяется по известной формуле
Qл = РЛУ, (1)
где Qл -объемная подача, м/с;
¥л - площадь рабочей части лопатки, м ;
V - линейная скорость, м/С.
Объем поступающего через дефлектор воздуха зависит от массовой расходной концентрации воздуха, линейной скорости и ряда других факторов. Величина подачи определяется по выражению
Qq = (2)
где Qq - объемная подача воздуха через дефлектор, м /с;
^ - площадь окна дефлектора, м2 ; ¡л - массовая расходная концентрация, ¡л = 3.. .6.
Основное назначение дефлектора - устранить «вентиляторный эффект» по оси диска. Исключить образование разряжения по оси диска обеспечивается при условии равенства м поступающего и исходящего воздуха из
междискового пространства. Учитывая выражение 2 равенство поступающего и исходящего воздуха составляет
Н • I• vV = Нч • ^ , (3)
где Н, НС1 - высота лопатки и высота корпуса дефлектора соответствен но,
м;
I, ^ - длина лопатки и длина дефлектора соответственно, м ;
Учитывая, что равенство длины лопаток и длины дефлектора
I = ¡г
Тогда из рисунка 1 следует, что
\ = вq • tgaq , (4)
где в(1 - ширина окна дефлектора;
ач - угол наклона корпуса дефлектора.
Проведя несложное преобразование, определяем ширину окна дефлектора. Принимаем высоту лопатки пропорциональной ширине окна
в.Т^- (5)
Из выражения (5) установим зависимость ширины окна от высоты лопатки и угла наклона корпуса дефлектора.
Пусть отношение
Н
— = а
вя
Значение коэффициента пропорциональности а = 0,5 .1,0.
tgа = . (6)
Из выражения 6 следует, что угол наклона корпуса дефлектора к горизонтали составляет ач = 9...320. При этом увеличение ширины окна приводит к пропорциональному уменьшению угла наклона корпуса дефлектора.
На основе изучения воздушных потоков нами обоснована необходимость постановки дефлекторов на активном кольце. Дефлектор обеспечивает подачу воздуха в пространство перед лопатками. Конструктивные
размеры дефлектора, как установлено теоретически, зависят от высоты лопатки и, в конечном итоге, определяют объем поступающего воздуха. Однако из-за неполного знания механизма движения воздушных потоков внутри диска описать аналитически процесс подачи воздуха через дефлектор оказалось невозможным. Причина - действие различных факторов. Среди них - частота вращения диска, температура и влажность воздуха, размеры частиц сыпучего материала и многие другие. Решение данной задачи, как известно, возможно методом экстремальных экспериментов. На этапе поисковых исследований определялись уровни варьирования факторов с учетом аналитических зависимостей, полученных ранее. За параметр оптимизации принят коэффициент относительной кучности ук [2].
— Lmax Lmin
к L '
max
где Lmax , Lmin - максимальная и минимальная дальность полета частиц, м.
Уменьшение этого показателя свидетельствует о повышение равномерности распределения сыпучего материала.
Основные факторы, влияющие на кучность распределения сыпучего материала - угол наклона корпуса дефлектора и ширина окна. Выбраны уровни варьирования и интервал варьирования и представлен в таблице 1.
Таблица 1- Уровни варьирования факторов относительной кучности
Факторы Угол наклона корпуса дефлектора к горизонтали, х^, град Ширина окна дефлектора, x2, мм
Основной 22,0 30,0
уровень х0
Верхний уровень хег=+1 32,0 40,0
Нижний уровень хн=-1 12,0 20,0
Интервал варьирования 10,0 10,0
Принят симметричный композиционный план типа Вк состоящий из
2
ядра плана полно факторного эксперимента 2 , звездных точек. Матрица планирования эксперимента представлена в таблице 2.
Таблица 2- Матрица планирования трех уровневый Вк плана второго порядка
№ Фактор Параметр Примечание
п/п Х1 Х2 оптимизации У
1 +1 +1 0,87 ПФЭ 22
2 +1 -1 0,67 ядро плана
3 +1 0,69
4 -1 0,60
5 +1 0 0,85 Звездные точки
6 0 0,73
7 0 +1 0,87 Звездные точки
8 0 -1 0,71
9 0 0 0,88 Нулевые точки
10 0 0 0,88
Х1 Х2
11 0 0 0,87
12 0 0 0,83
Перевод натуральных значений факторов в кодированное проводили по известной формуле [6]
_ Хг Хй
Хг
ЛХг
где хг - кодированное значение фактора; хвг - натуральное значение фактора; хо - натуральное значение фактора на нулевом уровне; хг - интервал варьирования факторов. Коэффициенты регрессии определяем по формулам [7]
Ьо = СЕА - с2ЕЕ
К N
2
Ь= С3 Е XX
и ' ^6,
п=1 п=1 п=1
(7)
(8)
(9)
Ъ,=СЕ ади (10)
п=1
N K N N /1 1\
Ь = C5 Ех 2шУи + с6 Е Е^ X - с и (11)
в
п=1
п=1
где С1, С2, С3, С4, С5 - вспомогательные константы; К - число факторов; N - число опытов. После реализации плана определялись коэффициента регрессии. Вспомогательные константы С; , входящие в формулы для рассматриваемого плана представлены в таблице 3. Таблица 3 - Вспомогательные константы для плана Вк [7]
bii Вспомогательные константы
С1 С2 С3 С4 С5 С6 С7 С8 С9
2 1,250 0,750 0,1667 0,250 0,250 1,118 0,408 0,500 0,660
Планы Вк неортогональны и для оценки коэффициентов регрессии Ь0 и Ьн ковариации соуЬ0Ьи и ковариации соуЬнЬу не равны нулю. Поэтому после расчетов доверительных интервалов и проверки их статистической значимости проводилось исключение незначимых коэффициентов.
Дисперсии, среднеквадратичные ошибки и ковариации определяли по формулам [7]
о 2 _О 2. о _Q . О 2 _о 2 . о 2 _о
S bo = C1Sy ; Sb = C7Sy;S bi = C3Sy ;S bi = C5Sy,
S 2boij = C4Sy2'; Sbij = C2Sy ;COVbobii =-C22Sy
(12)
Проведя расчеты, получили нижеследующие значения величины коэффициентов регрессии
bo=0,874; b1= 0,0614; b2=0,075. b12=0,0275; b11=-0,0804; b22=-0,0877. Дисперсия ошибки опыта [5]
NM _
X X (Yu - Y)2
S2 = -, (13)
y N(m -1)
где m - число повторностей опыта.
'у
Дисперсия ошибки опыта составила Sy = 0,00043. Адекватность модели второго порядка проводили с помощью F - критерия Фишера [6].
Р = , (14)
У
где - дисперсия неадекватности математической модели.
м _
т(Е А - А )2 ^ = ^ . (15)
N - (к +1)
Число степеней свободы числителя и знаменателя /2 определяется по формулам [6]
/ =К - (к + 1), (16)
/ =К(Ш - 1).
Табличное значение Б 0,о5=3,0 при/ =3 /2 =6 [2].
Уравнение регрессии имеет вид
22
У = 0,8744 + 0,0614^7 + 0,075х2 + 0,0275х7х2 - 0,084х7 + 08776х22. (17)
Для анализа и систематизации полученного уравнения необходимо привести к типовой канонической форме. Для преобразования к каноническому виду проведен перенос начала координат в новую точку и поворот старых осей на некоторый угол а в факторном пространстве, в результате чего исчезают линейные части и изменяется значение свободного члена.
Для определения координат нового центра необходимо взять частные производные по независимым переменным и приравнять их к нулю.
^ = 0,0614 + 0,0275х2 - 0,1608х7 = 0
дх
^ = 0,075 + 0,0275х7 - 0,1755х2 = 0
дх2
Координаты нового центра х1з=0,0698, х2з= 0,0983.
Значения отклика в новом центре координат У1з=0,900. Проверка правильности вычислений осуществляется сравнением сумм коэффициента уравнения регрессии при квадратичных членах [6].
Угол поворота осей определяется по формуле [6]
Ви
tg 2а =-—
В11 - В12 tg 2а = -0,267
а = -14,95°
Так как величина угла отрицательна, то оси необходимо повернуть почасовой стрелке на угол а = -14,95°
Уравнение регрессии в каноническом виде
Ук = 0,92 - 0,0698х12 - 0,0983х22
Поскольку коэффициенты имеют одинаковые знаки, поверхность отклика - эллипсоид, а ее центр - экстремум, причем, максимум, так как канонические коэффициенты отрицательны.
пилона, Гвм
Ii
S
20 . 25 30 35 Ширина. ми
Рисунок 2 - Двумерные сечения поверхности отклика
На основании полученных данных построено двумерное сечение поверхности отклика рисунок 2, где представлена система контурных кривых показателя относительной кучности, представляющих эллипсы. Анализ кривых позволил установить оптимальные конструктивные параметры дефлектора. Область оптимума находится в пределах х1= 22...25 мм, угол наклона корпуса х2 = 8...12°. Дальнейшее увеличение конструктивных параметров сопровождается уменьшением показателя относительной кучности.
Для комплектования посевного агрегата следует использовать методику разработанную авторами [8,9]. Снижение расхода минеральных удобрений достигается использованием грубых кормов в качестве мульчи. Для измельчения грубых кормов разработана конструкция измельчителя [10].
Авторами разработана технология и подобраны технические средства для приготовления концентрированных кормов [11].
На урожайность сельскохозяйственных культур влияет качество посевного материала. Повысить качество семенного материала , как установлено авторами, обеспечивается за счет выбора рациональной технологии послеуборочной обработки семян, а параметры семяочистительных машин настраиваются с учетом скорости витания материала [12,13]. Авторами установлено, что включение в состав технологической линии фотосепараторов позволяет достигнуть степени очистки 99,9 %.
Список литературы
1. Устройство для поверхностного рассева минеральных удобрений и других сыпучих материалов Якимов Ю.И., Иванов В.П., Припоров Е.В., Заярский В.П., Волков Г.И., Селивановский О.Б. патент на изобретение RUS 2177216 14.03.2000.
2. Центробежный аппарат с подачей материала вдоль лопаток Припоров Е.В., Кар-тохин С.Н. // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2015. № 112. С. 1499-1511.
3. Центробежный рабочий орган для рассева сыпучего материала Якимов Ю.И., Припоров Е.В., Иванов В.П., Заярский В.П., Волков Г.И., Селивановский О.Б. патент на изобретение RUS 2177217. заявл.14.03.2000.
4. Патент на изобретение RUS 2197807. Центробежный разбрасыватель сыпучих материалов. Якимов Ю.И., Припоров Е.В., Заярский В.П., Волков Г.И., Селивановский О.Б. заявл. 20.04.2001.
5. Прибор для исследования центробежных аппаратов разбрасывателей сыпучих материалов Якимов Ю.И., Припоров Е.В., Карабаницкий А.П., Ткаченко В.Т., Якушев А.А. патент на изобретение RUS 2201059. заявл. 20.04.2001.
6. Адлер Ю.П., Марков Е.В., Грановский Ю.В. Планирование эксперимента при поиске оптимальных условий.- М.: Наука, 1971 - С. 283.
7. Новиков Ф.С., Арсов Я.Б. Оптимизация процессов технологии металлов методами планирования экспериментов. - М.: Машиностроение, 1980, 304с.
8. Сошники зерновых сеялок ресурсосберегающих технологий Припоров Е.В. // В сборнике: Связь теории и практики научных исследований Сборник статей Международной научно-практической конференции. Ответственный редактор: Сукиасян Асатур Альбертович. 2016. С. 63-66.
9. Определение энергосберегающего режима работы тягового агрегата Припоров Е.В. // Инновации в сельском хозяйстве. 2015. № 5 (15). С. 92-95.
10. Измельчитель грубых кормов Маслов Г.Г., Артемов В.Е., Припоров Е.В., Небавский В. А. патент на изобретение RUS 2222175. Заявл. 11.06.2002
11. Технологии и технические средства приготовления концентрированных кормов с использованием соевого белка. Фролов В.Ю., Сысоев Д.П., Припоров И.Е., Горб С.С.// International Scientific and Practical Conference World science. 2016. Т. 1. № 3 (7). С. 53-58.
12. Припоров И.Е., Лазебных Д.В. Рациональная технология послеуборочной обработки семян подсолнечника // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. 2015. № 112. С. 1475-1485.
13. Влияние толщины, ширины и индивидуальной массы семян подсолнечника на скорость их витания Шафоростов В.Д., Припоров И.Е. // Масличные культуры. Научно-технический бюллетень Всероссийского научно-исследовательского института масличных культур. 2010. № 1 (142-143). С. 7б-80.
14. Обоснование применения оптического фотоэлектронного сепаратора в составе универсального семяочистительного комплекса. Припоров И.Е. // В сборнике: Конкурентная способность отечественных гибридов, сортов и технологии возделывания масличных культур Сборник материалов 8-й международной конференции молодых учёных и специалистов. 2015. С. 138-141.
15. Классификация оптических фотосепараторов для сортирования семян подсолнечника. Припоров И. Е., Шафоростов В. Д. // Вестник Казанского государственного аграрного университета. 2015. Т. 10. № 1. С. б8-70.
16. Моделирование процесса сепарирования семян подсолнечника в вертикальном пневмоканале ветро-решетных зерноочистительных машин Шафоростов В.Д., Припоров И.Е.// Масличные культуры. Научно-технический бюллетень Всероссийского научно-исследовательского института масличных культур. 2011. № 1 (14б-147). С. 113-118.
References
1. Ustrojstvo dlja poverhnostnogo rasseva mineral'nyh udobrenij i drugih sypuchih materialov Jakimov Ju.I., Ivanov V.P., Priporov E.V., Zajarskij V.P., Volkov G.I., Selivanovskij O.B. patent na izobretenie RUS 217721б 14.03.2000.
2. Centrobezhnyj apparat s podachej materiala vdol' lopatok Priporov E.V., Kartohin S.N. // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta. 2015. № 112. S. 1499-1511.
3. Centrobezhnyj rabochij organ dlja rasseva sypuchego materiala Jakimov Ju.I., Priporov E.V., Ivanov V.P., Zajarskij V.P., Volkov G.I., Selivanovskij O.B. patent na izobretenie RUS 2177217. zajavl.14.03.2000.
4. Patent na izobretenie RUS 2197807. Centrobezhnyj razbrasyvatel' sypuchih materialov. Jakimov Ju.I., Priporov E.V., Zajarskij V.P., Volkov G.I., Selivanovskij O.B. zajavl. 20.04.2001.
5. Pribor dlja issledovanija centrobezhnyh apparatov razbrasyvatelej sypuchih materialov Jakimov Ju.I., Priporov E.V., Karabanickij A.P., Tkachenko V.T., Jakushev A.A. patent na izobretenie RUS 2201059. zajavl. 20.04.2001.
6. Adler Ju.P., Markov E.V., Granovskij Ju.V. Planirovanie jeksperimenta pri poiske op-timal'nyh uslovij.- M.: Nauka, 1971 - S. 283.
7. Novikov F.S., Arsov Ja.B. Optimizacija processov tehnologii metallov metodami pla-nirovanija jeksperimentov. - M.: Mashinostroenie, 1980, 304s.
8. Soshniki zernovyh sejalok resursosberegajushhih tehnologij Priporov E.V. // V sbornike: Svjaz' teorii i praktiki nauchnyh issledovanij Sbornik statej Mezhdunarodnoj nauch-no-prakticheskoj konferencii. Otvetstvennyj redaktor: Sukiasjan Asatur Al'bertovich. 201б. S. б3-бб.
9. Opredelenie jenergosberegajushhego rezhima raboty tjagovogo agregata Priporov E.V. // Innovacii v sel'skom hozjajstve. 2015. № 5 (15). S. 92-95.
10. Izmel'chitel' grubyh kormov Maslov G.G., Artemov V.E., Priporov E.V., Nebavskij V.A. patent na izobretenie RUS 2222175. Zajavl. 11.0б.2002
11. Tehnologii i tehnicheskie sredstva prigotovlenija koncentrirovannyh kormov s ispol'zovaniem soevogo belka. Frolov V.Ju., Sysoev D.P., Priporov I.E., Gorb S.S.// International Scientific and Practical Conference World science. 2016. T. 1. № 3 (7). S. 53-58.
12. Priporov I.E., Lazebnyh D.V. Racional'naja tehnologija posleuborochnoj obrabotki semjan podsolnechnika // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Ku-banskogo gosudarstvennogo agrarnogo universiteta. 2015. № 112. S. 1475-1485.
13. Vlijanie tolshhiny, shiriny i individual'noj massy semjan podsolnechnika na skorost' ih vitanija Shaforostov V.D., Priporov I.E. // Maslichnye kul'tury. Nauchno-tehnicheskij bjul-leten' Vserossijskogo nauchno-issledovatel'skogo instituta maslichnyh kul'tur. 2010. № 1 (142-143). S. 76-80.
14. Obosnovanie primenenija opticheskogo fotojelektronnogo separatora v sostave uni-versal'nogo semjaochistitel'nogo kompleksa. Priporov I.E. // V sbornike: Konkurentnaja sposobnost' otechestvennyh gibridov, sortov i tehnologii vozdelyvanija maslichnyh kul'tur Sbornik materialov 8-j mezhdunarodnoj konferencii molodyh uchjonyh i specialistov. 2015. S. 138-141.
15. Klassifikacija opticheskih fotoseparatorov dlja sortirovanija semjan podsolnechnika. Priporov I.E., Shaforostov V.D. // Vestnik Kazanskogo gosudarstvennogo agrarnogo universiteta. 2015. T. 10. № 1. S. 68-70.
16. Modelirovanie processa separirovanija semjan podsolnechnika v vertikal'nom pnevmokanale vetro-reshetnyh zernoochistitel'nyh mashin Shaforostov V.D., Priporov I.E.// Maslichnye kul'tury. Nauchno-tehnicheskij bjulleten' Vserossijskogo nauchno-issledovatel'skogo instituta maslichnyh kul'tur. 2011. № 1 (146-147). S. 113-118.
