Научная статья на тему 'Компьютерное моделирование процессов кристаллизации как резерв повышения качества поршней ДВС'

Компьютерное моделирование процессов кристаллизации как резерв повышения качества поршней ДВС Текст научной статьи по специальности «Прочие технологии»

CC BY
108
28
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЛИТАЯ ДЕТАЛЬ / КОКИЛЬ / КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ / ТЕРМОИЗОЛИРУЮЩЕЕ ПОКРЫТИЕ / MOLDED PIECE / METAL MOLD / COMPUTER MODELING / THERMAL INSULATING COATING

Аннотация научной статьи по прочим технологиям, автор научной работы — Пономаренко О. И., Тренёв Н. С.

В статье предложено использование процедуры распознавания образов для описания локализации дефектов в литых деталях «поршень» для ДВС. Приведены результаты компьютерного моделирования процесса заполнения кокиля и кристаллизации сплава для усовершенствованной технологии литья, предполагающей применение термоизолирующих покрытий кокиля. Результаты компьютерного моделирования могут быть использованы для формализации процесса описания локализации дефектов в поршнях.

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

Computer modeling of crystallization processes as a reserve of improving the quality of pistons of ICE

The use of pattern recognition procedure for the description of the localization of defects in the molded pieces «piston» for ICE (internal combustion engine) is proposed in the paper. The relevance of this research is explained by the need of rapid development of new high quality molded pieces for the automotive industry, or optimization of the existing equipment design to improve the product quality. The results of computer modeling of the process of metal mold casting and alloy crystallization for the advanced casting technology, involving the use of thermal insulating coatings of metal mold are given. The results of computer modeling can be used for formalization of the process of describing the localization of defects in pistons and, based on this, for development of measures on resource saving and improving the quality of molded pieces. It was proposed to use the results of statistical classification on the basis of such formalized description.

Текст научной работы на тему «Компьютерное моделирование процессов кристаллизации как резерв повышения качества поршней ДВС»

16. Chien-Ching, Ma. Analytical exact solutions of heat conduction problems for anisotropic multi-layered media [Текст] / Chien-Ching Ma, Shin-Wen Chang // Int. Journal of Heat and Mass Transfer. - 2004. - Vol. 47. - P. 1643-1655.

17. Ozisik, M. N. Heat Conduction Wiley [Текст] / M. N. Ozisik // New York, 1993.

18. Bejan, A. Heat Transfer [Текст] / A. Bejan. — John Wiley & Sons, New York, 1993.

19. Dryden, I. G. C. The Efficient Use of Energy [Текст] / I. G. C. Dryden. - Ed. 2. - Butterworth Scientific, Oxford, 1982.

ОБОСНОВАНИЕ ПАРАМЕТРОВ И ЧИСЛЕННЫЕ ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ ЭЛЕКТРООБОГРЕВАЕМОГО ПОЛА

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

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

Романченко Микола Анастасійович, кандидат технічних наук, професор, кафедра інтегрованої електротехнології та процесів,

Харківський національний технічний університет сільського господарства ім. Петра Василенка, Україна.

Слесаренко Анатолій Павлович, доктор фізико-математич-них наук, професор, ведучий науковий співробітник, Лауреат Державної премії України, Інститут проблем машинобудування ім. А. Н. Підгорного НАН України, Україна.

Сорока Олександр Степанович, кандидат фізико-матема-тичних наук, доцент, кафедра мікроелектроніки, електронних приладів і пристроїв, Харківський Національний університет радіоелектроніки, Україна.

Романченко Николай Анастасиевич, кандидат технических наук, профессор, кафедра интегрированных электротехнологий и процессов, Харьковский Национальный технический университет сельского хозяйства им. Петра Василенка, Украина. Слесаренко Анатолий Павлович, доктор физико-математических наук, профессор, Лауреат Государственной премии Украины, Ведущий научный сотрудник, Институт проблем машиностроения им. А. Н. Подгорного НАН Украины, Украина. Сорока Александр Степанович, кандидат физико-математических наук, доцент, кафедра микроэлектроники, электронных приборов и устройств, Харьковский Национальный университет радиоэлектроники, Украина.

