Научная статья на тему 'Влияние влажности сыпучих сред на их распределение в рабочем объеме смесителя с подвижной лентой'

Влияние влажности сыпучих сред на их распределение в рабочем объеме смесителя с подвижной лентой Текст научной статьи по специальности «Физика»

CC BY
206
34
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
БАРАБАННО-ЛЕНТОЧНЫЙ СМЕСИТЕЛЬ / БИЛА / ВЛАЖНОСТЬ / СЫПУЧИЕ КОМПОНЕНТЫ / УГОЛ РАСПЫЛИВАНИЯ / ФАЗОВЫЙ ОБЪЕМ / ФУНКЦИЯ РАСПРЕДЕЛЕНИЯ / HUMIDITY / LOOSE COMPONENTS / SPRAYING ANGLE / PHASE VOLUME / DISTRIBUTION FUNCTION / BEATERS / DRUM-RIBBON BLENDER / MIXING

Аннотация научной статьи по физике, автор научной работы — Капранова Анна Борисовна, Бакин Михаил Николаевич, Верлока Иван Игоревич, Зайцев Анатолий Иванович

С учетом влажности смешиваемых компонентов предложен стохастический метод определения дифференциальных функций распределения для частиц потоков сыпучих материалов при срыве с бил, закрепленных по винтовой линии на поверхности барабана, по углу распыливания. Согласно стохастическому подходу для описания указанной технологической операции применяется теория ударных взаимодействий дисперсных систем. Выполнен анализ указанных функций с примерами расчета их параметров при смешивании сыпучих сред двух сортов манной крупы и природного песка. Теоретически установлено, что с учетом вращательной части движения частиц смешиваемых компонентов при оценке их стохастических энергий в разреженных потоках максимальный угол рассеивания частиц существенным образом зависит от шага винтовой линии, вдоль которой закреплены била на поверхности барабана. Полученные результаты могут быть использованы для оценки качества получаемой смеси, в том числе при порционном внесении в смесь одного из компонентов, например c объемно-весовым соотношением 1:10 и более.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Капранова Анна Борисовна, Бакин Михаил Николаевич, Верлока Иван Игоревич, Зайцев Анатолий Иванович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Influence of Particles Humidity on their Distribution in the Working Volume of a Blender with a Moving Ribbon

Given the humidity of the blend components we proposed a stochastic method for determining the differential distribution functions of particles mixed flows of bulk materials at the disruption from beaters enshrined in a helical line on the drum surface, the angle of atomization. According to the stochastic approach for the description of the specified process step, we applied the theory of shock interaction of disperse systems. The analysis of the mentioned functions was performed and illustrated by the examples of calculation of the parameters during mixing of the bulk solids of two sorts farina and natural sand. Taking into account the rotational motion of the particles of the components to be mixed, and assessing their energies in stochastic flow dilution, we found that the maximal angle of particles dispersion depended on twist of the line along which the beaters were fastened to the drum’s surface. The results can be used to assess the quality of the granular mixture, including when making a batch in a mixture of one of the components, for example, with volume-weight correlation 1:10 and more.

Текст научной работы на тему «Влияние влажности сыпучих сред на их распределение в рабочем объеме смесителя с подвижной лентой»

УДК 66.022.54

Б01: 10.17277/уе81шк.2015.01.рр.097-104

ВЛИЯНИЕ ВЛАЖНОСТИ СЫПУЧИХ СРЕД НА ИХ РАСПРЕДЕЛЕНИЕ В РАБОЧЕМ ОБЪЕМЕ СМЕСИТЕЛЯ С ПОДВИЖНОЙ ЛЕНТОЙ

А. Б. Капранова, М. Н. Бакин, И. И. Верлока, А. И. Зайцев

Кафедра «Теоретическая механика», ФГБОУ ВПО «Ярославский государственный технический университет», г. Ярославль; [email protected]

Ключевые слова: барабанно-ленточный смеситель; била; влажность; сыпучие компоненты; угол распыливания; фазовый объем; функция распределения.

Аннотация: С учетом влажности смешиваемых компонентов предложен стохастический метод определения дифференциальных функций распределения для частиц потоков сыпучих материалов при срыве с бил, закрепленных по винтовой линии на поверхности барабана, по углу распыливания. Согласно стохастическому подходу для описания указанной технологической операции применяется теория ударных взаимодействий дисперсных систем. Выполнен анализ указанных функций с примерами расчета их параметров при смешивании сыпучих сред двух сортов - манной крупы и природного песка. Теоретически установлено, что с учетом вращательной части движения частиц смешиваемых компонентов при оценке их стохастических энергий в разреженных потоках максимальный угол рассеивания частиц существенным образом зависит от шага винтовой линии, вдоль которой закреплены била на поверхности барабана. Полученные результаты могут быть использованы для оценки качества получаемой смеси, в том числе при порционном внесении в смесь одного из компонентов, например с объемно-весовым соотношением 1:10 и более.

