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

АЛГОРИТМ РАСЧЕТА ДОВЕРИТЕЛЬНЫХ ИНТЕРВАЛОВ КРИВОЙ МАЛОЦИКЛОВОЙ УСТАЛОСТИ

Сызранцев В.Н. 1 Сызранцева К.В. 1 Ильиных В.Н. 1
1 ФГБОУ ВПО «Тюменский государственный нефтегазовый университет»
Для прогнозирования долговечности деталей с заданной вероятностью неразрушения, подвергающихся в условиях эксплуатации малоцикловому циклическому деформированию, использована разработанная в кинетической теории усталости математическая модель, на основе которой выполняется обработка результатов усталостных испытаний образцов. Показано, что два параметра модели: предел прочности и число циклов до верхней точки перегиба кривой малоцикловой усталости, связанное с пределом прочности нелинейной зависимостью, – являются величинами случайными, в общем случае с неизвестными законами распределения. В результате левые границы доверительных интервалов кривой малоцикловой усталости, соответствующие заданной вероятности (1–5 %) неразрушения образцов и используемые для прогнозирования их долговечности, рассчитать возможным не представляется. В работе предложен оригинальный алгоритм определения параметров модели кривой усталости, проходящей через любую экспериментальную точку; с использованием алгоритма который в процессе компьютерного моделирования рассчитывается серия кривых усталости длиной, равной общему числу экспериментальных точек. Установленная серия кривых усталости позволяет для любого фиксированного уровня действующих напряжений сформировать выборку случайной величины, – числа циклов до разрушения образцов, необходимую для расчета ее квантильных оценок при построении границы доверительного интервала. Для восстановления неизвестной функции плотности распределения этой случайной величины в работе использован математический аппарат непараметрической статистики, обеспечивающий решение задачи независимо от сложности закона распределения исследуемой случайной величины. Результаты работы проиллюстрированы на примере определения границ доверительных интервалов данных малоцикловых испытаний образцов из гибких труб.
циклическое деформирование
малоцикловая усталость
предел прочности
вероятность разрушения образцов
методы непараметрической статистики
границы доверительных интервалов
1. Почтенный Е.К. Кинетическая теория механической усталости и ее приложения. – Минск: Наука и техника, 1973. – 213 с.
2. Степнов М.Н. Статистические методы обработки результатов механических испытаний: Справочник. – М.: Машиностроение, 1985. – 232 с.
3. Сызранцев В.Н. Оценка безопасности и прочностной надежности магистральных трубопроводов методами непараметрической статистики / В.Н. Сызранцев, В.В. Новоселов, П.М. Созонов, С.Л. Голофаст. – Новосибирск: Наука, 2013 – 172 с.
4. Сызранцев В.Н. Расчет прочностной надежности изделий на основе методов непараметрической статистики / В.Н. Сызранцев, Я.П. Невелев, С.Л. Голофаст. – Новосибирск: Наука, 2008 – 218 с.
5. Сызранцева К.В. Расчет прочностной надежности деталей машин при случайном характере внешних нагрузок. – Тюмень: ТюмГНГУ, 2011. – 88 с.
6. Сызранцев В.Н., Сызранцева К.В., В.Н.Ильиных. Обработка данных малоцикловых испытаний на основе кинетической теории усталости // Фундаментальные исследования. – 2015. – № 3. – С. 139–143.

Основой решения задач прогнозирования долговечности работы деталей и конструкций машин, оценки их прочностной надежности является кривая усталости (кривая Велера) [2], связывающая число циклов деформирования (N) с уровнем действующего в опасном месте детали напряжения (σ). Среди математических моделей, используемых для описания кривых усталости, наиболее перспективной является полуэмпирическая модель, разработанная в рамках кинетической теории механической усталости [1, 6], учитывающая процессы накопления усталостных повреждений в материале при циклическом деформировании деталей. В частности, для кривой усталости в малоцикловой области математическая модель имеет следующий вид [1, 6]:

syzran01.wmf (1)

где σВ – предел прочности материала; ϑ – угол наклона кривой усталости в системе координат lgN – σ; N и H соответственно число циклов нагружения и число циклов деформирования до верхней точки перегиба кривой малоцикловой усталости, рассчитываемые по зависимостям

syzran02.wmf

syzran03.wmf (2)

