МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ
Расчеты и испытания на прочность МЕТОДЫ
СХЕМАТИЗАЦИИ СЛУЧАЙНЫХ ПРОЦЕССОВ Strength
calculation and testing. Representation |
ГОСТ |
Постановлением Государственного комитета СССР по стандартам от 12 июля 1983 г. № 3066 дата введения установлена
Ограничение срока действия снято по протоколу № 3-93 Межгосударственного совета по стандартизации, метрологии и сертификации (ИУС 5-6-93)
Настоящий стандарт устанавливает методы схематизации случайных процессов нагружения, возникающих в условиях эксплуатации или испытаний элементов машин и конструкций, и статистического представления результатов схематизации.
Результаты схематизации используют для:
расчетной оценки усталостной долговечности элементов по критериям накопления усталостных повреждений на стадиях до появления макроскопической усталостной трещины и ее развития;
сравнительной оценки процессов нагружения однотипных элементов с целью выявления наиболее нагруженных элементов;
количественной оценки эксплуатационных режимов машин и конструкций по их повреждающему воздействию;
моделирования реального нагружения элементов при испытаниях на усталость и расчетном определении характеристик их сопротивления усталости.
Настоящий стандарт не распространяется на методы вычисления корреляционных, автокорреляционных функций и спектральных плотностей процессов реального нагружения, а также методы выбора теоретического закона распределения и проверки согласия с эмпирическим распределением.
1.1. Объектом исследования является реализация случайного процесса нагружения элементов машин или конструкции (изменение во времени нормальных и касательных напряжений, сил, моментов, линейных и угловых деформаций и др.) x(t), далее - процесса нагружения.
1.2. Схематизация процессов нагружения элементов машин и конструкций и статистическое представление результатов схематизации включает следующие этапы:
предварительная подготовка процесса нагружения к схематизации;
дискретизация процесса нагружения (ГОСТ 23207-78);
вычисление статистических характеристик дискретной последовательности ординат процесса нагружения;
выделение экстремумов процесса по дискретной последовательности его ординат;
замена реального процесса нагружения схематизированным по выбранному методу схематизации;
получение эмпирических распределений нагрузок схематизированного процесса и вычисление основных статистических характеристик распределений.
Примечание. В зависимости от применяемых алгоритмов и технических средств допускается изменение последовательности и исключение этапов.
1.3. Схематизацию процессов нагружения проводят применительно к блоку нагружения (ГОСТ 23207-78), который соответствует совокупности последовательных значений переменных нагрузок, возникающих в элементе конструкции за какой-либо характерный период эксплуатации (один полет самолета, lб километров пробега автомобиля по дороге определенного профиля, lб часов работы и т.д.).
1.4. Длительность записи или длина реализации процесса нагружения, соответствующая блоку нагружения, должна содержать не менее 1000 экстремумов пэ процесса.
Примечание. Допускается проводить схематизацию для коротких процессов с числом экстремумов менее 1000. В этом случае рекомендуется оценивать погрешность, вносимую использованием недостаточно представительных по числу экстремумов реализаций.
1.5. При схематизации не учитывают колебания, амплитуда которых меньше половины ширины класса.
1.6. Термины и определения, применяемые в стандарте, соответствуют ГОСТ 23207-78, ГОСТ 21878-76. Условные обозначения приведены в приложении 1.
1.7. Теоретические основы стандарта приведены в приложении 5.
2.1. Подготовка процесса нагружения к схематизации
2.1.1 Разбиение процесса нагружения на классы
Диапазон изменения нагрузок процесса
разбивается равноотстоящими уровнями нагрузок на классы шириной
(2)
Величина D отражает масштаб при схематизации и показывает, какая нагрузка соответствует одному классу. Количество классов т должно удовлетворять неравенству
14 £ т £ 32.
Нумерация классов i производится снизу вверх в направлении возрастания нагрузок процесса, начиная с первого класса в области Xmin. Пример разбиения на классы показан на черт. 1.
Разбиение на классы
Класс, в который попадает значение средней нагрузки , называется классом средней нагрузки и обозначается . Класс, в котором располагается значение нагрузки, равное нулю, обозначается i0.
Если величина нагрузки попадает на границу класса, то ее следует относить к классу с бóльшим номером.
2.1.2. Определение максимальной частоты процесса Fmax
Для определения Fmax разбивают процесс на 20 - 50 равных фрагментов длительностью tф, измеряемой в секундах. В каждом фрагменте подсчитывают количество экстремумов. Соседние экстремумы, образующие размах меньше ширины класса, не подсчитывают. Значение Fmax определяют по формуле
(3)
где пэ - максимальное для рассмотренных фрагментов число экстремумов.
2.1.3. Определение минимальной длины процесса нагружения
Минимальную длину процесса нагружения, необходимую для схематизации, выбирают с учетом п. 1.4, что обеспечивает достаточную статистическую представительность результатов.
Примечание. Допускается выбирать из общей длины процесса нагружения несколько отдельно взятых участков, соответствующих каждому режиму нагружения, входящему в блок нагружения, при условии выполнения п. 1.4 для каждого участка (например, запись процесса нагружения элементов самолета за один полет может быть представлена участками руления, взлета, горизонтального полета и посадки).
2.2. Дискретизация процесса нагружения
Проведению схематизации предшествует дискретное представление процесса нагружения в виде последовательности ординат процесса для ввода цифровой информации в ЭВМ и вычисления статистических характеристик.
Дискретизация может осуществляться методом случайных ординат или методом пересечений.
2.2.1. Дискретизация по методу случайных ординат
На черт. 2 показана дискретизация по методу случайных ординат. Значения процесса нагружения x(f) определяют через равные промежутки времени Dt, в результате получают случайную последовательность ординат процесса нагружения.
Дискретизация по методу случайных ординат
Интервал дискретизации Dt определяют по формуле
(4)
2.2.2. Дискретизация по методу пересечений
При дискретизации по методу пересечений ординаты процесса нагружения определяют в момент пересечения верхней границы класса восходящей или нисходящей ветвями процесса. На черт. 3 иллюстрируется дискретизация по восходящим ветвям процесса.
Примечание. Допускается использовать метод пересечений для схематизации процесса нагружения (приложение 2).
Дискретизация по методу пересечений
2.3. Оценка статистических характеристик процесса нагружения
Для выбора метода схематизации, ее проведения и общей характеристики процесса нагружения вычисляют ряд статистических характеристик по последовательности ординат процесса: среднеарифметическое значение процесса нагружения (средняя нагрузка) , вычисляют как выборочное среднее;
дисперсию процесса нагружения S2, вычисляют как выборочную дисперсию;
среднеквадратическое отклонение процесса нагружения S, вычисляют как выборочное среднеквадратическое отклонение;
коэффициент нерегулярности процесса нагружения k вычисляют согласно ГОСТ 23207-78 по формуле
(5)
где n0 - число пересечений процессом уровня средней нагрузки;
пэ - число экстремумов того же процесса.
Методы схематизации устанавливают правила выделения цикла или полуцикла регулярного нагружения из исследуемого процесса нагружения.
Полуцикл нагружения составляет часть цикла регулярного нагружения (ГОСТ 23207-78), которая соответствует половине периода изменения нагрузки.
Схематизированная реализация случайного процесса нагружения (ГОСТ 23207-78), далее - схематизированный процесс, представляет собой совокупность выделенных по одному из методов схематизации циклов или полуциклов нагружения.
Последовательность нагрузок в схематизированном процессе может быть представлена:
при однопараметрической схематизации - последовательностью амплитуд или размахов выделенных циклов aj;
при двухпараметрической схематизации - корреляционной таблицей (п. 3.3.2) или последовательностью приведенных к симметричному циклу нагружения амплитуд . Полученную по методу схематизации последовательность нагрузок группируют по интервалам, образуя эмпирическое распределение частот повторения нагрузок схематизированного процесса h(xs) (п. 5.1).
3.1. Выделение экстремумов
Непосредственно перед схематизацией по любому из методов должно проводиться выделение экстремумов. Соседние экстремумы, образующие размах меньше ширины класса, выделению не подлежат.
3.2. Однопараметрические методы схематизации
3.2.1. Метод экстремумов
По данному методу рассматривают все положительные максимумы и все отрицательные минимумы (ГОСТ 23207-78) процесса нагружения. На черт. 4 иллюстрируется выделение амплитуд по методу экстремумов. За амплитуды нагружения принимают модуль разности каждого из указанных экстремумов и медианой экстремумов процесса нагружения x50 (ГОСТ 23207-78)
(6)
Схематизация по методу экстремумов
По последовательности значений определяют распределение частот повторения амплитуд полуциклов нагружения h(xa).
По данному методу рассматривают все положительные максимумы. На черт. 5 иллюстрируется схематизация по методу максимумов. За амплитуды нагружения принимают разность
(7)
По последовательности значений определяют распределение частот повторения амплитуд циклов h(xa).
Схематизация по методу максимумов
По данному методу рассматривают все отрицательные минимумы. На черт. 6 иллюстрируется схематизация по методу минимумов. За амплитуду принимают разность
(8)
По последовательности значений определяют распределение частот повторения амплитуд циклов h(xa).
Схематизация по методу минимумов
3.2.4. Однопараметрический метод размахов
На черт. 7 иллюстрируется схематизация по однопараметрическому методу размахов. За размах принимают абсолютное значение разности следующих друг за другом экстремумов
(9)
Схематизация по однопараметрическому методу размахов
Учитывают размахи как на восходящих, так и на нисходящих ветвях процесса. Амплитуды полуцикла нагружения определяют как половину размаха
(10)
Примечание. Допускается проводить схематизацию только по восходящим или по нисходящим ветвям процесса.
3.3. Двухпараметрические методы схематизации
При двухпараметрической схематизации результат представляют либо в виде корреляционной таблицы, либо в виде распределения приведенных амплитуд.
3.3.1. Двухпараметрический метод размахов
Схематизацию по этому методу проводят аналогично однопараметрическому методу размахов, однако учитывают не только величину размаха, но и положение размаха относительно нагрузки, равной нулю. На черт. 8 показана схематизация по двухпараметрическому методу размахов. В результате схематизации заполняют корреляционную таблицу (п. 3.3.2). Схематизацию проводят по восходящим и нисходящим ветвям процесса нагружения.
Примечание. Допускается проводить схематизацию только по восходящим или по нисходящим ветвям процесса.
Схематизация по двухпараметрическому методу размахов
Корреляционная таблица является одной из форм представления результата при двухпараметрической схематизации. В каждой ее клетке содержится накопленное число выделенных циклов или полуциклов нагружения. Величина интервала разбиения каждого параметра корреляционной таблицы равна или кратна ширине класса D.
Допускается использовать две формы корреляционной таблицы:
корреляционная таблица максимумов и минимумов (параметрами таблицы являются максимумы и минимумы выделенных циклов или полуциклов);
корреляционная таблица амплитуд и средних значений (параметрами таблицы являются амплитуда и среднее выделенных циклов или полуциклов).
При заполнении корреляционной таблицы указывают, построена она для циклов или для полуциклов.
3.3.2.1. Корреляционная таблица максимумов и минимумов.
Корреляционная таблица максимумов и минимумов показана в табл. 1. По вертикали сверху вниз и по горизонтали слева направо указывают номера классов.
Регистрацию j-го цикла (полуцикла) проводят в клетке корреляционной таблицы с координатами: по вертикали - номер класса (), в которой попадает максимум; по горизонтали - номер класса (), в который попадает минимум. Табл. 1 заполнена для полуциклов процесса нагружения, изображенного на черт. 8 при схематизации его по двухпараметрическому методу размахов. Так, например, полуцикл с размахом а1 попадает в заштрихованную в табл. 1 клетку, поскольку лежит в классе № 6, а - в классе № 4. В данном примере ширина класса и интервала совпадает.
Таблица 1
Корреляционная таблица максимумов и минимумов
На основании корреляционной таблицы максимумов и минимумов могут быть получены три однопараметрических распределения:
суммирование содержимого клеток в таблице в направлении строк дает распределение минимумов h(xmin);
суммирование в направлении столбцов дает распределение максимумов h(xmax);
суммирование в направлении диагоналей, указанных в табл. 1, дает распределение амплитуд h(xa).
3.3.2.2. Корреляционная таблица амплитуд и средних значений.
Амплитуду и среднее значение для j-го цикла (полуцикла) в единицах класса определяют по формулам:
(11)
(12)
Регистрацию j-го цикла (полуцикла) проводят в клетке корреляционной таблицы с координатами: по вертикали ; по горизонтали .
Корреляционная таблица амплитуд и средних значений показана в табл. 2. Табл. 2 заполнена для полуциклов процесса нагружения, изображенного на черт. 8 при схематизации его по двухпараметрическому методу размахов. Например, полуцикл с размахом a10 попадает в заштрихованную клетку таблицы, поскольку амплитуда полуцикла , а среднее цикла .
На основании корреляционной таблицы амплитуда и средних значений могут быть получены два однопараметрических распределения:
суммирование по строкам таблицы дает распределение амплитуд h(ха);
суммирование по столбцам таблицы дает распределение средних значений циклов (полуциклов) h(xm).
Таблица 2
Корреляционная таблица амплитуд и средних значений
3.3.3. Методы выделения полных циклов
К методам выделения полных циклов относят:
метод полных циклов;
метод «дождя».
Оба метода приводят к практически одинаковым результатам.
Допускается применять методы выделения полных циклов и при однопараметрической схематизации.
3.3.3.1. Метод полных циклов
По данному методу осуществляют выделение циклов постепенно в несколько просмотров, начиная с циклов с наименьшим размахом. На черт. 9 приведена схематизация по методу полных циклов. Исходный процесс нагружения показан на черт. 9а. Здесь вначале выделяют заштрихованные циклы с размахом D £ а < 2D, которые регистрируют и из дальнейшего рассмотрения исключают. Процесс x1(t), изображенный на черт. 9б, получен из исходного х(t) после исключения циклов с D £ а < 2D. Затем выделяют циклы с 2D £ а < 3D (черт. 9б). Процесс х2(t), изображенный на черт. 9в, получен из процесса х1(t) после исключения циклов с размахом 2D £ а < 3D. После исключения циклов с размахами а £ 7D будет получен процесс х3(t), изображенный на черт. 9г. Процесс х3(t) представляет собой один цикл, образованный Хmax и Хmin с наибольшим размахом аmax.
Схематизация по методу полных циклов
3.3.3.2. Метод «дождя»
Схематизация по методу «дождя» приведена на черт. 10. При схематизации удобно представить, что ось времени направлена вертикально вниз. Пусть линии, соединяющие соседние экстремумы - это последовательность крыш, по которым стекают потоки дождя. Номерам максимумов соответствуют нечетные числа, номерам минимумов - четные.
Схематизация по методу «дождя»
Траектории потоков определяют в соответствии со следующими правилами:
1. Потоки начинаются с внутренней стороны экстремумов последовательно. Каждый поток определяет полуцикл нагружения. Величину размаха определяют проекцией траектории потока на ось нагрузки.
2. Поток, начавшийся в точке максимума, прерывается в тот момент, когда встретится максимум больший, чем исходный. Например, поток из максимума 1 стекает по направлению к максимуму 2 и прерывается напротив максимума 3, поскольку он больше исходного.
3. Поток, начавшийся в точке минимума, прерывается, когда встретится минимум меньший, чем исходный. Например, поток из точки 10 прерывается напротив 16, поскольку минимум 16 меньше исходного.
4. При встрече на одной из крыш нескольких потоков движение продолжает тот, который берет начало в экстремуме с меньшим номером, а остальные прерываются. Например, поток из точки 5 продолжает свой путь, а потоки из точек 7 и 9 прерываются.
5. Поток, не встретивший препятствий, падает на землю, например поток из максимума 19.
Метод «дождя» позволяет ввести схематизацию в режиме реального времени. При этом запоминают лишь траектории потоков дождя и обработку проводят по мере ввода экстремумов.
3.4. Приведение асимметричных циклов нагружения к эквивалентным симметричным
Для учета влияния асимметрии цикла нагружения при определении характеристик сопротивления усталости элементов машин и конструкций полученные в результате схематизации корреляционные таблицы преобразуют в однопараметрические таблицы распределения эквивалентных (в смысле вносимого усталостного повреждения) амплитуд. Приведение осуществляют по формуле
(13)
где j - коэффициент чувствительности к асимметрии цикла нагружения.
Коэффициент чувствительности к асимметрии цикла нагружения для сталей выбирают по ГОСТ 25.504-82. Для других материалов значение коэффициента j должно быть подтверждено экспериментально или теоретически.
В ряде случаев более удобным является получение однопараметрических распределений приведенных амплитуд по последовательности амплитуд и средних нагрузок каждого выделенного цикла нагружения , без использования корреляционных таблиц.
Каждый из методов схематизации, изложенных в разд. 3, имеет свои особенности применения в зависимости от характера подлежащего схематизации процесса нагружения. Для обоснованного выбора метода схематизации в пп. 4.1 - 4.7 кратко изложены особенности каждого из них.
4.1. При использовании методов экстремумов, максимумов и минимумов (пп. 3.2.1; 3.2.2; 3.2.3) не учитывают отрицательные максимумы и положительные минимумы. Это приводит к заниженному количеству зарегистрированных циклов нагружения по сравнению с количеством циклов, образуемым всеми максимумами и минимумами процесса. В целом схематизированный процесс по методу экстремумов будет более повреждающим, чем реальный, поскольку мелкие наложенные циклы нагружения учитывают с завышенной амплитудой, например амплитуда процесса нагружения, изображенного на черт. 4.
4.2. Методы максимумов и минимумов дают удовлетворительные результаты для процессов с симметричным относительно средней нагрузки распределением экстремумов.
4.3. При использовании однопараметрического метода размахов (п. 3.2.2) не учитывают средние значения каждого выделенного цикла. Поэтому метод может быть использован только для процессов с незначительным изменением среднего значения циклов нагружения или для оценки нагружения элементов машин и конструкций, сопротивление усталости материала которых зависит практически только от переменной составляющей цикла нагружения.
4.4. Методы размахов как однопараметрический, так и двухпараметрический (п. 3.3.1) приводят к схематизированному процессу менее повреждающему, чем реальный, так как не учитывают циклы с большой амплитудой, на которые наложены мелкие циклы нагружения.
4.5. Методы выделения полных циклов (п. 3.3.3) имеют преимущество по сравнению с перечисленными методами, поскольку учитывают как основные, так и наложенные циклы (см. приложение 3).
4.6. Все методы схематизации, изложенные в настоящем стандарте, в равной степени применимы к процессам нагружения с коэффициентом нерегулярности k = 1.
Для меньших значений k возрастает различие в результатах схематизации, полученных с применением разных методов.
4.7. Область применения методов схематизации, изложенных в разд. 3, для процессов нагружения с различными коэффициентами нерегулярности k приведена в табл. 3. Критерием выбора метода схематизации является коэффициент нерегулярности исследуемого процесса нагружения k.
Таблица 3
Область применения методов схематизации
Пункт стандарта |
k процесса нагружения |
|
1. Метод пересечения границ классов |
0,8 £ v £ 1,0 |
|
2. Метод экстремумов |
0,5 £ k £ 1,0 |
|
3. Метод максимумов |
0,5 £ k £ 1,0 |
|
4. Метод минимумов |
0,5 £ k £ 1,0 |
|
5. Однопараметрический метод размахов |
0,8 £ k £ 1,0 |
|
6. Двухпараметрический метод размахов |
0,5 £ k £ 1,0 |
|
7. Методы выделения полных циклов |
0 < k £ 1,0 |
5.1. Функции распределения
Для оценки нагруженности и сопротивления усталости элементов машин и конструкций используют следующие функции распределения:
h(xs) - распределение частот повторения xs;
H(xs) - распределение накопленных частот повторения xs;
Fэ(xs) - эмпирическая функция распределения;
F(xs) - сглаженная эмпирическая функция распределения;
f(хs) - эмпирическая функция плотности распределения.
Перечисленные функции распределения могут строиться для амплитуд ха, приведенных амплитуд, размахов a, максимумов хmax и минимумов хmin выделенных циклов. Функции распределения представляют в виде таблиц с числом интервалов т. Ширину интервала группирования случайной величины принимают равной или кратной ширине класса D. В таблицах распределений h(хs) и H(xs) указывают в циклах или полуциклах полученные частоты распределения. Примеры построения рядов распределений, графическая интерпретация функции Fэ(xs) и функции плотности распределения fэ(xs) приведены в приложении 4.
5.2. Расчет статистических характеристик распределений
Расчет статистических параметров распределений случайных величин, полученных в результате схематизации, проводят по формулам:
среднее выборочное значение случайной величины
(14)
выборочная дисперсия
(15)
выборочное среднеквадратическое отклонение
(16)
коэффициент вариации
(17)
где hsk - значение h(xs) в k-м интервале.
x(t) - функция, описывающая реализацию случайного процесса нагружения;
t - время;
lб - длительность периода работы машины или конструкции, соответствующая одному блоку нагружения, измеряемая в часах работы, километрах пробега, количестве полетов и т.п.;
пэ - количество экстремумов в процессе нагружения;
amax - диапазон изменения нагрузок процесса;
Xmin - абсолютный минимум процесса нагружения;
Xmax - абсолютный максимум процесса нагружения;
D - ширина класса;
i - номер класса, i = 1, 2, ..., т;
т - количество классов или интервалов;
- среднеарифметическое значение процесса нагружения (средняя нагрузка);
- класс средней нагрузки;
х50 - медиана экстремумов процесса нагружения;
- класс нагрузки, соответствующий медиане экстремумов;
i0 - класс, в котором располагается значение нагрузки, равное нулю;
Fmax - наибольшая частота процесса нагружения;
tф - длительность одного фрагмента, при определении Fmax;
Dt - интервал дискретизации процесса нагружения по времени;
xl - ордината процесса нагружения, получаемая при дискретизации;
S2 - дисперсия нагрузок процесса нагружения;
S - среднеквадратическое отклонение нагрузок процесса нагружения;
k - коэффициент нерегулярности;
n0 - число пересечений процессом уровня средней нагрузки;
- значение нагрузок схематизированного процесса. Под xs подразумевается одна из характеристик выделенных циклов: ха, хт, хmax, хmin, хпр;
- амплитуда схематизированного процесса;
аj - размах схематизированного процесса;
- амплитуда схематизированного процесса, приведенная к симметричному циклу;
h(xs) - распределение частот повторения xs;
- экстремум процесса нагружения;
- минимальное значение нагрузки в j-м цикле нагружения;
- максимальное значение нагрузки в j-м цикле нагружения;
- среднее значение нагрузки в j-м цикле нагружения;
- номер интервала корреляционной таблицы, в который попадает максимум;
- номер интервала корреляционной таблицы, в который попадает минимум;
- номер интервала корреляционной таблицы, в который попадает амплитуда;
- номер интервала корреляционной таблицы, в который попадает среднее значение цикла;
y - коэффициент чувствительности к асимметрии цикла нагружения;
Н(xs) - распределение накопленных частот повторения xs;
Fэ(xs) - эмпирическая функция распределения xs;
F(xs) - сглаженная эмпирическая функция распределения xs;
f(xs) - эмпирическая функция плотности распределения xs;
hsk, Hsk, Fэk, Fk, fsk - значения функций h(xs), H(xs), Fэ(xs), F(xs), f(xs) соответственно в интервале k;
- выборочное среднее значение нагрузки xs;
- выборочная дисперсия нагрузок xs;
Ss - выборочное среднеквадратическое отклонение xs;
vs - коэффициент вариации;
hвi - число пересечений верхней границы i-го класса;
xвsk - верхняя граница k-го интервала нагрузок схематизированного процесса;
хнsk - нижняя граница k-го интервала нагрузок схематизированного процесса;
x¢sk - середина k-го интервала нагрузок схематизированного процесса;
nб - общее число циклов в блоке нагружения;
e - деформация.
Для числа пересечений отдельных уровней и числа максимумов и минимумов в классе i справедливо соотношение
(1)
где hвi - число пересечений верхней границы i-го класса восходящими или нисходящими ветвями процесса.
Для того, чтобы воспользоваться уравнением (1), принимают следующие допущения:
если
Зная распределения hвi с использованием соотношений (1) и (2), строят распределение h(хmax) и h(xmin), используя для выделения амплитуд принцип метода экстремумов (п. 3.2.1), и на основании распределения h(хmax) и h(xmin) строят распределение амплитуд h(xa)
Экстремумы, попадающие в класс , не учитывают при построении распределения h(хa), так как для них амплитуда меньше половины ширины класса D.
Для процесса, изображенного на чертеже, распределения hвi, h(хmax), h(xmin) и h(xa) приведены в таблице.
Схематизация процесса нагружения по методу пересечений
Распределения, полученные по методу пересечений
hвi |
hвi - hвi-1 |
h(хmaxi) |
h(xmini) |
|
|
1 |
1 |
1 |
- |
1 |
1 |
2 |
2 |
1 |
- |
1 |
2 |
3 |
4 |
2 |
- |
2 |
4 |
4 |
5 |
1 |
- |
1 |
2 |
5 |
5 |
0 |
- |
0 |
2 |
|
6 |
1 |
- |
- |
- |
7 |
5 |
1 |
1 |
- |
- |
8 |
4 |
2 |
1 |
- |
- |
9 |
2 |
1 |
2 |
- |
- |
10 |
1 |
1 |
1 |
- |
- |
11 |
0 |
- |
1 |
- |
- |
12 |
- |
- |
- |
- |
- |
Для процесса, изображенного на черт. 1, коэффициент нерегулярности k = 0,6. Для процессов, имеющих такой коэффициент нерегулярности, допускается применять любой из изложенных методов схематизации (см. табл. 3 настоящего стандарта).
В таблице приведены распределения, полученные при однопараметрической схематизации процесса разными методами. На черт. 2 приведена графическая интерпретация распределений. Следует отметить, что для данного процесса имеет место небольшое расхождение результатов схематизации по разным методам.
Распределения амплитуд, полученных по разным методам схематизации
k |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
Метод экстремумов |
h(xs) |
9 |
12 |
9 |
7 |
4 |
3 |
1 |
2 |
|
H(xs) |
9 |
21 |
30 |
37 |
41 |
44 |
45 |
47 |
|
Fэ(xs) |
0,181 |
0,436 |
0,628 |
0,777 |
0,862 |
0,925 |
0,947 |
0,989 |
Метод пересечений |
h(xs) |
9 |
8 |
6 |
5 |
3 |
2 |
2 |
1 |
|
H(xs) |
9 |
17 |
23 |
28 |
31 |
33 |
35 |
36 |
|
Fэ(xs) |
0,236 |
0,458 |
0,625 |
0,764 |
0,847 |
0,903 |
0,958 |
0,986 |
Метод размахов |
h(xs) |
19 |
23 |
15 |
1 |
1 |
- |
- |
- |
|
H(xs) |
19 |
42 |
57 |
58 |
59 |
- |
- |
- |
|
Fэ(xs) |
0,314 |
0,703 |
0,958 |
0,974 |
0,9915 |
- |
- |
- |
Метод «дождя» |
h(xs) |
31 |
12 |
7 |
- |
5 |
2 |
1 |
1 |
|
H(xs) |
31 |
43 |
50 |
50 |
55 |
57 |
58 |
59 |
|
Fэ(xs) |
0,517 |
0,720 |
0,839 |
0,839 |
0,924 |
0,958 |
0,974 |
0,991 |
Примечание. Распределения h(xs) и H(xs) приведены для полуциклов.
1 - метод размахов; 2 - метод «дождя»; 3 - метод экстремумов; 4 - метод пересечений
1. Распределение частот повторений h(xs)
Функция h(xs) представлена в табл. 1.
1 |
2 |
3 |
4 |
5 |
|
хнsk + хвsk |
(0 ¸ 1)D |
(1 ¸ 2)D |
(2 ¸ 3)D |
(3 ¸ 4)D |
(4 ¸ 5)D |
х¢sk |
0,5D |
1,5D |
2,5D |
3,5D |
4,5D |
hsk |
13 |
3 |
1 |
2 |
- |
6 |
7 |
8 |
9 |
10 |
|
хнsk + хвsk |
(5 ¸ 6)D |
(6 ¸ 7)D |
(7 ¸ 8)D |
(8 ¸ 9)D |
(9 ¸ 10)D |
х¢sk |
5,5D |
6,5D |
7,5D |
8,5D |
9,5D |
hsk |
1 |
1 |
- |
2 |
1 |
В первой строке таблицы указывают номера интервалов k, во второй - границы интервалов: xвsk - верхняя граница, xнsk - нижняя граница k-го интервала; в третьей строке - середины интервалов - x¢sk; в четвертой - частоты попадания случайной величины в k-й интервал - hsk. Если случайная величина попадает на границу интервала, суммирование осуществляют в интервале с большим номером.
Табл. 1 заполнена для амплитуд полуцикла процесса нагружения, изображенного на черт. 10 настоящего стандарта при схематизации по методу «дождя». Величины интервалов и классов совпадают.
2. Распределение накопленных частот повторений Н(xs)
Вычисление функции Н(xs) по данным табл. 1 представлено в табл. 2. Накопленные частоты вычисляют по формуле
(1)
Общее число циклов в блоке nб вычисляется по формуле
(2)
и составляет 24 цикла (по данным табл. 2).
Таблица 2
Функция накопленных частот
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
х¢sk |
0,5 |
1,5 |
2,5 |
3,5 |
4,5 |
5,5 |
6,5 |
7,5 |
8,5 |
9,5 |
Hsk |
13 |
16 |
17 |
19 |
19 |
20 |
21 |
21 |
23 |
24 |
3. Эмпирическая функция распределения Fэ(xs)
Функция Fэ(xs) по тем же цифровым данным представлена в табл. 3. Значение функции Fэ(xs) определяет относительную частоту события xs < Xs. Для вычисления Fэk используют формулу
(3)
Вычисление эмпирической функции распределения Fэ(xs)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
х¢sk |
0,5D |
1,5D |
2,5D |
3,5D |
4,5D |
5,5D |
6,5D |
7,5D |
8,5D |
9,5D |
Fэk |
0,5210 |
0,6460 |
0,6870 |
0,7710 |
0,7710 |
0,8120 |
0,8540 |
0,8540 |
0,9370 |
0,9792 |
По эмпирической функции распределения с использованием методов математической статистики, при необходимости, проводят подбор теоретического закона распределения.
4. Сглаженная эмпирическая функция распределения F(xs), эмпирическая функция плотности распределения f(xs)
Построение F(xs) показано на чертеже. На нормальную вероятностную бумагу наносят точки (xsk, Fэk) из табл. 3, которые затем соединяют плавной кривой линией, изображающей функцию F(xs). Значения F(xs) приведены в табл. 4. В табл. 4 приведены значения функции плотности распределения f(xs), которую строят как ступенчатую аппроксимацию функции F(xs).
fk = (Fk - Fk-1).
Таблица 4
Сглаженная функция распределения F(xs) и функция плотности f(xs)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
х¢sk |
0,5D |
1,5D |
2,5D |
3,5D |
4,5D |
5,5D |
6,5D |
7,5D |
8,5D |
9,5D |
Fk |
0,5210 |
0,6064 |
0,6627 |
0,7257 |
0,7710 |
0,8023 |
0,8340 |
0,8850 |
0,9370 |
0,9790 |
Fэk |
0,5210 |
0,0854 |
0,0563 |
0,0630 |
0,0453 |
0,0313 |
0,0317 |
0,0510 |
0,0520 |
0,0420 |
Функция распределения F(xs), представленная на нормальной вероятностной бумаге
+ - функция Fэ(xs), ○ - функция F(xs)
Настоящий стандарт разработан на основе стандартов ГДР (12) и СССР (1, 2), а также с учетом новых разработок и опыта промышленности двух стран (3 - 11, 13 - 15). Стандарт является унифицированным стандартом ГДР и СССР.
Рекомендация длительности реализации, как содержащей не менее 1000 экстремумов, основана на опыте ряда отраслей промышленности в СССР и согласуется с соответствующей рекомендацией ГДР, как обеспечивающей надлежащую точность оценок параметров функций распределения амплитуд напряжений.
Из многочисленных известных методов схематизации случайных процессов в стандарте основное внимание уделено группе методов, основанных на выделении полных циклов, а именно, методу полных циклов и методу «дождя». Эти методы дают практически одинаковые результаты и по мнению большинства специалистов являются наиболее обоснованными.
Если метод экстремумов приводит к схематизированному процессу более повреждающему, чем реальный, метод размахов - к менее повреждающему, то методы полных циклов и «дождя» занимают промежуточное положение и дают расчетные оценки долговечности, более соответствующие экспериментальным данным. Сущность методов заключается в том, что меньшие циклы рассматривают как наложенные на плавный ход нагрузки в одном направлении. На чертеже иллюстрируется явление наложения циклов.
Наложение циклов
а - фрагмент процесса нагружения; б - диаграмма «напряжение-деформация»
Для наглядности зависимость «напряжение-деформация» показана непропорциональной. Данный фрагмент образует две замкнутые петли гистерезиса или два цикла. Цикл 2-3-2 наложен на цикл с большей амплитудой 1-4-5. В свою очередь цикл 2-4-5 может быть наложен на еще больший цикл.
Особенность методов такова, что как бы далеко не располагались в реализации Xmax и Хmin, наибольший размах аmах будет выделен.
При приведении асимметричных циклов к симметричным используется коэффициент влияния асимметрии цикла ys. В стандарте ГДР рекомендовалось для всех случаев принимать ys = 0,33. Однако, как показывают опытные данные для сварных и болтовых соединений, для других деталей с резкой концентрацией напряжений это значение является завышенным. Поэтому в настоящем стандарте рекомендовано принимать ys по рекомендации ГОСТ 25.504-82 для стальных деталей и по соответствующим справочникам для других материалов.
Выбор метода схематизации проводят по величине коэффициента нерегулярности k, равного отношению числа нулей процесса к числу экстремумов. Для процессов простой структуры, для которых k » 1, все методы схематизации приводят практически к одинаковым результатам, поэтому для использования может быть рекомендован любой из них.
Для процессов сложной структуры, для которых k « 1, рекомендуются методы полных циклов или «дождя», как обеспечивающие наиболее достоверную оценку уровня нагруженности.
Литература
1. ГОСТ 23207-78 «Сопротивление усталости. Основные термины, определения и обозначения».
2. ГОСТ 23604-79* «Надежность в технике. Статистическая оценка нагруженности машин и механизмов. Методы обработки данных о нагруженности. Общие положения.
* С 1 июля 1986 г. действует методический документ (ИУС 8-86).
3. Когаев В.П. Расчеты на прочность при напряжениях, переменных во времени. М., Машиностроение, 1977, 232 с.
4. Дмитриченко С.С., Полев В.А., Боровик А.П. Автоматизация расчета на ЭВМ долговечности элементов машин при случайном нагружении. «Вестник машиностроения», № 1, 1982, стр. 7 - 11.
5. Трофимов О.Ф., Злобин Б.Н. Статистический анализ измерений случайной нагруженности для оценки накопления усталостных повреждений. «Вестник машиностроения», № 10, 1969, стр. 3 - 5.
6. К.Н. Вöhmе. Beanspruchungskollektive und ihre Erwartenswahrscheinlichkeit als Grundlage zur Berechnung der Ausfallwahrscheinlichkeit von Bauteilen. «JFL-Mitteilungen», H 1/2, 1976.
7. F. Pfeiffer. Untersuchungen zum EinfluB regelloser Beanspruchungsfolgen auf die Bauteilermüdung «JFL- Mitteilungen», H, 1/2, 1976.
8. J. Dathe, D. Meisel und W. Weygandt. Die Schädigungsgrenze bei Kollektivbeanspruchung und ihre Beracksichtigung beim experimentellen Betriebsfestigkeitsnachweis. «JFL-Mitteilungen», H. 1/2, 1976.
9. F. Pfeiffer und K. Frintert. Zangzeitklassiergerät zur Ermittlung von repräsentativen Belastungskollek-tiven. «JFL-Mitteilungen», H. 1/2, 1976.
10. DDR-Standard TGL 36766. Schwingfestigkeit Ermüdungsprüfung von Werkstoffproben.
11. Hänel, G. Wirthgen. Die Berechnung der Dauerfestigkeit nach dem verfahren von Kogaev und Serensen. «JFL-Mitteilungen», H. 3, 1981.
12. DDR-Standard TGL 33787. Schwingfestigkeit Regellose Zeitfunktion Statidtische Auswertung.
13. Frank Lange und Fritz Pfeiffer. Statistisches Verfahren zur Auswertung von Betriebsbeanspruchungen auf der Basis der Bewertung geschlossener Hystereseschleifen («Rain-Flow») «IFL-MITT» (21), 1982, № 1.
14. Karl-Heinz Bohme. Auswertung von regellosen Zeitfunktionen nach geschlossenen Hystereseschleifen («Rain-Flow» mittels Prozeßrechners «IFL-MITT» (21), 1982, № 1.
15. F. Pfeiffer. Vergleichende Schwingfestigkeitsversuche an geschweißten Proben bei Random-und-Block-beanspruchung «IFL-MITT», 1981, № 2.
СОДЕРЖАНИЕ