Архив рубрики ‘Надежность’
Стыд-то какой.
Читал ДБН В.1.2-14-2009, “Принципы обеспечения надежности зданий и сооружений”, типа наш отечественный аналог еврокода EN 1990:2002.
Большего позорища видеть в жизни не приходилось.
Если вы не видели, как выглядит надгробная плита на могиле технических нормативов двадцать первого века, то вот она, натурально.
А вот так выглядит ссылка на еврокод EN 1990. Девять-ноль, бляха-муха. В государственном нормативном документе. Это ад и позорище. Внимание на правописание ещё.
О наполнении можно не говорить. Нет, хорошо, конечно, когда есть куда сослаться, если тебе заявляют, что термины, которые ты используешь, непонятны. Но содрать две трети содержания из наших атомных НД и делать вид, что так и надо – это слишком.
Я пошел и посмотрел список авторов. Две трети фамилий знакомые, половина – люди, с которыми я вживую общался.
Ну, понятно, у нас все занятые стали. Написание нормативов можно отдать на растерзание секретаршам.
P.s. Ну его в баню, два доменных имени и два блога. Буду писать на bausk.wordpress.com и на английском, и на русском, буржуи потерпят как-нибудь.
Очередной англоязычный сайт…
…с материалами по инженерному делу.
www.thestructuralengineer.info
Там есть замечательная подборка ссылок. При всем моем скепсисе по отношению к подобным ссылкосборникам, этот – удивительно качественный. Особенно интересны разделы про иностранные стандарты и (как обычно) надежность.
Matlab: статистические штучки для получения реализаций случайной величины
Так-с. Исходные данные. Предполагается, что читатель знаком с:
- методом Монте-Карло (вкратце: случайные параметры -> много-много расчетов посредством произвольной расчетной модели -> вероятностный результат, надежность и прочие радости)
- и с тем, что такое функция распределения и функция плотности случайной величины.
Предположим, что у нас есть случайная прочность R с параметрами m=350, s=30. Это выглядит так:
cdf('norm', x, Mean1, StdDev1)pdf('norm', x, Mean1, StdDev1)
Но для задачи вероятностного расчета методом Монте-Карло нам нужно средство генерации множества реализаций нашей случайной прочности.
Для равномерного распределения такое средство одинаково в любом языке программирования и называется RND. И в Екселе тоже есть.
Для любого другого распределения для генерации случайных чисел нам нужна функция, обратная CDF, отвечающая на вопрос “какое число соответствует вероятности непревышения такой-то?”. Такой способ получения случайных чисел называется методом обратного преобразования.
Откуда нам взять такого зверя? На помощь приходит замечательный статистический пакет Stixbox, о котором я обещал написать :
>> MonteCarloData1 = rnorm(1000, 350, 30); >> Mean1 = mean(MonteCarloData1) Mean1 =349.9133 >> StdDev1 = std(MonteCarloData1) StdDev1 =30.7880
Аналогичные функции есть в этом пакете и для других известных распределений. Для задач вероятностного анализа расчетных моделей конструкций более чем достаточно. Спасибо за этот замечательный инструмент Андерсу Хольтсбергу.
Про 0,99865
Как показало былинное обсуждение «коэффициента спокойного сна» на dwg.ru, насчет критериев обеспеченности расчетных значений существует значительная неразбериха. Не говоря уже про вероятностный анализ.
Итак. Откуда берутся упомянутые в первой ссылке величины обеспеченности 0,99865 и 0,997. (Предполагается, что читатель знает, что такое среднее значение случайной величины и её стандартное отклонение.)
- 0,997 – вероятность того, что реализация случайной величины с нормальным распределением окажется в промежутке между значением “среднее минус три стандартных отклонения” и “среднее плюс три стандартных отклонения”. Этот промежуток будет называться доверительным интервалом.
- 0,99865 – вероятность того, что эта реализация будет больше, чем значение “среднее минус три стандартных отклонения”.
Первое можно получить с помощью так называемой функции ошибки (используем Matlab/Octave, как обычно):
>> erf(3/sqrt(2)) ans = 0.9973
Второе, как уже обсуждалось в более ранней записи,
>> 1 - cdf('norm', 1000-3*50, 1000, 50)
ans = 0.998650101968370
Как видим, 0,997 пришло к нам из статистики и еретического учения маркетинга сотоварищи. 0,99865, наоборот, есть правильная инженерная цифирь и, попросту говоря, является обеспеченностью (вероятностью), например, расчетной прочности бетона.
В нормах, понятно, все это отражено крайне мутными словами, и расчетное значение вычисляется на основе нормативного (с обеспеченностью 0,95) при помощи частных коэффициентов. Все это приводит к тому, что доценты нашего института не знают, что такое “обеспеченность” и в чем разница между коэффициентами перегрузки и коэффициентами безопасности (правильнее – частные коэффициенты надежности).
На картинке показана:
темно-красным – плотность распределения некоторой прочности с матожиданием 350, стандартным отклонением 30;
красным – гистограмма дискретных реализаций этой величины (1000 штук);
синим – расчетная и нормативная прочности с обеспеченностью 0,99865 и 0,95.
