О некоторых особенностях планирования и статистической обработки результатов усталостных испытаний образцов Текст научной статьи по специальности «Строительство и архитектура»

О некоторых особенностях планирования и статистической обработки результатов усталостных испытаний образцов

А.В. Кулагин Удмуртский государственный университет, г. Ижевск

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

Некоторые общие теоретические закономерности оценки прочности и надежности в условиях статического и циклического нагружения элементов конструкций и деталей машин достаточно хорошо согласованы с практикой эксплуатации, что нашло отражение в ряде нормативных документов (ГОСТ 25.502-79. Расчеты и испытания на прочность в машиностроении. Методы механических испытаний металлов. Методы испытаний на усталость. РД 50705-91. Методические указания. Расчеты и испытания на прочность в машиностроении. Планирование и статистическая обработка результатов статических испытаний и испытаний на усталость.) и научно-технических изданиях [1,2]. Однако, если при статическом нагружении проведение эксперимента и обработка его результатов не вызывает больших трудностей, то оценка усталостной прочности в силу многофакторности особенностей работы многих технических систем бывает довольно сложной и может привести к достаточно большому разбросу механических характеристик материалов. Поэтому попытаемся систематизировать статистику обработки результатов усталостных испытаний образцов в виде проведения 4-х основных этапов.

¡.Выбор числовых характеристик параметра статической прочности Механические свойства материала определяются по результатам испытаний ограниченного числа образцов (по выборке). Поэтому количественные оценки свойств материала являются статистическими оценками, так называемых, генеральных характеристик и определяются по формулам математической статистики [1-3]. Обозначим:

л - значение характеристики (параметра прочности), полученное в результате испытания /-ого образца;

среднее арифметическое параметра прочности; - выборочное среднее квадратичное отклонение; 5- выборочная дисперсия параметра прочности; 1; - выборочный коэффициент вариации; п -объем выборки (количество образцов в партии).

Для расчета выборочных характеристик применяются следующие

формулы:

Уа= *

< й>

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

\1Ь Л-'Гь

где M(R)- математическое ожидание параметра прочности;

^-уровень значимости;

у=п-1- число степеней свободы;

Ьдда- критерий Стьюдеита для уровня значимости а и числа степеней свободы

V.

При определении доверительного интервала для генерального среднего квадратичного отклонения параметра прочности можно пользоваться выражением

в котором <7(Я) - генеральное среднее квадратичное отклонение параметра прочности;

q -величина, зависящая от уровня доверительной вероятности Р и числа степеней свободы V, р = 1-й.

2. Обработка результатов усталостных испытаний. Методика обработки зависит от того, какова была цель испытаний

[4,5].

2.1. Определение долговечности (ресурса) образцов. Выборочные характеристики долговечности N определяется по формулам (1) с учетом равенства й, = ]£АГС. При этом принимается гипотеза о нормальном распределении случайной величины lgN.

2.2. Определение порога чувствительности по циклам. Считается, что долговечность (или, точнее ее количественная оценка, ресурс) распределена непрерывно в интервале от N = до Л? = Нижняя граница ЛГ0 этого интервала и называется порогом чувствительности по циклам. При определении порога чувствительности исходят из предположения о нормальном законе распределения случайной величины ж = — ЛРц). В соответствии с методом максимума правдоподобия расчет производится с использованием уравнения

;

В основе регрессивного анализа лежит нормальный закон распределения [4]. Целью его является установление количественной связи между максимальным напряжением и долговечностью образца. Уравнение эмпирической линии регрессии имеет вид

где X -оценка условного математического ожидания величины х; х = нормально распределенная случайная величина; у= 1§<гП1ая, - независимая величина; <?шяяг - максимальное напряжение (параметр цикла); <х>,<у>,т -параметры уравнения линии регрессии.

Параметры уравнения (2) рассчитываются по формулам:

где п-число уравнений напряжений.

Полученное таким образом уравнение соответствует осредненной кривой усталости (для вероятности разрушения P=0,5). Для построения доверительной области теоретической линии регрессии вычисляют: дисперсию вокруг эмпирической лини регрессии

X =<х > < у >1

(2)

т =

;

;

,

;

дисперсию параметров m:

дисперсию оценки математического ожидания для уровня у

I <у>)3 (3)

Границы области для доверительной вероятности рассчитывают по формуле

, (4)

