Введение В настоящее время основной причиной аварий на строительных объектах, возникающих в процессе эксплуатации, является значительное снижение показателей их надежности [1]. Данный факт обусловлен увеличением смещений в узловых точках каркаса как результатом разболтанности конструктивных элементов и появлением различного рода дефектов при длительной работе кранового технологического оборудования. Для определения периода безопасной эксплуатации здания в современной инженерной практике большое распространение получили детерминированные методы оценки его надежности [2]. Существующие методики нельзя отнести к тем, что отличает высокая степень достоверности получаемых результатов, поскольку они не учитывают значительное количество условий. Вероятностные методики не получили большого распространения ввиду большой сложности вычислительного характера, особенно в части расчета пространственных моделей зданий и сооружений. Таким образом, актуальной задачей становится разработка автоматизированной системы расчета и прогнозирования остаточного ресурса строительных объектов, реализующей алгоритм объединения комплекса методик последовательного решения прямой, обратной и прогнозной задач вероятностной оценки изменений деформации каркаса строительного объекта в процессе его эксплуатации с учетом вариативности воздействий. Анализ предметной области, разработка алгоритма поиска остаточного срока службы промышленного здания Развитие промышленности послужило толчком к изучению принципов проектирования различных одноэтажных каркасных сооружений промышленного назначения с использованием грузоподъемных механизмов. Однако довольно быстро стало понятно, что работу конструкций каркасных одноэтажных сооружений промышленного типа необходимо рассматривать в непрерывной взаимосвязи друг с другом. Таким образом, возникает потребность учета пространственной работы каркаса, своевременного контроля запаса прочности и роста обобщенной величины нагрузки вследствие технологического перевооружения производственного процесса. Основным фактором, который влияет на изменение прочностных характеристик промышленного здания, является степень восприятия его отдельными конструктивными элементами обобщенной нагрузки, включающей воздействие мостовых кранов, сейсмическую и ветровую составляющие кратковременной нагрузки, собственный вес конструкций. Обеспечение надежности прочностных характеристик каркаса строительного объекта, оборудованного мостовыми кранами, является одной из актуальных проблем при их проектировании, возведении и эксплуатации. Следует выделить ряд составных частей, которые являются неотъемлемыми элементами исследования расчетной схемы строительного объекта промышленного типа: геометрические, жесткостные и нагрузочные. Стоит отметить, что работа конструкций каркаса напрямую зависит от диапазона смещений в соответствующих расчетных точках [3]. В настоящее время используется несколько методик оценки надежности промышленных зданий, реализация которых зависит от выбора модели представления и обработки данных, полученных при обследовании объектов. Результаты анализа достоинств и недостатков этих моделей, а также возможности их численной реализации для условий функционирования конкретно взятого объекта представлены в табл. 1. Таблица 1 Модели оценки надежности промышленных зданий Класс модели Модель учитывает Достоинства Недостатки Вероятностная модель Изменение надежности строительных конструкций. Вероятностное распознавание категорий технического состояния строительных конструкций. Изменение вероятности отказа промышленного объекта. Обоснованность. Детальная проработка. Построение прогноза динамики изменения показателей надежности на длительной период времени эксплуатации объекта. Отсутствие критериев по дифференциации объектов. Изменение поведения конструктивных элементов каркаса сводится к сложной системе нелинейных дифференциальных уравнений. Сложность вычисления. Детерминированная модель Изменение параметров технического состояния промышленного здания. Ухудшение характеристики несущей способности. Утвержденные нормативы сроков эксплуатации промышленного здания до капитального ремонта. Степень физического износа промышленного объекта. Обоснованность. Детальная проработка. Возможность численного решения прямой задачи строительной механики. Получение результата на момент проведения технического обследования эксплуатируемого объекта. Отсутствие критериев по дифференциации объектов. Невозможность учета динамики изменений физических параметров материалов конструкций. Ограниченность объема имеющейся статистической информации об изменении свойств объекта исследования. Невозможность использования для определения о статочного ресурса промышленного здания. Анализ результатов изучения предметной области позволил сформировать схему, отображающую последовательность действий исследователя и направленную на получение конечного результата в виде сроков достижения объектом предельных состояний в работе его конструктивных элементов. Данный алгоритм сводится к оценке величины резерва прочности конструкций каркаса здания, который определяется из соотношения несущей способности и максимума обобщенной нагрузки при заданных условиях функционирования (рис. 1). Рис. 1. Схема проведения исследования напряженно-деформированного состояния промышленного здания: МКЭ - метод конечных элементов Для оценки остаточного ресурса промышленного здания организуется поиск изменений амплитудных значений смещений в расчетных точках построенной физической модели. Данные изменения определяются по итогам сравнительного анализа натурных измерений соответствующих параметров в разные моменты времени эксплуатации данного объекта при неизменной величине обобщенной нагрузки. Однако для определения изменений напряженно-дефор-мированного состояния каркаса здания намеренно, заведомо ложно, предполагаем, что зафиксированное увеличение смещений в расчетных точках вызвано возрастанием нагрузки. А потому функцию нагрузки, построенную в виде регрессионной зависимости по фактору времени, считаем фиктивной. Зафиксировав новое состояние объекта на момент проведения обследования, возвращаясь к исходному загружению расчетной схемы, корректируем предыдущую матрицу жесткости каркаса. Организовав поиск расчетных значений математических ожиданий резерва прочности, определяем среднегодовую скорость износа. Согласно представленному комплексному подходу данная система показателей и является основанием для реализации алгоритма по определению остаточного срока службы промышленного здания. Реализация данного алгоритма в вероятностной постановке с получением конкретного результата возможна только при задействовании аппаратных средств электронно-вычис-лительных машин. Известно несколько систем автоматизированного проектирования работ, позволяющих решить поставленную задачу. В табл. 2 проведен сравнительный анализ их функциональных возможностей. Таблица 2 Автоматизированные системы проектирования Программный продукт Преимущества Недостатки SCAD Высокая скорость расчета. Автоматическая генерация произвольной сетки конечных элементов на плоскости. Возможность формирования сложных расчетных моделей путем сборки из различных схем. Развитые графические средства формирования и корректировки геометрии расчетных схем, описания физико-механических свойств материалов, задания условий опирания и примыкания, а также нагрузок. Эффективные методы оптимизации матрицы жесткости. Высокая стоимость. Невозможность определения остаточного ресурса здания. Расчет ориентирован на жилое строительство. Отсутствие явных инструментов для отображения на расчетной схеме динамических крановых нагрузок. Невозможность корректировки матрицы жесткости по истечении фиксированного времени эксплуатации здания. ANSYS Широкая специфика применения, возможность проектирования самолетов, космических кораблей и других технологических продуктов. Моделирование и анализ в некоторых областях промышленности, работа на основе геометрического ядра Parasolid. Высокая стоимость. Невозможность определения остаточного ресурса здания. ЛИРА В расчете применяется метод конечных элементов. Возможность выполнения расчетов на статические (силовые и деформационные) и динамические воздействия. Возможность проведения подбора или проверки сечений стальных конструкций и(или) армирование сечений железобетонных конструкций. Высокая стоимость. Невозможность определения остаточного ресурса здания. Расчет ориентирован на жилое строительство. На основе анализа информации, изложенной в табл. 2, можно сделать вывод о необходимости разработки программно-расчетного комплекса, специализирующегося на проведении расчетов работы конструкций промышленного здания, оборудованного мостовыми кранами, при целом комплексе возмущающих воздействий и расширяющего спектр решаемых вопросов. В ходе исследования комплексной схемы, представленной на рис. 1, был разработан сводный алгоритм расчета промышленного здания, который в дальнейшем может быть реализован при разработке автоматизированной системы проектирования промышленных зданий. Данный алгоритм изображен на рис. 2. Рис. 2. Алгоритм оценки изменения напряженно-деформированного состояния каркаса промышленного здания при действии комплекса нагрузок: P - количество воздействий; V - уровень модели в зависимости от изменения матрицы жесткости; L - вид кранового воздействия на строительные конструкции сооружения; n - количество повторений при указании параметра L; h - бинарная переменная вариантов расчета на крановые нагрузки (0 - нормативный; 1 - вероятностный); k - параметр цикла Алгоритм включает поиск значений изгибающих моментов и напряжений, возникающих под действием конкретного загружения конструкций каркаса. На представленной блок-схеме продемонстрирован выбор варианта комплекса нагрузок, действующих на объект в фиксированный момент времени. Предложенный алгоритм реализует вероятностный подход, который наиболее соответствует случайному характеру техногенных [4] и природных воздействий [5], а также условиям их восприятия [6]. Техногенные воздействия на каркас промышленного здания определяются нагрузками, вызванными работой грузоподъемного оборудования. Алгоритм учитывает расчет трех составляющих крановой нагрузки: - L = 1 - вертикальное давление крана с грузом; - L = 2 - торможение крановой тележки; - L = 3 - боковые силы, возникающие при биении моста крана, перпендикулярно рельсовому пути. Согласно предложенному алгоритму значения крановых нагрузок генерируются в виде непрерывных случайных величин в соответствии с нормальным законом распределения при заданных параметрах математического ожидания и среднеквадратического отклонения. Для повышения точности итогового результата определение величины нагрузки и расчет последствий от их восприятия каркасом здания повторяется многократно. Алгоритм также позволяет выполнить сравнительный анализ полученных вероятностных значений с результатами расчета по нормативным значениям нагрузок при неизменных условиях функционирования системы. Основным преимуществом предложенной схемы проведения расчета является возможность решения последовательно сразу трех задач: прямой, обратной и прогнозной. Прямая сводится к поиску смещений в узловых точках расчетной схемы здания при известных геометрических и жесткостных его характеристиках, а также случайной величины обобщенной нагрузки. Решение обратной задачи позволяет скорректировать матрицу жесткости каркаса по истечении некоторого времени эксплуатации здания на основе зафиксированного увеличения смещений в контролируемых точках конструкции при очередном обследовании. Коррелируя значения частот собственных колебаний технической системы, полученных по итогам проведения нескольких обследований при неизменной величине обобщенной нагрузки, можно получить численный результат прогнозной задачи, что и является основанием для определения величины остаточного ресурса и других показателей надежности здания. Представленный алгоритм положен в основу разработки программного комплекса «АСР-2018», выступающего в качестве инструмента реализации всего спектра методик расчета. Данный программно-расчетный продукт предоставляет удобный интерфейс, что расширяет круг потенциальных пользователей системы. Заключение Главное достоинство разработанной схемы проведения расчета - возможность решения последовательно сразу прямой, обратной и прогнозной задач для поиска остаточного срока службы промышленного здания. Прямая задача представляет собой поиск смещений в узловых точках расчетной схемы промышленного объекта при известных геометрических и жесткостных характеристиках, а также учитывает случайную величину обобщенной нагрузки.