Введение

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

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

Учет влажности смешиваемых компонентов в стохастическом описании поведения их потоков

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

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

т1У= 6~1икм,1рТ1ё11, (1)

описывается в фазовом пространстве с элементом вида

ё = П ^ ¿уХ1]ауу1] = П г]=1 , (2)

где индекс 1 соответствует условному номеру смешиваемого материала, 1 = 1,..., Пк; индекс ] - номеру деформированного гибкого элемента в первой четверти отсчета угловой координаты полярной системы, ] = 1,...,пь; к№1 - параметр, учитывающий отношение масс влажной и невлажной частицы, (кг-м)/рад; Рт - истинная плотность вещества данного компонента, кг/м3; ё^ - диаметр сферической частицы, м; ю - угловая скорость вращения барабана, рад/с.

Здесь с учетом приближения о сравнимости кинематических характеристик (в частности, уХу] или ущ проекций на поперечную плоскость вращающегося

барабана скоростей движения частиц -го компонента смеси при их отрыве от конца ]-го била в декартовой или полярной системах координат) выражение (2) может быть записано в форме

2

ёО.] = ёё^ёТ ] =-ю ёё^1Г]ёГ]ёд]. (3)

Аналогично развитому в работах [4 - 6] подходу стохастическая энергия Е]

смешиваемых зернистых компонентов рассчитывается согласно энергиям ее поступательного и вращательного движений [5], случайного характера момента импульсов частицы при срыве с била [6], а также энергии этой частицы при упругом взаимодействии с билом, угловой коэффициент жесткости которого равен ки [7]

Е] = т^щ!2 + 5т 1 (к/угШ] ))(2ё1) + ки&]2/2, (4)

Заметим, что поиск Е] зависит от вида моделируемого уравнения винтовой линии г8 (9]) в указанной поперечной плоскости, как геометрическое местоположение проекций концов бил на эту плоскость, когда кратчайшее расстояние от вектора скорости частицы у^ в полярной системе координат до оси вращения

барабана равно

hJ (e J ) {b (e J )]2 + R0 } . (5)

где Щ = 2к,(I- -к0)Ь-1; к, - шаг винтовой линии на цилиндрической поверхности барабана, м; Ь- - длина барабана, м; I- - длина била, м; к0) - высота валково-ленточного зазора, м.

Следовательно, распределение частиц разреженного потока, описываемое [4]

с помощью выражения йЫ1- = Л1-ехр (-Еу!^]- )О-- с параметром нормировки Л1-, по углу распыливания частиц а,- в поперечной плоскости вращению барабана определяется выражением

Р1 (а1, а2, -,ап- ) = П-=1 ^ (а] ), (6)

когда

(а- ) = ] йЫ1]/йа], (7) или согласно (3), (4) соответствует функциям

fJ (aJ ) = w2iJ { erf I w2j [l + g3 (aJ + Ф J )]2 ) - erf ( )

:exp

-ku (Eoj ) 1 (aJ + фJ)

(8)

где коэффициенты gз = ^ (г- + И)) 1.

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

Анализ дифференциальных функций распределения для частиц смешиваемых потоков с примером расчета их параметров

В частности, приведем примеры расчета Е1 (а!,а2,аз) для двух сыпучих

компонентов (щ = 2) при следующих значениях (рис. 1, 2) параметров: конструктивных - гь = 3,0 ■ 10-2 м; 1ь = 4,5 ■ 10-2 м; к, = (1,6 - 3,0) ■ 10-2 м; Ьь = 1,85 ■ 10-1 м;

—4 / \

^ = 5,0-10 (кгм)/рад при п- = 3 - числе билов (] = 1,..., п-) с поворотом бара-

—2

бана на угол п / 2 рад; режимных - в диапазонах изменения к = (2,0 . 3,0) 10 м; ю = (47 ... 53) рад/с. Физико-механические характеристики каждого из смешиваемых материалов соответственно равны - для манной крупы ГОСТ 7022-97 (1 = 1): = 4,0-10—4 м; рп = 1,440-103 кг/м3; = (1,0 ... 1,4) и для природного песка

ГОСТ 8736-93 (1 = 2): ё52 = 1,5-10—4 м; рТ2 = 1,525-103 кг/м3; ^ = (1,0 ... 1,4). Анализ графиков для функций Ф^а -, ю) и Ф2(а -, ю) в случае -2