в которых обозначено: Q – коэффициент выносливости; σr – предел выносливости детали при коэффициенте асимметрии цикла r; σrT – циклический предел текучести (ниже его уровня следы пластической деформации даже после нескольких миллионов циклов нагружения отсутствуют).

Алгоритмы и процедуры расчета значений параметров Q, σr, σrT, ϑ модели (1) на основе совокупности данных σi, Ni, syzran04.wmf, полученных в процессе испытаний n образцов на долговечность, рассмотрены в работе [6]. В результате реализации этих алгоритмов определяется кривая усталости в виде регрессионной зависимости N = N(σ), соответствующая 50 % вероятности разрушения (неразрушения) образцов. Естественно, что использовать такую зависимость для прогнозирования долговечности деталей не имеет смысла. В связи с изложенным важнейшим практическим приложением результатов обработки данных усталостных испытаний образцов является не сама кривая усталости, а ее левые (нижние) границы доверительных интервалов, соответствующие, например, вероятности разрушения образцов 1 или 5 %. Цель настоящей статьи состоит в разработке методики расчета точек границ доверительного интервала для кривой малоцикловой усталости в виде (1).

Расчет границ доверительного интервала кривой малоцикловой усталости

Расчет границ доверительного интервала кривой усталости требует знания функции плотности распределения f(N) или f(lgN) при σ = const. Традиционно [2] принимают, что эти функции соответствуют нормальному (логнормальному) закону распределения случайных величин. Однако, как показано в работах [3, 4, 5], функции f(N) и f(lgN) являются существенно более сложными, не описываемыми с помощью законов, исследованных в теории параметрической статистики.

Для решения задачи в настоящей статье реализована идея статистического моделирования, впервые предложенная в работе [4] для получения выборки предела выносливости с целью восстановления его функции плотности распределения методами непараметрической статистики [3, 5].

Из анализа зависимости (1) следует, что она содержит два параметра, природа которых случайная, – это предел прочности (σВ) и число циклов до верхней точки перегиба кривой малоцикловой усталости (H), которое относительно σВ описывается нелинейной зависимостью (2). На момент решения рассматриваемой задачи параметры syzran05.wmf, syzran06.wmf, Q = Q*, ϑ = ϑ* известны. Воспользуемся данными разрушения образцов на разрывной машине, на основе которых определим статистические характеристики σВ: среднее значение syzran07.wmf и границы его доверительного интервала, например, для вероятности 99 % syzran08.wmf и syzran09.wmf.

Подставляя в выражение (2) syzran10.wmf, syzran11.wmf, syzran12.wmf, рассчитаем математическое ожидание syzran13.wmf числа циклов H, а также границы его 99 % доверительного интервала: syzran14.wmf и syzran15.wmf.

Поставим задачу получения зависимости для кривой усталости, которая проходила через любую экспериментальную точку σi = const, Ni = const syzran16.wmf и находилась в пределах установленных границ доверительного интервала. Введем безразмерную величину χ, используя которую зададим текущие значения σВ и H:

syzran17.wmf

syzran18.wmf (3)

Подставляя функции σВ(χ) и H(χ) в выражение (1), получим

syzran19.wmf (4)

Входя в это выражение при σ = σi = const и N = Ni = const, имеем трансцендентное уравнение относительно одной переменной χ:

syzran20.wmf (5)

Решая данное уравнение любым численным методом для каждой пары значений σi, Ni, syzran21.wmf экспериментальных данных, определим выборку χi, syzran22.wmf, которая позволяет получить серию длиной n кривых малоцикловой усталости:

syzran23.wmf (6)

Данная совокупность кривых усталости позволяет сформировать выборку syzran24.wmf, syzran25.wmf при σ = σ* = const:

syzran26.wmf (7)

необходимую для решения задачи восстановления неизвестной функции плотности распределения случайной величины N или lgN при σ = σ* = const.

Пусть требуется восстановить функцию плотности распределения числа циклов до разрушения fN(N) при фиксированной величине напряжения σ = σ* = const. Поскольку априори закон распределения случайной величины N неизвестен, воспользуемся математическим аппаратом непараметрической статистики [4], успешно применяемым в последнее время для решения подобных задач [3, 4, 5]. Исходной информацией для определения функции fN(N) является совокупность значений syzran27.wmf, syzran28.wmf, рассчитанная по зависимости (7).

