Научная статья на тему 'Аналитическое Исследование локальных температур при шлифовании'

Аналитическое Исследование локальных температур при шлифовании Текст научной статьи по специальности «Технологии материалов»

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

Аннотация научной статьи по технологиям материалов, автор научной работы — Унянин А. Н.

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

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

Похожие темы научных работ по технологиям материалов , автор научной работы — Унянин А. Н.

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

Текст научной работы на тему «Аналитическое Исследование локальных температур при шлифовании»

№ б 2006

621.923

АНАЛИТИЧЕСКОЕ ИССЛЕДОВАНИЕ ЛОКАЛЬНЫХ ТЕМПЕРАТУР

ПРИ ШЛИФОВАНИИ

Канд. техн. наук, доц. А. Н. УНЯНИН

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

Physical and mathematical models for calculation of temperatures at the grinding, allowing in a complex to consider the shape of contacting bodies i,e. abrasive grain, workpiece, a circle and a swaft, relative migration of these bodies, an arrangement of separate sources of a thermal emission, superposition of thermal impulses from separate grains, dependence of thermophysical properties of bodies and mechanical performances of work material from temperature« Adequacy of the selected models to real conditions is

proved. Results of a numerical modeling of temperatures are tested.

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

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

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

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

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

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

можно не принимать во внимание.

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

Схематизируя объекты, участвующие в теплообмене, представим а.з. в виде усеченной пирамиды с углом при вершине 2у площадкой затупления размером /2 [5], отдающей в массу круга часть теплоты. Форму заготовки примем в виде прямоугольного параллелепипеда: такую форму имеют многие детали, заготовки которых обрабатывают по схеме плоского маятникового шлифования периферией круга (корпусные детали, плиты, шпонки, линейки, измерительные плитки). Стружку, имеющую переменную толщину, представим в виде совокупности прямоугольных параллелепипедов различной толщины: такое допущение не внесет значительных погрешностей в расчеты, поскольку изменением толщины стружки на малой длине можно пренебречь. Часть объема круга, в котором расположено рассматриваемое зерно, представим в виде прямоугольного параллелепипеда. Пограничный слой между а.з., стружкой и заготовкой и его термическое сопротивление можно не учитывать, так как контакт между этими телами в условиях пластической деформации приближается к идеальному. Поэтому считаем, что температура в любой точке на поверхности соприкосновения контактирующих тел одна и та же.

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

Общая мощность тепловыделения при работе режущего и пластически деформирующего а.з.

где Ц?2Т — соответственно мощность тепловыделения источников в результате

