CC BY
270
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
УДК 691.322:004.421 Дата подачи статьи: 29.08.2014
DOI: 10.15827/0236-235X.108.254-259
РАЗРАБОТКА ПРОГРАММНО-АЛГОРИТМИЧЕСКОГО СРЕДСТВА ОБРАБОТКИ ДАННЫХ ТРЕХФАКТОРНОГО ПЛАНИРОВАННОГО ЭКСПЕРИМЕНТА ДЛЯ РАСЧЕТА МАТЕМАТИЧЕСКОЙ МОДЕЛИ
ПРОЧНОСТИ БЕТОНА
В.В. Белов, д.т.н., профессор, проректор, зав. кафедрой, [email protected]; И.В. Образцов, инженер, [email protected];
Ю.Ю. Курятников, к.т.н., доцент, yuriy-k@)yandex.ru (Тверской государственный технический университет, наб. Аф. Никитина, 22, г. Тверь, 1 70026, Россия)
В технологии строительных композиционных материалов наряду с субстанционным (изготовление физического образца материала) и структурно-имитационным (имитация взаимодействия структурных элементов системы) моделированием широко применяется функциональное моделирование, результатом которого является получение некой математической функции, описывающей поведение объекта исследования, абстрагируясь от внутренней структуры вещественного субстрата.
Авторами разработана программа обработки данных трехфакторного планированного эксперимента B-D13, позволяющая рассчитывать план и обрабатывать результаты эксперимента. Алгоритм программы включает процедуру расчета коэффициентов функции отклика, статистической обработки и визуализации математической модели. Все основные вычисления производятся циклично, что позволяет моментально перестраивать математическую модель, изменяя входные данные. Кроме того, алгоритм включает вспомогательную процедуру, обеспечивающую проверку синтаксической правильности вводимых данных. При допущении ошибок ввода данных программа корректирует действия пользователя по средствам текстового оповещения. Программа оснащена дополнительными функциями загрузки и сохранения данных, а также функцией экспорта результатов расчетов в математический пакет Microsoft Excel, где пользователь может произвести дополнительные построения. Интерфейс программного продукта реализован в виде логических блоков, позволяющих вводить исходные данные и изменять параметры вывода математической модели в интерактивном режиме.
С применением разработанного программного продукта проведен планированный эксперимент по исследованию влияния рецептурных параметров бетонной смеси на свойства наполненного мелкозернистого бетона, получена статистически адекватная математическая модель прочности бетона. Разработанное программное средство можно применять в любых научно-прикладных задачах по оптимизации свойств объекта исследования, подбора рецептуры и технологических параметров, где используется математическое моделирование методом ортогонального планирования экспериментов.
Ключевые слова: функциональная модель, планированный эксперимент, трехфакторный план, объект исследования, функция отклика.
В настоящее время в научных исследованиях и образовании, в производственной и других сферах деятельности человека определяющее значение имеют информационно-вычислительные системы. Развитие информатики и применение программных продуктов в научных исследованиях ставят вопрос о пересмотре основных концепций представления научных знаний и выдвигают на первый план задачу структурирования этих знаний [1]. Электронные образовательные ресурсы на основе современной компьютерной трехмерной симуляции физических процессов и явлений реализуются в форме мультимедийных учебно-научных лабораторий, или виртуальных тренажеров. Новизна технологии виртуальных тренажеров аргументируется использованием современных средств компьютерного моделирования и активным внедрением информационных технологий в сферу образования как нового трансдисциплинарного направления [2]. Именно с этой точки зрения внедрение информационных технологий способствует оптимальному решению вышеназванных задач и устранению ряда недостатков традиционного
способа обучения. Эти вопросы во всей полноте можно решать с помощью мультимедийных учебно-научных лабораторий [3, 4].
Мультимедийная учебно-научная лаборатория, как правило, сочетает в себе имитационную динамическую модель оборудования и программную оболочку с методическим сопровождением лабораторной работы [5]. Программные средства - это математически обоснованная виртуальная модель, включающая систему графической визуализации, звуковое сопровождение и текстовую информацию [6]. Ввод и вывод информации осуществляются согласно разработанному алгоритму -программному коду виртуальной модели. Современные виртуальные тренажеры выполняются в трехмерной графике с максимальной имитацией материалов и освещения, что существенно повышает качество работы. Главной задачей здесь является приближение модели к реальному объекту за счет соблюдения правильных пропорций, размеров, цветовых решений и освещения [7, 8].
Моделирование является быстро развивающейся областью науки и техники. В технологии
254
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
строительных композиционных материалов наряду с субстанционным (изготовление физического образца материала) и структурно-имитационным (имитация взаимодействия структурных элементов системы) моделированием широко применяется функциональное моделирование, результатом которого является получение некой математической функции, описывающей поведение объекта исследования, абстрагируясь от внутренней структуры вещественного субстрата. Функциональная модель работает по принципу «черного ящика», при этом известны параметры входа - переменные или постоянные факторы, а также параметры выхода - критерий эффективности, отклик и т.д. [9, 10].
Построение функциональных моделей экспериментальных зависимостей свойств бетона от его состава включает следующие этапы:
- уточнение оптимизируемых параметров (прочности бетона, удобоукладываемости бетонной смеси и др.) в зависимости от конкретной задачи;
- выбор факторов, определяющих изменчивость оптимизируемых параметров;
- определение основного исходного состава бетонной смеси;
- выбор интервалов варьирования факторов;
- выбор плана и условий проведения эксперимента;
- обработка результатов эксперимента с построением математических моделей зависимостей свойств бетонной смеси и бетона от выбранных факторов.
В данном исследовании выбран трехфакторный план эксперимента B-D13, который позволяет получать нелинейные квадратичные модели и обладает хорошими статистическими характеристиками. Основная цель разработки программно -алгоритмического средства обработки данных трехфакторного планированного эксперимента -изготовление программного продукта, позволяющего моментально рассчитать план и обработать выходные данные эксперимента. Разработка программы осуществлялась в среде визуального объектно-ориентированного программирования Microsoft Visual Basic. В качестве физического объекта исследования выбран мелкозернистый бетон, наполненный тонкомолотым известняком. Применялись сырьевые материалы: песок природный (мелкий, средний и повышенной крупности), портландцемент ЦЕМ I 42,5 Н (Лц28=53,9 МПа), молотый известняк с удельной поверхностью 450 м2/кг, гиперпластификатор Melflux 1641 F.
Алгоритм программы включает основные процедуры: расчета коэффициентов функции отклика, статистической обработки и процедуру визуализации математической модели. Все основные вычисления производятся циклично, что позволяет моментально перестраивать математическую мо-
дель, изменяя входные данные. Кроме того, алгоритм включает вспомогательную процедуру, обеспечивающую проверку синтаксической правильности вводимых данных. При допущении ошибок ввода данных программа корректирует действия пользователя по средствам текстового оповещения.
Интерфейс программного продукта реализован в виде логических блоков, позволяющих вводить исходные данные и изменять параметры вывода математической модели в интерактивном режиме (http://www.swsys.ru/uploaded/image/2014-4-dop/10.jpg).
В первом логическом блоке устанавливаются входные факторы эксперимента. В эксперименте варьировались количество вяжущей части, содержание молотого известняка в вяжущей части и количество добавки гиперпластификатора по массе вяжущей части. Блок ввода значений входных факторов представлен в виде таблицы (табл. 1), значения факторов задавались в натуральном ви-
де.
В расчете факторного плана значения уровней входных факторов принимаются в кодированном виде, при этом основной уровень (центр плана) каждого фактора обозначается как «0», а нижний и верхний уровни - «-1» и «+1» соответственно. Пересчет заданных пользователем натуральных значений факторов производится путем линейной интерполяции значений по формуле
х. = -
Xi - Х0.
AX
(1)
где x, - значение i-го фактора в кодированном виде; X, - значение ,-го фактора в натуральном виде; ДХ, - интервал варьирования ,-го фактора.
В эксперименте контролировалась величина предела прочности бетона на сжатие в возрасте 28 суток нормального твердения Ясж, МПа. Для определения воспроизводимости измерений выходного параметра фиксировалось по три параллельных измерения прочности бетона на сжатие. Согласно плану эксперимента было заформовано 10 опытных серий образцов по 3 образца в каждой серии. Рассчитанный трехфакторный план и поля ввода выходных параметров отображаются в виде таблицы (табл. 2) во втором логическом блоке.
После автоматической проверки введенных данных программа рассчитывает коэффициенты математической модели и выводит функцию отклика в третьем логическом блоке в следующем виде:
Гср = b0+b1*x1+b2*x2+b3*x3+M1*x1A2+b22*x2A2+ +b33*x3A2+b12*x1 *x2+b 13*x1*x3+b23*x2*x3 y = (47,128)+(14,77)*x1+(1,734)*x2+(5,624)*x3+ +(0,016)*x1A2+(0,471)*x2A2+(0,141)*x3A2+ +2,28*x1*x2+6,01*x1*x3+4,05*x2*x3
После получения математической модели проверяются значимость (отличие от нуля) коэффициентов модели и ее адекватность. Проверка коэффициентов на значимость производится с по-
255
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
Факторы, определяющие изменчивость оптимизируемых параметров Factors defying optimized parameter variability
Таблица 1 Table 1
Матрица планирования Уровень Интервал варьирования Наименование фактора
нижний (-1) основной (0) верхний (+1)
xl: 15 20 25 5 Вяжущая часть, %
x2: 0 5 10 5 Известняк, %
x3: 0 0,75 1,5 0,75 Гиперпластификатор,%
План трехфакторного эксперимента типа B-D13 A plan for a “-D13” three-factor experiment
Таблица 2 Table 2
№ опыта (u) Матрица планирования Натуральные значения переменных Выходной параметр
x1 x2 x3 Вяжущая часть, % Известняк, % Гиперпластификатор, % Y(u, 1) Y(u, 2) Y(u, 3)
1 -1 -1 -1 15 0 0 38,7 37,2 38,1
2 +1 -1 -1 25 0 0 51,5 50,1 51,3
3 -1 +1 -1 15 10 0 28,6 28,7 28,9
4 -1 -1 +1 15 0 1,5 28,8 29,3 29,1
5 -1 0,19 0,19 15 5,95 0,8925 32,7 31,6 32,8
6 0,19 -1 0,19 20,95 0 0,8925 48,5 49,3 48,7
7 0,19 0,19 -1 20,95 5,95 0 43,6 42,9 42,56
8 -0,29 +1 +1 18,55 10 1,5 52,9 52,4 53,1
9 +1 -0,29 +1 25 3,55 1,5 71,2 71,5 72,6
10 +1 +1 -0,29 25 10 0,5325 60,5 60,3 65,8
мощью критерия Стьюдента (t-критерия), который рассчитывается по формуле
t =
1L
S {Ъ,} ’
(2)
где bi - i-й коэффициент математической модели; S{bj} - среднеквадратическое отклонение в определении коэффициентов.
Среднеквадратическое отклонение в определении коэффициентов функции отклика определяется по формуле
s {Ъ }=Jc~Si, (3)
где Ci - величины, приведенные для плана B-Di3 в таблице 3; SE2 - дисперсия воспроизводимости в параллельных опытах.
Таблица 3
Величины С, для плана B-D13
Table 3
The Ct values for plan B-D13
*0 *1 *2 *3 *11 *12 *13 *22 *23 *33
0,868 0,159 0,159 0,159 0,594 0,226 0,226 0,594 0,226 0,594
Дисперсия воспроизводимости в параллельных опытах рассчитывается по формуле
1 N m , --ч2
-ZZ(>V-Уп) > (4)
St = -
N (m -1)
и=1 j=1
где N - количество опытов в плане; m - количество параллельных измерений в каждом опыте; yuj- -значение выходного параметра в u-м опыте, j-м параллельном замере; у - среднее значение выходного параметра в u-м опыте.
Расчетное значение t-критерия сравнивается с табличным ^абл для выбранного уровня значимости (как правило, 5 %) и данного числа степеней свободы N(m-l). При ti < ^абл коэффициент bi считается значимым.
Проверка адекватности математической модели производится по критерию Фишера (F-критерий). Для этого вычисляется дисперсия адекватности по формуле
S2 = ~Т7~ £(У - Уп )г, (5)
где n - количество значимых коэффициентов; уи - значение отклика, предсказанное по уравнению математической модели.
В свою очередь критерий Фишера рассчитывается как отношение
F = S2JS2 . (6)
Расчетное значение F-критерия сравнивается с табличным ^табл для выбранного уровня значимости (как правило, 5 %) и чисел степеней свободы
256
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
N(m-1) и (N-n3). При F<Fts5r уравнение математической модели считается адекватным. Результаты статистической обработки модели отображаются в четвертом логическом блоке (рис. 1).
по квадратичным функциям x2 = f (x3), x1 = f (x3) и x1 = f (x2) соответственно.
В шестом логическом блоке программы представлена интерактивная диаграмма изолиний, по-
Статистическая обработка данных планирования эксперимента:
Проверка значимости коэффициентов модели (при наличии параллельных измерений):
Дисперсия воспроизводимости в параллельных опытах...............:
Число степеней свободы..........................................:
Табличное значение критерия Стьюдента (;_та6л)..................:
1,256
20
■2U9"
Коэфф-т: b0 b1 b2 b3 b11 b12 b13 b22 b23 b33
К. Стьюд. (t) 45,136 33,051 3,88 12,585 0,019 4,272 11,279 0,545 7,6 0,163
Значимость: 1 1 1 1 0 1 1 0 1 0
Примечание: 1 - коэффициент значим, 0 - коэффициент незначим.
Проверка адекватности математической модели:
Дисперсия адекватности математической модели...
Число степеней свободы при значимых коэффициентах
Табличное значение критерия Фишера (Б_табл)....
Расчетное значение критерия Фишера (F).........
3,858
3
3,1
3,07
Вывод: по критерию Фишера уравнение математической модели является адекватным.
Модель применима для решения производственных задач.
Рис. 1. Блок статистической обработки математической модели Fig. 1. Mathematical model aggregation block
Математическая модель прочности мелкозернистого бетона признана адекватной по критерию Фишера (F=3,07 < Fj.a=3,1) и применима для решения рецептурно-технологических задач. Уравнение математической модели представляет собой квадратичную функцию трех переменных: у = 47,13 +14,77x + 1,73x2 + 5,62x3 + 0,02X +
+ 0,47 x2 + 0,14x32 + 2, 28XjX2 + 6,01x3x3 + 4,05x2x3. Графическая интерпретация функции трех переменных требует четырехмерного пространства, поэтому для визуального упрощения и удобства работы с математической моделью функцию трех переменных необходимо преобразовать в функцию двух переменных, поочередно принимая константой один из факторов. В пятом логическом блоке программы представлены средства для преобразования уравнения регрессии в функцию двух переменных. Пользователь может установить постоянный фактор и задать его значение (в пределах интервала варьирования) в кодированном и натуральном видах.
В результате преобразования получаются три варианта математической модели: y = f (x2, x3) при x1 = const, y = f (x1, x3) при x2=const и y = f (x1, x2) при x3=const. Для визуализации каждого из трех видов уравнений строится диаграмма линий равного уровня (изолиний), представляющая собой проекции трехмерных поверхностей на плоскости (x2 © x3), (x1 © x3) и (x1 © x2). Таким образом, кривая каждой изолинии строится в координатах (x2, x3), (x1, x3) и (x1, x2), а ее построение производится
зволяющая пользователю снимать координаты факторного поля и значения выходного параметра в режиме реального времени (рис. 2).
Обработка данных планированного эксперимента завершается процедурой обнаружения экстремума функции отклика. Для определения координат точки экстремума производится автоматическое вычисление первой производной по каждому из значений факторов. Корни полученной системы уравнений представляют собой координаты точки экстремума исследуемого уравнения регрессии:
dy ( x,, x2) dxl
<
dy ( xt, x2)
dx2
= 0, = 0.
(7)
а) б) в)
Рис. 2. Диаграмма изолиний математической модели прочности бетона: x1=const (а), x2=const (б), x3=const (в)
Fig. 2. A contour chart for a concrete strength mathematical model: x1=const (а), x2=const (б), x3=const (в)
257
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
46 -i
44 4
42 1
< -fc. о
38 \
36 -
34
32 1 CO o' 'Я 1 о
0,8 " 0,5 " 0,2 -0,1 -0,4
X,
-0,7 3
2
а)
X2
б)
Рис. 3. Поверхность отклика (а) при x1=const и ее сечение (б) при x1=const и x3=const
Fig. 3. The response surface (a) when x1=const and its cross section (б) when x1=const and x2=const
Экстремум функции отклика математической модели прочности мелкозернистого наполненного бетона обнаружен в пределах варьирования переменных факторов. Значение экстремума составляет 32,7 МПа. Экстремуму функции отклика соответствуют значения факторов: x2 = 5,3 % (процентное содержание известняка в вяжущей части) и x3 = 0,8 % (процентное содержание гиперпластификатора) при xj = 15 % (общее количество вяжущей части в смеси).
Программа оснащена дополнительными функциями загрузки/сохранения данных, а также
функцией экспорта результатов расчетов в математический пакет Microsoft Excel, где пользователь может осуществить дополнительные построения, например, построить поверхность отклика в трехмерной системе координат (рис. 3а), а также сечения поверхности отклика, позволяющие проанализировать изменения выходного параметра в зависимости от одного переменного фактора (рис. 3б).
Разработанное программное средство можно применять в любых научно-прикладных задачах по оптимизации свойств объекта исследования, подбора рецептуры и технологических параметров, где используется математическое моделирование методом ортогонального планирования экспериментов.
Литература
1. Белов М.А., Антипов О.Е. Принципы проектирования виртуальной компьютерной лаборатории на основе технологии облачных вычислений // Современные проблемы и пути их решения в науке, транспорте, производстве и образовании: сб. докл. Междунар. конф. Одесса: Изд-во УкрНИИМФ, 2010. С. 27-30.
2. Лесовик В.С. Геоника (геомиметика) как трансдисциплинарное направление исследований // Высшее образование в России. 2014. № 3. С. 77-83.
3. Соловов А.В. Виртуальные учебные лаборатории в инженерном образовании // Индустрия образования: сб. статей. Вып. 2. М.: Изд-во МГИУ, 2002. С. 386-392.
4. Норенков И.П., Зимин А.М. Информационные технологии в образовании. М.: Изд-во МГТУ им. Н.Э. Баумана, 2004. 352 с.
5. Белов В.В., Образцов И.В. Виртуализация физических процессов в теории и практике строительного образования // Теория и практика повышения эффективности строительных материалов: матер. V Всерос. конф. студ., аспирантов и молодых ученых. Пенза: Изд-во ПГУАС, 2010. С. 186-189.
6. Афанасьев В.О., Бровкин А.Г. Исследования и разработка системы интерактивного наблюдения индуцированной виртуальной среды (системы виртуального присутствия) // Космонавтика и ракетостроение. 2001. N° 20. С. 19-2 J.
7. Колганов Д.А. Нереальная физика. Тестирование NVIDIA PhysX на конфигурации SLI Multi-Card // Игромания. 2010. Февраль. С. 162-164.
8. Zhang G., Torquato S. Precise algorithm to generate random sequential addition of hard hyperspheres at saturation. Physical review, E 88, 20j3, pp. 0533j2-j-9.
9. Баженов Ю.М., Демьянова В.С., Калашников В.И. Модифицированные высококачественные бетоны: науч. издание. М.: Изд-во Ассоциации строительных вузов, 2006. 368 с.
10. Белов В.В., Образцов И.В., Реунов А.Г., Смирнов М.А. Компьютерный метод расчета зернового состава заполнителя строительного композита // Актуальные проблемы строительства и строительной индустрии: сб. матер. Х Междунар. науч.-технич. конф. Тула, 2009. С. 6-7.
DOI: 10.15827/0236-235X.108.254-259
Received 29.08.2014
DEVELOPMENT OF ALGORITHMIC DATA PROCESSING TOOL OF PLANNED THREE-FACTOR EXPERIMENT TO CALCULATE A STATISTICALLY ADEQUATE MATHEMATICAL MODEL OF CONCRETE STRENGTH Belov V.V., Dr.Sc. (Engineering), Professor, Head of Chair, [email protected];
258
Программные продукты и системы /Software & Systems
№ 4 (108), 2014
Obraztsov I.V., Engineer, [email protected];
Kuryatnikov Yu.Yu., Ph.D. (Engineering), Associate Professor, [email protected] (Tver State Technical University, Quay Nikitin 22, Tver, 170026, Russian Federation)
Abstract. Functional modeling is widely applied in technology of building composite materials alongside with physical (manufacturing of the material physical sample) and structurally-imitating (imitation of system structural elements interoperability) modeling. The result is reception of the certain mathematical function that describes research object behavior abstracting from internal structure of a material substratum.
The authors develop the program of B-D13 three-factor planned experiment data processing. It allows counting the plan and processing experiment results. The program algorithm includes a procedure of calculation of response function factors, statistical processing and mathematical model visualization. All basic calculations are made cyclically that allows reconstructing a mathematical model abruptly, changing input data. Besides, the algorithm includes the auxiliary procedure providing input data syntactic check. The program corrects user actions using text notification when there are input data errors. The program has additional functions of data loading and conservation, as well as the function of exporting calculations results to Microsoft Excel where the user can make additional developments. The software product interface is implemented in the form of logical blocks that allow entering initial data and changing parameters of mathematical model output in interactive mode.
The developed software product is applied to a planned experiment on research of concrete mixture parameters influence on properties of the filled fine-grained concrete. The authors received a statistically adequate mathematical model of concrete strength. The developed software can be used in any scientifically-applied problems to optimize research object properties, select compounding and technological parameters where mathematical modeling is used by an experiments orthogonal planning method.
Keywords: functional model, planned experiment, three-factor plan, research object, response function.
References
1. Belov M.A., Antipov O.E. Design principles for virtual computer laboratory on the basis of cloud computing technology. Sovremennye problemy i puti ikh resheniya v nauke, transporte, proizvodstve i obrazovanii: sb. dokl. Mezhdunar. konf. [Proc. of the Int. Conf. Modern Problems and their Solving in Science, Transport, Production and Education]. Odessa, UkrNUMF, 2010, pp. 27-30 (in Russ.).
2. Lesovik V.S. Geonica (geodaetica) as a transdisciplinary field of research. Vyshee obrazovanie vRossii [Higher Education in Russia]. 2014, no. 3, pp. 77-83 (in Russ.).
3. Solovov A.V. Virtual training laboratories in engineering education. Industriya obrazovaniya: sb. statey [Education Industry: Proc.]. Iss. 2. Moscow, Moscow St. Industrial Univ. Publ., 2002, pp. 386-392 (in Russ.).
4. Norenkov I.P., Zimin A.M. Informatsionnye tekhnologii v obrazovanii [Information Technologies in Education]. Moscow, MSTU Publ., 2004, 352 p.
5. Belov V.V., Obraztsov I.V. Virtualization of physical processes in the theory and practice of construction education. Teoriya ipraktikapovysheniya effektivnosti stroitelnykh materialov: materialy V Vseros. konf. stud., aspirantov i molodykh uchenykh [Proc. 5th All-Russian Conf. The Theory and Practice of Increasing Efficiency of Constructional Materials]. Penza, Penza St. Univ. Architecture and Construction Publ., 2010, pp. 186-189 (in Russ.).
6. Afanasyev V.O., Brovkin A.G. Research and development of a system for interactive monitoring induced virtual environment (virtual presence). Kosmonavtika i raketostroenie [Astronautics and Rocket Production]. 2001, no. 20, pp. 19-21 (in Russ.).
7. Kolganov D.A. Unrealistic physics. Testing NVIDIA PhysX in SLI configuration Multi-Card. Igromaniya. 2010, February, pp. 162-164 (in Russ.).
8. Zhang G., Torquato S. Precise algorithm to generate random sequential addition of hard hyperspheres at saturation. Physical Review. E 88. 2013, pp. 053312-1-9.
9. Bazhenov Yu.M., Demyanova V.S., Kalashnikov V.I. Modifitsirovannye vysokokachestvennye betony: nauchnoe izdanie [Modified High-Quality Concretes: Scientific Publication]. Moscow, Assotsiatsiya Stroitelnykh Vuzov Publ., 2006, 368 p.
10. Belov V.V., Obraztsov I.V., Reunov A.G., Smirnov M.A. Computer calculation method of the grain structure of a filler composite construction. Aktualnye problemy stroitelstva i stroit. industrii: sb. materialov X Mezhdunar. nauch.-tekhnich. konf. [Proc. 10th Int. Science and Technical Conf. Important Building and Building Industry Problems]. Tula, 2009, pp. 6-7 (in Russ.).
259
