Научная статья на тему 'Учебно-консультационная компьютерная программа для изучения теплообменных процессов'

Учебно-консультационная компьютерная программа для изучения теплообменных процессов Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
258
70
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МНОГОСЛОЙНАЯ СТЕНКА / ГАРНИСАЖ / ТЕПЛОВОЙ ПОТОК / ТЕПЛОПРОВОДНОСТЬ / ГРАНИЧНЫЕ УСЛОВИЯ

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Портянкин Артем Александрович, Тинькова Светлана Михайловна, Пискажова Татьяна Валериевна

Специальные компьютерные программы повышают качество обучения за счет проведения большего количества операций за короткий промежуток времени, структурирования учебной информации на разных уровнях, систематизации процесса представления информации и организации интерактивного общения. Целью работы является создание учебно-консультационного комплекса для преподавания курса теплотехники при подготовке специалистов горно-металлургической отрасли и решения типовых инженерных задач. Программа имеет следующие функциональные возможности: выполнение статических расчетов: теплового потока и распределения температур по толщине плоской стенки, расчет толщины слоя стенки (при задании различных граничных условий); выполнение динамических расчетов: перераспределения температур, плавление и кристаллизация гарнисажа как одного из слоев стенки, при подаче различных воздействий; использование встроенных баз данных, необходимых для проведения расчетов и проверки их правильности; визуализация результатов в виде графиков и таблиц. Компьютерная программа позволяет с помощью цветовых анимационных иллюстраций и графического изображения наглядно представить полученные результаты, предусматривает возможность интерактивного взаимодействия между пользователем и элементами программы, при этом обеспечивает мгновенный контроль за усвоением материала. Данная программа может быть использована также и для дистанционного обучения. Программное обеспечение было реализовано в среде объектно-ориентированного программирования C++ Builder. При создании программы были использованы: теория теплопередачи, численные методы для решения обыкновенных дифференциальных уравнений. Представленная программа опробована в процессе обучения бакалавров и магистрантов института цветных металлов и материаловедения СФУ по направлениям «Металлургия», «Автоматизация технологических процессов и производств», «Управление в технических системах».

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

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Портянкин Артем Александрович, Тинькова Светлана Михайловна, Пискажова Татьяна Валериевна

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

Training-support software for studying heat exchange processes

Special software improve the quality of training by processing many transactions in a short period of time, structuring learning information at different levels, organizing presentation of information, and managing interactive communication. The objective is to develop a training-support complex for teaching a course on Thermal Engineering when training specialists for mining and metallurgical industry and solving common engineering problems. The software has the following features: static calculation of heat flow and temperature distribution over the thickness of a flat side-wall, calculation of the side-wall layer thickness (with various boundary conditions); dynamic calculations of temperature redistribution, melting and side-ledge crystallization as one of the side-wall layers, when applying various impacts; use of embedded databases required to perform and verify calculations; Graphical and tabular rendering of the results. The software makes it possible to render the results obtained using color illustrations and animated graphic, provides for conversational interaction between the user and software components, and at the same time, ensures instant control of mastering the material delivered. It can be used for distant learning as well. The software has been implemented in the C ++ Builder environment. The theory of heat transfer and numerical methods of solving ordinary differential equations were used to design the software. The proposed software was trialed during training undergraduate and graduate students of the Institute of Nonferrous Metals and Materials Science of the Siberian Federal University majoring in Metallurgy, Automation of Technological Processes and Production, Management in Technical Systems.

Текст научной работы на тему «Учебно-консультационная компьютерная программа для изучения теплообменных процессов»

УДК 681.5 : 536 DOI:10.18503/1995-2732-2016-14-1-116-123

УЧЕБНО-КОНСУЛЬТАЦИОННАЯ КОМПЬЮТЕРНАЯ ПРОГРАММА ДЛЯ ИЗУЧЕНИЯ ТЕПЛООБМЕННЫХ ПРОЦЕССОВ

Портянкин А.А., Тинькова С.М., Пискажова Т.В.

Сибирский федеральный университет, Красноярск, Россия

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

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

Программа имеет следующие функциональные возможности:

- выполнение статических расчетов: теплового потока и распределения температур по толщине плоской стенки, расчет толщины слоя стенки (при задании различных граничных условий);

- выполнение динамических расчетов: перераспределения температур, плавление и кристаллизация гарнисажа как одного из слоев стенки, при подаче различных воздействий;

- использование встроенных баз данных, необходимых для проведения расчетов и проверки их правильности;

- визуализация результатов в виде графиков и таблиц.

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

Программное обеспечение было реализовано в среде объектно-ориентированного программирования C++ Builder. При создании программы были использованы: теория теплопередачи, численные методы для решения обыкновенных дифференциальных уравнений.

Представленная программа опробована в процессе обучения бакалавров и магистрантов института цветных металлов и материаловедения СФУ по направлениям «Металлургия», «Автоматизация технологических процессов и производств», «Управление в технических системах».

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

Введение

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

© Портянкин А.А., Тинькова С.М., Пискажова Т.В., 2016

Эту задачу можно решить путем создания специализированного программного обеспечения и его внедрения в процесс подготовки будущих специалистов горно-металлургического комплекса.

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

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

программы ЗБ-моделирования SolidWorks [1], Лп8у8 [2], в которых можно выполнить большую часть интересующих нас теплотехнических расчетов, но для овладения этими программами нужно долго учиться, и расчеты в них длительные, что не всегда дает возможность быстро оценить результаты вводимого в объект изменения. Существует еще определенный класс учебно-консультационных программ, которые имитируют производственные процессы и объекты уже непосредственно для нужд производства - повышения квалификации, переподготовки, тренировки специалистов, а также для выполнения производственных расчетов [3, 4].

Целью нашей программы является:

- повышение уровня подготовки молодых специалистов;

- сокращение временных затрат при обучении на повторяющиеся вычислительные операции;

- визуализация полученных результатов;

- создание у обучающихся четкого представления о работе исследуемого объекта и понимания выполняемых расчетов.

Описание программного обеспечения

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

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

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

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

лабораторных работ с использованием описываемого ПО:

• расчеты теплового потока и распределения температур по толщине плоской многослойной стенки при задании граничных условий 1-го и 3-го рода;

• расчет толщины третьего слоя стенки при задании граничных условий 1-го рода -выбор теплоизоляции;

• расчет толщины гарнисажа на внутренней стенке при задании граничных условий 3-го рода;

• динамический расчет изменения распределения температур слоев стенки при задании граничных условий 3-го рода;

• динамический расчет толщины гарнисажа на многослойной стенке при задании граничных условий 3-го рода.

Программа состоит из трех основных частей:

1. База данных по огнеупорным и теплоизоляционным материалам.

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

3. Интерфейс пользователя, позволяющий вводить условия и видеть результаты расчетов.

Базы данных

Для задания теплофизических характеристик создана расширяемая база данных, которая позволяет осуществлять подбор материала для конструкции при задании различных исходных данных. Фрагмент базы представлен на рис. 1, где отражены теплофизические характеристики огнеупорных и теплоизоляционных материалов, такие как коэффициент теплопроводности X, удельная теплоемкость с и плотность материалов Я.

Коэффициент теплопроводности определяется по уравнению

Л = А + В ■ Т,

где А и В - константы для данного материала; Т -температура материала, °С.

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

с = С1 + С2 • Т,

где С1 и С2 - коэффициенты для определения теплоемкости; Т- температура слоя, °С.

Рис. 1. Фрагмент базы данных

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

Математическое обеспечение

Одной из задач является определение тепловых потоков через стенку. Если рассматривать металлургические печи, то величина этих потоков влияет на энергозатраты. Расчет тепловых потерь основан на уравнении теплопередачи [5]

д = к-ДТ,

где АТ - разница температур; к = ^^Я - коп

эффициент теплопередачи; ^ Я - сумма тепло-

1=1

вых сопротивлений системы.

Для определения распределения температур в многослойной стенке при стационарных условиях разработана программа, которая базируется

на методе «холодных сопротивлений» при стационарных условиях [6].

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