где а=1-р - уровень значимости;

v=n-2 - число степеней свободы (здесь п-объем выборки); - критерий Стьюдента;

Расчеты по формулам (3) и (4) производят отдельно для каждого уровня напряжений.

4. Определение параметров функции распределения предела выносливости (метод ступенчатого изменения нагрузки).

Метод предусматривает испытание серии образцов на нескольких уровнях напряжения (ГОСТ 25.502-79. Расчеты и испытания на прочность в машиностроении. Методы механических испытаний металлов. Методы испытаний на усталость. РД 50-705-91. Методические указания. Расчеты и испытания на прочность в машиностроении. Планирование и статистическая обработка результатов статических испытаний и испытаний на усталость.)

Первый образец испытывают при напряжении, равном ожидаемому пределу выносливости, а уровень напряжения для каждого последующего образца повышается по сравнению с предыдущим образцом, если тот не разрушился, и понижается, если предыдущий образец разрушился. Перепад а между уровнями напряжений выбирается постоянным и не превышающим двух значений среднего квадратичного отклонения предела выносливости при любом выбранном методе схематизации процессов нагружения, то есть на рисунке 1а) это соответствует разбросу по экстремумам напряжений d4+d7, а на рисунке 1б) для этих же напряжений-это размах 2 А2 и так далее,

:

сравнивая попарно все пиковые напряжения при проведении однотипных испытаний и 8 измерениях.

Рис.1- Оценка уровня усталостных напряжений во времени После испытания подсчитывают общее число разрушений и общее число не разрушений и выбирают Д^^- наименьшее из них. Пусть менее частое событие происходило при напряжениях <73,0^,... с

соответствующими частотами п^,п1г...гпг, ... • Тогда оценка среднего значения предела выносливости производится по формуле

в которой знак плюс (+) ставится, когда число не разрушений, а знак минус (-) когда число разрушений.

Среднее квадратичное отклонение предела выносливости определяется в виде

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

Основная задача планирования испытания заключается в определении объема выборки. Если целью испытания является оценка характеристик статической прочности или усталостной долговечности при уровне

значимости а, то необходимое количество образцов вычисляют по формулам:

= —

' -1 „с : _ ^

'ч1 к- -

где % число образцов среднего значения характеристики прочности с предельной относительной ошибкой

У - коэффициент вариации генеральной совокупности (задаются на

основании уже имеющихся данных по аналогичным материалам);

ир - квантиль нормального распределения вероятностей Р=1-а/2;

■г.. - предельная относительная ошибка среднего квадратичного отклонения

характеристики прочности;

^■й/гг^оа" критерий Пирсона для вероятностей Р= ^ и Р=0,5, соответственно,

и числа степеней свободы у=п-1;

- число образцов, необходимое для определения среднего квадратичного

отклонения с ошибкой Еа, находится в результате решения второго

уравнения методом последовательных приближений.

Из двух значений объема испытаний пм и и0 выбирается наибольшее

значение. Если в результате испытания дается оценка новой технологии или

I Инженерный вестник Дона, №4 (2019)

ivdon.ru/ru/magazine/archive/n4y2019/5858

другого фактора, изменяющего прочностные свойства материала, то формулы имеют вид:

где а - уровень значимости, представляющий собой достаточно малую вероятность отрицательной оценки более эффективной технологии (ошибка первого рода);

в - достаточно малая вероятность положительной оценки менее эффективной технологии (ошибка второго рода).

Значения ошибок рекомендуется выбирать в пределах ем = 0,01 - 0,05;

■Г- = 0,1 - 0,5; а вероятности ошибок первого и второго рода - а=0,05 - 0,1; в

Для построения осредненной кривой усталости рекомендуется проводить испытания на 4 - 5 уровнях напряжения. Максимальный уровень выбирается с учетом требований на протяженность левой ветви кривой усталости, а минимальный подбирается опытным путем так, чтобы долговечность меньше базового числа циклов показали не более, чем 50% образцов, разрушенных при напряжении этого уровня. Образцы по уровням напряжений следует распределять равномерно. Общее количество образцов, необходимое для построения осредненной кривой усталости при относительной предельной ошибке е предела выносливости и уровне значимости а, находится по формуле

;

а

0,05 - 0,1.

,

где К- число параметров уравнения линии регрессии.

