Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,685

МЕТОД ОЦЕНКИ ВЕРОЯТНОСТИ КАТАСТРОФЫ СБОРКИ ДЛЯ СЛУЧАЯ, КОГДА УПРАВЛЯЮЩИЕ ПАРАМЕТРЫ ЯВЛЯЮТСЯ СЛУЧАЙНЫМИ ФУНКЦИЯМИ

Питухин А.В. 1 Скобцов И.Г. 1
1 ГОУ ВПО «Петрозаводский государственный университет»
В работе приведен метод оценки показателей надежности с позиций теории катастроф. В первой части статьи представлено введение в теорию катастроф, рассмотрена катастрофа сборки и проведен анализ критических значений управляющих параметров. Параметры управления рассматривались как стационарные случайные процессы (функции) и ставилась задача о выбросах. На следующем этапе исследований определялось среднее число пересечений нулевого уровня стационарным случайным процессом. Полученные аналитические зависимости дают возможность оценки вероятности безотказной работы, вероятности отказа (вероятности катастрофы сборки), средней наработки до отказа. Для определения характеристик случайного процесса (математического ожидания и дисперсии) использовался метод статистической линеаризации. Предложенный подход позволяет произвести статистический анализ положений равновесия вблизи критических точек. Полученные зависимости могут быть применены для оценки вероятности безотказной работы различных элементов конструкций машин и оборудования с позиций теории катастроф.
теория катастроф
катастрофа сборки
случайная функция
показатели надежности
1. Арнольд В.И. Теория катастроф. – М.: Наука, 1990. – 128 с.
2. Питухин А.В. Вероятностно-статистические методы механики разрушения и теории катастроф в инженерном проектировании. – Петрозаводск: Изд-во ПетрГУ, 1998. – 304 с.
3. Питухин А.В. Методы теории катастроф при проектировании элементов конструкций машин и оборудования лесного комплекса // Известия вузов «Лесной журнал». – 2007. – № 2. – С. 58–65.
4. Постон Т. Теория катастроф и ее приложения / Т. Постон, И. Стюарт; пер. с англ. – М.: Мир, 1980. – 608 с.
5. Свешников А.А. Прикладные методы теории случайных функций. – М.: Мир, 1962. – 463 с.
6. Томпсон Дж. М. Т. Неустойчивости и катастрофы в науке и технике; пер. с англ. – М.: Мир, 1985. – 254 с.
7. Gilmore R. Structural stability of the phase transition in Dicke-like models // J. Math. Phys. A. – 1977. – Vol. 18. – P. 17–22.
8. Gilmore R. Catastrophe Theory for Scientists and Engineers. – New York: Dover, 1993.
9. Rise S.O. Mathematical analysis of random noise // Bell System Tech. J. – 1945. – Vol. 24. – P. 46–156.
10. Thom R. Structural Stability and Morphogenesis: An Outline of a General Theory of Models. – Reading, MA: Addison-Wesley, 1989.
11. Whitney H. On singularities of Mappings of Euclidean Spaces I. Mappings of the Plane into the Plane // Ann. Math. – 1955. – Vol. 62. – P. 374–410.
12. Zeeman E.C. Catastrophe theory: Selected Papers 1972–1977. – Reading, MA: Adisson-Wesley, 1977.

Теория катастроф как раздел математики начала формироваться еще в середине ХХ века на основе теории особенностей гладких отображений, заложенной математиком Хасслером Уитни, и теории бифуркаций динамических систем А. Пуанкаре и А.А. Андронова. Теория особенностей представляет собой обобщение исследования функций на максимум и минимум, при этом в теории Уитни [11] функции заменены отображениями – наборами нескольких функций нескольких переменных. Теория бифуркаций (раздвоений) изучает качественные перестройки различных объектов при изменении параметров, от которых они зависят. Катастрофами называются скачкообразные изменения, возникающие в виде внезапного ответа системы на плавное изменение внешних условий [1].