При определении коэффициентов теплопроводности в зависимости от температуры материала используется итеративная процедура.

На рис. 3 представлен интерфейс программного продукта. Пользователь имеет возможность задавать теплофизические параметры стенки путем выбора материала из универсальной базы данных, различные граничные условия и размеры слоев стенки. Также интерфейс дает возможность графической визуализации и табличного представления результатов расчетов. На графике изображен температурный профиль в многослойной стенке, по горизонтальной оси представлена заданная толщина слоя стенки, по вертикальной - температура, °С, Т2 и Т3 - значения температур на границах слоев.

Материалы 1-3 слоев; Температуры: Внешней поверхности стенки (14), Внутренней поверхности стенки (11);

Толщина слоев (б^^з); Подключение базы данных (А^АгАз);

сК=М4;

Р=51/А1 +52/А2+5з/Аз=Ро1+^2+ ^з; А^сИ-б^-А!);

Задаются материалы слоев; Задается температура внешней поверхности стенки; Задается температура внутренней поверхности стенки; Задается толщина слоев; Определяются коэффициенты теплопроводности;

Определяется предполагаемая температура на границах слоев и проверяется температура внешней __поверхности стенки

Л2=А2+В2-(12+13)/2 Хз=Аз+Вз-аз+14)/2

Определяются коэффициенты теплопроводности

q=(t1-t4)/(5Дl+A2/)l2+AЗУЗlз)=(t1-t4)/(R1+R2+Rз);

Определяется тепловой поток

Ч; Хь Ь', Ь ▼ ^

Конец

Вывод результатов

Рис. 2. Блок-схема для расчета теплового потока и распределения температуры

по толщине плоской стенки

Рис. 3. Главная страница программы с исходными данными и результатами расчета

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

Изучение нестационарных процессов в курсе теплотехники требует особой наглядности и визуализации, является предметом для научной работы студентов и углубленного изучения курса.

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

МТ п п т •с •—=а* - &ых.

м

При этом средние температуры в слоях рассчитываются следующим образом:

МТ

т • с •—1 = 5• к • (Т -Т)-5• к_2 • (Т -т2);

МТ

т2 ■ с2--2 = 5 • к12 • (Т -Т2) - 5 • к2 3 • (Т2 -Т3);

МТ

