Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,749

ALGORITHM FOR CONFIDENCE INTERVALS CALCULATION OF LOW-CYCLE FATIGUE CURVE

Syzrantsev V.N. 1 Syzrantseva K.V. 1 Ilinykh V.N. 1
1 Federal State Budget Higher Professional Educational Institution «Tyumen State Oil and Gas University»
For forecasting of stress-cycled machine parts lifetime with specified probability of non-destruction authors used developed in kinetic theory of fatigue mathematical model. The results processing of samples fatigue tests is carrying out on a base of this model. Paper illustrates, that two parameters of these model: tensile strength and number of cycles before upper inflection point of low-cycle fatigue curve, related with tensile strength by nonlinear dependence, are random values with unknown distribution laws. In that way left limits of confidence intervals of low-cycle fatigue curve, corresponding to specified probability (1–5 %) of samples non-destruction used for forecasting of its lifetime are impossible to calculate. Authors propose the original algorithm for model parameters determination of fatigue curve passing through any experimental point. Using this algorithm it is possible to calculate during computer modeling the set of fatigue curves. Obtained set of fatigue curves allows to generate the random sample – the number of cycles before destruction, for any fixed level of actual stress. This random sample is necessary for calculation of its quantile estimation at plotting of limits of confidence intervals. For regeneration of unknown density function of distribution this random value authors use mathematical apparatus of nonparametric statistics ensuring the task decision independently of complexity distribution law of researched random value. The results of developed algorithm realization are illustrated on example of limits of confidence intervals determination for results of low-cycle fatigue test of flexible pipes samples.
cyclic deforming
low-cycle fatigue
tensile strength
probability of samples destruction
methods of nonparametric statistics
limits of confidence intervals
1. Stepnov M.N. Statisticheskie metody obrabotki rezultatov mekhanicheskikh ispytaniy [Statistical methods of mechanical test data processing] Moscow, Mashinostroenie, 1985, 232 p.
2. Pochtenny E.K. Kineticheskaya teoriya mekhanicheskoy ustalosti i ee prilozheniya [Kinetic theory of mechanical fatigue and its application], Minsk, Science and engineering, 1973, 213 p.
3. Syzrantsev V.N., Novoselov V.V., Sozonov P.M., Golofast S.L. Ocenka bezopasnosti i prochnostnoy nadezhnosti magistralnykh truboprovodov metodami neparametricheskoy statistiki [Estimation of safety and durability of trunk pipeline by the methods of nonparametric statistic] Novosibirsk, Nauka, 2013, 172 p.
4. Syzrantsev V.N., Nevelev Ja.P., Golofast S.L. Raschet prochnostnoy nadezhnosti izdeliy na osnove metodov neparametricheskoy statistiki [Durability estimation of production based on methods of nonparametric statistic] Novosibirsk, Nauka, 2008, 218 p.
5. Syzrantseva K.V. Raschet prochnostnoy nadezhnosti detaley mashin pri sluchaynom kharaktere vneshnikh nagruzok [Calculation of machine parts durability at random loading] Tyumen: TyumGNGU, 2011, 88 p.
6. Syzrantsev V. N., Syzrantseva K.V., Ilinykh V.N. Fundamentalnye issledovaniya, 2015, no3, pp. 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.