перехода в теплоту работы деформирования и работы сил трения стружки о зерно (№{Т) и зерна о заготовку (УУ2Т), Вт.

Мб

2006

Все источники тепловыделения считаем двумерными, что не вносит существенных погрешностей в результаты расчета [4]. Мощности источников определим по следующим зависимостям:

2Т гггр у кр5

где Рг- — касательная составляющая силы диспергирования единичным а.з.; Ру и Р — силы трения соответственно между стружкой и а.з. и между а.з. и заготовкой; Укр — рабочая скорость крута; У1 = Укр/г| — скорость перемещения стружки относительно а.з. [4]; г| — коэффициент усадки стружки.

Принимая во внимание известные зависимости для расчета сил диспергирования и трения [5], получим для зерен, осуществляющих микрорезание,

W

кч т

1,5т, cos(90° - у - Р) ViM-oKp

W.

3

27 = Зц5 ;

и л m к « п кр 1 I 7

sin р Sin (pv

где xs — напряжение сдвига; ат — глубина внедрения зерна в заготовку; bk — ширина среза; (3 — угол сдвига; ф — угол внутреннего трения; |х — коэффициент внутреннего трения в плоскости сдвига; ц0 - коэффициент внешнего трения;/тр - площадь поверх-ности трения зерна о заготовку; 1{ — длина контакта стружки с поверхностью а,з,; у — угол при вершине зерна.

Для зерен, осуществляющих пластическую деформацию, имеем

W

27*

W ^hlLlL ° 2

3

f . cosv^ siny +-L

f V

Jкд кр>

где Нт — твердость материала заготовки при скорости и температуре деформации шлифования;— площадь поверхности трения пластически-деформирующего а.з. о заготовку; /кд — площадь поверхности зерна, осуществляющей пластическое деформирование материала заготовки.

Учитывали, что средние по плоскости сдвига касательные напряжения связаны с интенсивностью напряжений при скорости и температуре деформации процесса шлифования <7/ш зависимостью [5]: х5 =

На схеме теплообмена, показанной на рис. 1, действительная траектория движения режущего зерна представлена в виде прямой линии ГВ{. Очевидно, это допущение не внесет ощутимой погрешности в результаты расчета, поскольку длина дуги контакта кру-

№6

2006

Рис. I. Схема теплообмена в зоне контакта режущего а.з, с заготовкой: 1 — заготовка; 2 — а.з.; 3 ■— стружка;

4 — круг

га с заготовкой /к при реализации большинства схем шлифования ничтожно мала в сравнении с его размерами,

Считаем, что источник плотностью возникающий на плоскости сдвига, распределен равномерно [4], а соответствующая плотность тепловыделения равна

Я

о

0

о

апА

(1)

Плотность тепловыделения источника, действующего на поверхности контакта между стружкой и зерном, примем распределенной по несимметричному нормальному закону распределения [4] с максимальной плотностью потока, определяемой по

зависимости

С1

I т

Ък ^тсегГ

где к

01

коэффициент; егГ

функция, называемая интегралом вероятности.

№6

2006

Плотность тепловыделения в точке, расположенной на расстоянии х от точки с

максимальной плотностью тепловыделения,

<?17(*н) = 417-ехР[-Ал2]'

(3)

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

I

где I = 2 .— егГ

Ч

Ж

2 Т

2 Т

(4)

/

#2 Т

Ь2 л/тсегГ

(5)

Ьт и„) = 42Г ехр [~к0 хи ],

(б)

где к0 — коэффициент; 12 — размер площадки затупления на а.з.

Источник плотностью ци:, возникающий при работе а.з., пластически деформирующего материал заготовки (рис. 2), считаем распределенным равномерно

Чо

\¥0 СОБ у

7

Л КД

(7)

Плотность тепловыделения д2Т на участке контакта пластически деформирующего

зерна с заготовкой примем распределенной по экспоненциальному закону аналогично закону распределения нормальных напряжений. Для ее расчета используем зависимости (4)—(6), как и для режущих а.з.

Считаем, что заготовка 1 неподвижна, а режущее а.з. 2, шлифовальный круг стружка 3 и плоскость сдвига ОЕ перемещаются относительно заготовки со скоростью, равной сумме скоростей Укр и скорости подачи (рис. 1). Систему координат хоу свяжем с зерном и учтем, что она перемещается относительно заготовки в направлении оси ох со скоростью Укр. Такая схематизация относительных перемещений контактирующих объектов, когда один из них считают неподвижным, широко используется при теплофизичес-ком анализе шлифования и не влияет на точность расчета.

Уравнения теплопроводности для заготовки 7, режущего зерна 2, стружки 3 и шлифовального круга 4 имеют вид

Э7:

1

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

1

Эх схрх

Э

дх

/

V

Эл:

+

э

/

/

ду

X

1

V

Э7;

Эу

\

У

V.

кр

К,

дх

(8)

дТ] дх

1

С2Р

—(х ^

Эл' дх

V

Э

/

/

Эу

'дт,

ду

\

№6

2006

1

Рис. 2. Схема контакта заготовки с пластически деформирующим а.з.: 1 — заготовка; 2 — а.з.;

3 — шлифовальный круг

дТ, дх

1

дх дх

д

у

ду

к

\

дТ\

ду

\

ВТ, дТ,

V, э1п у

ду

дх '

(10)

К

дх

1

С4р4

Э7;4

4 Эх

дх

Э

ч

у

Эу

\

ду

(11)

где , Х2, — коэффициенты теплопроводности; сх, с2, с4 — теплоемкости; рр р2, р4 ■— плотности; Г2> Т4 — температуры заготовки, а.з., стружки и круга соответственно.

Уравнения теплопроводности для пластически деформирующего а.з. 2, заготовки 1 и круга 4 (рис. 2) аналогичны уравнениям (8), (9) и (11). Начальное условие

Цх, у, О) = Т()

означает равенство в начальный момент времени температур любой точки заготовки, а.з., стружки и круга температуре окружающей среды.

Условия взаимодействия поверхностей обрабатываемой заготовки, режущего а.з., стружки, круга и окружающей среды описываются следующими граничными условиями: в пределах зоны контакта заготовки 1 с а.з. 2

дТ

1

ду

<?2т дГ2 X, ' ду

Яг т

Ал

,7] — Т2 , ^27 Ягт

д2Т,

в пределах зоны контакта заготовки 1 со стружкой 3

дп

Чо

дт

X, дп

К

№6

2006

в пределах зоны контакта стружки 3 с а.з, 2

_ Яхт ___ / _

— Г — Г ' 2 ~~ 3' Я\г*

(Ж Л 2 о?! А1

Примем, что сопряжение а.з. 2 со связкой круга является идеальным, а это означает равенство температур сопряженных поверхностей зерна 2 и круга 4, т.е. = Г4.

Теплообмен на границах объектов, контактирующих со смазочно-охлаждающей жидкостью (СОЖ) или воздухом, зададим в форме закона Ньютона—Рихмана, описывающего процесс конвективного теплообмена, в виде

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

где г — номер объекта (г = 1, 2, 3, 4 для заготовки, а.з., стружки и круга соответственно); ; — номер поверхности; осу< — коэффициент теплоотдачи оту-ой поверхности.

Границы Р(2 и области Р0№§ шлифовального круга считаем теплоизолированными.

Условия взаимодействия при работе пластически деформирующего а,з, (рис. 2) имеют особенности в пределах зоны контакта заготовки 1 с «передней» поверхностью а.з. 2

дТх qD ЭГ2 q"D , „

дп Xj дп X

2

Время, в течение которого происходит теплообмен, разобьем на конечные малые промежутки Дт, Через к обозначим число положений а,з. относительно дуги контакта с

заготовкой

*=—+ 1. КрАт

Глубина внедрения зерна в заготовку при его /-ом положении

а . = (i -f a

mi кр 3i>

где а — критическая глубина микрорезания; аз/ — подача на зерно при г-м положении

кр

последнего:

a3j = —— sin

max кр

O'-i)- ,

V к J

где /х — среднее расстояние между диспергирующими зернами; пк — частота вращения

шлифовального круга; Vs — скорость продольной подачи; DK — диаметр наружной поверхности круга; атах — угол, определяемый как

Ас-2гф «n«x = arCC0S-

где tih — глубина шлифования.

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

№6

2006

I V

к кр

П =-—

VV<>

Таким образом, чтобы учесть взаимное влияние тепловых импульсов, необходимо моделировать движение по траектории РВХ зерен, число которых равно п. Начальную температуру заготовки при воздействии на нее ;-го зерна (2 < ] < п) принимаем равной температуре, до которой она успевает остыть после прохождения (/' - 1)-го зерна. При моделировании температурного поля от первого зерна считаем, что начальная температура равна температуре окружающей среды.

Для расчета температур необходимо определить коэффициенты теплоотдачи от заготовки, а.з., стружки и выделенного объема шлифовального круга к СОЖ и воздуху. Поверхности ЫК и УЬ а.з., МО и МЫ стружки, СБ и РТ7 заготовки, а также РК и ¿5 круга обмениваются теплотой с СОЖ, проникающей в зону контакта круга с заготовкой; поверхности В^С и заготовки омываются жидкостью, подаваемой поливом (рис. 1). Остальные поверхности заготовки, а также поверхность шлифовального круга обмениваются теплотой с окружающим воздухом. При определении коэффициентов теплоотдачи от поверхностей объектов, находящихся в зоне контакта круга с заготовкой, учитывали, что в эту зону проникает газожидкостная среда (смесь СОЖ с воздухом), а в зависимости от температуры поверхности устанавливается либо режим конвективного теплообмена, либо теплообмена при кипении [6].

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

механических свойств материала заготовки в зависимости от температуры создают серьезные трудности для аналитического определения локальных температур. Однако эту задачу можно решить численным методом (использовали метод конечных разностей) [7].

Для каждого типа расчетных элементов а.з., заготовки, стружки и круга методом баланса [7] получали дискретные аналоги соответствующих дифференциальных уравнений теплопроводности. Значение устойчивого шага интегрирования по времени определяли до начала расчета теплового состояния контактирующих объектов при двух значениях температуры: минимальной Г0 и максимально возможной. Из всех полученных значений выбирали наименьшее, которое и принимали в качестве шага интегрирования по времени Ат.

Тепловое состояние системы рассчитывали на основе дискретных аналогов уравнений теплопроводности в соответствии с разработанной методикой. Методика реализована в оригинальных программах для персонального компьютера на языке Microsoft Visual Basic.

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

В качестве примера приведем некоторые результаты численного моделирования при следующих исходных данных: материал заготовки — сталь 3ОХ; материал а.з. — элект-

Nq 6

2006

рокорунд белый; Укр = 35 м/с; Уу = 10 м/мин; глубина шлифования Гф = 0,01 мм; СОЖ эмульсия Укринол-1 (3 %); варьировали коэффициентом, определяющим долю СОЖ в объеме воздушно-жидкостной смеси, проникающей в зону шлифования [б], в пределах 0,25 ,..., 1,0, а также коэффициентом трения зерна о заготовку от 0,3 до 0,6.

Установлено, что вблизи от вершины а.з. (точки Е на рис. 1) на площадке его контакта с заготовкой локальные температуры достигают температуры плавления материала зерна» На расстоянии 2 мкм вглубь материала а.з. максимальная температура не превышает 2000 К (температура плавления электрокорунда — 2120. .,2170 К), т. е. до температуры плавления разогревается очень тонкий поверхностный слой зерна. На участке, размер которого составляет около 70 % от размера площадки затупления а.з., температура превышает температуру плавления материала заготовки и, следовательно, создаются наиболее благоприятные условия для адгезионного взаимодействия а.з. с заготовкой [1]. Максимальная температура на поверхности контакта а.з. со стружкой ниже и не превышает 1800 К, поэтому интенсивность адгезионного взаимодействия на этом участке соответственно также ниже [2]. Локальные температуры а.з., пластически деформирующих материал заготовки, существенно ниже температуры режущих а.з.

Охлаждающее действие СОЖ практически не влияет на локальные температуры. Увеличение коэффициента внешнего трения от 0,3 до 0,6 привело к увеличени. как локальных температур на а.з., так и температур в поверхностных слоях заготовки почти в 2 раза, что свидетельствует о недопустимости значительного увеличения этого коэффициента. Поэтому важно принимать меры, способствующие предотвращению или замедлению процесса налипания частиц материала заготовки на а.з, или следует удалять отходы с рабочей поверхности круга, используя, например, механическую очистку [8].

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

СПИСОК ЛИТЕРАТУРЫ

L X у д о б и н Л. В., У н я н и н А. Н. Исследование кинетики процесса налипания частиц материала заготовки на абразивные зерна// Фундаментальные и прикладные проблемы технологии машиностроения; Материалы междун. научно-техн. конф. — Орел: Орловский госуд. технич. университет, 2002. -— С. 53— 56.

2. Худобин J1. В,, У н я н и н А. Н. Исследование процесса микрорезания заготовок из пластичных материалов единичным абразивным зерном // Процессы абразивной обработки, абразивные инструменты и материалы: Сб. статей международной научно-технической конференции. — Волжский ин-т строит, и технологии, Волжский, 2004. — С. 108—111.

3. БокучаваГ. В. Трибология процесса шлифования — Тбилиси: Сабчота Сакартвело, 1984. — 238 с.

4. Резников А.Н. Тепловые процессы в технологических системах. — М: Машиностроение, 1990. — 288 с.

5. К о р ч а к С. Н. Производительность процесса шлифования стальных деталей. — М.: Машиностроение, 1974. —280 с,

6. Унянин А. Н. К вопросу о теплообмене контактирующих при шлифовании объектов с окружающей средой // Физические и компьютерные технологии в народном хозяйстве: Сб. трудов международной научно-технической конференции. —• Харьков; ХНПК «ФЭД», 2002. — С. 79—82.

7. Ш и Д. Численные методы в задачах теплообмена: пер. с англ. — М.: Мир, 1988. — 544 с.

Мб

2006

8. У н я н и н А. Н. Стабилизация режущей способности шлифовального круга при обработке заготовок из пластичных материалов//Динамика технологических систем: Сб. трудов 7-ой междун. научно-техн. конференции. — Саратов: Саратовский госуд. технич. университет, 2004. — С. 353—356.

МЕТОДИКА ОЦЕНКИ СРОКА СЛУЖБЫ ГИДРОЦИЛИНДРОВ УРАВНОВЕШИВАНИЯ РАБОЧИХ ВАЛКОВ ЛИСТОВЫХ СТАНОВ*

Д-р тгхн наук, проф. В, П. АНЦУПОВ, ст. препод. А.А. БАЖЕНОВ, канд. техн. наук, доц. А.В. АНЦУПОВ, асп. Ал. В. АНЦУПОВ

Одной из основных причин остановки стана листовой прокатки является износ элементов механизма уравновешивания рабочих валков, а именно, износ бронзовых втулок и резинотканевых уплотнений. Анализ литературных данных показывает отсутствие в настоящее время какой-либо инженерной методики расчета ресурса работы узлов гидр о- и пневмосистем типа « плунжер-втул ка». Поэтому в данной работе сделана попытка оценки межремонтного срока службы гидрог^илиндров на основе математической модели процесса изнашивания пар трения скольжения гидроцилиндра («плунжер-бронзовая втулка», «плунжер-уплотнение»).

One of the principal causes of a contour rolling mill stopping is component deterioration of the work roll balance mechanism, namely, deterioration of bronze sleeves and rubber-fabric seals. The analysis of published research data displays absence of any engineering design procedure as for resource of operation of hydraulics-and pneumatic systems of a «plunger-plug» type. Therefore in this article we attempted to estimate a turnaround time of hydraidic actuators on the basis of mathematical model of deterioration

process of sliding couplers («plunger-bronze sleeve», «plunger-condensation»).

Модель построена с использованием энергетического подхода В.Д. Кузнецова [1], согласно которому изношенный объем детали АV пропорционален работе сил трения, совершенной в контакте пары трения при ее эксплуатации А .

где — показатель энергетической интенсивности изнашивания (энергетическая интенсивность изнашивания).

В основу расчета /н, и А^ (7) положена базовая схема И.В. Крагельского [2], согласно которой рассматривается стационарный (установившийся) режим граничного трения, определяемый стабильностью во времени фрикционных характеристик сопряжения, температуры и физико-механических свойств поверхностных слоев трущихся тел.

Рассмотрим поочередно пары скольжения гидроцилиндра.

Объемный износ бронзовой втулки из геометрических соображений (рис. 2)

669.1.002.5-192

А^(0 = /,Л(0 >

(1)

(2)

Sadc = $аосв - Sacs = 'Л ~ {rm ~ АЯВТ - (гвт - У)) гпл sin (ср,),

*При финансовой поддержке Правительства Челябинской области.

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