к, = 3,0 ■ 10 м; ю = 47,12 рад/с при различных влажностях (рис. 1) показывает, что для песка с ростом массы его частиц за счет увлажнения их поверхности в 1,4 раза (кривые 1 и 3, см. рис. 1, б) наблюдается уменьшение максимального значения угла распыливания в 1,2 раза и некоторое сужение дисперсионного интервала.

F1

0,2

0,1

2 3 а)

ау, рад

F2 3

/jT 2

0,5 и А,;

0,4

0,3

0,2

а], рад

0 0,2 0,4 0,6 0,8

б)

Рис. 1. Зависимости ^(оо/) и ^(а/) для манной крупы ГОСТ 7022-97 (а) и природного песка ГОСТ 8736-93 (б) при режимных параметрах = 3,0-10 м; ю = 47,12 рад/с и различных степенях их влажности:

1 - кц1 = 1,0; 2 - кц1 = 1,2; 3 - кЦг- = 1,4

F1 0,18

0,14

0,10

0,06

0,02

F2 0,2

0,1

а], рад

а], рад

0 1

2 3 а)

4

0

1

4

2 3

б)

Рис. 2. Зависимости ^(оо,) и ^2(а,) для манной крупы ГОСТ 7022-97 (а)

«—2

и природного песка ГОСТ 8736-93 (б) при режимных параметрах ^ = 1,6-10 м; ю = 52,36 рад/с и различных степенях их влажности:

1 - кц1 = 1,0; 2 - кщ = 1,142; кц2 = 1,074; 3 - кц1 = 1,2

Аналогичная картина прослеживается при = 1,6 -10 м; ю = 52,36 рад/с; (рис. 2), когда увеличение массы частиц песка в 1,2 раза (кривые 1 и 3, см. рис. 2, б) приводит к снижению максимального значения угла распыливания для второго

компонента (природного песка ГОСТ 8736-93) а^ах в 1,15 раза.

Для манной крупы ГОСТ 7022-97 характер поведения частиц несколько другой - с ростом увлажнения составляющих этот сыпучий компонент приобретает отклонение на больший угол по сравнению с менее влажным вне зависимости от значений шага навивки гибких элементов И8 (кривые 1 и 3, см. рис. 1, а и 2, а). Это объясняется двумя причинами: 1) слипаемость изначально более крупных частиц манной крупы приводит к достижению большего углового разлета их со-

0

1

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

Таким образом, построенная модель смешивания в его рабочем объеме бара-банно-ленточного аппарата с расположением бил по винтовой линии в поперечной плоскости вращения барабана имеет следующие особенности: применима при смешивании нескольких сыпучих компонентов без традиционного деления на ключевой и транспортирующий материалы; определяет дифференциальные функции распределения смешиваемых компонентов и смеси по углам распространения потоков с учетом формы твердых частиц и их возможной влажности. Заметим, что описание поведения одиночных частиц при входе в слой сыпучей смеси рассматривается в работах [8, 9].

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

Заключение

Предложенное стохастическое описание распределения смешиваемых компонентов в рабочем объеме барабанно-ленточного аппарата выполнено с учетом влажности зернистых материалов. При этом теоретически установлено, что с учетом вращательной части движения частиц смешиваемых компонентов при оценке их стохастических энергий в разреженных потоках максимальный угол рассеивания частиц существенным образом зависит от шага винтовой линии, вдоль которой закреплены била на поверхности барабана. В частности при снижении указанного шага в 1,88 раза наблюдается эффект сближения дифференциальных функций распределения для усредненных потоков компонентов с размерами частиц, отличающимися в 2,7 раза, при изменении отношения значений их максимальных углов распыливания а]^/а]^ от 5 до 1,1 ед., что позволяет теоретически выявить условия наиболее качественного смешивания указанных сыпучих компонентов.

Полученные результаты в соответствии с выражениями (6) - (8) для полной дифференциальной функции распределения сыпучих компонентов по углам рас-пыливания для каждого ]-го деформированного била (] = 1,..., п-), установленного по винтовой линии на поверхности вращающихся барабанов аппарата, позволяют оценить один из основных показателей качества смеси - коэффициент ее неоднородности [10].

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

1. Бакин, М. Н. Современные аппараты с подвижной лентой для смешивания сыпучих материалов [Электронный ресурс] / М. Н. Бакин, А. Б. Капранова, И. И. Верлока // Фундаментальные исследования. - 2014. - № 5(ч. 4). - С. 687 - 691. -Режим доступа : www.rae.ru/fs/?section=content&op=show_article&article_id=10003 246 (дата обращения: 25.02.2015).