Создателями теории катастроф являются французский математик Р. Том [10] и российский математик В.И. Арнольд [1]. Большой вклад в ее развитие внесли Т. Постон, И. Стюарт [4], Дж.М.Т. Томпсон [6], Р. Гилмор [7, 8], К. Зиман [12].

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

pityh01.wmf, (1)

и многообразие M катастрофы задается уравнением

pityh02.wmf (2)

Рис. 1 иллюстрирует поведение поверхности равновесия (или многообразие катастрофы) с ее проекцией на плоскость ab, определяющую управляющие параметры a и b.

Точка (a,b), медленно изменяя свои координаты a и b, описывает некоторую траекторию на плоскости ab. Тогда наблюдаемое положение равновесия пройдет путь на поверхности M, лежащий над путем в плоскости ab. Из-за сборки на поверхности равновесия M этому пути, возможно, придется перескакивать с одного листа поверхности на другой. Это объясняется тем, что точки многообразия M, расположенные на внутренней поверхности складки, соответствуют неустойчивому состоянию системы. Такой очень быстрый скачок системы (катастрофа) наступает только при выходе из области I, так как у системы не остается другого выбора (принцип максимального промедления (perfect delay) Р. Тома [10]). Таким образом, гладкие изменения переменных управления a и b могут повлечь за собой разрывные изменения переменной состояния x, вызывая катастрофические прыжки.

pic_3.tif

Рис. 1. Катастрофа сборки

Кубическое уравнение (2) имеет от одного до трех вещественных корней. Природа этих корней зависит от дискриминанта

D = 4a3 + 27b2. (3)

Пространство управляющих параметров, образующих плоскость ab, разбивается на пять подмножеств: B1, B2, P, I, E, изображенных на рис. 2. Линия, образуемая множествами B1, B2, и P, соответствует условию D = 0, т.е. 4a3 + 27b2 = 0.

При этом, если a ≠ 0 и b ≠ 0, то точка (a, b) лежит в B1 или в B2, и уравнение (2) имеет три вещественных корня, два из которых совпадают между собой. Точка P соответствует a = b = 0 и имеются три совпадающих вещественных корня, равных нулю. Точка P называется точкой возврата.

Если (a, b) ∈ I, то имеется три различных вещественных корня и D < 0.

Если (a, b) ∈ E, то имеются только один вещественный корень (D > 0).

Катастрофа происходит, когда траектория точки (a, b) покидает область I и при этом D меняет знак с отрицательного на положительный.

pic_4.tif

Рис. 2. Пространство управляющих параметров

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

Рассмотрим катастрофу сборки (1) со стохастических позиций. Переменные управления a и b в общем случае изменяются во времени, и характеристика состояния будет определяться случайным процессом D(t). Таким образом, необходимо решать задачу о выбросах случайного процесса из области I. При этом вероятность возникновения катастрофы

P(t) = P{D(t) > 0}.

Пусть D(t) – дифференцируемый случайный процесс,

pityh03.wmf

Нас интересует вероятность того, что реализация случайного процесса D(t) пересечет нулевой уровень (рис. 3). Для того чтобы выброс произошел, необходимо осуществление двух событий:

1) в момент времени t, предшествующий выбросу, ордината D(t) должна быть меньше нуля (D(t) < 0);

2) в момент времени t + dt, следующий за моментом t (в тот момент, когда произойдет выброс) ордината случайной функции должна быть положительной, т.е. D(t + dt) > 0.

pic_5.tif

Рис. 3. Реализация случайного процесса D(t)

Таким образом, вероятность того, что выброс произойдет в бесконечно малом интервале времени (t, t + dt) [9]

pityh04.wmf (4)

где p0(t) – временная плотность вероятности.

Выразим в выражении (4) D(t + dt) через D(t). Для этого разложим D(t + dt) в ряд Тейлора в окрестности точки t

pityh05.wmf

Отсюда

