ГОСТ 25.101-83 Расчеты и испытания на прочность. Методы схематизации случайных процессов нагружения элементов машин и конструкций и статистического представления результатов.
ГОСТ 25.101-83
Группа Т59
МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ
Расчеты и испытания на прочность
МЕТОДЫ СХЕМАТИЗАЦИИ СЛУЧАЙНЫХ ПРОЦЕССОВ НАГРУЖЕНИЯ ЭЛЕМЕНТОВ МАШИН И КОНСТРУКЦИЙ И СТАТИСТИЧЕСКОГО ПРЕДСТАВЛЕНИЯ РЕЗУЛЬТАТОВ
Strength calculation and testing. Representation of random loading of machine elements and structures and statistical evaluation of results
МКС 03.120.30
19.060
ОКСТУ 0025
Дата введения 1984-07-01
Постановлением Государственного комитета СССР по стандартам от 12 июля 1983 г. N 3066 дата введения установлена 01.07.84
Ограничение срока действия снято по протоколу N 3-93 Межгосударственного совета по стандартизации, метрологии и сертификации (ИУС 5-6-93)
Настоящий стандарт устанавливает методы схематизации случайных процессов нагружения, возникающих в условиях эксплуатации или испытаний элементов машин и конструкций, и статистического представления результатов схематизации.
Результаты схематизации используют для:
расчетной оценки усталостной долговечности элементов по критериям накопления усталостных повреждений на стадиях до появления макроскопической усталостной трещины и ее развития;
сравнительной оценки процессов нагружения однотипных элементов с целью выявления наиболее нагруженных элементов;
количественной оценки эксплуатационных режимов машин и конструкций по их повреждающему воздействию;
моделирования реального нагружения элементов при испытаниях на усталость и расчетном определении характеристик их сопротивления усталости.
Настоящий стандарт не распространяется на методы вычисления корреляционных, автокорреляционных функций и спектральных плотностей процессов реального нагружения, а также методы выбора теоретического закона распределения и проверки согласия с эмпирическим распределением.
1. ОБЩИЕ ПОЛОЖЕНИЯ
1.2. Схематизация процессов нагружения элементов машин и конструкций и статистическое представление результатов схематизации включает следующие этапы:
предварительная подготовка процесса нагружения к схематизации;
дискретизация процесса нагружения (ГОСТ 23207-78);
вычисление статистических характеристик дискретной последовательности ординат процесса нагружения;
выделение экстремумов процесса по дискретной последовательности его ординат;
замена реального процесса нагружения схематизированным по выбранному методу схематизации;
получение эмпирических распределений нагрузок схематизированного процесса и вычисление основных статистических характеристик распределений.
Примечание. В зависимости от применяемых алгоритмов и технических средств допускается изменение последовательности и исключение этапов.
Примечание. Допускается проводить схематизацию для коротких процессов с числом экстремумов менее 1000. В этом случае рекомендуется оценивать погрешность, вносимую использованием недостаточно представительных по числу экстремумов реализаций.
1.5. При схематизации не учитывают колебания, амплитуда которых меньше половины ширины класса.
1.6. Термины и определения, применяемые в стандарте, соответствуют ГОСТ 23207-78, ГОСТ 21878-76. Условные обозначения приведены в приложении 1.
1.7. Теоретические основы стандарта приведены в приложении 5.
2. ПОДГОТОВКА, ДИСКРЕТИЗАЦИЯ И ОЦЕНКА СТАТИСТИЧЕСКИХ ХАРАКТЕРИСТИК ПРОЦЕССА НАГРУЖЕНИЯ
2.1. Подготовка процесса нагружения к схематизации
2.1.1. Разбиение процесса нагружения на классы
Диапазон изменения нагрузок процесса
разбивается равноотстоящими уровнями нагрузок на классы шириной
Разбиение на классы
Черт.1
Если величина нагрузки попадает на границу класса, то ее следует относить к классу с большим номером.
2.1.3 Определение минимальной длины процесса нагружения
Минимальную длину процесса нагружения, необходимую для схематизации, выбирают с учетом п.1.4, что обеспечивает достаточную статистическую представительность результатов.
Примечание. Допускается выбирать из общей длины процесса нагружения несколько отдельно взятых участков, соответствующих каждому режиму нагружения, входящему в блок нагружения, при условии выполнения п.1.4 для каждого участка (например, запись процесса нагружения элементов самолета за один полет может быть представлена участками руления, взлета, горизонтального полета и посадки).
2.2. Дискретизация процесса нагружения
Проведению схематизации предшествует дискретное представление процесса нагружения в виде последовательности ординат процесса для ввода цифровой информации в ЭВМ и вычисления статистических характеристик.
Дискретизация может осуществляться методом случайных ординат или методом пересечений.
2.2.1. Дискретизация по методу случайных ординат
Дискретизация по методу случайных ординат
Черт.2
2.2.2. Дискретизация по методу пересечений
При дискретизации по методу пересечений ординаты процесса нагружения определяют в момент пересечения верхней границы класса восходящей или нисходящей ветвями процесса. На черт.3 иллюстрируется дискретизация по восходящим ветвям процесса.
Примечание. Допускается использовать метод пересечений для схематизации процесса нагружения (приложение 2).
Дискретизация по методу пересечений
Черт.3
2.3. Оценка статистических характеристик процесса нагружения
3. МЕТОДЫ СХЕМАТИЗАЦИИ ПРОЦЕССОВ НАГРУЖЕНИЯ
Методы схематизации устанавливают правила выделения цикла или полуцикла регулярного нагружения из исследуемого процесса нагружения.
Полуцикл нагружения составляет часть цикла регулярного нагружения (ГОСТ 23207-78), которая соответствует половине периода изменения нагрузки.
Схематизированная реализация случайного процесса нагружения (ГОСТ 23207-78), далее - схематизированный процесс, представляет собой совокупность выделенных по одному из методов схематизации циклов или полуциклов нагружения.
3.1. Выделение экстремумов
Непосредственно перед схематизацией по любому из методов должно проводиться выделение экстремумов. Соседние экстремумы, образующие размах меньше ширины класса, выделению не подлежат.
3.2. Однопараметрические методы схематизации
3.2.1. Метод экстремумов
Схематизация по методу экстремумов
Черт.4
3.2.2. Метод максимумов
По данному методу рассматривают все положительные максимумы. На черт.5 иллюстрируется схематизация по методу максимумов. За амплитуды нагружения принимают разность
Схематизация по методу максимумов
Черт.5
3.2.3. Метод минимумов
По данному методу рассматривают все отрицательные минимумы. На черт.6 иллюстрируется схематизация по методу минимумов. За амплитуду принимают разность
Схематизация по методу минимумов
Черт.6
3.2.4. Однопараметрический метод размахов
На черт.7 иллюстрируется схематизация по однопараметрическому методу размахов. За размах принимают абсолютное значение разности следующих друг за другом экстремумов
Схематизация по однопараметрическому методу размахов
Черт.7
Учитывают размахи как на восходящих, так и на нисходящих ветвях процесса. Амплитуды полуцикла нагружения определяют как половину размаха
Примечание. Допускается проводить схематизацию только по восходящим или по нисходящим ветвям процесса.
3.3. Двухпараметрические методы схематизации
При двухпараметрической схематизации результат представляют либо в виде корреляционной таблицы, либо в виде распределения приведенных амплитуд.
3.3.1. Двухпараметрический метод размахов
Схематизацию по этому методу проводят аналогично однопараметрическому методу размахов, однако учитывают не только величину размаха, но и положение размаха относительно нагрузки, равной нулю. На черт.8 показана схематизация по двухпараметрическому методу размахов. В результате схематизации заполняют корреляционную таблицу (п.3.3.2). Схематизацию проводят по восходящим и нисходящим ветвям процесса нагружения.
Схематизация по двухпараметрическому методу размахов
Черт.8
Примечание. Допускается проводить схематизацию только по восходящим или по нисходящим ветвям процесса.
3.3.2. Корреляционная таблица
Допускается использовать две формы корреляционной таблицы:
корреляционная таблица максимумов и минимумов (параметрами таблицы являются максимумы и минимумы выделенных циклов или полуциклов);
корреляционная таблица амплитуд и средних значений (параметрами таблицы являются амплитуда и среднее выделенных циклов или полуциклов).
При заполнении корреляционной таблицы указывают, построена она для циклов или для полуциклов.
3.3.2.1. Корреляционная таблица максимумов и минимумов.
Корреляционная таблица максимумов и минимумов показана в табл.1. По вертикали сверху вниз и по горизонтали слева направо указывают номера классов.
Таблица 1
Корреляционная таблица максимумов и минимумов
На основании корреляционной таблицы максимумов и минимумов могут быть получены три однопараметрических распределения:
суммирование в направлении диагоналей, указанных в табл.1, дает распределение амплитуд
3.3.2.2. Корреляционная таблица амплитуд и средних значений.
Таблица 2
Корреляционная таблица амплитуд и средних значений
На основании корреляционной таблицы амплитуд и средних значений могут быть получены два однопараметрических распределения:
3.3.3. Методы выделения полных циклов
К методам выделения полных циклов относят:
метод полных циклов;
метод "дождя".
Оба метода приводят к практически одинаковым результатам.
Допускается применять методы выделения полных циклов и при однопараметрической схематизации.
3.3.3.1. Метод полных циклов
Схематизация по методу полных циклов
Черт.9
3.3.3.2. Метод "дождя"
Схематизация по методу "дождя" приведена на черт.10. При схематизации удобно представить, что ось времени направлена вертикально вниз. Пусть линии, соединяющие соседние экстремумы - это последовательность крыш, по которым стекают потоки дождя. Номерам максимумов соответствуют нечетные числа, номерам минимумов - четные.
Схематизация по методу "дождя"
Черт.10
Траектории потоков определяют в соответствии со следующими правилами:
1. Потоки начинаются с внутренней стороны экстремумов последовательно. Каждый поток определяет полуцикл нагружения. Величину размаха определяют проекцией траектории потока на ось нагрузки.
2. Поток, начавшийся в точке максимума, прерывается в тот момент, когда встретится максимум больший, чем исходный. Например, поток из максимума 1 стекает по направлению к максимуму 2 и прерывается напротив максимума 3, поскольку он больше исходного.
3. Поток, начавшийся в точке минимума, прерывается, когда встретится минимум меньший, чем исходный. Например, поток из точки 10 прерывается напротив 16, поскольку минимум 16 меньше исходного.
4. При встрече на одной из крыш нескольких потоков движение продолжает тот, который берет начало в экстремуме с меньшим номером, а остальные прерываются. Например, поток из точки 5 продолжает свой путь, а потоки из точек 7 и 9 прерываются.
5. Поток, не встретивший препятствий, падает на землю, например, поток из максимума 19.
Метод "дождя" позволяет ввести схематизацию в режиме реального времени. При этом запоминают лишь траектории потоков дождя и обработку проводят по мере ввода экстремумов.
3.4. Приведение асимметричных циклов нагружения к эквивалентным симметричным
Для учета влияния асимметрии цикла нагружения при определении характеристик сопротивления усталости элементов машин и конструкций полученные в результате схематизации корреляционные таблицы преобразуют в однопараметрические таблицы распределения эквивалентных (в смысле вносимого усталостного повреждения) амплитуд. Приведение осуществляют по формуле
4. ВЫБОР МЕТОДА СХЕМАТИЗАЦИИ. КРИТЕРИИ ВЫБОРА
Каждый из методов схематизации, изложенных в разд.3, имеет свои особенности применения в зависимости от характера подлежащего схематизации процесса нагружения. Для обоснованного выбора метода схематизации в пп.4.1-4.7 кратко изложены особенности каждого из них.
4.2. Методы максимумов и минимумов дают удовлетворительные результаты для процессов с симметричным относительно средней нагрузки распределением экстремумов.
4.3. При использовании однопараметрического метода размахов (п.3.2.2) не учитывают средние значения каждого выделенного цикла. Поэтому метод может быть использован только для процессов с незначительным изменением среднего значения циклов нагружения или для оценки нагружения элементов машин и конструкций, сопротивление усталости материала которых зависит практически только от переменной составляющей цикла нагружения.
4.4. Методы размахов как однопараметрический, так и двухпараметрический (п.3.3.1) приводят к схематизированному процессу менее повреждающему, чем реальный, так как не учитывают циклы с большой амплитудой, на которые наложены мелкие циклы нагружения.
4.5. Методы выделения полных циклов (п.3.3.3) имеют преимущество по сравнению с перечисленными методами, поскольку учитывают как основные, так и наложенные циклы (см. приложение 3).
Таблица 3
Область применения методов схематизации
|
|
|
Метод схематизации | Пункт стандарта | процесса нагружения |
1. Метод пересечения границ классов | 2.2.2 | 0,8 1,0 |
2. Метод экстремумов | 3.2.1 | 0,5 1,0 |
3. Метод максимумов | 3.2.2 | 0,5 1,0 |
4. Метод минимумов | 3.2.3 | 0,5 1,0 |
5. Однопараметрический метод размахов | 3.2.4 | 0,8 1,0 |
6. Двухпараметрический метод размахов | 3.3.1 | 0,5 1,0 |
7. Meтоды выделения полных циклов | 3.3.3 | 0 1,0 |
5. СТАТИСТИЧЕСКОЕ ПРЕДСТАВЛЕНИЕ РЕЗУЛЬТАТОВ СХЕМАТИЗАЦИИ
5.1. Функции распределения
Для оценки нагруженности и сопротивления усталости элементов машин и конструкций используют следующие функции распределения:
иложении 4.
5.2. Расчет статистических характеристик распределений
Расчет статистических параметров распределений случайных величин, полученных в результате схематизации, производят по формулам:
среднее выборочное значение случайной величины
выборочная дисперсия
выборочное среднекрадратическое отклонение
коэффициент вариации
ПРИЛОЖЕНИЕ 1
Обязательное
УСЛОВНЫЕ ОБОЗНАЧЕНИЯ, ПРИМЕНЯЕМЫЕ В СТАНДАРТЕ
|
|
- функция, описывающая реализацию случайного процесса нагружения; | |
- время; | |
- длительность периода работы машины или конструкции, соответствующая одному блоку нагружения, измеряемая в часах работы, километрах пробега, количестве полетов и т.п.; | |
- количество экстремумов в процессе нагружения; | |
- диапазон изменения нагрузок процесса; | |
- абсолютный минимум процесса нагружения; | |
- абсолютный максимум процесса нагружения; | |
- ширина класса; | |
- номер класса, =1, 2, ..., ; | |
- количество классов или интервалов; | |
- среднеарифметическое значение процесса нагружения (средняя нагрузка); | |
- класс средней нагрузки; | |
- медиана экстремумов процесса нагружения; | |
- класс нагрузки, соответствующий медиане экстремумов; | |
- класс, в котором располагается значение нагрузки, равное нулю; | |
- наибольшая частота процесса нагружения; | |
- длительность одного фрагмента, при определении ; | |
- интервал дискретизации процесса нагружения по времени; | |
- ордината процесса нагружения, получаемая при дискретизации; | |
- дисперсия нагрузок процесса нагружения; | |
- среднеквадратическое отклонение нагрузок процесса нагружения; | |
- коэффициент нерегулярности; | |
- число пересечений процессом уровня средней нагрузки; | |
- значение нагрузок схематизированного процесса. Под подразумевается одна из характеристик выделенных циклов: , , , , ; | |
- амплитуда схематизированного процесса; | |
- размах схематизированного процесса; | |
- амплитуда схематизированного процесса, приведенная к симметричному циклу; | |
- распределение частот повторения ; | |
- экстремум процесса нагружения; | |
- минимальное значение нагрузки в -м цикле нагружения; | |
- максимальное значение нагрузки в -м цикле нагружения; | |
- среднее значение нагрузки в -м цикле нагружения; | |
- номер интервала корреляционной таблицы, в который попадает максимум; | |
- номер интервала корреляционной таблицы, в который попадает минимум; | |
- номер интервала корреляционной таблицы, в который попадает амплитуда; | |
- номер интервала корреляционной таблицы, в который попадает среднее значение цикла; | |
- коэффициент чувствительности к асимметрии цикла нагружения; | |
- распределение накопленных частот повторения ; | |
- эмпирическая функция распределения ; | |
- сглаженная эмпирическая функция распределения ; | |
- эмпирическая функция плотности распределения ; | |
, , , ,
| - значения функций , , , , соответственно в интервале ; |
- выборочное среднее значение нагрузки ; | |
- выборочная дисперсия нагрузок ; | |
- выборочное среднеквадратическое отклонение ; | |
- коэффициент вариации; | |
- число пересечений верхней границы -го класса; | |
- верхняя граница -го интервала нагрузок схематизированного процесса; | |
- нижняя граница -го интервала нагрузок схематизированного процесса; | |
- середина -го интервала нагрузок схематизированного процесса; | |
- общее число циклов в блоке нагружения; | |
- деформация. |
ПРИЛОЖЕНИЕ 2
Справочное
МЕТОД ПЕРЕСЕЧЕНИЙ ДЛЯ СХЕМАТИЗАЦИИ ПРОЦЕССОВ НАГРУЖЕНИЯ
Для того, чтобы воспользоваться уравнением (1), принимают следующие допущения:
(2)
Схематизация процесса нагружения по методу пересечений
Распределения, полученные по методу пересечений
|
|
|
|
|
|
| |||||
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 | 6 | 1 | - | - | - |
7 | 5 | 1 | 1 | - | - |
8 | 4 | 2 | 1 | - | - |
9 | 2 | 1 | 2 | - | - |
10 | 1 | 1 | 1 | - | - |
11 | 0 | - | 1 | - | - |
12 | - | - | - | - | - |
ПРИЛОЖЕНИЕ 3
Справочное
ПРИМЕР СХЕМАТИЗАЦИИ ПРОЦЕССА НАГРУЖЕНИЯ
В таблице приведены распределения, полученные при однопараметрической схематизации процесса разными методами. На черт.2 приведена графическая интерпретация распределений. Следует отметить, что для данного процесса имеет место небольшое расхождение результатов схематизации по разным методам.
Распределения амплитуд, полученных по разным методам схематизации
|
|
|
|
|
|
|
|
|
|
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | |
Метод экстремумов | 9 | 12 | 9 | 7 | 4 | 3 | 1 | 2 | |
| 9 | 21 | 30 | 37 | 41 | 44 | 45 | 47 | |
| 0,181 | 0,436 | 0,628 | 0,777 | 0,862 | 0,925 | 0,947 | 0,989 | |
Метод пересечений | 9 | 8 | 6 | 5 | 3 | 2 | 2 | 1 | |
| 9 | 17 | 23 | 28 | 31 | 33 | 35 | 36 | |
| 0,236 | 0,458 | 0,625 | 0,764 | 0,847 | 0,903 | 0,958 | 0,986 | |
Метод размахов | 19 | 23 | 15 | 1 | 1 | - | - | - | |
| 19 | 42 | 57 | 58 | 59 | - | - | - | |
| 0,314 | 0,703 | 0,958 | 0,974 | 0,9915 | - | - | - | |
Метод "дождя" | 31 | 12 | 7 | - | 5 | 2 | 1 | 1 | |
| 31 | 43 | 50 | 50 | 55 | 57 | 58 | 59 | |
| 0,517 | 0,720 | 0,839 | 0,839 | 0,924 | 0,958 | 0,974 | 0,991 |
Черт.1
1 - метод размахов; 2 - метод "дождя"; 3 - метод экстремумов; 4 - метод пересечений
Черт.2
ПРИЛОЖЕНИЕ 4
Справочное
ПРИМЕРЫ СТАТИСТИЧЕСКОГО ПРЕДСТАВЛЕНИЯ РЕЗУЛЬТАТОВ СХЕМАТИЗАЦИИ
Таблица 1
|
|
|
|
|
|
1 | 2 | 3 | 4 | 5 | |
(0 1) | (1 2) | (2 3) | (3 4) | (4 5) | |
0,5 | 1,5 | 2,5 | 3,5 | 4,5 | |
13 | 3 | 1 | 2 | - |
Продолжение табл.1
|
|
|
|
|
|
6 | 7 | 8 | 9 | 10 | |
(5 6) | (6 7) | (7 8) | (8 9) | (9 10) | |
5,5 | 6,5 | 7,5 | 8,5 | 9,5 | |
1 | 1 | - | 2 | 1 |
Табл.1 заполнена для амплитуд полуцикла процесса нагружения, изображенного на черт.10 настоящего стандарта при схематизации по методу "дождя". Величины интервалов и классов совпадают.
и составляет 24 цикла (по данным табл.2).
Таблица 2
Функция накопленных частот
|
|
|
|
|
|
|
|
|
|
|
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
0,5 | 1,5 | 2,5 | 3,5 | 4,5 | 5,5 | 6,5 | 7,5 | 8,5 | 9,5 | |
13 | 16 | 17 | 19 | 19 | 20 | 21 | 21 | 23 | 24 |
Таблица 3
|
|
|
|
|
|
|
|
|
|
|
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
0,5 | 1,5 | 2,5 | 3,5 | 4,5 | 5,5 | 6,5 | 7,5 | 8,5 | 9,5 | |
0,5210 | 0,6460 | 0,6870 | 0,7710 | 0,7710 | 0,8120 | 0,8540 | 0,8540 | 0,9370 | 0,9792 |
По эмпирической функции распределения с использованием методов математической статистики, при необходимости, проводят подбор теоретического закона распределения.
Таблица 4
|
|
|
|
|
|
|
|
|
|
|
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | |
0,5 | 1,5 | 2,5 | 3,5 | 4,5 | 5,5 | 6,5 | 7,5 | 8,5 | 9,5 | |
0,5210 | 0,6064 | 0,6627 | 0,7257 | 0,7710 | 0,8023 | 0,8340 | 0,8850 | 0,9370 | 0,9790 | |
0,5210 | 0,0854 | 0,0563 | 0,0630 | 0,0453 | 0,0313 | 0,0317 | 0,0510 | 0,0520 | 0,0420 |
ПРИЛОЖЕНИЕ 5
Справочное
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ СТАНДАРТА
Настоящий стандарт разработан на основе стандартов ГДР (12) и СССР (1, 2), а также с учетом новых разработок и опыта промышленности двух стран (3-11, 13-15). Стандарт является унифицированным стандартом ГДР и СССР.
Рекомендация длительности реализации, как содержащей не менее 1000 экстремумов, основана на опыте ряда отраслей промышленности в СССР и согласуется с соответствующей рекомендацией ГДР, как обеспечивающей надлежащую точность оценок параметров функций распределения амплитуд напряжений.
Из многочисленных известных методов схематизации случайных процессов в стандарте основное внимание уделено группе методов, основанных на выделении полных циклов, а именно, методу полных циклов и методу "дождя". Эти методы дают практически одинаковые результаты и по мнению большинства специалистов являются наиболее обоснованными.
Если метод экстремумов приводит к схематизированному процессу более повреждающему, чем реальный, метод размахов - к менее повреждающему, то методы полных циклов и "дождя" занимают промежуточное положение и дают расчетные оценки долговечности, более соответствующие экспериментальным данным.
Сущность методов заключается в том, что меньшие циклы рассматривают как наложенные на плавный ход нагрузки в одном направлении. На чертеже иллюстрируется явление наложения циклов.
Наложение циклов
а - фрагмент процесса нагружения; б - диаграмма "напряжение-деформация"
Для наглядности зависимость "напряжение-деформация" показана непропорциональной. Данный фрагмент образует две замкнутые петли гистерезиса или два цикла. Цикл 2-3-2 наложен на цикл с большей амплитудой 1-4-5. В свою очередь цикл 2-4-5 может быть наложен на еще больший цикл.
Литература
1. ГОСТ 23207-78 "Сопротивление усталости. Основные термины, определения и обозначения".
2. ГОСТ 23604-79* "Надежность в технике. Статистическая оценка нагруженности машин и механизмов. Методы обработки данных о нагруженности. Общие положения".
________________
* С 1 июля 1986 г. действует методический документ (ИУС 8-86).
3. Когаев В.П. Расчеты на прочность при напряжениях, переменных во времени. М., Машиностроение, 1977, 232 с.
4. Дмитриченко С.С., Полев В.А., Боровик А.П. Автоматизация расчета на ЭВМ долговечности элементов машин при случайном нагружении. "Вестник машиностроения", N 1, 1982, стр. 7-11.
5. Трофимов О.Ф., 3лобин Б.Н. Статистический анализ измерений случайной нагруженности для оценки накопления усталостных повреждений. "Вестник машиностроения", N 10, 1969, стр.3-5.
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, N 1.