Для восстановления функции fN(N) воспользуемся методом Парзена – Розенблатта [3, 4]. Следуя этому методу, неизвестная функция плотности fN(N) записывается в виде

syzran29.wmf (8)

где K(N) – ядерная функция (ядро); hN – параметр «размытости».

На основании работ [3, 4, 5] воспользуемся ядерной функцией с нормальным ядром. В этом случае функция fN(N) описывается выражением

syzran30.wmf (9)

В работе [4] показано, что оптимальное значение syzran31.wmf рассчитывается по формуле

syzran32.wmf (10)

где syzran33.wmf syzran34.wmf

Имея функцию fN(N), требуемая по условиям обработки данных усталостных испытаний при σ = σ* = const величина квантиля syzran35.wmf определяется в результате решения относительно syzran36.wmf численным методом следующего уравнения:

syzran37.wmf (11)

Установленное в процессе реализации изложенного алгоритма значение syzran38.wmf при заданной величине напряжения σ* = const определяет точку границы доверительного интервала кривой усталости (1), соответствующего вероятности разрушения образцов 1 %. Для расчета других точек этого доверительного интервала достаточно повторить рассмотренную процедуру для напряжений σ* = const в требуемом диапазоне изменения.

Реализация разработанных алгоритмов на примере обработки данных малоцикловых испытаний образцов гибкой трубы HS-80

В работе [6] представлены результаты испытаний на долговечность прямоугольных образцов, вырезанных из гибкой трубы HS-80. Для этих данных осуществим расчет границ доверительных интервалов. На основе обработки данных растяжения образцов на разрывной машине получено

syzran39.wmf syzran40.wmf

syzran41.wmf syzran42.wmf

syzran43.wmf syzran44.wmf

По формулам (2) рассчитаем значения

syzran45.wmf syzran46.wmf

syzran47.wmf syzran48.wmf

syzran49.wmf syzran50.wmf

Обратимся к уравнению (5). Решая его численным методом n раз для каждой пары значений σi, Ni, определим выборку безразмерной величины χi, syzran51.wmf. После чего для любой фиксированной величины напряжения σ = σ* = const по выражению (7), используя массив χi, syzran52.wmf, рассчитаем выборку syzran53.wmf, syzran54.wmf. Воспользовавшись математическим аппаратом непараметрической статистики на основе выборки syzran55.wmf, syzran56.wmf, восстановим неизвестную функцию плотности распределения fN(N*), описываемую зависимостью (9). В качестве примера на рис. 1 показана гистограмма распределения syzran57.wmf, syzran58.wmf и функция fN(N*), – в виде (9), при σ = σ* = 250 МПа. Для расчета границ доверительных интервалов, например при вероятности разрушения образцов 0,5; 2,5; 5 и 50 %, достаточно для ряда фиксированных величин напряжений реализовать рассмотренную процедуру восстановления fN(N*) и рассчитать путем решения интегрального уравнения (11) соответствующие квантильные оценки чисел циклов. Результаты таких выполненных расчетов представлены на рис. 2 в системе координат lgN – σ.

pic_22.wmf

Рис. 1. Функция плотности распределения fN(N*)

pic_23.wmf

Рис. 2. Границы доверительных интервалов в системе координат lgN – σ

Основные результаты работы

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

2. Разработанные методики и вычислительные алгоритмы проиллюстрированы на примере обработки данных малоцикловых усталостных испытаний образцов гибких труб HS-80.

Рецензенты:

Быков И.Ю., д.т.н., профессор кафедры «Машины и оборудование нефтяной и газовой промышленности», ФГОУ ВПО «Ухтинский государственный технический университет», г. Ухта;

Плеханов Ф.И., д.т.н., профессор, директор Глазовского инженерно-экономического института (филиала), ФГОУ ВПО «Ижевский государственный технический университет им. М.Т. Калашникова», г. Глазов.

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


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

Сызранцев В.Н., Сызранцева К.В., Ильиных В.Н. АЛГОРИТМ РАСЧЕТА ДОВЕРИТЕЛЬНЫХ ИНТЕРВАЛОВ КРИВОЙ МАЛОЦИКЛОВОЙ УСТАЛОСТИ // Фундаментальные исследования. – 2015. – № 2-15. – С. 3289-3293;
URL: https://fundamental-research.ru/ru/article/view?id=37770 (дата обращения: 28.03.2024).

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

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