Формулой можно воспользоваться и при расчете количества образцов, необходимого для определения медианного значения предела выносливости с заданной точностью и надежностью. При этом кривая усталости описывается двухпараметрическим уравнением (К = 2).

Далее, используя результаты 4-х этапов можно, например, перейти к решению задачи по единому критерию усталостной прочности, уточнению параметров внешних циклических воздействий и рационального выбора материалов, что потребует экспериментальной проверки [7-10].

Литература

1.Степнов М.Н., Шаврин А.В. Статистические методы обработки результатов механических испытаний. Справочник. - М.: Машиностроение, 2005.- 399 с. -2-е издание исправленное и дополненное.

2.Александров А.В., Потапов В.Д., Державин Б.П. Сопротивление материалов: Учеб. для вузов/ Под ред. Александрова А.В. — 3-е издание исправленное. М.: Высшая школа, 2003. 560 с: ил.

3.Хальд А. Математическая статистика с техническим приложением. М.: Изд. Иностр. лит., 1956. 664 с.

4.Дрейпер Н., Смит Г. Прикладной регрессионный анализ. - М.: Вильямс, 2016. - 912 с.

5.Протодьяконов М.М., Тедер Р.И. Методика рационального планирования экспериментов. - Л.: Наука, 2013. - 338 c.

6. Болотин, Владимир Васильевич. Ресурс машин и конструкций. М.: Машиностроение, 1990. - 446 с

7. Davidson D.L., Lankford J. Fatigue crack growth in metals and alloys: mechanisms and micromechanics // International Materials Reviews. - 1992 - V. 37- № 2 - pp. 45-76.

8. Mughrabi H., Christ H.-J. Cyclic deformation and fatigue of selected ferritic and austeniticsteels; specific aspects // ISIJ International. - 1997 - V. 37 -№ 12 - pp. 1154-1169.

9.Кулагин А.В., Дородов П.В. О запасе прочности и оценке надежности узлов металлоконструкций // Инженерный вестник Дона, 2012. - № 2 URL: ivdon.ru/ru/magazine/archive/n2y2012/810.

10.Кулагин А.В., Дородов П.В. Исследование напряжений возле плоского горизонтального выреза // Инженерный вестник Дона. - 2012. - № 2 URL: ivdon.ru/ru/magazine/archive/n2y2012/813.

References

1. Stepnov M. N., Shavrin A.V. Statisticheskie metody obrabotki rezul'tatov mekhanicheskikh ispytaniy [Statistical methods of processing the results of mechanical tests]. HandbookMoscow: Mechanical Engineering, 2005. 399 p. 2nd edition, revised and supplemented.

2. Aleksandrov A.V., Potapov V. D., Derzhavin B. P. Soprotivlenie materialov [The mechanics of materials]. Proc. for universities. Aleksandrov A.V., Potapov, Derzhavin; under the editorship of Aleksandrova A. V. 3rd ed. ISPR. M.: Higher. SHK, 2003. 560 p: Il.

3. Hald A. Matematicheskaya statistika s tekhnicheskim prilozheniem [Mathematical statistics with technical application]. M.: Izd. Foreign. lit., 1956. 664 p.

4. Draper N., Smith G. Prikladnoy regressionnyy analiz [Applied regression analysis]. M.: Williams, 2016. 912 p.

5. Protodyakonov M. M., Teder R.I. Metodika ratsional'nogo proektirovaniya eksperimentov [Methods of rational planning of experiments]. L.: Science, 2013. 338 p.

6. Bolotin V.V. Resurs mashin i konstrukcij [Resource of machines and structures]. M.: Mechanical Engineering, 1990. 446 p.

II Инженерный вестник Дона, №4 (2019) ivdon.ru/ru/magazine/archive/n4y2019/5858

7. Davidson D.L., Lankford J. International Materials Reviews. 1992. V. 37. № 2. pp. 45-76.

8. Mughrabi H., Christ H.-J. ISIJ International. 1997 V. 37. № 12. pp. 11541169.

9. Kulagin A.V., Dorodov P.V. Inzenernyj vestnik Dona (Rus), 2012. № 2. URL: ivdon.ru/ru/magazine/archive/n2y2012/810.

10. Kulagin A.V. Dorodov, P. V. Inzenernyj vestnik Dona (Rus), 2012. № 2 URL: ivdon.ru/ru/magazine/archive/n2y2012/813.