Romanchenko Nickolay, Kharkiv Petro Vasylenko National Technical University of Agriculture, Ukraine.

Slesarenko Anatoliy, A. N. Podgorny Institute for Mechanical Engineering Problems of NAS of Ukraine, Ukraine.

Soroka Alexander, Kharkiv National University of Radioelectronics, Ukraine

УДК 621.74

Пономаренко 0. И., КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ Тренёв Н'с ПРОЦЕССОВ КРИСТАЛЛИЗАЦИИ КАК РЕЗЕРВ ПОВЫШЕНИЯ КАЧЕСТВА ПОРШНЕЙ ДВС

В статье предложено использование процедуры распознавания образов для описания локализации дефектов в литых деталях «поршень» для ДВС. Приведены результаты компьютерного моделирования процесса заполнения кокиля и кристаллизации сплава для усовершенствованной технологии литья, предполагающей применение термоизолирующих покрытий кокиля. Результаты компьютерного моделирования могут быть использованы для формализации процесса описания локализации дефектов в поршнях.

Ключевые слова: литая деталь, кокиль, компьютерное моделирование, термоизолирующее покрытие.

1. Введение

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

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

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

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

2. Постановка задачи исследования

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

3. Технологический процесс изготовления поршней

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

Поршни, изготавливаемые литьем, в основном получают литьем в кокиль. Этим методом производится большинство автомобильных и легких промышленных поршней. Внутренняя поверхность поршня формируется стержнем, который может быть цельным (моноклин) или состоять из 3—7 частей (клиньев). Наружная поверхность поршня формируется матрицей, образованной несколькими частями (обычно двумя), иногда еще и верхним стержнем. Собранные вместе все вышеперечисленные компоненты образуют литейную форму — кокиль.

После литья материал заготовки не обладает достаточной прочностью, твердостью и стабильностью при высоких температурах. Для улучшения структуры материала и придания необходимых механических свойств заготовку подвергают термообработке. После остывания отливки проходят этап так называемого «старения». Их помещают в печь, где нагревают до температуры примерно 520 °С, затем отливки медленно остывают вместе с печью. В зависимости от требований, предъявляемых к изделию, далее может иметь место один из двух вариантов:

1. Поршням дают остыть до комнатной температуры, отрезают литники и выпоры (прибыли) и загружают в печь при температуре около 515 °С на 4 часа. Затем поршни закаляют в кипящей воде, которая их быстро

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

2. После удаления из пресс-формы отливки немедленно обдуваются воздухом (воздушная закалка). Это достаточно быстро их охлаждает, чтобы исключить выделение осадка и связать (заморозить) медь и магний в твердом растворе. Этот процесс называется «Воздушная закалка после пресс-формы» и позволяет снизить энергозатраты и подготовительно-наладочное время, благодаря устранению обработки на твердый раствор. Таким методом сейчас производят 80 % поршней.

При обоих вариантах обработки поршни затем подвергаются процессу старения (примерно 4 часа при температуре около 220 °С), в результате которого образуется контролируемая структура материала по всему поршню.

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

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

4. Материалы для изготовления поршня