pityh06.wmf

После переноса V(t)dt в правую часть неравенства получим

D(t) > 0 – V(t)dt.

С учетом зависимости (4)

– V(t)dt < D(t) < 0.

вероятность выброса запишется в виде

pityh07.wmf

Здесь f(d, v; t) – двумерный закон распределения ординат D и скоростей V.

После применения теоремы о среднем

pityh08.wmf

получим

pityh09.wmf

Здесь pityh10.wmf – временная плотность для вероятности выброса.

Временную плотность p0(t) можно трактовать как среднее число пересечений случайным процессом D(t) нулевого уровня в единицу времени.

В общем случае среднее число пересечений уровня 0 за время τ (математическое ожидание числа выбросов)

pityh15.wmf (5)

Для стационарного процесса плотность распределения ординат и скоростей не зависит от времени, т.е. f(d, v; t) = f(d, v) и, следовательно, p0(t) = p0. Отсюда

pityh16.wmf (6)

Во многих задачах практический интерес представляет вариант, при котором среднее число выбросов за данный промежуток времени достаточно мало и можно считать появления последовательных выбросов независимыми «редкими» событиями. В этом случае число появлений выбросов можно считать приближенно подчиняющимся закону распределения Пуассона [5], при этом единственным параметром, входящим в закон распределения, является математическое ожидание числа выбросов, определяемое формулами (5) или (6). Тогда в общем случае вероятность безотказной работы определится

pityh17.wmf (7)

для стационарного процесса

pityh18.wmf (8)

Имеет место и другой подход, изложенный в работах А.В. Питухина [2, 3]. При этом для оценки вероятности катастрофы сборки предложены аналитический метод и метод статистической линеаризации для варианта, когда управляющие параметры a и b являются случайными величинами. Воспользуемся методом статистической линеаризации для оценки математического ожидания и дисперсии случайного процесса D(t) в случае, когда переменные управления являются стационарными случайными функциями (процессами). Исходя из формулы (3), получим

pityh11.wmf (9)

pityh12.wmf (10)

где pityh13.wmf – математические ожидания стационарных случайных процессов a(t) и b(t); pityh14.wmf – дисперсии стационарных случайных процессов a(t) и b(t).

Имея значения оценок (9) и (10) и полагая вид закона распределения дискриминанта D известным, применяя метод моментов Пирсона, можно определить неизвестные значения параметров θ1, θ2 предполагаемого распределения F(d, θ1, θ2).

Тогда вероятность безотказной работы определится

pityh19.wmf (11)

Вероятность отказа (вероятность катастрофы сборки)

pityh20.wmf (12)

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

Работа выполнена при поддержке Программы стратегического развития (ПСР) Петрозаводского государственного университета в рамках реализации комплекса мероприятий по развитию научно-исследовательской деятельности на 2012–2016 гг.

Рецензенты:

Григорьев И.В., д.т.н., профессор кафедры технологии лесозаготовительных производств Санкт-Петербургского государственного лесотехнического университета, г. Санкт-Петербург;

Заика Ю.В., д.ф.-м.н., профессор, заведующий лабораторией моделирования природно-технических систем, Институт прикладных математических исследований Карельского научного цента РАН, г. Петрозаводск.

Работа поступила в редакцию 27.01.2014.


Библиографическая ссылка

Питухин А.В., Скобцов И.Г. МЕТОД ОЦЕНКИ ВЕРОЯТНОСТИ КАТАСТРОФЫ СБОРКИ ДЛЯ СЛУЧАЯ, КОГДА УПРАВЛЯЮЩИЕ ПАРАМЕТРЫ ЯВЛЯЮТСЯ СЛУЧАЙНЫМИ ФУНКЦИЯМИ // Фундаментальные исследования. – 2014. – № 1. – С. 24-27;
URL: https://fundamental-research.ru/ru/article/view?id=33496 (дата обращения: 09.12.2022).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.685