CC BY
13УДК 624.012.36:311.214
ПОЛНЫЙ ВЕРОЯТНОСТНЫЙ РАСЧЕТ НАДЕЖНОСТИ КОНСТРУКТИВНОГО ЖЕЛЕЗОБЕТОННОГО ЭЛЕМЕНТА
канд. техн. наук, доц. Д.О. ГЛУХОВ; канд. техн. наук, доц. Р.П. БОГУШ; канд. техн. наук, доц. Е.Д. ЛАЗОВСКИЙ; Т.М. ГЛУХОВА (Полоцкий государственный университет)
Представлен комплексный подход оценки надежности строительных конструктивных элементов на основе нелинейных вероятностных моделей наиболее опасных участков и расчетных ситуаций. Сделан вывод о необходимости проверки гипотезы о том, что функция предельного состояния является гауссовой смесью вероятностных распределений, а в случае состоятельности данной гипотезы необходимо выполнение процедуры разделения смеси на компоненты. Предложена интерпретация данного факта как эффекта влияния нелинейных свойств и сингулярности моделей на распределение функции предельного состояния.
Ключевые слова: конструктивный железобетонный элемент, оценка надежности, нелинейные вероятностные модели.
При рассмотрении воздействия Q на конструктивный элемент как случайную величину, учитывая случайный характер функции предельного состояния Я, задачу оценки надежности конструктивного элемента рассматривают как задачу анализа характеристик случайной величины М = Я - 5".
Если в результате статистического моделирования в рамках полного вероятностного расчета получены оценки плотности вероятности для Q и Я, то вероятность разрушения элемента Ру определяют как область перекрытия распределений, что показано на рисунке 1.
Плотность
вероятности
Вероятность i
разрушения, Ffj
R, Q
Рисунок 1. - Интерпретация области перекрытия распределений Q и R
Вероятность разрушения в этом случае определяется формулой
Pf = P(P < Q) = Р(М < 0),
где Pf - вероятность разрушения элемента.
Определение точного значения данного интеграла является в большинстве случаев невыполнимой задачей, поэтому применяются упрощенные методы оценки значения интеграла, которые носят название FORM (First Order Reliability Method) и SORM (Second Order Reliability Method) [18].
Для оценки вероятности Pf широкое распространение получила процедура определения индекса надежности в (индекс безопасности по С. А. Корнеллу [5], характеристика безопасности по А.Р. Ржаницыну [6], индекс надежности по Раквицу и Фислеру [3; 4]).
Индекс надежности рассчитывается исходя из предположения о нормальном законе распределения функции воздействия, функции предельного состояния и независимости случайных величин Q и R:
р=Ф-1 (-P )=
4
S 2 + Sq
где mR и mQ - математические ожидания распределений Q и R, sR и Sq - соответствующие среднеквадра-тические отклонения.
R
Q
mR - mQ
Применению метода Монте-Карло для статистического моделирования, разработке методов уменьшения размерности пространства оценки, фокусировке процесса моделирования на пограничной подобласти и т.п. посвящено множество публикаций [4; 20-23].
В соответствии с ИСО СТБ 2394 критерий обеспечения надежности задается неравенством:
Р^,
где вщ - минимальный целевой индекс надежности, численные значения которого задаются в зависимости от класса ответственности конструкции.
Для оценки индекса надежности выполняется полный вероятностный расчет, в котором с помощью генераторов псевдослучайных последовательностей (ГПСП) моделируются прочностные характеристики материалов, геометрические параметры, площадь арматуры. Строится вероятностная модель воздействия на конструктивный элемент. Причем модель воздействия не влияет на оценку распределения функции предельного состояния, но оказывает существенное влияние на распределение оценок по деформациям.
Для реализации полного вероятностного расчета нелинейных моделей строительных элементов нами были построены нелинейные модели основных расчетных ситуаций. Общая схема полного вероятностного нелинейного расчета приведена на рисунке 2.
Оценка надежности
строительной
конструкции
С
Вероятностная модель функции предельного состояния
Проверка гипотезы о смеси ЕМ или GEM алгоритм разделения
Рисунок 2. - Общая схема полного вероятностного нелинейного расчета
Проектирование строительных конструкций опирается на полувероятностные методы расчета, основанные на частных коэффициентах безопасности. Калибровка коэффициентов безопасности выполняется с учетом особенностей производства и систем контроля качества, принятых в той либо иной стране [8]. Действующие нормы Республики Беларусь, определяя параметрический ряд классов бетонов и арматуры, определяют семейство вероятностных моделей материалов. Параметры распределения случайной величины восстанавливаются по нормативным и средним значениям прочностных характеристик. Однако, в сравнении с анализом результатов при обследовании, для моделирования тех или иных величин применяются распределения, отличные от нормального распределения.
Общие принципы обеспечения надежности строительных и иных конструкций изложены в базовом документе EN 1990 [8] и ISO 2394 [9]. Еврокоды предлагают готовую систему теоретических распределений для анализа и моделирования случайных свойств конструкций. В частности, в Eurocode 3 «Стальные конструкции» [10] при определении прочности (предела текучести) стали по результатам испытаний применяется логарифмически нормальный закон распределения для моделирования как геометрических, так и прочностных параметров.
Объединенный комитет по безопасности конструкций JCSS, выполняющий работы по калибровке частных коэффициентов безопасности, закладываемых в европейские нормы, предлагает вероятностную модель PMC (Probabilistic Model Code), где также определяет вид распределения для каждого учитываемого параметра:
- прочностные характеристики - логарифмически нормальный закон распределения;
- плотность бетона - нормальное распределение;
- длительно действующие нагрузки - гамма-распределение;
- кратковременные нагрузки - экспоненциальное распределение;
- сопротивление конструкции - нормальный закон распределения.
Отдельные параметры допускается рассматривать как подчиненные распределению Вейбула и Гумбеля [7]. Нами в работе [11] систематизированы методы идентификации параметров применяемых теоретических распределений по результатам детального обследования строительных конструкций, в работе [12] выполнен выбор и разработано программное обеспечение для реализации генераторов псевдослучайных последовательностей для статистического моделирования строительных конструкций.
Нелинейные модели строительных элементов. Системы нелинейных сингулярных уравнений являются основой целого ряда математических моделей. В частности, такие системы возникают при попытке описать напряженно-деформированное состояние железобетонных элементов, работающих под нагрузкой при внецентренном сжатии и изгибе, при работе такого элементы в предопорной зоне на срез с изгибом. Сингулярные свойства таких моделей связаны с тем, что те либо иные параметры, зависимости, описывающие свойства модели, имеют точки разрыва первого и второго рода. Такие точки разрыва фактически разделяют всю область определения входных параметров на подобласти непрерывного поведения. Такие точки разрыва позволяют описать скачкообразное изменение напряжения в бетоне в момент трещинообразования.
Опыт нескольких лет изучения подобного рода математических моделей с целью разработки эффективных алгоритмов решения сингулярных систем нелинейных уравнений показал, что наиболее верным подходом является детальное изучение всех особенностей многомерной поверхности решения и применение различных стратегий поиска решения на поверхности решений в зависимости от области непрерывного поведения системы, в которой находятся значения входных параметров [13-17].
Из многолетнего опыта обследования технического состояния строительных конструкций, накопленного в Полоцком государственном университете, наибольший интерес представляют 3 нелинейные модели:
1) нелинейная деформационная модель сечения, нормального к продольной оси железобетонного элемента, на совместное действие изгибающего момента в двух осях и продольной сжимающей или растягивающей силы (данная модель применяется для расчета наиболее нагруженных средних частей пролетов плит, балок ригелей, а также расчета колонн);
2) нелинейная деформационная модель наклонного сечения на совместное действие изгибающего момента продольной и перерезывающей силы (применяется для моделирования предопорных зон балок, ригелей);
3) нелинейная блочная модель растянутого элемента с изгибом (является новой моделью для описания процесса деформирования растянутого элемента, описывает процесс перераспределения напряжений в блоках между трещинами нормального отрыва, в частности для моделирования процесса трещино-образования в многопустотных плитах безопалубочного формования в составе платформенных стыков сборно-монолитных перекрытий).
Рассматриваемые нелинейные модели описывают наиболее опасные участки строительных конструкций, как правило, испытывающие максимальные усилия и приводящие к разрушению конструкций. Нами разработаны:
1) эффективные методы решения прямых и обратных задач в рамках рассматриваемых нелинейных моделей;
2) алгоритмы и программные модули расчета рассматриваемых моделей;
3) проведена апробация моделей в рамках сопоставления результатов расчета и экспериментальных данных на широком классе конструктивных элементов и видов нагружения.
Рассмотрим далее в кратком изложении основные системы уравнений изучаемых моделей.
Нелинейная деформационная модель сечения, нормального к продольной оси железобетонного элемента
Нелинейная деформационная модель напряженно-деформированного состояния сечения, нормального к продольной оси, представляет собой систему сингулярных нелинейных уравнений интегрального вида:
М* = "Я(X, У )) (Х " X» )ёхёУ + N(Хое -х»);
с
Му =-Ц о(е( X, У)) (У - У» )ёхёу + Ы(уое - у»); (1)
с
N = Л о(е( х, у )) ёхёу,
с
где о(е(х, у)) - нормальные напряжения в точке (х, у) сечения в бетоне или арматуре; е(х, у) - продольные деформации в точке (х, у) сечения; у0 - расстояние от центра изгиба сечения, нейтральной точки, до оси у; и Уое - расстояние от места приложения продольного усилия до оси у; х0 - расстояние от центра изгиба сечения, нейтральной точки, до оси х; и х0е - расстояние от места приложения продольного усилия до оси х; N - продольная сила от действия внешней нагрузки; Мх, Му - изгибающие моменты в плоскости XX и УХ соответственно.
Уравнения совместности деформаций для точки (х, у) выглядит следующим образом: е(х,у) = е2 -—(у -уо) -—(х- хо) + дй(х, у) + дДх,у) + д2(х, у),
Гу г,
1 1
где--кривизна продольной оси элемента относительно оси у;--кривизна продольной оси элемента
г г
У х
относительно оси х; д0, д1, д2 - остаточные деформации от предварительного напряжения элемента до усиления, добавочные деформации от нагрузки до усиления и остаточные деформации от процесса динамического обжатия после усиления.
Математическая модель устанавливает отображение множества допустимых значений деформаций, определенное в метрическом пространстве деформаций, во множество значений усилий, определенное в метрическом пространстве усилий М: Я ® Я3. Сама модель лишена предположений относительно способа расчета напряженно-деформированного состояния. Из деформационной модели следует формулировка задачи расчета как задачи поиска такого расположения плоскости распределения деформаций
1111
в 3-х мерном пространстве е = <!—,—,е2 >, которое бы соответствовало действию на элемент вектора
I гх гу
нагрузки, заданного в трехмерном пространстве нагрузок Р = {Мг, Му, N1. Таким образом, решение задачи ищем в 3-х мерном пространстве.
Нелинейная деформационная модель наклонного сечения. Система уравнений исследуемой нами модели представлена в следующем виде:
Му = -1 о (е (у)) х(у)(у - у» )йу + N(у»е - у»);
N = | о(е( у)) х(у)йу; (2)
V = | Т( х, у) х( у)ёу,
где о(е(у)) - нормальные напряжения в слое на уровне у, сечения в бетоне или арматуре; е(у) - продольные деформации элементарной площадки бетона или арматуры (проекция на ось X суперпозиции главных сжимающих и главных растягивающих деформаций, полученная с учетом поворота элементарных площадок под действием перерезывающей силы); х(у) - функция формы; у0 - расстояние от центра изгиба сечения, нейтральной точки, до начала координат; у0е - расстояние от места приложения продольного усилия до начала координат; N - продольная сила от действия внешней нагрузки; V - перерезывающая сила; Му - изгибающий момент, действующий в плоскости И.
В модели постоянного (фиксированного) угла принимается предпосылка о том, что трещины, возникающие в железобетонном элементе в условиях плоского напряженного состояния, ориентированы по от-
ношению к продольной оси под некоторым постоянным углом 9, не зависящим от уровня нагружения и наличия новых трещин и определяемым углом наклона осей главных напряжений в элементе без трещин.
В модели переменного угла принимается допущение о том, что угол наклона 9 трещин и сжатых бетонных полос между наклонными трещинами на произвольном этапе после их образования совпадает с углом наклона главных напряжений а, значение которого изменяется при нагружении.
Уравнения равновесия и совместности деформаций для одноосного случая выглядят следующим образом (ось Z - совпадает с продольной осью элемента):
s z = ос2 cos2 (a) + acl sin2 (a) + pz aíZ, sy = sC2 sin2 (a) + Gcicos2 (a) + p y ov, ty = (-sC2 + sci) sin(a) • cos(a), e z = ec2 cos2 (a) + ec1 sin2 (a), e y = ec2 sin2 (a) + ec1 cos2 (a),
yz
= (-ec2 + ec1) sin(a) • cos(a),
где oz, oy, oc1, oc2, Tzy - соответственно проекции главных сжимающих и растягивающих напряжений на оси Z и Y, главные сжимающие и растягивающие напряжения, тангенциальные напряжения в бетоне; osz, osy - напряжения в арматуре; Sz, ey, yzy - относительные деформации; a - угол наклона осей главных напряжений по отношению к продольной оси.
Постановка задачи расчета данной модели методом двух сечений может быть представлена отображением:
< My1, My2, 5, N > ® < Sx, Sz, Sy, Yzy > ® V,
где 5 - расстояние между сечениями.
Нелинейная блочная модель растянутого элемента с изгибом
Данная нелинейная модель была построена, реализована программно и верифицирована экспериментально в рамках диссертационной работы А.Д. Лазовского [17].
Для любого поперечного сечения блока суммарное усилие в растянутом бетоне Nc и арматурном стержне Ns постоянно на всей длине блока. Тогда условие равновесия усилий с учетом правила знаков на элементарном участке в пределах бетонного блока имеет вид:
dSs,ad • As.ad + dsc • A,eff,red = 0
Исходя из условия равенства усилия в арматурном стержне диаметром 0, как результирующего от действия растягивающих нормальных напряжений as,ad, распределенных по площади его поперечного сечения Asad, и усилия от действия касательных напряжений ть по площади контакта арматурного стержня с бетоном, а также разности относительных деформаций арматуры Ss,ad и бетона Sct как взаимного смещения s, закона сцепления ть = f(s), связывающего касательные напряжения по площади контакта арматурного стержня с бетоном и их взаимное смещение, учитывающее поперечное сжимающее давление от вертикального прижатия, можно записать:
dss,ad = p0 •Wp,r • tb(s)dx;
(es,ad -ec )dx = ds;
Wp,r -tb = f (s).
Коэффициент Q.p,tr к ординатам принятого закона сцепления арматуры с бетоном для учета влияния давления вертикального прижатия ptr в окружающем арматуру бетоне средней прочности fcm,ad
W p,r = 1 - tanh
(
0,2-
0,1 f
> J cm
,ad
\
График, иллюстрирующий принятый закон сцепления с учетом влияния поперечного давления прижатия р(г, представлен на рисунке 3.
В случае достижения относительными деформациями арматуры значений, соответствующих пределу текучести, в,,асС > 8,У,аС, к ординатам принятого закона сцепления арматуры с бетоном вводится понижающий коэффициент йу:
« ь
а = 1,0-0,85• (1 -в~5а ),
где а =
е - е
,,аС ,у ,аС
ь =
/
^ уп
- коэффициенты; 8,М;аС - относительные деформации, соответствую-
щие пределу прочности (предельной растяжимости) дополнительной надопорной арматуры; /ут,асС - среднее значение предела текучести арматуры при растяжении; /т = к/утаС - среднее значение предела прочности арматуры при растяжении.
При 8,,аС < 8,у,аС коэффициент Qy = 1.
График, иллюстрирующий принятый закон сцепления с учетом влияния текучести арматуры, представлен на рисунке 4.
Рисунок 3. - График, иллюстрирующий закон сцепления с учетом влияния давления поперечного прижатия рг
Рисунок 4. - График, иллюстрирующий закон сцепления с учетом влияния относительных деформаций, соответствующих пределу текучести арматуры
Для решения задачи получения распределения относительных деформаций арматуры и бетона центральный бетонный блок в растянутой зоне поперечного сечения высотой Нф выделенный соседними трещинами нормального отрыва, разделяем п поперечными сечениями на (п — 1) интервалов длиной Дх таким образом, чтобы границы интервала Дх^ совпадали с границами вертикального шва между торцами многопустотных плит (рисунок 5). Учитывая при этом низкую адгезию бетона заполнения вертикального шва с бетоном плит, работу бетона на растяжение поперек шва в сечениях, совпадающих с его границами, в расчет не берем. Поэтому на интервале, равном длине вертикального шва Дх^, принимаем эффективную площадь поперечного сечения растянутой зоны в бетоне замоноличивания Асф,аС и диаграмму деформирования бетона замоноличивания.
Рисунок 5. - Схема разделения центрального бетонного блока на интервалы Ах
2
е
Таким образом, рассматриваемую модель можно описать следующей системой уравнений:
d _ ( )_ (N-a.-'A
dx
s,ad s,ad
A
c,eff ,red
(3)
d 4
—a ,, _ — dx s,ad 0
a ,, _ — W -W , -t(s).
s,ad sy, y p,tr V у
Как и для пролетных бетонных блоков, для центрального бетонного блока процесс трещинообра-зования проходит две стадии. Особенностью для центрального блока является уменьшение площади бетона растянутой зоны в месте вертикального шва (в работе на растяжение участвует только бетон замо-ноличивания площадью Acfd, что способствует появлению в конце первой стадии трещинообразования в месте контакта бетона плит и бетона замоноличивания поперечной трещины нормального отрыва). Для поперечных сечений с трещинами нормального отрыва (по краям выделенного блока) значения напряжений Gs,ad,(i) и os,ad,(r) известны, а значения смещений и S(r) - неизвестны.
Анализ результатов статистического моделирования. В силу нелинейности рассматриваемых моделей результатом статистического моделирования являются распределения выходных переменных модели (предельные значения внутренних усилий, момент трещинообразования, ширина раскрытия трещины, оценка деформаций: прогибы, перемещения), представляющие собой смеси вероятностных распределений. В частности, помимо сингулярности, связанной с появлением трещин, на изменение характера поведения модели под нагрузкой сказываются такие переходы, как достижение зоны текучести и предела текучести растянутой арматуры, причем отдельно для основного сечения и элементов усиления конструкции.
При нелинейном моделировании предопорных зон влияние на изменение характера выходных переменных оказывает также состояние поперечной арматуры и процесс возникновения наклонных трещин под воздействием перерезывающих сил.
Модель, представляющая собой смесь нормальных распределений, задается в виде:
k
p(x) _ XwiP(x|i),
где p(x | i) - нормальное распределение для i-го кластера; Wj - вес i-го кластера в смеси.
Анализируемая смесь вероятностных распределений характеризуется неопределенностью относительно количества распределений в смеси и ограниченностью объема выборки. Кроме этого, распределения в смеси сильно перекрываются, что делает неприменимым алгоритм разделения смеси на основе кластеризации методом k-средних.
Для решения задачи разделения смеси вероятностных распределений, обладающей вышеуказанными особенностями, не рекомендуется применение классического EM-алгоритма, основными недостатками которого является, во-первых, высокая чувствительность к выбору первого приближения, во-вторых, неустойчивость [1; 2]. Также следует отметить, что EM-алгоритм и его известные модификации требуют значительных вычислительных ресурсов, в частности, адресного пространства для хранения nxk скрытых переменных и необходимости пересчета значений всех скрытых переменных на каждой итерации алгоритма.
Имеются медианные и стохастические модификации EM-алгоритма, которые делают алгоритм более устойчивым к выбору первого приближения, однако не снимают проблемы, связанной с неопределенностью по количеству кластеров [1]. Нами предложен генетический алгоритм разделения смеси, варьирующий количество и параметры компонентов (GEM - Genetic EM).
Общие положения теории надежности. Целью анализа результатов статистического моделирования является поиск вероятностных характеристик функции случайных переменных Y = g(X).
В общем случае, когда речь идет о надежности строительных конструкций, характеристическая функция представляет собой разность случайной величины (в некоторых случаях - суперпозиции случайных величин), описывающей прочность конструкции, и случайной величины (суперпозиции случайных величин в силу линейного характера сочетания усилий), описывающей воздействие на конструкцию:
Y = g(X) = XR - Xs.
Тогда g(X) = 0 представляет собой определение предельного состояния, g(X) > 0 описывает область безопасной эксплуатации и g(X) < 0 - область разрушения конструкции. Вероятности этих событий можно интерпретировать как
P[g(X) > 0] = R - надежность, P[g(X) < 0] = 1 - R = Pf - вероятность разрушения.
В случае если модели входных переменных меняются со временем, имеем дело с вероятностным процессом, в рамках которого надежность оценивается как функция времени
R(t) = P[g(X(t) > 0)].
л
Переход к стандартному нормальному распределению. В результате статистического моделирования мы получаем распределение несущей способности и воздействия. Если исходить из предположения, что закон распределения этих случайных величин является нормальным, то разность этих случайных величин есть нормально-распределенная случайная величина с параметрами: Мя-я = Мл -
Дс
о
/о2 + о2
JR ' "S '
Тогда, перейдя к стандартному нормальному распределению, имеем
Pf = P[R - S < 0] = J /r_s (x)dx = Ф
0 mR-S l = F(-p),
где в - интерпретируется как индекс надежности.
Рассмотрение зависимости интеграла вероятности от базовых переменных модели. Для более детального описания вероятностной модели применяют переход от одномерной модели воздействия-сопротивления к многомерной модели с и-мерным вектором входных вероятностных параметров (прочностные характеристики материала, геометрические параметры изделий и т.п.).
В этом случае функцию g(X) рассматривают как функцию в и-мерном пространстве, а вектор X = (Х1, Х2,..., Хи)т - как и-мерный вектор входных случайных параметров.
В таком случае интеграл вероятности определяется как и-мерный интеграл:
Pf = J fx (x)dx.
д(х)£0
В случае если все параметры представлены независимыми случайными величинами, объединенная функция распределения плотности вероятности в и-мерном пространстве будет представлять собой произведение функций распределения плотностей вероятности параметров:
л (х)=п/- (х).
/=1
Переход и-мерного представления к стандартному нормальному распределению. Для упрощения процедуры вычисления искомого интеграла вероятности используют переход к стандартным нормальным распределениям. Такое допущение возможно, если указанный переход корректировать на каждом шаге итерационных процессов поиска наиболее вероятной точки области ошибки (разрушения).
В таком случае и-мерный интеграл вероятности будет определен как
Pf = J Ф„ (u)du = Ф(-Р).
q(u)£0
Для определения вероятности разрушения в пространстве стандартизованных нормальных распределений получили широкое распространение методы, основанные на разложении в ряд Тейлора функции g(u). В зависимости от того, сколько членов ряда Тейлора берется во внимание для построения итерационного процесса поиска наиболее вероятной точки границы g(u) = 0, различают FORM (First order reliability method) и SORM (Second order reliability method).
q(U) = „<„•) + - »>)' + \(U - »>) H ^ - u.)<U - „.f,
где H(u ) - Гессиан, V(u*) - градиент функции g(u) в точке u .
V(u*) =
дq dq
dU, dU2 д2 q д2 q
Э2 q 'dU1U„
H (u*) =
dUj2 дU1U2 д2 q д2 q
д 2 q dUiU„ д 2 q
дU„Ul дUяU,
R-S
f
\
u=u*
Поиск наиболее вероятной точки границы выживаемости или проектной точки (в зарубежных источниках MPP - Most probable point) осуществляется итерационно. В большинстве публикаций предлагается поиск методом простой итерации с критерием остановки по условию малости изменения оценки индекса надежности. В пространстве стандартизованных нормальных распределений индекс надежности определяется как расстояние от начала координат до проектной точки:
Р = minj |u||, q(u) = 0.
Методы FORM и SORM применимы в том случае, когда мы имеем возможность построить дифференцируемое аналитическое выражение функции g(u). Однако при использовании метода Монте-Карло для нелинейных моделей такое условие становится практически не выполнимым.
Для случая применения метода Монте-Карло оценка интеграла надежности будет представлять собой отношение количества численных экспериментов, в которых обеспечена прочность элемента, к общему количеству выполненных модельных расчетов:
nf
Pf .
f N
С учетом малости вероятности ошибки (10-4... 10-6) для оценки вероятности разрушения требуется объем эксперимента и обработки больших объемов данных.
Нами применен метод построения теоретических распределений по данным меньших объемов, и анализ результатов распределений как гауссовой смеси.
Учет нелинейных эффектов в процедуре оценки интеграла вероятности
Поскольку, как было показано, результирующее распределение несущей способности в силу влияния нелинейных эффектов является смесью нормальных распределений, то решением нелинейной задачи оценки индекса надежности будет являться оценка суммарной вероятности разрушения от каждой компоненты смеси.
Учитывая линейность преобразования Лапласа, сумма вероятностей отказа будет представлена следующим образом:
У1Ф(-р1) + У2Ф(-р2) + ... + УтФ(-рт) = Ф(^1р1 - V2p2 ... -Vmpm),
m
p=E v,p,,
1=1
где V1, ., Vm - веса компонентов смеси.
Заключение. Предложен комплексный метод оценки надежности строительных конструктивных элементов на основе нелинейных вероятностных моделей наиболее опасных участков и расчетных ситуаций. Сделан вывод о необходимости проверки гипотезы о том, что функция предельного состояния является гауссовой смесью вероятностных распределений, а в случае состоятельности данной гипотезы необходимо выполнение процедуры разделения смеси на компоненты. Предложена интерпретация данного факта как эффекта влияния нелинейных свойств и сингулярности моделей на распределение функции предельного состояния.
ЛИТЕРАТУРА
1. Королёв, В.Ю. ЕМ-алгоритм, его модификации и их применение к задаче разделения смесей вероятностных распределений. Теоретический обзор / В.Ю. Королёв. - М. : Изд-во ИПИРАН, 2007. - 94 с.
2. Сеточные методы разделения смесей вероятностных распределений и их применение к декомпозиции волатильности финансовых индексов / В.Ю. Королёв [и др.] // Информатика и её применения. -2008. - Т. 2. - Вып. 2. - С. 3-18.
3. Rackwitz, R. Structural Reliability under Combined Random Load Sequences / R. Rackwitz and B. Fiessler // Computers and Structures. - 1978. - 9. - P. 489-494.
4. Allen, T.M. Calibration to Determine Load and Resistance Factors for Geotechnical and Structural Design / T.M. Allen, A.S. Nowak, R.J. Bathurst // Transportation Research Circular E-C079, Transportation Research Board. - Washington, D.C., 2005.
5. Cornell, C.A. Bounds on the Reliability of Structural Systems / C.A. Cornell // American Society of Civil Engineers : Journal of the Structural Division ASCE. - 1967. - Vol. 93, No. ST. - February. - P. 171-200.
6. Ржаницын, А.Р. Теория расчета строительных конструкций на надежность / А.Р. Ржаницын. - М. : Стройиздат, 1978. - 239 с.
7. Vrouwenvelder, Ton Reliability Based Code calibration. The use of the JCSS Probabilistic Model Code / Vrouwenvelder Ton // Joint Committee of Structural Safety Workshop on Code Calibration, March 21/22, Zurich.
8. EN 1990 Eurocode - Basis of structural design. CEN 2002.
9. ISO 2394 General principles on reliability for structures, ISO 1998.
10. CEN, «ENV 1993-1-1 : Eurocode 3 Tiel 1-1 : Annex Z - Determination of design resistance from tests» // European Committee for Standardisation, Brussels, 1993.
11. Глухов, Д.О. Оценка параметров теоретических распределений случайных величин в расчете надежности строительных конструкций / Д.О. Глухов, Д.Н. Лазовский // Вестник Полоцкого государственного университета. Сер. F, Строительство. Прикладные науки. - 2008. - № 6. - С. 2-12.
12. Глухов, Д.О. Генераторы случайных чисел для расчета надежности строительных конструкций / Д.О. Глухов, Т.М. Глухова // Вестник Полоцкого государственного университета. Сер. С, Фундаментальные науки. - 2008. - № 3. - С. 60-68.
13. Глухов, Д.О. Численный анализ расчетной модели нормального сечения проекта СНБ 2.03.01-98 / Д.О. Глухов, Д.Н. Лазовский // Совершенствование железобетонных конструкций, оценка их состояния и усиление : сб. материалов респ. науч.-техн. конф. ; под ред. Т.М. Пецольда. - Минск : Техно-принт, 2001. - С. 97.
14. Глухов, Д.О. Программы БЕТА и RADUGA для расчета строительных конструкций / Д.О. Глухов, Д.Н. Лазовский, О.Н. Лешкевич // Архитектура и строительство. - Минск : Полипринт, 2001. -№ 6 (148). - С. 28-30.
15. Глухов, Д.О. Оценка вероятностных характеристик распределения эффектов от внешнего воздействия и функции предельного состояния в расчетах железобетонных элементов / Д.О. Глухов, Т.М. Глухова, Д.Н. Лазовский // Вестник Полоцкого государственного университета. Сер. F, Строительство. Прикладные науки. - 2009. - № 6. - С. 35-40.
16. Применение мягких вычислений для сглаживающей аппроксимации сложных сингулярных зависимостей / Д.О. Глухов [и др.] // Вестник Полоцкого государственного университета. Сер. С, Фундаментальные науки. - 2012. - № 12. C. 2-5.
17. Лазовский, А.Д. Алгоритм расчета сопротивления изгибу многопустотных плит перекрытий безопалубочного формования в составе платформенных стыков зданий / А.Д. Лазовский, Д.О. Глухов, Т.М. Глухова // Вестник Полоцкого государственного университета. Сер. F, Строительство. Прикладные науки. - 2017. - № 8. - С. 82-89.
18. Haldar, A. Probability, Reliability, and Statistical Methods in Engineering Design / A. Haldar and S. Mahadevan, John Wiley. - New York, 2000.
19. Xiaoping Du. Second-Order Reliability Method with First-Order Efficiency [Electronic resource] / Xiaoping Du and Junfu Zhang // Proceedings of the ASME Design Engineering Technical Conference (2010). -[Electronic resource]. - Mode of access: http://works.bepress.com/xiaoping-du/75.
20. Au, S.K. Estimation of small failure probabilities in high dimensions by subset simulation / S.K. Au, J.L. Beck // Probabilistic Engineering Mechanics. - 2001. - 16(4). - P. 263-277.
21. Au, S.K. Importance sampling in high dimensions / S.K. Au, J.L. Beck // Structural Safety. - (2003а). -25(2). - P. 139-163.
22. Au, S.K. Subset Simulation and its application to seismic risk based on dynamic analysis / S.K. Au, J.L. Beck // Journal of Engineering Mechanics. - (2003b). - 129(8). - P. 1-17.
23. Schueller, G.I. A critical appraisal of reliability estimation procedures for high dimensions / G.I. Schueller, H.J. Pradlwarter, P.S. Koutsourelakis // Probabilistic Engineering Mechanics. - (2004). - 19. - P. 463-474.
Поступила 01.12.2017
PROBABILITY ANALYSIS OF THE STRUCTURAL REINFORCED CONCRETE ELEMENT
D. GLUKHOV, R. BOGUSH, Е. LAZOVSKY, T. GLUKHOVA
The paper presents a comprehensive approach to assessing the reliability of structural components based on nonlinear probability models of the most dangerous areas and design situations. It is concluded that it is necessary to test the hypothesis that the resistance probability distribution function is a Gaussian mixture of probability distributions, and in the case of the validity of this hypothesis, the procedure for separating of the mixture into components is performed. An interpretation of this fact is proposed as an effect of the influence of nonlinear properties and the singularity of models.
Keywords: structural reinforced concrete element, evaluation of reliability, non-linear probabilistic models.