Сплавы алюминия с некоторым количеством легирующих элементов (кремний, медь, магний, никель и др.) Кремний ^) в количестве 5—22 % придает сплаву жаропрочность, уменьшает коэффициент теплового расширения, улучшает литейные свойства, но при этом охрупчивает сплав и ухудшает обрабатываемость. Медь (Си) 1—8 % улучшает механические свойства материала. Магний (Мg) 0,8—1,3 % придает сплаву высокиемеханические характеристики, повышая предел прочности. Никель (№) 0,5—2 % придает сплаву жаропрочность и жаростойкость. Титан (Т^ 0,05—2 % улучшает механические свойства материала, за счет создания с алюминием интер-металлидных соединений уменьшает ползучесть. Различают три типа алюминиевых сплавов, применяемых в производстве поршней, которые отличаются главным образом содержанием кремния.

Доэвтектические сплавы. Содержание кремния ниже, чем в точке эвтектики — обычно около 9 % (для сплавов алюминий-кремний точка эвтектики 11—13 %).

Эвтектические сплавы. Содержание кремния составляет 11—13 % (в зависимости от количества других легирующих элементов).

Заэвтектические сплавы. В этих сплавах содержится значительное количество нерастворенного кремния в алюминии, что обеспечивает повышение износостойкости, но ухудшает технологические свойства. Содержание кремния может колебаться от 13 до 25 %.

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

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

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

5. Модернизация технологии изготовления литого поршня 51-03-23 на основе результатов компьютерного моделирования

В качестве объекта исследования была выбрана конструкция кокиля и технология изготовления литых поршней на ОАО «АВТРАМАТ» (г. Харьков, Украина). Некоторые результаты моделирования термонапряженного состояния литого поршня, зависящего от конструкции кокиля и принятых технологических режимов литья, были рассмотрены в работах [7, 8]. После анализа приведенных в этих работах результатов для литого поршня модификации 51-03-23 и определения вероятности формирования в его теле усадочной пористости, было установлено, что существующая технология изготовления поршнем на ОАО «АВТРАМАТ» имеет ряд недостатков. Эти недостатки существующего технологического процесса приводят к формированию усадочной пористости и недостаточной плотности структуры сплава в ответственных местах отливки в силу слишком большой скорости кристаллизации. На основе этого было выбрано направление проведения дальнейших исследований на основе компьютерного моделирования процесса кристаллизации литого поршня. В качестве технического решения, обеспечивающего по предварительным выводам возможность регулирования процесса кристаллизации, было принято использование теплоизолирующего кожуха, свойства которого приведены в табл. 1, и предназначенного для поддержания сплава в жидком состоянии в течение более длительного времени.

Теплоизолирующий кожух, укрывающий кокиль, был смоделирован с помощью программы Solid Works и поршень был «пролит» с помощью пакета прикладных программ LVMFlow.

Результаты моделирования термонапряженного состояния поршня и вероятности образования усадочной пористости представлены на рис. 1.

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

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

Рис. 1. Зоны внутренних дефектов: а, б — зоны вероятного образования усадочной пористости; в, г — зоны пониженной плотности сплава

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

W

Рис. 2. Результаты компьютерного моделирования мест формирования внутренних дефектов

Таким образом, оказалось возможным поставить задачу о формализации процесса локализации внутренних дефектов.

Таблица 1

Свойства материала теплоизолирующего кожуха

в

г

Характеристика теплоизолирующего материала Природа материала Плотность, кг/м3 Коэффициент теплопроводности, Вт/м*К Прочность на сжатие кгс/см2 Коэффициент теплового линейного расширения

Блоки и плиты из пеностекла FOAMGLAS® Неорганическое алю-мосиликатное стекло 120—160 0,045—0,06 5—20 9 * 10-6 (Ді = 0,27 мм)

6. Формализация процесса локализации внутренних дефектов

Идея применения метода распознавания образов для формализации процедуры локализации иллюстрирует рис. 3.

б

Рис. 3. Схема, поясняющая идею применения метода распознавания образов для формализации процедуры локализации внутренних дефектов в детали: а — двухмерное представление; б — трехмерное представление

f-2 Pa(X ) = ke^

f--

Pa(X) = ke^ 2

X - АЛ I cov 1(X )l A - A

-1 (x -(x

(1)

где k — постоянный множитель, СОУА(Х) и СОУВ(Х) — ковариационные матрицы классов А и В соответственно.

3. Расчет условных вероятностей принадлежности j-го объекта (координата точки отливки) к классу А и В:

P{([xi ) =

P(A) Pa(Xі )

P (A)pa( Xі ) + P (A) Pa ( Xі )’

P Axi ) =_________P (A)pB(Xi)_____-.

^ 1 > P (A)pa( x! ) + P (A) Pa (x!)

(2)

4. Получение классифицирующего правила в виде (3) или (4):

х1 є А если Р[А\х])> Р[-В|хі), х1 єВ если Р(АХ)<Р(х1).

х1 єА если F(х1) у0, х1 єВ если F{АХ)<у0.

(З)

(4)

где F(xj) — дискриминантная функция, определяемая по уравнению (5), у0 — пороговое значение дискриминантной функции.

F (x ) = x' cov 1 (x )(mA - AB ) =

1A = 21 AA + A

(mA + mA) cov 1 (x)(mA -mA)-ln^(5)

В формуле (5) Р(А) и Р(В) — вероятности наблюдения классов А и В соответственно, а координатами вектора Х выбираются значения координат XYZ, снимаемые с результатов компьютерного моделирования.

Пример применения данного метода для одного из описаний плоскости для рассматриваемой литой детали «поршень» приведен в форме классифицирующего правила (6):

Разделяющие прямые — это границы, описываемые математически, и отделяющие область существования внутренних дефектов от бездефектной области и эти границы могут быть получены на основе применения методов распознавания образов, в частности методов статистической классификации [6, 9, 10]. Согласно этим работам, для решения задачи классификации должны быть выполнены следующие процедуры.

1. Расчет математических ожиданий тА и тв классов А и В, определяющих области существования дефектов в отливке от бездефектной области.

2. Расчет плотности распределения вероятностей классов А и В — рА(Х) и рВ(Х) по формулам:

Если тА и тв — математические ожидания вектора Х для классов А и В, и ковариационные матрицы распределения вектора Х для классов А и В равны (соуа(Х) = соув(Х)), плотности распределения вероятностей рА(Х) и рВ(Х) могут быть представлены в виде:

х* е А если [-7,803х1 -51,025х2 -7,575х3 > (-3510,79)], х* еВ если [-7,803[ -51,025х2 -7,575х3 < (-3510,79)]. (6)

На рис. 4 показано распределение для классов А и В в задаче построения поверхности, разделяющей дефектную и бездефектную области.

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

8. Савченко, Ю. Применение компьютерно-интегрированных систем и технологий в производстве поршней [Текст] / Ю. Савченко // Технологический аудит и резервы производства. — 2012. — Т. 1, № 1(3). — С. 8—13.

9. Дёмин, Д. А. Применение параметрических методов распознавания образов для исследования длительности модифицирующего эффекта [Текст] : материалы Научно-технической конференции / Д. А. Дёмин // Информационные технологии: наука, техника, технология, образование, здоровье. — 1996. — С. 56—58.

10. Дёмин, Д. А. Идентификация чугуна для определения рациональных режимов легирования [Текст] / Д. А. Дёмин, А. Б. Божко, А. В. Зрайченко, А. Г. Некрасов // ВосточноЕвропейский журнал передовых технологий. — 2006. — № 4/1(22). — С. 29—32.

11. Дёмин, Д. А. Обработка экспериментальных данных и построение математической модели технологического процесса методом наименьших квадратов (МНК) [Текст] / Д. А. Дёмин // Восточно-Европейский журнал передовых технологий. — 2006. — № 3/1(26). — С. 47—50.

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

12. Дёмин, Д. А. Оптимизация технологического процесса в цехе предприятия [Текст] / Д. А. Дёмин // Восточно-Европейский журнал передовых технологий. — 2005. — № 6/1(18). — С. 48—59.

13. Дёмин, Д. А. Оптимизация технологических режимов [Текст] / Д. А. Дёмин // Восточно-Европейский журнал передовых технологий. — 2006. — № 2/1(20). — С. 32—35.

КОМП'ЮТЕРНЕ МОДЕЛЮВАННЯ ПРОЦЕСІВ КРИСТАЛІЗАЦІЇ ЯК РЕЗЕРВ ПІДВИЩЕННЯ ЯКОСТІ ПОРШНІВ ДВЗ

У статті запропоновано використання процедури розпізнавання образів для опису локалізації дефектів у литих деталях «поршень» для ДВС. Наведено результати комп’ютерного моделювання процесу заповнення кокілю і кристалізації сплаву для вдосконаленої технології лиття, що передбачає застосування термоізолюючих покриттів кокіль. Результати комп’ютерного моделювання можуть бути використані для формалізації процесу опису локалізації дефектів у поршнях.

Ключові слова: лита деталь, кокіль, комп’ютерне моделювання, термоізолююче покриття.

Пономаренко Ольга Ивановна, доктор технических наук, профессор, кафедра литейного производства, Национальный технический университет «Харьковский политехнический институт», Украина, e-mail: [email protected].

Тренёв Николай Сергеевич, кафедра литейного производства, Национальный технический университет «Харьковский политехнический институт», Украина, e-mail: [email protected].

Пономаренко Ольга Іванівна, доктор технічних наук, професор, кафедра ливарного виробництва, Національний технічний університет «Харківський політехнічний інститут», Україна. Треньов Микола Сергійович, кафедра ливарного виробництва, Національний технічний університет «Харківський політехнічний інститут», Україна.

Ponomarenko Olga, National Technical University «Kharkiv Polytechnic Institute», Ukraine, e-mail: [email protected]. Trenev Nicholay, National Technical University «Kharkiv Polytechnic Institute», Ukraine, e-mail: [email protected]

ства. - 2012. - T 1, № 1(3). - С. 3-8.

□ Распределение для класса А (бездефектная область)

□ Распределение для класса В (область с наличием дефектов)

Значение критерия - величины y

Рис. 4. Зависимость частоты попадания в соответствующий интервал от значения критерия — величины у

7. Выводы

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

Литература

1. Пономаренко, О. И. Оптимизация технологических решений в условиях работы литейных цехов [Текст] / О. И. Пономаренко. — Харьков: НТУ «ХПИ», 2007. — 320 с.

2. Алёхин, В. И. Компьютерное моделирование процессов при производстве литых деталей двигателя [Текст] / В. И. Алёхин, О. В. Акимов, А. П. Марченко // Литейное производство. — Москва, 2010. — № 9. — С. 31—33.

3. Алёхин, В. И. Компьютерно-интегрированное моделирование литейных процессов в автомобильных поршнях на основе конструкторско-технологической методики проектирования деталей ДВС [Текст] / В. И. Алёхин, А. В. Белогуб, А. П. Марченко, О. В. Акимов // Двигатели внутреннего сгорания. — Харьков: НТУ «ХПИ», 2009. — № 2. — С. 101 — 104.

4. Алёхин, В. И. Моделирование литейных процессов при изготовлении автомобильных поршней [Текст] / В. И. Алёхин, А. В. Белогуб, А. П. Марченко, О. В. Акимов // Цветные металлы. — Москва, 2010. — № 8. — С. 81—83.

5. Алёхин, В. И. Моделирование мест проявления дефектов усадочного характера при проектировании литых деталей ДВС [Текст] / В. И. Алёхин, А. В. Белогуб, А. П. Марченко, О. В. Акимов // Металл и литье Украины. — Киев, 2010. — № 12. — С. 27—30.

6. Дёмин, Д. А. Диагностика технологического процесса. Руководство для технолога [Текст] / Д. А. Дёмин. — ВосточноЕвропейский журнал передовых технологий. — 2005. — № 5/1(17). — С. 29—40.

7. Краснокутский, Е. Компьютерное моделирование процессов кристаллизации литой детали в кокиле [Текст] / Е. Крас-нокутский // Технологический аудит и резервы производ-

7

6

5

з

2 -

1 -

о

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