Введение. На данный момент проектирование стальных рамно-консольных стационарных конструкций покрытий над трибунами стадионов базируется на методе предельных состояний. Данный метод лежит в основе таких нормативных документов в области строительства как ДБН (Украина), СП (Россия), EUROCODE (Евросоюз), а также многих других. При использовании методики предложенной в этих документах возникает ряд спорных вопросов в области определения надежности конструкции. Метод предельных состояний требует сопоставления расчетного напряжения в элементе (Qp) с расчетным значением несущей способности (Ry). При выполнении условия Qp< Ry, считается, что надежность обеспечена. Таким образом, надежность конструкций во всех расчетных случаях оказывается приблизительно одинаковой и зависит только от конструктивной схемы и материала конструкции. Пытаясь учесть этот недостаток норм, ДБН [14)]и EUROCODE [18] ранжируют конструкции по классу последствий аварий (экономических и людских потерь) и категориям ответственности, а также приводят максимально допустимые вероятности отказа для этой классификации. Причем методика вычисления вероятности отказа приведена крайне обобщенной и не имеет четкой структуры или алгоритма действий. Поэтому значение вероятности отказа запроектированной конструкции, как правило, остается неизвестным для инженера. Учитывая все выше сказанное можно сказать, что современные нормы проектирования для уникальных и конструкций повышенной ответственности нуждаются в доработке, для таких конструкций предлагается выполнять численный расчет надежности и живучести на стадии проектирования, с учетом таких случайных величин и факторов как расчетное сопротивление материала, монтажные несовершенства и дефекты, величины временных нагрузок (в частности снеговой), просадки основания.Выполняя анализ литературы по этому вопросу можно сказать, что ряд авторов рассматривает проблему обеспечения надежности строительных систем в основном как конструктивную (требуемый уровень обеспечивается использованием более прочных материалов либо наиболее надежных конструктивных схем) [22]. Данный подход не всегда может быть применим с точки зрения экономических требований, выдвигаемых к конструкции, к тому же открытым остается вопрос численного определения показателей надежности. Он рассмотрен в ряде работ [1, 2, 3, 19, 20, 21]. Данные авторы предлагают различные методики вычисления показателей надежности (в некоторых случаях даже с приведением примеров расчета), но для сложных много раз статически неопределимых систем они тяжело реализуемы, а в некоторых случаях нереализуемы вообще. Так же данная тема обширно освещена в работах А.В. Перельмутера [6] и С.Ф. Пичугина [7] под руководством, которых была проведена большая исследовательская работа в этом направлении, использованная в дальнейшем в ДБН [13, 14]. К сожалению, авторы не привели единого метода вычисления показателей надежности. Выше перечисленные работы базируются на классических трудах (в этой области) Г. Шпете, В.Д. Райзера, А. Р. Ржаницына, В. В. Болотина [8, 9, 10, 11]. В них детально разобрана методика вычислений вероятности отказа и характеристики безопасности (как показателей надежности системы), но также на примере простых конструкций, что вызывает ряд спорных моментов при расчете данных показателей для сложных многоэлементных систем.В данной статье предложена методика подбора поперечных сечений стержневых элементов стальных рамно-консольных стационарных покрытий над трибунами стадионов с учетом нижней границы надежности наиболее ответственных стержней.Методы исследования. При расчете конструкций по методу предельных состояний задача обеспечения заданного уровня надежности ложится на ряд коэффициентов надежности. Для уникальных конструкций (какими являются стационарные покрытия над трибунами стадионов) их значения, как правило, указаны некорректно, либо не рассматриваются вообще. Таким образом, требования надежности для таких конструкций могут быть не выполненными в полной мере.Как известно из работ Ржаницына А.Р., Шпете Г., Райзера В.Д, основное уравнение безотказности работы конструкции имеет вид:Ŷ(t)=Ř(t)-Ŝ(t)>0;                      (1)где Ř(t) – обобщенная несущая способность конструкции или элемента; Ŝ(t) – обобщенная нагрузка на конструкцию (напряжения в элементе); Ŷ(t) – характеристика резерва прочности или резерва несущей способности;Вероятность отказа может быть представлена в виде интеграла:                 (2)где  – функция распределения вероятности случайной величины ;  – плотность распределения вероятностей случайной величины.В качестве случайных величин (для задачи, приведенной в статье) рассматриваются расчетное сопротивление материала, геометрические характеристики плоских сечений элементов конструкции, величины временных нагрузок (в частности снеговой), просадки основания [4, 10].Так как система много раз статически неопределима, наиболее целесообразно вычислять показатели надежности для наиболее ответственных элементов. С этой целью производится итерационный геометрически нелинейный расчет конструкции. Итогом, которого, является список элементов, поочередно вышедших из строя. Критерием отказа предлагается рассматривать наступление первого предельного состояния в каждом стержне в частности, так как для рассматриваемых конструкций переход через второе предельное состояние не повлечет за собой таких масштабных экономических и человеческих последствий.Для выбранных в предыдущем этапе расчете элементов выполняется итерационный расчет с учетом указанных случайных величин. Таким образом, формируется выборка напряжений для рассматриваемой группы стержней (Ŝ(t)). Выборка должна иметь значительный объем. Исследуя эту проблему ряд авторов (Райзер, Шпетте, Пичугин) приводят числа в интервале 104-108.В качестве второй составляющая уравнения №1 выступает случайная величина характеристики сопротивления материала конструкции (предел текучести) Ř(t)). Которое сформирована на основе анализа статистических данных полученных на металлургических заводах.Полученные две обобщенных случайных величины обрабатываются методами математической статистики. Определяются их законы и плотности распределения. В случае невозможности или трудности выполнения данной задачи рекомендуется подобрать наиболее соответствующее распределение и аппроксимировать случайную величину с точность удовлетворяющей расчет. Зная виды распределения случайных величин, вычисляются их характеристик (математическое ожидание, стандартное отклонение).Оперируя плотностью, законами и характеристиками распределений вычисляются вероятности отказа группы выбранных элементов. Определенная вероятность отказа является нижней границей отказа [11]. То есть полученная для системы, у которой все элементы включены в работу последовательно и выход из строя одного стержня повлечет за собой выход из строя всей системы. Существует еще верхняя граница отказа, в ней предполагается последовательное включение элементов в работу, но вероятность отказа в этом случае значительно выше и ее мы рассматривать не будем.Насчитанные характеристики надежности сравниваются с нормативными, с учетом ответственности рассматриваемой конструкции. В случае удовлетворения требований надежности расчет заканчивается. Если надежность не обеспеченна, пересчитываются частные коэффициенты надежности и решается обратная задача, целью которой, становится определить требуемое значение математического ожидания напряжений в элементах конструкций, а от него перейти к требуемым площадям стержней. Конструкция с переопределенными размерами сечений пересчитывается заново по выше описанному алгоритму. Итерации ведутся до момента удовлетворения условий прочности рассматриваемой группой элементов.Подобранные по приведенной методике сечения стержней стальных рамно-консольных покрытий над трибунами стадионов удовлетворяют не только требованиям расчета по предельным состояниям, но и соответствуют заданному уровню надежностиОбъект исследования. Рассмотрим методику подбора сечений стальных рамно-консольных стационарных покрытий над трибунами стадионов на примере уже эксплуатируемого деформационного блока покрытия южной трибуны главного игрового поля ФК Олимпик в городе Донецке.Конструктивная схема стационарного покрытия представляет собой стержневую систему, выполненную из «Г»-образных плоских консольных рам (в качестве основных несущих конструкций) установленных с шагом 6 м., связанных между собой ферменными связями и прогоном на консоли. Все узлы принимаются жесткими; следовательно, конструкция является много раз статически не определимой системой. Крепление к фундаменту конструктивно выполнено таким образом, что в расчетной схеме целесообразно рассматривать его как жесткую заделку (рис. 1).        Рис. 1. Конструкция стационарного покрытия над трибунами стадиона ФК «Олимпик» Рис. 2. Конструктивная схема рамно-консольного покрытия  Рис. 3. Геометрическая схема плоской рама и ферменной связи с указанием номеров жесткостей элементов конструкции Таблица 1 Сечения конструкции покрытия над трибунами стадиона СК «Олимпик»№НазваниеРазмер1Двутавр][ №162Швеллер[ №163Гнуто-сварнаяпрямоугольная труба□ 100×60×34Гнуто-сварная квадратная труба□ 60×2,55Гнуто-сварнаяпрямоугольная труба(прогон на консоли)□ 200×60×4,5 Алгоритм исследования. Универсального программного продукта для решения поставленной задачи на данный момент авторами не найдено. Поэтому, был создан ряд программ, написанных на высокоуровневом интерпретируемом языке программирования MATLAB. Все ниже приведенные расчеты были реализованы в них. Проверку и отладку программного комплекса выполняли с использованием программного комплекса ЛИРА-САПР 2013.Весь созданный программный комплекс целесообразно разбить на 2 блока:Определение наиболее ответственных элементов в конструкции, которые будут определять ее надежность в целом.Формирование выборки напряжений группы наиболее ответственных элементов системы. Обработка выборок напряжений и расчётного сопротивления методами математической статистики. И вычисление критериев надежности.Рассмотрим каждый из них более детально на примере описанной конструкции.Блок 1Нелинейный расчет конструкции основывается методе конечных элементов в пространственной постановке. Геометрическая нелинейность работы конструкции учтена следующим образом:На каждой итерации обновляются координаты узлов с учетом их перемещения в предыдущей.Напряжения, возникшие в предыдущей итерации от смещения узлов конструкции, учитываются на этапе формирования матрицы жесткости элемента путем ввода функций Новожилова В.В.К рассмотренному покрытию прикладывают снеговую нагрузку (приведя ее к узловой) и производится итерационный расчет до момента достижения заданного количества элементов, вышедших из строя. На каждом цикле расчета дополнительно прикладывается нагрузка 1 кН. Если на предыдущем цикле какой-либо элемент вышел из строя, он исключается из системы и она пересчитывается заново без приложения дополнительной нагрузки. Таким образом, учитывается явление лавинообразного разрушения (рис. 4).У рассматриваемого покрытия на четвертой итерации вышли из строя элементы 32, 33, 60, 61,88, 89, 116,117 (рис. 5–8)   Рис. 4. Блок схема определения наиболее ответственных элементов конструкции  Рис. 5. Элементы конструкции, вышедшие из строя на 4 цикле расчета  Рис. 6. Элементы конструкции, вышедшие из строя на 5 цикле расчета  Рис. 7. График выхода из строя элементов конструкции на каждом шаге нагружения   Рис. 8. Положение стержней, определяющих надежность конструкции  Блок 2Моделирование снеговой нагрузки.Снеговая нагрузка как случайная величина принималась с двойным экспоненциальным законом распределения Гумбеля. Методом математической инверсии, было получено выражение насчитывающее выборку значений снеговой нагрузки для Донецкого региона.  (3)где R – генератор равномерной случайной величины в пределах [0…0,935].Необходимые параметры распределения Гумбеля получены путем анализа статистики снеговых максимумов за год в течении 40 лет (11).  Анализируя выражение «3» видно что должно выполнятся условие <0, чтобы значение снеговой нагрузки не перешли в область комплексных чисел. При использовании максимального значения R=0,935 этого не произойдет. Данное выражение насчитывает снеговую нагрузку в пределах от 200 до 1500 Па, что соответствует характеру нагрузки для рассматриваемого региона.Моделирование прочностных характеристик материала и геометрических характеристик сечений.Геометрические характеристики плоских сечений и прочностные характеристики   генерируются на основании анализа металлопроката двух металлургических заводов: «Луганский трубный завод» и «ПАО Днепропетровский металлургический завод им. Коминтерна» (5). Таблица 2Стохастические характеристики размеров сечений гнутосварных труб (ГОСТ 8645-68), полученные опытным путем№СечениеКол-возамеровТолщина стенки (мм)Высота сечения (мм)Ширина сечения (мм)Мат. ожид.Ст. откл.Мат. ожид.Ст. откл.Мат. ожид.Ст. откл.160×40×32502,960,2558,981,4838,772,152120×60×32502,980,19118,073,9858,392,793100×100×32502,930,2699,135,6098,345,854100×60×32502,890,2397,515,6357,791,71   Выборка прочностных свойств сечения и выборки геометрических замеров всех сечений были подвергнуты χ2- анализу, на предмет возможности аппроксимации распределения случайных величин «нормальным законом распределения», по критерию Пирсона. Анализ показал, что распределение каждой из случайных величин отличается от нормального не более чем на 5 % (был задан уровень значимости 0,05). Следовательно, можно считать данные случайные величины подчиняющимися нормальному закону распределения.Моделирование геометрических несовершенств и просадок основания.Просадки фундаментов и геометрическую неточности конструкции учли путем анализа геодезической сьемки южной трибуны главного игрового поля ФК «Олимпик». Сьемка геометрии конструкции проходила в два цикла и дала результаты, отображенные в гистограмме (рис. 11). Данные учтены на этапе формировании локальных матриц жесткости расчетной схемы конструкции. Таблица 3Значения прочностных характеристик сечения гнутосварных труб (ГОСТ 8645-68) полученные опытным путем№СечениеСтальМат. ожид.СтандартЧисл. станд. (μ)Эксп. знач. (мПа)Нормативное значение (мПа)              160×40×3Ст3пс482,58317,4422,3835,153,15206,72411,502253702120×60×3Ст3пс443,67332,2131,0722,073,15262,69345,802253703100×100×3Ст1пс428,85323,6328,5331,883,15223,21338,98--4100×60×3Ст3пс418,54379,0325,8429,463,15286,23337,14225370   Рис. 9. Схема геодезической сети для мониторинга конструкций трибун стадиона  Рис. 10. Схема контрольных точек на конструкциях трибун    Рис. 11. Гистограмма распределения отклонений координат узлов конструкции с аппроксимирующимикривыми  Рис.12. Блок-схема алгоритма определения критериев надежности стационарных рамно-консольных покрытий над трибунами стадионов Зная законы распределения всех рассматриваемых стохастических величин, используя генератор случайных чисел, формируем выборки в объёме 10000. Далее, используя детерминированный расчет на основе метода конечных элементов, насчитываются выборки напряжений в рассматриваемых восьми элементах.Полученные случайные величины напряжений и предела текучести анализировались методами мат статистики (ХИ2-анализ, критерий Шапиро-Уилка). Так ка полученные выборки напряжений в стержнях и предела текучести аппроксимируются нормальным законом распределения, в качестве вероятностной характеристики вычислялась характеристика безопасности β (дальность отказа). Блок схема расчета представлена на рисунке 12.Результаты исследования. Полученные значения β сопоставляем с нормативными (14). Требования норм не выполняются ( ). Вычисляем требуемое значение математического ожидания площади элементов. Принимаем новые сечения и полностью пересчитываем второй блок программы с учетом вычисленных коэффициентов условий работы. Результаты расчета сведены в таблицу 4. Анализируя полученные результаты можно заметить: 1) явно прослеживается симметрия конструкции; 2) насчитанные значения γс в значительной мере отличаются от рекомендованных нормами ДБН для конструкций повышенной ответственности (класс последствий “СС3”, категория ответственности “А”) Таблица 4Результаты подбора площади сечения элементов конструкции с учетом требований надежностиИтерация№1№2Нормативные значения№ элементаСеченияβрасчγсСеченияβрасчγсγс(ДБН) (ДБН)32□ 100×60×34,42550,9729□ 100×60×54,98270,97340,94,7633□ 100×60×34,64040,9742□ 100×60×55,15730,974860□ 100×60×34,07750,9708□ 100×60×54,76340,971261□ 100×60×34,29660,9722□ 100×60×54,81350,972588□ 100×60×34,07750,9708□ 100×60×54,76340,971289□ 100×60×34,29660,9722□ 100×60×54,81350,9725116□ 100×60×34,42550,9729□ 100×60×54,98270,9734117□ 100×60×34,64040,9742□ 100×60×55,15730,9749  Выводы.Методика расчета по предельным состояниям не обеспечивает требуемого уровня надежности уникальных конструкций, каковыми являются стационарные стальные рамно-консольные покрытия над трибунами стадионов. Следовательно, можно сделать вывод о необходимости применения прямых методов теории надежности для подбора сечений элементов конструкции, определения показателей надежности и последующего уточнения значений γс.Рекомендуемые нормами общие меры учета ответственности сооружения (учет коэффициента надежности по ответственности γn, класс последствий здания и сооружения, категории ответственности конструкций и их элементов) не всегда обеспечивают требуемый уровень надежности  для сооружений повышенной ответственности. Для отдельных ключевых элементов конструкции, определяющих вероятность ее отказа, наиболее просто это сделать корректировкой значений коэффициента условий работы γc, но в существующей нормативной литературе их значения предлагаются без такого учета степени ответственности сооружения и нуждаются в уточнении.Так как нормативные документы в области строительства не дают четкого алгоритма расчета показателей надежности уникальных (да и каких-либо других) строительных конструкций, для стержневых систем стационарных покрытий над трибунами возможно применения предлагаемой методики как в части прямого вероятностного расчета конструкций, так и в части нормирования коэффициентов условий работы наиболее ответственных элементов конструкций, обеспечивающих требуемый уровень проектной надежности.