2. Явник, П. М. Моделирование процесса непрерывного весового дозирования / П. М. Явник, С. В. Першина, В. Ф. Першин // Вестн. Тамб. гос. техн. университета. - 2012. - Т. 18, № 4. - С. 912 - 916.

3. Бакин, М. Н. Современные методы математического описания процесса смешивания сыпучих материалов [Электронный ресурс] / М. Н. Бакин, А. Б. Капранова, И. И. Верлока // Фундаментальные исследования. - 2014. - № 5(ч. 5). -С. 923 - 927. - Режим доступа : www.rae.ru/fs/?section=content&op=show_article& агйс1е_1ё=10003287 (дата обращения: 25.02.2015).

4. Бакин, М. Н. Исследование распределения сыпучих компонентов в рабочем объеме барабанно-ленточного смесителя [Электронный ресурс] / М. Н. Бакин, А. Б. Капранова, И. И. Верлока // Фундаментальные исследования. - 2014. -№ 5(ч. 5). - С. 928 - 933. - Режим доступа : www.rae.ru/fs/?section=content&op= show_artic1e&artic1e_id= 10003288 (дата обращения: 25.02.2015).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5. Зайцев, А. И. Ударные процессы в дисперсно-пленочных системах / А. И. Зайцев, Д. О. Бытев. - М. : Химия, 1994. - 176 с.

6. Исследование ударного смешивания твердых дисперсных сред при вторичных столкновениях частиц / А. Б. Капранова [и др.] // Изв. ВУЗов. Химия и хим. технология. - 2013. - Т. 56, № 6. - С. 83 - 85.

7. Лебедев, А. Е. Математическая модель механики движения сыпучих материалов в разреженных потоках аппаратов с эластичными рабочими элементами / А. Е. Лебедев, А. И. Зайцев, А. Б. Капранова, И.О. Кузьмин // Изв. ВУЗов. Химия и химическая технология. - 2009. - Т. 52, № 5. - С. 111 - 113.

8. Капранова, А. Б. Исследование этапа ударного входа одиночной частицы в неподвижный слой сыпучей среды / А. Б. Капранова, А. А. Петров, А. И. Зайцев // Изв. ВУЗов. Химия и хим. технология. - 2013. - Т. 56, № 9. - С. 109 - 111.

9. Моделирование процесса образования кратера в слое сыпучей среды после удара одиночной частицы / А. Б. Капранова [и др.] // Изв. ВУЗов. Химия и хим. технология. - 2013. - Т. 56, № 10. - С. 119 - 121.

10. Оценка коэффициента неоднородности зернистой смеси в объеме бара-банно-ленточного устройства / А. Б. Капранова [и др.] // Изв. ВУЗов. Химия и хим. технология. - 2014. - Т. 57, № 9. - С. 104 - 106.

Influence of Particles Humidity on their Distribution in the Working Volume of a Blender with a Moving Ribbon

A. B. Kapranova, M. N. Bakin, I I. Verloka, A. I. Zaitsev

Department "TheoreticalMechanics", Yaroslavl State Technical University; [email protected]

Keywords: beaters; distribution function; drum-ribbon blender; humidity; loose components; mixing; phase volume; spraying angle.

Abstract: Given the humidity of the blend components we proposed a stochastic method for determining the differential distribution functions of particles mixed flows of bulk materials at the disruption from beaters enshrined in a helical line on the drum surface, the angle of atomization. According to the stochastic approach for the description of the specified process step, we applied the theory of shock interaction of disperse systems. The analysis of the mentioned functions was performed and illustrated by the examples of calculation of the parameters during mixing of the bulk solids of two

sorts - farina and natural sand. Taking into account the rotational motion of the particles of the components to be mixed, and assessing their energies in stochastic flow dilution, we found that the maximal angle of particles dispersion depended on twist of the line along which the beaters were fastened to the drum's surface. The results can be used to assess the quality of the granular mixture, including when making a batch in a mixture of one of the components, for example, with volume-weight correlation 1:10 and more.

References

1. www.rae.ru/fs/?section=content&op=show_article&article_id=10003246 (assecced 25 February 2015).

2. Yavnik P.M., Pershina S.V., Pershin V.F. Transactions of the Tambov State Technical University, 2012, vol. 18, no. 4, pp. 912-916.

3. www.rae.ru/fs/?section=content&op=show_article&article_id=10003287 (assecced 25 February 2015).

4. www.rae.ru/fs/?section=content&op=show_article&article_id=10003288 (assecced 25 February 2015).

