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

Размещено на http://www.allbest.ru/

Северо-Кавказский филиал Белгородского Государственного Университета им. В.Г. Шухова

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

Любимов С.В., аспирант

1. Постановка задачи

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

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

2. Генетическое программирование

Генетическое программирование представляет один из видов эволюционных алгоритмов. С помощью методов генетического программирования решаются комбинаторные задачи, в которых необходимо подобрать не столько значение параметров, а скорее их расположение. Например, с помощью методов генетического программирования решаются задачи подбора расписания, подбор электрической схемы, подбор топологии нейронной сети [1]. В генетическом программировании, эволюционирует популяция программ. Каждая программа может быть представлена различными способами: деревом, стековым языком, языком списков (например, Lisp).

При представлении программы в виде дерева, вершинами являются функции, а их потомки-аргументами. Например, выражение будет иметь такое представление (Рис. 1):

Рис. 1. Дерево арифметического выражения

Оператор кроссовера. Оператор скрещивания - предназначен для получения новых решений, на основе тех, что находятся в данный момент в популяции. Получает на вход две или более хромосом, на выходе выдает комбинированное решение, которое построено на основе входных решений. Как и в генетических алгоритмах, операция кроссовера в генетическом программировании служит для получения потомка из двух и/или более предков. Самый распространенный вид кроссовера - subtreecrossover. На вход, получаем две программы-дерева и выбираем случайно (и независимо для каждого дерева) точку скрещивания. Затем создается потомок, в котором поддерево с корнем в точке скрещивания в качестве корня заменено на такое же поддерево из второго дерева (рис. 2). Оператор скрещивания должен контролировать размер получившегося потомка, что глубина дерева не превышала заранее установленную.

Оператор мутации. Предназначен для вывода популяции из локального оптимума. Он применяется с некоторой вероятностью к каждому решению полученному в результате операции скрещивания. Самый популярный вид мутации - subtreemutation. Алгоритм работы аналогичен алгоритму кроссовер аsubtreecrossover. Случайно выбирается точка мутации в дереве и заменяется на случайно сгенерированное дерево. Чаще всего это реализуют, как операцию скрещивания мутирующего дерева, с деревом сгенерированным оператором инициализации.

Рис. 2. Оператор кроссовера

Рис. 3. Оператор мутации

Операция вычисления. Данная операция применяется к решению для получения закодированного в решении значения. Например, для программы, представленной деревом на рисунке 4 является 7 (max(1+2,3+4)).

Рис. 4. Арифметическое выражение max (1+2,3+4)

Методика тестирование. В наличии имеются два журнала испытаний [4]. В первом содержится 18 испытаний, во втором 50. Каждое испытание содержит следующие параметры:

Процентное содержание C₃S, C₂S,C₃A, C₄AF.

Прочность на 7-е сутки

Прочность на 28 сутки

Задачей алгоритма является подбор функции, аппроксимирующей прочность на 28 сутки

Параметры алгоритма:

Размер популяции - 300

Вероятность мутации - 90

Максимальная глубина дерева - 4

Процент данных, участвующих в отборе - 75%

Максимальное количество поколений - 200

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

Аппроксимация с учетом всех параметров. Результаты представлены на рис. 5-10.

Рис. 5.Результат аппроксимации (среднее значение ошибки 17%)

Рис. 6. Результат аппроксимации (среднее значение ошибки 6%)

Аппроксимация с учетом всех параметров и применением bootstrap. Для увеличения количества данных можно применить такой метод, как bootstrap. Метод формирования нескольких выборок данных того же размера, что и исходная генеральная совокупность, но с разными распределениями интересующей величины. Бутстрэп-метод был предложен в 1977 году Б. Эфроном. Он представляет разновидность рандомизированной обработки данных.

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

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

Рис. 7. Результат аппроксимации (среднее значение ошибки составляет 15%)

Рис. 8. Результат аппроксимации (среднее значение ошибки составляет 5%)

Рис. 9. Результат аппроксимации (среднее значение ошибки составляет 22%)

Рис. 10. Результат аппроксимации (среднее значение ошибки составляет 15%)

Из диаграмм видно, что существенное влияние на точность аппроксимации оказывает информация о прочности на 7-е сутки, что подтверждается и эмпирическими формулами [4]:

,

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

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

Литература

строительный материал кроссовер программирование

1. Poli R., William B. Langdon, Nicholas F. McPhee J., Koza R. A Field Guide to Genetic Programming.

2. Сайт компании Base Group Labs: http://www.basegroup.ru/glossary/definitions/bootstrap/.

3. Яковис Л.М. Многокомпонентные смеси для строительства. Расчетные методы оптимизации состава - Л.: Стрройиздат, 1988.

4. Акиева Е.А. Прогнозирование марочной прочности цементных систем по результатам краткосрочных испытаний и минералогическому составу: Дисс. канд. техн. наук. спец. 05.23.05 / Е.А. Акиева. Научн. рук. Ш.М. Рахимбаев. Белгород: БГТУ им. Шухова, 2006.

5. Борисов В.В., Круглов В.В., Федулов А.С. Нечеткие модели и сети. - М.: Горячая линия - Телеком, 2007.

Размещено на Allbest.ru