т3 ■ с3 ■ ~т3 = 5 • к2-3 • (Т2 - Т3) - 5 • аенеш • (Т3 - Теоэд X

где Qвx- входные тепловые потоки; Qвыx - выходные тепловые потоки; тп - масса слоев стенки; с12,3 - удельная теплоемкость соответствующих

мТ

материалов стенки; - - скорость изменения

М

средней температуры слоя; 5 - площадь поверхности стенки; Т - температура внутренней поверхности стенки; Т12 3 - средние температуры соответствующих слоев; Овнеш- коэффициент теплоотдачи во внешнюю среду; Т3 - температура границы между вторым и третьим слоем; Твозд - температура внешней среды; к!,2,3 - коэффициенты теплопередачи, которые рассчитываются по следующим формулам [7]:

к = 2 • ^

к1-2 - /

2

5к+

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

V ^ 2

к -_

2-3 / ,

2

Л'

V ^2 ^3 У

Здесь ^1,2,3 - коэффициент теплопроводности материала соответствующего слоя; 51,2,3 - толщина соответствующего слоя.

Значения температур на границах слоев определяются из условий равенства температур и потоков на границе раздела слоев.

На рис. 4 представлены результаты расчета динамического изменения средних температур слоев и температуры границ слоев стенки в результате резкого увеличения температуры внутренней поверхности стенки (например, с температуры 1700 до 1900°С), что реально может наблюдаться при работе ряда металлургических печей.

Рис. 4. Интерфейс для отображения расчетов динамического изменения температур стенки

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

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

При управлении некоторыми металлургическими печами актуальной задачей является определение толщины гарнисажа [8]. Гарнисаж - твёрдый защитный слой, образующийся при плавке на внутренней (рабочей) поверхности стенок некоторых металлургических агрегатов, подвергающихся интенсивному охлаждению. При этом образуются тонкие корки застывших продуктов, толщина которых зависит от условий плавки и, особенно, от скорости охлаждения печи. Нами был реализован динамический расчет толщины гарнисажа при задании граничных условий третьего рода, с использованием следующего уравнения [7]:

(

8п+1 = 0,5 •

5„ --

' авнут • (Т0 - Т)

рн • К

0,25 •

М • а

^ внут

• (Та - Т)

Р« • К

2-Л • X н • (Т - Т) Рн •К

где 5п - толщина гарнисажа в данный момент времени; авнут - коэффициент внутренней теплоотдачи; Хв, рм - коэффициенты теплопроводности и плотность гарнисажа соответственно; Кн - удельная теплота плавления гарнисажа; Т0 -температура внутренней среды; Т

- температура ликвидуса расплава; Т1 - средняя температура слоя гарнисажа; М - шаг по времени.

На рис. 5 представлен результат по определению динамического изменения толщины гарнисажа при задании условий 3-го рода и подаче воздействия путем увеличения температуры ликвидуса расплава на 5 град.

Выводы

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

Рис. 5. Интерфейс для отображения результатов расчета изменения толщины гарнисажа

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

1. Официальный сайт SolidWorks [Электронный ресурс]. URL: http://www.solidworks.ru (дата обращения: 06.03.2014).

2. Официальный сайт Ansys [Электронный ресурс]. URL: http://www.ansys.com (дата обращения: 06.03.2014).

3. Piskazhova T.V., Mann V.C. The Use of a Dynamic Aluminum Cell Model // JOM. 2006. Vol. 58, №2. P. 48-52.

4. Информационно-аналитический журнал Rational Enterprise Management // Имитационное моделирование как инструмент оптимизации производственных процессов в металлургии [Электронный ресурс]. URL: http://www.anylogic.ru/upload/iblock/e56/ e56ccf70ee38f9080c9bb7f69f2b5908.pdf (дата обращения 07.04.2014).

5. Теплотехника металлургического производства. Т. 1: Теоретические основы: учеб. пособие для вузов / В.А. Криван-дин [и др.]. М.: МИСиС, 2002.

6. Тинькова С.М. Металлургическая теплотехника. Красноярск: Государственный университет цветных металлов и золота, 2005. 143 с.

7. Белолипецкий В.М., Пискажова T.B. Математическое моделирование процесса электролитического получения алюминия. Решение задач управления технологией. Красноярск: Сибирский федеральный университет, 2013. 271 с.

8. Chuck C., Chen J.J., Barry J. Modeling of Dynamic Ledge Heat Transfeer // Light Metals. 1997. P. 309-317.

9. Арутюнов B.A., Бухмиров B.B., Крупенников C.A. Математическое моделирование тепловой работы промышленных печей. М.: Наука, 1990. 319 с.

10. Портянкин А.А., Тинькова С.М., Пискажова T.B. // Конференция «Молодёжь и наука» [Molodezh i nauka] 2013г. [Электронный ресурс]. URL: http://conf.sfu-kras.ru/sites/mn2013/thesis/ s008/s008-007.pdf (дата обращения: 06.03.2014).

Материал поступил в редакцию 20.07.15.

INFORMATION ABOUT THE PAPER IN ENGLISH

D0I:10.18503/1995-2732-2016-14-1-116-123

TRAINING-SUPPORT SOFTWARE FOR STUDYING HEAT EXCHANGE PROCESSES

Portyankin Artem Aleksandrovich - Postgraduate Student, Siberian Federal University, Krasnoyarsk, Russia. E-mail: [email protected].

Tinkova Svetlana Mikhailovna - Ph.D. (Eng.), Associate Professor, Siberian Federal University, Krasnoyarsk, Russia.

Piskazhova Tatiana Valerievna - D.Sc. (Eng.), Siberian Federal University, Krasnoyarsk, Russia. E-mail: [email protected].

Abstract. Special software improve the quality of training by processing many transactions in a short period of time, structuring learning information at different levels, organizing presentation of information, and managing interactive communication.

The objective is to develop a training-support complex for teaching a course on Thermal Engineering when training specialists for mining and metallurgical industry and solving common engineering problems.

The software has the following features:

- static calculation of heat flow and temperature distribution over the thickness of a flat side-wall, calculation of the side-wall layer thickness (with various boundary conditions);

- dynamic calculations of temperature redistribution, melting and side-ledge crystallization as one of the side-wall layers, when applying various impacts;

- use of embedded databases required to perform and verify calculations;

- Graphical and tabular rendering of the results.

The software makes it possible to render the results obtained using color illustrations and animated graphic, provides for conversational interaction between the user and software components, and at the same time, ensures instant control of mastering the material delivered. It can be used for distant learning as well.

The software has been implemented in the C ++ Builder environment. The theory of heat transfer and numerical methods of solving ordinary differential equations were used to design the software.

The proposed software was trialed during training undergraduate and graduate students of the Institute of Nonferrous Metals and Materials Science of the Siberian Federal University majoring in Metallurgy, Automation of Technological Processes and Production, Management in Technical Systems.

Keywords: Multi-layered wall, ledge, heat flow, thermal conductivity, boundary conditions.

References

1. The SolidWorks official website [electronic resource]. URL: http://www.solidworks.ru (accessed date: 03.06.2014).

2. The Ansys official website [electronic resource]. URL: http://www.ansys.com (accessed date: 03.06.2014).

3. Piskazhova T.V., Mann V.C. The Use of a Dynamic Aluminum Cell Model. JOM. 2006, vol. 58, no. 2, pp. 48-52.

4. Simulation modeling as a tool for optimization of production processes in the steel industry [electronic resource]. URL: http://www.anylogic.ru/upload/iblock/e56/e56ccf70ee38f9080c9bb 7f69f2b5908.pdf (accessed date: 04.07.2014).

5. Krivandin V.A. Teplotekhnika metallurgicheskogo proizvodstva. Tom 1: Teoreticheskie osnovy: uchebnoe posobie dlya vuzov [Heat engineering of metals practice. Vol. 1. Theoretical Foundations: Textbook for Universities]. Moscow: MISiS, 2002.

Tinkova S.M. Metallurgicheskaya teplotekhnika [Metallurgical heat engineering]. Krasnoyarsk: State University of Non-ferrous Metals and Gold, 2005. 143 p.

Belolipetsky V.M., Piskazhova T.V. Matematicheskoe modeliro-vanie protsessa ehlektroliticheskogo polucheniya alyuminiya. Reshenie zadach upravleniya tekhnologiej [Mathematical modeling of electrolytic production of aluminum. Meeting the challenges of technology management]. Krasnoyarsk: Siberian Federal University, 2013. 271 p.

10.

Chuck C., John J.J. Chen, Barry J. Modeling of Dynamic Ledge Heat Transfeer. Legkie metally [Light Metals]. 1997, pp. 309-317. Arutyunov V.A., Bukhmirov V.V., Krupennikov S.A. Mathematical modeling of thermal performance of industrial furnaces. Moscow: Nauka, 1990. 319 p.

Portyankin A.A., Tinkova S.M., Piskazhova T.V. "Youth and Science" Conference [Molodiozh i nauka] 2013. [Electronic resource]. URL: http://conf.sfu-kras.ru/sites/mn2013/thesis/s008/s008-007.pdf (accessed date 06.03.2014).

Портянкин A.A., Тинькова C.M., Пискажова T.B. Учебно-консультационная компьютерная программа для изучения теплообменных процессов // Вестник Магнитогорского государственного технического университета им. Г.И. Носова. 2016. Т. 14. №1. С. 116-123. doi:10.18503/1995-2732-2016-14-1-116-123

Portyankin A.A., Piskazhova T.V., Tinkova S.M. Training-support software for studying heat exchange processes. Vestnik Magnitogorskogo Gosudarstvennogo Tekhnicheskogo Universiteta im. G.I. Nosova [Vestnik of Nosov Magnitogorsk State Technical University]. 2016, vol. 14, no. 1, pp. 116-123. doi:10.18503/1995-2732-2016-14-1-116-123

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