5. Zaitsev A.I., Bytev D.O. Udarnye protsessy v dispersno-plenochnykh sistemakh (Shock processes in dispersed-film systems), Moscow: Khimiya, 1994, 176 p.

6. Kapranova A.B., Bakin M.N., Lebedev A.E., Zaitsev A.I. Izvestiya VUZov. Khimiya i khimicheskaya tekhnologiya, 2013, vol. 56, no. 6, pp. 83-85.

7. Lebedev A.E., Zajcev A.I., Kapranova A.B., , Kuz'min I.O. Izv. VUZov. Himija i himicheskaja tehnologija, 2009, vol. 52, no. 5, pp. 111-113.

8. Kapranova A.B., Petrov A.A., Zaitsev A.I. Izvestiya VUZov. Khimiya i khimicheskaya tekhnologiya, 2013, vol. 56, no. 9, pp. 109-111.

9. Kapranova A.B., Petrov A.A., Bakin M.N., Zaitsev A.I. Izvestiya VUZov. Khimiya i khimicheskaya tekhnologiya, 2013, vol. 56, no. 10, pp. 119-121.

10. Kapranova A.B., Bakin M.N., Lebedev A.E., Zaitsev A.I. Izvestiya VUZov. Khimiya i khimicheskaya tekhnologiya, 2014, vol. 57, no. 9, pp. 104-106.

Einfluss der Feuchtigkeit der streubaren Umgebungen auf ihre Verteilung im Arbeitsumfang des Mischers mit dem beweglichen Band

Zusammenfassung: Unter Berücksichtigung der Feuchtigkeit der gemischten Komponenten ist die stochastische Methode der Bestimmung der Differentialfunktionen der Verteilung für die Teilchen der gemischten Ströme der streubaren Materialien bei der Vereitelung von den Schläger, die t nach der Schraubenlinie auf der Oberfläche der Trommel gefestigt sind, nach dem Winkel des Sprühens. Laut dem stochastischen Herangehen für die Beschreibung der angegebenen technologischen Operation die Theorie der Stosswechselwirkungen der dispersen Systeme verwendet. Es ist die Analyse der angegebenen Funktionen mit den Beispielen der Berechnung ihrer Parameter bei der Vermischung der streubaren Umgebungen zwei Sorten - des Grießes und des natürlichen Sandes erfüllt. Es ist theoretisch bestimmt, dass unter Berücksichtigung des drehenden Teiles der Bewegung der Teilchen der gemischten Komponenten bei ihrer Einschätzung der stochastischen Energien in den entladenen Ströme der maximale Winkel der Streuung der Teilchen vom Schritt der Schraubenlinie in der wesentlichen Weise abhängt, entlang der gefestigt sind schlug auf der Oberfläche der Trommel. Die bekommenen Ergebnisse können für die Einschätzung der Qualität der bekommenen Mischung verwendet sein, Einschließlich, bei der Portionseintragung in die Mischung einen der Komponenten, zum Beispiel, vom räumlich-Waageverhältnis 1:10 und mehr.

Influence de l'humidité des milieux pulvérilents sur leur répartition

dans un volume de travail de mélangeur avec une bande mobile

Résumé: Est proposée une méthode stochastique de la définition des fonctions différentielles de la répartition pour les particules des courants mélangés des matériaux pulvérilents compte tenu de l'humidité des composants. Suivant une approche stochastique pour la définition de l'opération mentionnée est employée la théorie des interactions de choc des systèmes de dispersion. Est effectuée l'analyse des fonctions indiqueées avec les exemples du calcul de leurs paramètres. Les résultats reçus peuvent être utilisés pour l'estimation de la qualité du mélange y compris lors de l'addition par portions dans le mélange d'un des composants avec la relation volume-poids 1:10 et plus.

Авторы: Капранова Анна Борисовна - доктор физико-математических наук, доцент, профессор кафедры «Теоретическая механика»; Бакин Михаил Николаевич - соискатель кафедры «Теоретическая механика»; Верлока Иван Игоревич - аспирант кафедры «Теоретическая механика»; Зайцев Анатолий Иванович -доктор технических наук, профессор, заведующий кафедрой «Теоретическая механика», ФГБОУ ВПО «Ярославский государственный технический университет», г. Ярославль.

Рецензент: Петерсон Станислав Антонович - доктор технических наук, профессор, заведующий кафедрой «Сопротивление материалов», ФГБОУ ВПО «Ярославский государственный технический университет», г. Ярославль.

i Надоели баннеры? Вы всегда можете отключить рекламу.